平曲线视距横净距地计算Word格式.docx

平曲线视距横净距地计算Word格式.docx

《平曲线视距横净距地计算Word格式.docx》由会员分享，可在线阅读，更多相关《平曲线视距横净距地计算Word格式.docx（10页珍藏版）》请在冰豆网上搜索。

　　对于半径为R1、缓和曲线长度为LS1的平曲线，行车轨迹线上对应的曲线半径R2＝R1－△R，对应的缓和曲线长度LS２应是多少呢?

　　我们分别用p、q、E、T表示圆曲线设置缓和曲线后的内移值、切垂距、外距、切线长，并分别用下标1、2对路中线、行车轨迹线上有关的量进行的区分。

如图2所示，因路中线与行车轨迹线的径向间距为△R，则应有：

图二

　　将有关各量代入并整理后可得:

　　在R1R2LS1是定值的情况下，使△p、△q同时等于零是不可能的（LS1＝0时除外）。

也就是说，无论LS２取何值，由R2LS２确定的轨迹线（下文中称之为计算行车轨迹线）与理论行车轨迹线（距未加宽时的行车道内侧边缘1.5m）不可能完全重合，两轨迹线间的偏离程度可表示为：

将有关各量代入并整理后可得：

上式中，△E表示计算行车轨迹线与理论行车轨迹线的曲线中点间的距离，△T表示两轨迹线上的曲线起（终）点沿切线方向的距离。

　　既然计算行车轨迹线与理论行车轨迹线不能完全重合，那么我们能否选取适当的LS２值，使两轨迹线尽量接近，从而满足横净距计算的精度要求呢?

由

（1）、

（2）式可以看出，△p对两轨迹线间偏离程度的影响远比△q的影响大，如果我们通过LS２的取值使△p等于零，即用牛顿迭代法按下式解算LS２

此外，两迹线的偏离程度可表示为：

　　△E＝0△T＝△q

　两轨迹线在圆曲线段重合，下面我们来看看两轨迹线在缓和曲线段的偏离情况。

　如图3所示，A为路中线上缓和曲线段内的任一点，距缓和曲线起点ZH1的长度为L1，在xoy坐标系中的坐标为（XA，YA），A点的法线AF与x轴的夹角为：

沿法线AF由A点向曲线内侧移动△R得B点，则B点的坐标为：

由路中线与理论行车轨迹线的关系知：

B点在理论行车轨迹线上。

　　令C点为计算行车轨迹线上缓和曲线段内一点，距缓和曲线段起点ZH２的长度为L２，在x'

o'

y'

坐标系中的坐标为（x'

c，y'

c），经坐标轴平移后可得C点在xoy坐标系中的坐标：

　　若Xc＝Xb，则（Yc－Yb）cos（γ－π/2）即可视为该位置处计算行车轨迹线与理论行车轨迹线间的径向偏离值。

　　按上述方法，我们用计算机对《公路曲线测设用表》（第二册）中“缓和曲线长度总表”列出的R1、Ls1的所有组合进行了有关计算。

计算结果表明：

计算行车轨迹线与理论行车轨迹线在缓和曲线段的径向偏离值有以下特点：

（1）两轨迹线在缓和曲线段中部的偏离值最大，向两端逐渐趋于零；

（2）Ls1值不变时，R1值愈小，两轨迹线在缓和曲线段的偏离值就愈大；

（3）R1值不变时，Ls1愈大，两轨迹线在缓和曲线段的偏离值就愈大；

（4）表列各种组合中，两轨迹线在缓和曲线段的最大偏离值一般均小于0.01m，最大偏离值超过0.01m的组合只有以下4个：

①R1＝65、Ls1＝60、b＝7.5时，最大偏离值为0.015m；

②R1＝60、Ls1＝60、b＝7.5时，最大偏离值为0.018m；

③R1＝35、Ls1＝35、b＝6时，最大偏离值为0.012m；

④R1＝30、Ls1＝35、b＝6时，最大偏离值为0.016m。

通过对两轨迹线间偏离值的计算可知：

由R２＝R1－△R和按（3）式解算的Ls２确定的计算行车轨迹线与理论行车轨迹线非常接近（Ls1＝0时两轨迹线重合），若用计算行车轨迹线代替理论行车轨迹线来计算横净距，最大误差不会超过0.02m，完全能够满足横净距计算的精度要求。

2.3横净距的计算

为方便计算，我们也用桩号表示行车轨迹线上各点间的相对位置关系，并假定行车轨迹线上曲线起点ZH2（ZY2）与路中线上曲线起点ZH1（ZY1）的桩号相同（两线上的桩号是相互独立的）。

这样对于任一给定的桩号，能且只能在行车轨迹线上找到一点与之对应。

