数值分析论文打印上交.doc

基于方程求根求解巷道间最短距离的研究分析 杨波 班级:研硕2012-3班学号：36 摘要：本文应用方程求根数值分析思想，结合mat lab程序算法，解算了煤矿两和邻巷道间 的最短距离，对巷道支护、缓解顶板来压起到重要的作用，也为巷道施工的安全提供了保障。 关键词：方程求根；数值分析；最短距离 ?引言 煤矿井下巷道数量众多，相互Z间关系错踪复杂。巷道之问有的直接相交连通，有的空 间交错。为了保证巷道施工和煤矿生产的安全，常常需要知道空间两个相邻巷道之间的距离 [h2|o尤其当相邻巷道之间距离较短时，准确地计算两巷Z间的最短距离对生产和安全更是 意义重大。 空间交错的两巷道之间，有且只有一?个最短距离，该最短距离所在的线段应是两巷道的 公垂线。计算的关键和难点在于找出公垂线的位置。过去，测量人员往往采用简易的方法来 粗略估算两巷交巷道之间的最短距离；有的先计算投影到平而上的交点，再把投影交点处的 两巷道高程之差当作最短距离；有的先在采掘工程平面图上粗略地确定出两巷道之间大致最 近的位置，然后再近似地计算两巷道最近点的坐标和高程，最后反算出两点间的空间距离。 这些方法确定的两巷道Z间的最短距离不是“最短”，很不准确，有的与准确值相去茯远。 所以，为了确保安全，两空间交错巷道之间的最短距离必须准确地计算岀来Q引。为此，本 论文推导出两交叉巷道Z间的最短距离公式，很大程度的提高了计算速度和准确度。 ?问题分析及模型 如图1所示,巷道11湘巷道22空间交错（在平面内并不相交）,其中，“A”和“B”分 别是2巷道中的任意点。 2 图I相邻巷道间的故短距离示慝图 在巷道ir中，1点的坐标高程为XI、Yl、HI,同时，还知道该巷道的方位角cd和倾 角61；在巷道22,中，2点的坐标高程为X2、Y2、H2,同时，也知道该巷道的方位角a2和 倾角62。现要计算出2巷道之间的最短距离D,即AB2点间的最短倾斜长度。公式推导如 下： 已知数据为： XI、Yl、H1 巷道起始点1的坐标和高程； al、81 巷道1J的方位角和倾角； X2、Y2、H2 巷道起始点2的坐标和高程； 血、62 巷道22的方位角和倾角。 设:S1代农巷道1』中1点A点的倾斜距离，S2代表巷道22中2点B点的倾斜距离， 则计算A、B两点三维坐标的公式如下： X a = X] +S] cos§ cose Ya =7(+5] cos J] sinof] I (1)Ha =H^S} sin J, J XB = X2+ S2 cos COS(72 「 Yh =Y2 +52cos J2sincr2 (2) Hb = H2 4-52sin(72 J 故AB两点间的距离: D2 = (X.a — X J +(乙-Yb)2+(Ha-Hb)2 即： D2 =(X| +S|Cos§ coseZ] 一 X2 -S2cosy2cos?2)2 +(乙 +S} cos5)s\na} -Y2 -S2cos J2sin?2)2 + (77i +S] s\n3} -H2-S2 sin^2)2 (4) 先对(4)式求一阶导数,并令其等于0,即组成1个以Si和S2为未知数的二元一次方程组: 6 —= 2cos爲 cosQ](X[ +S] cos/〕cosa} - X2 -S2 cosS2 cosa2) + 2cosSx sina}(Y} + 5 S] cos爲 sintZ] -Y2 - S2 cos^2 sina2) + 2sin(//( +S] sin5〕-H2 -S2sin^2) = 0 (5) d 2 = -2cos J2 cos6z2(Xj +S] cos J] cosof] - X2 -S2 cos J2 cos?2)-2cos J2 sina2(Y} + 久2 S] cos sinc^j - Y2 - S2 cos S2 sin a2) - 2sin S2 +S】sin -H2-S2 sin J2) = 0 (6) 令： k = cos § cos cos J2 cos a2 + cos § sin a} cos sin a2 +sin sin J2 [ m = cos J] cosG|(X| — X2)+ cos?] sine(Y| -y^ + sin/KHi — H?) r (7) Z?=COs52COSG2(X| -X2)+ COs52Si】]Q2(Y| — Y2)+sin》2(a —丹2) J 因： (cos J] COS6Z])2 +(cos5] sin?])2 +(sin^)2 = I 「 (8) (cosS2 cos?2)2 +(cos岛 sin?2)2 +(sin J2)2 = 1 」 解(5)式和(6)式纽成的方程纽，则有： m 一 nk a (9) ink 一 n S2=KF J