基于cad和vba地表移动观测站数据处理系统.doc

安徽理工大学毕业论文

本科毕业论文

基于CAD和VBA的地表移动观测站数据处理系统
DATA PROCESSING SYSTEM ON MINING AREA SURFACE MOVEMENT STATION BASED ON AUTOCAD VBA

学院（部）： 测绘学院
专业班级： 测绘09-1班
学生姓名： 刘石磊
指导教师： 王磊讲师

2013 年 5 月 31 日

55

基于CAD和VBA的地表移动观测站数据处理系统

摘 要

采用AutoCAD VBA二次开发技术，实现了地表移动观测数据处理、成图自动一体化。克服了传统方法成图只能查看、不易定量分析的缺点。在输出下沉、倾斜、曲率、水平移动、水平变形计算结果的同时，将其各曲线以一定比例在CAD中成图，可从曲线上量取任意点的对应值。且图中按原比例绘出对应地下煤层开采情况，可很方便地从图中直接量取各角量参数，距离参数。生成的曲线图更有利于直观现时地分析开采对地表的影响。

关键词：AutoCAD VBA；二次开发；地表移动变形；数据处理

DATA PROCESSING SYSTEM ON MINING AREA SURFACE MOVEMENT STATION BASED ON AUTOCAD VBA

ABSTRACT

AutoCAD VBA secondary development technique is used to realize the mining area surface movement observation data fast processing and automatic mapping in this paper. This way can overcome traditional methods’ shortcomings which only can make the map for watching, but not for quantitative analysis. The mining surface’s subside, lean, curvature, horizontal moving and horizontal deformation is calculated. And the all results are put out. At the same time each curve is drawn in AutoCAD with certain proportion. Any point’s value can be taken from the AutoCAD form. By the way, the underground coal mining condition is painted in the map with original proportion. So angle parameters and distance parameters can be easily measured and marked. All of the above can make the visual description on the surface deformation of the underground mining convenient and immediate.

Key words: AutoCAD VBA; secondary development; surface movement and deformation; data processing

目录

摘 要 II
ABSTRACT III
绪论 1
1研究背景概述 1
1.1研究意义 1
1.2研究现状 2
1.2.1 国外研究现状 2
1.2.2 国内研究现状 3
1.3研究内容 4
1.4拟采用的研究思路与方法 4
1.4.1研究思路与方法 4
1.4.2技术路线图 5
2观测数据处理原理 6
2.1坐标规划 6
2.2移动和变形的计算 6
2.2.1 下沉 7
2.2.2倾斜 7
2.2.3曲率 7
2.2.4水平移动 7
2.2.5水平变形 8
3软件实现 9
3.1 AutoCAD VBA简介 9
3.2实现步骤 9
3.2.1 确定煤层采动信息 9
3.2.2 确定观测线基准点坐标 10
3.2.3读入观测数据进行坐标规划 11
3.2.4文件的打开，文件的导出与程序的退出 11
3.2.4 计算移动和变形 12
3.2.5 曲线图的绘制 12
4实例应用 31
4.1矿区概况 31
4.1.1地理交通 31
4.1.2生产情况 31
4.1.3地形地貌 32
4.2矿井地质采矿条件 32
4.2.1煤层赋存情况 32
4.2.2 东坪煤矿15412工作面地质采矿条件 33
4.3地表移动观测站设计 34
4.3.1观测站设计依据和原则 34
4.3.2观测线长度及位置的确定 35
4.3.3 观测站点位布置情况 36
4.4地表移动观测站的观测工作 37
4.4.1连接测量 37
4.4.2全面观测 38
4.4.3日常观测工作 40
4.5观测站数据处理 41
4.6地表移动参数的求取 45
4.6.1角量参数的求取 45
4.6.2 动态参数的求取 46
结论 51
参考文献： 52
致 谢 53

绪论

