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

请使用WORD的Times New Roman字体，若使用其它字体，字母O与数字0是区别不出的.切记切记!

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

第十一版

1  (主程序) (使用时请运行这个)

L1  LbI 0:{ QDV}

L2  U<4=> Prog “4”△

L3  U=5=> Prog “2”△

L4  U=6=> Prog “2”△ 

L5  Int U=7=> D=0:V=0△

L6  U=8=>{M}: PSW:D=(B+M“01”P)÷Cos Abs(90-AbsV): Q=W+S√(1-M2)△

L7  U>5=> 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):J<0=>J=J+360:△A=J: Prog“9”:Goto 3△

L10  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 I÷1000: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=B÷2÷Sin Abs V:J=R:U=7.1=>D=B÷2: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”: V≤0=>V=V+180: ≠>V=V-180: △D=S:Goto 1△

L14 U<5=>Q=Q+B: △Goto 0

10（子程序）（为共用 计算中桩/加转角/用正算计算边桩/反算/显示 的过程）

L1  Prog “F”:U=11=>R=0:△J= R+V: Prog “6” : Prog “8” :  Prog “5”

9 （子程序）（为点垂桩垂距/涵交桩交角的共同运行部分）

L1  Prog“7”: Prog“8”: Abs(J-R)<1=>Q=Q-I:≠> Q=Q+I:△I<0.001=>T=P:K=S: Prog “8”: “OK”: Q: Pause 0:I: Pause 0: J=J-R:J>180=>J=-J+180:△J→DMS◢

U=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)+S

5  (子程序)( 显示部分)

L1  T: Pause 0:K: Pause 0:I: Pause 0:J→DMS◢

F  (子程序) (中桩计算部分)

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):Q>G =>R=R+π:△DEG:R=180R/π

F2  (子程序) (中桩计算部分)(使用时把F2改为F;同时把F改为F1)(这个子程序计算速度快,但在缓曲长到2或3百米且半径小到6或7十米时,计算误差有几厘米)

L1  Defm 7

L1  I=Q-H

L2  J= (Z-C)I÷(G-H)

L3  Z[1]=A+90(J÷8+2C)(I÷8)÷π

L4  Z[2]=A+90(J÷4+2C)(I÷4)÷π

L5  Z[3]=A+90(3J÷8+2C)(3I÷8)÷π

L6  Z[4]=A+90(J÷2+2C)(I÷2)÷π

L7  Z[5]=A+90(5J÷8+2C)(5I÷8)÷π

L8  Z[6]=A+90(6J÷8+2C)(6I÷8)÷π

L9  Z[7]=A+90(7J÷8+2C)(7I÷8)÷π

L10  R=A+90(J+2C) I÷π

L11  T=E+ AbsI(Cos A+4(CosZ[1]+ CosZ[3]+ CosZ[5]+ CosZ[7])+2(CosZ[2]+ CosZ[4]+ CosZ[6])+ CosR)÷24

L12  K=F+ AbsI(Sin A+4(SinZ[1]+ SinZ[3]+ SinZ[5]+ SinZ[7])+2(SinZ[2]+ SinZ[4]+ SinZ[6])+ SinR)÷24

8 (子程序) (座标反算部分)

(也可作为座标反算主程序独立运行, 方位角按RCL R SHIFT °′″ 显示)

L1  J=TK: Pol ( T-Y,K-O): Pause 0:J<0=>J=J+360△

6 (子程序) (边桩和座标正算计算部分) 

(也可作为座标正算主程序独立运行) 

L1  TK:T=T+ Rec (D,J): Pause 0:K=K+J

2 (子程序) (推算开挖/填筑桩部分)

{平原微丘一般推算1-2次,山岭重丘一般4-5次,能达到5厘米以内的精度,看你的熟练程度}

L1  {WSPB}：I= -W+ S◢    

L2  D= Abs IP+B◢

L3  {D}:D

3 (导线点数据库)(数据库中的数字请按你所在公路的数据替换)

L1  Y=21=>Y=36640.666:O=65111.666△

