【正文】 下沉盆地內(nèi)兩點(diǎn)間單位長(zhǎng)度的水平移動(dòng)差稱(chēng)之為水平變形，以 ε 表示，單位為 mm/m。煤層信息有平均采深、采厚、傾角；采動(dòng)信息有沿走向線(xiàn)上開(kāi)切眼坐標(biāo)、推進(jìn)距離，沿傾向線(xiàn)上一邊界坐標(biāo)、傾向?qū)挾?。 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ì)算移動(dòng)和變形 各點(diǎn)規(guī)劃至觀測(cè)線(xiàn)方向后，根據(jù)公 式 （ 23） ~（ 26） 計(jì)算出移動(dòng)和變形。煤層走向推進(jìn)距離 ,開(kāi)切眼坐標(biāo)，基準(zhǔn)點(diǎn)坐標(biāo) ,開(kāi)切眼至基準(zhǔn)點(diǎn)距離 HZ = Val() CH = Val() ZD = Val() KX = Val() KY = Val() JX1 = Val() JY1 = Val() 39。煤層走向推進(jìn)距離 ,開(kāi)切眼坐標(biāo)，基準(zhǔn)點(diǎn)坐標(biāo) ,開(kāi)切眼至基準(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() 安徽理工大學(xué)畢業(yè)論文 17 JX1 = Val() JY1 = Val() 39。 , 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 CommandButton7_Click() Dim layer4 As AcadLayer Set layer4 = (曲率曲線(xiàn) ) = acGreen = layer4 Call qlqx Dim wn() As Double Dim wk() As Double Dim wi() As Double Dim AQ() As Double Dim ds() As Double Dim gx() As Double Dim gy() As Double Dim n As Double Dim j As Integer 安徽理工大學(xué)畢業(yè)論文 23 Dim m As Integer 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。 , amp。 AQ(3 * k + 2) amp。 AQ(3 * k + 1) amp。 AQ(3 * k) amp。 AQ(3 * k + 1) amp。 , amp。 , amp。 , 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。沿走向 安徽理工大學(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。 strtext amp。 程序界面如圖 2 。 曲率 觀測(cè)線(xiàn)下沉的彎曲程度叫曲率，用 K表示，單位 mm/m2，則 n 號(hào)點(diǎn)附近的曲率，即n 1 點(diǎn)至 n+1 號(hào)點(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 號(hào)點(diǎn)至 n 號(hào)點(diǎn)和 n 號(hào)點(diǎn)至 n1 號(hào)點(diǎn)的傾斜； ln+1~n， ln~n1——分別為 n+1 號(hào)點(diǎn) 至 n 號(hào)點(diǎn)和 n 號(hào)點(diǎn)至 n1 號(hào)點(diǎn)的水平距離。sinα12 （ 22） 其中 ， S’=S 圖 1 中 JZ1 和 JZ2 是為確定觀測(cè)線(xiàn)方向，在觀測(cè)線(xiàn)兩端選取的兩基準(zhǔn)點(diǎn)。此外，還有很多學(xué)者對(duì)開(kāi)采沉陷數(shù)據(jù)處理、管理和預(yù)算等這一方面的研究。 隨著計(jì)算機(jī)技術(shù)的發(fā)展和生活的需要許多和能源相關(guān)的高校和研究機(jī)構(gòu)都開(kāi)始對(duì)開(kāi)采沉陷監(jiān)測(cè)數(shù)據(jù)的計(jì)算機(jī)處理方面做過(guò)很多的研究，特別是中國(guó)礦業(yè)大學(xué)、山東科技大學(xué)、安徽理工大學(xué)等工礦類(lèi)院校的研究已經(jīng)取得一些重要理論和實(shí)際效果。楊倫提出了巖層二次壓縮理論 ,將地表下沉直接與巖體的物理力學(xué)性質(zhì)聯(lián)系起來(lái)。如國(guó)際著名的有限元分析軟件 ADINA、 ANSYS 以及專(zhuān)門(mén)用于巖土工程問(wèn)題數(shù)值分析的 FLAC、 UDEC 等 [10]。 90 年代至今，隨著科學(xué)的發(fā)展和研究手段的進(jìn)步 ,又掀起了礦山開(kāi)采沉陷的再一次研究 熱潮。 1907 年 , Korten 發(fā)表了他自己的觀測(cè)成果 ,提出了水平移動(dòng)和水平變形的分布規(guī)律 [4]。 文中 AutoCAD 進(jìn)行二次開(kāi)發(fā)，應(yīng)用 VBA 計(jì)算出各觀測(cè)點(diǎn)的各期移動(dòng)變形值，各種變形曲線(xiàn)依不同比例在 CAD 中成圖，且對(duì)應(yīng)繪出地下煤層采動(dòng)情況。 在礦區(qū)，為保護(hù)地面建筑物、水體及鐵路，免受地下開(kāi)采 活動(dòng)的影響，減少地下資源損失，就必須研究地下開(kāi)采所引起的巖層與地表移動(dòng)規(guī)律。在確保野外觀測(cè)數(shù)據(jù)的正確性后，計(jì)算出 各觀測(cè)點(diǎn)的移動(dòng)變形信息，最終將計(jì)算結(jié)果繪制成曲線(xiàn)圖，以描述沿觀測(cè)線(xiàn)的地表移動(dòng)和變形的分布特征及其發(fā)展過(guò)程，以便對(duì)開(kāi)采損害做出評(píng)價(jià)。 secondary development?？朔藗鹘y(tǒng)方法成圖只能 查 看、不易定量分析的缺點(diǎn)。 傳統(tǒng)的人為管理模式會(huì)導(dǎo)致開(kāi)采沉陷數(shù)據(jù)管理上的混亂，人力和物力的過(guò)多浪費(fèi)，開(kāi)采沉陷信息管理的費(fèi)用也會(huì)增加，從而造成該企業(yè)單位的負(fù)擔(dān)過(guò)重，影響整個(gè)工程的進(jìn)展速度。造成地表建筑物損害和塌陷區(qū)生態(tài)環(huán)境破壞，給周邊人民的生產(chǎn)、生活帶來(lái)了不同程度的影響。 采用 AutoCAD VBA 二次開(kāi)發(fā)技術(shù)，實(shí)現(xiàn)了地表移動(dòng)觀測(cè)數(shù)據(jù) 處理、成圖自動(dòng)一體化 。 據(jù)參考文獻(xiàn)，早在 1825 年，比利時(shí)的工程師們?cè)谡{(diào)查 Liege 市地表裂縫的原因時(shí)，證實(shí)了礦山巷道和地表之間 300 英尺的距離足以阻止地表沉陷，開(kāi)始了最早的開(kāi)采沉陷工作的研究 [2]。 1931 年德國(guó)高等學(xué)校的礦業(yè)學(xué)院開(kāi)始講授 “ 開(kāi)采沉陷學(xué) ” 這門(mén)課程。澳大利亞的 Mine TER 以及 Mine Map,用于礦山設(shè)計(jì)和繪制工程圖 ,三維開(kāi)采模擬。周?chē)?guó)銓、虞 萬(wàn)波、崔繼憲等提出了負(fù)指數(shù)函數(shù)法計(jì)算地表移動(dòng)。楊碩建立采動(dòng)損害空間的力學(xué)模式。 曹化平等研究者采用基于概率積分法開(kāi)采沉陷預(yù)計(jì)結(jié)果 , 利用 ArcGIS的 ArcScene模塊實(shí)現(xiàn)地表變形預(yù)計(jì)數(shù)據(jù)的三維可視化表達(dá) , 并以預(yù)計(jì)的下沉量為例進(jìn)行可視化結(jié)果的等值線(xiàn)繪制、任意剖面圖制作、任意點(diǎn)位形變信息提取應(yīng)用分析 [14]。 AutoCAD 自 版本開(kāi)始就支持 VBA 作為二次開(kāi)發(fā)工具。 α 12為觀測(cè)線(xiàn)方位角， S 為 N 到基準(zhǔn)點(diǎn) JZ1 的距離。該方法 使用坐標(biāo)觀測(cè)量進(jìn)行計(jì)算， 避免了復(fù)雜的支距改正，提高了計(jì)算精度。則 n 號(hào)點(diǎn)至 n+1 號(hào)點(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 號(hào)點(diǎn)至 n 號(hào)點(diǎn)在首次和 m 次觀測(cè)時(shí)的水平距離。 安徽理工大學(xué)畢業(yè)論文 10 圖 2 軟件界面 確定觀測(cè)線(xiàn)基準(zhǔn)點(diǎn)坐標(biāo) 在觀測(cè)線(xiàn)兩端選取兩點(diǎn)，坐標(biāo)可在 cad 觀測(cè)線(xiàn)布置圖中直接獲取。 從界面窗口中可以看到移動(dòng)和變形的計(jì)算信息，同時(shí) 各 計(jì)算結(jié)果會(huì)自動(dòng)保存在 txt 文檔里。沿走向 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。沿走向 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。 wn(i) amp。MsgBox read brr = Split(read, vbCrLf) ReDim wn(UBound(brr))