刘 胜,李 鹏,杨翊仁
(西南交通大学力学与工程学院,成都 610031)
图1为一种新型水力缓冲器结构模型,主要由套筒、活塞、碰撞部件、缓冲头和缓冲弹簧组成。套筒与活塞杆之间安装有预紧缓冲弹簧,使得相对静止状态下套筒底部与活塞头部紧密贴合。运动部件安装在缓冲器上部,并置于静水中。工作时,运动部件随缓冲器一起从一定高度落下。当进入下落行程末段时,缓冲头与碰撞部件接触,并以一定速度和加速度冲击碰撞部件。冲击过程中,缓冲头压缩碰撞部件弹性件,同时自身速度不断减小,直至弹性件被完全压缩,而后缓冲头由刚性件提供支撑;另一方面,冲击时缓冲头压缩活塞杆,从而挤压套筒内的冷却剂和缓冲弹簧,起到对运动部件缓冲的效果。运动部件完全静止后,缓冲弹簧推动活塞杆和缓冲头复位至安装初始状态。
本文依据水力缓冲器的结构和工作原理,建立下落过程和缓冲过程的结构控制方程和流体控制方程,计算缓冲过程中缓冲头对碰撞部件的冲击力、套筒内部压强、套筒和活塞位移等参数的动态变化,分析各节流孔对冲击力和压强影响,对缓冲器结构进行优化。
图1 水力缓冲器结构
为建立缓冲过程的控制方程,首先给出以下基本假设:①在下落过程中,套筒与活塞为无相对运动,可视为一个整体;②各部件之间的碰撞为完全弹性碰撞,忽略碰撞阻尼,碰撞等效为弹簧作用,组件之间的碰撞刚度为组件拉压刚度的串、并联值[1-2];③ 流体不可压缩,其能量损失仅为局部能量损失;④仅考虑外部流场对缓冲器的流体形状阻力。
缓冲器在下落和缓冲过程中的结构示意图如图2所示。M1,M2分别为套筒和活塞杆的质量;Mf1,Mf2分别为二者排开水的质量;Cf,Cs分别为黏性阻力系数和形状阻力系数;D1,D2,H1,H2分别为套筒与活塞杆的外直径和高度;Sht为缓冲头截面积;K,F0分别为缓冲弹簧刚度和预紧力;K1,K2,K3分别为组件各部分的碰撞刚度;H为活塞有效工作高度;He为碰撞部件弹性件最大压缩行程;P1,P'1,P2,P'2分别为作用与套筒和活塞上的流体压力。
图2 缓冲器结构简化示意图
下落过程结构的运动方程可写为[3]:
缓冲过程结构控制方程为
其中:
缓冲头与碰撞部件之间的冲击力为
其中δ为Dirichlet函数。
图3给出了缓冲过程中流量的计算模型。套筒与活塞之间的相对运动会挤出套筒主腔内的流体或吸入外部流体,因此各节流孔的流体流向及流量与此相对运动相关。
图3 流量计算示意图
首先假定当主腔内流体被挤出,且挤出流体均是从活塞头部流出,而后经过各节流孔流出套筒;当流体被吸入主腔时,认为所有流体均流向活塞头部。各节流孔流量分别记为αQ,βQ,γQ,ηQ,ξQ。忽略非定常项和位置水头对总水头的影响,可由Bernoulli方程[4]沿流线建立能量方程:
流体从节流孔流出主腔时
流体从节流孔流进主腔时
其中:i=1,2,…,5,分别对应5个不同节流孔;p,U分别为活塞头部压强与流速;pi,Ui分别为流出i号节流孔后的流体压强与流速;∑hm,∑h'm分别为沿流线的所有局部能量损失,流出及流进二者不同。由于套筒侧壁开有节流孔,各条流线上的能量方程与下落过程中活塞头与各节流孔的相对位置相关。
忽略外部流场的流动,外部流场满足
流动中的流量方程为
其中:λi= ±1,i=1,2,3,4,5为流向迭代指标;ΔV为副腔体积变化;λV=±1为副腔体积增大和缩小指标。
联立式(5)~(8),可得流量控制方程:
流量控制方程中,活塞头部流体流速应等于套筒与活塞的相对运动流速,即
采用 Newton迭代方法[5]求解式(9)。由式(9)获得节流孔流量后,代入式(5)或(6)即可得到缓冲器内部压强分布。
采用Newmark积分格式[6]对结构控制方程求解。基于Fortran程序[7]编制整个求解过程数值计算程序。图4给出了程序计算流程。图5给出了当下落高度为1.75 m时冲击力和活塞头部压力的时程曲线。显然,在第1次冲击之后,冲击力和压强都迅速降低。此外,当活塞经过套筒侧壁节流孔时,压强会出现显著突变。
图4 缓冲分析的数值模拟流程
图6给出了不同下落高度时最大冲击力和最大压强的变化曲线,可见随着下落高度的增加,最大冲击力和最大压强随之增大。
图7给出了节流孔尺寸对最大冲击力的影响曲线,节流孔编号与图3中对应。由图7可以看出:随着#4、5节流孔直径的增大,最大冲击力减小,#5节流孔直径的影响更加明显。#1节流孔直径较小时对最大冲击力影响也很小,当直径大于3 mm后,冲击力随之先增大后减小。
图5 冲击力、压强时程曲线
图6 最大冲击力和最大压强随下落高度变化
图7 最大冲击力随节流孔直径变化
针对一新型水力缓冲器建立了理论模型。该模型包括结构控制方程、流量控制方程和结构-流量耦合控制方程。并根据该模型,进行了缓冲器的程序实现与数值模拟,得到以下结论:
1)最大冲击力与最大压强出现在冲击开始的短暂时间内,之后,二者均显著下降。
2)随着初始下落高度的增加,最大冲击力和最大压强随着增大。
3)套筒顶部节流孔在尺寸达到一定程度后对缓冲影响较大,在侧壁节流孔中,上方节流孔对缓冲影响更大。
[1]孙训方,方孝淑,关来泰.材料力学(I)[M].4版.北京:高等教育出版社,2002.
[2]倪振华.振动力学[M].西安:西安交通大学,1989.
[3]金栋平,胡海岩.碰撞振动与控制[M].北京:科技出版社,2005.
[4]张鸣远.流体力学[M].北京:高等教育出版社,2010.
[5]John H M,Kurtis D F.数值方法(MATLAB 版)[M].4版.陈渝,钱方,译.北京:电子工业出版社,2005.
[6]王勖成.有限单元法[M].北京:清华大学出版社,2003.
[7]彭国伦.Fortran95程序设计[M].北京:中国电力出版社,2006.