本程序以高斯四节点法为计算内核.docx

1、本程序以高斯四节点法为计算内核本程序以高斯四节点法为计算内核,多条线路正反算,中桩,边桩高程计算,超高计算,超高缓和计算,加宽计算,加宽缓和计算,边坡开口线计算,挡土墙坡脚线计算，考虑了中间绿化带的影响，适用与国家高速公路至乡村四级公路计算和放样，路基路面工程可以直接得出中边桩的设计三维坐标，去掉加宽和超高影响的计算困难，边坡和坡脚线计算可直接在边坡上提取坐标带入程序，经过所有设计因素的综合，得出更改边距并显示修改偏差后的坐标反算速度明显提高,只需要3秒.程序无错，可放心输入，另本程序可以增加隧道超欠挖计算子程序,非常方便,因为每个隧道的断面数据不一样,所以在此没有明确写出,有需要的可以联系我

2、.本人QQ76805071,只为交友.计算器主程序：ZHU-CHENG-XULbi0:“1，ZS=FS，2ZS，3FS，4XYSG，5。”？U:U=1Prog”ZS-FS”:进入公路三维程序U=2Porg”ZS”:进入坐标正算程序U=3Porg”FS”:进入坐标反算程序U=4Porg”XY SG”: 进入大地坐标转施工坐标程序Goto0:说明：计算器总的主程序，进入选择各种分支计算程序。1为公路三维计算，2为普通正算，3为普通反算，4为大地坐标转施工坐标。选择错误重新选择。此程序可以不用输入，只为给大家一个思路,可以把计算器所有程序集中到一个主程序内管理.公路三维部分正反算选择程序：ZS-FS

3、Deg:/设置角度模式20DimZ:/扩展变量“1LZXY，2XYLZ,3BIANPO-FY”？U： /正反算选择，正算选1，反算选2,坡口坡脚选3If U=1: ThenProg”ZS-XH”:IfEnd: /进入正算循环主程序IfU=2:Then Prog”FS-XH”:IfEnd:/进入反算循环主程序IfU=3:ThenProg”BIANPO-XH”:IfEnd: /进入边坡开挖主程序正算循环主体程序ZS-XH“1PT-2SJ”？W： /普通计算和设计边距计算选择“XL-XZ“？U： /选择线路 1NIf W=1: Then /运行普通计算模式 边踞自由输入Lbi0:”L=”?L:”Z=

4、”?Z:/输入桩号和边距Porg”ZS-XH-1“：Goto0:IfEnd:运行正算循环子程序If W=2:Then /运行设计计算模式 边距按设计输入Lbi1:”L=”? L：/输入桩号Porg”0。5LUKUAN”: /录入公路加宽后，左右侧设计边距Porg”ZS-XH-1“：Goto1:IfEnd:/运行正算循环子程序正算循环子程序ZS-XH-1Prog”PM-XL-XZ”:/进入平面线路选择程序对线路线元进行选择Prog”CHAOXIAN-CL”:/里程超限处理Prog”GL-ZS”:/进入坐标正算程序Prog”GL-BZ-ZS”:/进入公路边桩正算程序Prog”GC-XL-XZ”:/

5、进入高程线路选择程序对线路线元进行选择Prog”GC”: /进入高程部分提取中桩高程ProgCG-XL-XZ: /进入超高数据库提取超高数据If W=1:ThenProg”B-H”: /普通计算时进入边桩高程部分提取边桩高程ElseZ3 Z: Porg“B-H”：/按设计计算时候代入减去半边绿化带后的路基边距，提取边桩高程Ifend:Prog”ZS-XS”:/进入正算显示程序，提取正算三维结果反算循环主体程序FS-XH0Z:“XL-XZ”？U： /选择线路1NLbi0：“L”？L：“X=”？T：“Y=”？S：/输入起算假定桩号，实测X YProg”GL-FS”:/进入公路反算程序Prog”FS

