车线动力学(耦合振动3)方案.ppt

?第五节　半车━轨道空间耦合振动模型的有限单元解法 车辆━轨道空间耦合振动求解，涉及到系统中众多的： 耦合关系 非线性因素 时变特性 求解过程十分复杂，因而最有效的方法为有限单元法。 一、模型特点及变量选择 耦合特性： 轮轨垂向振动耦合：轮轨接触力或弹簧 轮轨横向振动耦合：蠕滑力和轮缘力 车辆自由度 车体：沉浮(车体在曲线上受到的未被平衡离心力则以等效离心力和力矩施加于转向架二系悬挂点上) 转向架构架：横移、沉浮、摇头、侧滚及点头 轮对：横移、沉浮、摇头、侧滚 车辆自由度共14个。 轨道自由度 钢轨：8N(每一枕跨为一个单元，每一钢轨节点有垂向位移、绕y轴的转角、横向位移、绕z轴的转角4个未知变量，两股钢轨) 轨枕：3N(沉浮、横移、侧滚) 轨道自由度数：11N 轮轨间设置4个轮缘力未知变量。 (1)当轮缘贴靠钢轨时，轮缘力未知量列入方程，同时引入轮轨横向位移协调条件； (2)而当轮轨不贴靠时，总刚中与轮缘力对应的主元赋大值。 振动模型中总的未知变量数为11N+18个。 变量位置 二、系统动能、势能、虚功及振动方程的组建 １、振动系统的总动能(质量矩阵) ２. 振动系统的总应变能(刚度矩阵与荷载列阵)　　（未包含轮轨接触弹簧的应变能） ３、振动系统中的阻尼力虚功(阻尼矩阵) 4、蠕滑力虚功(荷载列阵) 5、轮缘力及轮缘摩擦力虚功(刚度矩阵) 6、曲线上未被平衡的离心力所作的虚功(荷载列阵) 7、轮轨接触力虚功(刚度矩阵及荷载列阵) 8、重力所作的虚功(荷载列阵) 三、振动方程的求解方法及过程 振动系统的非线性振动方程组： 循环顺序 时间及车轮位置 轮缘与钢轨的贴靠判断 接触力迭代 蠕滑力饱和判断及蠕滑力迭代 　　1、选择积分方法 根据方程组的强非线性和时变特性，选取无条件稳定且无超调现象、对非线性方程有很强适应性的Park积分方法，对微分方程进行求解。 PARK方法的表达式： 代入振动方程可得到： ２、处理非线性因素 振动系统中主要包含非线性接触刚度和非线性蠕滑力两种非线性因素。根据这两种非线性因素不同的特点，采用不同的迭代方法进行处理。 （1） 用直接迭代法处理非线性蠕滑力 因蠕滑力与蠕滑的关系曲线是凸曲线，蠕滑力的直接迭代过程是稳定和收敛的。 对于每一次蠕滑力迭代，振动方程中的荷载列阵变为常量。 （2） 用牛顿迭代法处理非线性接触力 轮轨非线性接触力与位移的关系曲线是凹曲线，直接迭代方法是发散的。因迭代中非线性刚度取为切线刚度，改善迭代收敛性的最简化途径是采用牛顿迭代法。 ３、启步方法 Park方法是线性多步法，不能自动启步，需借助其它方法启步。自动启步的方法很多，可选取最常用的Wilson－θ法。 启步将轮轨接触简化为线性弹簧，可考虑非线性蠕滑力。使启步运算比较简便，避免对于非线性方程不收敛。 振动方程组变为： 4、程序编制 计算程序包含一个主程序和多个子程序（1）钢轨有限和无限单元刚度矩阵；（2）钢轨有限和无限单元质量矩阵；（3）钢轨梁无限单元阻尼矩阵； 　（4）组集总质量矩阵； 　　（5）计算各车轮位置、车轮所在钢轨单元形函数值、车轮下轨面不平顺值；  （6）处理各部弹簧及相关的常数项； （7）处理各部阻尼及相关的常数项； （8）处理曲线上求被平衡的离心力、重力、轮轨间垂向力的水平分力、轨面不平顺； （9）计算轮轨间的蠕滑力。 （10）判断轮缘与钢轨的贴靠状态，处理轮缘力及轮缘摩擦力；  （11）进行蠕滑力饱和判断及修正； （12）用全主元高斯方法求解代数方程； （13）时间为0时计算各单元初始位置； （14）第一、二步Wilson－θ方法启步； （15）从第三步开始用PARK方法处理振动方程。 程序外层结构框图 PARK方法程序框图 程序功能 程序中线路曲线半径、轨道垂向和横向不平顺可任意选择，且内外轨不平顺可以不一样。 运算结果输出为轮轨间蠕滑功、蠕滑力、轮轨垂向力、轮缘力，轨道振动的力、位移、速度和振动加速度。 5、计算方法的参数选择 正式进行计算之前，除对车辆、轨道系统中各有关结构参数要进行研究和选择外， 须对对程序计算结果进行校验， 并进行数值试验，选择与计算方法有关的合适的计算参数，如积分步长、收敛精度、模型长度等。 （1）积分步长 Park方法是无条件稳定的，任意积分步长都能计算出结果。 但步长增加时算法阻尼增大，步长太小耗费机时且增加舍入误差，因此必须根据求解的频率范围选择合适的步长。 根据Park方法算法阻尼与积分步长的关系，对应振幅误差率小于2.5％的最高振型频率为0.1/Δt（Hz），对应振幅误差率小于15％的最高振型频率为0.15