一种计算各类线路任意线元坐标的方法与流程

1.本发明属于线路工程施工技术领域，，具体是一种计算各类线路任意线元坐标的方法。


背景技术：

2.由于受地形、地质、技术条件的限制和经济发展的需要，各类线路工程施工的方向要不断改变。为了保持线路的圆顺，在改变方向的两相邻直线间须用曲线连接起来，这种曲线称为平面曲线。常见的线路平面曲线分为两种类型：一种是圆曲线与前后直线直接连接，主要用于铁路专用线和行车速度不高的线路上；另一钟是前后带缓和曲线的圆曲线，通过缓和曲线的过渡，使圆曲线与前后直线平顺衔接。上述两种曲线类型属常规曲线，通过切线支距法、偏角法进行曲线坐标的计算已比较成熟，作业精度和效率均可满足施工需求。
3.在匝道、立交桥等工程应用中，为实现线路在短距离内调整方向、位置和高度的目的，常设置由直线、圆曲线、缓和曲线、卵形曲线等各种线元任意组合的不规则平面线元布置，如公式选用不当就会出现较大计算误差，即便是能对切线支距公式进行多项展开，导致该类线路计算放样的精度、效率层次不齐。复化simpson公式是该领域应用较多的一种计算方式，但该公式涉及需将单一线元平分为n等分分别计算，在n的取值上，取值太小会影响线元坐标的计算精度，取值太大则导致计算步骤多，工作量大，数据可靠性和作业效率无法得到有效保证。
4.因此，亟需一种能适应各种线元类型，且数据计算精度和效率都能满足施工需求的计算方法以解决上述至少一种问题。


技术实现要素：

