Design and Matching Test of Two-Stage Sequential Supercharging System for Diesel Engines
-
摘要:
为改善低工况条件下因增压压力不足而造成柴油机缸内燃烧恶化问题,将两级增压(two-stage turbocharging,TST)与相继增压(sequential turbocharging,STC)技术相结合. 首先,通过推进特性试验进行两级相继增压与柴油机匹配计算,确定增压器型号;其次,利用GT-power软件得到增压器与柴油机联合运行曲线,验证方案的可行性;最后,进行台架试验,分析推进特性和负荷特性下两级相继增压对柴油机性能影响. 试验结果表明:推进特性条件下,在0~50%负荷区间,采用单涡轮增压(1TC)可以改善燃烧恶化问题;在50%~100%负荷区间,采用双涡轮增压(2TC)能够进一步提高空燃比,改善缸内燃烧,尤其在60%负荷下,燃油消耗率相比原机降低8.08%;两级相继增压模式在低转速、全扭矩及中高转速、高扭矩运行区域,可实现增压器的良好匹配.
Abstract:To improve the deterioration of combustion in the cylinder of diesel engines due to insufficient supercharging pressure under low operating conditions, two-stage turbocharging is combined with sequential turbocharging technology. First, the matching calculation of two-stage sequential turbocharging and diesel engine is carried out through propulsion characteristic tests, and the turbocharger model is determined. Secondly, the combined operation curve of the supercharger and diesel engine is obtained by using GT-power software to verify the feasibility of the scheme. Finally, the bench test is carried out to analyze the influence of two-stage sequential turbocharging on diesel engine performance under propulsion and load characteristics. The test results show that with propulsion characteristics, single turbocharging (1TC mode) can effectively improve the combustion deterioration under the engine load of 0–50%. Under the engine load of 50%–100%, the twin turbocharging (2TC mode) can further increase the air-fuel ratio and improve the in-cylinder combustion; especially at 60% engine load, the fuel consumption rate is 8.08% lower than that of the original engine. The two-stage sequential turbocharging mode can achieve good matching of turbochargers at low speed and full torque as well as medium-high speed and high torque.
-
Key words:
- diesel engine /
- two-stage sequential turbocharging /
- GT-power /
- supercharger matching
-
城市轨道交通具有运量大、节约环保、安全准时等特点. 截至2020年底,全国(不含港澳台)共有44个城市开通运营城市轨道交通线路233条,运营里程达7545.5公里,整体呈逐年上涨趋势. 目前城市轨道交通轨道结构主要采用现浇式,包括短枕式整体道床和长枕式整体道床:短轨枕式几何形位保持能力较差,施工精度较低;长枕式新旧混凝土接触面多,界面裂纹相对较多. 同时,现浇轨道结构均存在施工进度缓慢、混凝土质量控制难、作业环境恶劣、道床裂纹较多、维修困难、减振效果难以保证等诸多弊端[1-2].
为适应城市轨道交通的大规模建设,促进轨道结构技术创新,可结合国内高铁各类预制板式无砟轨道应用经验,引进并优化为适用于城市轨道交通的新型预制装配板式无砟轨道系统,解决目前现浇道床结构的各类缺点. 装配板式轨道相较于现浇式轨道具有七大优势:1) 工厂化预制质量更好、精度更高,在结构信息化方面优势明显;2) 机械化施工、现场拼装铺设轨道效率高、工期短;3) 施工占地小、城市交通压力小;4) 运营期间维修量少、可维修性好;5) 有减振需求时可快速实现减振升级改造;6) 能够应对我国即将到来的人力资源紧缺及人工成本持续提高的难题;7) 节能减排、降低能耗,能够有效降低施工对环境的污染. 此外,在轨道交通建设与运营中必须考虑减振问题,为最大限度减少振动造成对周围环境的影响,研发减振效果更好、结构更稳定的装配板式轨道十分必要[3-4].
目前已有部分学者对装配式轨道进行了相关研究. 李浩宇等[5]对城市轨道交通新型装配式无砟轨道的施工技术进行研究;冉丹丹等[6]对地铁装配式轨道板铺设的施工技术进行研究;李奇等[7]对装配式长型浮置板的静动力性能进行了研究;刘伟斌等[8-9]对城市轨道交通板式无砟轨道的设计应用进行相关研究;叶军等[10]对装配式轨道板的结构动力及减振特性进行研究;杨秀仁等[11]对城市轨道交通智能装配式减振轨道系统进行研究. 综上可知,既有针对城市轨道交通装配式轨道的研究主要集中在施工技术方面,对减振型装配式轨道设计研究相对较少,整体尚不够系统.
本文在系统总结分析装配式轨道设计需考虑因素的基础上,提出主要设计原则,并结合城市轨道交通减振需求,创新性提出一种减振型装配式无砟轨道设计方案,建立轨道动力分析模型和室内实尺试验模型,从理论和试验的角度对减振型装配式无砟轨道力学性能进行研究,以期为我国城市轨道交通减振型装配式无砟轨道结构设计提供参考.
1. 减振型装配式无砟轨道设计方案
1.1 设计原则
从既有城市轨道交通轨道结构运营实践来看,装配式轨道结构设计需考虑以下因素:
1) 小半径曲线适应性
城市轨道交通多运行在市区,受选线综合影响,小曲线半径较多,目前预制板式轨道结构在小半径曲线多利用扣件系统的调整量进行钢轨的圆矢度控制. 设计时需充分考虑针对小曲线应用需求,采用预制短板进行装配来适应线路调整,同时还需考虑外轨偏差、缓和曲线设置、轨道超高适应性等.
2) 施工吊装运输的便利性
城市轨道交通隧道地段铺设施工占比较大,预制轨道板尺寸及重量偏大时,在隧道内运输时需要专用吊装设备,洞内组装难度较大,施工效率偏低,平均铺设进度工效需进一步提升. 可在高速铁路轨道板设计尺寸参数等基础上进行优化调整,提高施工便利性和施工效率,兼顾施工质量与工效.
3) 桥上结构适应性
桥上铺设时需控制二期恒载,轨道结构长度需考虑温度效应. 对于露天路基地段及桥梁地段,轨道板运输及施工条件较好区域,轨道板长度可适当加长,进一步提高铺轨效率和轨道系统整体性能.
4) 减振区段适应性
城市轨道交通市区穿行,减振地段应用范围广泛,需充分考虑线路减振性能要求.
5) 排水适应性
受地下线盾构限界影响,现阶段预制板式轨道结构普遍采用双侧水沟,利用板两侧与隧道空隙设置排水沟,在小半径曲线地段内侧易造成排水沟尺寸较小,导致区间排水不畅,诱发相关病害. 可进一步优化轨道板结构尺寸,采用平板或框架板结构,适应中间排水与两侧排水应用需求.
6) 结构标准化
为进一步提升轨道建造经济性、降低轨道板制造费用,结合项目应用需求,考虑不同区段轨道板布置,通过优化结构尺寸及轨道板布置减少轨道板长度型号,形成标准化轨道板.
7) 运维便利性
城市轨道交通由于施工控制不佳或线路沉降超标,部分地段无法满足线路运营状态要求,需进行调整. 为降低对运营影响,需在天窗期实现单元板的切除、调运、更换、安装.
综合考虑上述因素,本文在城市轨道交通装配式无砟轨道设计时,主要遵循以下原则:1) 轨道结构简单、传力清晰,便于施工;2) 轨道结构各层匹配,无明显薄弱环节;3) 轨道结构采用工厂化制造、机械化施工,提高工程质量及施工工效;4) 具备较好的可维修性;5) 轨道结构减振实现方便且减振效果良好;6) 考虑一定的经济性;7) 充分考虑与相关工程的接口技术要求,系统设计.
1.2 设计方案
减振型装配式无砟轨道结构主要组成部件包括钢轨、扣件、预制轨道板、自密实混凝土层、减振垫、板间连接装配件、限位构件和基座等,如图1所示.
1) 钢轨及扣件
采用60 kg/m钢轨,材质为U75V. 采用无挡肩的DTVI2型扣件,如图2所示,扣件间距600 mm.
2) 轨道板
轨道板为工厂预制非预应力、单元分块式框架板结构. 板面设4对承轨台,板底设2排共14根门型筋,混凝土强度等级为C60. 轨道板分中间板和端板两种型式,中间标准轨道板宽2100 mm、厚200 mm、长2100 mm,端板长2200 mm. 板中间框架沿线路纵向长900 mm,沿横向宽800 mm,倒角半径200 mm. 框架板中间部分可作为排水沟检查使用,设置盖板作为疏散通道使用. 轨道板端部纵向预留2组共8根连接钢筋. 两块端板加上N块中间板通过板端连接钢筋和板间连接装配件装配形成较长的装配单元,如图3所示. 装配单元之间设100 mm伸缩缝,板端预埋接地端子.
3) 板间连接装配件
板间连接装配件为钢筋套筒,如图4所示. 板间连接装配件连接轨道板端外伸的连接钢筋装配形成一定长度的装配单元,形成的接头需符合《钢筋机械连接技术规程》[12]的Ⅰ级接头指标要求.
4) 自密实混凝土垫层
板下自密实混凝土垫层厚100 mm,强度等级为C40,起到调平、填充轨道板与下部基础的作用,如图5所示. 垫层宽度650 mm,与框架板两侧宽度相同,纵向与板长相同. 垫层内配置单层ϕ10钢筋网片,垫层混凝土浇筑施工时通过板一端向另一端进行灌注.
5) 减振垫
减振垫设置在自密实混凝土垫层与基座之间,可一定程度提高减振垫上部轨道结构的参振质量,提高减振效果. 减振垫采用聚氨酯发泡材料或者具有锥形突起结构的橡胶材料,减振层厚度为30 mm,减振垫刚度0.019 N/mm3.
6) 限位构件
限位构件现浇于装配单元中的中间板与中间板之间或中间板与端板之间,与自密实混凝土垫层共同浇筑形成整体. 限位构件分左右两块,如图6(a)所示,单块宽300 mm、长650 mm、高300 mm,起纵横向限位作用. 限位构件混凝土等级为C40,设钢筋笼,通过8根ϕ 20长540 mm的限位钢筋与下部混凝土基座相连,限位钢筋外包裹一层套筒,如图6(b)所示,钢筋可沿套筒长度方向自由伸缩,以此适应垂向减振的变形需求.
7) 基座
基座为现浇钢筋混凝土结构,混凝土强度等级为C40. 基座在限位构件位置植有限位钢筋,钢筋植入深度250 mm. 考虑到城市轨道交通曲线半径较小,超高较大,轨道结构超高通过基座实现.
总体来说,减振型装配式无砟轨道采用预制框架板型式,轨道板重量较轻,施工运输方便;板中框架可作为排水沟的观察孔和清理孔,盖上盖板还可作为疏散通道;单块轨道板在线路纵向上长度较短,对小曲线半径适应性较强;板与板之间通过装配形成较长的装配单元,结构施工精度高、整体性好、可靠性强,同时增加了轨道单元的整体长度,减振效果进一步加强;板端限位构件与自密实混凝土浇筑形成整体,限位能力强、稳定性能好;后期养护维修时,可通过解锁装配单元完成对轨道板单元的更换调整,可维修性较强.
2. 减振型装配式无砟轨道动力特性
结合轨道结构设计方案特点,建立可考虑新型装配式无砟轨道各结构层相互作用关系的车辆-轨道耦合动力分析模型,研究轨道结构在列车荷载作用下的动力特性.
2.1 新型装配式无砟轨道动力分析模型
基于车辆-轨道动力学耦合理论[13],建立车辆-装配式无砟轨道空间耦合动力分析模型. 模型包括车辆、装配式无砟轨道、轮轨耦合、轨道不平顺共4个部分. 车辆采用整车模型,考虑车体和转向架的沉浮、横摆、点头、摇头、侧滚以及轮对的沉浮、摇头、横摆和侧滚共计31个自由度,车辆参数参照地铁B型车辆[14]. 新型装配式无砟轨道各组成部件的几何尺寸按设计取值,材料参数见表1.
表 1 轨道结构模型参数Table 1. Track structure model parameters部件 弹性模量/GPa 泊松比 密度/(kg·m−3) 钢轨 206.0 0.3 7800 轨道板 36.5 0.2 2500 混凝土垫层 34.0 0.2 2400 基座 34.0 0.2 2500 限位构件 34.0 0.2 2500 为更好反映各部件受力情况,钢轨、轨道板、板下自密实混凝土垫层、基座和限位构件均采用实体单元模拟. DTVI2型扣件采用三向非线性弹簧-阻尼单元模拟,考虑纵、横、垂3个方向的刚度,垂向刚度取30 kN/mm,切向静刚度取40 kN/mm[15]. 轨道板通过门型筋与自密实混凝土层相连,两界面采用绑定连接. 减振垫与自密实混凝土界面切向采用摩擦接触模拟,摩擦系数取0.6,法向采用硬接触模拟[16],减振垫与基座界面采用绑定连接,减振垫刚度取0.019 N/mm3,以隧道地段为例,支撑刚度取1200 MPa/m[16-17]. 钢轨两端采用对称约束,轮轨间接触遵循Hertz非线性接触理论,轨道不平顺采用美国普通铁路的美国五级谱[18]. 综上,新型装配式无砟轨道车辆-轨道空间耦合动力有限元分析模型如图7所示.
2.2 轨道结构动力特性
行车速度100 km/h条件下,新型装配式无砟轨道的部分典型动力指标时程曲线如图8所示,系统动力响应汇总如表2、3所示.
表 2 车-轨动力响应统计Table 2. Statistics of vehicle-rail dynamic response动力响应指标 峰值 轮轨横向力/kN 5.457 轮轨垂向力/kN 73.115 脱轨系数 0.073 轮重减载率 0.108 表 3 系统动力响应统计Table 3. Statistics of system dynamic response项目 动力响应指标 峰值 车体 垂向加速度/(m·s−2) 0.429 横向加速度/(m·s−2) 0.285 钢轨 垂向加速度/(m·s−2) 190.380 垂向位移/mm 2.141 轨道板 垂向加速度/(m·s−2) 31.200 垂向位移/mm 0.812 混凝土垫层 垂向加速度/(m·s−2) 26.020 垂向位移/mm 0.761 基座 垂向加速度/(m·s−2) 1.160 垂向位移/mm 0.122 由图8和表2、3可知:新型装配式无砟轨道在行车速度100 km/h条件下,脱轨系数和轮重减载率最大分别为0.073、0.108,轮轨垂向力和横向力峰值分别为73.115 kN和5.457 kN,结果均未超限,满足行车安全要求;车体垂、横向加速度分别为0.429 m/s2和0.285 m/s2,满足《地铁设计规范》(GB 50157—2013)[19]垂向加速度az≤0.13g和横向加速度ay≤0.10g的行车安全性和乘坐舒适性指标要求;钢轨最大垂向位移为2.141 mm,轨道板垂向位移和自密实混凝土层相差不大,分别为0.812 mm和0.761 mm. 基座位移很小,约0.122 mm;钢轨、轨道板、自密实混凝土层和基座的振动加速度从上往下均有不同程度的衰减,垂向加速度峰值分别为190.380、31.200、26.020、1.160 m/s2.
综上可知,新型装配式无砟轨道结构在列车荷载作用下动力性能良好,满足行车安全性和舒适性要求.
3. 减振型装配式无砟轨道试验研究
借助试验室静动力加载装置以及落轴试验装置,对轨道结构进行静载试验、疲劳试验和落轴冲击试验,研究在静载、疲劳动载和落轴冲击载荷作用下,减振型装配式无砟轨道的力学性能,为新型装配板式无砟轨道的设计以及实际应用提供试验支撑. 减振型装配式无砟轨道实尺模型由3块轨道板(2块端板 + 1块中间板)组成,总长7.1 m. 现场部分测点布置如图9所示.
3.1 静载试验
轨道系统静载试验考虑的最大载荷为3Pj (Pj为列车静轮载,轴重17 t),包含规范推荐的设计垂向荷载2Pj. 从0加载至510 kN,在400 kN前每一级荷载为50 kN,在400 kN后每一级荷载为20 kN,记录每次载荷作用下轨道系统各测点的应变情况. 各级荷载作用下典型测点的混凝土应变变化曲线如图10所示,其中测点1、2分别位于轨道板横向两侧截面中间位置,测点3~6等间距分布于轨道板上表面板中位置.
由图10可知:各测点混凝土应变绝对值随荷载的增加整体呈线性增加趋势;在3Pj的静载作用下,轨道混凝土最大压应变为 −129.6 με,远小于混凝土极限应变值;在静载试验过程中轨道结构未出现开裂现象,减振型装配式轨道结构强度和刚度满足承载要求.
图11为各级荷载作用下轨道结构变形曲线. 由图11可知:各测点垂向位移均随着荷载的增大而增大,左、右两边钢轨最大垂向位移分别为2.44 mm和2.40 mm,轨道板两侧的最大垂向位移分别为1.32 mm和1.34 mm,数值相差不大,轨道结构两侧对称受力. 图12为加、卸载过程中钢轨位移变化曲线,由图12可知:在加载卸载过程中曲线总体平稳,未见位移的突变;卸载后基本没有残余位移,说明加载、卸载过程中轨道板主要处于弹性受力阶段,未出现较大的塑性变形,受力性能良好.
3.2 疲劳试验
采用MTS疲劳试验机进行疲劳试验,疲劳荷载的上限按规范要求取为1.5Pj,疲劳荷载下限值按疲劳机最大加载能力的3%确定,即15 kN. 通过MTS作动器施加正弦荷载,加载频率为4 Hz,累计疲劳作用次数为200万次.
图13为钢轨和轨道板位移随疲劳试验次数变化曲线. 由图12可知:在1.5Pj的200万次疲劳荷载作用下,钢轨最大垂向位移为1.303 mm,轨道板最大垂向位移为0.843 mm;疲劳试验过程中,钢轨和轨道板位移在疲劳荷载作用下随疲劳次数的增加在一定范围内上下波动,200万次疲劳试验后钢轨最大垂向位移增加了0.180 mm,轨道板最大垂向位移增加了0.130 mm,整体变化量较小,表明装配式无砟轨道在经200万次疲劳荷载作用下的工作性能良好.
在疲劳试验过程中,轨道结构表面及层间界面未见肉眼可见裂缝,减振型装配式无砟轨道抵抗疲劳性能良好,新型装配式无砟轨道能够经受1.5Pj的疲劳荷载作用200万次的检验.
3.3 落锤冲击试验
采用力锤锤击法,分别测试一般地段和减振地段装配式无砟轨道结构在力锤脉冲荷载作用下的动力响应,对比分析轨道各部件在脉冲荷载作用下的Z振级,掌握减振型装配式无砟轨道振动传递特性和减振性能. 试验采用DFC-2型力锤,重量8 kg,力锤垂直下落高度为120 mm. 为避免轨道板纵向边界效应的影响,并且充分激发轨道板的垂向振动,将落锤点设置在轨道板板中截面.
力锤作用下一般地段和减振地段装配式无砟轨道基座垂向加速度1/3倍频程如图14所示.
由图14可知:减振型装配式无砟轨道具有显著的减振效果,频率为100~125 Hz,最大减振效果可达12.4 dB.
将力锤锤击试验工况下的轨道结构系统简化为单一激励输入的质量弹簧系统,分析力锤锤击激励输入与基座振动响应输出间的传递函数,对比一般地段与减振地段新型装配式无砟轨道传递函数幅值,如图15所示.
由图15可知:相比一般地段,减振地段的新型装配式无砟轨道可使力锤锤击激励与基座振动响应间的传递函数幅值大幅削减,减振型装配式无砟轨道减振效果较好;疲劳加载前、疲劳加载至100万次、200万次工况下的传递函数幅值曲线峰值起伏频率范围基本一致,表明轨道结构的减振特性经过疲劳加载后的稳定性较好.
4. 结 论
1) 减振型装配式无砟轨道结构在100 km/h的行车条件下,系统动力响应均在安全限值范围内,满足行车安全性和舒适性要求.
2) 轨道结构在静载试验过程中混凝土应变随荷载增加整体呈线性增加趋势,最大压应变为−129.6 με,远小于混凝土极限应变值,结构受力在强度允许范围内;卸载后残余变形很小,轨道弹性较好;静载试验过程中结构未见肉眼可见裂纹,结构承载能力和抗裂性能满足要求.
3) 在200万次疲劳荷载作用下,轨道结构整体受力变形较小,疲劳试验过程中变化幅值不大,轨道结构表面及层间界面未见裂缝,轨道结构抗疲劳性能良好.
4) 减振型装配式无砟轨道减振效果显著且减振稳定性较好,频率为100~125 Hz,最大减振效果可达12.4 dB.
致谢:中铁第四勘察设计院集团有限公司科技研发课题(2019K011,2020K155).
-
表 1 柴油机主要技术参数
Table 1. Main technical parameters of diesel engine
参数 参考值 标定功率/kW 444 标定转速/(r·min−1) 1800 气缸排列方式 V 型,60° 夹角 缸径/mm 128 冲程/mm 140 排量/L 28.8 压缩比 15∶1 发火次序 参考文献[6] 增压器型号 J100 -
[1] WINTRUFF I, BÜCHELER O, RALL H, et al. The new MTU series 4000 rail engines certified for EU IIIB[J]. MTZ Worldwide, 2011, 72(5): 30-34. doi: 10.1365/s38313-011-0050-0 [2] 王银燕,何清林,王贺春,等. 柴油机采用4台增压器相继增压性能的试验研究[J]. 内燃机工程,2012,33(1): 6-10,17. doi: 10.3969/j.issn.1000-0925.2012.01.002WANG Yinyan, HE Qinglin, WANG Hechun, et al. Experimental research on performance of four-turbocharger sequential turbocharged diesel engine[J]. Chinese Internal Combustion Engine Engineering, 2012, 33(1): 6-10,17. doi: 10.3969/j.issn.1000-0925.2012.01.002 [3] SCHMITT F. Powerful turbocharging system for passenger car diesel engines[J]. MTZ Worldwide, 2014, 75(3): 12-19. doi: 10.1007/s38313-014-0028-9 [4] 魏名山,张志,方金莉,等. 带有放气阀的二级增压系统的设计与试验[J]. 内燃机学报,2009,27(2): 166-170. doi: 10.3321/j.issn:1000-0909.2009.02.012WEI Mingshan, ZHANG Zhi, FANG Jinli, et al. Design and experiments for a regulated two-stage turbocharging system with a wastegate valve[J]. Transactions of CSICE, 2009, 27(2): 166-170. doi: 10.3321/j.issn:1000-0909.2009.02.012 [5] 高占斌,王彬彬,高敬博,等. 掺混生物柴油对相继增压柴油机燃烧性能的影响[J]. 西南交通大学学报,2022,57(4): 791-796. doi: 10.3969/j.issn.0258-2724.20200544GAO Zhanbin, WANG Binbin, GAO Jingbo, et al. Emission characteristics of successively supercharged diesel engine blended with biodiesel[J]. Journal of Southwest Jiaotong University, 2022, 57(4): 791-796. doi: 10.3969/j.issn.0258-2724.20200544 [6] 李先南. 船用V型柴油机大小增压器相继增压系统设计研究[D]. 哈尔滨: 哈尔滨工程大学, 2012. [7] 吴鹏. 柴油机米勒循环与两级涡轮增压的研究[D]. 哈尔滨: 哈尔滨工程大学, 2013. [8] 刘博. 柴油机可调二级涡轮增压系统研究[D]. 上海: 上海交通大学, 2011. [9] 钱跃华. 双涡轮增压系统匹配方法和优化控制的研究[D]. 上海: 上海交通大学, 2015. [10] 顾宏中. 涡轮增压柴油机性能研究[M]. 上海: 上海交通大学出版社, 1998: 7-14. [11] 顾宏中,郭中朝. 16VPA6STC柴油机改用MIXPC涡轮增压系统研究[J]. 中国造船,2008,49(1): 90-96. doi: 10.3969/j.issn.1000-4882.2008.01.014GU Hongzhong, GUO Zhongchao. A study of MIXPC turbocharging system on 16VPA6STC diesel engine[J]. Ship Building of China, 2008, 49(1): 90-96. doi: 10.3969/j.issn.1000-4882.2008.01.014 [12] LENG L, QIU H J, LI X N, Et al. Effects on the transient energy distribution of turbocharging mode switching for marine diesel engines[J]. Energy, 2022, 249: 123746.1-123761.11. [13] LIU R L, YANG C H, ZHANG Z J, et al. Experimental study on thermal balance of regulated two-stage turbocharged diesel engine at variable altitudes[J]. Journal of Thermal Science, 2019, 28(4): 682-694. doi: 10.1007/s11630-019-1151-3 [14] WANG J, SHEN L Z, BI Y H, et al. Power recovery of a variable nozzle turbocharged diesel engine at high altitude by response surface methodology and sequential quadratic programming[J]. Proceedings of the Institution of Mechanical Engineers, Part D: Journal of Automobile Engineering, 2019, 233(4): 810-823. doi: 10.1177/0954407017753913 -