开采沉陷的研究对国民经济的发展和人民生活水平的提高有着重要的意义，然而随着开采沉陷研究的规模越大，有关开采沉陷的信息量也相应的增加，对这些信息的存储和管理也就相当繁琐。 传统的人为管理模式会导致开采沉陷数据管理上的混乱，人力和物力的过多浪费，开采沉陷信息管理的费用也会增加，从而造成该企业单位的负担过重，影响整个工程的进展速度。目前很多的开采沉陷信息数据都是用Word、Excel或者其他软件进行数据管理和处理并存储在计算机中，但再对这些信息进行查看和处理时很不方便和直观。如何有效地对开采沉陷信息管理和利用起来，如何使得开采沉陷的信息能够在更大范围内安全地进行传输、共享和管理，一直是很多企业单位急需解决的问题。
1研究背景概述
1.1研究意义
在矿区，为保护地面建筑物、水体及铁路，免受地下开采活动的影响，减少地下资源损失，就必须研究地下开采所引起的岩层与地表移动规律。目前是最直接有效的方法是建立地表移动观测站，定期对观测站上的监测点及控制点进行观测。在确保野外观测数据的正确性后，计算出各观测点的移动变形信息，最终将计算结果绘制成曲线图，以描述沿观测线的地表移动和变形的分布特征及其发展过程，以便对开采损害做出评价。
地下煤层被开采出来以后，岩体内部形成一个采空区，导致周围岩体的原有应力平衡状态受到破坏，引起应力的重新分布，使岩体产生移动和破坏，直到达到新的平衡。此过程中，上覆岩层将产生移动、变形与破坏，当开采面积达到一定范围时，移动变形将波及地表，形成地表移动盆地[1]。造成地表建筑物损害和塌陷区生态环境破坏，给周边人民的生产、生活带来了不同程度的影响。为了有效预防和治理这类沉陷灾害，需要获取地表移动规律。
建立地表移动观测站是研究开采沉陷和建筑物破坏规律的最可靠手段。鉴于煤矿从未有过地表移动观测资料，使得矿井在合理留设保护煤柱和开展“三下”采煤设计方面缺乏科学依据。为此，根据《煤矿安全规程》和《建筑物、水体、铁路及主要井巷煤柱留设与压煤开采规程》中的有关规定，在煤矿工作面上方建立地表移动观测站，开展地表移动规律研究，以期获得本区的地表移动盆地角值参数、动态变形参数和概率积分法预计参数，为煤矿各类保护煤柱留设和开采沉陷预测研究提供科学依据。
在矿区，为保护地面建筑物、水体及铁路，免受地下开采活动的影响，减少地下资源损失，就必须研究地下开采所引起的岩层与地表移动规律。目前是最直接有效的方法是建立地表移动观测站，定期对观测站上的监测点及控制点进行观测。在确保野外观测数据的正确性后，计算出各观测点的移动变形信息，最终将计算结果绘制成曲线图，以描述沿观测线的地表移动和变形的分布特征及其发展过程，以便对开采损害做出评价。
采用AutoCAD VBA二次开发技术，实现了地表移动观测数据处理、成图自动一体化。克服了传统方法成图只能查看、不易定量分析的缺点。在输出下沉、倾斜、曲率、水平移动、水平变形计算结果的同时，将其各曲线以一定比例在CAD中成图，可从曲线上量取任意点的对应值。且图中按原比例绘出对应地下煤层开采情况，可很方便地从图中直接量取各角量参数，距离参数。生成的曲线图更有利于直观现时地分析开采对地表的影响。文中AutoCAD进行二次开发，应用VBA计算出各观测点的各期移动变形值，各种变形曲线依不同比例在CAD中成图，且对应绘出地下煤层采动情况。需要求取的地表移动参数，可在图中直接精确量取和标注。
1.2研究现状
1.2.1 国外研究现状
国外早开始了开采沉陷方面的研究工作，特别是比利时、前苏联、波兰、德国、澳大利亚、加拿大、美国等采矿业比较发达的国家对开采沉陷理论和实践都进行了较深入研究。 据参考文献，早在1825年，比利时的工程师们在调查Liege市地表裂缝的原因时，证实了矿山巷道和地表之间300英尺的距离足以阻止地表沉陷，开始了最早的开采沉陷工作的研究[2]。在1838年,多里斯提出了开采沉陷的第一个理论“垂线理论”。1885年法国矿长H.Fay总结了开采沉陷时间的矛盾意见[3]。1903年, Halbaum将采空区上方岩层作为悬臂梁,提出了地表应变与曲率半径成反比的理论。1900年Wachsmann首次提出了开采沉陷的现代观点，即煤层开采时底部煤层最先崩塌，中部煤层下沉且破裂，上部煤层下沉但不破裂。1907年, Korten发表了他自己的观测成果,提出了水平移动和水平变形的分布规律[4]。本世纪20至30年代,斯奇米茨、凯因斯特和巴斯研究了开采影响的作用面积及分带,可作为影响函数法的初级阶段;凯因斯特又提出了水平移动的算式,即U =Wtg (地表点到开采中心连线与铅垂线的夹角)。1931年德国高等学校的矿业学院开始讲授“开采沉陷学”这门课程。英国自1930年开始了移动变形的观测并于1950年发现了观测和地表变形之间关系的重要性，建立了不同采动程度下的下沉系数修正体系。学者K.Wardell改进了沉陷观测的方法，并对开采沉陷的的理论方法也做出了巨大贡献[5]。 另外, 60年代初,英国学者Berry和Sales将岩体视为均质弹性体,分为平面各向同性、横观各向同性与空间问题3类进行分析,提出了计算岩体下沉的方法[6]。南非的Salamon应用弹性理论提出了面元原理,将连续介质力学与影响函数法相结合,为现在的边界单元法奠定了基础。90年代至今，随着科学的发展和研究手段的进步,又掀起了矿山开采沉陷的再一次研究热潮。在西方发达国家,利用信息技术开发出了一批适用于矿山的应用软件如加拿大的LYNX系统,用于地质学研究从勘探资料分析地质特征;英国的Datamine用于矿山开采设计及数据处理;澳大利亚的Mine TER以及Mine Map,用于矿山设计和绘制工程图,三维开采模拟。国外很多国家用基于概率基本法[7]预计计算来指导建筑物下、水体下和铁路下采煤。 随着计算机和3S技术的广泛应用，将新技术应用于开采沉陷观测和数据处理也成为国际上普遍研究的的一个重要方面。B.N.Whittaker、D.J.Reddish和D.Fitzpatrick等采用计算程序计算了长壁式工作面采煤沉陷的地表应力模式[8]。德国将数字摄影测量技术和遥感技术手段应用到鲁尔矿区的开采沉陷数据采集和处理[9]。如国际著名的有限元分析软件ADINA、ANSYS以及专门用于岩土工程问题数值分析的FLAC、UDEC等[10]。
1.2.2 国内研究现状
我国是以煤炭为能源的国家，自新中国成立以来就开始开采沉陷方面的研究，80、90年代至今，我国开采沉陷理论和实践研究出现了日新月异的发展。何国淸、马伟民、王金庄建立了碎块体理论—地表沉陷的威布尔分布。周国铨、虞万波、崔继宪等提出了负指数函数法计算地表移动。邹友峰对地表下沉预计方法进行研究。王金庄等建立了用于预计主断面地表下沉的典型曲线法。谢和平应用有限元分析岩层移动规律。刘天泉提出和发展了矿山岩体采动影响与控制工程学。杨伦提出了岩层二次压缩理论,将地表下沉直接与岩体的物理力学性质联系起来。李增琪建立了岩层移动的三维层状模型。张玉卓提出了岩层移动的位错理论。杨硕建立采动损害空间的力学模式。邓喀中提出了岩体开采沉陷的结构效应。吴立新、王金庄提出了条带开采覆岩破坏的托板理论。于广明从非线性科学角度认识开采沉陷的复杂性。徐乃忠、郭惟菇研究了采动覆岩离层形成的基本规律和离层注浆控制地表下沉的理论与实践[11]。 随着计算机技术的发展和生活的需要许多和能源相关的高校和研究机构都开始对开采沉陷监测数据的计算机处理方面做过很多的研究，特别是中国矿业大学、山东科技大学、安徽理工大学等工矿类院校的研究已经取得一些重要理论和实际效果。 陈宜金，黄绍东对开采沉陷信息处理软件系统各模块间关系和设计原则进行了探讨[12]。王卷乐，吴侃把空间数据可视化思想引入开采沉陷的防治中，对数据的可视化的视觉表达和分析进行了分析[13]。 曹化平等研究者采用基于概率积分法开采沉陷预计结果, 利用ArcGIS的ArcScene模块实现地表变形预计数据的三维可视化表达, 并以预计的下沉量为例进行可视化结果的等值线绘制、任意剖面图制作、任意点位形变信息提取应用分析[14]。郑艳等人提出智能信息处理系统在开采沉陷预测中的应用，通过建立基于遗传算法的改进BP神经网络开采沉陷预测模型，利用矿区大量的地表沉陷实际观测数据样本对该神经网络进行训练和学习，并用该网络对几组数据进行采动地表沉陷预测，改进的BP神经网路具有收敛速度快、预测参数精度高的优点等[15]。 郝长胜等研究者利用VB作为系统开发前台，基于GIS桌面软件MapInfo及MapX组件来进行系统二次开发设计：根据模块各功能应尽量共享、数据相互调用、减少资源内存储量等原则，该系统煤矿开采沉陷预计分析整体流程分为数据管理、变形预计参数计算、图形输出、三维显示及预计结果分析五大模块。马天勤利用Surfer软件来实现开采沉陷可视化的效果，在Surfer软件平台基础上，实现地表沉陷观测和预计数据的可视化管理，实现等值线图绘制、立体图绘制和简单空间分析（图解分析、量算分析、数理统计分析等），提高了用户对信息的管理效率，可更直观、有效地对沉陷数据进行分析应用[16]。王兆亮等人利用GIS技术设计和开发了土地管理信息系统，系统利用Visual Studio.NET和ArcGIS Server为开发平台，综合运用Web GIS技术、Internet技术和数据库技术（ArcSDE+SQLServer）,并利用ASP.NET2.0和ArcgisserverADFfor .NET进行界面设计实现该信息管理系统。此外，还有很多学者对开采沉陷数据处理、管理和预算等这一方面的研究。例如贾小敏、余学祥等讨论了应用Visual Basic和Auto CAD二次开发技术进行矿区地表移动变形值计算和变形曲线绘制的方法，实现了在此环境下变形值的计算和曲线的
1.3研究内容
VBA ( visual basic for application) 是由Microsoft创建的，用来自动执行任务的一个编程环境，它提供了一些用来创建图形用户界面的可拖拉工具和用来与AutoCAD对象交换的编程语言。其最大特点在于两个方面：第一，与VB有着几乎相同的开发环境和语法，具备功能强大和易于掌握的特点；第二，在于它的for Application功能，针对性强，它驻留在主程序的内部，使其结构精简，且代码运行效率高，使语言的引擎在技术上与开发环境分离。AutoCAD自R14.01版本开始就支持VBA作为二次开发工具。
对AutoCAD进行二次开发，应用VBA计算出各观测点的各期移动变形值，各种变形曲线依不同比例在CAD中成图，且对应绘出地下煤层采动情况。需要求取的地表移动参数，可在图中直接精确量取和标注
1.4拟采用的研究思路与方法
1.4.1研究思路与方法
以矿区实测数据为例,应用VBA 与AutoCAD 二次开发技术进行矿区地表移动变形值的解
算, 分析绘制地表移动变形曲线图的理论和方法, 并以VBA 为基础, 实现了矿区地表移动变形值的快速求解, 从而达到了对地下开采造成的地表移动变形的直观性描述。
1.4.2技术路线图
收集资料，学习理论知识
收集整理数据
检验数据质量
输出相关曲线
数据处理求取移动变形值
求取预计参数
地表移动特征分析
求取角量参数
求取动态参数
结论与建议
下沉
水平移动
倾斜
水平变形
曲率
边界角
移动角
充分采动角
超前影响距
超前影响角
下沉速度曲线
下沉系数
主要影响角正切
拐点偏距
水平移动系数

