缓和曲线与原曲线任意点坐标计算程序

Z=90S÷(R∏):U=S2÷(24R)-S^3÷(2688R^3):Q=S÷2-S^3÷(240R2):T=(R＋U)tan(N÷2)＋Q▲　　

L=R(N-2Z)∏÷180▲　　

Lbl0:{G}:K=Abs(G-F):S=0=＞Goto2: ≠＞Goto1⊿　　

Lbl1:K≤S=＞C=K-K^5÷(40R2S2) +K^9÷(3456R^4S^4):V=K^3÷(6RS)-K^7÷(336R^3S^3)+K^11÷(42240R^5S^5):O=tan-1(V÷C): Goto3: ≠＞Goto2⊿　　

Lbl2: O=180(K-S)÷(R∏)+Z:C=RsinO+Q:V=R(1-cosO)+U　　

Lbl3:Rec(C,B):X=D+I:Y=P+J: Rec(V,B):X=X+I:Y=Y+J▲　　

{H,W}:A=B+MO:ReC(H,A+90W):X=X+I▲Y=Y+J▲　　

Goto0　　

S：缓和段长度　　

R：圆曲线　　

∏：圆周率　　

T：缓和曲线切线长　　

N：缓和曲线偏角（输入绝对值）　　

L：圆曲线长　　

K：距离缓直点或圆起点距离　　

G：计算里程桩号F：缓直点（直缓点）或直圆点（圆直点）里程桩号　　

D：缓直点或圆起点X坐标P：缓直点（直缓点）或圆起点（圆终点）Y坐标　　

B：前切线方位角（计算缓直至圆缓或直圆至圆直用）或后切线方位角（计算直缓至缓圆或圆直至直圆用）　　

M：缓和曲线偏角方向控制（左输“　-1”　，右输“　1”　）（反向计算时相应反号）　　

X：计算点X值　　

Y：计算点Y值　　

H：边桩至中桩距离W：计算左边桩坐标时，输“　-1”　。　　

适用于计算：1、从缓直点至圆缓点或直缓点至缓圆点坐标　　

2、从直圆点至圆直点或圆直点至直圆点坐标，此时S=0。　　

本人编的4800计算器程序，有什么不到之处，恳请各位朋友给出宝贵的意见，本人一定予以改进。在此多谢各位！　　

本文来自：中国路桥网（www.9to.com），原文地址：http://www.9to.com/Article/Technique/daolu/200902/29227.html