如图1，欲求横净距的断面M的桩号为L0，我们不妨先假设行车轨迹线上视距起终点A、B的桩号分别为L0－S／2、L0＋S／2（在应用范围内，A、B不会落在法线MN的同一则），此时沿行车轨迹线A、B间的长度等于设计视距S，计算AB与MN的交点到M点的距离值H；

然后将A、B两点沿行车轨迹线向前移动0.01m，即A、B点的桩号同时增大0.01m，再计算MN与新直线AB的交点到M点的距离值H。

根据两次计算的距离值H的大小，使A、B向距离值H增大的方向同步等量逐渐移动（以保证沿行车轨迹线A、B间的长度恒等于设计视距S）直到距离值H开始减小为止，此时的距离值H与△R之差即为M点的横净距。

按上述方法，我们编写了计算横净距的程序。

3实例应用

鹤大公路鸡牡段林口支线为山岭区三级公路，行车道宽6m，设计视距（二倍停车视距）S＝60m；

JD2转角α＝76°

32′，曲线半径100m，缓和曲线长Ls＝25m，曲线起点桩号为K2＋838.68，计算各有关断面的横净距。

程序运行时的操作及计算结果如下（带下画线者为输入的数据）：

曲线起点桩号?

838.68

转角：

?

76°

32′?

0〃

平曲线半径?

100

缓和曲线长?

25

行车道宽度?

6

设计视距?

60

曲线主点桩号：

ZH＝838.68

HY＝863.68

QZ＝917.97

YH＝972.26

HZ＝997.26

断面桩号（＜0结束）?

800横净距：

0.00

820横净距：

0.34

838.68横净距：

1.41

840横净距：

1.53

860横净距：

3.52

863.68横净距：

3.80

880横净距：

4.45

900横净距：

4.53

917.97横净距：

920横净距：

940横净距：

960横净距：

4.37

980横净距：

3.14

997.26横净距：

1000横净距：

1.19

1020横净距：

0.23

1040横净距：

－1

OK

4源程序

10CLEAR：

DEFDBLA－H，K－Z：

PI＝3.141592653589793#

20CLS：

LOCATE11，1：

INPUT〃曲线起点桩号：

〃；

ZH1

30PRINT〃转角：

40LOCATE12，6：

INPUTZJ1：

LOCATE12，14：

INPUTZJ2：

LOCATE12，22：

INPUTZJ3

50INPUT〃平曲线半径＝〃；

R1：

INPUT〃缓和曲线长＝〃；

LS1

60INPUT〃平曲线半径＝〃；

LMK：

INPUT〃设计视距＝〃；

S

70ZJ＝（ZJ1＋ZJ2／60＋ZJ3／3600）*PI／180：

DR＝LMK／2－1.5：

R2＝R1－DR

80IFLS1＞0THENGOSUB640：

REM计算LS2

90R＝R1：

LS＝LS1：

GOSUB710：

LY1＝LY：

P1＝P：

Q1＝Q：

T1＝T

100R＝R2：

LS＝LS2：

LY2＝LY：

P2＝P：

Q2＝Q：

T2＝T

110DT＝Q1－Q2：

ZH2＝ZH1：

SS＝S／2／R2*R1

120HY1＝ZH1＋LS1：

YH1＝HY1＋LY1：

HZ1＝YH1＋LS1：

QZ1＝（ZH1＋HZ1）／2

130HY2＝ZH2＋LS2：

YH1＝HY2＋LY2：

HZ2＝YH2＋LS2：

QZ2＝（ZH2＋HZ2）／2

140PRINT〃曲线主点桩号：

〃：

A＄＝〃######.##〃：

IFLS1＝0THEN180

150PRINT〃ZH＝〃；

USINGA＄；

ZH1：

PRINT〃HY＝〃；

HY1

160PRINT〃QZ＝〃；

QZ1：

PRINT〃YH＝〃；

YH1

170PRINT〃HZ＝〃；

HZ1：

GOTO200

180PRINT〃ZY＝〃；

PRINT〃QZ＝〃；

QZ1

190PRINT〃YZ＝〃；

UXINGA＄；

HZ1

200INPUT〃断面桩号（＜0结束）〃；

L0：

IFL0＜0THEN450

210IFL0＜＝ZH1－SORL0＞＝HZ1＋STHENH＝0：

GOTO440

220IFL0＞QZ1THENL0＝ZH1＋（HZ1－L0）：

REM若欲求横净距断面在曲线后半部，则计算曲线前半部对称位置的横净距

230IFL0－SS＞＝HY1ANDL0＋SS＜＝YH1THENH＝R2*（1－COS（S／2／R2））：

240R＝R1：

P＝P1：

Q＝Q1：

X0＝T1＋T1*COS（ZJ）：