5.本发明为了解决上述问题，提供一种计算各类线路任意线元坐标的方法。
6.本发明采取以下技术方案：一种计算各类线路任意线元坐标的方法，包括以下步骤，s100～设置初始计算参数；s200～输入线元相关参数；s300～输入待求点相关参数；s400～计算待求点切线方位角；s500～分别计算n等分点切线方位角；s600～分别计算2n等分点切线方位角；s700～计算待求点对应中桩坐标；s800～计算待求点坐标；s900～输出计算结果并返回下一个待求点计算。
7.步骤s100包括以下步骤，
8.s101～设置计算数值保留位数；s102～设定角度计算单位；s103～使用拓展变量，并设定拓展变量最大数为30；s104～设定线元方向默认右转。
9.步骤s300包括以下步骤，s301～在待求点相关参数输入前设置定位标记，包括待求点里程、夹角以及偏距，使下个待求点计算时无需重复输入线元相关参数；s302～输入边桩任意夹角和偏距；s303～默认夹角为90度。
10.步骤s500包括以下步骤，s501～设定n＝6；s502～以n
‑
1为循环步数，利用公式分别计算各n等分点切线方位角；
[0011][0012][0013]
α
i
＝α
a
+β1[0014]
式中：ρi—待求点曲率；
[0015]
ρ
a
—曲线元起点曲率；
[0016]
ρ
b
—曲线元终点曲率；
[0017]
k
i
—待求点里程；
[0018]
k
a
—曲线元起点里程；
[0019]
k
b
—曲线元终点里程；
[0020]
β1—切线方位角增量；
[0021]
α
i
—待求点切线方位角；
[0022]
α
a
—曲线元起点切线方位角；
[0023]
s503～通过变量n自动计算循环次数。
[0024]
步骤s600包括以下步骤，s601～设定n＝6；s602～以n为循环步数，利用公式计算各2n等分点切线方位角；
[0025][0026][0027]
α
i
＝α
a
+β1[0028]
s603～通过变量n自动计算循环次数。
[0029]
步骤s700中，计算待求点对应中桩坐标，使用以下公式计算，
[0030][0031][0032]
式中：x
z
—待求点中桩x坐标；
[0033]
y
z
—待求点中桩y坐标；
[0034]
x
a
—曲线元起点x坐标；
[0035]
y
a
—曲线元起点y坐标；
[0036]
α
k
—曲线元n等分点处的切线方位角；
[0037]
—曲线元2n等分点处的切线方位角；
[0038]
h—等分段长度；
[0039]
步骤s800包括以下步骤，判断偏距输入值的正负，为正值时，方位角不变，为负值
时，方位角加180度，计算反方向偏距；s802～使用以下公式将极坐标数据转换为直角坐标增量及待求点坐标；
[0040]
δx＝d.cosα
[0041]
δy＝d.sinα
[0042]
x＝x
z
+δx
[0043]
y＝y
z
+δy。
[0044]
式中：δ
x
—待求点较中桩x坐标增量；
[0045]
δ
y
—待求点较中桩y坐标增量；
[0046]
d—待求点较中桩距离；
[0047]
α—中桩至待求点方位角；
[0048]
x—待求点x坐标；
[0049]
y—待求点y坐标。
[0050]
步骤s900包括以下步骤，s901～切线方位角直接显示六十进制数值；s902～通过矩阵显示排列，将输出结果输出。
[0051]
与现有技术相比，本发明对计算设备的初始参数进行设置后，将线元基本参数赋值于变量，完成线元基本数据库的建立；设置定位标记后，输入待求点对应中桩的里程、夹角及偏距信息，计算待求点对应中桩切线方位角；通过循环计算，依次计算待求点对应中桩至线元起点段落各n等分点的切线方位角及各2n等分点的切线方位角；根据线元基础数据及前述计算参数，使用复化simpson公式计算待求点对应中桩的坐标；将待求点与对应中桩的极坐标关系转换为直角坐标增量，计算待求点坐标，并将待求点坐标及其对应的切线方位角予以输出显示。
[0052]
本发明提供的simpson公式计算任意线元坐标的方法，适用于各种常规或非常规线元布置，适用范围广；通过fx5800语言的编程，减少人工干预步骤，增强了线元坐标计算结果精度，提高了工作效率，解决了非常规曲线布置数据计算精度与效率的矛盾冲突。
附图说明
[0053]
图1为本发明流程图。
具体实施方式
[0054]
以下结合附图对本发明的具体实施方式进行详细说明。应当理解的是，此处所描述的具体实施方式仅用于说明和解释本发明，并不用于限制本发明。
[0055]
需要说明的是，在不冲突的情况下，本发明中的实施例及实施例中的特征可以相互组合。
[0056]
在本发明中，在未作相反说明的情况下，使用的方位词如“上、下、顶、底”通常是针对附图所示的方向而言的或者是针对竖直、垂直或重力方向上而言的各部件相互位置关系描述用词。
[0057]
下面将参考附图并结合实施例来详细说明本发明。
[0058]
在一个实施例中，一种计算各类线路任意线元坐标的方法，包括以下步骤，
[0059]
s100～设置初始计算参数；曲线方向分为左转向与右转向，对应的切线方位角增
量符号会发生相应改变。鉴于此，对代表线元转向的变量赋值为1，当实施例线元方向为右转时，不用单独输入变量。反之线元为左向时，重新输入
‑
1。s101～设置计算数值保留位数；
[0060]
s102～设定角度计算单位；s103～使用拓展变量，并设定拓展变量最大数为30；s104～设定线元方向默认右转。
[0061]
fx5800语言如下：
[0062]
fix 4：deg：设定保留位数，设定角度单位；
[0063]
30
→
dimz：1
→
l：设定扩展变量数量，默认右转曲线。
[0064]
s200～输入线元相关参数；依次输入线元转向参数、线元起点坐标、线元起点方位角、线元起终点里程及半径。同时，将线元起终点半径转化为曲率。
[0065]
fx5800语言如下：
[0066]”r＝+1,l＝
‑
1”？l：”xq＝”？r：”yq＝”？o：输入曲线转向及曲线起点坐标；
[0067]”aq＝”？k：”k1＝”？e：”k2＝”？f：输入曲线起点方位角及起终点里程；
[0068]”r1＝”？g：”r2＝”？h：输入曲线起点和终点的半径；
[0069]
g≠0＝>1/g
→
g：h≠0＝>1/h
→
h：计算曲线起点和终点的曲率。
[0070]
s300～输入待求点相关参数；s301～在待求点相关参数输入前设置定位标记，使下个待求点计算时无需重复输入线元相关参数；s302～输入边桩任意夹角和偏距；s303～默认夹角为90度。
[0071]
待求点位置通常用线路中桩的相对关系来表示。偏距夹角指中桩沿切线前进方向右偏的角度，大多数情况下待求点垂直于对应中桩，因此该变量默认为90度，即默认待求点位于中桩右侧垂直方向，也可根据实际情况输入；偏距为正值时，代表中桩沿设置夹角方向至待求点的距离，为负值时，则代表待求点位于该方向的反向。在工程习惯上，偏距符号通常设置为左负右正，如默认90度夹角不发生改变时，正偏距代表中桩右侧垂直偏距，负偏距代表中桩左侧垂直偏距。
[0072]
设置定位标记，每个待求点计算无需重复输入线元基本参数；每个待求点依次输入对应中线里程、偏距夹角及偏距。
[0073]
fx5800语言如下：
[0074]
lbi 0：待求点定位标记；
[0075]”k＝”？m：90
→
p：”jj＝”？p：”d＝”？q：输入待求点里程、夹角、偏距。
[0076]
s400～计算待求点切线方位角；通过曲率内插、偏角计算等过程，计算待求点对应中桩切线方位角。该步骤为常规计算。
[0077]
fx5800语言如下：
[0078]
l(g+(h
‑
g)/(f
‑
e)*(m
‑
e)+g)(m
‑
e)90/π+k
→
z[25]：计算待求点对应中桩切线方位角。
[0079]
s500～分别计算n等分点切线方位角；
[0080]
s501～设定n＝6；
[0081]
s502～以n
‑
1为循环步数，利用公式分别计算各n等分点切线方位角；
[0082]
[0083][0084]
α
i
＝α
a
+β1[0085]
式中：ρi—待求点曲率；
[0086]
ρ
a
—曲线元起点曲率；
[0087]
ρ
b
—曲线元终点曲率；
[0088]
k
i
—待求点里程；
[0089]
k
a
—曲线起点里程；
[0090]
k
b
—曲线终点里程；
[0091]
β1—切线方位角增量；
[0092]
α
i
—待求点切线方位角；
[0093]
α
a
—曲线起点切线方位角。
[0094]
s503～通过变量n自动计算循环次数。
[0095]
待求点对应中桩所处线元平分的n等分越多，数据计算的精度越高，但同时也意味着计算量的增大，手工计算难以快速统计结果，且极容易出错。通过使用fx5800语言编程，运用do lp.w等循环函数功能，设定相应参数后，由计算设备自动执行相应次数的计算，可解决上述问题。经实践，将n赋值为6时，运用fx5800语言编程数据精度可优于1mm，计算用时4.5秒，可满足施工需求。n的取值可动态调整，满足对计算精度和速度有不同要求的场景应用。设置测段长度变量及计数变量，运用do lp.w循环函数，以n
‑
1为循环次数，分别计算各n等分点的切线方位角，赋值变量为z[2]、z[4]、z[6]、z[8]、z[10]。
[0096]
fx5800语言如下：
[0097]0→
i：0
→
n：等分段距离及计数常数归零；
[0098]
do：(m
‑
e)/6+i
→
i：循环起点，设n＝6,分别计算n等分点长度；
[0099]
n+1
→
n：计数常数增加次数；
[0100]
l(g+(h
‑
g)/(f
‑
e)*i+g)i90/π+k
→
z[2n]：分别计算n等分点切线方位角；
[0101]
lp.w n<5：循环计算5次。
[0102]
s600～分别计算2n等分点切线方位角；
[0103]
s601～设定n＝6；
[0104]
s602～以n为循环步数，利用公式计算各2n等分点切线方位角；
[0105][0106][0107]
α
i
＝α
a
+β1[0108]
s603～通过变量n自动计算循环次数。
[0109]
设置测段长度变量及计数变量，运用do lp.w循环函数，以n为循环次数，分别计算各2n等分点的切线方位角，赋值变量为z[1]、z[3]、z[5]、z[7]、z[9]、z[11]。
[0110]
fx5800语言如下：
[0111]
‑
(m
‑
e)/12
→
i：0
→
n：设置等分起点，计数常数归零；
[0112]
do：(m
‑
e)/6+i
→
i：n+1
→
n：循环起点，分别计算2n等分点长度；
[0113]
l(g+(h
‑
g)/(f
‑
e)*i+g)i90/π+k
→
z[2n
‑
1]分别计算2n等分点切线方位角；
[0114]
lp.w n<6：循环计算6次。
[0115]
s700～计算待求点对应中桩坐标；使用以下公式计算，
[0116][0117][0118]
式中：x
z
—待求点中桩x坐标；
[0119]
y
z
—待求点中桩y坐标；
[0120]
x
a
—曲线元起点x坐标；
[0121]
y
a
—曲线元起点y坐标；
[0122]
α
k
—曲线元n等分点处的切线方位角；
[0123]
—曲线元2n等分点处的切线方位角；
[0124]
h—等分段长度。
[0125]
运用复化simpson公式，引用前述步骤计算得出的各类参数，计算待求点对应中桩相对于线元起点的直角坐标增量，计算待求点对应中桩坐标。
[0126]
fx5800语言如下：
[0127]
r+(m
‑
e)/36(cos(k)+4(cos(z[1])+cos(z[3])+cos(z[5])+cos(z[7])+cos(z[9])+cos(z[11]))+2(cos(z[2])+cos(z[4])+cos(z[6])+cos(z[8])+cos(z[10]))+cos(z[25])
→
z[21]：计算待求点对应中桩x值；
[0128]
o+(m
‑
e)/36(sin(k)+4(sin(z[1])+sin(z[3])+sin(z[5])+sin(z[7])+sin(z[9])+sin(z[11]))+2(sin(z[2])+sin(z[4])+sin(z[6])+sin(z[8])+sin(z[10]))+sin(z[25])
→
z[22]：计算待求点对应中桩y值。
[0129]
s800～计算待求点坐标；
[0130]
s801～判断偏距输入值的正负，为正值时，方位角不变，为负值时，方位角加180度，计算反方向偏距；
[0131]
s802～使用以下公式将极坐标数据转换为直角坐标增量及待求点坐标；
[0132]
δx＝d.cosα
[0133]
δy＝d.sinα
[0134]
x＝x
z
+δx
[0135]
y＝y
z
+δy。
[0136]
式中：δ
x
—待求点较中桩x坐标增量；
[0137]
δ
y
—待求点较中桩y坐标增量；
[0138]
d—待求点较中桩距离；
[0139]
α—中桩至待求点方位角；
[0140]
x—待求点x坐标；
[0141]
y—待求点y坐标。
[0142]
在步骤s700中，运用复化simpson公式，引用步骤s400
‑
s600计算得出的各类参数，计算待求点对应中桩相对于线元起点的直角坐标增量，计算待求点对应中桩坐标。
[0143]
fx5800语言如下：
[0144]
q<0＝>p+180
→
p：偏距小于0时，方位角加180度；
[0145]
rec(abs(q),z[25]+p)：计算待求点与对应中桩的坐标增量；
[0146]
z[21]+i
→
z[23]：z[22]+j
→
z[24]：计算待求点坐标。
[0147]
s900～输出计算结果并返回下一个待求点计算。通过
→
dms函数将待求点对应中桩切线方位角显示为六十进制；将显示数据矩阵排列，将输出结果显示在同一屏上。输出结束后，返回至步骤s300，重复步骤s300
‑
s900，在同一线元范围内循环计算各待求点坐标，无需重复输入线元参数。
[0148]
fx5800语言如下：
[0149]“qx＝”：z[25]
→
dms显示待求点切线方位角；
[0150]
cls：locate 1,1,”xb＝”：locate 1,3,”yb＝”：
[0151]
locate 4,2,z[23]：locate 4,4,z[24]显示待求点坐标；
[0152]
goto 0返回待求点输入界面。
[0153]
本发明提供的fx5800语言复化simpson公式计算任意线元坐标的方法，通过fx5800语言编程将计算步骤进行有机整合。对计算设备的初始参数进行设置后，将线元基本参数赋值于变量，完成线元基本数据库的建立；设置定位标记后，输入待求点对应中桩的里程、夹角及偏距信息，计算待求点对应中桩切线方位角；通过循环计算，依次计算待求点对应中桩至线元起点段落各n等分点的切线方位角及各2n等分点的切线方位角；根据线元基础数据及前述计算参数，使用复化simpson公式计算待求点对应中桩的坐标；将待求点与对应中桩的极坐标关系转换为直角坐标增量，计算待求点坐标，并将待求点坐标及其对应的切线方位角予以输出显示。本发明提供的fx5800语言复化simpson公式计算任意线元坐标的方法，适用于各种常规或非常规曲线布置，适用范围广；通过fx5800语言的编程，减少人工干预步骤，增强了线元坐标计算结果精度，提高了工作效率，解决了非常规曲线布置数据计算精度与效率的矛盾冲突。
[0154]
(1)经济效益
[0155]
在匝道、立交桥等工程的不规则卵形曲线平面布置中，线路坐标的计算步骤繁琐，难以通过手工计算实现，常见的做法是购置的商业软件和其对应的平台进行计算。商业软件通常运行于计算机、智能全站仪、智能手机等平台，软件费用从5000元
‑
20000元不等，平台购置费用3000元
‑
200000元不等，单套平均综合费用约82337元。
[0156]
本发明应用于施工现场普及应用的casiofx5800平台，计算器购置成本低，无需额外增加软件成本，在满足施工现场需求的前提下，单套系统配置可节约软件平台购置成本81837元，增效99％以上。
[0157]
(2)技术经济分析
[0158]
采用复化simpson公式进行一个精度优于1mm的线路坐标点计算，需完成计算式43步。如采用手工计算，熟练情况下计算加检核过程约用时70分钟，基本无法满足施工现场需求；通过fx5800语言对计算步骤进行编译，完成相应计算和和检核约用时10秒，增效99％以
上。
[0159]
(3)社会效益
[0160]
工法适用范围广，推广成本低，对提高工程质量、节约施工成本及加强人才培养具有积极的意义，具有较大的普及性和广泛的推广应用价值。
[0161]
以上结合附图详细描述了本发明的优选实施方式，但是，本发明并不限于上述实施方式中的具体细节，在本发明的技术构思范围内，可以对本发明的技术方案进行多种简单变型，这些简单变型均属于本发明的保护范围。
[0162]
另外需要说明的是，在上述具体实施方式中所描述的各个具体技术特征，在不矛盾的情况下，可以通过任何合适的方式进行组合，为了避免不必要的重复，本发明对各种可能的组合方式不再另行说明。
[0163]
此外，本发明的各种不同的实施方式之间也可以进行任意组合，只要其不违背本发明的思想，其同样应当视为本发明所公开的内容。