图1.1 技术路线图

2观测数据处理原理
2.1坐标规划
目前，已广泛应用全站仪、GPS对地表移动观测站进行坐标观测，而不再是相邻两点间的水平距离和各测点偏离观测线方向的支距[1]。由于各点或多或少要偏离观测线，而要计算的则是各点沿观测线方向的移动和变形，因此本文中采用首先将各点规划到观测线方向上，之后再进行变形量的计算。
图1中JZ1和JZ2是为确定观测线方向，在观测线两端选取的两基准点。N为观测线上的一观测点，N’为N受地下开采影响后的点。n、n’为N、 N’在观测线方向上的投影，即规划到观测线上的点位。α12为观测线方位角，S为N 到基准点JZ1的距离。

图1 坐标规划示意图
设点JZ1、JZ2的坐标为（x1，y1），（x2，y2），N、n坐标为（x，y）、（x’，y’）。方位角α12、α，S的求取在此就不给出了。则n的坐标为
x’=x1+S’·cosα12 （2-1）
y’=y1+S’·sinα12 （2-2）
其中，
S’=S·cosβ，β=α12—α
同理可求出n’的坐标。然后由n、n’的坐标计算该点的移动变形量。该方法使用坐标观测量进行计算，避免了复杂的支距改正，提高了计算精度。
2.2移动和变形的计算
观测数据经规划后，便可计算观测线上各测点和测点间的移动和变形。移动和变形的计算主要包括：各测点的下沉和水平移动，相邻两测点间的倾斜和水平变形，相邻三点的曲率变形[2-4]。
2.2.1 下沉
各点的下沉以本次与首次的标高差表示，单位mm，则m次观测n点的下沉Wn为
Wn = Hn 0 — Hn m （2-3）
式中，
Hn 0，Hn m——分别为首次和m次观测时n点的高程。
2.2.2倾斜
地表下沉盆地沿某一方向的坡度叫倾斜，用i表示，单位mm/m，相邻两点间的倾斜为
（2-4）
式中，
ln~n+1——n号点至n+1号点的水平距离；
Wn，Wn+1——分别为n号点和n+1号点的下沉量。

