汽车用橡胶超弹性材料参数反求方法研究
21世纪以来我国汽车工业飞速发展,汽车的安全舒适性能已经是消费者关注的重要性能。为了达到隔振减振的效果,汽车结构中采用许多橡胶柔性连接,主要应用在发动机悬置系统、驱动系统、车身系统等;操纵稳定性影响汽车的安全性,为了提高车辆的操纵稳定性,汽车前后悬架、操纵机构中也采用很多橡胶柔性连接。汽车研发技术不断的提高要求汽车子系统也要向高性能方向发展,汽车用橡胶结构件自身设计过程中的计算机仿真分析越来越受到重视,而在工程实际中橡胶超弹性材料参数的缺少给仿真分析带来不便。描述橡胶材料力学性能的本构模型有多种形式,橡胶超弹性材料参数获得的传统方法来自橡胶材料基础试验,需要制备专门的试件,由于某些原因无法得到标准试件;试验标准不完善、国内试验条件的不足和试验设备的缺乏导致这些标准试验很难完成。在这种情况下,本文应用仿真分析、优化算法结合的反求方法的思路,提出了橡胶超弹性材料参数的反求方法。在反求过程中引入响应面优化的方法,对响应面建立过程中两种&
(本文共68页)
权威出处：
近日,马来西亚种植业和商品部表示,今年1据了解,2016年马来西亚橡胶及橡胶制品出~7月橡胶及橡胶制品出口额比2016年同期的口额为181亿令吉。销售额的增长主要来自美138亿令吉增长38.5%至191亿令吉,预计今年国、欧洲、中国和印度的购买支持,部分是由于一将达到200亿令吉(47....&
(本文共1页)
权威出处：
据外媒报道,泰国天然橡胶管理局和泰国国道局在橡胶整合利用方面的广泛合作有了进展,首个合作项目就是从2018年开始,利用橡胶作为原料生产公路警示桩,预计消耗橡胶13万t。作为全球天然橡胶最重要的产区和出口国,泰国橡胶绝大部分用于出口,而国内利用和使用的比例相对较低。因此,泰国政府极力希望国内各行业加大橡胶作为原料或配料的使用...&
(本文共1页)
权威出处：
为发展家庭式橡胶种植业,柬埔寨农业部日前计划向法国开发署(AFD)申请1300万美元贷款。这笔贷款将用于支持拉达那基里省和马德望省农民种植橡胶。目前,柬埔寨农业部正在研究拉达那基里省和马德望省共6000公顷橡胶种植土地。这2个省已有2500户家庭表示有意投...&
(本文共1页)
权威出处：
泰国橡胶协会、橡胶种植农联盟、泰国天然橡家橡胶机构暂缓出清31万t橡胶,并将出清库存胶委员会主席邬泰近期表示,与南部16府橡胶棕橡胶行动推迟至3个月后,即日榈种值农团体负责人拓沙蓬、全国橡胶农联盟代以后再执行,并冻结上述库存...&
(本文共1页)
权威出处：
10.14028/j.cnki.17.01.010引言高分子物理是高分子材料与工程专业的一门专业基础课,其教学目标是以高分子的结构特点及结构与性能的关系为主线,在充分了解高分子结构的基础上,从分子运动的角度,讨论高分子结构与外界条件对高分子材料宏观性能的影响,建立高分子结构与性能之间的关系[1~3]。该课程内容多、概念多、头绪多、关系多、数学推导多[4,5],对学生来说是一门难度较大的专业基础课。橡胶是一种具有优异性能的高分子材料。早期的橡胶是取自橡胶树、橡胶草等植物的胶乳,加工后制成具有弹性和绝缘性、不透水与空气的材料。合成橡胶是化学工业中占有重要地位的三大合成材料之一,它与天然橡胶一起构成了橡胶加工行业的最基本原材料,广泛应用于制造轮胎、胶管、胶带、胶鞋、机械配件等橡胶制品。合成橡胶已经成为我国国民经济发展中不可缺少的材料之一。橡胶材料的结构与性能在很多方面综合了高分子物理的知识。在高分子物理的教学过...&
(本文共4页)
权威出处：
扩展阅读：
CNKI手机学问
有学问，才够权威！
xuewen.cnki.net
出版：《中国学术期刊（光盘版）》电子杂志社有限公司
地址：北京清华大学 84-48信箱 大众知识服务
京ICP证040431号&
服务咨询：400-810--9993
订购咨询：400-819-9993
传真：010-超弹性分析的概括和执行【元计算吧】_百度贴吧
&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&签到排名：今日本吧第个签到，本吧因你更精彩，明天继续来努力！
本吧签到人数：0可签7级以上的吧50个
本月漏签0次！成为超级会员，赠送8张补签卡连续签到：天&&累计签到：天超级会员单次开通12个月以上，赠送连续签到卡3张
关注：29贴子：
超弹性分析的概括和执行
(1)概要　　超弹性的性质包括应力和应变的关系如图 1 的那样不是直线，并且还有大应变，卸载时沿着加载路径的反向返回，载荷回到 0，则应变(变形)也为 0。　　变形返回原来的样子这一点来说是弹性的，而纵弹性模量所依赖的应变这一点上却是 非线性。　　这样性质的材料的分析用超弹性分析，多数场合伴随着大变形或大应变。橡胶材料是超弹性材料中的一种，这种场合泊松比为 0.5(非压缩材料)或者在其附近。 用橡胶材料制成的 O 型环、垫圈、衬套、密封垫、轮胎等，在处理用眼能看得见大的　　变形场合，一起利用大变形，大应变等的性能。 即使金属由于应力的有或无内部结晶结构会发生变化，由它的变态大应变度可达到超　　过 10%，而仍持有超弹性性质。被实用化了的有眼镜架等。(2)执行　　分析的执行顺序与线性分析相同，但必须用材料物理性质输入超弹性材料的特性。 根据需要指定考虑大变形或大应变。由分析得到应力、位移、应变等。附：FELAC 2.0软件简介　　FELAC 2.0采用自定义的有限元语言作为脚本代码语言，它可以使用户以一种类似于数学公式书写和推导的方式，非常自然和简单的表达待解问题的微分方程表达式和算法表达式，并由生成器解释产生完整的并行有限元计算C程序。　　FELAC 2.0的目标是通过输入微分方程表达式和算法之后，就可以得到所有有限元计算的程序代码，包含串行程序和并行程序。该系统采用一种语言(有限元语言)和四种技术(对象技术、组件技术、公式库技术生成器技术)开发而成。并且基于FELAC 1.0的用户界面，新版本扩充了工作目录中右键编译功能、命令终端输入功能，并且丰富了文本编辑功能，改善了用户的视觉体验，方便用户快速便捷的对脚本或程序进行编辑、编译与调试。其中并行版在前后处理上进行了相应的改进。
贴吧热议榜
使用签名档&&
保存至快速回贴基于Mooney_Rivlin和h模型的超弹性橡胶材料有限元分析_黄建龙_百度文库
赠送免券下载特权
10W篇文档免费专享
部分付费文档8折起
每天抽奖多种福利
两大类热门资源免费畅读
续费一年阅读会员，立省24元！
基于Mooney_Rivlin和h模型的超弹性橡胶材料有限元分析_黄建龙
阅读已结束，下载本文需要
想免费下载本文？
定制HR最喜欢的简历
下载文档到电脑，同时保存到云知识，更方便管理
加入VIP
还剩2页未读，
定制HR最喜欢的简历
你可能喜欢论文发表、论文指导
周一至周五
9：00&22：00
基于Creo Simulate的橡胶减震器超弹性材料分析
　　本文介绍了Creo Simulate超弹性材料模型的本构方程和材料系数拟合方法，在Creo Simulate中对橡胶减震器进行了超弹性材料、大变形和接触三个方面的组合非线性分析，模拟了减震器的变形情况和应力状态。 中国论文网 /1/view-5370700.htm　　本文介绍了Creo Simulate超弹性材料模型的本构方程和材料系数拟合方法，在Creo Simulate中对橡胶减震器进行了超弹性材料、大变形和接触三个方面的组合非线性分析，模拟了减震器的变形情况和应力状态。 　　一、超弹性材料的力学性能 　　超弹性材料（例如橡胶）是指可对大应变进行瞬时弹性响应的非线性材料，如图1所示的典型固体橡胶材料单轴拉伸应力―应变曲线。在低应变区的弹性系数为1MPa左右，仅为钢（弹性系数约为2×105M P a）的二十万分之一。 　　通常橡胶可以拉长到原长的600%，最长可达1000%，而钢仅在伸长1%时才保持弹性。橡胶的热学属性表现为受热缩短，和受热膨胀的其他固体相反。拉伸状态下，材料先软化再硬化，而压缩时材料急剧硬化。 　　橡胶是链状高分子聚合物，添加硫或其他无机物使链状分子相互搭桥形成网状结构，具有高度弹性，而且几乎不发生体积变化。自然状态这些链状高分子处于无规则蜷缩状态，受拉时卷曲链状分子通过内部旋转被拉直，但被拉直的链状分子的无规则运动力图使其恢复卷曲状态。在整个变形过程中，橡胶材料的应力应变关系是非线性的，但是在一个很小的增量段范围内仍可看成是线性的，因此可以使用增量形式创建橡胶材料的应力应变关系。 　　二、Creo Simulate中选用的超弹性材料模型 　　Creo Simulate包括下列6种超弹性材料模型（表1）。 　　下面对这些材料模型进行简要说明。 　　橡胶材料的应力应变关系表达有两种理论，一种是统计热力学，有高斯模型和非高斯模型；另一种是基于连续介质力学，把橡胶作为一个连续统一体的唯象理论。前者认为橡胶弹性恢复力主要来自熵的减少，橡胶的伸长使得橡胶结构由高度无序变得有序，由对橡胶中分子链的长度、方向以及结构的同价得到橡胶的本构关系，软件中也提供了Arruda-Boyce模型；后者假设在未变形状态下橡胶为各向同性材料，长分子链方向在橡胶中随机分布，该假设用单位体积弹性应变能密度描述橡胶特性，软件中也提供了多项式形式模型特例。 　　1.统计热力学模型 　　（1）高斯统计模型。 　　（1） 　　式中，n 为平均单位体积的网链数，k 为Boltzmann常数，T 为绝对温度，λ i为主伸长率。 　　（2）非高斯统计模型。 　　有单链、满链、P链等模型形式，P链模型有3链、4链、8链等。典型的8链模型是Arruda-Boyce模型，Arruda-Boyce模型可以在较少的试验数据下得到较好的结果，但是当材料发生大变形时，计算结果不精确，原因是其中的朗之万反函数在展开过程中仅取了前几项（例如5项）。该模型需要的实验数据很少，应变可达300%。 　　2.多项式形式 　　对于各向同性材料，应变能密度函数分解成应变偏量能和体积应变能两部分，形式为： 　　（2） 　　（1）多项式中的项说明如下。 　　W为应变能密度，I 1、I 2、I 3是变形张量不变量，l 1、l 2、l 3是主伸长率。体积比J =l 1l 2l 3=V /V 0，热膨胀体积变形J t h=（1+ε t h）3。J e是弹性体积比，对于不可压缩材料，J e=1。 　　弹性体积变形、总体积变形与热体积变形的关系为： 　　J e=J=J total/J th （3） 　　偏差主伸长率和偏差不变量被定义为： 　　其中p =1，2，3。 　　（2）将完全多项式的I 2项略掉，C i j=0（j ≠0），则可以得到缩减多项式，也称之为简化多项式。 　　（4） 　　（3）对于完全多项式。 　　N =1即为Mooney-Rivlin模型二项形式，该模型在工程上广泛应用于弹性体变形研究。但是该模型处理双轴拉伸和平面拉伸（纯剪切）问题时，效果非常差。为此高阶项的Mooney-Rivlin模型相继被提出。该模型是比较常用的本构模型，对于没有碳黑的橡胶来说，能得到比较准确的效果。对密封圈、轮胎、O型圈等材料进行分析时常用，但是模拟加了碳黑的橡胶不适用。M-R模型的适用范围：无填料，平衡态。无填料是指没有炭黑、白炭黑等增强粒子，平衡态是指橡胶中的分子链要在外力作用下达到平衡，充分松弛，所以施加外力要非常缓慢，对于拉力机上的橡胶拉伸是不合适的，拉伸应变可达90%～100%。 　　N =2即为2阶多项式模型，也就是Mooney-Rivlin模型五项形式。实际应用的应变能函数不仅仅只是取Rivlin级数方程的所有低阶项，而是要通过观察实验中的应力应变数据，然后根据经验选择所需级数项。拉伸应变可达100%～200%。 　　（4）对于减缩多项式。 　　N =1即为Neo-Hookean模型。当C10=0.5nkT 时，它等价于高斯统计模型。neo-Hookean模型一般只适用于近似预测30%～40%的单轴拉伸和80%～90%的纯剪切橡胶力学行为，也就是说，该模型适用于小应变到中等应变。 　　N =2即为2阶缩减多项式模型。 　　N =3即为Yeoh模型，Yeoh模型和有着常剪切模量的neo-Hookean与Mooney-Rivlin模型不一样，它能描述有随变形而变化的剪切模型的填料橡胶，而且由某种简单变形实验数据拟合的参数可以用来预测其他变形的力学行为，描述的变形范围也较宽，一般适合于模拟大变形，但是它不能很好解释等双轴拉伸实验。Yeoh模型是模拟加了碳黑填料的橡胶的大变形行为，并且可以用简单试验的数据模拟其他变形力学行为，但是不能很好地解释双轴试验数据，在小变形（伸长比1.5）时不适用。
　　表2是6种超弹性材料模型的适用说明。 　　对于完全不可压缩问题，D =0或D i=0。 　　当给材料模型系数赋值不合适时，会出现如下提示。 　　◎对于Neo-Hookean、2阶缩减多项式模型，如图2。 　　◎对于Mooney-Rivlin、2阶完全多项式、Yeoh模型，如图3。 　　◎对于Arruda-Boyce模型，如图4。 　　三、Creo Simulate超弹性材料模型的材料系数确定方法 　　1.通过实验数据拟合材料模型 　　（1）实验类型。 　　对于超弹性材料， 可以创建以下测试类型： 单轴（Uniaxial）（默认）、等双轴（Equibiaxial）、平面（Planar）和体积（Volumetric），如图5。 　　（2）通过实验数据拟合材料。 　　在Creo Simulate用于从测试数据确定超弹性材料属性的过程中，假定此材料类似于橡胶，几乎无法压缩。输入测试数据来定义超弹性材料时，应输入公称或工程应变，以及公称或工程应力，必须指定应力单位，先前为模型设置的主单位制决定了所显示的默认单位。Creo Simulate显示RMS错误和选定材料模型对指定测试的有效性。如果材料模型无效，则Creo Simulate会为该模型显示红色感叹号，如图6。 　　Creo Simulate采用在测试中指定的数据来最佳拟合曲线，以计算任何材料模型的系数值。可以更改这些值，清除“使用最佳拟合系数”（Use Best Fit Coefficients）复选框，然后指定该材料模型的系数值。 　　（3）材料模型的材料系数拟合方法。 　　CreoSimulate通过使用最小二乘拟合算法来最小化实验应力值与拟合应力值间的平方误差之和，从而使材料模型曲线拟合至测试数据。用于最小二乘拟合算法的应力值有归一化和非归一化两种，受config.pro选项“sim_hyperelastic_material_fit”控制，缺省值是normalized（归一化），也可以修改为non-normalized（非归一化）。 　　对于归一化应力值，Creo Simulate将按以下方式计算RMS误差： 　　（5） 　　对于非归一化，Creo Simulate将按以下方式计算RMS误差： 　　（6） 　　其中n 为测试数据点的总个数； 为第i 个测试数据应力值；Ti（Cj）为根据系数Cj从拟合的材料模型曲线得出的应力值。如果选择非归一化选项，则会在较高应变值处获得较好的拟合，而归一化选项会在较低应变值处提供较好的拟合。 　　2.直接设置材料系数 　　对于已经成型的橡胶元件，通常不容易通过上述试验来确定其材料常数。经验公式是通过橡胶的IRHD硬度指标来确定材料的弹性模量和切变模量，再由材料常数和弹性模量的关系来确定材料常数。基本公式为（小应变条件）： 　　（7） 　　将得到的材料常数代入Mooney-Rivlin模型进行计算。通过其他工具已经获得材料模型常数或者已经在Creo Simulate中保存的超弹性材料模型可以直接调用。 　　使用该方式指定材料模型的系数值以定义超弹性材料的方法：清除“由测试定义”（Define By Tests）复选框；在“材料定义”（Ma t e r ia l De f init ion）对话框的“材料模型”（Material Model）区域中选择6种材料模型之一，并指定模型的系数值，如图7。 　　3.超弹性材料数据信息的保存 　　（1）使用“由测试定义”（Define By Tests）选项创建新材料并将材料保存到材料库时，Creo Simulate使用以下指导方针：仅存储材料模型的系数；不保存用于选取材料模型的测试数据，并显示警告消息；会清除“由测试定义”（Define By Tests）复选框。 　　类似地，当编辑材料库中的材料时，无法选择“材料定义”（Material Definition）对话框中的“由测试定义”（Define By Tests）复选框，可在“材料定义”（Material Definition）对话框中选择材料模型并指定材料模型系数。 　　（2）在当前模型中创建的实验数据信息可以保存到该模型中。 　　（3）输入的实验数据最多10组。 　　四、Creo Simulate超弹性材料分析过程 　　将超弹性材料分配给模型中的任何部分后，可以运行大变形静态结构分析（L D A）或任意类型的小应变分析。运行小应变分析时，Creo Simulate会使用材料的小应变属性。 　　非线性选项：如果模型使用超弹性材料，则在打开“静态分析定义”对话框时，“非线性/使用载荷历史”复选框处于选择状态（默认情况下），“超弹性”在非线性选项区域中突出显示。分析具有超弹性材料的模型时，Creo Simulate始终计算大变形，无法清除此复选框。如果模型具有接触界面，“接触”也将突出显示。 　　超弹性分析策略：①分析材料的具体变形行为：组件、超弹性、大变形、摩擦接触、过盈配合和微小载荷，使用位移约束定义预留自由度的连接；②变形情况已经假设，分析材料的应力行为：通过隔离零件、超弹性、大变形和可行性研究获得合适外力或者强制位移。 　　五、橡胶减震器分析 　　1.前处理 　　（1）减震器装配模型。 　　箱体零件通过4个螺栓固定在机架上，每个螺栓连接处包括一个缸套和2个减震器组成的器件组，共4组，刚体的重力和其他作用力通过减震器传递到缸套，如图8。 　　（2）模型过盈配合。
　　箱体与缸套之间通过减震器接触方式约束，减震器和箱体之间存在过盈配合，假设减震器和箱体之间的重叠距离如图9所示。 　　（3）减震器模型几何简化。 　　这里主要研究减震器的超弹性材料性能，因此将箱体零件的大部分切除，仅保留与减震器零件相关的体积部分，如图10。同时仅保留一个减震器组合进行研究。 　　（4）模型载荷设置。 　　假设箱体的重力和其他作用力对4组器件（一个缸套和2个减震器）的作用效果相同，这样可以简化器件组分析时的受力情况，将箱体的综合作用力简化为作用在箱体侧面上的力，如图11。 　　（5）模型约束设置。 　　箱体侧面轴向自由，缸套顶面、内侧柱面、底面三个平面全部约束，底层的减震器底面约束轴向，如图12。因为是实体模型，旋转自由度不能设置。 　　（6）接触定义。 　　该装配体中有10处界面，需要合理设置连接、接触界面方式。如果一个零件两侧各有一个接触或有一个环形接触，则可能约束不足，发生移动。该模型的一个方案是：缸套和减震器之间1个接触，减震器和箱体之间4个接触，减震器之间1个接触，共6个接触，如图13。没有勾选“分割曲面”和“无穷大摩擦”项。 　　（7）材料设置。 　　箱体、缸套使用steel材料，减震器使用Mooney-Rivlin超弹性材料模型。对于大变形分析，Creo Simulate支持线性弹性、弹塑性和超弹性材料。将线性弹性材料用于大变形静态分析时，Creo Simulate会根据线性弹性的自然推广来解释材料属性。使用与线弹性中所用的相同公式，可将E 和ν 转换为拉梅常数λ 和μ 。应力使用Neo-Hookean材料法则进行计算，它与λ 和μ 线性相关。 　　（8）网格划分。 　　大变形分析中不能使用平面、销钉和球连接，使用位移约束定义自由度。大变形分析包括：3D、2D平面应力、2D轴对称或2D平面应变模型，只有实体和质量元素。接触分析不支持壳，网格信息如图14。 　　2.定义和运行分析 　　创建静态分析，勾选“非线性”选项，由于有超弹性材料，“大变形”选项自动勾选，选择接触和超弹性材料，设置“局部网格细化”、“检查接触力”。设置输出步数为6。在定义分析前一定要设置好工作目录，保证结果文件的存放位置。设置使用内存和手动方式，如图15，也可以使用Config.pro选项“sim_solver_memory_allocation”，其值一般设为物理内存的1/4。 　　Creo Simulate以若干个迭代步骤计算结果。当分析达到收敛或执行了最大迭代次数时，数值迭代将停止。默认情况下，最小迭代次数为1，最大迭代次数为200。使用配置选项“sim_max_contact_iterations”可指定接触分析的最大迭代次数。 　　3.查看并评估分析结果 　　（1）Von Misess应力分布。 　　查看整个组件的应力结果（图16），由于橡胶材料和金属材料应力相差太大，橡胶材料结果显示不明显。采用下列方式单独显示减震器的应力情况：显示位置设置“元件/层类型”，仅显示减震器。确保减震器全部在窗口内，“信息→视图最大值/视图最小值”，查看窗口中显示的减震器元件应力的上下限度值（从信息栏读取），设置标签。 　　（2）位移分布。 　　过盈配合部位的位移最大，可以看出减震器既有上下方向的移动，也有环向的移动，两个减震器之间的间隙出现中间倾斜角度，两个减震器的压缩状态分布也不一样，这主要是受减震器与缸套零件的不同部位接触的影响（图17）。 　　（3）接触压力。 　　箱体与减震器的交线位置接触压力最大，这也反应了减震器最容易破坏的位置。由于简化原则舍弃了缸套内侧端面和减震器外端面之间的接触定义，图18中没有显示该部分接触面。 　　（4）单位体积应变能。 　　箱体、缸套金属材料的应变能密度非常小，应变能密度最大的位置在减震器与缸套、箱体的内侧交界线处，如图19。 　　六、注意事项 　　和通用的有限元分析软件一样，由于受材料模型数据获取方法的限制，使用Creo Simulate超弹性材料分析同样需要注意下列几点。 　　（1）通过实验拟合获取材料模型系数的影响。 　　好的实验数据拟合是决定采用哪种超弹性模型的最好方式。确保实验数据包括预期的应变范围，若数据针对50%应变，不要指望它能很好地与200%应变相关。 　　确保实验数据包括预期的变形模式。对于复杂的材料响应，仅单轴拉伸的数据拟合不能满足，使实验数据与关心的应变范围和变形模式相关联是保证选择了合适的超弹性模型的最好方法。 　　（2）合理选择材料模型。 　　Neo-Hookean和2项Mooney-Rivlin模型是简单的模型，可以用于分析的出发点。注意这些模型不可能很好地预测大的拉伸应变（材料硬化效应）或压缩模式，仅基于单轴数据，Arruda-Boyce模型可很好地预测变形的多重模式。
转载请注明来源。原文地址：
【xzbu】郑重声明：本网站资源、信息来源于网络，完全免费共享，仅供学习和研究使用，版权和著作权归原作者所有，如有不愿意被转载的情况，请通知我们删除已转载的信息。
xzbu发布此信息目的在于传播更多信息，与本网站立场无关。xzbu不保证该信息（包括但不限于文字、数据及图表）准确性、真实性、完整性等。