橡胶减振器参数化有限元法优化设计_黄祖宇

橡胶减振器参数化有限元法优化设计黄祖宇焦作市高级技工学校实习工厂(焦作机床厂) (454100)摘 要 为获得与理想减振特性相吻合的橡胶减振器，建立了1/8立体参数化有限元模型。

通过试验确定有限元仿真的材料模型参数并对分析结果进行试验验证。

基于灵敏度分析选择橡胶减振器的外倾角α、半径r 2、高度h 2作为优化变量，以理想载荷-变形特性曲线为优化目标，运用多约束非线性二次规划算法优化得到相应的尺寸值并对橡胶减振器的强度进行了校核。

借助于理想载荷-变形特性曲线并结合有限元法优化得到橡胶减振器的最优外形尺寸，为橡胶元件性能最优化设计提供了新的思路和方法。

关键词 橡胶减振器 参数化建模 有限元 优化设计橡胶减振器是汽车上常见而重要的减振元件，由于橡胶减振器外形结构的复杂性以及它的弹性模量随形状系数而变化的时变特性，使形状不规则的橡胶减振器很难通过解析公式预测其刚度、强度及疲劳寿命，因此早期的橡胶产品开发中大多采用反复试验修正的方法。

20世纪70年代中期以来，随着计算机技术的发展，有限元仿真已成为橡胶件各种性能研究的有力工具：Morman [1]等人用有限元法分析了发动机橡胶减振件的静态和动态特性，Seong Beom Lee [2]借助有限元法预测了汽车橡胶衬套扭转特性，J Pelc [3]研究了充气橡胶轮胎在外力作用下的变形和分层应力特性，Lee [4]用有限元法对发动机橡胶支撑件进行了形状参数优化设计，但未见由理想载荷-变形特性曲线优化得到橡胶减振器外形参数的相关文献。

本文以某工程自卸车橡胶减振器为研究对象，对该减振器进行参数化建模、灵敏度研究以及试验验证，优化得到与理想载荷—变形特性曲线相吻合的工程自卸车橡胶减振器。

1 橡胶大变形数值模拟基本方程橡胶材料在外载作用下的大变形，同时具有材料非线性、几何大变形非线性以及接触非线性的特点，解析求解困难，目前一般借助有限元方法对其研究，其基本求解过程如下：根据虚功原理，外力在虚位移上所作的虚功等于内力在虚应变上所作的虚功，将虚位移原理应用于橡胶元件初始构形的每一个单元： 0()()ee TT e e T e V a B SdV W a F δδ≡=∫(1)式中 B ——应变-位移矩阵S ——克希荷夫应力张量V 0e ——单元初始构形下的体积 F e ——作用在单元上的外力。

对于橡胶类材料，不可压缩条件在每个单元上平均地满足，因此有如下补充方程03(1)0eV I dV −=∫(2)上式中3I 是格林应变张量的第三不变量。

按照整体拉格朗日描述，由式(1) ～ (2)可得t t +Δ时刻第n+1次单元平衡方程的迭代近似解如式(3) ～ (4)所示。

000()*()**ee e e t te T T n n n n V K a T p F B S dV λλ+ΔΔ+Δ=−∫(3)00()*ee T n n n V T B L dVλλ=∫(4) 00000()[***()*]e e T T T n n n T n n V e e en Ln nK G M G B D B dVK K K λλλλσ=+=++∫(5)由式(5)可以看出，橡胶变形时切线刚度矩阵()e T n K 由三部分组成：一是小位移矩阵0enK ；二是大位移矩阵eLnK ；三是初始应力矩阵和几何刚度矩阵enK σ。

