首页/文章/ 详情

基于有限元与预变形的橡胶脚垫动态刚度建模研究

7月前浏览343

在汽车工程中,结构传播噪声(structure‑borne noise)是座舱内重要的噪声来源,而发动机振动通过机脚垫传递至车身是其中关键路径之一。为了有效抑制高频振动,提高整车 NVH 性能,必须准确研究和建模橡胶机脚垫的动态传递刚度(dynamic transfer stiffness)特性。

本文介绍一种结合有限元分析(FEA)的橡胶机脚垫动力学建模方法,重点讨论**预压缩(static preload)**对机脚垫动态刚度的影响。


1. 引言

橡胶机脚垫承担发动机支撑与隔振的双重任务,其主要功能是阻隔发动机高频振动向车身的传递。橡胶材料本身具有高度非线性与不可压缩性,因此为了评估机脚垫在宽频范围内的隔振性能,必须建立一个足够真实的数值模型以模拟其在预载荷条件下的动态行为。

研究采用了商用有限元软件 Abaqus,通过非线性静态分析与叠加线性谐波分析,构建一套用于预测机脚垫动态刚度随频率变化规律的方法。


2. 方法与分析流程

2.1 橡胶材料特性

橡胶具有以下典型力学特征:

大变形特性

:变形过程中应变可达数十甚至上百个百分点;

非线性本构关系

:载荷—位移曲线明显非线性;

近似不可压缩

:体积变化极小;

黏弹性

:动态加载下表现出频率相关性。

因此传统线性材料模型难以描述其真实行为,需要采用非线性材料模型并结合静态预变形。


2.2 分析步骤

分析过程分为两个主要阶段:

由于发动机重量会使机脚垫产生显著的轴向预压缩,因此需要首先计算这一初始变形状态(图1)。该步骤使用非线性求解器,采用真实橡胶材料模型,得到受载后的稳定形状。

图1:机脚垫在轴向载荷下的静态预变形 u₁。

在静态预变形基础上,叠加小幅振动并进行线性化谐波分析(图2)。假设发动机运行产生的振动幅值相对较小,可将其视为在预变形平衡状态周围的线性扰动。


3. 动态传递刚度的计算

机脚垫的动态行为可用一个分块动态传递刚度矩阵描述,其中包含三个平动方向与三个旋转方向的耦合项。

对于轴向方向,动态传递刚度定义为:

k12=F2/u1

其中:

F2:作用在机脚垫另一侧的反力;

u1:在上端施加的轴向振动位移。

图3与图4分别展示了不同预变形水平下,机脚垫轴向动态传递刚度的幅值与相位随频率变化的典型结果。

图3:不同预变形水平下轴向动态传递刚度的幅值



图4:不同预变形水平下动态传递刚度的相位曲线

从图中可观察到:预压缩会显著改变刚度曲线,特别是在中高频段。


4. 结果与讨论

研究发现:

预压缩显著影响机脚垫的动态刚度;

某些频率范围内,随着预变形程度增加,动态传递刚度可能急剧变化;

动态刚度幅值的增加会使机脚垫隔振性能下降(因为更高刚度 → 更强传递路径),从而导致结构传播噪声加剧;

相位变化反映了机脚垫在不同频率下的能量传递效率与滞后行为,对动力系统 NVH 调校具有工程意义。

这意味着在实际整车 NVH 仿真与标定过程中,忽略静态预载荷会导致动态刚度评估显著偏差,从而影响隔振设计。


5. 结论

本文基于有限元建模展示了橡胶机脚垫在静态预变形与动态加载共同作用下的动态刚度特征,主要结论如下:

  • 静态预压缩对动态刚度影响显著,必须在模型中正确体现。

  • 机脚垫在宽频范围内表现出明显的频率相关性,尤其是高频段。

  • 随着预压缩程度增加,某些频率范围的动态刚度会明显上升,使隔振性能下降。

基于预变形 + 小幅线性扰动的建模方法适用于机脚垫等橡胶构件的 NVH 数值分析。

