HXD2型电力机车轮对压装过程仿真分析.pdf

4 机 械设 计 与制 造 Ma c hi n e r y De s i g n&Ma n u f a c t u r e 第 5期 2 0 1 4年 5月 H X D 2 型电力机车轮对压装过程仿真分析 蔡慧 ， 赵永翔 ， 梁红琴 ( 1 ． 西南交通大学 牵引动力国家重点实验室， 四川 成都6 1 0 0 3 1 ； 2 ． 西南交通大学 机械工程学院， 四川 成都6 1 0 0 3 1 ) 摘要： 采用大型有限元分析软件 A N S Y S对轮对压装过程进行接触问题的动态有限元分析， 模型采用 H X D ， 型电力 机车的轮对模型， 轮轴材料采用课题组得到的轮轴应力应变本构曲线。通过采用不同的过盈量， 摩擦系数进行仿真分 析， 得到其对轮毂孔、 轮座面的Mi s e s 应力沿轮轴接触面轴向位置的分布的影响， 以及其对压装压力曲线的影响。同时 还将根据得到的压装压力曲线及实际图纸要求的过盈量来为本轮对压装时采用的摩擦系数得出一定的范围，供实际 应用参考。 关键词： 轮对； AN S YS ； 摩擦系数； 过盈量： 压装压力曲线 中图分类号： T H1 6 ； U 3 6 4 ． 8 文献标识码： A 文章编号： 1 0 0 1 — 3 9 9 7 ( 2 0 1 4 ) 0 5 — 0 0 0 4 — 0 4 Th e W h e e l s e t F i ~ i n g Si mu l a t i o n o f HXD2 E l e c t r i c L o c o mo t i v e ， CAI Hui ，ZHAO Yo n g - x i a n g ，L I ANG Ho n g - q i n 。 ( 1 ．T r a c t i o n P o w e r S t a t e K e y L a b o r a t o r y o f S o u t h w e s t J i a o t o n g U n i v e r s i t y ， S i c h u a n C h e n g d u 6 1 0 0 3 1 ， C h i n a ； 2 ． S c h o o l o f Me c h a n i c M E n g i n e e r i n g ， S o u t h w e s t J i a o t o n g U n i v e r s i t y ， S i c h u a n C h e n g d u 6 1 0 0 3 1 ， C h i n a ) A b s t r a c t ： T h e fin i t e e l e m e n t a n a l y s i s s o fi w a F e A N S Y S W O $u s e d t o m a k e t h e d y n a m i c fi n i t e e l e me n t a n a l y s i s o f c o n t a c t p r o b l e m d u r i n g w h e e l s e t fit t i n g ． HX D 2 e l e c t r i c l o c o m o t i v e w h e e l s e t m o d e l w i l l b e u s e d ． U s i n g t h e s t r e s s - s t r a i n c u r v e s o f t h e m a t e r i a l s w h i c h i s o b t a i n e d b y s t u d y i n g t e a m． U s i n g d if f e r e n t a mo u n t o f i n t e d e r e n c e ， f r i c t i o n c o e f fic i e n t t o ma k e S i mu l at i o n a n a l y s is t o g e t t h e a c c u r a t e s t r e s s d i s t r i b u t i o n o ft h e s u o r ace o fw h e e l - h o l e a n d w h e e l - s e a t d o n g t h e a x i a l p o s i t i o n， ∞w e l l a s t h e p r e s s u r e c u r v e ． A t t h e s a m e t i me w e w i l l use t h e p r e s s u r e c u rve and t h e act u a l a mo u n t of i n t e de r e n c e t o g e t a c e r t ai n r ang e o f f r i c t i o n c o e f f wi e n t ， f o r t h e p r act i c a l a p p l i c ati o n o fr e f e r e n c e ． Ke y W o r d s ： W he e l s e t ； ANS YS ； Fr i c t i o n Co e f fic i e n t ； I n t e r f e r e n c e ； Pr e s s u r e Cu r v e 1引言 随着铁道车辆速度以及轴重的不断提高，对轮轴安全性的 要求也越来越高。研究表明， 轮轴配合部位的应力状态对轮轴的 疲劳强度具有重要的影响。 因此对轮轴的压装过程的研究具有相 当重要的意义。 文献f l从配合面的光洁度、 几何形状误差、 尺寸误差、 压入速 度、 温度、 润滑、 材料的机械性能和压装的操作过程等多个因素探 讨了其对压装质量的影响_1 ]。 文献 轮对压装中出现的压装曲线 不合格等问题作了相应的分析和探讨。 文献嗍 A B A Q U S软件建 立了转 K 6型转向架轮对压装力学模型， 通过对轮对压装过程进 行仿真模拟计算，分析了压装参数对轮对压装质量的影响规律。 文献吲 针对动车组轮对压装出现的不合格 F — S曲线， 分析了E N 标准动车组轮对压装 F — S曲线的正确判定方法 ，并提出了建议。 文献圈 对 H X N 型机车轮对压装过程进行弹塑性数值模拟， 分析 了摩擦系数及过盈量对压装应力的影响， 以及油槽的存在对应力 分布的影响。 文献 阐述了轮对突悬组装的形成、 突悬量的确定及 压装压力曲线理论长度计算方法， 提出了在生产实际中解决突悬 组装对压装压力曲线影响的方法。 文献[7 1介绍了轮对裂纹新的非破坏性检测技术感应电流的 集中潜能下沉( I C F P D) 技术， 并进行了铁路车轴轴向和径向的裂 纹检测。 文献【田 介绍了在压装曲线的预测中， 简化分析模型和有限 元分析模型的各自的优缺点， 同时也详细的介绍了如何利用试样 进行摩擦测量的方法， 以及压装参数对压装质量的影响。 采用大型有限元分析软件 A N S Y S 对 H X D ： 型机车轮对压 装过程进行仿真分析， 轴毂过盈配合考虑了接触非线性， 同时轮 对材料采用的是弹塑性本构关系曲线，压装时车轮一端固定， 在 车轴端面采用位移控制加载方式， 计算过程同轮对实际压装的条 件接近 。 2轮对压装基本分析方案 2 ． 1基本数据的确定 车轴材料为 A1 N ， 车轮材料为R 7 T ， 由于进行弹塑胜分析， 材 料采用课题组得到的应力应变本构曲线。 车轴与车轮采用过盈配合 来稿 日期 ： 2 0 1 3 — 1 1 - 0 5 基金项目： 国家自 然科学基金( 5 1 1 7 5 4 3 9 ) ； 教育部博士点基金( 2 0 1 1 0 1 8 4 1 1 0 0 0 3 ) 作考简介： 蔡慧， ( 1 9 9 0 一 ) ， 男， 浙江丽 水人， 硕士 研究生， 主 要研 究方向： 强度可 靠性研究 赵永翔， ( 1 9 6 3 一 ) ， 男， ~ 1 ) l l成都人， 博士生导师， 主要研究方向： 强度可靠性研究 第5 期 蔡慧等： H X D ： 型电力机车轮对压装过程仿真分析 5 连接 ， 车轮与车轴配合面的直径为 2 5 0 m m， 由T B f F 1 7 1 8 铁道车 辆轮对组装技术条件可知： 过盈量为轮座直径的( 0 ． 0 8 1 ．5 ) ％， 计 算得过盈量为 ( O ． 2 0 — 0 ． 3 7 5 ) mm，取以下 7 个过盈量进行研究： 0 ． 2 0 、 0 ． 2 2 、 0 ． 2 5 、 0 ． 2 8 、 0 ． 3 1 、 0 ． 3 4 、 0 ． 3 7 。摩擦系数取 ： 0 ． 1 、 0 ． 1 3 、 0 ． 1 6 、 0 ． 1 9 、 0 ． 2 2 、 0 ． 2 5 、 0 ． 2 8 。由于压装速度的改变将涉及到摩擦传热等 许多其它关于热方面的东西， 需要考虑的东西相当复杂， 且由先 前国内外的研究可知， 压装速度对结果的影响并不大， 就固定一 个相对较低的压装速度， 这样就可以忽略有摩擦引起的热应力的 影响， 压装速度的详细的研究有待以后进一步探讨。 2 _ 2有限元模型及压装分析过程 实体采用 3 I ) " - S o l i d 1 8 5网格单元模型， 轮轴的接触， 套筒与 轮毂端面的接触采用面面接触单元 C O N T A C T 1 7 4和T A R G E — T 1 7 0模拟。 有限元模型加载时采用固定轮毂一端面的轴向位移，具体 处理采用一个套筒， 套简一端固定， 另一端与轮毂端面接触 ， 限制 其轴向位移， 此接触的摩擦系数取为 0 ． 2 ， 在轴的右端面施加一个 移动位移， 模拟轴的装入过程。如图 1 所示。 图 1轮对有限元模型图 Fi g ． 1 W h e e l s e t F i n i t e El e me n t Mo de l 3相关参数对轮对压装接触应力及压装 力的影响 3 ． 1过盈量的影响 通过模拟不同的过盈量下轮对的压装过程，分析得到以下 结果。固定摩擦系数为0 ． 1 9 。轮毂孔和车轴轮座面表面Mi s e s 应 力沿轮轴接触面轴向位置的分布情况， 如图 2 、 图 3 所示。由图 2 和图 3中应力曲线可知过盈量的大小直接影响轮轴接触配合面 应力的大小， 过盈量越大， V o n M i s e s 应力越大， 但是过盈量对应 力分布的趋势基本无影响。 兰 0-2 2 0 ． 2 5 0 ． 2 8 O _ 3 1 0． 3 4 0 _ 3 7 0． 2 图2不同过盈量下轮毂孑 L 内表面应力的分布情况 Fi g ．2 S t r e s s Di s t ri b u t i o n o f t h e Wh e e l Bo r e S u rfa c e i n Di f f e r e nt I n t e rfe r e n c e 三 U U儿H J l 5 U 2 Uu L ( mm ) 图 3不同过盈量下轮座表面应力的分布情况 F i g ． 3 S t r e s s Di s t rib u t i o n o f t h e W h e e l S e a t S u rfa c e i n Di ffe r e nt I n t e rfe r e n c e 通过提取不同过盈量下的压装过程中多个位置的压装力， 得到以下不同过盈量下的压装压力曲线， 如图 4 所示。为了让读 者能更加清楚的了解不同过盈量下最大压装力的大小比较， 列出 其数值 ，如表 1 所示。由表 1中的数据可知 ：当过盈量增加 0 ．0 3 mm时， 最大压装力大约提高了2 3 0 k N， 但对压装压力曲线的 形状基本上无影 响。 L( mm) 图 4不同过盈量下的压装压力曲线 F i g ．4 Th e P r e s s u r e Cu r v e i n Di f f e r e n t I n t e rfe r e n c e 表 1不同过盈量下最大压装力 T a b ． 1 Ma l l[Pr e s s -F i t F o r c e i n Di ff e r e n t I n t er f e r e n c e 3 ． 2摩擦系数的影响 。 研究时固定过盈量为 0 ． 2 8 ，对 6 种不同的摩擦系数进行仿 真分析。 轮毂孔和车轴轮座面表面 M i s e s 应力沿轮轴接触面轴向 位置的分布情况， 如图5 、 图 6 所示。 兰 +0． 1 —0． 1 3 +0 ． 1 6 一 0 ． 1 9 — 一 0 ． 2 2 — 一0 ． 2 5 —0 ．2 8 图5不同摩擦系数下轮毂子 L 内表面应力的分布情况 F i g ．5 S t r e s s Di s t r i b u t i o n o f t h e W h e e l B o r e S u rfa c e i n Di f f e r e n t Fr i c t i o n Co e f f i c i e n t ∞ ∞ 加 ∞ ∞ ∞ 加 0 6 机 械设 计 与制 造 No ． 5 Ma v ． 2 01 4 +0 ． 1 一O． 1 3 +O ． 1 6 —0 ． 1 9 —* 一 0 ． 2 2 +0 ． 2 5 — 。0 ． 2 8 ￡ ( ram) 图6不同摩擦系数下轮座表面应力的分布情况 Fi g ．6 S t r e s s Di s t r i b u t i o n o f t h e Wh e e l S e a t S u r f a c e i n Di f f e r e n t F r i c t i o n Co e f f i c i e n t 由图5和图6可知，不同摩擦系数下应力分布趋势基本相 同，且摩擦系数对轮毂孔内表面及轮座表面的应力影响较小， 因 此压装时在轴毂配合面涂上润滑油尽管对降低表面接触应力的 作用不明显 ， 但却能起到防止擦伤的作用。 不同摩擦系数对压装压力曲线的影响， 如图7所示。 0． 1 +0 ． 1 3 ’ 一0 ． 1 6 一* 0 ． 1 9 — *’0 ． 2 2 * ． O ． 2 5 O ． 2 8 L ( m m) 图 7不同摩擦系数下的压装压力曲线 F i g ．7 T h e Pr e s s u r e Cu r v e i n Di f f e r e n t F r i c t i o n Co e ff i c i e n t 为了让读者能更加清楚的了解不同摩擦系数下最大压装力 的大小比较 ， 列出其数值， 如表 2 所示。 表2不同摩擦系数下最大压装力 Ta b ．2 Ma x P r e s s —F i t F o r c e i n Di ff e r e n t F r ic t io n Co e ff i c i e n t 由表 2可知， 当摩擦系数增加 0 ．0 3时， 最大压装力大约增加 了 3 0 0 k N， 说明其对压装压力的影响很大， 但对压装压力曲线的 形状基本无影响。 4对研究的 H X D 轮对的压装的探讨 由于实际压装时， 无法观察到轮轴结合面的应力情况 ， 因此 采用压装曲线来判断压装的质量，由于影响压装曲线的因素很 多， 且由上文的研究及相关文献可知， 摩擦系数以及过盈量是其 中最重要的两个，下面将在图纸上给出的过盈量范围的基础上， 仿真研究采用合适的摩擦系数使压装曲线及最大压装力满足 T B ／ T 1 7 1 8 — 2 0 0 3 ( 铁道车辆轮对组装技术条件》 的要求。 4 ． 1过盈量的确定 在进行轮对过盈量的设计时， 主要考虑以下两个方面： ( 1 ) 为了保证行车的安全， 最小过盈量应该能保证车轮与车 轴过盈配合传递牵引力和横向力的要求。 ( 2 ) 最大过盈量应保证在轮对组装时不破坏轴毂结合处应 力分布的均匀性 ， 以避免轮座早期产生微小的裂纹。 基于以上两点以及铁道标准的规定，大同电力机车有限公 司给出了压装时轮对过盈量的范围为( O - 3 l 1 — 0 ． 3 6 9 ) mm。 4 _2根据标准的规定确定压装时合适的摩擦系数 由于T B ／ T 1 7 1 8 — 2 0 0 3 ( 铁道车辆轮对组装技术条件》 规定， 压 装后检查验收是否合格的标准是压装曲线。 根据其关于最终压装 压力和压装压力曲线的规定，车轴按轮毂孔直径每 1 0 0 ra m的最 终压装力应为( 3 4 3 ～ 5 8 8 ) k N ， 计算得所使用车轴的最终压装力为 ( 8 5 7 ． 5 ~ 1 4 7 0 ) k N 。 首先采用过盈量为 0 ．3 1 1 进行模拟仿真， 得到合适的摩擦系 数， 使其得到的最大压装力小于 1 4 7 0 k N； 再采用过盈量为 0 ． 3 6 9 进行模拟仿真， 得到合适的摩擦系数， 使其得到的最大压装力大 于8 5 7 ．5 k N， 这样就可以得到压装时摩擦系数的合适范围。 4 ． 2 ． 1过 盈量 为 0 - 3 1 1 时 取过盈量为0 ． 3 1 1 ，采用了多种不同的摩擦系数进行实验， 当采用摩擦系数为 0 ． 1 2 5时， 得到的最大压装力为 1 4 3 5 k N ， 满足 要求且接近于 1 4 7 0 k N， 得到的压装曲线， 如图 8 所示。 2 0 o O l 5 o o 至l O O O 5 0 o No ． 5 Ma v ． 2 0 1 4 机 械 设 计 与制 造 7 经测量，末端下降曲线长为 9 ． 3 m m，小于曲线投影长度的 1 0 ％( 1 8 ．2 7 mm) ， 满足要求。 压力下降值为 1 9 k N， 小于最大压力的 5 ％( 7 3 ． 5 k N) 。经校核， 压装压力曲线符合标准。 Z ￡( m m) 图 9压装压力 曲线 Fi g ． 9 P r e s s - Fit Pr e s s u r e C ur v e 通过上文的分析 ，建议实际压装时 ，摩擦系数的范围为 ( 0 ． 0 6 5 — 0 ． 1 2 5 ) 。 5结论 根据图纸要求建立了轮轴的模型，然后采用大型有限元分 析软件 A N S Y S 对轮对建立了有限元模型， 进行仿真分析， 对轮轴 进行了强度分析， 论文研究结论归纳如下： ( 1 ) 对轮对压装过程中不同过盈量及摩擦系数下轮轴应力 的结果分析表明：过盈量对轮轴装配面的应力有很大的影响， 摩 擦系数对其应力影响很小。 ( 2 ) 对轮对压装过程中不同过盈量及摩擦系数下最大压装 力及压装压力曲线的结果分析表明： 摩擦系数及过盈量对最大压 装压力及压装压力曲线的影响都比较大， 其中摩擦系数的影响相 对更大。 ( 3 ) 对 HX D ： 轮对在实际过盈量下， 分析得到了压装时合适 的摩擦系数范围为( O ． 0 6 5 ～ 0 ． 1 2 5 ) ， 供实际参考。 参考文献 [ 1 ] 徐国明．影响轮对压装质量的因素[ J ] ．机车车辆工艺， 1 9 8 4 ( 2 ) ： 4 5 _ _ 4 8 ． ( X u G u o — m i n g ． F a c t o r s a f f e c t i n g t h e q u a l i t y o f t h e p r e s s - fi t w h e e l [ J ] ． L o c o m o t i v e ＆R o l l i n g S t o c k T e c h n o l o g y ， 1 9 8 4 ( 2 ) ： 4 5 — 4 8 ． ) [ 2 ] 周兰英． 轮对压装曲线不合格的原因及对策[ J ] ． 铁道机车车辆工人， 2 0 0 1 ( 1 ) ： 2 - 5 ． ( Z h o u L a n- y i n ~T h e r e a s o n f o r t h e p r e s s - f i t c u r v e r o u n d a n d c o u n t e r me - a s u r e s f a i l[ J ] ．R a i l w a y Loc o m o t i v e ＆R o l l i n g S t o c k Wo r k e r s ， 2 0 0 1 ( 1 ) ： 2 — 5 ． ) [ 3 ] 李晓艳．铁路货车轮对压装仿真计算[ J ] ． 铁道车辆， 2 0 1 3 ， 5 1 ( 4 ) ： 8 - 1 1 ． ( L i X i a o — y a n ．T h e w h e e l s e t p r e s s s i m u l a t i o n o n r a i l w a y w a g o n s[ J ] ． R a i l w a y V e h i c l e s ， 2 0 1 3 。 5 1 ( 4 ) ： 8 - 1 1 ． ) [ 4 ] 景彪．关于动车组轮对压装曲线分析判定及建议[ J ] _铁道车辆， 2 0 1 1 ， 4 9 ( 4 ) ： 4 O — 4 2 ． ( j i n g B i a o ．A b o u t E MU w h e e l m o u n t e d o n t h e p r e s s u r e c u rve a n a l y s i s a n d r e c o m m e n d a t i o n s j u d g e m e n t [ J ] ． R a i l w a y V e h i c l e s ， 2 0 1 1 ， 4 9 ( 4 ) ： 4 I — 4 2 ． ) [ 5 ] 张利强． H X N 型机车轮对压装工艺仿真研究[ J ] ． 内燃机车， 2 0 1 3 ( 5 ) ： 1 9 — 2 2 ． ( Z h a n g L i - q i a n g ．T h e p r e s s - fi t S i m u l a t i o n R e s e a r c h i n g o f H XN 5 t y p e o f l o c o m o t i v e [ J ] ．D i e s e l l o c o m o t i v e ， 2 0 1 3 ( 5 ) ： 1 9 — 2 2 ． ) [ 6 ] 陈艳伶． 轮对突悬组装对压装压力曲线的影响[ J ] ． 铁道车辆， 2 0 0 8 ， 4 6 ( 7 ) ： 3 2 — 3 4 ． ( C h e n Y a n- l i n g ．T h e e f f e c t o f Wh e e l o v e r h ang p r e s s - fi t a s s e m b l y o n p r e s s u r e c u rve[ J ] _ R a i l w a y Ve h i c l e s ， 2 0 0 8 ， 4 6 ( 7 ) ： 3 2 - 3 4 ． ) [ 7 ]K w o n S e o k J i n ， S h o j i T e t s u o ． C r a c k d e t e c t i o n i n p r e s s fi t r a i l w a y a x l e u s i n g i n d u c e d c u r r e n t f o c u s i n g p o t e n t i al d r o p t e c h n i q u e [ J ] ． K e y E n g i n e e r i n g Ma t e ria l s ， 2 0 0 4， 2 7 0 - 2 7 3： 1 0 02 —1 0 0 7． [ 8 ]B e n u z z i ． D ， D o n z e l l a ． G ． P r e d i c t i o n o f t h e p r e s s - fi t c u rve i n t h e a s s e m b l y o f a r a i l w a y axl e a n d w h e e l [ J ] J o u r n al o f R a i l a n d R a p i d T r ans i t ， 2 0 0 4 ， 2 1 8 ( 1 ) ： 5 l 一 5 6 ． [ 9 ]T B ／ T 1 7 1 8 — 2 0 0 3 ． 铁道车辆轮对组装技术条件[ s ] ． ( T B ／ r 1 7 1 8 - 2 0 0 3 ．T e c h n i c al c o n d i t i o n s f o r r ai l w a y v e h i c l e w h e e l ass e mb l y [ S ] ． ) ( 上接 第 3页 ) [ 5 ] S j r n D i t t m e r ， K l a u s K t t i g e r ， D e t l e f R i e g e r ， T h o m as M a t s c h u l l a t a n d A r n o Db b e l e r ．As y mme t r i c al p o we r c o n t r o l o f AC- EArs b y s t r u c t u r e - bo rne s o u n d e v a l u a t i o n [ J ] ． A I S T e c h P r o c e e d i n g s ． 2 0 1 0 ( 1 ) ： 9 3 7 — 9 4 5 ． [ 6 ] R o n G e r h a n ， Y u r y K r o t o v and N i c o l d s L u g o ． S e c o n d a r y s y s t e m m e c h a n i c a l r e s o n a n c e i n n e w t a l l — s h e l l S D I — B u t l e r a r c f u r n a c e [ J ] ． A I S T e c h P r o c e e d — i n g s ， 2 0 1 0 ( 1 ) ： 9 4 7 - 9 5 7 ． [ 7 ] N a o t a k a E k i n a g a ． S t r e s s a n a l y s i s o f n i p p l e i n t h e j o i n t o f a r t i fi c i al g r a p h i t e e l e c t r o d e s [ J ] J o u rna l o f t h e C e r a m i c S o c i e t y o f J a p an， 1 9 9 2 ， 1 0 0 ( 6 ) ： 8 0 3 - 8 o 7 ． [ 8 ] 余杨， 赵渭康． 电炉导电横臂技术的发展趋势[ J ] ．冶金设备 ， 1 9 9 9 ， 3 ( 1 1 5 ) ： 4 6 — 4 7 ． ( Y u Yang ， Z h a o We i - k a n g ．T e c h n i c al d e v e l o p me n t t e n d e n c y o f f u r n a c e e l e c t r o d e c o n d u c t i v e a r m [ J 3 ．M e t al l u r - gi c a l E q u i p m e n t ， 1 9 9 9 ， 3 ( 1 1 5 ) ： 4 6 — 4 7 ． ) [ 9 ] 陈文来， 吕红兵， 贾庆远． 电弧炉冶炼中石墨电极折断原因分析[ J ] _ 炭 素技术 ， 2 0 1 0 ， 2 9 ( 3 ) ： 3 6 — 3 8 ． ( C h e n We n - l a i ， L v Ho n g— b i n g ， J i a Q i n g- y u an．C a u s e anal y s i s for t h e f r a c t u r e o f u l t r a — h i g h p o w e r e l e c t r ode[ J ] ．C ar b o n T e c h n i q u e s ， 2 0 1 0 ， 2 9 ( 3 ) ： 3 6 - 3 8) [ 1 0 ] 王勇勤， 刘峻亦， 严兴春． 1 2 ． 5 M N双柱下拉式快锻液压机整机模态分 析[ J ] ．机械设计与制造， 2 0 1 3 ( 2 ) ： 1 6 6 - - 1 6 8 ． ( Wang Y o n g - q i n ， “u J u n - y i ， Y a n X i n g - - c h u rr Mo d e l a n a l y s i s o f t h e w h o l e m a c h i n e o f 1 2 ． 5 MN t w i n c o l u mn s p u l l - - d o w n h i g h s p e e d for gi n g h y d r a u l i c m a c h i n e [ J ] ．M a c h i n e r y D e s i g n & M anu f a c t u r e ． 2 0 1 3 ( 2 ) ： 1 6 6 - 1 6 8 ． )