单圆曲线中边桩坐标计算程序PM4-3

单圆曲线中边桩坐标计算程序PM4-3

单交点单圆曲线中边桩坐标计算程序(PM4-3)

(2) 主程序——PM4-3，占用内存1070字节。 "SINGLE CIRCLE CURVE" 显示程序标题1 "METHOD OF COORDINATE PM4-3" 显示坐标法标题2

Deg:ClrStat:FreqOn:Fix 3 基本设置 14

DimZ 定义额外变量

"JD MILEAGE PEG(m)="?Z 输入以m为单位的交点桩号

"JD X(m)="?U 输入交点x坐标

"JD Y(m)="?V 输入交点y坐标 "ZD X(m)="?M 输入转点x坐标

"ZD Y(m)="?N 输入转点y坐标 Pol(U-M,V-N):Cls 计算ZDJD的边长与方位角 If J<0:Then J+360A:Else J

A:IfEnd 判断ZDJD的方位角A "TURNING ANGLE -L,+R(Deg)="?Q 输入交点路线转角,左偏为负,右偏为正 If Q<0:Then -1Z[10]:Else 1Z[10]:IfEnd 确定转角数值 QZ[4]:Abs(Q)D 保存转角及符号 "R(m)="?R 输入圆曲线半径 Rtan(0.5D)

T 计算切线长 RDπ÷180

L 计算曲线长 R(cos(0.5D)-1-1)E 计算外距 2T-L

J 计算切曲差

"T(m)=":T 显示切线长

"L(m)=":L 显示曲线长 "E(m)=":E 显示外距 "J(m)=":J 显示切曲差

Z-TZ[1]:Z[1]+0.5LZ[2]:Z[2]+0.5LZ[3] 计算ZY，QZ，YZ点桩号 Z[1]List X[1] 存储ZY点的桩号到统计串列

U-Tcos(A)List Y[1]:V-Tsin(A)List Freq[1] 计算并存储ZY点的测量坐标到统计串列 "ZY PEG(m)=":List X[1] 显示ZY点桩号及坐标 "ZY-X(m)=":List Y[1] "ZY-Y(m)=":List Freq[1]

输入路面左半幅宽度 "WL(m),0 NO="?B

If B>0:Then List Y[1]+Bcos(A-90)Z[11] 计算并显示ZY点的左边桩坐标 List Freq[1]+Bsin(A-90)Z[12] "XL(m)=":Z[11] "YL(m)=":Z[12]

IfEnd

输入路面右半幅宽度 "WR(m),0 NO="?C

If C>0:Then List Y[1]+Ccos(A+90)Z[13] 计算并显示ZY点的右边桩坐标 List Freq[1]+Csin(A+90)Z[14] "XR(m)=":Z[13] "YR(m)=":Z[14]

IfEnd

存储ZY点的桩号到统计串列 Z[2]List X[2]:2

K

Prog "SUB4-31" 调子程序计算并显示QZ点坐标 Z[3]List X[3]:3

K 存储YZ点桩号 Prog "SUB4-31" 调子程序计算并显示YZ点坐标 If R<30:Then 5I:Else If R<60:Then 10I:Else 20

I:IfEnd:IfEnd 确定整桩间距

单圆曲线中边桩坐标计算程序PM4-3

"INT DIST(m)="?O 输入整桩间距 If O>0:Then O

I:IfEnd 输入的整桩间距大于零时，使用输入值 "ONLY CALC +PEG(1)"?G 输入1为只计算加桩点的坐标 If G=1:Then Goto 1:IfEnd

计数变量计数 K+1

K

Int(Z[1]÷I)I+IList X[K] 计算并存储从ZY点开始的第一个整桩号 Prog "SUB4-31" 调子程序计算并显示坐标 Do 计算ZY点至YZ点 K+1

K 计数变量计数 List X[K-1]+IList X[K] 计算并存储整桩号 Prog "SUB4-31" 调子程序计算并显示坐标 LpWhile List X[K]+I<Z[3] 没计算到YZ点时继续循环

Lbl 1 Do

"+PEG(m)="?

F 输入加桩号 If F

List X[1]:Then Break:IfEnd 加桩号小于ZH点桩号时结束程序运行 If F

List X[3]:Then Break:IfEnd 加桩号大于HZ点桩号时结束程序运行 K+1

K 计数变量计数 FList X[K] 存储加桩号到统计串列 Prog "SUB4-31" 调子程序计算并显示偏角弦长

LpWhile F>0 没有计算到YZ点时继续循环 "PM4-3END"

(3) 子程序——SUB4-31，占用内存498字节。