将单元方程(3) ～ (5)组装成橡胶元件整体有限元方程，得大变形状态下橡胶元件的整体求解方程如式(6)所示：0T T n n nK T a R T p r Δ⎡⎤⎧⎫⎧⎫=−⎨⎬⎨⎬⎢⎥Δ⎣⎦⎩⎭⎩⎭ (6) 适当约束刚体位移，由方程组(6)可解出n a Δ 和n p Δ，在每一个载荷步内重复应用修正的Newton-Raphson 方法直至全部载荷收敛即可求解出橡胶元件大变形时的位移值。

2 参数化建模与试验验证由橡胶大变形求解公式(3) ～ (6)可见对橡胶元件尤其是不规则橡胶件直接进行数值计算相当复杂而难以完成有效设计。

本文以ANSYS 为平台，采用有限元法对图1所示的橡胶减振器进行数值模拟和优化，降低了橡胶元件大变形分析的难度，为 橡胶元件结构分析提供切实可行的方法。

整个橡胶减振器组合体由上钢套1、橡胶体2、下钢套3组成，如图1所示。

完整的橡胶减振器参数化有限元优化过程包括合理确定材料模型、参数化建模、灵敏度分析以及优化求解等几部分。

1.上钢套2. 橡胶体3.下钢套图1 橡胶减振器组合体2.1 合理选择橡胶材料模型橡胶材料的本构模型主要分为现象学描述的多项式模型以及Ogden 模型[5]和基于热力学统计描述的Arruda-Boyce 模型两大类。

相对低阶多项式模型，Ogden 模型除了在非常大的应变下(λ>7)才会出现显著的偏差，与试验的吻合程度较好，高阶多项式模型尽管与试验吻合一致，但由于累计误差的原因将之用于有限元计算时会导致收敛困难。

因此，Ogden 模型具有较广泛的适用性，Ogden 本构模型的应变能表达式为：()212321121(3)1i i i N Ni i i i i iU J D αααμλλλα===++−+−∑∑ (7)Ogden 应变能以三个主伸长率λ1、λ2、λ3为变量，其应变能函数的第一部分只与12I I ,有关，需要由试验确定的参数有µi 和αi 。

应力单轴拉伸应变图2 单轴拉伸试验数据拟合应力双轴拉伸应变图3 双轴拉伸试验数据拟合图2和图3分别是单轴拉伸和双轴拉伸的拟合结果与试验结果的对比图。

由两图可以看出，单轴拉伸试验与拟合结果吻合较好，由于双轴拉伸的试验数据由单轴压缩转化而来，因此其拟合效果稍差。

从图上看Ogden 三阶模型与橡胶试件的应力—应变曲线更为吻合且其残差为所有拟合形式中最小者，因此Ogden 三阶六参数模型是最能描述本文橡胶减振器的材料模型。

采用Marquard-Levenberg 法[6]拟合得到的橡胶材料参数如下：µ1=0.002 1185，a 1=8.314，µ2=1.392 41，a 2=1.519 6，µ3=1.434 1，a 3=1.519 6。

2.2 参数化模型与试验验证为提高运算效率并真实模拟橡胶减振器的变形状态，取橡胶减振器的一半建模并增加一刚性支撑板，橡胶减振器1/8立体有限元模型截面尺寸如图4所示，图中h 1表示减振器橡胶体高度；α是外部圆锥台与支撑钢板之间的外倾角；r 2是上部圆锥台的挖空半径；h 2表示上部圆锥台的挖空高度；r 3表示橡胶减振器中部过渡圆角半径。

各尺寸参数的基准值分别为h 2=39 mm ，r 2=38.5 mm ， r 3=10 mm ，α=58°，1/8立体有限元模型由图4所示的参数化平面模型旋转而成。

图4 橡胶减振器半剖面图本文采用3-D 8节点位移-压力杂交单元Solid185模拟橡胶大变形，材料模型选用由试验确定的Ogden 三阶模型，划分完网格后的立体有限元模型共2 778个节点、2 248个单元，如图5所示。

图5 橡胶减振器参数化有限元模型为正确模拟橡胶体与钢套之间的接触状态，采用拉格朗日和罚函数混合法进行接触分析，并在橡胶体与钢套以及橡胶体与下支撑板之间分别定义刚-柔接触面，如图6所示。