2.2.3曲率
观测线下沉的弯曲程度叫曲率，用K表示，单位mm/m2，则n号点附近的曲率，即n -1点至n+1号点之间的曲率为
（2-5）
式中，
in+1~n，in~n-1——分别为n+1号点至n号点和n号点至 n-1号点的倾斜；
ln+1~n，ln~n-1——分别为n+1号点至n号点和n号点至 n-1号点的水平距离。
2.2.4水平移动
下沉盆地内，某点沿某一水平方向的位移叫水平移动，用U表示，单位mm，则m次观测n点的水平移动为
Un= Ln m — Ln 0 （2-6）
式中，
Ln m，Ln 0——分别为m 次观测时和首次观测时 n 号点至观测线控制点间的水平距离，直接用坐标求得。
2.2.5水平变形
下沉盆地内两点间单位长度的水平移动差称之为水平变形，以ε表示，单位为 mm/m。则 n 号点至 n+1 号点间的水平变形为
（2-7）
式中，
（l n+1 ~ n）m，（l n+1 ~ n）0——分别为n+1号点至n号点在首次和m 次观测时的水平距离。

3软件实现
3.1 AutoCAD VBA简介
VBA ( visual basic for application) 是由Microsoft创建的，用来自动执行任务的一个编程环境，它提供了一些用来创建图形用户界面的可拖拉工具和用来与AutoCAD对象交换的编程语言[5]。其最大特点在于两个方面：第一，与VB有着几乎相同的开发环境和语法，具备功能强大和易于掌握的特点；第二，在于它的for Application功能，针对性强，它驻留在主程序的内部，使其结构精简，且代码运行效率高，使语言的引擎在技术上与开发环境分离[6]。
AutoCAD自R14.01版本开始就支持VBA作为二次开发工具。
3.2实现步骤
基于上文所述计算原理，在AutoCAD 2008下使用VBA实现数据计算及成图。程序界面如图2 。
3.2.1 确定煤层采动信息
为反应地下开采与地表移动和变形规律，在进行各项计算之前应当首先确定观测数据所对应的地下开采情况。煤层信息有平均采深、采厚、倾角；采动信息有沿走向线上开切眼坐标、推进距离，沿倾向线上一边界坐标、倾向宽度。

