【正文】 動變形值，各種變形曲線依不同比例在 CAD 中成圖，且對應(yīng)繪出地下煤層采動情況。鄧喀中提出了巖體開采沉陷的結(jié)構(gòu)效應(yīng)。國外很多國家用基于概率基本法 [7]預(yù)計計算來指導(dǎo)建筑物下、水體下和鐵路下采煤。在 1838 年 ,多里斯提出了開采沉陷的第一個理論 “ 垂線理論 ” 。為了有效預(yù)防和治理這類沉陷災(zāi)害，需要獲取地表移動規(guī)律。在輸出下沉、傾斜、曲率、水平移動、水平變形計算結(jié)果的同時，將其各曲線以一定比例在 CAD 中成圖，可從曲線上量取任意點的對應(yīng)值。目前是最直接有效的方法是建立地表移動觀測站，定期對觀測站上的監(jiān)測點及控制點進行觀測。 生成的曲線圖更 有利于直觀現(xiàn)時地分析開采對地表的影響。南非的 Salamon應(yīng)用彈性理論提出了面元原理 ,將連續(xù)介質(zhì)力學(xué)與影響函數(shù)法相結(jié)合 ,為現(xiàn)在安徽理工大學(xué)畢業(yè)論文 3 的邊界單元法奠定了基礎(chǔ)。劉天泉提出和發(fā)展了礦山巖體采動影響與控制工程學(xué)。王兆亮等人利用 GIS 技術(shù)設(shè)計和開發(fā)了土地管理信息系統(tǒng)，系統(tǒng)利用 Visual ArcGIS Server為開發(fā)平臺，綜合運用 Web GIS技術(shù)、 Inter技 術(shù) 和 數(shù) 據(jù) 庫 技 術(shù) （ ArcSDE+SQLServer ） , 并 利 用 和ArcgisserverADFfor .NET 進行界面設(shè)計實現(xiàn)該信息管理系統(tǒng)。cosα12 （ 21） y’=y1+S’ 實現(xiàn)步驟 基于上文所述計算原理，在 AutoCAD 20xx 下使用 VBA 實現(xiàn)數(shù)據(jù)計算及成圖。煤層走向推進距離 ,開切眼坐標(biāo)，基準(zhǔn)點坐標(biāo) ,開切眼至基準(zhǔn)點距離 Dim plineObj As AcadLWPolyline Dim points(0 To 19) As Double HZ = Val() CH = Val() ZD = Val() KX = Val() KY = Val() JX1 = Val() JY1 = Val() 39。 dn amp。 AQ(3 * k) amp。 vbCrLf 安徽理工大學(xué)畢業(yè)論文 30 Next For k = 0 To UBound(brr) 1 AQ(3 * k) = ds(k): AQ(3 * k + 1) = we(k) / n + 20xx: AQ(3 * k + 2) = 0 mystr = mystr amp。 vbCrLf Next = mystr Dim splineObj As AcadSpline Dim startTan(0 To 2) As Double Dim endTan(0 To 2) As Double startTan(0) = 0: startTan(1) = 0: startTan(2) = 0 endTan(0) = 1500: endTan(1) = 0: endTan(2) = 0 Set splineObj = (AQ, startTan, endTan) End Sub ： Private Sub CommandButton9_Click() Dim layer6 As AcadLayer Set layer6 = (水平變形曲線 ) = 20 = layer6 Call bx Dim ws() As Double Dim AQ() As Double Dim ds() As Double 安徽理工大學(xué)畢業(yè)論文 28 Dim hs() As Double Dim gx() As Double Dim gy() As Double Dim hx() As Double Dim hy() As Double Dim n As Double Dim k As Integer Dim j As Integer Dim we() As Double n = Val() JX1 = Val() JY1 = Val() JX2 = Val() JY2 = Val() a = Atn((JY2 JY1) / (JX2 JX1)) Dim ss As String ss = Set fso = CreateObject() Set ttfile = (ss, 1) read = read = Replace(read, , ) Do While InStr(1, read, String(2, vbTab)) 0 read = Replace(read, String(2, vbTab), vbTab) Loop 39。MsgBox read brr = Split(read, vbCrLf) ReDim wn(UBound(brr)) ReDim AQ(3 * UBound(brr) 7) ReDim ds(UBound(brr)) ReDim gx(UBound(brr)) ReDim gy(UBound(brr)) ReDim wi(UBound(brr) 1) ReDim wk(UBound(brr) 2) For i = LBound(brr) To UBound(brr) If brr(i) Then arr = Split(brr(i), vbTab) 安徽理工大學(xué)畢業(yè)論文 24 If UBound(arr) 5 Then dn = arr(0) wn(i) = Format(CDbl(arr(6)) CDbl(arr(3)), ) gx(i) = Format(Sqr((CDbl(arr(1)) JX1) ^ 2 + (CDbl(arr(2)) JY1) ^ 2) * Cos(a) + JX1, ) gy(i) = Format(Sqr((CDbl(arr(1)) JX1) ^ 2 + (CDbl(arr(2)) JY1) ^ 2) * Sin(a) + JY1, ) ds(i) = Format(Sqr((gx(i) JX1) ^ 2 + (gy(i) JY1) ^ 2), ) End If End If Next For j = LBound(brr) To UBound(brr) 2 wi(j) = (wn(j + 1) wn(j)) / (ds(j + 1) ds(j)) 39。沿走向 Q = Sqr((QBX JX1) ^ 2 + (QBY JY1) ^ 2) Dim solidObj As AcadSolid Dim point1(0 To 2) As Double Dim point2(0 To 2) As Double Dim point3(0 To 2) As Double Dim point4(0 To 2) As Double 39。 從界面窗口中可以看到移動和變形的計算信息，同時 各 計算結(jié)果會自動保存在 txt 文檔里。則 n 號點至 n+1 號點 間的水平變形為 n 1 ~ n 1 ~ 01~n 1 ~ 0( ) ( )()n m nnnnlll? ?????? （ 27） 式中， （ l n+1 ~ n） m，（ l n+1 ~ n） 0——分別為 n+1 號點至 n 號點在首次和 m 次觀測時的水平距離。 α 12為觀測線方位角， S 為 N 到基準(zhǔn)點 JZ1 的距離。 曹化平等研究者采用基于概率積分法開采沉陷預(yù)計結(jié)果 , 利用 ArcGIS的 ArcScene模塊實現(xiàn)地表變形預(yù)計數(shù)據(jù)的三維可視化表達(dá) , 并以預(yù)計的下沉量為例進行可視化結(jié)果的等值線繪制、任意剖面圖制作、任意點位形變信息提取應(yīng)用分析 [14]。周國銓、虞 萬波、崔繼憲等提出了負(fù)指數(shù)函數(shù)法計算地表移動。 1931 年德國高等學(xué)校的礦業(yè)學(xué)院開始講授 “ 開采沉陷學(xué) ” 這門課程。 采用 AutoCAD VBA 二次開發(fā)技術(shù)，實現(xiàn)了地表移動觀測數(shù)據(jù) 處理、成圖自動一體化 。 傳統(tǒng)的人為管理模式會導(dǎo)致開采沉陷數(shù)據(jù)管理上的混亂，人力和物力的過多浪費，開采沉陷信息管理的費用也會增加，從而造成該企業(yè)單位的負(fù)擔(dān)過重，影響整個工程的進展速度。 secondary development。 在礦區(qū)，為保護地面建筑物、水體及鐵路，免受地下開采 活動的影響，減少地下資源損失，就必須研究地下開采所引起的巖層與地表移動規(guī)律。 1907 年 , Korten 發(fā)表了他自己的觀測成果 ,提出了水平移動和水平變形的分布規(guī)律 [4]。如國際著名的有限元分析軟件 ADINA、 ANSYS 以及專門用于巖土工程問題數(shù)值分析的 FLAC、 UDEC 等 [10]。 隨著計算機技術(shù)的發(fā)展和生活的需要許多和能源相關(guān)的高校和研究機構(gòu)都開始對開采沉陷監(jiān)測數(shù)據(jù)的計算機處理方面做過很多的研究，特別是中國礦業(yè)大學(xué)、山東科技大學(xué)、安徽理工大學(xué)等工礦類院校的研究已經(jīng)取得一些重要理論和實際效果。 圖 1 中 JZ1 和 JZ2 是為確定觀測線方向，在觀測線兩端選取的兩基準(zhǔn)點。 曲率 觀測線下沉的彎曲程度叫曲率，用 K表示，單位 mm/m2，則 n 號點附近的曲率，即n 1 點至 n+1 號點之間的曲率為 1 ~ ~ 1 1 ~ ~ 11 ~ ~ 1n 1 ~ n ~ 1 n 1 ~ n ~ 12 ( )( ) / 2n n n n n n n nn n nn n n ni i i iK l l l l? ? ? ???? ? ? ??????? （ 25） 式中， in+1~n， in~n1——分別為 n+1 號點至 n 號點和 n 號點至 n1 號點的傾斜； ln+1~n， ln~n1——分別為 n+1 號點 至 n 號點和 n 號點至 n1 號點的水平距離。 strtext amp。沿走向 安徽理工大學(xué)畢業(yè)論文 16 Q = Sqr((QBX JX1) ^ 2 + (QBY JY1) ^ 2) points(0) = 0: points(1) = HZ points(2) = 1000: points(3) = HZ + 1000 * Tan(j * / 180) points(4) = 1000: points(5) = HZ CH / Cos(j * / 180) + 1000 * Tan(j * / 180) points(6) = 0: points(7) = HZ CH / Cos(j * / 180) points(8) = 0: points(9) = HZ points(10) = Q: points(11) = HZ + Q * Tan(j * / 180) points(12) = Q + D * Cos(j * / 180): points(13) = HZ + Q * Tan(j * / 180) + D * Sin(j * / 180) points(14) = Q + D * Cos(j * / 180): points(15) = HZ + Q * Tan(j * / 180) + D * Sin(j * / 180) CH / Cos(j * / 180) points(16) = Q: points(17) = HZ + Q * Tan(j * / 180) CH / Cos(j * / 180) points(18) = Q: points(19) = HZ + Q * Tan(j * / 180) Set plineObj = (points) ZoomAll End Sub 煤層傾向采空區(qū)填充 Public Sub QXTC() Dim HZ, D, QBX, QBY, JX1, JY1, Q, CH, j 39。 , amp。 AQ(3 * k + 1) amp。 AQ(3 * k + 1) amp。 , amp。 , amp。煤層走向推進距離 ,開切眼坐標(biāo)，基準(zhǔn)點坐標(biāo) ,開切眼至基準(zhǔn)點距離 HZ = Val() CH = Val() ZD = Val() KX = Val() KY = Val() JX1 = Val() JY1 = Val() 39。煤層信息有平均采深、采厚、傾角；采動信息有沿走向線上開切眼坐標(biāo)、推進距離，沿傾向線上一邊界坐標(biāo)、傾向?qū)挾?。然后由 n、 n’ 的坐標(biāo)計算該點的移動變形量。其最大特點在于兩個方面：第一，與 VB 有著幾乎相同的開發(fā)環(huán)境和語法，具備功能強大和易于掌握的特點；第二，在于它的 for Application功能，針對性強，它駐留在主程序的 內(nèi)部，使其結(jié)構(gòu)精簡，且代碼運行效率高，使語言的引擎在技術(shù)上與開發(fā)環(huán)境分離。張玉卓提出了巖層移動的位錯理論。英國的 Datamine 用于礦山開采設(shè)計及數(shù)據(jù)處理 。 研究現(xiàn)狀 國外研究現(xiàn)狀 國外早開始了開采沉陷方面的研究工作，特別是比利時、前蘇聯(lián)、波蘭、 德國、澳大利亞、加拿大、美國等采礦業(yè)比較發(fā)達(dá)的國家對開采沉陷理論和實踐都進行了較深入研究。此過程中，上