《车线动力学耦合振动课件.ppt》由会员分享，可在线阅读，更多相关《车线动力学耦合振动课件.ppt（52页珍藏版）》请在taowenge.com淘文阁网|工程机械CAD图纸|机械工程制图|CAD装配图下载|SolidWorks_CaTia_CAD_UG_PROE_设计图分享下载上搜索。
1、车线动力学耦合振动第1页，此课件共52页哦n车辆轨道空间耦合振动求解，涉及到系统中众多的：耦合关系 非线性因素 时变特性n求解过程十分复杂，因而最有效的方法为有限单元法。第2页，此课件共52页哦一、模型特点及变量选择一、模型特点及变量选择第3页，此课件共52页哦半车轨道的“垂向横向非线性空间耦合振动时变模型”Mc,Jc第4页，此课件共52页哦n耦合特性：轮轨垂向振动耦合：轮轨接触力或弹簧 轮轨横向振动耦合：蠕滑力和轮缘力n非线性特性(坐标原点位于系统不受任何力时的质心位置)：轮轨间垂向力或非线性赫兹接触弹簧 轮轨间非线性蠕滑力n时变特性：半车模型在轨道移动，车辆每前进一步，振动方程均需重新组建
2、。第5页，此课件共52页哦车辆自由度n车体：沉浮(车体在曲线上受到的未被平衡离心力则以等效离心力和力矩施加于转向架二系悬挂点上)n转向架构架：横移、沉浮、摇头、侧滚及点头n轮对：横移、沉浮、摇头、侧滚n车辆自由度共14个。第6页，此课件共52页哦轨道自由度n钢轨：8N(每一枕跨为一个单元，每一钢轨节点有垂向位移、绕y轴的转角、横向位移、绕z轴的转角4个未知变量，两股钢轨)n轨枕：3N(沉浮、横移、侧滚)n轨道自由度数：11N第7页，此课件共52页哦n轮轨间设置4个轮缘力未知变量。(1)当轮缘贴靠钢轨时，轮缘力未知量列入方程，同时引入轮轨横向位移协调条件；(2)而当轮轨不贴靠时，总刚中与轮缘力对
3、应的主元赋大值。n振动模型中总的未知变量数为11N+18个。第8页，此课件共52页哦变量位置部件位置曲线外轨(4N)14N曲线内轨(4N)4N+18N轨枕(3N)8N+111N第1轮对(4)11N+111N+4第2轮对(4)11N+511N+8构架(5)11N+911N+13车体(1)11N+14轮缘力(4)11N+1511N+18第9页，此课件共52页哦二、系统动能、势能、虚功二、系统动能、势能、虚功及振动方程的组建及振动方程的组建第10页，此课件共52页哦、振动系统的总动能、振动系统的总动能(质量矩阵质量矩阵)左右钢轨轨枕轮对构架车体第11页，此课件共52页哦.振动系统的总应变能振动系统的
4、总应变能(刚度矩阵与荷载列阵刚度矩阵与荷载列阵)（未包含轮轨接触弹簧的应变能）外轨弯曲内轨弯曲外轨扣件横向刚度内轨扣件横向刚度第12页，此课件共52页哦外轨扣件垂向刚度内轨扣件垂向刚度道床横向刚度道床外侧垂向弹簧刚度道床内侧垂向弹簧刚度第13页，此课件共52页哦一系悬挂横向刚度一系悬挂纵向刚度曲线半径的影响，进入荷载列阵第14页，此课件共52页哦一系悬挂垂向刚度二系悬挂摇头角刚度二系悬挂侧滚角刚度二系悬挂点头角刚度二系悬挂垂向刚度二系悬挂横向刚度第15页，此课件共52页哦、振动系统中的阻尼力虚功、振动系统中的阻尼力虚功(阻尼矩阵阻尼矩阵)外轨扣件横向阻尼内轨扣件横向阻尼外轨扣件垂向阻尼内轨扣件
5、垂向阻尼第16页，此课件共52页哦道床横向阻尼道床外侧垂向阻尼道床内侧垂向阻尼一系悬挂垂向阻尼第17页，此课件共52页哦二系悬挂垂向阻尼二系悬挂摇头阻尼二系悬挂侧滚阻尼二系悬挂点头阻尼二系悬挂横向阻尼第18页，此课件共52页哦4 4、蠕滑力虚功、蠕滑力虚功(荷载列阵荷载列阵)纵向蠕滑力横向蠕滑力自旋蠕滑力矩第19页，此课件共52页哦 5 5、轮缘力及轮缘摩擦力虚功、轮缘力及轮缘摩擦力虚功(刚度矩阵刚度矩阵)轮缘横向力轮缘摩擦力第20页，此课件共52页哦6 6、曲线上未被平衡的离心力所作的虚功、曲线上未被平衡的离心力所作的虚功(荷载列阵荷载列阵)轮对构架车体第21页，此课件共52页哦7 7、轮轨
6、接触力虚功、轮轨接触力虚功(刚度矩阵及荷载列阵)处理成弹簧的应变能第22页，此课件共52页哦处理成力的虚功垂向力垂向力的水平分力第23页，此课件共52页哦第24页，此课件共52页哦8 8、重力所作的虚功、重力所作的虚功(荷载列阵)第25页，此课件共52页哦三、振动方程的求解方法及过程三、振动方程的求解方法及过程第26页，此课件共52页哦n振动系统的非线性振动方程组：不同时刻车轮处于轨道的不同位置，钢轨不平顺值也不相同，荷载列阵发生变化。荷载列阵中包含蠕滑力，隐含振动位移、速度。轮轨接触力是波动的，轮轨接触弹簧瞬时刚度依赖于接触力。第27页，此课件共52页哦n循环顺序 时间及车轮位置 轮缘与钢轨
7、的贴靠判断 接触力迭代 蠕滑力饱和判断及蠕滑力迭代第28页，此课件共52页哦1 1、选择积分方法、选择积分方法n根据方程组的强非线性和时变特性，选取无条件稳定且无超调现象、对非线性方程有很强适应性的Park积分方法，对微分方程进行求解。第29页，此课件共52页哦nPARK方法的表达式：第30页，此课件共52页哦n代入振动方程可得到：第31页，此课件共52页哦、处理非线性因素、处理非线性因素n振动系统中主要包含非线性接触刚度和非线性蠕滑力两种非线性因素。根据这两种非线性因素不同的特点，采用不同的迭代方法进行处理。第32页，此课件共52页哦（1 1）用直接迭代法处理非线性蠕滑力用直接迭代法处理非线
8、性蠕滑力n因蠕滑力与蠕滑的关系曲线是凸曲线，蠕滑力的直接迭代过程是稳定和收敛的。0第33页，此课件共52页哦n对于每一次蠕滑力迭代，振动方程中的荷载列阵变为常量。第34页，此课件共52页哦（2 2）用牛顿迭代法处理非线性接触力用牛顿迭代法处理非线性接触力n轮轨非线性接触力与位移的关系曲线是凹曲线，直接迭代方法是发散的。因迭代中非线性刚度取为切线刚度，改善迭代收敛性的最简化途径是采用牛顿迭代法。第35页，此课件共52页哦第36页，此课件共52页哦、启步方法、启步方法nPark方法是线性多步法，不能自动启步，需借助其它方法启步。自动启步的方法很多，可选取最常用的Wilson法。第37页，此课件共5
9、2页哦n启步将轮轨接触简化为线性弹簧，可考虑非线性蠕滑力。使启步运算比较简便，避免对于非线性方程不收敛。第38页，此课件共52页哦n振动方程组变为：第39页，此课件共52页哦4 4、程序编制、程序编制n计算程序包含一个主程序和多个子程序（1）钢轨有限和无限单元刚度矩阵；（2）钢轨有限和无限单元质量矩阵；（3）钢轨梁无限单元阻尼矩阵；（4）组集总质量矩阵；（5）计算各车轮位置、车轮所在钢轨单元形函数值、车轮下轨面不平顺值；第40页，此课件共52页哦（6）处理各部弹簧及相关的常数项；（7）处理各部阻尼及相关的常数项；（8）处理曲线上求被平衡的离心力、重力、轮轨间垂向力的水平分力、轨面不平顺；（9）
10、计算轮轨间的蠕滑力。（10）判断轮缘与钢轨的贴靠状态，处理轮缘力及轮缘摩擦力；第41页，此课件共52页哦（11）进行蠕滑力饱和判断及修正；（12）用全主元高斯方法求解代数方程；（13）时间为0时计算各单元初始位置；（14）第一、二步Wilson方法启步；（15）从第三步开始用PARK方法处理振动方程。第42页，此课件共52页哦n程序外层结构框图第43页，此课件共52页哦nPARK方法程序框图第44页，此课件共52页哦程序功能n程序中线路曲线半径、轨道垂向和横向不平顺可任意选择，且内外轨不平顺可以不一样。n运算结果输出为轮轨间蠕滑功、蠕滑力、轮轨垂向力、轮缘力，轨道振动的力、位移、速度和振动加速
11、度。第45页，此课件共52页哦5 5、计算方法的参数选择、计算方法的参数选择n正式进行计算之前，除对车辆、轨道系统中各有关结构参数要进行研究和选择外，n须对对程序计算结果进行校验，n并进行数值试验，选择与计算方法有关的合适的计算参数，如积分步长、收敛精度、模型长度等。第46页，此课件共52页哦（1 1）积分步长）积分步长nPark方法是无条件稳定的，任意积分步长都能计算出结果。n但步长增加时算法阻尼增大，步长太小耗费机时且增加舍入误差，因此必须根据求解的频率范围选择合适的步长。第47页，此课件共52页哦n根据Park方法算法阻尼与积分步长的关系，对应振幅误差率小于2.5的最高振型频率为0.1/
12、t（Hz），对应振幅误差率小于15的最高振型频率为0.15/t（Hz）。n如积分步长取为0.001秒，则对应100Hz振型的振幅误差率为2.5，对应150Hz的振型的振幅误率为15。n如要确保400Hz以内的振幅误差率小于2.5，则积分步长应小于0.00025秒。第48页，此课件共52页哦n实际计算中，为节省机时，根据不平顺的长度和列车速度进行时间步长数值试验。步长从0.001秒开始反复折半，计算出23组结果，当前后两组结果的误差在510以内时，认为步长是合适的。第49页，此课件共52页哦（2 2）迭代误差限）迭代误差限n迭代误差限依据节省机时且又满足最终结果精度的原则而确定。n在轮轨非线性接触力的迭代过程中，首先控制广义位移列阵（含位移、转角）中各元素的精度，使位移增量列阵的2-范数小于10-8，即：n相当于位移、转角的精度为10-610-7，力的精度为10-410-5。第50页，此课件共52页哦n当位移增量满足精度时，判断轮轨接触力的精度，要求前后两次接触力值的相对误差小于1，而后轮轨接触力按常力处理。第51页，此课件共52页哦n蠕滑力迭代循环中，控制轮轨间纵向和横向蠕滑力的精度，如前后两次蠕滑力值的相对误差小于1。第52页，此课件共52页哦