积分公式计算的线元法CASIO4850计算器程序.doc

请使用WORD的TimesNewRoman字体，若使用其它字体，字母O与数字0是区别不出的.切记切记!积分公式计算的线元法CASIO4850计算器程序第十一版1 (主程序) (使用时请运行这个)L1 LbI 0: QDVL2 U Prog “4”L3 U=5= Prog “2”L4 U=6= Prog “2” L5 Int U=7= D=0:V=0L6 LbI 4:U=8=W: D=BSin AbsV+W“01”P:SL7 U5= Prog “4” L8 LbI2:IntU=9= PS: Prog“F”:A=R+90: Prog“9”:Goto 2 L9 LbI3:U=10=PSWB: Prog“F”: Pol ( P-W,S-B):JJ=J+360:A=J: Prog“9”:Goto 3L10 QDVHEFACZG: Frac Y=0= Prog “3”: O: Prog “10” L11 U=2=T=T+ Int (-5+10 Ran# )1000: K=K+ Int (-5+10Ran# )1000: Prog “8”: Prog “5” L12 U=3= I= -S+2S Ran# :J=(S2-I2): I= Int I1000:J= Int (-J+2J Ran# )1000:T+I: Pause 0:K+J: Pause 0:I: Pause 0:J(I2+J2) L13 LbI 1: IntU=7= DV: Prog“10”: S=D:D=B2Sin Abs V:J=R:U=7.1=D=B2:J=R+V-90: Prog “6”: Prog “8”: Prog “5”:D=2D:J=R-180: U=7.1= J=R+V+90: Prog “6”: Prog “8”: Prog “5”: V0=V=V+180: V=V-180: D=S:Goto 1L14 U=8= Goto 4L15 UQ=Q+B: Goto 010（子程序）（为共用 计算中桩/加转角/用正算计算边桩/反算/显示 的过程）L1 Prog “F”:U=11=R=0:J= R+V: Prog “6” :U=8=D=(1-W2) Abs S :SJ=180+R: J=R:Prog “6”: Prog “8” : Prog “5”9 （子程序）（为点垂桩垂距/涵交桩交角的共同运行部分）L1 Prog“7”: Prog“8”: Abs(J-R)Q=Q-I: Q=Q+I:IT=P:K=S: Prog “8”: “OK”: Q: Pause 0:I: Pause 0: J=J-R:J180=J=-J+180:JDMSU=9.1= Prog “11”: Prog “4”7 （子程序）（两点斜直线交点部分）(也可作为两点斜直线交点主程序独立运行)L1 E=Tan A:F= Tan R:Y=(EP-S-FT+K) (E-F): Pause 0:O=E(Y-P)+S5 (子程序)( 显示部分)L1 T: Pause 0:K: Pause 0:I: Pause 0:JDMSF (子程序) (中桩计算部分)L1 RAD:I= Abs(Q-H):J= Abs(G-H):R=A/180+CI+I2(Z-C)/2/J:T=E+(Cos (A/180+CX+X2(Z-C)/2/J),0,I,4):K= F+(Sin (A/180+CX+X2(Z-C)/2/J),0,I,4):QG =R=R+:DEG:R=180R/F2 (子程序) (中桩计算部分)(使用时把F2改为F;同时把F改为F1)(这个子程序计算速度快,但在缓曲长到2或3百米且半径小到6或7十米时,计算误差有几厘米)L1 Defm 7L1 I=Q-HL2 J= (Z-C)I(G-H)L3 Z1=A+90(J8+2C)(I8)L4 Z2=A+90(J4+2C)(I4)L5 Z3=A+90(3J8+2C)(3I8)L6 Z4=A+90(J2+2C)(I2)L7 Z5=A+90(5J8+2C)(5I8)L8 Z6=A+90(6J8+2C)(6I8)L9 Z7=A+90(7J8+2C)(7I8)L10 R=A+90(J+2C) IL11 T=E+ AbsI(Cos A+4(CosZ1+ CosZ3+ CosZ5+ CosZ7)+2(CosZ2+ CosZ4+ CosZ6)+ CosR)24L12 K=F+ AbsI(Sin A+4(SinZ1+ SinZ3+ SinZ5+ SinZ7)+2(SinZ2+ SinZ4+ SinZ6)+ SinR)248 (子程序) (座标反算部分)(也可作为座标反算主程序独立运行, 方位角按RCL R SHIFT 显示)L1 J=TK: Pol ( T-Y,K-O): Pause 0:JJ=J+3606 (子程序) (边桩和座标正算计算部分) (也可作为座标正算主程序独立运行) L1 TK:T=T+ Rec (D,J): Pause 0:K=K+J2 (子程序) (推算开挖/填筑桩部分)平原微丘一般推算1-2次,山岭重丘一般4-5次,能达到5厘米以内的精度,看你的熟练程度L1 WSPB：I= -W+ S L2 D= Abs IP+BL3 D:D3 (导线点数据库)(数据库中的数字请按你所在公路的数据替换)L1 Y=21=Y=36640.666:O=65111.666L2 (省略号表示类同上行, 并非真在程序中输入省略号. 注意: 第一个Y为导线点编号, 第二个Y为导线点X座标,O为导线点Y座标, 数据库中的数字请按你所在公路的数据替换. 但注意: 导线点编号必须没字母也没小数;导线点X座标则必须有小数,若刚好没有小数,请加上0.0001并不影响精度)4 (平曲线数据库) (数据库中的数字请按你所在公路的数据替换)L1 Q444.978=H=374.978:E=36771.643:F=64929.942:A=1512938:C=0: G=444.978:Z=-600-1: Goto 1L2 Q811.819=H=444.978: E=36710.800:F=64964.534:A=1480906:C=-600-1: G= 811.819:Z=-600-1: Goto 1L3 Q881.819= H=811.819:E=36475.596:F=65238.599:A= 1130715:C=-600-1: G=881.819:Z=0: Goto 1L4 (省略号表示类同上行,并非真在程序中输入省略号)L5 LbI 1 (注意:这行始终是平曲线数据库的最后一行, 之前行都是相同的) (为更好的理解这个平曲线数据库的结构,请对照上面数据看下图)11 （子程序）（双圆垂轴隧道 断面开挖 超欠测量部分）A:实测高程H:设计高程 I:实测宽度B:超欠值 (设棱镜半径0.03)L1 AHI:AH=B=110:AH+2.75+1.056=B=(A-(H+1.905)2+(I+1.813)2) +0.03-4.5:B = (A-(H+2.75)2+I2)+0.03-2.5: B平曲线数据库的数据图示:(图画得不好,现丑了.图示是对称缓曲,但线元法程序对任意公路线形都适用的)ZHHZ0+374.9780+881.819X=36771.6430Y=64929.942HY15129380+444.9780X=36710.800Y=64964.534YH14809060+811.819-600-1X=36475.596Y=65238.5991130715-600-1使用说明:本程序是针对CASIO4850计算器的，其它的计算器请作相应修改有方框的字为一个字符，不是几个字母，为各种计算状态名或函数名，如RAD与Cos等等，为积分符号，L1L2L3L4等为程序行号不要在程序中输入, 省略号表示类同上行,并非真在程序中输入省略号如有其它疑问或进行技术交流，请打詹德军的电话015918402837欢迎交流!运行主程序 1 (按1 FILE进入)U?(判断进入何种计算状态)1为数据库输入原始数据计算中边桩，2为数据库输入计算放线测量记录表，3为数据库输入计算平面位置检查记录表，4为人工输入原始数据计算中边桩，5为人工输入推算开挖/填筑桩状态，6为数据库输入推算开挖/填筑桩状态，7为数据库输入计算斜交斜做涵洞/单跨桥七点放样状态 ，7.1为数据库输入计算斜交正做涵洞七点放样状态 ，(斜交斜做涵洞的伸缩缝处的放样也可用)8为数据库输入计算锥坡放样，9为数据库输入反算线外点垂桩垂距，9.1为数据库输入双圆垂轴隧道 断面开挖 超欠测量，10为数据库输入求线外两点交路线的桩号与夹角， 11为计算曲线桥平行布桩的桩/柱放样状态(因为D=0时就是中桩,D0时就是边桩,所以不再单独设立中桩或边桩的计算状态,方便程序进一步开发)（U=1）数据库输入原始数据计算中边桩1. U？1(判断进进入何种计算状态)2. Q？所求桩号（计算第二个桩开始Q已经自动加了步长，不用手工输入，但仍然会显示你看，并可修改成别的桩号）3. D？中桩到边桩的距离（求中桩输入0）4. V？中桩到边桩方位角与中桩前进切线方位角间的夹角(注意左转为负值,角度按需要以度分秒输入任意输入)5. Y?置仪点X座标 （也可以输入导线点编号，因为此时使用了导线点数据库）6. O?置仪点Y座标（当X座标输入了导线点编号，则此行不显示，因为数据库已经装入了导线点的X与Y座标）7. =所求桩号X座标 （注:四个结果同时显示于屏幕）=所求桩号Y座标=所求桩号放样距离=所求桩号放样方位角8B？步长（步长就是间隔长度，当要计算的是大量有相同距离间隔 且 相同中边距或都是中桩 时，步长就是好东西了，计算完第一个桩后所有桩，都只需闭着眼睛连按四次EXE，就是下个桩要计算的所有结果了）9以后循环第2/3/4/7步（U=2）计算放线测量记录表1 U？2(判断进进入何种计算状态)2 Q？所求桩号（计算第二个桩开始Q已经自动加了步长，不用手工输入，但仍然会显示你看，并可修改成别的桩号）3 D？中桩到边桩的距离（求中桩输入0）4 V？中桩到边桩方位角与中桩前进切线方位角间的夹角(注意左转为负值,角度按需要以度分秒输入任意输入)5 Y?置仪点X座标 （也可以输入导线点编号，因为此时使用了导线点数据库）6 O?置仪点Y座标（当X座标输入了导线点编号，则此行不显示，因为数据库已经装入了导线点的X与Y座标）7 =所求桩号计算X座标 （注:四个结果同时显示于屏幕）=所求桩号计算Y座标=所求桩号计算放样距离=所求桩号计算放样方位角8 =所求桩号复测检查X座标 （注:四个结果同时显示于屏幕）=所求桩号复测检查Y座标=所求桩号复测检查放样距离=所求桩号复测检查放样方位角9 B？步长（步长就是间隔长度）10 以后循环第2/3/4/7/8步（U=3）计算平面位置检查记录表1 U？3(判断进进入何种计算状态)2 Q？所求桩号（计算第二个桩开始Q已经自动加了步长，不用手工输入，但仍然会显示你看，并可修改成别的桩号）3 D？中桩到边桩的距离（求中桩输入0）4 V？中桩到边桩方位角与中桩前进切线方位角间的夹角(注意左转为负值,角度按需要以度分秒输入任意输入)5 Y?置仪点X座标 （也可以输入导线点编号，因为此时使用了导线点数据库）6 O?置仪点Y座标（当X座标输入了导线点编号，则此行不显示，因为数据库已经装入了导线点的X与Y座标）7=所求桩号设计X座标 （注:四个结果同时显示于屏幕，但后两个数据是用不上的）=所求桩号设计Y座标=所求桩号放样距离=所求桩号放样方位角8 ？所允许的最大误差（以毫米为单位，比如高速公路的中线允许的最大误差为20MM，就输入20）9=所求桩号检测X座标 （注:四个结果同时显示于屏幕）=所求桩号检测Y座标=所求桩号X轴偏差=所求桩号Y轴偏差10=所求桩号双轴偏差11 B？步长（步长就是间隔长度）12 以后循环第2/3/4/7/9/10步（U=4）人工输入原始数据计算中边桩1 Q？所求桩号（计算第二个桩开始Q已经自动加了步长，不用手工输入，但仍然会显示你看，并可修改成别的桩号）2 D？中桩到边桩的距离（求中桩输入0）3 V？中桩到边桩方位角与中桩前进切线方位角间的夹角(注意左转为负值,角度按需要以度分秒输入任意输入)4 H?起算点桩号5 E?起算点X座标6 F?起算点Y座标7 A?起算点方位角(以度分秒输入)8 C?起算点曲率(即半径分之一，直线为0，注意左转为负值)9 Z?终算点曲率(即半径分之一，直线为0，注意左转为负值)10. G?终算点桩号11. Y?置仪点X座标12. O?置仪点Y座标13. =所求桩号X座标 （注:四个结果同时显示于屏幕）=所求桩号Y座标=所求桩号放样距离=所求桩号放样方位角14. ？步长（步长就是间隔长度）15. 以后循环第1/2/3/13步（U=5）人工输入原始数据推算开挖/填筑桩状态(平原微丘一般推1-2次，山岭重丘4-5次达到5厘米的精度，要看你的熟练程度)1W?实测大地高程2S?变坡点设计高程3=变坡点设计高程 与 实测大地高程 之高差(正为填方,负为挖方)4P?填方或挖方的坡比(若1:1.5侧输1.5,无论填挖坡比全为正)5B?中桩到变坡点宽度6=推算的中桩到开挖/填筑桩距离7. D?采用的中桩到开挖/填筑桩距离7 Q?所求开挖/填筑桩桩号8 V?中桩到边桩的方位角与中桩的前进切线方位角间的夹角(注意左转为负值,只能输入右边90或左边-90，因为开挖/填筑桩总是垂直的)9 H?起算点桩号10.E?起算点X座标11.F?起算点Y座标12.A?起算点方位角(以度分秒输入)13.C?起算点曲率(即半径分之一，直线为0，注意左转为负值)14.Z?终算点曲率(即半径分之一，直线为0，注意左转为负值)15.G?终算点桩号16. Y?置仪点X座标17. O?置仪点Y座标18.=推算开挖/填筑桩桩号X座标 （注:四个结果同时显示于屏幕）=推算开挖/填筑桩桩号Y座标=推算开挖/填筑桩桩号放样距离=推算开挖/填筑桩桩号放样方位角12以后循环第1/2/3/4/5/6/7/8/18步（U=6）数据库输入原始数据推算开挖/填筑桩状态(平原微丘一般推1-2次，山岭重丘4-5次达到5厘米的精度，要看你的熟练程度)1 W?实测大地高程2 S?变坡点设计高程3 =变坡点设计高程 与 实测大地高程 之高差(正为填方,负为挖方)4 P?填方或挖方的坡比(若1:1.5侧输1.5,无论填挖坡比全为正)5 B?中桩到变坡点宽度6 =推算的中桩到开挖/填筑桩距离7 D?采用的中桩到开挖/填筑桩距离8 Q?所求开挖/填筑桩桩号9 V?中桩到边桩的方位角与中桩的前进切线方位角间的夹角(注意左转为负值,只能输入右边90或左边-90，因为开挖/填筑桩总是垂直的)10 Y?置仪点X座标 （也可以输入导线点编号，因为此时使用了导线点数据库）11 O?置仪点Y座标（当Y输入了导线点编号，则此行不显示，因为数据库已经装入了导线点的X与Y座标）12 =推算开挖/填筑桩桩号X座标 （注:四个结果同时显示于屏幕）=推算开挖/填筑桩桩号Y座标=推算开挖/填筑桩桩号放样距离=推算开挖/填筑桩桩号放样方位角12以后循环第1/2/3/4/5/6/7/8/11步推算开挖/填筑桩的图示(以有二级填方边坡的填筑线为例,其它类推)中桩S变坡点高程P本级边坡坡比 W大地点高程 中桩至该变坡点距离 B地面线D中桩至该大地点距离（U=7/7.1）为计算斜交斜做涵洞单跨桥 / 斜交正做涵洞七点放样状态(U=7/7.1并能避免曲线内单跨桥涵进出水口长度不一至)(U=7同时适用于,曲线桥扇形布桩的盖梁放样; 曲线桥扇形布桩的支座放样)1.U? 7/7.1进入计算斜交斜做涵洞单跨桥/ 斜交正做涵洞七点放样状态2.Q? 涵洞桩号3.Y? 置仪点X座标 （也可以输入导线点编号，因为此时使用了导线点数据库）4. O?置仪点Y座标（当Y输入了导线点编号，则此行不显示，因为数据库已经装入了导线点的X与Y座标）5.=涵洞中桩X座标 （注:四个结果同时显示于屏幕）=涵洞中桩Y座标=涵洞中桩放样距离=涵洞中桩放样方位角6.D?右侧涵斜长7.V? 右侧转角(进流水轴 与 中桩的前进切线方位角间的夹角,左负右正)8.= 右侧桩X座标 （注:四个结果同时显示于屏幕）=右侧桩Y座标=右侧桩放样距离=右侧桩放样方位角9.B?正的跨径长度(程序自算斜跨径长度)10. =右侧前进侧台身脚X座标 （注:四个结果同时显示于屏幕）=右侧前进侧台身脚Y座标=右侧前进侧台身脚放样距离=右侧前进侧台身脚放样方位角11. =右侧后退侧台身脚X座标 （注:四个结果同时显示于屏幕）=右侧后退侧台身脚Y座标=右侧后退侧台身脚放样距离=右侧后退侧台身脚放样方位角12.D?左侧涵斜长13.V? 左侧转角(出流水轴 与 中桩的前进切线方位角间的夹角,左负右正,程序自动按180度转换了,其实不用输入,只是显示观看)14.= 左侧桩X座标 （注:四个结果同时显示于屏幕）=左侧桩Y座标=左侧桩放样距离=左侧桩放样方位角15. =左侧前进侧台身脚X座标 （注:U=7.1时为后退侧结果，四个结果同时显示于屏幕）=左侧前进侧台身脚Y座标=左侧前进侧台身脚放样距离=左侧前进侧台身脚放样方位角16. =左侧后退侧台身脚X座标 （注: U=7.1时为前进侧结果，四个结果同时显示于屏幕）=左侧后退侧台身脚Y座标=左侧后退侧台身脚放样距离=左侧后退侧台身脚放样方位角17.以后循环第6-8/10-17步(当然这种循环只是重新检查而已,按两次AC/ON退出)计算涵洞单跨桥七点放样图示（注: U=7.1时左侧的15步与16步顺序互换）前进方向10步15步 8步 右侧左侧14步5步11步16步切线路线（U=8）为数据库输入原始数据计算锥坡放样1 U?8为进入数据库输入原始数据计算锥坡放样2 B？正交半路宽，包括耳墙的正宽度3 V？转角（上/下流水轴 与 中桩的前进切线方位角间的夹角,左负右正）4 01？等分比(必须是0/0.1/0.2至0.8/0.9/0.95/0.985/1,否则在0到1之外会显示程序行出错)(0为最靠近桥的点,1为最远离桥的点)5 P？椭圆单侧斜长轴6 S？椭圆单侧斜短轴(当从远桥1点往近桥0点推桩号增加时，S输入正值。当从远桥1点往近桥0点推桩号减少时，S输入负值)7 Q？左右两椭圆锥顶连线对应的中桩桩号8 Y? 置仪点X座标 （也可以输入导线点编号，因为此时使用了导线点数据库）9 O?置仪点Y座标（当Y输入了导线点编号，则此行不显示，因为数据库已经装入了导线点的X与Y座标）10.=在椭圆上对应等分点X座标 （注:四个结果同时显示于屏幕）=在椭圆上对应等分点Y座标=在椭圆上对应等分点放样距离=在椭圆上对应等分点放样方位角计算锥坡放样图示(本图是以圆上桥锥坡画的,所以两桥台不平行)Q点(此台S为正值)后退方向桥台1等分点0等分点1等分点Q点(此台S为负值)前进方向桥台(U=9)为数据库输入原始数据反算线外点垂桩垂距状态(300米内推算时间直线3-6秒,曲线12-15秒)1. U?9为进入数据库输入原始数据反算线外点垂桩垂距状态2. Q?开始推算的桩号(离真正的垂直桩号不要太远,在同一平曲线内