Influence of Nose Rail Reconstruction Profile Parameters Based on B-Spline Curve
-
摘要:
针对固定辙叉钢轨廓形不规则变化特点,基于B样条理论提出一种考虑固定辙叉钢轨断面廓形特征的参数化设计方法,并提出拟合钢轨廓形贴合度评价指标和关键控制参数对廓形贴合度影响权重指标;以60 kg/m钢轨12号固定辙叉为例,采用DOE (design of experiments)方法分析轨顶横坡、侧面纵坡和复合圆弧段比例系数对心轨轨头整体廓形、轨顶面廓形和侧面廓形的影响. 结果表明:1) 在心轨20 mm断面,轨顶横坡、侧面纵坡和复合圆弧段比例系数对心轨整体廓形变化的影响权重分别为21.08%、56.89%和22.02%,对心轨顶面廓形变化影响权重分别为8.42%、61.95%、29.63%;在心轨50 mm断面,各关键控制参数对心轨轨头整体廓形变化影响权重分别为55.9%、33.38%、10.72%;侧面纵坡对心轨20 mm和50 mm断面的侧面廓形变化影响权重分别为76.82%、66.04%. 2) 当心轨轨头宽度为20 mm时,侧面纵坡对各部分廓形变化影响均达到50%以上;随着心轨轨头宽度的增大,心轨轨顶横坡对整体廓形变化影响权重由21.8%增至55.9%,侧面纵坡和复合圆弧段比例系数影响权重则分别减小41.3%和51.3%.
Abstract:According to irregular changes in rail profiles in a fixed-nose crossing, a parametric design method was proposed based on the B-spline theory by considering the section characteristics of rail profiles in a fixed-nose crossing. In addition, the rail profile fitting evaluation indexes and the influence weights of key control parameters on the rail profile fitting were proposed. With No.12 turnout of fixed-nose crossing of 60 kg/m rail as an example, the design of experiment (DOE) was employed to analyze the influence of railhead slope, rail side slope, and proportion coefficient of composite circular arc on the changes in the full railhead profile, rail top profile, and rail side profile of the nose rail. The results show that 1) for nose rail with a section of 20 mm, the influence weights of railhead slope, rail side slope, and proportion coefficient of composite circular arc on the full railhead profile change are 21.08%, 56.89%, and 22.02%, respectively, and the influence weights on the rail top profile change of the nose rail are 8.42%, 61.95%, and 29.63%, respectively. For nose rail with a section of 50 mm, the influence weights of the key control parameters on the full railhead profile change are 55.9%, 33.38%, and 10.72%, respectively. For nose rails with a section of 20 mm and 50 mm, the influence weights of the rail side slope on the rail side profile change are 76.82% and 66.04%, respectively. 2) When the nose railhead width is 20 mm, the influence of the rail side slope on the nose profiles is more than 50%. As the nose railhead width increases, the influence weight of railhead slope on full railhead profile gradually increases from 21.8% to 55.9%, while that of rail side profile and proportion coefficient of composite circular arc is reduced by 41.3% and 51.3%, respectively.
-
Key words:
- turnout /
- fixed-nose crossing /
- B-spline curve /
- rail profile /
- parametric design /
- design of experiments
-
固定辙叉是道岔重要组成部分,广泛应用于我国重载和提速铁路. 固定辙叉区轮轨动力相互作用异常剧烈,关键钢轨部件磨耗磨损等病害较区间线路更加频发[1-6]. 我国固定辙叉平均使用寿命仅为美国的40%,欧洲的60%,造成其服役寿命相对较低的主要原因之一是轮载过渡区辙叉钢轨总承载面积相对欧美较小[7].
辙叉区不规则的钢轨廓形是影响轮载过渡区钢轨承载面积的直接因素. 为此,国内外学者从钢轨廓形参数优化角度进行了大量研究:王树国等[7]从理论分析、动力仿真和现场试验角度对辙叉心轨顶宽、查照间隔和护轨轮缘槽宽度等参数进行优化;徐井芒等[8-9]提出基于轮轨廓形净差值的心轨廓形优化方法;曹洋等[10]针对固定辙叉轮载过渡段轮轨匹配性能对心轨关键断面廓形进行优化;沈钢等[11]采用粒子群算法求解优化模型对心轨顶面降低值参数进行优化;张鹏飞等[12]研究翼轨加高值对列车过岔动力特性的影响;Wan等[13-14]采用数值模拟和现场试验分析辙叉区各轮轨动力响应的权重大小,并结合多目标优化方法提出固定辙叉心轨廓形设计方法.
上述文献研究主要针对辙叉廓形单一几何参数,通过调整局部廓形进行优化设计,并未从固定辙叉几何特征控制参数角度探寻廓形参数化设计方法并分析各参数对辙叉钢轨廓形的影响. 因此,本文基于B样条理论,采用DOE (design of experiments)方法,分析各关键控制参数对心轨廓形变化的影响权重,对实现固定辙叉区钢轨廓形参数化设计具有重要意义.
1. 固定辙叉廓形特征参数
在固定辙叉中,有害空间引起轨距线不连续,且翼轨和心轨廓形呈连续不规则变化. 然而,固定辙叉平面布置、心轨各关键断面廓形以及翼轨断面几何特征具有如下规律:对于任意心轨和翼轨断面,轨顶横坡和轨头侧面纵坡不变;沿辙叉走行方向,心轨和翼轨高度分段线性变化,心轨断面宽度线性增大;心轨和翼轨轮缘槽宽度一定;距轨距测量点标准轨头半宽35.5 mm处的翼轨顶面高度以及翼轨宽度可精确确定. 以60 kg/m钢轨12号固定辙叉[15]为例,心轨轨顶横坡为1∶20,轨头侧面纵坡为1∶5;翼轨工作边纵坡为1∶5,轨顶横坡为1∶20,非工作边侧面纵坡1∶20,如图1所示,W为翼轨轨头宽度.
由此可见,虽然辙叉钢轨轨头宽度和高度沿辙叉方向不断变化,但仍存在保持不变的设计参数,分别为各断面处心轨和翼轨轨头宽度和高度、轨顶横坡、轨头侧面坡度. 本文将以上参数统称为固定辙叉钢轨全断面廓形特征参数,为后文中不断重构更新辙叉钢轨廓形提供了设计依据.
2. 辙叉钢轨廓形参数化拟合方法
2.1 B样条曲线理论
辙叉区钢轨各断面廓形是由多段不同曲率圆弧组成的复杂曲线,整个廓形需保持连续平滑. B样条曲线在计算机辅助设计中具有广泛应用[16],可通过较少的轮廓关键控制顶点实现达到期望逼近精度的廓形拟合,且具有连续平滑性和强凸包性. 因此,通过B样条曲线拟合辙叉钢轨廓形,对利用辙叉廓形特征参数进行廓形参数化设计具有明显优势. 已有相关研究将其应用于钢轨廓形设计:Yang等[17]采用B样条理论,选取实测廓形33个控制点建立磨耗钢轨廓形曲线,对不同打磨方式钢轨廓形进行磨耗预测;温士明[18]运用B样条曲线对钢轨廓形进行参数化描述,提出钢轨打磨目标廓形优化方法;刘建桥[19]结合遗传算法和B样条曲线,以降低钢轨磨耗为目标确定钢轨优化廓形;唐彦玲等[20]通过B样条曲线对钢轨廓形进行参数化描述,结合遗传算法对曲线地段钢轨廓形进行设计;王亮等[21]采用B样条曲线和NSGA-Ⅱ算法对控制点权因子寻优计算,对钢轨打磨区域轨头廓形进行设计;林凤涛等[22]采用B样条理论对高速铁路辙叉钢轨廓形拟合,基于标准辙叉廓形设置16个控制点设计了钢轨经济性打磨廓形;史振帅[23]以道岔转辙器区域直尖轨为研究对象,采用NURBS理论对道岔直尖轨打磨廓形进行设计重构. 已有文献虽实现了钢轨廓形重构,但仍是通过标准或实测钢轨廓形的稠密离散点数据进行廓形拟合,且大多针对区间线路和可动心辙叉,对固定辙叉钢轨廓形研究较少.
在B样条方法中,任意曲线可用式(1)表示.
C(t)=n−1∑v=0Bv,d(t)Qv, (1) 式中:Qv为序号v的廓形关键控制点;n为控制点Qv个数,依次连接控制点可以构成一个特征多边形;t为节点,t0≤t≤tm-1,T=(t0,t1,…,tv,…tm-1),T为节点向量,m=n + d + 1,m为节点个数,d为B样条曲线的阶数;Bv,d(t)为由递归方程表示的d−1次B样条基函数,如式(2)所示.
{Bv,d(t)=t−tvtv+d−tvBv,d−1(t)+tv+d+1−ttv+d+1−tv+1Bv+1,d−1(t),Bv,0={1,tv⩽t⩽tv+1,0,t<tv或t>tv+1. (2) 根据B样条曲线理论与相关工程应用情况可知,采用3次样条可保证拟合曲线2阶导数连续,满足工程实际需要,因此,本文采用3次非均匀B样条曲线拟合辙叉钢轨廓形.
2.2 关键控制点确定
依据心轨和翼轨廓形关键几何特征参数可确定B样条曲线拟合廓形所需关键控制点,由各控制点组成心轨和翼轨廓形拟合控制多边形,如图2所示. 由于心轨断面廓形轴对称,因此只选取心轨中心轴一侧进行廓形参数化说明.
在图2(a)中,心轨廓形拟合的关键控制点为{N1,N2…,Ni},Ni=(xi,yi),i =1,2,…,8. 其中:N1为心轨轨顶高度控制点;N2、N3用来调节心轨顶面平直段廓形范围;N4、N6、N8用来调节心轨轨头复合圆弧段廓形;N5为轨顶横坡和侧面纵坡交点;N7为心轨宽度控制点;P为△N4N5N6中垂线与N4N6的交点,为确定控制点N8的辅助点,N8位于N5P的连线上.
在图2(b)中,翼轨廓形拟合的关键控制点为{W1,W2,…,Wj},j=1,2,…,9. 其中:W1和W9为翼轨轨头宽度控制点;W3为轨顶横坡和翼轨工作边纵坡交点;W5为翼轨高度控制点;W7为轨顶横坡和翼轨非工作边纵坡交点;W2、W4、W6、W8用来调节翼轨轨头两侧复合圆弧段廓形;kW1、kW2、kW3分别为翼轨工作边纵坡、轨顶横坡和非工作边纵坡斜率;hW为翼轨高度;wW为标准轨头宽度的一半,即35.5 mm.
由于心轨和翼轨廓形参数化方法一致,且心轨廓形变化较翼轨复杂,因此,以心轨为例详细说明辙叉钢轨廓形参数化过程. 采用向量法确定各控制点间关系,如式(3).
\left\{\begin{array}{l} \dfrac{{{\boldsymbol{N}}}_{4}{{\boldsymbol{N}}}_{5}}{\left|{{\boldsymbol{N}}}_{4}{{\boldsymbol{N}}}_{5}\right|}=(1\text{,}{k}_{1}),\\ \dfrac{{{\boldsymbol{N}}}_{5}{{\boldsymbol{N}}}_{6}}{\left|{{\boldsymbol{N}}}_{5}{{\boldsymbol{N}}}_{6}\right|}=(1\text{,}{k}_{2}),\\ {{\boldsymbol{N}}}_{4}{{\boldsymbol{N}}}_{5}={\lambda }_{1}{{\boldsymbol{N}}}_{1}{{\boldsymbol{N}}}_{5},\\ {{\boldsymbol{N}}}_{1}{{\boldsymbol{N}}}_{2}={\lambda }_{2}{{\boldsymbol{N}}}_{1}{{\boldsymbol{N}}}_{4},\\ {{\boldsymbol{N}}}_{5}{{\boldsymbol{N}}}_{6}={\lambda }_{3}{{\boldsymbol{N}}}_{5}{{\boldsymbol{N}}}_{7},\\ {{\boldsymbol{N}}}_{2}{{\boldsymbol{N}}}_{3}={\lambda }_{4}{{\boldsymbol{N}}}_{2}{{\boldsymbol{N}}}_{4},\\ {{\boldsymbol{N}}}_{5}{{\boldsymbol{N}}}_{8}={\lambda }_{5}{{\boldsymbol{N}}}_{5}{\boldsymbol{P}}, \end{array}\right. (3) 式中:λ1、λ2、λ3、λ4、λ5分别为调节控制点N4、N2、N6、N3、N8位置的比例系数,为可变参数,记为λs∈[0,1],其中,s=1,2,…,5,λ5用来调节心轨轨头复合圆弧段廓形,简称心轨复合圆弧段比例系数;k1、k2分别为心轨轨顶横坡和侧面纵坡斜率.
当其他控制点一定时,心轨轨头复合圆弧段廓形主要由控制点N8位置决定,为了保持心轨轨头廓形凸曲线的特点,N8取值需在N5P范围内. 依据向量法求三角形垂点P,如式(4)所示.
{\boldsymbol{P}} = {{\boldsymbol{N}}_4} + \frac{{{{\boldsymbol{N}}_4}{{\boldsymbol{N}}_5} {\text{•}} {{\boldsymbol{N}}_4}{{\boldsymbol{N}}_6}}}{{\left| {{{\boldsymbol{N}}_4}{{\boldsymbol{N}}_6}} \right| {\text{•}} \left| {{{\boldsymbol{N}}_4}{{\boldsymbol{N}}_6}} \right|}} {\text{•}} ({{\boldsymbol{N}}_6} - {{\boldsymbol{N}}_4}). (4) 结合式(3)、(4),根据定比分点原理确定各关键控制点坐标,如式(5)所示.
\left\{ \begin{array}{l}\left({x}_{1},{y}_{1}\right)=(0,h),\\ \left({x}_{2},{y}_{2}\right)=(1-{\lambda }_{2})({x}_{1},{y}_{1}) + {\lambda }_{2}({x}_{4},{y}_{4}) + (0\text{,}\varDelta ),\\ \left({x}_{3},{y}_{3}\right)=(1-{\lambda }_{4})({x}_{2},{y}_{2}) + {\lambda }_{4}({x}_{4},{y}_{4}),\\ \left({x}_{4},{y}_{4}\right)={\lambda }_{1}({x}_{1},{y}_{1}) + (1-{\lambda }_{1})({x}_{5},{y}_{5}),\\ \left({x}_{5},{y}_{5}\right)=\left(\dfrac{h + {k}_{2}w}{{k}_{2}-{k}_{1}},h + \dfrac{{k}_{1}(h + {k}_{2}w)}{{k}_{2}-{k}_{1}}\right),\\ \left({x}_{6},{y}_{6}\right)=(1-{\lambda }_{3})({x}_{5},{y}_{5}) + {\lambda }_{3}({x}_{7},{y}_{7}),\\ \left({x}_{7},{y}_{7}\right)=(w,0),\\ \left({x}_{8},{y}_{8}\right)=(1-{\lambda }_{5})({x}_{5},{y}_{5}) + {\lambda }_{5}({x}_{\text{p}},{y}_{\text{p}}), \end{array} \right. (5) 式中: (xp,yp)为点P坐标;Δ为心轨平直段廓形修正系数,通过调整Δ使控制点N2与N1在同一水平线上;h为x处心轨高度,x为心轨任意断面距辙叉理论尖端的距离;w为x处心轨轨头半宽.
由式(3)~(5)可知,心轨任意断面廓形各控制点坐标可表示为各特征参数和比例系数的函数. 翼轨同理,则固定辙叉全断面廓形控制点可表示为
\left\{ \begin{gathered} {{{\boldsymbol{N}}_i}} = {f_{{\text{cp}}}}(h,w,{k_1},{k_2},\varDelta ,[{\lambda _s}]), \\ {{{\boldsymbol{W}}_j}} = {g_{{\text{wp}}}}({h_{\text{W}}},{w_{\text{W}}},{k_{\text{W}}}_1,{k_{\text{W}}}_2,[{\lambda _t}]), \\ \end{gathered} \right. (6) 式中:fcp(·)、gwp(·)分别为求解心轨和翼轨任意断面所有控制点的函数;λt为翼轨控制点调节比例系数,记为λt∈[0,1],其中,t为其对应的比例系数编号.
综上所述,关键控制点确定所涉及的参数具体包括两类,一类是与廓形相关的有实际物理意义的特征参数,另一类是B样条曲线几何控制点调节比例系数. 为方便后续分析,将上述两类参数统称为辙叉钢轨廓形关键控制参数.
2.3 廓形拟合
以60 kg/m钢轨12号固定辙叉钢轨为例,采用B样条曲线进行廓形构建,心轨和翼轨控制点分布及拟合廓形分别如图3所示.
由图3可知:各控制点位置随心轨轨顶高度和宽度的变化而改变,从而实现心轨廓形连续变化;本文不考虑翼轨高度和轨头宽度的变化,因此翼轨拟合廓形唯一.
在此基础上,结合辙叉区心轨和翼轨走行方向变化,可进一步确定辙叉三维廓形,考虑一侧翼轨的固定辙叉三维廓形如图4所示.
由图4可知,虽然固定辙叉钢轨廓形复杂多变,心轨和翼轨沿辙叉走行方向廓形不断变化,但结合辙叉结构设计中各关键断面的关键控制参数,基于B样条方法可快速确定沿辙叉走行方向的任意截面廓形,该方法可极大拓展以固定辙叉廓形优化设计为基础的轮轨动力仿真以及轮轨接触关系等的研究.
3. 心轨廓形参数试验设计分析方法(DOE)分析
3.1 评价指标
心轨廓形关键控制参数调整会改变心轨廓形,廓形面积差可以直观量化地表达心轨整体廓形的变化情况;与此同时,不同参数对心轨轨头廓形和心轨侧面廓形变化影响程度不同,而心轨轨头廓形变化量和侧面廓形变化量可以更细致地反映各参数对心轨轨头和侧面廓形的单独影响. 因此,为分析各关键控制参数对心轨廓形的影响,选取ΔS以及D1、D2评价各方案廓形与标准廓形的贴合度,如式(7)所示. 本文借鉴《铁路线路修理规则》中对钢轨垂磨和侧磨的定义,选取心轨轨头顶面和侧面某一位置处廓形距离偏差量来定义D1和D2.
\left\{ \begin{gathered} \Delta S = \left| {{S_{{\text{case}}}} - {S_{{\text{std}}}}} \right|, \\ {D_1} = \left| {{Y_{{\text{case}}}} - {Y_{{\text{std}}}}} \right|, \\ {D_2} = \left| {{X_{{\text{case}}}} - {X_{{\text{std}}}}} \right|, \\ \end{gathered} \right. (7) 式中:Scase为各方案廓形曲线与X坐标轴围成的面积;Sstd为采取表1参数确定的标准廓形与X坐标轴围成的面积;ΔS为各方案廓形与标准廓形面积差,反映心轨整体廓形的变化程度;Ycase和Ystd分别为距心轨中心轴2w/3处,各方案廓形曲线与标准廓形曲线纵坐标值;Xcase和Xstd分别为距心轨顶面h/2处各方案廓形曲线与标准廓形曲线横坐标值;D1和D2分别为各方案廓形与标准廓形在垂向和横向上的廓形偏差.
表 1 DOE全阶乘方案设计Table 1. DOE full factorial scheme design方案 k1 k2 λ5 方案 1 1/30 1/6 1/4 方案 2 1/30 1/6 1/2 方案 3 1/30 1/4 1/4 方案 4 1/30 1/4 1/2 方案 5 1/10 1/6 1/4 方案 6 1/10 1/6 1/2 方案 7 1/10 1/4 1/4 方案 8 1/10 1/4 1/2 标准 1/20 1/5 1/3 在DOE分析中,为明确各关键控制参数在所选合理取值范围内对心轨廓形变化的影响权重E,除了对D1、D2及ΔS各自的均值进行分析,还需对E进行分析,如式(8)所示.
E = \frac{{\dfrac{{I^+ - I^ - }}{{V_b^ + - V_b^{{ - }}}}}}{{\displaystyle\sum\limits_{b = 1}^3 {\dfrac{{I^ + - I^ - }}{{V_b^ + - V_b^{{ - }}}}} }} \times 100{\text{%}}, (8) 式中:$I^ -$和$I^ +$分别为各控制参数分别取2个水平(−和 + )时的各贴合度指标均值,$ V_b^ + $和$ V_b^ - $分别为各关键参数的2个水平值(−和 + ).
3.2 方案设置
在工程实践中,心轨轨头顶面和轨头侧面复合圆弧部分更易发生严重磨耗[24],k1和k2分别对心轨顶面廓形和轨头侧面廓形起控制作用,比例系数λ5能够调整控制点N8的位置,从而实现轨头侧面复合圆弧段廓形的调整. 因此,采用DOE全阶乘设计方法对心轨廓形关键控制参数k1、k2及λ5进行方案设计. 依据《道岔设计手册》60 kg/m钢轨12号道岔[15]中的心轨标准设计参数,在一定范围内对各参数调整取值进行方案设计,以满足心轨整体廓形符合工程实际. 各方案参数取值如表1.
3.3 重构廓形轮轨匹配特征
轮轨接触点分布规律是轮轨廓形匹配的最直观反映[25]. 选取LMA型车轮踏面,以心轨20 mm断面和50 mm断面为例,分析轮对不同横移量下各方案轮轨接触点分布规律,进一步验证该参数化设计方法的合理性. 各方案心轨拟合廓形如图5所示,各方案轮轨接触点分布情况如图6所示.
由图5可知:各方案心轨廓形与标准廓形贴合度良好,说明取得了较好的钢轨廓形重构效果. 所建立的钢轨廓形B样条曲线构造方法为后文参数分析中不断重构更新钢轨廓形奠定了基础.
由图6可知:各方案拟合廓形轮轨接触点分布规律均与现有标准廓形一致;在心轨顶宽20 mm断面,各方案轮轨接触点位置基本一致,仅接触点发生突变对应的横移量有极小差别,说明在心轨顶宽较小时,各参数对轮轨接触点分布影响较小;在心轨顶宽50 mm断面,与标准廓形方案相比,方案1~4中轮轨接触点分布偏向于车轮踏面外侧,变化规律差别较小;而方案5~8则偏向于车轮踏面轮缘侧,变化规律也基本相近;随着心轨顶宽增大,参数k1对轮轨接触点分布的影响较k2和λ5大.
3.4 DOE分析
固定辙叉钢轨廓形影响轮轨匹配性能,因此,明确各控制参数对辙叉钢轨廓形影响权重,对辙叉钢轨廓形参数化设计中的参数调整优化具有重要意义. 以心轨20 mm断面和50 mm断面为例,通过DOE设计方法分析各控制参数对心轨廓形变化的影响效果及影响权重. 图7、8所示为心轨各控制参数对D1、D2和ΔS变化范围的影响. 图9为反映各控制参数对心轨廓形绝对影响效应的Pareto图. 图10为各控制参数对心轨廓形变化的影响权重.
由图7、8可知:在心轨20 mm断面,k2和λ5对ΔS、D1和D2变化范围均影响较大,k1对三者影响较小;在心轨50 mm断面,k1对ΔS、D1影响较大,对D2影响较小,k2和λ5对ΔS、D1和D2影响相当. 说明当心轨轨头宽度较小时,心轨整体廓形变化、心轨顶面和侧面廓形均主要受k2和λ5影响,k1影响较小;随着心轨轨头宽度增大,k1对心轨整体廓形和顶面廓形影响增大,对侧面廓形影响较小;k2和λ5对心轨各部位廓形影响效果相当.
由图9可知:在心轨20 mm断面,k1和λ5对ΔS影响效应较大,各参数对D1和D2影响效应相当且均较小;在心轨50 mm断面,各参数对ΔS影响效应最大,D1次之,对D2影响很小,其中,k1对ΔS影响效应最大,k2和λ5对ΔS影响效应相当;各参数对D1影响效应相当.
由图10可知:1) 在心轨20 mm断面,k1和λ5对ΔS和D1影响权重较小,而k2影响较大. 其中,ΔS对应的E分别为21.08%、56.89%和22.02%,D1对应的E分别为8.42%、61.95%、29.63%. 说明在心轨20 mm断面,k1和λ5对心轨整体廓形和顶面廓形的影响权重远不及k2,主要原因是在心轨轨头宽度较小时,侧面高度相对较大,心轨侧面廓形对心轨廓形变化起主导作用. 2) 在心轨50 mm断面,k1、k2和λ5对ΔS和D1的影响权重依次减小. 其中,ΔS对应的E分别为55.90%、33.38%、10.72%,D1对应的E分别为64.81%、26.71%、8.48%,说明随着心轨轨头宽度增大,各关键控制参数在方案所选取值范围内时,k1对心轨顶面和心轨整体廓形的影响权重逐渐明显,k2和λ5的影响逐渐减弱. 3) 在心轨20 mm断面,k1、k2和λ5对心轨侧面廓形偏差D2对应的E分别为9.93%、76.82%、13.25%;在心轨50 mm断面,k1、k2和λ5对心轨侧面廓形偏差D2对应的E分别为14.09%、66.04%、19.87%;由此可见,由于不同心轨轨头宽度对应的高度变化不大,因此,无论心轨宽度大小,心轨侧面廓形主要受k2影响较大,受λ5影响次之,受k1影响最小.
4. 结 论
1) 基于B样条理论,提出了一种考虑固定辙叉全断面钢轨廓形几何特征的参数化设计方法,该方法利用B样条曲线定义工业产品几何形状的优势,通过固定辙叉钢轨全断面廓形关键特征参数确定任意断面廓形控制点,进而实现钢轨廓形参数化重构;提出了钢轨重构廓形贴合度评价指标和关键控制参数对廓形贴合度影响权重指标,可直观量化地反映廓形参数变化对心轨廓形的影响.
2) 在心轨20 mm和50 mm断面,各方案心轨重构廓形与标准廓形贴合度良好,不同横移量下轮轨接触点分布规律符合实际. 各参数对轮轨接触点变化规律与心轨顶面宽度有关,在心轨顶宽较小时,各参数对轮轨接触点变化影响较小,在心轨顶宽较大时,k1对轮轨接触点分布的影响较k2和 λ5大.
3) 当心轨轨头宽度较小时,心轨整体廓形、心轨顶面和侧面廓形变化范围主要受k2和λ5影响,k1影响较小;k2对各部分廓形变化影响权重较大,k1和λ5影响较小;随着心轨轨头宽度的增大,k1对心轨整体廓形和顶面廓形变化范围和影响权重均逐渐增大,k2和λ5的影响逐渐减弱. 由此可见,各控制参数对心轨各部分廓形影响权重不同,通过考虑各控制参数对钢轨各部分廓形的影响权重,有针对性地重点考虑相应的控制参数进行廓形方案设计,避免了经验—性能计算—修正等反复试凑的传统设计过程,使辙叉廓形设计更具有针对性和可预测性.
-
表 1 DOE全阶乘方案设计
Table 1. DOE full factorial scheme design
方案 k1 k2 λ5 方案 1 1/30 1/6 1/4 方案 2 1/30 1/6 1/2 方案 3 1/30 1/4 1/4 方案 4 1/30 1/4 1/2 方案 5 1/10 1/6 1/4 方案 6 1/10 1/6 1/2 方案 7 1/10 1/4 1/4 方案 8 1/10 1/4 1/2 标准 1/20 1/5 1/3 -
[1] 高原,王平,陈嵘,等. 重载铁路固定辙叉区轮轨瞬态滚动接触行为分析[J]. 铁道学报,2020,42(6): 35-43.GAO Yuan, WANG Ping, CHEN Rong, et al. Analysis of 3D transient wheel-rail rolling contact behavior for fixed frog in heavy haul railways[J]. Journal of the China Railway Society, 2020, 42(6): 35-43. [2] 杨新文,张昭,孟玮,等. 道岔区固定辙叉心轨垂磨对轮轨动态接触的影响[J]. 同济大学学报(自然科学版),2020,48(11): 1595-1604.YANG Xinwen, ZHANG Zhao, MENG Wei, et al. Effect of vertical wear of unmovable frog nose rail on dynamical wheel-rail contact in turnout zone[J]. Journal of Tongji University (Natural Science), 2020, 48(11): 1595-1604. [3] MA Y W, MASHAL A A, MARKINE V L. Modelling and experimental validation of dynamic impact in 1: 9 railway crossing panel[J]. Tribology International, 2018, 118: 208-226. doi: 10.1016/j.triboint.2017.09.036 [4] WIEDORN J, DAVES W, OSSBERGER U, et al. Simplified explicit finite element model for the impact of a wheel on a crossing-Validation and parameter study[J]. Tribology International, 2017, 111: 254-264. doi: 10.1016/j.triboint.2017.03.023 [5] SKRYPNYK R, NIELSEN J C O, EKH M, et al. Metamodelling of wheel-rail normal contact in railway crossings with elasto-plastic material behaviour[J]. Engineering With Computers, 2019, 35(1): 139-155. doi: 10.1007/s00366-018-0589-3 [6] 于淼,张军,鹿广清,等. 重载货车作用下固定辙叉心轨磨耗分析[J]. 铁道学报,2014,36(7): 30-35.YU Miao, ZHANG Jun, LU Guangqing, et al. Analysis on wear of unmovable frog nose rail under the action of heavy-haul freight car[J]. Journal of the China Railway Society, 2014, 36(7): 30-35. [7] 王树国,葛晶,司道林,等. 固定辙叉查照间隔及心轨加宽研究[J]. 中国铁道科学,2014,35(1): 7-12.WANG Shuguo, GE Jing, SI Daolin, et al. Study on guard check gauge of fixed frog and design on widening nose rail[J]. China Railway Science, 2014, 35(1): 7-12. [8] 徐井芒,王平. 基于轮轨廓型的固定辙叉优化设计方法[J]. 中国铁道科学,2014,35(2): 1-6.XU Jingmang, WANG Ping. Optimization design method for rigid frog based on wheel/rail profile type[J]. China Railway Science, 2014, 35(2): 1-6. [9] 徐井芒,郑兆光,赖军,等. 轨道参数对高速道岔轮轨接触行为的影响[J]. 西南交通大学学报,2022,57(5): 990-999.XU Jingmang, ZHENG Zhaoguang, LAI Jun, et al. Influence of track parameters on wheel/rail contact behavior of high-speed turnout[J]. Journal of Southwest Jiaotong University, 2022, 57(5): 990-999. [10] 曹洋,王平. 固定辙叉心轨轨顶降低值优化[J]. 西南交通大学学报,2015,50(6): 1067-1073.CAO Yang, WANG Ping. Optimization of nose depth for rigid frog[J]. Journal of Southwest Jiaotong University, 2015, 50(6): 1067-1073. [11] 沈钢,汪徐江. 固定辙叉区心轨轨顶降低值优化设计方法[J]. 同济大学学报(自然科学版),2020,48(11): 1605-1611.SHEN Gang, WANG Xujiang. An optimization method for height of nose rail at turnout crossing[J]. Journal of Tongji University (Natural Science), 2020, 48(11): 1605-1611. [12] 张鹏飞,朱旭东,雷晓燕. 提速道岔辙叉翼轨的加高值方案优化[J]. 西南交通大学学报,2021,56(3): 602-610.ZHANG Pengfei, ZHU Xudong, LEI Xiaoyan. Optimization of wing rail lifting value for rigid frog of speed-up turnout[J]. Journal of Southwest Jiaotong University, 2021, 56(3): 602-610. [13] WAN C, MARKINE V L, SHEVTSOV I Y. Improvement of vehicle−turnout interaction by optimising the shape of crossing nose[J]. Vehicle System Dynamics, 2014, 52(11): 1517-1540. doi: 10.1080/00423114.2014.944870 [14] WAN C, MARKINE V L, SHEVTSOV I. Optimisation of the elastic track properties of turnout crossings[J]. Proceedings of the Institution of Mechanical Engineers Part F Journal of Rail and Rapid Transit, 2016, 230(2): 360-373. doi: 10.1177/0954409714542478 [15] 铁道部第三设计院. 道岔设计手册[M]. 北京: 人民铁道出版社, 1975: 413-503. [16] 韩江,江本赤,夏链,等. 基于轮廓关键点的B样条曲线拟合算法[J]. 应用数学和力学,2015,36(4): 423-431.HAN Jiang, JIANG Benchi, XIA Lian, et al. A B-spline curve fitting algorithm based on contour key points[J]. Applied Mathematics and Mechanics, 2015, 36(4): 423-431. [17] YANG Y, QIU W S, ZENG W, et al. A prediction method of rail grinding profile using non-uniform rational B-spline curves and Kriging model[J]. Journal of Central South University, 2018, 25(1): 230-240. doi: 10.1007/s11771-018-3732-9 [18] 温士明. 地铁线路小半径曲线钢轨打磨廓形研究[D]. 成都: 西南交通大学, 2018: 11-25. [19] 刘建桥. 地铁曲线轨道钢轨廓型优化研究[D]. 成都: 西南交通大学, 2019: 24-29. [20] 唐彦玲,吴磊,董勇,等. 重载铁路曲线钢轨廓形多目标优化设计[J]. 机械,2020,47(12): 1-9.TANG Yanling, WU Lei, DONG Yong, et al. Multi-objective optimization design of rail profiles on curve tracks of heavy haul railway[J]. Machinery, 2020, 47(12): 1-9. [21] 王亮,向伟彬,刘宏达,等. 小半径曲线钢轨非对称打磨目标型面优化[J]. 机械科学与技术,2021,40(6): 949-954.WANG Liang, XIANG Weibin, LIU Hongda, et al. Target profile optimization of asymmetrical grinding for rail with sharp-radius curve[J]. Mechanical Science and Technology for Aerospace Engineering, 2021, 40(6): 949-954. [22] 林凤涛,吴涛,杨洋,等. 高速铁路辙叉区钢轨打磨廓形设计方法[J]. 交通运输工程学报,2021,21(6): 124-135.LIN Fengtao, WU Tao, YANG Yang, et al. Design method of rail grinding profile in frog area of high-speed railway[J]. Journal of Traffic and Transportation Engineering, 2021, 21(6): 124-135. [23] 史振帅. 基于NURBS曲面拟合的道岔打磨廓形研究[D]. 南昌: 华东交通大学, 2021: 22-27. [24] 陈雨,安博洋,潘自立,等. 考虑尖轨变截面廓形的轮轨接触与磨耗分析[J]. 西南交通大学学报,2022,57(6): 1250-1258.CHEN Yu, AN Boyang, PAN Zili, et al. Analysis of wheel-rail contact and wear considering variable cross-sections of switch rail[J]. Journal of Southwest Jiaotong University, 2022, 57(6): 1250-1258. [25] 王璞,王树国,赵振华,等. 基于改善轮轨共形度的60 kg/m钢轨廓形优化[J]. 西南交通大学学报,2022,57(6): 1233-1238.WANG Pu, WANG Shuguo, ZHAO Zhenhua, et al. Optimization of 60 kg/m rail profile based on improving wheel-rail conformal degree[J]. Journal of Southwest Jiaotong University, 2022, 57(6): 1233-1238. -