图7是橡胶减振器垂向变形特性仿真与试验对比，从图上可以看出有限元仿真与橡胶减振器试验结果基本一致，变化趋势一致，只是在很大的压力载荷下相对误差稍大，但基本都小于工程应用允许的5%，有限元模拟有效预测了橡胶减振器的非线性变形。

因此，有限元法能预测并实现橡胶减振器的优化设计。

垂向载荷F V /k N垂向变形D V /mm图6 橡胶减振器垂向载荷-变形曲线仿真与试验对比2.3 灵敏度分析为保持结构的稳定性并保证有限元优化的效率，本文只对α、r 2、h 2、r 3这四个参数对载荷-变形尺寸的灵敏度进行研究。

橡胶减振器参数一阶百分比灵敏度计算公式如式(8)所示：()/r r r r d x s x x Δ=Δ (r =1…7) (8)式中 x r ——第r 个设计设计参数Δx r —— 第r 个参数的变化值 x r 0—————各参数的基准值Δd (x r )——橡胶减振器垂向位移变化值 S r ——橡胶悬架第r 个参数的灵敏度。

据式(8)计算得到的橡胶减振器空载工况下参数灵敏度如图7所示满载工况下得到的灵敏度曲线与此基本一致。

灵敏度值参数相对基准值的百分比(%)图7 空载下参数灵敏度分析由图7可以看出，α、h 2、r 2的灵敏度值都为正值，表明上述参数的增加能使垂向位移迅速增大。

外倾角α的改变对橡胶减振器垂向变形量的影响最为明显，挖空高度h 2、挖空半径r 2的影响也较为显著，但减振器过渡圆角半径r 3的影响很不明显。

因此，橡胶减振器优化设计时以参数α、h 2、r 2作为系统的设计变量，不考虑过渡圆角半径r 3的影响。

3 橡胶减振器有限元优化研究本文利用工程自卸车整车数学模型优化得到了橡胶减振器最佳载荷-变形特性曲线(另文撰述)。

根据该曲线，采用有限元法优化得到减振器的外形参数，从而设计出具有理想变刚度特性的橡胶减振器。

3.1橡胶减振器优化数学模型表1是根据遗传算法优化得到的理想载荷-变形特性，本文利用有限元一阶优化方法和橡胶减振器参数模型对挖空半径r 2，挖空高度h 2以及橡胶减振器与支撑钢板的外倾角α三个参数进行优化。

优化目标定义为各载荷下有限元仿真结果与橡胶减振器优化结果之差的平方和最小，计算公式如(9)所示。

表1 橡胶减振器理想载荷变形特性曲线 载荷/kN 5 10 15 2025变形量/mm26.1 44.5 57.8 68.276.8 h (1/2变形量)/mm 13.05 22.25 28.9 34.138.4载荷/kN 30 35 40 45—变形量/mm84.3 90.8 96.7 102.— h (1/2变形量)/mm 42.15 45.4 48.35 51.0—2122m a x 2m in 22m a x 2m in 22m a x m in m a x m in ()()[]..[]n i i i T f X x h X r h s t r r r h h h ασσααα=⎧=−⎪⎪⎪=⎪≤⎨⎪≤≤⎪≤≤⎪⎪≤≤⎩∑ (9) 式中x i ——1/8立体有限元模型各载荷下有限元仿真结果；h —— 橡胶减振器一半变形量，如表1所示；σmax —— 满载载荷下橡胶减振器橡胶体的V on mises 应力，是性能约束函数； i —— 载荷步(i =1…9)；X = [r 2 h 2 α]T 表示待优化的设计变量，是边界约束函数。

3.2 优化算法优化过程中采用了多约束非线性二次规划法，假设经过k 次优化迭代，得到设计点xk ，现在从该点出发继续搜索，最终逼近问题的最优点。