L2  ……  (省略号表示类同上行, 并非真在程序中输入省略号. 注意: 第一个Y为导线点编号, 第二个Y为导线点X座标,O为导线点Y座标, 数据库中的数字请按你所在公路的数据替换. 但注意: 导线点编号必须没字母也没小数;导线点X座标则必须有小数,若刚好没有小数,请加上0.0001并不影响精度)

4 (平曲线数据库) (数据库中的数字请按你所在公路的数据替换)

L1 Q≤444.978=>H=374.978:E=36771.643:F=64929.942:A=151°29′38″:C=0: G=444.978:Z=-600-1:  Goto 1△

L2  Q≤811.819=>H=444.978: E=36710.800:F=64964.534:A=148°09′06″:C=-600-1: G= 811.819:Z=-600-1: Goto 1△

L3  Q≤881.819=> H=811.819:E=36475.596:F=65238.599:A= 113°07′15″:C=-600-1: G=881.819:Z=0: Goto 1△

L4  ……  (省略号表示类同上行,并非真在程序中输入省略号)

L5  LbI 1  (注意:这行始终是平曲线数据库的最后一行, 之前行都是相同的)

   (为更好的理解这个平曲线数据库的结构,请对照上面数据看下图)

11 （子程序）（双圆垂轴隧道 断面开挖 超欠测量部分）

A:实测高程    H:设计高程     I:实测宽度    B:超欠值  (设棱镜半径0.03)
L1  {AHI}:A≤H=>B=110:≠>A≤H+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◢

平曲线数据库的数据图示:

(图画得不好,现丑了.图示是对称缓曲,但线元法程序对任意公路线形都适用的)

ZH HZ

0+374.978 0+881.819

X=36771.643 0

Y=64929.942 HY

151°29′38″ 0+444.978

0 X=36710.800

Y=64964.534 YH

148°09′06″ 0+811.819

-600-1 X=36475.596

Y=65238.599

113°07′15″

-600-1

使用说明:

{本程序是针对CASIO4850计算器的，其它的计算器请作相应修改}
{有方框的字为一个字符，不是几个字母，为各种计算状态名或函数名，如RAD与Cos等等，∫为积分符号，L1\ L2\ L3\L4等为程序行号不要在程序中输入, ……省略号表示类同上行,并非真在程序中输入省略号}
{如有其它疑问或进行技术交流，请打詹德军的电话015918402837欢迎交流!}

运行主程序 1   (按1 FILE进入)

U? (判断进入何种计算状态)

1为数据库输入原始数据计算中边桩，

2为数据库输入计算放线测量记录表，

3为数据库输入计算平面位置检查记录表，

4为人工输入原始数据计算中边桩，

5为人工输入推算开挖/填筑桩状态，

6为数据库输入推算开挖/填筑桩状态，

7为数据库输入计算斜交斜做涵洞/单跨桥七点放样状态 ，

7.1为数据库输入计算斜交正做涵洞七点放样状态 ，(斜交斜做涵洞的伸缩缝处的放样也可用)

8为数据库输入计算锥坡放样，

9为数据库输入反算线外点垂桩垂距，

9.1为数据库输入双圆垂轴隧道 断面开挖 超欠测量，

10为数据库输入求线外两点交路线的桩号与夹角， 

11为计算曲线桥平行布桩的桩/柱放样状态

(因为D=0时就是中桩,D≠0时就是边桩,所以不再单独设立中桩或边桩的计算状态,方便程序进一步开发)

（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座标

=所求桩号放样距离

=所求桩号放样方位角

8．B？步长（步长就是间隔长度，当要计算的是大量有相同距离间隔 且 相同中边距或都是中桩 时，步长就是好东西了，计算完第一个桩后所有桩，都只需闭着眼睛连按四次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厘米的精度，要看你的熟练程度)

1．W?实测大地高程

2．S?变坡点设计高程

3．=变坡点设计高程 与 实测大地高程 之高差(正为填方,负为挖方)

4．P?填方或挖方的坡比(若1:1.5侧输1.5,无论填挖坡比全为正)

5．B?中桩到变坡点宽度

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． P？椭圆单侧斜长轴

3． S？椭圆单侧斜短轴(当从近桥0点往远桥1点推桩号增加时，S输入负值。当从近桥0点往远桥1点推桩号减少时，S输入正值)

