有缓和曲线或无缓和曲线任意桩号中、边桩坐标4800计算程序
有缓和曲线或无缓和曲线任意桩号中、边桩坐标4800计算程序
四、QXZB(曲线坐标)
A”ZH(ZH)=〃:B”X(ZH)=〃:C”Y(ZH)=〃 ”X(HZ)=〃:E”Y(HZ)=〃:F”FWJ-1=〃:G”FWJ-2=〃:H”R=〃:I”L1=〃:J”L=〃:K”L2=〃:L”ZJ(Z=1,Y=-1)=〃
O=180/pi (pi为圆周率,电脑中找不到此符号以此代替)
LbI A
{MN}
M”ZHUANG HAO=〃
P=M-A
MM>(A+I+J+K)?GOTO A⊿
N”ZHUANG JU(Z+,Y-)=〃
M>(A+I+J)?GOTO D⊿
M>(A+I)?GOTO C⊿
LbI B(第一缓和曲线)
X=B+(P-(P5÷40÷H2÷I2))COSF+(P3÷6÷H÷I-P7÷336÷H3÷I3)LSINF
Y=C+(P-(P5÷40÷H2÷I2))SINF-(P3÷6÷H÷I-P7÷336÷H3÷I3)LCOSF
X=X+NCOS(F-(P2÷2÷H÷I)OL-90)◢
Y=Y+NSIN(F-(P2÷2÷H÷I)OL-90)◢
GOTO A
LbI C(圆曲线)
X=B+(HSIN(((P-I)÷H+I÷2÷H)O)+(I÷2-I3÷240÷H2))COSF-(H(1-COS(((P-I)÷H+I÷2÷H)O))+ I2÷24÷H)(-1)LSINF
Y=C+(HSIN(((P-I)÷H+I÷2÷H)O)+(I÷2-I3÷240÷H2))SINF+(H(1-COS(((P-I)÷H+I÷2÷H)O))+I2÷24÷H)(-1)LCOSF
X=X+NCOS(F-((P-I)÷H+I÷2÷H)OL-90)◢
Y=Y+NSIN(F-((P-I)÷H+I÷2÷H)OL-90)◢
GOTO A
LbI D(第二缓和曲线)
Q=A+I+J+K-M
X=D-(Q-Q5÷40÷H2÷K2)COSG+(Q3÷6÷H÷K-Q7÷336÷H3÷K3)LSING
Y=E-(Q-Q5÷40÷H2÷K2)SING-(Q3÷6÷H÷K-Q7÷336÷H3÷K3)LCOSG
X=X+NCOS(F-(I÷2÷H+J÷H+K÷2÷H-Q2÷2÷H÷K)OL-90)◢
Y=Y+NSIN(F-(I÷2÷H+J÷H+K÷2÷H-Q2÷2÷H÷K)OL-90)◢
GOTO A
变量说明:
以上变量依次为:
ZH(ZH)=直缓点桩号,无缓和曲线则为直圆点;
X(ZH)=直缓点X坐标;无缓和曲线则为直圆点;
Y(ZH)=直缓点Y坐标,无缓和曲线则为直圆点;
X(HZ)=缓直点X坐标,无缓和曲线则为圆直点;
Y(HZ)=缓直点Y坐标,无缓和曲线则为圆直点;
FWJ-1=第一切线方位角;
FWJ-2=第二切线方位角;
R=半径;
L1=第一缓和曲线长;
L=圆曲线长;
L2=第二缓和曲线长;
ZJ(Z=1,Y=-1)=转角形式(左=1,右=-1)
ZHUANG HAO=所求点桩号
ZHUANG JU(Z+,Y-)=边桩桩距,左边桩为正,右为负,中桩0
我已经用过了,没发现什么问题,请大家提意见。
MM>(A+I+J+K)?GOTO A⊿ 里面的“?”是什么?