图2 软件界面
3.2.2 确定观测线基准点坐标
在观测线两端选取两点，坐标可在cad观测线布置图中直接获取。计算观测线方向代码如下：
'读入基点坐标，计算规划方向
JX1 = Val(UserForm1.TextJZ1x.Text)
JY1 = Val(UserForm1.TextJZ1y.Text)
JX2 = Val(UserForm1.TextJZ2x.Text)
JY2 = Val(UserForm1.TextJZ2y.Text)
a = Atn((JY2 - JY1) / (JX2 - JX1))
3.2.3读入观测数据进行坐标规划
数据文件编辑成txt文档，首次观测平面坐标及高程在前，第m次观测在后。首次和第m次观测的平面坐标规划实现代码如下：
' 数组gX1、gY1为首次观测的规划后平面坐标，数组hX2、hY2为第m次观测的规划后平面坐标
gx(i) = Format(Sqr((CDbl(arr(1)) - JX1) ^ 2 + (CDbl(arr(2)) - JY1) ^ 2) * Cos(a) + JX1, "0.000")
gy(i) = Format(Sqr((CDbl(arr(1)) - JX1) ^ 2 + (CDbl(arr(2)) - JY1) ^ 2) * Sin(a) + JY1, "0.000")
ds(i) = Format(Sqr((gx(i) - JX1) ^ 2 + (gy(i) - JY1) ^ 2), "0.000")
hx(i) = Format(Sqr((CDbl(arr(4)) - JX1) ^ 2 + (CDbl(arr(5)) - JY1) ^ 2) * Cos(a) + JX1, "0.000")
hy(i) = Format(Sqr((CDbl(arr(4)) - JX1) ^ 2 + (CDbl(arr(5)) - JY1) ^ 2) * Sin(a) + JY1, "0.000")
hs(i) = Format(Sqr((hx(i) - JX1) ^ 2 + (hy(i) - JY1) ^ 2), "0.000")
ws(i) = hs(i) - ds(i)
3.2.4文件的打开，文件的导出与程序的退出
1）文件的打开代码
Private Sub CommandButton1_Click()
Dim strtext As String
CommonDialog1.Filter = "文本文件(*.txt)|*.txt"
CommonDialog1.ShowOpen
If CommonDialog1.FileName <> "" Then
TextBox17.Text = ""
Open CommonDialog1.FileName For Input As #1
Do While Not EOF(1)
Line Input #1, strtext
TextBox17.Text = TextBox17.Text & strtext & Chr(13) & Chr(10)
Loop
Close #1
End If
End Sub
2）文件的导出代码：
Private Sub CommandButton2_Click()
CommonDialog1.Filter = "文本文件(*.txt)|*.txt"
CommonDialog1.ShowSave
If CommonDialog1.FileName <> "" Then
Open CommonDialog1.FileName For Output As #1
Print #1, TextBox17.Text
Close #1
End If
End Sub
3)程序的退出代码：
Private Sub CommandButton4_Click()
End
End Sub
3.2.4 计算移动和变形
各点规划至观测线方向后，根据公式（2-3）~（2-6）计算出移动和变形。从界面窗口中可以看到移动和变形的计算信息，同时各计算结果会自动保存在txt文档里。
3.2.5 曲线图的绘制
首先依原比例绘出煤层采动信息，及煤层与测点的空间关系。煤层绘制的部分代码如下：
1) 煤层走向信息绘制
Private Sub OptionButton1_Click()
Dim layer1 As AcadLayer
Set layer1 = ThisDrawing.Layers.Add("坐标")
layer1.color = 255
ThisDrawing.ActiveLayer = layer1
Call MC
Call CKQTC
End Sub
2）煤层倾向信息绘制
Private Sub OptionButton2_Click()
Dim layer1 As AcadLayer
Set layer1 = ThisDrawing.Layers.Add("坐标")
layer1.color = 255
ThisDrawing.ActiveLayer = layer1
Call mcq
Call QXTC
End Sub
Public Sub MC()
Dim HZ#, ZD#, KX#, KY#, JX1#, JY1#, SKQ#, CH# '煤层中间点采深，平均采厚
'煤层走向推进距离,开切眼坐标，基准点坐标,开切眼至基准点距离
Dim plineObj As AcadLWPolyline
Dim points(0 To 19) As Double
HZ = Val(UserForm1.TextCS.Text)
CH = Val(UserForm1.TextCH.Text)
ZD = Val(UserForm1.textJL.Text)
KX = Val(UserForm1.textKX1.Text)
KY = Val(UserForm1.textKY1.Text)
JX1 = Val(UserForm1.TextJZ1x.Text)
JY1 = Val(UserForm1.TextJZ1y.Text)
'沿走向
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 = ThisDrawing.ModelSpace.AddLightWeightPolyline(points)
ZoomAll
End Sub
煤层走向代码如下：
Public Sub CKQTC()
Dim HZ#, ZD#, KX#, KY#, JX1#, JY1#, SKQ#, CH# '煤层中间点采深，平均采厚
'煤层走向推进距离,开切眼坐标，基准点坐标,开切眼至基准点距离
HZ = Val(UserForm1.TextCS.Text)
CH = Val(UserForm1.TextCH.Text)
ZD = Val(UserForm1.textJL.Text)
KX = Val(UserForm1.textKX1.Text)
KY = Val(UserForm1.textKY1.Text)
JX1 = Val(UserForm1.TextJZ1x.Text)
JY1 = Val(UserForm1.TextJZ1y.Text)
'沿走向
SKQ = Sqr((KX - JX1) ^ 2 + (KY - 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
' 定义实体
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
' 在模型空间中创建实体对象
Set solidObj = ThisDrawing.ModelSpace.AddSolid(point1, point2, point3, point4)
ZoomAll
End Sub
煤层倾向代码如下：
Public Sub mcq()
Dim HZ#, D#, QBX#, QBY#, JX1#, JY1#, Q#, CH#, j# '煤层中间点采深，平均采厚
'煤层走向推进距离,开切眼坐标，基准点坐标,开切眼至基准点距离
Dim plineObj As AcadLWPolyline
Dim points(0 To 19) As Double
HZ = Val(UserForm1.TextCS.Text)
CH = Val(UserForm1.TextCH.Text)
D = Val(UserForm1.QD.Text)
j = Val(UserForm1.QJ.Text)
QBX = Val(UserForm1.QX.Text)
QBY = Val(UserForm1.QY.Text)
JX1 = Val(UserForm1.TextJZ1x.Text)
JY1 = Val(UserForm1.TextJZ1y.Text)
'沿走向
Q = Sqr((QBX - JX1) ^ 2 + (QBY - JY1) ^ 2)
points(0) = 0: points(1) = -HZ
points(2) = 1000: points(3) = -HZ + 1000 * Tan(j * 3.1415926 / 180)
points(4) = 1000: points(5) = -HZ - CH / Cos(j * 3.1415926 / 180) + 1000 * Tan(j * 3.1415926 / 180)
points(6) = 0: points(7) = -HZ - CH / Cos(j * 3.1415926 / 180)
points(8) = 0: points(9) = -HZ
points(10) = Q: points(11) = -HZ + Q * Tan(j * 3.1415926 / 180)
points(12) = Q + D * Cos(j * 3.1415926 / 180): points(13) = -HZ + Q * Tan(j * 3.1415926 / 180) + D * Sin(j * 3.1415926 / 180)
points(14) = Q + D * Cos(j * 3.1415926 / 180): points(15) = -HZ + Q * Tan(j * 3.1415926 / 180) + D * Sin(j * 3.1415926 / 180) - CH / Cos(j * 3.1415926 / 180)
points(16) = Q: points(17) = -HZ + Q * Tan(j * 3.1415926 / 180) - CH / Cos(j * 3.1415926 / 180)
points(18) = Q: points(19) = -HZ + Q * Tan(j * 3.1415926 / 180)
Set plineObj = ThisDrawing.ModelSpace.AddLightWeightPolyline(points)
ZoomAll
End Sub
煤层倾向采空区填充
Public Sub QXT