4． W？椭圆中心对应的中桩桩号

5． B？正交半路宽，包括耳墙的宽度

6． 01？等分比(必须是0/0.1/0.2至0.8/0.9/0.95/0.985/1,否则在0到1之外会显示程序行出错)(0为最靠近桥的点,1为最远离桥的点)

7． V？转角（上/下流水轴 与 中桩的前进切线方位角间的夹角,左负右正）

8． Y?  置仪点X座标 （也可以输入导线点编号，因为此时使用了导线点数据库）

9． O? 置仪点Y座标（当Y输入了导线点编号，则此行不显示，因为数据库已经装入了导线点的X与Y座标）

10.=在椭圆上对应等分点X座标 （注:四个结果同时显示于屏幕）

=在椭圆上对应等分点Y座标

=在椭圆上对应等分点放样距离

=在椭圆上对应等分点放样方位角

11.循环第5/6/9步可放完一个桥台的两侧锥坡, 按AC/ON中断后输入另一桥台数据

计算锥坡放样图示

(本图是以圆上桥锥坡画的,所以两桥台不平行)

W点(此台S为正值)后退方向桥台

1等分点

0等分点

1等分点

W点(此台S为负值)

前进方向桥台

(U=9)为数据库输入原始数据反算线外点垂桩垂距状态

(300米内推算时间直线3-6秒,曲线12-15秒)

1. U?9为进入数据库输入原始数据反算线外点垂桩垂距状态

2. Q?开始推算的桩号(离真正的垂直桩号不要太远,在同一平曲线内或平曲线两边的直线段最好)

3. P?线外点的X座标

4. S?线外点的Y座标

5. =OK(OK表示计算完成)（注:三个结果同时显示于屏幕）

=线外点对应路线的垂直桩号

=线外点对应路线的垂直距离

=垂直桩号切线前进方向 与 垂直桩号到线外点方向的转角(左转为负右转为正)

(U=9.1)为数据库输入双圆垂轴隧道 断面开挖 超欠测量

1. U?9.1为进入数据库输入双圆垂轴隧道 断面开挖 超欠测量状态

2. Q?开始推算的桩号(离真正的垂直桩号不要太远,在同一平曲线内或平曲线两边的直线段最好)

3. P?遂壁任意点的X座标

4. S?遂壁任意点的Y座标

5. =OK(OK表示计算完成)（注:三个结果同时显示于屏幕）

=遂壁任意点对应路线的垂直桩号

=遂壁任意点对应路线的垂直距离

=垂直桩号切线前进方向 与 垂直桩号到线外点方向的转角(左转为负右转为正)

6. A? 遂壁任意点的Z座标(高程)

7. H?第5步求出的遂壁任意点对应路线的垂直桩号的中桩设计高程

8. I?第5步求出的遂壁任意点对应路线的垂直距离(这里不用输入,只是显示)

9. =遂壁任意点超欠值(正为超挖值,负为欠挖值)

(10)为数据库输入原始数据求线外两点交路线的桩号与夹角(精确涵位)

(300米内推算时间直线3秒,曲线12-15秒)

1. U?10为数据库输入原始数据求线外两点交路线的桩号与夹角

2. Q? 开始推算的桩号(离真正的垂直桩号不要太远,在同一平曲线内或平曲线两边的直线段最好)

3. P?线外第一点的X座标

4. S?线外第一点的Y座标

5. W?线外第二点的X座标

6. B?线外第二点的Y座标

7. =OK(OK表示计算完成)（注:三个结果同时显示于屏幕）

=线外两点交路线的桩号

=线外第一点交路线的距离

=交点桩切线前进方向 与 交点桩到线外第一点方向的转角(左转为负右转为正)

（U=11）为计算曲线桥平行布桩的桩/柱放样状态

(千万不能用于,曲线桥平行布桩的盖梁放样; 曲线桥平行布桩的支座放样)

1． Q？桩/柱对应的中桩号

2． D？中桩至桩/柱的距离

3． V？中桩至桩/柱的方位角（注意：这里不再是转角）

4． 以后循环第1/2/3步