【正文】 m wn(UBound(brr)) ReDim AQ(3 * UBound(brr) 1) ReDim ds(UBound(brr)) ReDim gx(UBound(brr)) ReDim gy(UBound(brr)) 安徽理工大學(xué)畢業(yè)論文 19 For i = LBound(brr) To UBound(brr) If brr(i) Then arr = Split(brr(i), vbTab) 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), ) mystr = mystr amp。煤層中間點(diǎn)采深，平均采厚 39。煤層中間點(diǎn)采深，平均采厚 39。首次和 第 m次觀測 的 平面坐標(biāo) 規(guī)劃實現(xiàn)代碼如下： 39。 AutoCAD 自 版本開始就支持 VBA 作為二次開發(fā)工具。 安徽理工大學(xué)畢業(yè)論文 7 下沉 各點(diǎn)的下沉以本次與首次的標(biāo)高差表示，單位 mm，則 m次觀測 n 點(diǎn)的下沉 Wn 為 Wn = Hn 0 — Hn m （ 23） 式中， Hn 0， Hn m—— 分別為首次和 m次觀測時 n點(diǎn)的高程。則 n 的坐標(biāo)為 x’=x1+S’ 技術(shù)路線圖 圖 技術(shù)路線圖 收集資料，學(xué)習(xí)理論知識 收集整理數(shù)據(jù) 檢驗數(shù)據(jù)質(zhì)量 輸出相關(guān)曲線 數(shù)據(jù)處理求取移動變形值 求取預(yù)計參數(shù) 地表移動特征分析 求取角量參數(shù) 求取動態(tài)參數(shù) 結(jié)論與建議 下沉 水平移動 傾斜 水平變 形 曲率 邊界角 移動角 充分采動角 超前影響距 超前影響角 下沉速度曲線 下沉系數(shù) 主要影響角正切 拐點(diǎn)偏距 水平移動系數(shù) 安徽理工大學(xué)畢業(yè)論文 6 2 觀測數(shù)據(jù)處理原理 坐標(biāo)規(guī)劃 目前，已廣泛應(yīng)用全站儀、 GPS 對地表移動觀測站進(jìn)行坐標(biāo)觀測，而不再是相鄰兩點(diǎn)間的水平距離和各測點(diǎn)偏離觀測線方向的支距 [1]。馬天勤利用 Surfer 軟件來實現(xiàn)開采沉陷可視化的效果，在 Surfer 軟件平臺基礎(chǔ)上，實現(xiàn)地表沉陷觀測和預(yù)計數(shù)據(jù)的可視化管理，實現(xiàn)等值線圖繪制、立體圖繪制和簡單空間分析（圖解分析、量算分析、數(shù)理統(tǒng)計分析等），提高了用戶對信息的管理效率，可更直觀 、有效地對沉陷數(shù)據(jù)進(jìn)行分析應(yīng)用 [16]。于廣明從非線性科學(xué)角度認(rèn)識開采沉陷的復(fù)雜性。謝和平應(yīng)用有限元分析巖層移動規(guī)律。 、 和 等采用計算程序計算了長壁式工作面采煤沉陷的地表應(yīng)力模式 [8]。 另外 , 60 年代初 ,英國學(xué)者 Berry 和 Sales 將巖體視為均質(zhì)彈性體 ,分為平面各向同性、橫觀各向同性與空間問題 3 類進(jìn)行分析 ,提出了計算巖體下沉的方法 [6]。 1903 年 , Halbaum 將采空區(qū)上方巖層作為懸臂梁 ,提出了地表應(yīng)變與曲率半徑成反比的理論。且圖中按原比例繪出對應(yīng)地下煤層開采情況，可很方便地從圖中直接量取各角量參數(shù)，距離參數(shù)。鑒于煤礦從未有過地表移動觀測資料，使得礦井在合理留設(shè)保護(hù)煤柱和開展“三下”采煤設(shè)計方面缺乏科學(xué)依據(jù)。 1 研究背景概述 研究意義 在礦區(qū)，為保護(hù)地面建筑物、水體及鐵路，免受地下開采活動的影響，減少地下資源損失，就必須研究地下開采所引起的巖層與地表移動規(guī)律。 生成的曲線圖更 有利于直觀現(xiàn)時地分析開采對地表的影響。且圖中按原比例繪出對應(yīng)地下煤層開采情況，可很方便地從圖中直接量取各角量參數(shù)，距離參數(shù)。如何有效地對開采沉陷信息管理和利用起來，如何使得開采沉陷的信息能夠在更大范圍內(nèi)安全地進(jìn)行傳輸、共享和管理，一直是很多企業(yè)單位急需解決的問題。 建立地表移動觀測站是研究開采沉陷和建筑物破壞規(guī)律的最可靠手段 。在輸出下沉、傾斜、曲率、水平移動、水平變形計算結(jié)果的 同時，將其各曲線以一定比例在 CAD 中成圖，可從曲線上量取任意點(diǎn)的對應(yīng)值。 1885年法國礦長 總結(jié)了開采沉陷時間的矛盾意見 [3]。學(xué)者 改進(jìn)了沉陷觀測的方法，并對開采沉陷的的理論方法也做出了巨大貢獻(xiàn) [5]。 隨著計算機(jī)和 3S 技術(shù)的廣泛應(yīng)用，將新技術(shù)應(yīng)用于開采沉陷觀測和數(shù)據(jù)處理也成為國際上普遍研究的的一個重要方面。王金莊等建立了用于預(yù)計主斷面地表下沉的典型曲線法。吳立新、王金莊提出了條帶開采覆巖破壞的托板理論。 郝長勝等研究 者利用 VB 作為系統(tǒng)開發(fā)前臺，基于 GIS 桌面軟件 MapInfo及 MapX 組件來進(jìn)行系統(tǒng)二次開發(fā)設(shè)計：根據(jù)模塊各功能應(yīng)盡量共享、數(shù)據(jù)相互調(diào)用、減少資源內(nèi)存儲量等原則，該系統(tǒng)煤礦開采沉陷預(yù)計分析整體流程分為數(shù)據(jù)管理、變形預(yù)計參數(shù)計算、圖形輸出、三維顯示及預(yù)計結(jié)果分析五大模塊。需要求取的地表移動參數(shù)，可在圖中直接精確量取和標(biāo)注 擬采用的研究思路與方法 研究思路與方法 以礦區(qū)實測數(shù)據(jù)為例 ,應(yīng)用 VBA 與 AutoCAD 二次開發(fā)技術(shù)進(jìn)行礦區(qū)地表移動變形值安徽理工大學(xué)畢業(yè)論文 5 的解 算 , 分析繪制地表移動變形曲線圖的 理論和方法 , 并以 VBA 為基礎(chǔ) , 實現(xiàn)了礦區(qū)地表移動變形值的快速求解 , 從而達(dá)到了對地下開采造成的地表移動變形的直觀性描述。方位角 α12 、 α ， S 的求取在此就不給出了。移動和變形的計算主要包括：各測點(diǎn)的下沉和水平移動，相鄰兩測點(diǎn)間的傾斜和水平變形，相鄰三點(diǎn)的曲率變形 [24]。其最大特點(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]。讀入基點(diǎn)坐標(biāo)，計算規(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次觀測在后。 煤層繪制的部分代碼如下： 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。 在 模型空間中創(chuàng)建實體對象 Set solidObj = (point1, point2, point3, point4) ZoomAll End Sub 煤層傾向代碼如下： Public Sub mcq() Dim HZ, D, QBX, QBY, JX1, JY1, Q, CH, j 39。 在模型空間中創(chuàng)建實體對象 Set solidObj = (point1, point2, point3, point4) ZoomAll End Sub 3）曲線繪制代碼 ： Private Sub CommandButton5_Click() Dim layer1 As AcadLayer Set layer1 = (下沉曲線 ) = acRed = layer1 Call ZS 安徽理工大學(xué)畢業(yè)論文 18 Dim wn() As Double Dim AQ() As Double Dim ds() As Double Dim gx() As Double Dim gy() As Double Dim n 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。 AQ(3 * k) amp。 wi(j) amp。MsgBox read brr = Split(read, vbCrLf) ReDim ws(UBound(brr)) ReDim AQ(3 * UBound(brr) 1) ReDim ds(UBound(brr)) ReDim hs(UBound(brr)) ReDim gx(UBound(brr)) ReDim gy(UBound(brr)) ReDim hx(UBound(brr)) ReDim hy(UBound(brr)) ReDim ws(UBound(brr)) For i = LBound(brr) To UBound(brr) If brr(i) Then arr = Split(brr(i), vbTab) If UBound(arr) 5 Then dn = arr(0) 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, ) 安徽理工大學(xué)畢業(yè)論文 27 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) End If End If Next For k = 0 To UBound(brr) 1 AQ(3 * k) = ds(k): AQ(3 * k + 1) = ws(k) / n: AQ(3 * k + 2) = 0 mystr = mystr amp。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 之后，將 單位一致的 變形 量繪在同一坐標(biāo)下 ，