【文章內(nèi)容簡介】 ln~n+1——n 號點(diǎn)至 n+1 號點(diǎn)的水平距離； Wn， Wn+1——分別為 n 號點(diǎn)和 n+1 號點(diǎn)的下沉量。 曲率 觀測線下沉的彎曲程度叫曲率，用 K表示，單位 mm/m2，則 n 號點(diǎn)附近的曲率，即n 1 點(diǎn)至 n+1 號點(diǎn)之間的曲率為 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 號點(diǎn)至 n 號點(diǎn)和 n 號點(diǎn)至 n1 號點(diǎn)的傾斜； ln+1~n， ln~n1——分別為 n+1 號點(diǎn) 至 n 號點(diǎn)和 n 號點(diǎn)至 n1 號點(diǎn)的水平距離。 水平移動 下沉盆地內(nèi)，某點(diǎn)沿某一水平方向的位移叫水平移動，用 U 表示，單位 mm，則 m次觀測 n 點(diǎn)的水平移動為 Un= Ln m — Ln 0 （ 26） 式中， Ln m， Ln 0——分別為 m 次觀測時和首次觀測時 n 號點(diǎn)至觀測線控制點(diǎn)間的水平距離，直接用坐標(biāo)求得。 安徽理工大學(xué)畢業(yè)論文 8 水平變形 下沉盆地內(nèi)兩點(diǎn)間單位長度的水平移動差稱之為水平變形，以 ε 表示，單位為 mm/m。則 n 號點(diǎn)至 n+1 號點(diǎn) 間的水平變形為 n 1 ~ n 1 ~ 01~n 1 ~ 0( ) ( )()n m nnnnlll? ?????? （ 27） 式中， （ l n+1 ~ n） m，（ l n+1 ~ n） 0——分別為 n+1 號點(diǎn)至 n 號點(diǎn)在首次和 m 次觀測時的水平距離。 安徽理工大學(xué)畢業(yè)論文 9 3 軟件實(shí)現(xiàn) AutoCAD VBA 簡介 VBA ( visual basic for application) 是由 Microsoft 創(chuàng)建的，用來自動執(zhí)行任務(wù)的一個編程環(huán)境，它提供了一些用來創(chuàng)建圖 形用戶界面的可拖拉工具和用來與AutoCAD 對象交換的編程語言 [5]。其最大特點(diǎn)在于兩個方面：第一，與 VB 有著幾乎相同的開發(fā)環(huán)境和語法，具備功能強(qiáng)大和易于掌握的特點(diǎn)；第二，在于它的 for Application功能，針對性強(qiáng)，它駐留在主程序的內(nèi)部，使其結(jié)構(gòu)精簡，且代碼運(yùn)行效率高，使語言的引擎在技術(shù)上與開發(fā)環(huán)境分離 [6]。 AutoCAD 自 版本開始就支持 VBA 作為二次開發(fā)工具。 實(shí)現(xiàn)步驟 基于上文所述計(jì)算原理，在 AutoCAD 20xx 下使用 VBA 實(shí)現(xiàn)數(shù)據(jù)計(jì)算及成圖。 程序界面如圖 2 。 確定煤層采動信息 為反應(yīng)地下開采與地表移動和變形規(guī)律，在進(jìn)行各項(xiàng)計(jì)算之前應(yīng)當(dāng)首先確定觀測數(shù)據(jù)所對應(yīng)的地下開采情況。煤層信息有平均采深、采厚、傾角；采動信息有沿走向線上開切眼坐標(biāo)、推進(jìn)距離，沿傾向線上一邊界坐標(biāo)、傾向?qū)挾?。 安徽理工大學(xué)畢業(yè)論文 10 圖 2 軟件界面 確定觀測線基準(zhǔn)點(diǎn)坐標(biāo) 在觀測線兩端選取兩點(diǎn)，坐標(biāo)可在 cad 觀測線布置圖中直接獲取。 計(jì)算 觀測線 方向代碼如下： 39。讀入基點(diǎn)坐標(biāo)，計(jì)算規(guī)劃方向 JX1 = Val() JY1 = Val() JX2 = Val() JY2 = Val() a = Atn((JY2 JY1) / (JX2 JX1)) 安徽理工大學(xué)畢業(yè)論文 11 讀入 觀測數(shù)據(jù)進(jìn)行坐標(biāo)規(guī)劃 數(shù)據(jù)文件編輯成 txt 文檔，首次觀測平面坐標(biāo)及高程在前，第 m次觀測在后。首次和 第 m次觀測 的 平面坐標(biāo) 規(guī)劃實(shí)現(xiàn)代碼如下： 39。 數(shù)組 gX gY1 為 首次觀測 的規(guī)劃后 平面坐標(biāo) ，數(shù)組 hX hY2為 第 m次觀測 的規(guī)劃后 平面坐標(biāo) 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), ) hx(i) = Format(Sqr((CDbl(arr(4)) JX1) ^ 2 + (CDbl(arr(5)) JY1) ^ 2) * Cos(a) + JX1, ) hy(i) = Format(Sqr((CDbl(arr(4)) JX1) ^ 2 + (CDbl(arr(5)) JY1) ^ 2) * Sin(a) + JY1, ) hs(i) = Format(Sqr((hx(i) JX1) ^ 2 + (hy(i) JY1) ^ 2), ) ws(i) = hs(i) ds(i) 文件的打開，文件的導(dǎo)出與程序的退出 1）文件的打開代碼 Private Sub CommandButton1_Click() Dim strtext As String = 文本文件 (*.txt)|*.txt If Then = Open For Input As 1 Do While Not EOF(1) Line Input 1, strtext = amp。 strtext amp。 Chr(13) amp。 Chr(10) 安徽理工大學(xué)畢業(yè)論文 12 Loop Close 1 End If End Sub 2）文件的導(dǎo)出代碼： Private Sub CommandButton2_Click() = 文本文件 (*.txt)|*.txt If Then Open For Output As 1 Print 1, Close 1 End If End Sub 3)程序的退出代碼： Private Sub CommandButton4_Click() End End Sub 計(jì)算移動和變形 各點(diǎn)規(guī)劃至觀測線方向后，根據(jù)公 式 （ 23） ~（ 26） 計(jì)算出移動和變形。 從界面窗口中可以看到移動和變形的計(jì)算信息，同時 各 計(jì)算結(jié)果會自動保存在 txt 文檔里。 曲線圖的繪制 首先依原比例繪出煤層采動信息，及煤層與測點(diǎn)的空間關(guān)系。 煤層繪制的部分代碼如下： 1) 煤層走向信息繪制 Private Sub OptionButton1_Click() Dim layer1 As AcadLayer 安徽理工大學(xué)畢業(yè)論文 13 Set layer1 = (坐標(biāo) ) = 255 = layer1 Call MC Call CKQTC End Sub 2）煤層傾向信息繪制 Private Sub OptionButton2_Click() Dim layer1 As AcadLayer Set layer1 = (坐標(biāo) ) = 255 = layer1 Call mcq Call QXTC End Sub Public Sub MC() Dim HZ, ZD, KX, KY, JX1, JY1, SKQ, CH 39。煤層中間點(diǎn)采深，平均采厚 39。煤層走向推進(jìn)距離 ,開切眼坐標(biāo)，基準(zhǔn)點(diǎn)坐標(biāo) ,開切眼至基準(zhǔn)點(diǎ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。沿走向 安徽理工大學(xué)畢業(yè)論文 14 SKQ = Sqr((KX JX1) ^ 2 + (KY JY1) ^ 2) points(0) = 0: points(1) = HZ points(2) = 1000: points(3) = HZ points(4) = 1000: points(5) = HZ CH points(6) = 0: points(7) = HZ CH points(8) = 0: points(9) = HZ points(10) = SKQ: points(11) = HZ points(12) = SKQ + ZD: points(13) = HZ points(14) = SKQ + ZD: points(15) = HZ CH points(16) = SKQ: points(17) = HZ CH points(18) = SKQ: points(19) = HZ Set plineObj = (points) ZoomAll End Sub 煤層走向代碼如下： Public Sub CKQTC() Dim HZ, ZD, KX, KY, JX1, JY1, SKQ, CH 39。煤層中間點(diǎn)采深，平均采厚 39。煤層走向推進(jìn)距離 ,開切眼坐標(biāo)，基準(zhǔn)點(diǎn)坐標(biāo) ,開切眼至基準(zhǔn)點(diǎn)距離 HZ = Val() CH = Val() ZD = Val() KX = Val() KY = Val() JX1 = Val() JY1 = Val() 39。沿走向 SKQ = Sqr((KX JX1) ^ 2 + (KY JY1) ^ 2) Dim solidObj As AcadSolid 安徽理工大學(xué)畢業(yè)論文 15 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。 定義實(shí)體 point1(0) = SKQ: point1(1) = HZ CH: point1(2) = 0 point2(0) = SKQ + ZD: point2(1) = HZ CH: point2(2) = 0 point3(0) = SKQ: point3(1) = HZ: point3(2) = 0 point4(0) = SKQ + ZD: point4(1) = HZ: point4(2) = 0 39。 在 模型空間中創(chuàng)建實(shí)體對象 Set solidObj = (point1, point2, point3, point4) ZoomAll End Sub 煤層傾向代碼如下： Public Sub mcq() Dim HZ, D, QBX, QBY, JX1, JY1, Q, CH, j 39。煤層中間點(diǎn)采深，平均采厚 39。煤層走向推進(jìn)距離 ,開切眼坐標(biāo)，基準(zhǔn)點(diǎn)坐標(biāo) ,開切眼至基準(zhǔn)點(diǎn)距離 Dim plineObj As AcadLWPolyline Dim points(0 To 19) As Double HZ = Val() CH = Val() D = Val() j = Val() QBX = Val() QBY = Val() JX1 = Val() JY1 = Val() 39。沿走向 安徽理工大學(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 *