卡西欧5800交点法计算缓和曲线内任一点.doc

CASIOfx-5800p交点法计算缓和曲线内任意点计算程序SpiraLcurve（缓和曲线任意点坐标计算）程序说明：本程序是基于计算对称性缓和曲线而编写的程序，是以交点法为基础，故输入数据多了点。程序思路：根据曲线的交点数据，进而推算曲线要素，并进行曲线内任意点的坐标计算。LbL0：Cls19→DimZ变量赋值“INTα(s)=”?F输入交点坐标方位角“JDPEG=”?K输入交点桩号“INTX(JD)=”?X输入交点X坐标“INTY(JD)=”?Y输入交点Y坐标“INTR(s)=”?R输入缓和曲线半径“INTL(s)=”?L输入缓和曲线长“TURNINGANGLE=”?A输入转角“TURNINGDIRECTION–L,+R=”?I输入路线转向：左负右正I赋值为1L2÷(24×R)→PL÷2-L3÷(240×R2)→Q90×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）→GZ[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坐标LbL2“XIANSHIQXYS0→YES,1→NO”→？Z:CLs显示曲线要素？输入：0→YES,1→NOIfZ=0:ThenGoTo3ELseIfZ=1:ThenGoTo4:IfEnd:IfEndLbL3<本段为曲线要素显示结果部分>“QIEXIANC=”:T▲“QUXIANC=”:C▲“YQUXIAN=”:O▲“WAIJU=”:E▲“QIEQUCHA=”:D▲“ZHPEG=”:Z[3]▲“X(ZH)=”:Z[16]▲“Y(ZH)=”:Z[17]▲“HYPEG=”:Z[4]▲“X(HY)=”:Z[8]▲“Y(HY)=”:Z[9]▲“QZPEG=”:Z[7]▲“X(QZ)=”:Z[14]▲“Y(QZ)=”:Z[15]▲“YHPEG=”:Z[5]▲“X(YH)=”:Z[10]▲“Y(YH)=”:Z[11]▲“HZPEG=”:Z[6]▲“X(HZ)=”:Z[12]▲“Y(HZ)=”:Z[13]▲GoTo0LbL4<本段为待求点桩号计算部分>“INTPPEG=”?H输入待求点桩号IfH<Z[3]:ThenGoTo4:IfEndIfH≥Z[3]AndH<Z[4]:ThenGoToAELseIfH≥Z[4]AndH<Z[5]:ThenGoToBELseIfH≥Z[5]AndH≤Z[6]:ThenGoToC:IfendIfend:IfendIfH>Z[6]:ThenGoTo4:IfendLbLA{待求点位于ZH到HY段上时}“PIANJIAO-L,+R=”?V输入边桩与路线中心线夹角：左负右正“PIANJU=”?W边桩与路中法线距离H–Z[3]→JJ–J5÷(90×R2×L2)→SZ[16]+Scos(F+30×I×J2÷(π×R×L)→MZ[17]+Ssin(F+30×I×J2÷(π×R×L)→NM+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]▲GoTo4LbLB{待求点位于HY到YH段上时}“PIANJIAO-L,+R=”?V“PIANJU=”?WH–Z[4]→JZ[8]+2Rsin(90×J÷(π×R)cos[F+90×I×J÷(π×R)+90×I×L÷(π×R)]→MZ[9]+2Rsin(90×J÷(π×R)sin[F+90×I×J÷(π×R)+90