计算橡胶隔振器静态特性的数值分析方法.docVIP

振 动与冲击 Vo. l 24N o . 32005 第24卷第3期J O URNAL OF VIB RAT I ON AND S HOCK 科研简报 计算橡胶隔振器静态特性的数值分析方法 陈 莲 周海亭 (上海交通大学振动、冲击、噪声国家重点实验室, 上海200030 摘 要 以橡胶隔振器为研究对象, 利用有限元ABAQ US 分析软件, 进行了橡胶隔振器有限元建模的数值试验分 析; 为了能更精确的描述橡胶的特性, 引用了多年来发展的几种弹性应变能函数方法, 选取本构方程来拟合试验数据以作为ABA QUS 软件对橡胶材料特性曲线的评估, 寻找一种描述橡胶特性的最佳的方法。预测的计算结果与试验吻合较好, 表明所建的数学模型和有限元计算方法能较理想的获得橡胶隔振器静态解。 关键词:橡胶隔振器, 应变能, 本构方程, 静态特性中图分类号:TB535. 1 文献标识码:A 0 引 言 常规的橡胶隔振器设计中, 一定形状的橡胶隔振器其应力与应变的关系不是弹性模量E , 而是由形状系数确定的。显然, 形状系数是相当复杂的, 不同的研究者提出不同的计算关系式。其次, 由于橡胶隔振器的性能随着橡胶的配方、硬度、温度、受力状态等诸多因素的不同而变化, 因此以往的橡胶隔振器的设计计算只能是粗略的。 随着橡胶材料机械特性研究的发展及商业有限元软件(如AB A Q US 、M ARC 、A NS YS 、ADI N A 等 逐步具备有对超弹性特性材料模型和体积近似不可压缩的大变形计算功能, 使橡胶制品的有限元仿真技术在工程设计中迅速发展起来。很多研究人员一直致力于有限元法对汽车、舰船、列车等橡胶产品的设计分析 [2, 3] [1] W -59W 1929=2(K (1 B -K B +K B 29K B 9I 19I 式中T B 为应力; W 为应变能函数; K B 为主拉伸T B = 比; I i (i =1, 2 为主应变不变量。确定橡胶材料应变能函数理论有两种。一种是基于统计热力学理论; 另一种是看作连续介质的唯象理论; 目前发展的研究方法都属于这两种理论, 用主应变不变量I i 或主拉伸比K i 表示的应变能函数, 常用的有Ogden 、A rruda Boy 2ce 、Van derWaals 、Yeoh 和Polyno m ial 方程等表示的应变能函数 [4, 5, 6] 。 1. 1 Ogden 方程 Ogden 放弃了应变能函数必须是拉伸比偶次幂的限制, 用级数形式描绘了不可压缩橡胶的应变能函数。W = 。 E n L n A n A n a n (K 1+K 2+K 3-3 A n (2 本文叙述了橡胶超弹性特性本构关系的确定, 应力-应变曲线的评估; 橡胶材料特性试验的测定及橡胶隔振器的有限元建模以及橡胶隔振器静态特性的有限元数值仿真计算方法, 其目的在于寻找描述橡胶特性的最普遍或最简单的途径和表征橡胶在工程应用或其它类似问题中最合适一般方法。 式中A n 可以取任何值, 正或负, 而且不一定是整数, 而L n 为常数。值得注意的是, 在Mooney 方程中包含两项, 分别对应于此公式在A , A 1=22=-2的特例。 与应变能函数(2 相对应的主应力为 t i = E n n L , 2, 3 n K 1-P (i =1 A (3 任意常数P 引进的不确定性, 是同不可压缩性假设有关的, 它在主应力差的公式中不出现, 此方程(3 改写为 t 1-t 2= 1 几种应变能函数 橡胶材料本构关系对于橡胶隔振器的有限元仿真至关重要。橡胶材料是一种各向同性的体积不可压缩或体积近似不可压缩材料, 它的非线性弹性特性 可以用超弹性模型描述, 其材料的本构关系可用弹性应变能函数描述为 收稿日期:200404 [4] E n n n L n (K 1-K 2 A A (4 如果考虑热膨胀效应及橡胶的压缩特性, 其更精 确的应变能函数为: N W = E i=1 2L i i A i A i (K 1+K 2+K 3-3 +A i N E i=1 12i (Jel -1 D i (5 。 09 第一作者 陈 莲 女, 硕士生, 1978年10月生 初始剪切模量及体积模量为: k 0=(6 D 1i=1 式中N 是材料参数, L i 、A i 和D i 是同温度有关的材料参数。 为了利用Ogden 公式, 首先要检验它能否用较少 L 0= N 式中C ij 和D i 12同上。初始剪切模量L 0=2(C01+C 10, k 0同式(13 。 E L i 2 算例 橡胶隔振器数值计算 本研究采用的橡胶隔振器以丁腈橡胶为基材。橡胶隔振器数值计算方法是采用AB AQUS 软件进行建模及有限单元网格划分, 并定义橡胶材料参数和约束材料的物理参数。选取