根据存储在List X[K]中的逐桩点桩号，计算中、左、右桩点坐标。 List X[K]-Z[1]Z[5] 计算ZY点至逐桩点的弧长 90Z[5]÷(Rπ)Z[6] 计算ZY点至逐桩点的弦切角 2Rsin(Z[6])Z[7] 计算ZY点至逐桩点的弦长 A+Z[10]Z[6]Z[8] 计算ZY点至逐桩点的方位角 A+2Z[10]Z[6]Z[9] 计算j点切线方位角 List Y[1]+Z[7]cos(Z[8])List Y[K] 存储X坐标到统计串列 List Freq[1]+Z[7]sin(Z[8])List Freq[K] 存储Y坐标到统计串列 If List X[K]=Z[2]:Then "QZ PEG(m)=" 判断桩号提示字符 Else If List X[K]=Z[3]:Then "YZ PEG(m)=" Else "INT PEG(m)=":IfEnd:IfEnd

显示桩号 List X[K]

"X(m)=":List Y[K] "Y(m)=":List Freq[K] "WL(m),0 NO="?B

计算左边桩坐标 If B>0:Then List Y[K]+Bcos(Z[9]-90)Z[11]

List Freq[K]+Bsin(Z[9]-90)Z[12] "XL(m)=":Z[11] "YL(m)=":Z[12]

IfEnd

"WR(m),0 NO="?C

计算右边桩坐标 If C>0:Then List Y[K]+Ccos(Z[9]+90)Z[13]

List Freq[K]+Csin(Z[9]+90)Z[14] "XR(m)=":Z[13] "YR(m)=":Z[14]

IfEnd

单圆曲线中边桩坐标计算程序PM4-3

(4) 程序说明

1) 程序能计算单交点单圆曲线的中边桩坐标，其中三个主点ZY，QZ，YZ的桩号顺序存储在统计串列List X[1]~List X[3]中，其中桩的x坐标顺序存储在统计串列List Y[1]~List Y[3]中，y坐标顺序存储在统计串列List Freq[1]~List Freq[3]中，其后紧接着存储整桩点及加桩点的桩号及其x，y坐标。

2) 每计算并显示完一个中桩点的坐标后，屏幕提示输入该中桩点的左、右幅半宽如下： WL(m),0 NO=? WR(m),0 NO=?

只有当用户输入的左幅半宽>0时，程序才计算并显示左边桩点坐标；用户输入的右幅半宽>0时，程序才计算并显示右边桩点坐标。程序计算出的左、右边桩点坐标只通过屏幕显示，没有存储到统计串列中，所以，运行程序时，请用户及时抄录计算程序结果。

3) 程序不能计算路面加宽的情形，如需要计算路面加宽的渐变值，请先使用程序PM10-3计

=算，求出逐桩点的加宽值后，在执行本程序。 =

=x，y坐标、4) 程序显示完计算标题后，顺序提示用户输入的已知数据分别为：交点桩号及其==转点桩号及其x，y坐标，路线转角(左转输入负数，右转输入正数)与圆曲线半径。完成已知数据=输入后，屏幕先显示曲线要素：切线长、曲线长、外距、切曲差及ZY，QZ，YZ三个主点的桩==号；然后按桩号大小顺序显示主点与逐桩点的桩号及其x，y坐标。

表4-3为左转角坐标计算案例，路面左半幅宽wL=wR=5m。如要计算超高加宽缓和段的加

宽值，需要先执行程序PM10-3计算加宽值。

表4-3 使用程序PM4-3进行单交点单圆曲线左转角中边桩坐标计算案例

y(m) xL(m) 序 桩号 x(m)

1 ZY K62+24.764 7639.376 3185.7317644.3122 QZ K62+117.656 7633.146 3278.2817638.1443 YZ K62+210.549 7644.119 3370.3897649.0084 40 7637.176 3200.8087642.1345 60 7634.988 3220.6867639.9686 80 7633.597 3240.6367638.5917 100 7633.004 3260.6267638.0048 120 7633.212 3280.6247638.2099 140 7634.219 3300.5977639.20610 160 7636.024 3320.5147640.99311 180 7638.623 3340.3437643.567

12 200 7642.014 3360.0527646.92413 70 7634.193 3230.6547639.18114 190 7640.220 3350.2157645.148

=

yL(m)

3186.5283278.1533369.3393201.4543221.1343240.8853260.6753280.4723300.2463319.9643339.594

3359.1063231.0023349.367

yR(m) 7634.4397628.1487639.2317632.2187630.0087628.6037628.0057628.2147629.2317631.0547633.680

7637.1057629.2057635.293

yR(m) 已知数据与要素 3184.934 JDK62+118.743278.410 △=-21°17′22″3371.438 R500m3200.161 T93.976m3220.239 L185.785m3240.388 E8.755m3260.578 J2.167m3280.776 ZY62024.764m3300.949 QZ62117.656m3321.065 YZ62210.549m

xJD=7624.3943341.092

yJD=3278.5063360.999

xZD=7643.2513230.307 yZD=3161.7323351.063

