CASIO5800P缓和曲线正反算
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
SpiraL curve（缓和曲线任意点坐标计算）
程序说明：本程序是基于计算对称性缓和曲线而编写的程序，是以交点法为基础，故输入数据多了点。程序思路：根据曲线的交点数据，进而推算曲线要素，并进行曲线内任意点的坐标计算。
LbL 0：Cls
19→DimZ变量赋值
“INTα(s)=”?F输入交点坐标方位角
“JD PEG=”?K输入交点桩号
“INT X(JD)=”?X输入交点X坐标
“INT Y(JD)=”?Y输入交点Y坐标
“INT R(s)=”?R输入缓和曲线半径
“INT L(s)=”?L输入缓和曲线长
“TURNING ANGLE=”?A输入转角
“TURNING DIRECTION –L,+R=”?I输入路线转向：左负右正 I赋值为1
L2÷(24×R)→P
L÷2-L3÷(240×R2)→Q
90×L÷（π×R）→B
（R+P）Tan（A÷2）+Q→T计算切线长
R（A-2B）×π÷180→O计算圆曲线长
R（A-2B）×π÷180+2L→C计算曲线长
（R+P）÷cos（A÷2）-R→E计算外距
2T-C→D计算切曲差
K-T→Z[3]ZH点
Z[3]+L→Z[4]HY点
Z[4]+O→Z[5]YH点
Z[5]+L→Z[6]HZ点
Z[6]-C÷2→Z[7]QZ点
X+Tcos（F+180）→Z[16] ZH点X坐标
Y+TsIn（F+180）→Z[17] ZH点Y坐标
L-L3÷（90R2）→G
Z[16]+Gcos(F+30×I×L÷(π×R) )→Z[8]HY点X坐标
Z[17]+Gsin(F+30×I×L÷(π×R)) →Z[9]HY点Y坐标
Z[8]+2Rsin(45×O÷(πR))cos(F+I×45×O÷(π×R)+I×90×L÷(π×R)) →Z[14] QZ点X坐标
Z[9]+2Rsin(45×O÷(πR))sin(F+I×45×O÷(π×R)+I×90×L÷(π×R)) →Z[15] QZ点Y坐标
Z[8]+2Rsin(90×O÷(πR))cos(F+I×90×O÷(π×R)+I×90×L÷(π×R)) →Z[10] YH点X坐标
Z[9]+2Rsin(90×O÷(πR))sin(F+I×90×O÷(π×R)+I×90×L÷(π×R)) →Z[11] YH点Y坐标
X+Tcos（F+I×A）→Z[12] HZ点X坐标
Y+TsIn（F+I×A）→Z[13] HZ点Y坐标
LbL 2
“XIAN SHI QXYS 0→YES,1→NO” →？Z:CLs显示曲线要素？输入：0→YES,1→NO If Z=0:Then GoTo 3
ELse If Z=1:Then GoTo 4:IfEnd:IfEnd
LbL 3 <本段为曲线要素显示结果部分>
“QIE XIAN C=”:T▲
“QU XIAN C=”:C▲
“Y QU XIAN=”:O▲
“WAI JU=”:E▲
“QIE QU CHA=”:D▲
“ZH PEG=”:Z[3] ▲
“X(ZH)=”:Z[16] ▲
“Y(ZH)=”:Z[17] ▲
“HY PEG=”:Z[4] ▲
“X(HY)=”:Z[8] ▲
“Y(HY)=”:Z[9] ▲
“QZ PEG=”:Z[7] ▲
“X(QZ)=”:Z[14] ▲
“Y(QZ)=”:Z[15] ▲
“YH PEG=”:Z[5] ▲
“X(YH)=”:Z[10] ▲
“Y(YH)=”:Z[11] ▲
“HZ PEG=”:Z[6] ▲
“X(HZ)=”:Z[12] ▲
“Y(HZ)=”:Z[13] ▲
GoTo 0
LbL 4 <本段为待求点桩号计算部分>
“INT P PEG=”?H输入待求点桩号
If H If H≥Z[3]And H ELse If H≥Z[4] And H ELse If H≥Z[5] And H≤Z[6]:Then GoTo C :Ifend Ifend : Ifend If H>Z[6]:Then GoTo 4:Ifend LbL A {待求点位于ZH到HY段上时} “PIAN JIAO -L ,+R = ” ? V输入边桩与路线中心线夹角：左负右正“PIAN JU = ”? W边桩与路中法线距离 H–Z[3]→J J–J5 ÷(90×R2×L2)→S Z[16]+S cos(F+30×I×J2÷(π×R×L)→M Z[17]+S sin(F+30×I×J2÷(π×R×L)→N M+W×cos(F+90×I×J2÷(π×R×L +V) →Z[18] N+W×sin(F+90×I×J2÷(π×R×L +v) →Z[19] “ X(p)= ”: Z[18] ▲ “ Y(p)= ”: Z[19] ▲ GoTo 4 LbL B {待求点位于HY到YH段上时} “PIAN JIAO -L ,+R = ” ? V “PIAN JU = ”? W H–Z[4]→J Z[8]+2Rsin(90×J÷(π×R)cos[F+90×I×J÷(π×R)+90×I×L÷(π×R)]→M Z[9]+2Rsin(90×J÷(π×R)sin[F+90×I×J÷(π×R)+90×I×L÷(π×R)]→N M+Wcos(F+90×I×L÷(π×R)+180×I×J÷(π×R)+V) →Z[18] N+Wsin(F+90×I×L÷(π×R)+180×I×J÷(π×R)+V) →Z[19] “X(P)= ”: Z[18] ▲ “Y(P)= ”: Z[19] ▲ GoTo 4 LbL C {待求点位于YH到HZ段时} “PIAN JIAO -L ,+R = ” ?V “PIAN JU = ”? W Z[6] –H →J