abaqus帮助文档中轮胎的例子 下载本文

内容发布更新时间 : 2024/3/29 23:16:43星期一 下面是文章的全部内容请认真阅读。

外胎是由胎体、缓冲层(或称带束层)、胎面、胎侧和胎圈组成

1、Bead:胎唇部;2、sidewall:胎侧;3、tread:胎面;4belt:缓冲层;5、carcass:胎体帘布层。

3.1.8 Treadwear simulation using adaptive meshing in ABAQUS/Standard

3.1.8使用自适应网格在Abaqus/Standard中进行轮胎磨损仿真分析

软件:Abaqus/Standard

这个例子在Abaqus/Standard中使用自适应网格技术对稳态滚动的轮胎进行建模。这次分析使用类似“Steady-state rolling analysis of a tire”Section 3.1.2来建立稳态滚动轮胎的接地印迹和状态。接着,进行稳态传输分析来计算和推测持续分析步,在稳态过程中产生一个近似瞬态磨损解。

问题描述和建模

轮胎描述和有限元建模和“Import of asteady-state rolling tire,”Section 3.1.6一样,但是有一些不一样,在这里需要指出。由于这次分析的中心是轮胎磨损,所以胎面建模需要更加精细。另外台面使用线性弹性材料模型来避免超弹性材料在网格自适应过程中不收敛。

图1所示的是轴对称175SR14轮胎的一半模型。橡胶层用CGAX4和 CGAX3单元建模。加强层使用带有rebar层的SFMGAX1单元模拟。橡胶层和加强层之间潜入单元约束。橡胶层的弹性模量为6Mpa,泊松比为0.49。剩下的轮胎部分用超弹性材料模型模拟。多应变能使用系数C10=10^6,C01=0和D1=2*10^8。用来模拟骨架纤维的刚性层和径向成0°,弹性模量为9.87Gpa。压缩系数设置成受拉系数的百分之一。名义应力应变数据用马洛超弹性模型定义材料本构关系。Belt fibers材料的拉伸弹性模量为172.2Gpa。压缩系数设置成拉伸系数的的百分之一。Belt的纤维走向在轴向±20°内。 旋转前面的轴对称一半模型可得到局部三位模型,如图2所示。我们关注轮胎印迹区域的网格。将局部模型镜像后可得到完整的三维模型。

自适应网格在轮胎磨损计算中的局限性

在这个例子中使用自适应网格必须严格遵守以下条件:

1、 圆柱网格不支持自适应网格并且在本例子也没有使用

2、 由于梯度状态变量的变形错误严重,自适应网格使用超弹性材料时表现很差。因此胎面用弹性材料定义

3、 在自适应网格的范围内不能用包含刚性层的嵌入网格。

4、 自适应网格通过网格几何特征来决定自适应网格在自由面光滑的方向,网格几何的特征通常不容易和描述的磨损方向一致。因此,下面将讨论到,通常你需要做额外的工作来明确地描述磨损的方向。

加载

分析分为5个阶段,用轴对称模型开始,以使用symmetric model generation生成的完整三维模型结束。前4阶段和―Steady-state rolling analysis ofa tire,‖ Section 3.1.2和类似。 1、 对称充气:轮胎内部施加200kpa的压力,中间平面使用对称条件。

2、 一半三维接地印迹分析:轴对称模型沿着对称轴旋转。

3、 完整的三维模型接地印迹分析:一半三维模型镜像生成完整的三维模型 4、 稳态滚动:在稳定速度32km/h车速进行完整模型分析,这是轮胎的滚动速度为25rad/s。这些条件符合制动工况,本次分析考虑惯性和迟滞作用。

5、 胎面磨损分析:胎面磨损分析在最后一步进行,本次分析轮胎速度保持为一定值,考虑轮胎表面的磨损,使用损耗的摩擦能来计算磨损。惯性和迟滞同样是本次分析中的考虑因素。本次分析在车速32km/h的情况下,持续进行3.6*10^6秒,轮胎前进32000公里。

最后的分析为磨损分析,来预测磨损或者面消融,根据稳态侧倾轮胎得到评估。我们关注由磨损评估结果得到的轮胎外形的改变;因此我们需要介绍在稳态过程中允许瞬态效果的建模假设。

基本的假设是用当前实时持续的滚动角速度来解释稳态前进分析步。我们认为在任何时候轮胎滚动时轮胎的磨损造成的轮胎外形变化仅仅有很小的效果。因此在整个分析步的每一步稳态的结果都是合适的。有了这些假设,我们就能同时考虑两个不同时间范围的效果:短的轮胎转动时间范围和长的轮胎寿命时间范围。

磨损模型

为了举例说明磨损的过程,假设磨损率是局部接触压力和滑移率的线性函数,进行一个简单的磨损例子。尽管我们能计算这些工程量,由于在稳态移动状态下使用欧拉公式,他们必须应用于胎面流线来模拟轮胎周长磨损。