Y0＝T1*SIN（ZJ）：

IFABS（ZJ－PI）＜1E－10THENX0＝0Y0＝2*（R1＋P1）

250IFL0＜＝ZH1THENA0＝PI／2：

X＝L0－ZH1：

Y＝0：

GOTO300

260IFL0＞＝HZ1THENA0＝ZJ＋PI／2：

X＝HZ1－L0：

GOSUB620：

270IFL0＞YH1THENL＝HZ1－L0：

A0＝ZJ＋PI／2－L^2／2／R1／LS1：

GOSUB600：

280IFL0＞＝HY1THENL＝L0－HY1＋LS1／2：

A0＝L／R1＋PI／2：

GOSUB590：

X＝X＋Q：

Y＝Y＋P：

290L＝L0－ZH1：

A0＝L^2／2／R1／LS1＋PI／2：

GOSUB600

300XM＝X－DT：

YM＝Y－DR

310IFA0＝PI／2THENA1＝1：

B1＝0：

C1＝－XMELSEK＝TAN（A0）：

A1＝K：

B1＝－1：

C1＝YM－K*XM

320R＝R2：

P＝P2：

Q＝Q2：

X0＝T2＋T2*COS（ZJ）：

Y0＝T2*SIN（ZJ）：

IFABS（ZJ－PI）＜1E－10THENX0＝0：

Y0＝2*（R2＋P2）

330A＝L0－S／2：

B＝L0＋S／2：

GOSUB460：

H1＝H：

A＝A＋.01：

B＝B＋.01：

A＝A－.01：

B＝B－.01

340REM分别以10m、1m、0.1m、0.01m的步距使A、B两点向H值增大的方向移动，以减少计算次数

350D＝100：

IFH＜H1THENF＝－1ELSEF＝1

360IFH＝H1ANDB＜＝ZH2THENF＝1

370IFH＝H1ANDA＞HZ2THENF＝－1

380FORI＝1TO3：

D＝D／10：

D1＝D*F

390FORJ＝1TO100：

A＝A＋D1：

B＝B＋D1：

GOSUB460

400IFH＜＝H1THENA＝A－D1：

B＝B－D1：

GOTO420

410H1＝H：

NEXTJ

420NEXTI

430H＝H1－DR

440U＝CSRLIN：

LOCATEU－1，48：

PRINT〃横净距：

USING〃###.##〃；

H：

GOTO200450END

460REM计算视线AB与法线MN的交点E至M点的距离H值

470G＝A：

GOSUB520：

XA＝X：

YA＝Y

480G＝B：

BOSUB520：

XB＝X：

YB＝Y

490A2＝YB－YA：

B2＝XA－XB：

C2＝（XB－XA）*YA－（YB－YA）*XA

500AB＝A1*B2－A2*A2*B2：

XE＝（B1*C2－B2*C1）／AB：

YE＝－（A1*C2－A2*C1）／AB

510H＝SQR（（XE－XM）^2＋（YE－YM）^2）：

RETURN

520REM计算行车轨迹线已知桩号的任一点的坐标

530IFG＜＝ZH2THENX＝G－ZH2：

540IFG＞＝HZ2THENX＝HZ2－G：

550IFLS2＝0THENL＝G－ZH2：

560IFG＞YH2THENL＝HZ2－G：

570IFG＞HY2THENL＝G－HY2＋LS2／2：

580L＝G－ZH2：

590X＝R*SIN（L／R）：

Y＝R*（1－COS（L／R））：

600C＝R*LS：

X＝L－L^5／40／C^2＋L^9／3456／C^4

610Y＝L^3／6／C－L^7／336／C^3＋L^11／42240!

／C^5：

620XI＝－X*COS（ZJ）－Y*SIN（ZJ）：

YI＝Y*COS（ZJ）－X*SIN（ZJ）

630X＝XI＋X0：

Y＝Y1＋Y0：

TETURN

640REM用牛顿迭代法计算LS2

650LS2＝LS1

660FORI＝1TO10

670FX＝LS2^2／24／R2－LS2^4／2688／R2^3＋LS2^6／506880!

／R2^5－LS1^2／24／R1＋LS1^4／2688／R1^3－LS1^6／506880!

／R1^5

680DX＝LS2／12／R2－LS2^3／672／R2^3＋LS2^5／84480!

／R2^5

690LS2＝LS2－FX／DX：

PRINTLS2：

NEXTI

700RETURN

710BT＝LS／2／R：

LY＝R*ZJ－LS：

IFLS＝0THEN730

720L＝LS：

Q＝X－R*SIN（BT）：

P＝Y－R*（1－COS（BT））

730IFZJ＝PITHENRETURN

740T＝（R＋P）*TAN（ZJ／2）＋Q：