6、-XS”:/进入反算显示程序Goto0:坐标正算程序GL-ZS.1739274226I:.3260725774J:.0694318442K:.3300094782E : /四节点法计算系数录入1-EF：1-KM：Q+KrO(C+KD) K: /数据库要用X、Y、Q作为起点坐标方位角Q+ErO(C+ED) E: /r是弧度,输入方法FUNCTION 选5Q+FrO(C+FD) F:Q+MrO(C+MD) M:X+O(ICos(K)+JCos(E)+JCos(F)+ICos(M) X: /计算出X坐标Y+O(ISin(K)+JSin(E)+JSin(F)+ISin(M) Y:/计算出Y坐标Q+Or

7、(C+D) Q：/计算出算点方位角,O是字母,r是弧度:边桩坐标计算程序GL-BZ-ZS:X+Zcos(Q+90) G:/计算出边桩X坐标Y+ZSin(Q+90) M:/计算出边桩Y坐标公路反算程序GL-FS20V:/V值设置为20，为第一次反算控制Do:Porg”PM-XL-XZ”: /进入平面线路选择程序Porg”CHAOXIAN-CL”: /里程超限处理程序If V10:/ V越小反算速度可能越慢，但是半径过小有可能算不出来，10为建议值Then Porg”GL-JY-ZS”: /进入公路简易正算二节点程序Else Porg”GL-ZS”:/进入公路正算4节点程序IfEnd:T-XI:S

8、-YJ:/试算点与测点X，Y距离差Pol(I,J):J-QJ: /方位角转换为接近试算点到测点的值：Rec(I,J): /计算出X差和Y差（X里程，J边）I/(1-CJ-2DJ) V:/计算出里程大致差距（此公式为加速反算）L+VL： /修改后的里程Abs(V) V: /修改值的绝对值LpWhile V0.001:/小于一毫米通过，大于一毫米循环计算中桩高程计算程序GCIf I-A0:Then1H:IfEnd: /开口向上部分赋值IfI-A=0:ThenGoto1: IfEnd: /无竖曲线部分赋值IfI-A0:Then-1H:IfEnd: /开口向下时赋值(L-C)I+K+(L-C)2/(2

9、E)HI： /计算得出竖曲线部分中桩高程Goto2:Lbi1:K+(L-C)II: /计算得出无竖曲线部分中桩高程Lbi2:边桩高程计算程序B-H（L-C）/（D-C）P:/计算点占本线元长比例3P(2(S)-2P(3(S) P:/计算得线性和高次抛物线系数F+P（T-F）E：/初步计算计算点横坡If V=1:Then If N= -1:Then/直线和直缓，缓直,曲线左偏时If Z0：ThenIfEZ4:ThenZ4E: IfEnd : /计算左侧时，横坡大于一般横坡取一般横坡I+Abs(Z)EO: IfEnd : /左侧高程IfZ0:Then I+ZFO:IfEnd: Ifend:/右侧高

10、程得数If N=1:Then IfZ0:Then I+Abs(Z)EO:IfEnd:/曲线右偏时，左侧高程得得IfZ0:Then If EZ4:Then Z4E:IfEnd/右侧高程横坡大于一般横坡，取一般横坡I+ZEO:IfEnd: IfEnd: IfEnd: /右侧高程得数IfV=2:ThenIf Z0:Then I+ZEO:IfEnd:/圆曲线全超高和S曲线超高渐变部分右侧高程得数IfZ0:ThenI-Abs(Z)EO:IfEnd:IfEnd: / 圆曲线全超高和S曲线超高渐变部分左侧高程得数备注:如果有更多的坡级,大家可以使用增加变量.完全满足大家边坡放样需求设计标准路面宽度子程序CH

11、ANGSHU12Z1：/录入12米为路面半宽度，根据设计修改1.5Z2：/录入1.5米为绿化带半宽，根据设计修改，无绿化带输入0-0.02 Z4: /录入-0。02为公路一般横坡，根据设计修改注：在以中轴线旋转的超高不考虑绿化带影响的，在常数程序中，绿化带修改为0路面宽度子程序0。5LUKUAN”Z-1,Y+1”?V:/左右边桩选择Porg”CHANG-SHU”:/设置标准路面半宽If V=1:Then Porg”GL-JK”/右侧时，进入加宽子程序Z1+ZZ: /计算右册正确公路半宽,Z-Z2 Z3：IfEnd：/计算右侧减去半边绿化带时计算边桩高程时用的边距IfV=-1: :Then Po

12、rg”GL-JK”/左侧时，进入加宽子程序-Z1-ZZ:/计算左册正确公路半宽Z+2 Z3:IfEnd:/计算左侧减去半边绿化带时计算边桩高程时用的边距公路加宽子程序GL-JK前点加宽里程A，加宽止点里程B，“加宽缓和段全长C，止点加宽宽度D本加宽程序欠缺S缓和曲线加宽程序，因为没有遇到过，所以暂时空缺Porg“JK-XL-XZ”: /进入加宽数据库线路选择程序If D=1:Then（4(L-A)/(B-A)3-3*(L-A)/(B-A)4）*CZ:IfEnd:/ 从小大到加宽If D=2:Then (4(B-A)-(L-A)/(B-A)3-3*（B-A)-(L-A)/(B-A)4)*CZ:I

13、fEnd:/从大到小加宽，以上两式为参考公式，可根据图纸更改公式IfD=3:ThenZ=Z：IfEnd: /圆曲线加宽直一样时IfZ0：then”JK=”:Z:IfEnd:/存在加宽值时显示加宽值正算显示ZS-XS“FWJ”：Q/显示计算里程方位角“ZHONG-X”:X/显示中线X坐标“ZHONG-Y”: Y /显示中线Y坐标“ZHONG-GC”:I /显示中桩高程“B-X”:G /显示边桩X坐标“B-Y”: M/显示边桩Y坐标“B-GC”:O /显示边桩高程反算显示FS-XS“L”：L/显示计算里程“Z”：J /显示计算边距里程超限CHAOXIAN-CLIf LA:Then “L-”；Sto

14、p:IfEnd:/小于本路线最小里程提示，退出程序If LB: Then”L+”:Stop:IfEnd:/大于本路线最大里程提示，退出程序公路简易正算程序GL-JY-ZS07886751346E：1-EF：Q+ErO(C+ED) E:Q+ErO(C+FD) F:X+O(Cos(E)+Cos(F)/2X:Y+O(Sin(E)+Sin(F)/2Y:Q+Or(C+D) Q:/Or,O为字母,r为弧度符号边坡放样程序BP-FY“XL-XZ“？U：/线路选择”S-L”?L:”X=”?T:”Y=”?S:/输入试算桩号和测点X，Y.TZ10:SZ11:Porg”PM-XL-XZ”: /对平面线路进行选择Po

15、rg”GL-FS“： /进入反算程序，反算测点里程偏距。If J0:Then”-YOU-“ /显示计算边坡为公路右侧Else”-ZUO-“Ifend: /显示计算边坡为公路左侧Porg”0.5lukuan”:/进入设计边距程序Porg”GC-XL-XZ”:/高程线路选择Porg”GC”:/对中桩高程进行计算PorgCG-XL-XZ /进入超高数据库提取超高数据Z3 Z：/录入减去半个隔离带加上加宽值后的边距Porg”B-H”: /计算该里程边桩高程Abs(Z)+Z2 Z:/录入设计绿化带Z10X：Z11Y： /坐标变量转换，抛弃中线点XY坐标，以第一测点XY为后步起算数据Porg”BP-XL-

16、XZ”:/进入边坡数据库,提取要素.“1，2，3，4-BP”T /选择所计算的边坡坡段Porg”BP-FY-1”: /进入边坡计算子程序1程序边坡计算子程序1BPFY-1O+CB+EFA: /计算出边坡底线高程Abs(Z)+B+D+EW:/计算中线点距离边坡底线距离 Abs(J)-WD :/试测点到起坡点平距VHC: /马道高差0B:/设置边坡底线平距为0If T=1:/计算第一段边坡Then Porg”BP-FY-2”:IfEnd: /进入计算子程序2A+G+CA: /第二段边坡底线高程H+GKB:/第二段边坡底距离第一段起线平距IfT=2：Then /进入计算第二段边坡MG：NK：/第二段

17、边坡数据变成子程序计算的第一边坡数据Porg”BP-FY-2“：IfEnd: /进入计算子程序2A+M+CA： /第三段边坡底线高程B+H+MNB/第三段边坡底距离第一段起线平距IfT=3：Then /进入计算第三段边坡IG：PK： /第三段边坡数据变成子程序计算的第一边坡数据Porg”BP-FY-2”:IfEnd: /进入计算子程序2A+I+CA:/第四段边坡底线高程B+H+IPB：/第四段边坡底距离第一段起线平距IfT=4：Then /进入计算第四段边坡RG：SK：/第四段边坡数据变成子程序计算的第一边坡数据Porg”BP-FY-2”:IfEnd: /进入计算子程序2边坡计算子程序3BP-

18、FY-3If G0:Then Return:IfEnd:IfNA: /如果测点高程小于该层边坡底线高程Then“+，SHANG-TIAO”：A-N /显示坡度错误，不在本坡段，上几米再次测量或计算低层边坡PorgBP-FY:IfEnd:/返回程序重新找点If N(A+G):/如果测点高程大于该段边坡定部高程Then”-，XIATIAO”: N-A-G/显示高程错误，不在本坡段，应在上面坡段，应忘上调整Porg: BP-FY”: IfEnd: /返回主程序边坡计算子程序2 BP-FY-2 “GC=”？N： /输入测点高程Porg”BP-FY-3”:/判断测点是否低于该坡段最低高程Abs(N-A)

19、*K+B-DE : /计算出正确偏差If Abs（E）0.05：Then/如果测点误差5厘米以内表示成功”CHENG-GONG,WUCHA=”:E/成功后，显示测点误差Prog”BP-FY”:IfEnd: /返回主程序测下个断面IfJ0:ThenEZ: /右侧录入误差Else -E Z: ZZ10 /左侧录入误差IfEnd: Porg”GL-BZ-ZS”: /计算下一点坐标“X=”：G“Y”：M /显示下一点坐标,进行放样GX：MY：/录入XY坐标,进行下次计算“GC=”N: /输入测点高程Abs(N-A)*K+B-(D+E)C: /计算出正确偏差If Abs(C)0.05：Then/如果测点

20、误差5厘米以内表示成功”CHENG-GONG,WUCHA=”:C/成功后，显示测点误差Prog”BP-FY”:IfEnd: /返回主程序测下个断面CE/(E-C) F:/由上次两点误差推算第三点IfJ0:ThenFZ:/右侧录入误差Else -FZ: IfEnd: / 左侧录入误差 Porg”GL-BZ-ZS”: /进行下一点坐标计算“X=”：G“Y=”：M /显示下一点坐标,进行放样GT：MS：/录入XY坐标,进行下次计算“NONO” :/本循环未计算出正确点Porg”BP-FY”: /，使用计算出的坐标，进入下个循环,程序自动把此坐标带入下个循环数据库部分1平曲线数据库平面线路选择程序PM-XL-XZU=1Prog”HUO-DA-XIAN”: /进入假定线路活大线。U=2Prog”X2”: /进入路线2U=3prog”X3”:/进入路线3（P-R）/（2HPR）D：L-OO：/正算程序需要补充部分ODD:1/PC： （数据库）活大线HUO-DA-XIAN 线路起点里程A：线路终点里程B：If L第一线元止点里程：Then 起点X