磨损率计算

磨损模型如下:

q是体积损失量或者磨损量;k是无量纲磨损系数;H是材料硬度;P是接触压力;A是接触面积;Y是接触滑移率。在这里我们可以认为用PAy描述摩擦耗损率。对于轮胎橡胶,我们假设磨损系数k=10^-3,材料硬度H=2GPa。

下面开发的目标是材料的磨损表达式能应用于磨损分析的节点上。首先,考虑用一条带状物围绕着轮胎,带状物的中心用包含胎面花纹的有序节点来定义。这条中心线是以和每个节点联系的辅助面的任意一边作为边界。这样的带状物包含轮胎与路面接触的所有面。我们认为发生在带状物上的磨损是均匀的;因此我们用下式表达整个带状物的磨损率,

其中t是时间,x是当前配置位置。因为我们使用欧拉稳态传输处理,现在表达式可以表示为只依赖于时间的方程,

其中S是沿着流线的位置,T(s)是带状物在S位置的宽度。我们也可以见表达式q写成局部材料衰减率的函数,

在整个带状物的离散化方程的处理的结果相等,得到

其中h是节点的消融速度,A是节点的接触区域。这个方程表明沿着带状物h是不均匀的,推导出的结论是带状物进入和离开接地印迹的宽度是不一样的。然而,因为我们为了维持一个合理的轮胎磨损后的结构,我们假定节点消融的速度是均匀的。设个假定使得如下表达式成立:

再次假定沿着带状物宽度方向的变化可以忽略,既Ti=T,同时认为节点接触区域,则表达式可以简化为没有接地面积的方程:

磨损过程实现

用曲面消融速度的磨损量方程,现在我们可以在稳态移动分析中应用磨损分析了。用户子程序UMESHMOTION用来指定轮胎外表面节点的磨损速度向量。UMESHMOTION用来定义自适应网格约束速度和用来连接自适应网格,每个收敛的增量步之后使用网格光滑技术。通过子程序指定胎面节点的消融速度,自适应网格技术用来调整橡胶层内部的节点来保持好的网格网格形状。

为了积累沿着每个胎面花纹的磨损量,必须在子程序中记录沿着花纹的节点编号。使用一个公共模块变量记录,公共模块记录属于集合NADAPT(图4)的节点和那些在整个模型横截面(0°)的参照点。普通模块变量同时也包括将模型旋转和镜像后节点的编号方式,和参考截面一起,完整地描述了轮胎表面的节点编号。如下的变量需要在外部公共模块中定义: 1、 nStreamlines:轮胎磨损分析中整个参考截面的节点数。 2、 nGenElem:在模型中沿着带状物体网格划分的数目。

3、 nRevOffset:通过*SYMMETRIC MODELGENERATION, REVOLVE指定节点偏移。

4、 nReflOffset:通过*SYMMETRIC MODELGENERATION, REFLECT指定节点偏移。(如果模型不用镜像,这个参数设置为0)。

5、 jslnodes:在参考截面下所有可能发生磨损的节点的节点信息数组。这个数组的大小为(2,nStreamlines)。每个流线的第一个分量是―根节点‖的节点编号(节点a在下面讨论),根节点是指在参考截面上细化流线部分的节点。第二个分量是指提供磨损方向的节点(节点b在下面讨论)。第二个分量仅仅在胎面的拐角处需要,将他设置成等于在参考截面的节点编号来定义磨损的方向。对于不在胎面拐角处的节点,第二个分量设置等于0。那些远离胎面拐角的节点将沿着局部坐标系的3方向进行磨损。

磨损表达式的变量通过函数GETVRN和GETVRMAVGATNODE从分析数据库中获得。

P从变量CSTRESS中获得;y从变量CDISP中获得;由带状物的节点坐标决定,从变量COORD中获得。 磨损运动的方向

磨损速率h是网格约束矢量变量ULOCAL的分量。这个变量通过在局部坐标系ALOCAL中定义的默认网格光滑运动传递到用户子程序中,这个局部坐标系测量出当前节点在曲面的法向。3的方向根据在节点附近网格表面法向的平均值定义为外法向的方向。在绝大部分情况下,这个法向方向对于描述和这个方向相反的消融和节点衰减的磨损结果是足够的。然而,在胎面拐角的区域,这个平均的法向方向不能准确的描述磨损的方向。这种情况下的法向应该像Figure3.1.8–5所示,

它计算如下:假设a是胎面拐角的节点。就有可能辨认出在胎面边上的节点b。在这种情况下,磨损的方向就是矢量ab。因为知道节点a和节点b的坐标,磨损方向就能通过整体坐标系计算得到,也就能转化为局部坐标系(ALOCAL)的方向。