该方法为整车NVH分析提供了更准确的橡胶隔振件动态特性输入,适合作为前期设计、材料选择与结构优化的依据。


来源:ABAQUS仿真世界
Abaqus振动非线性汽车材料NVH
著作权归作者所有,欢迎分享,未经许可,不得转载
首次发布时间:2026-01-13
最近编辑:7月前
yunduan082
硕士 | 仿真主任工程... Abaqus仿真世界
获赞 203粉丝 392文章 486课程 0
点赞
收藏
作者推荐

LS-DYNA 进行超弹性材料建模

目的:1)计算MooneyRivlin和Ogden材料参数,并使用dogbone试样拉伸试验的应力-应变数据(100%应变)对两者进行比较。2)给定的材料数据是工程应力-应变,单位为MPa/(mm/mm)。3)应该通过d3hsp文件和模拟来进行比较。理论:超弹性材料:超弹性材料模型常用于有限元分析(FEA)中,以表征材料的大变形行为。它们通常用于模拟未填充/填充弹性体的力学行为。除了弹性体之外,超弹性材料模型还用于近似生物组织、聚合物泡沫等的材料行为。线性弹性材料通过两个材料常数(例如杨氏模量和泊松比)来描述。相比之下,超弹性材料通过应变能密度函数来描述。应变能密度可用于推导非线性本构模型(即应力作为大应变变形量(如变形)的函数),这些模型由研究人员推导而来,其中一些被广泛接受的模型,例如Neo-Hookean模型、Mooney-Rivlin模型、Signorini模型、Ogden模型、Yeoh模型和Arruda-Boyce模型,都是常用的模型。Mooney-Rivlin:不可压缩Mooney-Rivlin材料的应变能密度函数定义为其中Cij是经验确定的材料常数,I1和I2分别是第一和第二不变量。通过改变i和j的指标值,我们可以改变本构方程,并使其适用于我们的组件。类似地,我们将有多个系数C。改变指标实际上就是改变能量密度方程多项式的阶数。我们需要找出不同阶数多项式对应的材料模型系数,并从D3HSp文件中提取最佳拟合曲线数据,并将其与材料曲线进行比较。程序:1.实验数据:1)下图所示的超弹性材料的工程应力和工程应变数据已在Excel中绘制出来。2.部件定义:1)为了验证和确认给定的材料数据与Mooney-Rivlin和Ogden的材料模型,取一个狗骨试样进行分析。2)为了找到Mooney-Rivlin常数和Ogden常数,N的值设定为1-3。密度和泊松比的值取为橡胶材料的通用值。对于N>0,采用单轴试验数据。SGL:试样标距长度SW:试样宽度ST:试样厚度3)将SGL、SW和ST的值设置为1(1.0),然后使用*DEFINE_CURVE工程应力与工程应变定义曲线LCID1。4)将狗骨试件的截面属性指定为厚度为2毫米的壳单元,ELFORM=2。5)将材料和截面分配给狗骨状试样部分。3.边界条件:1)由于拉伸试验过程中会发生横向膨胀,固定端的节点仅在X和Z方向上受到约束,而Y方向不受约束。2)由于中性轴沿X方向穿过试样的中间,因此中间的节点在Y方向上受到约束。3)使用位移载荷曲线LCID,为拉动端的节点分配X方向上的边界规定运动。4.控制卡:1)*CONTROL_IMPLICIT_GENERAL启用隐式分析并定义相关的控制参数。所有隐式分析都需要此关键字。IMFLAG:隐式/显式分析类型标志:1-隐式分析DT0:隐式分析的初始时间步长:0.0100000IMFORM:元素公式无缝回弹标志:2保留原始元素公式(默认)2)*CONTROL_IMPLICIT_SOLVER线性方程求解器执行CPU密集型的刚度矩阵求逆运算。LSOLVR:线性方程求解器方法:6(双精度求解器)3)*CONTROL_IMPLICIT_AUTO卡已定义,用于调整时间步长IAUTO:设置为1,自动调整时间步长。ITEOPT:优化操作计数步骤:11ITEWIN:是迭代窗口:54)*CONTROL_TERMINATION:此卡用于说明模拟的终止时间。ENDTIM:1(终止时间)5.材料模型验证的分析设置:1)在notepad++中打开每个材料模型的d3hsp文件,并针对每种情况获得Mooney-Rivlin和Ogden的材料常数以及拉伸和工程应力的最终拟合数据。2)从材料验证的输出文件中,用Notepad++打开d3hsp文件。得到的Mooney-Rivlin常数为:C10=c1=0.1768E+00C01=c2=0.1474E+00为了验证材料模型,材料卡中的N值被更改为0。这些值被输入到材料卡MAT_77_H中。从d3hsp文件获取的数据包括拉伸量和真实应力,真实应力被转换为工程应变和工程应力,以便与实验数据进行验证。3)从Ogden材料模型的d3hsp输出文件中获得的Ogden常数为:μ1=9.1076025585232E−01α1=1.3427238694706E+004)这些值被输入到材料卡MAT_77_O中。从d3hsp文件获取的数据是工程应力和拉伸比,将其转换为工程应变,以便与实验数据进行验证。5)用于材料模型验证的关键字与材料模型确认的关键字相同,区别在于材料卡(N=0)的更改以及用于计算弹性应变的数据库扩展二进制数据的添加。关键字文件保存后,可在LS-DYNA程序管理器中运行,以获得所需的输出文件。结果与讨论:1.冯·米塞斯应力动画轮廓:Mooney-Rivlin材料模型(N=1)的最大vm应力值为1.66801MPa,而Ogden材料模型(N=1)的最大vm应力值为1.65691MPa。因此,两种材料模型的最大vm应力值相等。2.有效塑性应变的动画轮廓:两种材料模型产生的最大有效塑性应变均为3.8。3.材料模型验证:从d3hsp文件中获取不同多项式值(N=1,2,3)的数据,在Excel中绘制工程应力与工程应变曲线,并与实验数据进行比较,结果如下图所示。从图中可以看出,不同多项式的曲线重合良好,与实验数据吻合度更高。与Mooney-Rivlin材料模型相比,Ogden材料模型拟合效果更佳。4.材料模型验证:为了验证结果,将狗骨形试件中间单元的真实应力-真实应变曲线保存为.csv文件。该文件在Excel中打开,并将真实应力和真实应变值转换为工程应力和工程应变,绘制出如下图所示的曲线。从图中可以看出,两种材料模型的模拟结果与实验数据均存在轻微偏差,这是因为材料具有高度非线性,且模拟中采用了通用密度和泊松比值。5.Mooney-Rivlin材料模型与Ogden材料模型的比较:从图中可以看出,对于给定的材料数据,两种材料模型表现相似。Ogden材料模型的多项式范围为1到8,易于使用,并且能更好地拟合实验数据。结论:1)Mooney-Rivlin和Ogden材料常数可通过在d3hsp文件中输入材料卡MAT_77_H和MAT_77_O中多项式N=1的值获得。2)通过d3hsp获得的两种超弹性本构模型的拉伸应力曲线数据与实验材料数据进行了验证,发现两者吻合良好。3)两种材料模型的拉伸试验模拟结果与实验数据相比略有偏差,这是由于所给超弹性材料具有高度非线性行为以及密度和泊松比的一般值所致。4)从模拟结果来看,两种超弹性本构材料模型均能有效描述给定材料数据的行为。5)Ogden材料模型具有更多1到8的多项式,易于使用,并且能够更好地拟合实验数据,但代价是计算时间较长。6)学会了验证和确认给定材料数据的超弹性材料模型。来源:ABAQUS仿真世界

未登录
还没有评论
课程
培训
服务
行家
VIP会员 学习计划 福利任务
下载APP
联系我们
帮助与反馈