2. 海南矿业股份有限公司, 海南 昌江 572700
2. Hainan Mining Co., Ltd., Changjiang 572700, China.
Corresponding author: DING Hang-xing, E-mail: hangxingding@163.com
在地下矿山生产中,开拓和采准等巷道掘进工作直接关系到三级矿量的平衡,因此占有很重要的地位.一般按照“采掘并举,掘进先行”的原则进行生产组织,掘进要领先于采矿回采,这样就要求掘进的设计工作必须快速、准确[1].而在掘进的设计工作中,为了指导施工放线,需给出巷道坐标参数表,该表格按照传统手工方法绘制是相当繁琐的,从而影响到掘进设计工作的效率.
20世纪90年代以来,以AutoCAD为代表的计算机辅助设计软件迅猛发展.AutoCAD是一种通用软件,很难满足各行业的具体应用需求,为了解决该问题,该软件提供了开放式结构和多种开发工具,使用户能完全摆脱AutoCAD的限制,根据自己的需要对其进行二次开发[2, 3, 4, 5].针对掘进设计工作出现的问题,可采用对AutoCAD二次开发的方法加以解决,从而实现巷道坐标参数表的自动化生成,使广大采矿工程师从繁杂的计算中解脱出来,把更多的精力投入到重要的专业问题上去.
本文基于AutoCAD VBA二次开发平台,针对传统掘进工程施工设计的耗时环节——巷道坐标参数表的绘制,提出了自动化生成巷道坐标参数表的构想,编写了程序FME(free mining engineer).该程序以其便捷、鲁棒和准确的优点,被多座矿山引入应用.
1 巷道坐标参数的传统计算方法为了便于测量人员放线指导现场的掘进工作,一般对开拓、采准等掘进设计,会要求给出巷道坐标参数表,该表包含巷道中心线各拐点的三维坐标(X,Y,Z),以及前后拐点之间的水平距离、线长、坡度和方位角等内容.
为了简化说明以上参数的传统计算方法,以图 1所示两点为例,图中点Pn为点Pn+1的前一点,施工方向为从点Pn到点Pn+1.在实际设计中,两点的坐标是已知的,在没有计算机辅助的条件下,可以用直尺在绘制的图纸上逐点量取各点坐标和两点的水平投影距离dn.坡度i根据工程具体情况给出,线长ln(两点的实际距离)可根据下式计算:
式中:dn为Pn到Pn+1的水平投影距离,单位为m;i为坡度,用百分数表示,正、负数分别表示上、下坡.
方位角α为从北方向顺时针到施工方向形成的夹角,可用量角器在设计图中量取.以上仅叙述了两个点的计算方法,而在实际设计中,通常点的数量是很多的.由此可见,在没有计算机辅助的条件下,传统的巷道坐标参数计算方法工作量繁重,准确度也很难得到保证.
在计算机辅助下,设计人员工作量得以减轻.只需要在AutoCAD环境下,将图纸坐标和系统空间坐标对应上,可以逐点获取各拐点的坐标,量取两点的距离和方位角,再将得到的数据导入到Excel中,编辑公式进行两点线长的计算,在Excel中形成完整的巷道坐标参数表,最终再导入到AutoCAD中.由以上也可以看出,虽然工作量得以减轻,但仍没有达到自动化的程度,还需要反复在AutoCAD中量取坐标、距离和方位角等数据,并要和Excel程序反复切换.因此,为了达到将采矿工程师从繁杂的设计中解脱出来的目的,需要实现各参数的自动化获取和计算,并在AutoCAD中生成巷道坐标参数表.
2 参数自动生成的实现原理AutoCAD有强大且丰富的二次开发功能,可基于该二次开发平台,利用程序自动获取并计算相关参数,从而实现巷道坐标参数的自动生成.整个实现过程人工参与量很小,并且均在AutoCAD中完成,避免了软件反复切换造成的繁琐操作,可以大大提高生成效率.
基于以上思路,首先需要研究各参数的获取和计算方法.在AutoCAD环境下,巷道的中心线一般设计为多段线,各拐点的X,Y坐标可以从多段线中获取,水平距离dn可通过两点平面坐标进行计算:
由于坡度i是人工给出,线长ln可根据式(1)进行计算,点Pn+1的Z坐标可通过下式进行计算:
由于AutoCAD默认的角度计算是以E方向为起始方向,逆时针到施工方向进行计算.而施工设计中,方位角的计算是以N方向为起始方向,顺时针到施工方向进行计算.因此,为了自动计算方位角,需要建立两者的转换关系.
根据巷道的施工方向,一般可将方位角分为图 2所示的4种情况αn~αn+3,这4种情况涵盖了除垂直和水平方向以外所有方向.
对于图 2中的方位角,可以用图 1中的两点Pn和Pn+1进行简化计算.假定点Pn固定的话,点Pn+1的位置共有8种可能情况,分别为位于第一~第四象限和N,E,S,W4个方向轴上.研究主要还是以4个象限内的情况为主,而4个方向轴上的情况可作为特例进行计算.
如图 3所示,图中(1)~(4)分别对应点Pn+1位于第一~第四象限的情况.当点Pn+1在第一象限和N,E方向轴时,方位角αn和AutoCAD系统计算角α的关系为:αn=90°-α;当点Pn+1在第二~第四象限和W,S方向轴上时,方位角αn和AutoCAD系统计算角α的关系为:αn=450°-α.
通过以上各参数的获取和计算,说明巷道坐标参数表的自动生成是可行的,下一步只需要根据具体的开发工具,对以上原理进行实现.
3 VBA程序开发AutoCAD提供了丰富的二次开发工具,主要有Object-ARX,VB/VBA,ADS,AutoLISP/VisualLISP以及.NET[6].在这些开发工具里,VBA以开发速度快、开发界面友好和易学易用的特点,被工程技术人员广泛使用.AutoCAD VBA (visual basic for application) 是自AutoCAD R14开始嵌套在AutoCAD之中的一个基于对象的编程环境,提供一系列用于定义新类的宏及对AutoCAD现有类的扩充功能,还提供内容丰富的类库,能直接在AutoCAD内部执行[2, 7, 8, 9].因此,针对本次巷道坐标参数表自动生成的开发,采用VBA工具进行.
在VBA的开发环境中,设计将整个系统分为三部分,即数据的输入、处理和输出(图 4).
在数据的输入部分,首先需要完成巷道中心线的选取操作,然后输入起始点编号、坡度和起始点坐标.对于起始点编号,根据设计的实际情况,需要提供纯数字编号和字母数字混合编号两种样式以供选择;对于巷道坡度,通常情况下,在整条巷道上为一固定值,不过在某些情况下,存在局部坡度变化,因此,需要在输入坡度的基础上,提供坡度的局部修改接口;对于起始点坐标,为了简化操作,默认为自动填入系统坐标,同时也需要提供手动输入和选取的输入方式.最后需要对输出表格的样式进行设置(表格行距、字体大小和小数点保留位数).
在数据的处理部分,对选取的巷道中心线,从起点到终点依次将各拐点X,Y坐标存入二维数组中.为了便于下一步的计算,需要在数组的前面加一与起点数值相同的虚拟点.之后进行循环计算,得出各点的Z坐标,以及前、后点之间的水平投影距离、线长和方位角,最终与点编号、坡度形成10列的二维数组.
在数据的输出部分,根据设置的表格样式和表格放置点,对得到的数组依次进行输出,并绘制表格线形成巷道坐标参数表.
根据以上系统设计,进行了程序FME的编制,由于VBA不能注册AutoCAD命令,针对ThisDrawing对象,采用菜单的方式调用各命令(见图 5),主体程序对应菜单中的“巷道坐标”项.考虑到巷道转弯的情况,增加了对弯道参数的自动生成.FME程序自动生成的巷道坐标参数表如表 1所示,经过在海南铁矿和谦比希铜矿等多座矿山的应用,验证了该程序操作的便捷性、输出数据的准确性和运行的稳定性.
传统人工绘制巷道坐标参数表的方法,不仅量取和计算操作的工作量大,而且准确度难以得到保证.在AutoCAD和Excel的辅助下,同样存在较繁琐的操作,因此需要编制自动化程度较高的巷道坐标参数表生成程序.通过对各参数的计算分析发现,在AutoCAD环境下,所有参数均可以用计算的方法进行自动获取,繁琐的操作步骤也可以用程序进行实现,说明巷道参数表的自动生成是可行的,并提出了实现原理,基于AutoCAD VBA二次开发平台,编制了FME辅助设计程序.该程序经过多座矿山的实际应用,表现出较强的稳定性、便捷性和准确性,可极大地提高巷道掘进设计的效率.
[1] | 任凤玉,丁航行,张东红,等.缓倾斜层状岩体巷道断面形状与支护方式[J].东北大学学报:自然科学版,2011,32(1):125-128. (Ren Feng-yu,Ding Hang-xing,Zhang Dong-hong,et al.Cross-sectional geometry of drifts in gently sloping layered rock and support method[J].Journal of Northeastern University:Natural Science,2011,32(1):125-128.)(1) |
[2] | 高成慧,李燕.关于AutoCAD二次开发工具的探讨[J].现代计算机(专业版),2002(2):31-33. (Gao Cheng-hui,Li Yan.The discussion about the tools in the second development of AutoCAD[J].Modern Computer,2002(2):31-33.)(2) |
[3] | Senthilkumar R,Ramesh P,Velmurugan P.Modeling of optimal screw jack design[C]//International Conference on Frontiers in Automobile and Mechanical Engineering (FAME 2010).Chennai,2010:313-316.(1) |
[4] | Alwisy A,Al-Hussein M,Al-Jibouri S H.BIM approach for automated drafting and design for modular construction manufacturing[C]//Computing in Civil Engineering—Proceedings of the 2012 ASCE International Conference on Computing in Civil Engineering.Clearwater Beach,2012:221-228.(1) |
[5] | Brito A C,Oliveira C D,Marques P S.A simulation study of a multi-site production plant using arena[C]//European Simulation and Modelling Conference:Modelling and Simulation.Guimaraes,2011:125-129.(1) |
[6] | 王永辉,胡青泥,李红彩.AutoCAD二次开发方法的研究[J].计算机系统应用,2007(3):94-96+100. (Wang Yong-hui,Hu Qing-ni,Li Hong-cai.Research on the methods of AutoCAD secondary development[J].Computer Systems & Applications,2007(3):94-96+100.)(1) |
[7] | 孙卫强.基于VBA的AutoCAD参数化绘图[J].机械研究与应用,2006(6):100-102. (Sun Wei-qiang.Parameterized drawing in AutoCAD based on VBA[J].Mechanical Research & Application,2006(6):100-102.)(1) |
[8] | Bowen B.Microstation VBA for AutoCAD users[J].Cadalyst,2003,20(5):50.(1) |
[9] | Li J.The method for obtaining elevation points from the contour lines using AutoCAD VBA[J].Research and Exploration in Laboratory,2009,28(10):172-175.(1) |