CN112199823A - 一种橡胶材料应力松弛及损伤效应的数值预测方法 - Google Patents
一种橡胶材料应力松弛及损伤效应的数值预测方法 Download PDFInfo
- Publication number
- CN112199823A CN112199823A CN202010978728.7A CN202010978728A CN112199823A CN 112199823 A CN112199823 A CN 112199823A CN 202010978728 A CN202010978728 A CN 202010978728A CN 112199823 A CN112199823 A CN 112199823A
- Authority
- CN
- China
- Prior art keywords
- relaxation
- damage
- time
- rubber
- stress
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 239000000463 material Substances 0.000 title claims abstract description 131
- 229920001971 elastomer Polymers 0.000 title claims abstract description 125
- 239000005060 rubber Substances 0.000 title claims abstract description 125
- 230000006378 damage Effects 0.000 title claims abstract description 45
- 238000000034 method Methods 0.000 title claims abstract description 35
- 230000000694 effects Effects 0.000 title claims abstract description 23
- 238000012360 testing method Methods 0.000 claims abstract description 76
- 238000004088 simulation Methods 0.000 claims abstract description 42
- 238000004364 calculation method Methods 0.000 claims abstract description 17
- 238000005516 engineering process Methods 0.000 claims abstract description 5
- 239000004033 plastic Substances 0.000 claims description 32
- 238000006073 displacement reaction Methods 0.000 claims description 24
- 230000006835 compression Effects 0.000 claims description 15
- 238000007906 compression Methods 0.000 claims description 15
- 230000006870 function Effects 0.000 claims description 15
- 238000005457 optimization Methods 0.000 claims description 12
- 230000002068 genetic effect Effects 0.000 claims description 10
- 238000004458 analytical method Methods 0.000 claims description 9
- 230000015556 catabolic process Effects 0.000 claims description 9
- 238000006731 degradation reaction Methods 0.000 claims description 9
- 238000004422 calculation algorithm Methods 0.000 claims description 7
- 238000013461 design Methods 0.000 claims description 6
- 238000011161 development Methods 0.000 claims description 4
- 238000005728 strengthening Methods 0.000 claims description 3
- 230000003902 lesion Effects 0.000 claims description 2
- 229920000459 Nitrile rubber Polymers 0.000 description 37
- 238000009864 tensile test Methods 0.000 description 8
- 230000008569 process Effects 0.000 description 5
- 238000012669 compression test Methods 0.000 description 3
- 230000003068 static effect Effects 0.000 description 3
- KVNRLNFWIYMESJ-UHFFFAOYSA-N butyronitrile Chemical compound CCCC#N KVNRLNFWIYMESJ-UHFFFAOYSA-N 0.000 description 2
- 230000005489 elastic deformation Effects 0.000 description 2
- 150000002825 nitriles Chemical class 0.000 description 2
- 230000004044 response Effects 0.000 description 2
- 230000009471 action Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 150000008363 butyronitriles Chemical class 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 230000001143 conditioned effect Effects 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 230000007797 corrosion Effects 0.000 description 1
- 238000005260 corrosion Methods 0.000 description 1
- 230000008878 coupling Effects 0.000 description 1
- 238000010168 coupling process Methods 0.000 description 1
- 238000005859 coupling reaction Methods 0.000 description 1
- 238000005520 cutting process Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000009826 distribution Methods 0.000 description 1
- 230000005284 excitation Effects 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 238000009396 hybridization Methods 0.000 description 1
- 230000007774 longterm Effects 0.000 description 1
- 230000005012 migration Effects 0.000 description 1
- 238000013508 migration Methods 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000035772 mutation Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000001550 time effect Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/10—Numerical modelling
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02P—CLIMATE CHANGE MITIGATION TECHNOLOGIES IN THE PRODUCTION OR PROCESSING OF GOODS
- Y02P90/00—Enabling technologies with a potential contribution to greenhouse gas [GHG] emissions mitigation
- Y02P90/30—Computing systems specially adapted for manufacturing
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
Abstract
本发明公开了一种橡胶材料应力松弛及损伤效应的数值预测方法，包括如下步骤：橡胶材料时变松弛‑损伤本构模型的建立；橡胶材料的力学试验测试；将时变松弛‑损伤本构模型应用于有限元软件，构建时变松弛‑损伤本构关系与有限元数值模拟技术的相容性仿真计算方法；根据第二步获得的橡胶材料的力学试验数据，利用第三步建立的相容性仿真计算方法对橡胶材料时变松弛‑损伤本构模型的关键参数进行拟合识别；根据建立的时变松弛‑损伤本构模型及识别的参数，通过第三步建立的相容性仿真计算方法，对橡胶材料的时变松弛行为及损伤变形数值进行模拟预测。本发明所公开的方法可对不同加载水平的橡胶材料应力松弛及损伤残余特性进行全面分析，预测精确。
Description
技术领域
本发明涉及材料分析技术领域，特别涉及一种橡胶材料应力松弛及损伤效应的数值预测方法。
背景技术
橡胶材料因其质量轻、耐腐蚀性强、非线性刚度及多向刚度等优点，广泛应用于汽车、船舶、电子、建筑及机械等工程领域，具有重要的实际应用价值。在实际作业中，橡胶材料呈现不可避免的时变应力松弛特性，影响其工作性能及使用寿命。橡胶材料的应力松弛是指，当外界加载位移保持某一恒定数值时，材料内部应力随时间缓慢衰减的现象。应力松弛严重影响橡胶材料及其制品的力学性能稳定性(如：内部应力衰减、静刚度退化等)。
橡胶材料的松弛变形可分为弹性变形与不可恢复残余变形。其中，橡胶材料的弹性变形在外界加载位移去除后即可恢复；但不可恢复的残余变形与加载时间、位移幅值有关，即使外界松弛载荷移除后，该残余变形仍然存在，整体表现为橡胶材料的宏观损伤效应。松弛卸载后的残余变形将影响橡胶材料的几何尺寸稳定性，进而影响橡胶制品与其它关联结构之间的配合关系。因此，准确预报橡胶材料的应力松弛行为及松弛-损伤变形效应对其工程应用具有十分重要的意义。
建立准确合理的本构模型是进行橡胶应力松弛-损伤特性分析的基础，所建立的本构模型将决定橡胶数值模拟能否全面表征材料的松弛行为及不可恢复损伤变形效应。此外，将所建立本构模型与现有数值模拟技术结合，构建准确的数值分析方法对橡胶制品的松弛性能预测尤为关键。
目前，数值模拟技术(有限元商业软件)中关于橡胶材料的内置本构模型主要集中在超弹性模型和粘弹性模型两类。其中，超弹性本构模型可描述橡胶的非线性大变形特性，但不可表征与时间相关的应力松弛行为；粘弹性本构模型可描述橡胶材料的时变松弛特性，但不可表征位移载荷相关的不可恢复损伤变形。由此可见，仅采用上述材料模型或者上述模型的组合叠加模型均不可对橡胶材料的应力松弛及损伤效应进行全面分析。如果模拟预测的橡胶材料松弛特性与实际情况不符，将会导致设计无效或者服役状态受损，影响橡胶制品乃至整体系统的正常使用，造成工程及经济损失。
发明内容
为解决上述技术问题，本发明提供了一种橡胶材料应力松弛及损伤效应的数值预测方法，以达到可对不同加载水平的橡胶材料应力松弛及损伤残余特性进行全面分析的目的。
为达到上述目的，本发明的技术方案如下：
一种橡胶材料应力松弛及损伤效应的数值预测方法，包括如下步骤：
第一步：橡胶材料时变松弛-损伤本构模型的建立；
第二步：橡胶材料的力学试验测试；
第三步：将时变松弛-损伤本构模型应用于有限元软件，构建时变松弛-损伤本构关系与有限元数值模拟技术的相容性仿真计算方法；
第四步：根据第二步获得的橡胶材料的力学试验数据，利用第三步建立的相容性仿真计算方法对橡胶材料时变松弛-损伤本构模型的关键参数进行拟合识别；
第五步：根据建立的时变松弛-损伤本构模型及识别的参数，通过第三步建立的相容性仿真计算方法，对橡胶材料的时变松弛行为及损伤变形数值进行模拟预测。
上述方案中，所述第一步的具体方法如下：
橡胶材料时变松弛-损伤本构模型在每一个时间步长内的增量dεij表示如下：
其中，
式中，Wr-e为橡胶材料时变松弛本构模型，t为橡胶材料应力松弛的持载时间；k为松弛幅值退化系数，表征松弛的幅值退化程度；r为松弛衰减系数，表征随时间变化的应力逐渐衰减程度；C10，C01，D10为橡胶材料的超弹参数；
为应变不变量；J是弹性体积比；
f是屈服函数，K为强化函数，
为等效塑性应变，
为等效塑性应变的增量，Sij为应力偏张量，
为等效屈服应力；橡胶材料塑性流动过程中的损伤可以分阶段表示，
为不同载荷阶段的等效塑性应变，
为不同载荷阶段的等效应力，当n＝0时，可得初始等效屈服应力
上述方案中，所述第二步的具体方法如下：
(1)开展橡胶试样在不同变形模式下的超弹力学试验，分别进行单轴拉伸力学测试、平面拉伸力学测试、等双轴拉伸力学测试及单轴压缩力学测试；
(2)采用与超弹力学试验测试相同批次的橡胶材料制成应力松弛试样，开展橡胶材料在不同加载位移水平的长时应力松弛试验；
(3)对橡胶试样在应力松弛卸载后的损伤残余变形进行测量，获取不同加载位移水平下的残余变形数据。
上述方案中，所述第三步的具体方法如下：
(1)在有限元软件ABAQUS中，建立橡胶试样的物理模型；
(2)对ABAQUS进行二次开发，基于Fortran语言、UHYPER子程序及ABAQUS脚本文件，将提出的时变松弛-损伤本构模型嵌入至商业有限元软件中，实现材料模型的编写入库；
(3)按照试验工况，对橡胶试样模型进行材料模型赋予、网格划分、边界条件和载荷条件施加、接触关系及求解计算设置，实现橡胶试样的有限元建模与分析。
上述方案中，所述第四步的具体方法如下：
(1)根据已获得的橡胶材料在单轴拉伸、平面拉伸、等双轴拉伸及单轴压缩变形模式下的超弹试验数据，采用最小二乘法对时变松弛-损伤本构模型中的超弹材料参数C01、C10和D10进行拟合识别；
(2)根据已获得的橡胶材料在不同加载位移水平的应力松弛试验数据，对时变松弛-损伤本构模型的损伤塑性应变材料参数
进行拟合识别；
(3)在已识别超弹材料参数和部分损伤塑性应变材料参数的基础上，根据已获得的橡胶材料在不同加载位移水平的应力松弛试验数据，对时变松弛-损伤本构模型的松弛材料参数及损伤塑性材料参数进行拟合识别；具体包括：利用软件Isight集成第三步中建立的橡胶材料相容性仿真计算方法，以时变松弛-损伤本构模型中的松弛材料参数k和r，和损伤塑性应力参数
作为优化设计变量，以松弛试验结果与仿真结果的差值作为优化目标函数，采用多岛遗传优化算法开展材料参数拟合识别。
通过上述技术方案，本发明提供的一种橡胶材料应力松弛及损伤效应的数值预测方法具有如下有益效果：
(1)本发明所提出的时变松弛-损伤本构模型不仅可描述时间尺度下橡胶材料的应力松弛行为，还可准确表征松弛卸载后的损伤残余变形效应，弥补了目前关于橡胶材料应力松弛行为及损伤效应研究的不足。
(2)本发明所提供的橡胶材料应力松弛数值模拟方法具有建模与分析过程明确、预测精度高等优点，可对不同加载水平的橡胶材料应力松弛及损伤残余特性进行全面分析.
(3)本发明所提出的时变松弛-损伤本构模型的参数识别方法，具有易理解、易操作、拟合精度高的优点。
(4)本发明为橡胶材料及制品的实际工程力学性能分析提供了良好的预测模拟方法。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍。
图1为本发明橡胶材料应力松弛及损伤效应的数值预测方法流程图；
图2a-2d为丁腈橡胶不同变形模式下的试验应力应变曲线图，其中，图2a为单轴拉伸试验的应力应变曲线图，图2b为平面拉伸试验的应力应变曲线图，图2c为等双轴拉伸试验应力应变曲线图，图2d为单轴压缩试验的应力应变曲线图；
图3a-3d为氢化丁腈橡胶不同变形模式下的试验应力应变曲线图，其中，图3a为单轴拉伸试验的应力应变曲线图，图3b为平面拉伸试验的应力应变曲线图，图3c为等双轴拉伸试验应力应变曲线图，图3d为单轴压缩试验的应力应变曲线图；
图4a-4b为丁腈橡胶及氢化丁腈橡胶的应力松弛曲线图，其中，图4a为丁腈橡胶的时间-应力松弛曲线图，图4b为氢化丁腈橡胶的时间-应力松弛曲线图；
图5为有限元模型图；
图6a-6h为不同拉伸试验和单轴压缩试验数据与数值拟合的结果对比图，其中，图6a，图6b，图6c，图6d分别为丁腈橡胶单轴拉伸、平面拉伸、等双轴拉伸、单轴压缩过程中的应力应变曲线对比图；图6e，图6f，图6g，图6h分别为氢化丁腈橡胶单轴拉伸、平面拉伸、等双轴拉伸、单轴压缩过程中的应力-应变曲线对比图；
图7为多岛遗传算法(MIGA)识别本构模型参数的流程图；
图8a-8b为应力松弛过程中应力衰减的数值模拟与试验测试结果对比图，其中，图8a为丁腈橡胶时变应力衰减曲线对比图，图8b为氢化丁腈橡胶时变应力衰减曲线对比图；
图9a-9b为应力松弛结束后损伤残余变形的数值模拟与试验测试结果对比图，其中，图9a为丁腈橡胶的损伤残余变形数值对比图，图9b为氢化丁腈橡胶的损伤残余变形数值对比图。
具体实施方式
下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述。
本发明提供了一种橡胶材料应力松弛及损伤效应的数值预测方法，如图1所示，具体实施例如下：
为验证本发明所构建的本构模型及数值模拟方法的可靠性，本实施例对两种不同类型的橡胶材料(丁腈橡胶和氢化丁腈橡胶)进行松弛特性数值预测。具体实施步骤如下：
第一步：时变松弛-损伤本构模型的提出
在工程应用中，橡胶往往处于大变形状态，其力学特征可采用基于应变能函数的超弹本构模型来模拟，基于应变不变量表示的连续介质唯象模型可写成，
式中，
为材料加载响应的能量密度偏量函数，W(J)应变能密度的体积能函数。
式中，
是偏延伸变量。J是橡胶变形后与变形前的弹性体积比，即：
式中，B为左柯西格林变形张量，F为变形梯度张量。
根据上述公式可推导出不同的超弹本构模型，其中，基于多项式形式的超弹本构模型被广泛应用于工业界，可预测橡胶在外荷载激励下的力学响应。为了方便工程应用，本实施例中选定N＝2的多项式形式超弹本构模型(即目前广泛应用的Mooney-Rivlin模型)，其表达式为，
超弹本构模型可用于分析橡胶材料的静力学特性，但其与时间效应无关，无法描述应力松弛特性。为了评估橡胶材料在常量荷载作用下的松弛退化效应，需对传统超弹本构模型进行修正。
本发明引入带有退化效应的本构模型，
将荷载引起的橡胶材料松弛与时间相连，修正的时变超弹松弛本构模型为：
式中，t为应力松弛的持载时间；k为松弛幅值退化系数，r为松弛衰减系数。
在每一个时间步长内，时变松弛本构模型的应变增量(超弹松弛变形)可表示成，
式中，σij为橡胶材料的应变能应力分量，
为超弹与松弛应变增量。
材料的损伤塑性准则反映了材料产生塑性流动时的应力水平，可以用屈服函数来表示，如下公式所示：
式中，f是屈服函数，K为强化函数，
为等效塑性应变，
为初始等效屈服应力。
当f＜K时，橡胶材料处于时变超弹松弛状态，此时只产生超弹松弛变形，此时，
式中，dεij为总应变增量。
当f＝K时，橡胶材料进入超弹松弛-损伤塑性状态，即初始屈服后，橡胶材料的变形由超弹松弛变形与损伤塑性变形共同组成，具体如下，
式中，
为损伤塑性应变增量。
进入超弹松弛-损伤塑性阶段后，随着加载位移的增加，损伤塑性变形也将进一步增加，而确定流动规则是确定残余变形变化趋势的前提条件。流动规则的增量形式如下，
式中，
为等效塑性应变的增量，Sij为应力偏张量，
为等效应力。橡胶材料塑性流动过程中的损伤可以分阶段表示，
为不同载荷阶段的等效塑性应变，
为不同载荷阶段的等效应力，当n＝0时，可得初始等效屈服应力
综上所述，橡胶材料时变松弛-损伤本构模型在每一个时间步长内的增量可以写成如下形式：
第二步：静态拉伸试验、应力松弛试验
(1)测试丁腈橡胶及氢化丁腈橡胶在单轴拉伸、等双轴拉伸、平面拉伸及单轴压缩变形下的应力-应变关系，测试方案为：首先将橡胶材料裁成标准哑铃型试样，进行准静态单轴拉伸力学测试；其次，对平面状橡胶材料进行平面拉伸力学测试；再次，对橡胶材料进行等双轴拉伸力学测试；最后，对橡胶圆柱试样进行单轴压缩力学测试。
针对不同变形模式下的超弹力学试验，分别将橡胶试样安装在相应的试验机上，调整好试验机参数，进行测试。丁腈橡胶在不同变形模式下的应力-应变曲线如图2a-图2d所示，氢化丁腈橡胶在不同变形模式下的应力-应变曲线如图3a-图3d所示。
(2)丁腈橡胶及氢化丁腈橡胶的应力松弛试验
对橡胶开展长时应力松弛试验，测试步骤如下：按照国家标准，首先对试样进行热调节和机械调节，然后进行不同加载水平的应力松弛试验。试验前，记录试样的初始直径和厚度，试验过程中保持压缩位移量不变。参照国家标准，应力松弛的加载压缩应变分别取为试样初始高度值的15％、20％及25％，试验周期为48小时。记录试验过程的时间、力、位移数据，并记录试样应力松弛试验结束时的直径和厚度。丁腈橡胶及氢化丁腈橡胶在不同加载水平下的时间-应力松弛曲线分别如图4a、图4b所示。
第三步：时变松弛-损伤本构关系与有限元数值模拟的相容性计算方法构建
通过二次开发及编译，将本发明所提供的时变松弛-损伤本构模型应用于有限元商业软件，构建时变松弛本构关系与数值模拟技术的相容性仿真计算方法。
(1)在商业有限元软件ABAQUS中建立橡胶试样的物理模型。参考试验样件为圆柱形试样(直径d为10.0mm，高度h为10.0mm)，建立橡胶试样的轴对称几何模型。
(2)对ABAQUS进行二次开发，基于Fortran语言及UHYPER子程序，将提出的时变松弛-损伤本构模型嵌入商业有限元分析软件中，实现材料模型的编写入库。将构建的应变能函数写入UHYPER子程序，定义超弹本构模型中的第一应变不变量、第二应变不变量、弹性体积比、不可压缩性、超弹参数、松弛参数以及各变量之间的数值关系。各变量之间的计算关系如下，
式中，u为UHYPER子程序定义的应变能函数，其中，松弛损伤部分
作适当简化，得到如下关系，
损伤本构模型可通过ABAQUS脚本文件进行输入。
(3)按照试验工况，对橡胶试样模型进行材料模型赋予、网格划分、边界条件和载荷条件施加、接触关系及求解计算设置等，实现橡胶试样的有限元建模与分析。
对于本发明所提出的时变松弛-损伤本构模型，通过ABAQUS脚本文件将塑性参数等材料属性赋值给有限元模型，通过UHYPER自定义程序将超弹和松弛等材料参数赋值给有限元模型。根据模型的对称性，采用四节点轴对称时变耦合四边形、双线性位移杂交单元来划分有限元模型的网格。橡胶试样模型的上、下表面分别与上、下解析刚体接触，位移载荷施加于解析刚体。构建的橡胶试样有限元模型如图5所示。
第四步：橡胶本构模型的关键参数识别
(1)辨识时变应力松弛-损伤本构模型的超弹材料参数
将单轴拉伸、平面拉伸、等双轴拉伸的试验数据导入ABAQUS，采用最小二乘法对时变松弛-损伤本构模型中的超弹材料参数(C01、C10和D10)进行数值拟合。对于每组应力应变试验数据，均使得相对误差Ee最小，
式中，Ee为相对误差，Ti test为试验应力结果，Ti th为按照Mooney-Rivilin超弹本构关系在与试验数据相同应变下对应的应力仿真结果。
将单轴压缩变形的超弹试验数据与上述数值拟合的结果进行对比，确定橡胶材料的超弹参数。丁腈材料参数：C10＝1414930Pa，C01＝-169289Pa，D10＝1.0×10-8；氢化丁腈材料参数：C10＝2030774Pa，C01＝-870827Pa，D10＝1.0×10-8。关于不同变形模式的超弹力学特性，丁腈试验与数值拟合的对比结果如图6a-6d所示，氢化丁腈试验与数值拟合的对比结果如图6e-6h所示。
(2)辨识时变应力松弛-损伤本构模型的损伤材料参数
第二步的试验中已记录橡胶试样松弛前后的尺寸数据，据此可确定时变松弛-损伤本构模型的损伤塑性应变材料参数
丁腈橡胶的试验数据见表1，氢化丁腈橡胶的试验数据见表2。
表1丁腈材料松弛前后的厚度变形结果
表2氢化丁腈橡胶松弛前后的厚度变形结果
可将平均残余变形与试样原始厚度平均值的比值视为等效损伤塑性应变，丁腈橡胶和氢化丁腈橡胶不同压缩位移条件下的等效塑性应变见表3。
表3丁腈橡胶和氢化丁腈橡胶等效塑性应变
(3)辨识时变松弛-损伤本构模型的损伤塑性材料参数(k，r和
)
通过在Isight优化软件中集成ABAUQS搭建橡胶样件本构参数识别的优化模型。以时变松弛-损伤本构模型参数
为优化设计变量，以试验数据和仿真数据之间的数值差异最小作为优化目标，优化模型的数学形式如下所示：
式中，ω为参数优化过程中的目标函数，α1,α2，α3,α4均为比例因子，
为第i组压缩位移下仿真模拟得到的残余变形，
为第i组压缩位移下试验测试得到的残余变形，j为压缩位移的组数，
为试验和仿真得到的残余变形之间差异的平均值，
为第i组压缩位移下仿真模拟得到的松弛曲线与横坐标轴所围面积，
为第i组压缩位移下试验测试得到的松弛曲线与横坐标轴所围面积，
为试验和仿真所获得的时间-应力松弛曲线面积差的平均值，pn,kd,rd为各优化设计变量对应的参数下限，pn+1,ku,ru为各优化设计变量对应的参数上限。
应用多岛遗传算法(MIGA)，以橡胶材料在不同压缩位移下的应力松弛试验数据为依据，识别本构模型的关键参数。多岛遗传算法是一种改进的并行分布遗传算法，与传统的遗传算法相比，它具有更好的全局求解能力、计算效率和收敛能力。结合优化模型，MIGA用于识别时变应力松弛-损伤本构模型参数的流程图如图7所示。具体识别流程如下：
首先，在寻优范围内随机生成整体种群，初始化之后将该种群分成若干个子种群，这些子种群在地理上分布在不同的岛内，子种群中个体携带的遗传信息中包含松弛参数k、r和不同阶段的等效应力
每一次迭代过程均可融合实时的本构模型并进行有限元分析，提取并计算分析结果中的目标变量值；其次，计算每个子种群内个体的适应度值，按照迁移规则选取可迁移个体，实现与各个种群的联姻，可保证遗传信息的多样性；再次，各子种群内个体则按照一定的概率进行突变和交叉，如此循环往复，不停繁衍，不断更新优秀子种群，进而将优良的遗传信息写入新的橡胶材料本构模型中，直至仿真结果满足精度条件。
丁腈橡胶和氢化丁腈橡胶时变松弛-损伤本构模型的参数识别结果见表4。
表4丁腈橡胶和氢化丁腈橡胶损伤塑性材料参数(k，r和
)
第五步：橡胶材料的时变松弛行为及损伤变形数值模拟预测
根据识别得到的时变松弛-损伤本构模型及识别的材料参数，通过第三步建立的有限元模型，对橡胶材料的时变松弛行为进行仿真模拟，应力随时间衰减的数值计算与试验测试结果对比如图8a和图8b所示(其中，图8a为丁腈橡胶结果对比，图8b为氢化丁腈橡胶结果对比)，应力松弛的损伤残余变形数值模拟与试验测试结果对比如图9a和图9b所示(其中，图9a为丁腈橡胶结果对比，图9b为氢化丁腈橡胶结果对比)。
对比结果显示：识别出的本构模型参数能够准确描述橡胶的应力松弛退化行为及损伤变形效应，仿真结果与试验结果误差均小于5％，这说明数值模拟的预测结果准确性较高。
对所公开的实施例的上述说明，使本领域专业技术人员能够实现或使用本发明。对这些实施例的多种修改对本领域的专业技术人员来说将是显而易见的，本文中所定义的一般原理可以在不脱离本发明的精神或范围的情况下，在其它实施例中实现。因此，本发明将不会被限制于本文所示的这些实施例，而是要符合与本文所公开的原理和新颖特点相一致的最宽的范围。
Claims (5)
1.一种橡胶材料应力松弛及损伤效应的数值预测方法，其特征在于，包括如下步骤：
第一步：橡胶材料时变松弛-损伤本构模型的建立；
第二步：橡胶材料的力学试验测试；
第三步：将时变松弛-损伤本构模型应用于有限元软件，构建时变松弛-损伤本构关系与有限元数值模拟技术的相容性仿真计算方法；
第四步：根据第二步获得的橡胶材料的力学试验数据，利用第三步建立的相容性仿真计算方法对橡胶材料时变松弛-损伤本构模型的关键参数进行拟合识别；
第五步：根据建立的时变松弛-损伤本构模型及识别的参数，通过第三步建立的相容性仿真计算方法，对橡胶材料的时变松弛行为及损伤变形数值进行模拟预测。
2.根据权利要求1所述的一种橡胶材料应力松弛及损伤效应的数值预测方法，其特征在于，所述第一步的具体方法如下：
橡胶材料时变松弛-损伤本构模型在每一个时间步长内的增量dεij表示如下：
其中，
式中，Wr-e为橡胶材料时变松弛本构模型，t为橡胶材料应力松弛的持载时间；k为松弛幅值退化系数，表征松弛的幅值退化程度；r为松弛衰减系数，表征随时间变化的应力逐渐衰减程度；C10，C01，D10为橡胶材料的超弹参数；
为应变不变量；J是弹性体积比；
f是屈服函数，K为强化函数，
为等效塑性应变，
为等效塑性应变的增量，Sij为应力偏张量，
为等效屈服应力；橡胶材料塑性流动过程中的损伤可以分阶段表示，
为不同载荷阶段的等效塑性应变，
为不同载荷阶段的等效应力，当n＝0时，可得初始等效屈服应力
3.根据权利要求1所述的一种橡胶材料应力松弛及损伤效应的数值预测方法，其特征在于，所述第二步的具体方法如下：
(1)开展橡胶试样在不同变形模式下的超弹力学试验，分别进行单轴拉伸力学测试、平面拉伸力学测试、等双轴拉伸力学测试及单轴压缩力学测试；
(2)采用与超弹力学试验测试相同批次的橡胶材料制成应力松弛试样，开展橡胶材料在不同加载位移水平的长时应力松弛试验；
(3)对橡胶试样在应力松弛卸载后的损伤残余变形进行测量，获取不同加载位移水平下的残余变形数据。
4.根据权利要求1所述的一种橡胶材料应力松弛及损伤效应的数值预测方法，其特征在于，所述第三步的具体方法如下：
(1)在有限元软件ABAQUS中，建立橡胶试样的物理模型；
(2)对ABAQUS进行二次开发，基于Fortran语言、UHYPER子程序及ABAQUS脚本文件，将提出的时变松弛-损伤本构模型嵌入至商业有限元软件中，实现材料模型的编写入库；
(3)按照试验工况，对橡胶试样模型进行材料模型赋予、网格划分、边界条件和载荷条件施加、接触关系及求解计算设置，实现橡胶试样的有限元建模与分析。
5.根据权利要求3所述的一种橡胶材料应力松弛及损伤效应的数值预测方法，其特征在于，所述第四步的具体方法如下：
(1)根据已获得的橡胶材料在单轴拉伸、平面拉伸、等双轴拉伸及单轴压缩变形模式下的超弹试验数据，采用最小二乘法对时变松弛-损伤本构模型中的超弹材料参数C01、C10和D10进行拟合识别；
(2)根据已获得的橡胶材料在不同加载位移水平的应力松弛试验数据，对时变松弛-损伤本构模型的损伤塑性应变材料参数
进行拟合识别；
(3)在已识别超弹材料参数和部分损伤塑性应变材料参数的基础上，根据已获得的橡胶材料在不同加载位移水平的应力松弛试验数据，对时变松弛-损伤本构模型的松弛材料参数及损伤塑性材料参数进行拟合识别；具体包括：利用软件Isight集成第三步中建立的橡胶材料相容性仿真计算方法，以时变松弛-损伤本构模型中的松弛材料参数k和r，和损伤塑性应力参数
作为优化设计变量，以松弛试验结果与仿真结果的差值作为优化目标函数，采用多岛遗传优化算法开展材料参数拟合识别。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010978728.7A CN112199823B (zh) | 2020-09-17 | 2020-09-17 | 一种橡胶材料应力松弛及损伤效应的数值预测方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010978728.7A CN112199823B (zh) | 2020-09-17 | 2020-09-17 | 一种橡胶材料应力松弛及损伤效应的数值预测方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN112199823A true CN112199823A (zh) | 2021-01-08 |
| CN112199823B CN112199823B (zh) | 2023-04-07 |
Family
ID=74015335
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202010978728.7A Active CN112199823B (zh) | 2020-09-17 | 2020-09-17 | 一种橡胶材料应力松弛及损伤效应的数值预测方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN112199823B (zh) |
Cited By (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112733415A (zh) * | 2021-01-14 | 2021-04-30 | 中国海洋大学 | 一种薄壁弹性体边界的无网格处理方法、装置、终端设备及计算介质 |
| CN112765734A (zh) * | 2021-01-27 | 2021-05-07 | 同济大学 | 一种胶粘异质车身门盖件固化变形及残余内应力预测方法 |
| CN113255134A (zh) * | 2021-05-27 | 2021-08-13 | 中国工程物理研究院总体工程研究所 | 基于匹配追踪的硅泡沫材料物理化学松弛自适应分离方法 |
| CN113627054A (zh) * | 2021-08-02 | 2021-11-09 | 东南大学 | 一种优化的延性金属损伤参数确定方法 |
| CN113808675A (zh) * | 2021-08-13 | 2021-12-17 | 中国人民解放军海军工程大学 | 发泡橡胶材料老化性能反映方法和装置 |
| CN114371076A (zh) * | 2022-01-06 | 2022-04-19 | 上海电气集团股份有限公司 | 工件应力值的测试方法、系统、电子设备及存储介质 |
| CN114566235A (zh) * | 2022-02-24 | 2022-05-31 | 南京理工大学 | 一种测定水凝胶软材料内部应力-应变场的方法 |
| WO2023201522A1 (zh) * | 2022-04-19 | 2023-10-26 | 石家庄铁道大学 | 海水老化的沉管隧道gina止水带时变本构模型构建方法 |
| CN118155770A (zh) * | 2024-03-11 | 2024-06-07 | 中南大学 | 双轴复杂应力状态下应力松弛过程中应力演化的预测方法 |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20150278413A1 (en) * | 2014-03-26 | 2015-10-01 | Livermore Software Technology Corporation | Numerical Model For Rubber-like Materials Suitable For Computer Aided Engineering Analysis |
| CN109408969A (zh) * | 2018-10-28 | 2019-03-01 | 北京工业大学 | 一种应用有限元软件辨识橡胶粘弹参数建立本构模型的方法 |
-
2020
- 2020-09-17 CN CN202010978728.7A patent/CN112199823B/zh active Active
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20150278413A1 (en) * | 2014-03-26 | 2015-10-01 | Livermore Software Technology Corporation | Numerical Model For Rubber-like Materials Suitable For Computer Aided Engineering Analysis |
| CN109408969A (zh) * | 2018-10-28 | 2019-03-01 | 北京工业大学 | 一种应用有限元软件辨识橡胶粘弹参数建立本构模型的方法 |
Non-Patent Citations (4)
| Title |
|---|
| DINGXIN LENG等: "A hyper-elastic creep approach and characterization analysis for rubber vibration systems", 《POLYMERS》 * |
| 张江涛等: "橡胶粘结颗粒材料粘弹性性能的试验研究与离散元模拟", 《固体力学学报》 * |
| 袁子豪等: "减振器橡胶动态剪切特性测试及黏-超弹模型仿真应用", 《航空科学技术》 * |
| 赵排航等: "采用多岛遗传算法的某型狙击榴弹发射器参数优化", 《火力与指挥控制》 * |
Cited By (14)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112733415B (zh) * | 2021-01-14 | 2022-07-01 | 中国海洋大学 | 一种薄壁弹性体边界的无网格处理方法、装置、终端设备及计算介质 |
| CN112733415A (zh) * | 2021-01-14 | 2021-04-30 | 中国海洋大学 | 一种薄壁弹性体边界的无网格处理方法、装置、终端设备及计算介质 |
| CN112765734A (zh) * | 2021-01-27 | 2021-05-07 | 同济大学 | 一种胶粘异质车身门盖件固化变形及残余内应力预测方法 |
| CN112765734B (zh) * | 2021-01-27 | 2023-06-02 | 同济大学 | 一种胶粘异质车身门盖件固化变形及残余内应力预测方法 |
| CN113255134A (zh) * | 2021-05-27 | 2021-08-13 | 中国工程物理研究院总体工程研究所 | 基于匹配追踪的硅泡沫材料物理化学松弛自适应分离方法 |
| CN113627054A (zh) * | 2021-08-02 | 2021-11-09 | 东南大学 | 一种优化的延性金属损伤参数确定方法 |
| CN113627054B (zh) * | 2021-08-02 | 2022-11-22 | 东南大学 | 一种优化的延性金属损伤参数确定方法 |
| CN113808675A (zh) * | 2021-08-13 | 2021-12-17 | 中国人民解放军海军工程大学 | 发泡橡胶材料老化性能反映方法和装置 |
| CN113808675B (zh) * | 2021-08-13 | 2023-12-12 | 中国人民解放军海军工程大学 | 发泡橡胶材料老化性能反映方法和装置 |
| CN114371076A (zh) * | 2022-01-06 | 2022-04-19 | 上海电气集团股份有限公司 | 工件应力值的测试方法、系统、电子设备及存储介质 |
| CN114566235A (zh) * | 2022-02-24 | 2022-05-31 | 南京理工大学 | 一种测定水凝胶软材料内部应力-应变场的方法 |
| CN114566235B (zh) * | 2022-02-24 | 2024-08-23 | 南京理工大学 | 一种测定水凝胶软材料内部应力-应变场的方法 |
| WO2023201522A1 (zh) * | 2022-04-19 | 2023-10-26 | 石家庄铁道大学 | 海水老化的沉管隧道gina止水带时变本构模型构建方法 |
| CN118155770A (zh) * | 2024-03-11 | 2024-06-07 | 中南大学 | 双轴复杂应力状态下应力松弛过程中应力演化的预测方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN112199823B (zh) | 2023-04-07 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN112199823B (zh) | 一种橡胶材料应力松弛及损伤效应的数值预测方法 | |
| Alliche | Damage model for fatigue loading of concrete | |
| Li et al. | Simulation of cyclic stress/strain evolutions for multiaxial fatigue life prediction | |
| Yang et al. | Structural damage localization and quantification using static test data | |
| Chalal et al. | Direct identification of the damage behaviour of composite materials using the virtual fields method | |
| CN110501126B (zh) | 基于支座反力和挠度曲率的梁结构损伤识别方法 | |
| CN112182935B (zh) | 基于兰姆波缺陷定量的结构可靠性敏感性分析方法 | |
| Meng et al. | Study on stress relaxation of membrane structures in the prestress state by considering viscoelastic properties of coated fabrics | |
| Liu et al. | Probabilistic fatigue life prediction model of natural rubber components based on the expanded sample data | |
| Lina et al. | FEM analysis of spring-backs in age forming of aluminum alloy plates | |
| Leng et al. | Experimental mechanics and numerical prediction on stress relaxation and unrecoverable damage characteristics of rubber materials | |
| KR101169393B1 (ko) | 고무 물성 평가 방법 | |
| CN108548720A (zh) | I型裂纹弹塑性理论公式获取延性材料j阻力曲线的方法 | |
| Peng et al. | Damage detection in beam-like structures using static shear energy redistribution | |
| KR101655058B1 (ko) | 모드 기여도를 고려한 건축 구조 부재의 휨강성 추정방법 및 이를 이용한 안전성 평가방법 | |
| Franulović et al. | Automation of LCF material model parameters’ identification | |
| CN116992708A (zh) | 一种基于壳单元的钢结构超低周疲劳断裂模拟方法 | |
| Wang | Wei-Qi | |
| Muliana et al. | Nested nonlinear viscoelastic and micromechanical models for the analysis of pultruded composite materials and structures | |
| CN115618676A (zh) | 一种基于连续累积损伤的低周疲劳裂纹的模拟方法及系统 | |
| CN114048646B (zh) | 基于有限元修正与人工智能的沥青面层损伤状态反演方法 | |
| Mohanty et al. | Effect of pressurized water reactor environment on material parameters of 316 stainless steel: a cyclic plasticity based evolutionary material modeling approach | |
| CN108984887A (zh) | 复合材料参数多阶段识别方法 | |
| CN115329644A (zh) | 一种大变形橡胶材料超弹复合行为精准仿真方法 | |
| Friswellt et al. | The location of damage from vibration data using genetic algorithms |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |