利用有限元法分轮胎胎面的弹塑性场.docVIP

利用有限元法分析轮胎胎面的弹塑性S K Deb Nath 日本东北大学纳米力学博士后研究员数值及实验研究利用得到轮胎弹塑性性能弹塑性应力应变关系三维实体单元Von Mises应力应变关系分别独立描述轮胎接触表面的弹性和塑性区域在不同区域Von Mises应力和不同细致地分析确定Von Mises应力的分布和不同的应力分量轮胎关键轮胎是由两种不同的材料,橡胶复合材料和各向同性橡胶根据其真正的几何学详细的弹塑性分析制造是如果应力分析考虑了轮胎胎面单元进行的和的材料性能这是不可能得到其精确的弹塑性行为所以确定真正的应力应变场、各向同性和各向异性的部分必须使用不同的材料特性和相关的边界条件分析一些研究人员假定轮胎橡胶的弹性特性如天然橡胶但在实际情况下,轮胎橡胶材料性质完全不同天然橡胶Nath获得轮胎橡胶弹塑性性能试验表明,压缩模量和张力系数是完全不同的弹塑性性能完全不同天然橡胶。Ahmed et al.定义了最优形状的轮胎由于侧向滑移考虑采用有限差分法二维平面应力的概念保证了轮胎接触表面最小磨损王,丹尼尔和黄Moire方法分析应力应变的试验方法，用多次翻新的轮胎胎肩和胎面做比较，其结果通过有限元方法分析。在橡胶复合材料结构载荷作用下的粘弹性特性一些作者有限元方法可以预测模态振型和固有频率装载轮胎形状轮胎钢帘线制品前大多数制造商使用专有软件无数的论文也处理涉及静态或者Padovan和Pardilok注重了黏弹性瞬态及不同速度稳态滚动接触最近值得Padova注意的是可以认定为轮胎各种驻波不同负载一些学者提出制定有限元模型计算应力状态由于因一辆车自重造成的轮胎接触平面负荷下刚性路面提供一个有用的调查考虑各向同性材料和相关的边界条件分析轮胎的弹塑性行为载荷作用下的弹塑性是的弹塑性行为这种分析将的指导生产轮胎的模型轮胎几何形状图1。轮胎模拟利用MARC有限元软件轮胎接触表面假定是自由的,因为侧向位移之间的摩擦力轮胎和道路有关边界条件显示在表1。实验获得的弹塑性应力应变关系图弹性模量是0.7 MPa、泊松比为0.45。塑性应变在接触表面的假定3MPa。 表1: 轮胎物理边界条件图1:图2:卡车轮胎橡胶压缩应力-应变关系图:卡车轮胎橡胶压缩作用下弹性模量的确定拉格朗日为解决以下轮胎接触问题的拉格朗日程序使用。在完全拉格朗日方法、平衡可以虚功原理表示为这里是对称第二Piola-Kirchhoff应力张量是Green-Lagrange应变是在参考是在参考牵引矢量是虚拟位移在原来的进行而进一步分解成线性的增量的线性部分增量是更新拉格朗日线性化方程(1)Δu和η分别是虚拟位移 ?s指?的对称部分，代表了当前梯度算子结果和讨论三维实体单元模型模拟轮胎的问题确定轮胎的弹塑性行为了解实际的弹塑性性能必要。分析弹塑性行为,我们需要得到弹性模量和塑性应力应变关系塑应力-应变关系方程转化为一个所有的弹塑性材料特性是用来解决当前轮胎问题。MARC有限元软件用来解决当前轮胎的问题从分析真正的弹塑性应力应变场的发现将有助于设计的轮胎图5描述了应用竖向荷载垂直位移的一个函数在两个不同节点的接触表面。图6中的阐明了Von Mises应力从这关系弹性和塑性区。Von Mises图7说明Von Mises应力区z / = -0.5作为x / b。Von Mises应力x / b比增加或减少Von Mises应力在两端，x / b=1或0处。在0≤x / b ≤0.9的区域，在x / b=0.5处观察到最小的应力。并且在第1个节点处压力下降的速率显著高于第2个节点处，对比第2个节点观察到相反的现象。在第3个节点0.15≤x / b≤0.85处，Von Mises应力分布是几乎不变的，但是在其他区域，应力值向两端减小。在4到6节点处应力分布几乎相同。在第1、2和3第节显示最高应力z /a = -0.5处相当于在轮胎不同部位x/b的功能。在第1节点的两端，正应力σxx集中，这是超过双外加压力。在0.15x/b0.85区域，节点1,2和3作为x/b显示常压。在3个节点中，第1个显示最高压力，且第2个显示最低压力，但是观察到在第2和第3之间有轻微的变化，其余部分是零压力。图9显示的正应力分量，σyy在z /a = -0.5处相当于x/b。在两端（x/b=0或1），而第1节点显示的最高压力几乎相当于双外加压力，节点2和3显示处微不足道的压力。在其他区域，三个节点中第3个节点显示出最大的压力，节点1为最低压力，并且节点2的应力值在他们之间。节点4，5和6显示出持续的压力，相当于在接触表面的应力。图10中的正压力分布，σzz在z /a = -0.5处相当于在轮胎不同部位x/b的功能。在第1节点的压力是双外加压力。在所有的部分，最大的压力是在第1节点观察到的。图11显示切应力，在z/a=-0.5