执行程序PM4-3，计算表4-3数据的屏幕提示与用户操作过程如下：

屏幕提示 按键

SINGLE CIRCLE CURVE METHOD OF COORDINATE PM4-3 JD MILEAGE PEG(m)=?

62118.74 JD X(m)=? 7624.394

JD Y(m)=? 3278.506

ZD X(m)=? 7643.251 ZD Y(m)=? 3161.732

17

22TURNING ANGLE -L,+R(Deg)=? -21

R(m)? 500 T(m)=93.976 L(m)=185.785 E(m)=8.755

说明

显示程序标题1

显示程序标题2

输入以m为单位的交点桩号 输入交点的x坐标 输入交点的y坐标 输入转点的x坐标 输入转点的y坐标 输入交点转角 输入圆曲线半径

显示切线长

显示圆曲线长

显示圆曲线外距

单圆曲线中边桩坐标计算程序PM4-3

J(m)=2.167 ZY PEG(m)=62024.764 ZY-X(m)=7639.376 ZY-Y(m)=3185.731

WL(m),0 NO=? 5 XL(m)=7644.312 YL(m)=3186.528

=

WR(m),0 NO=? 5

XR(m)=7634.439 YR(m)=3184.934 QZ PEG(m)=62117.656 X(m)=7633.146 Y(m)=3278.281 WL(m),0 NO=?5 XL(m)=7638.144 YL(m)=3278.153 WR(m),0 NO=?5 XR(m)=7628.148 YR(m)=3278.410 YZ PEG(m)=62210.549 X(m)=7644.119 Y(m)=3370.389 WL(m),0 NO=?5 XL(m)=7649.008 YL(m)=3369.339 WR(m),0 NO=?5 XR(m)=7639.231 YR(m)=3371.438 INT DIST(m)=?

0 ONLY CALC +PEG(1)?

0 INT PEG(m)=62040

X(m)=7637.176 Y(m)=3200.808 WL(m),0 NO=?5 XL(m)=7649.008 YL(m)=3369.339 WR(m),0 NO=?5 XR(m)=7639.231 YR(m)=3371.438

…… …… INT PEG(m)=62200

X(m)=7642.014 Y(m)=3360.052 WL(m),0 NO=?5 XL(m)=7646.924 YL(m)=3359.106 WR(m),0 NO=?5 XR(m)=7637.105 YR(m)=3360.999 +PEG(m)? 62070X(m)=7634.193 Y(m)=3230.654 WL(m),0 NO=?5 XL(m)=7639.181 YL(m)=3231.002 WR(m),0 NO=?5 XR(m)=7629.205

显示切曲差

显示ZY点桩号

显示ZY点的中桩坐标

输入左半幅路宽

显示ZY点的左边桩坐标

输入右半幅路宽

显示

ZY点的右边桩坐标

显示QZ点桩号

显示

QZ点的中桩坐标

按键使用原值

显示

QZ点的左边桩坐标

按键使用原值

显示

QZ点的右边桩坐标

显示YZ点桩号

显示

YZ点的中桩坐标

按键使用原值

显示

YZ点的左边桩坐标

按键使用原值

显示YZ

点的右边桩坐标

输入整桩间距0，程序自动选择 输入0为先计算整桩点

显示ZY点开始的第1个整桩号

显示第

1个整桩号的中桩坐标

按键使用原值

显示第

1个整桩号的左边桩坐标

按键使用原值

显示第1个整桩号的右边桩坐标

……

显示ZY点开始的第9个整桩号

显示第

9个整桩号的中桩坐标

键使用原值

按

显示第

9个整桩号的左边桩坐标

按键使用原值

显示第9个整桩号的右边桩坐标

输入第1个加桩号

显示第

1个加桩号的中桩坐标

按键使用原值

显示第

1个加桩的左边桩坐标

按键使用原值

显示第1个加桩的右边桩坐标

单圆曲线中边桩坐标计算程序PM4-3

YR(m)=3230.307 +PEG(m)? 62190X(m)=7640.220 Y(m)=3350.215 WL(m),0 NO=?5 XL(m)=7645.148 YL(m)=3349.367 WR(m),0 NO=?5 XR(m)=7635.293 YR(m)=3351.063 +PEG(m)? -2

PM4-3END

输入第2个加桩号

显示第2

个加桩号的坐标

按键使用原值

显示第2

个加桩的左边桩坐标

按键使用原值

显示第2个加桩的右边桩坐标

输入负数结束程序 程序运行结束显示

完成计算后，切线长、曲线长、外距、切曲差分别存储在字母变量T，L，E，J中，主点与逐桩点的桩号及其中桩坐标分别存储在List X，List Y，List Freq中，左、右边桩坐标没有保存，可按(REG)键进入双变量统计与回归模式查看中桩坐标。

=

=

本文来源：https://www.bwwdw.com/article/bgn4.html