类橡胶材料变形直到破坏的显式本构模型

本文通过直接、显式的方法提出一个多轴可压缩应变能弹性势来模拟类橡胶材料受载荷直到软化破坏的变形行为．首先，我们提出一个多轴可压缩应变能函数；其次，通过特定的不变量，该多轴应变能函数在单轴拉伸，平面应变和等双轴拉伸三个基准实验的情况下，可以退化为各自的单轴形函数形式；再次，我们显式给出带有软化破坏特性的形函数；最后，模型结果和试验数据可以精确匹配，同时可以预测材料临近破坏以后，接下来的变形行为．     

显式模拟类橡胶材料Mullins效应滞回圈

通过显式、直接的方法提出一个多轴可压缩应变能函数,用来模拟类橡胶材料在加载——卸载作用下,由于Mullins效应而产生的应力——应变滞回圈. 本文的创新点在于将表征能量耗散的变量引入到应变能函数.新的弹性势具有以下两个特点:第一,在加载情况下,新引入的变量不会对弹性势产生任何影响,因此,只要给出合适的形函数显式表达,3个基准实验,包括单轴拉伸和压缩,等双轴拉伸和压缩,以及平面应变,都可精确模拟;第二,新引入的变量在卸载情况下将被激活.在不同的卸载应力下,变量将发生改变,从而影响弹性势,使其最终产生不同的应力——应变关系卸载曲线,与对应的加载曲线共同构成应力——应变滞回圈.通过对Mullins效应实验数据进行分析和研究,得出了卸载形函数在不同卸载应力下变化的规律,并预测不同卸载应力下的应力——应变关系.最后,我们将得到精确匹配实验数据的数值模拟结果,从而证明本文方法不仅可以精确匹配至少3个基准实验,还可以模拟和预测类橡胶材料在加载——卸载作用下由于Mullins效应而产生的滞回圈.      

显式模拟类橡胶材料Mullins效应滞回圈

通过显式、直接的方法提出一个多轴可压缩应变能函数,用来模拟类橡胶材料在加载-卸载作用下,由于Mullins效应而产生的应力-应变滞回圈.本文的创新点在于将表征能量耗散的变量引入到应变能函数.新的弹性势具有以下两个特点:第一,在加载情况下,新引入的变量不会对弹性势产生任何影响,因此,只要给出合适的形函数显式表达,3个基准实验,包括单轴拉伸和压缩,等双轴拉伸和压缩,以及平面应变,都可精确模拟;第二,新引入的变量在卸载情况下将被激活.在不同的卸载应力下,变量将发生改变,从而影响弹性势,使其最终产生不同的应力-应变关系卸载曲线,与对应的加载曲线共同构成应力-应变滞回圈.通过对Mullins效应实验数据进行分析和研究,得出了卸载形函数在不同卸载应力下变化的规律,并预测不同卸载应力下的应力-应变关系.最后,我们将得到精确匹配实验数据的数值模拟结果,从而证明本文方法不仅可以精确匹配至少3个基准实验,还可以模拟和预测类橡胶材料在加载-卸载作用下由于Mullins效应而产生的滞回圈.     

帘线/橡胶复合材料各向异性黏-超弹性本构模型

帘线/橡胶复合材料广泛应用于轮胎等重要工程领域,为了描述其在服役条件下的大变形、非线性、各向异性和高应变率等材料力学行为,基于纤维增强复合材料连续介质力学理论,提出了一种考虑应变率效应的帘线/橡胶复合材料各向异性黏-超弹性本构模型. 该模型中单位体积的应变能被解耦为便于参数识别的基体等容变形能、帘线拉伸变形能、剪切应变能和黏性应变能四部分. 给出了模型参数的确定方法,并通过拟合文献中单轴拉伸、偏轴拉伸实验数据,得到了模型参数. 利用该模型预测了不同加载和变形条件下的力学行为,并将预测结果与实验结果对比分析. 结果表明, 考虑黏性模型和不考虑黏性模型对不同应变率变形条件下的预测结果相差很大,且考虑黏性模型的预测结果与实验结果吻合很好. 因此,与不考虑黏性模型相比,所提出的各向异性黏-超弹性本构模型能更好地表征帘线/橡胶复合材料在大变形、高应变率条件下的力学特性.      

基于横观各向同性超弹性理论的短纤维增强橡胶本构模型的建立与应用Establishment and application of short fiber reinforced rubber constitutive model based on transversely isotropic hyperelastic theory

提出了一种能够表征短纤维增强橡胶的横观各向同性超弹性本构模型,并结合试验体系,对其在数值分析中的应用方法和效果进行了研究。基于连续介质力学理论,建立了横观各向同性材料的应变能函数,推导得到不同变形形式下的应力应变关系,给出材料参数辨识试验方法,并成功应用于某短纤维增强橡胶测试中,得到表征其超弹性特性的相关材料参数。利用有限元软件ANSYS对不同纤维排布方向的单轴拉伸和平行纤维方向的平面拉伸进行仿真计算,并对比相应试验数据,以验证材料参数的可靠性。最后基于已验证的本构模型,建立了某铣槽装备减振环仿真模型,并对其进行了校核计算。研究结果表明,本文提出的本构模型能够有效表征短纤维增强橡胶的静态力学特性并且方便嵌入现有的有限元软件中,具有材料参数少、测试简便和结果准确等特点,工程实用性强。     

基于横观各向同性超弹性理论的短纤维增强橡胶本构模型的建立与应用

提出了一种能够表征短纤维增强橡胶的横观各向同性超弹性本构模型,并结合试验体系,对其在数值分析中的应用方法和效果进行了研究。基于连续介质力学理论,建立了横观各向同性材料的应变能函数,推导得到不同变形形式下的应力应变关系,给出材料参数辨识试验方法,并成功应用于某短纤维增强橡胶测试中,得到表征其超弹性特性的相关材料参数。利用有限元软件ANSYS对不同纤维排布方向的单轴拉伸和平行纤维方向的平面拉伸进行仿真计算,并对比相应试验数据,以验证材料参数的可靠性。最后基于已验证的本构模型,建立了某铣槽装备减振环仿真模型,并对其进行了校核计算。研究结果表明,本文提出的本构模型能够有效表征短纤维增强橡胶的静态力学特性并且方便嵌入现有的有限元软件中,具有材料参数少、测试简便和结果准确等特点,工程实用性强。     

短纤维增强EPDM包覆薄膜超弹性本构模型

短纤维增强三元乙丙橡胶(EPDM)包覆薄膜用于一种新型缠绕包覆工艺,主要解决复杂构型自由装填药柱外表面可靠性包覆问题.为了描述其在固体火箭发动机工作过程中产生的大变形、非线性和各向异性等力学行为,根据纤维增强复合材料连续介质力学理论,提出了各向异性超弹性本构模型.该模型中单位体积的应变能函数被解耦成两部分:表征各向同性的橡胶基体应变能和表征各向异性的纤维拉伸应变能,通过引入纤维方向对纤维应变能进行修正,给出了通过单轴拉伸、偏轴拉伸实验数据获取模型参数的具体方法.研究结果表明,该模型能够很好地预测材料在纤维方向0°~45°时的各向异性力学特性,并将预测结果与实验数据对比,误差在5%以下.所建立的各向异性超弹性本构模型准确性高、易于实现数值开发,在一定程度上能够为固体火箭发动机的装药结构完整性分析提供理论依据.     

高分子材料平面应变拉伸变形局部化

将基于应变软化玻璃状高分子材料微观特征建立的ＢＰＡ８－链分子网络模型引入ＵｐｄａｔｉｎｇＬａｇｒａｎｇｅ有限元方法，建立了适于变形局部化分析的大变形弹塑性有限元驱动应力法．在此基础上，数值模拟了初始各向同性高分子材料平面应变拉伸变形局部化的传播过程．探讨了ＢＰＡ模型对具有加工硬化特性的结晶性高分子材料变形分析的适应性；分析了局部化传播过程中颈缩截面的非均匀应力三轴效应；最后，讨论了网格尺寸以及初始几何不均匀性对颈缩扩散以及应力三轴效应的影响     

显式方法精确模拟形状记忆聚合物热力学行为

通过构建一个热耦合的多轴可压缩应变能函数,得到应力-应变、应力-温度和应变-温度之间的函数关系,建立形状记忆聚合物的本构方程.本文引入三个基于对数应变的不变量使得模型（i）可以模拟可压缩情况;（ii）适用于单轴拉伸和等双轴拉伸至少两个基准实验;（iii）多轴有效.通过显式方法（i）给出自由能和熵的具体表达,证明模型热力学定律;（ii）给出应变-应力,温度-应力以及,温度-应变的形函数具体表达.多轴模型在特定的情况下可以自动退化到各自的单轴情况. 通过调节形函数的参数,最终得到的模型结果和实验结果能够精确匹配.新方法建立的本构模型得到的结果能更加准确地指导形状记忆聚合物的工程设计。     

基于应变能的复杂应力场低周疲劳分析

塑性应变能使材料微观组织结构发生不可逆变化,从而引起等效宏观应力,该应力随循环加载而增大．假定材料疲劳源处破坏是由最大拉应力引起的,最大等效宏观应力与外加应力叠加达到材料本征断裂应力时形成微裂纹．微裂纹引起上述两部分应力变化,继续加载直至宏观裂纹出现,从而得到材料的疲劳寿命．本文所建立的多轴疲劳寿命公式包含材料参数、拉应力以及塑性应变能等,以上数据可通过单轴疲劳数据和有限元方法获得．通过对SM45C材料的计算验证,表明该模型对多轴随机应变加载低周疲劳寿命,具有良好的预测结果．     

Multi-axial strain-stiffening elastic potentials with energy bounds: explicit approach based on uniaxial data

According to the well-known models for rubberlike elasticity with strainstiffening effects, the unbounded strain energy is generated with the unlimitedly growing stress when the stretch approaches certain limits. Toward a solution to this issue, an explicit approach is proposed to derive the multi-axial elastic potentials directly from the uniaxial potentials. Then, a new multi-axial potential is presented to characterize the strain-stiffening effect by prescribing suitable forms of uniaxial potentials so that the strain energy is always bounded as the stress grows to infinity. Numerical examples show good agreement with a number of test data.     

Direct determination of multi-axial elastic potentials for incompressible elastomeric solids: an accurate,explicit approach based on rational interpolation

With a novel approach based on certain logarithmic invariants, we demonstrate that a multi-axial elastic potential for incompressible, isotropic rubber-like materials may be obtained directly from two one-dimensional elastic potentials for uniaxial case and simple shear case, in a sense of exactly matching finite strain data for four benchmark tests, including uniaxial extension, simple shear, bi-axial extension, and plane-strain extension. As such, determination of multi-axial elastic potentials may be reduced to that of two one-dimensional elastic potentials. We further demonstrate that the latter two may be obtained by means of rational interpolating procedures for uniaxial data and shear data displaying strain-stiffening effects. Numerical examples are presented in fitting Treloar’s data and other data.     

OBTAINING MULTI-AXIAL ELASTIC POTENTIALS FOR RUBBER-LIKE MATERIALS VIA AN EXPLICIT,EXACT APPROACH BASED ON SPLINE INTERPOLATION

An explicit, exact approach is proposed to obtain multi-axial elastic potentials for isotropic rubber-like materials undergoing large incompressible deformations. By means of two direct, explicit procedures, this approach reduces the problem of determining multi-axial poten- tials to that of determining one-dimensional elastic potentials. To this end, two one-dimensional potentials for uniaxial case and simple shear case are respectively determined via spline inter- polation and, then, the two potentials are extended to generate a multi-axial elastic potential using a novel method based on certain logarithmic invariants. Eventually, each of the multi-axial potentials will exactly match the finite strain data from four benchmark tests.     

Inhomogeneous shearing of strain-stiffening rubber-like hollow circular cylinders

The purpose of this paper is to investigate the effects of strain-stiffening for the classical problems of axial and azimuthal shearing of a hollow circular cylinder composed of an incompressible isotropic non-linearly elastic material. For some specific strain-energy densities that give rise to strain-stiffening in the stress–stretch response, the stresses and resultant axial forces are obtained in explicit closed form. While such results are well known for classical constitutive models such as the Mooney–Rivlin and neo-Hookean models, our main focus is on materials that undergo severe strain-stiffening in the stress–stretch response. In particular, we consider in detail two phenomenological constitutive models that reflect limiting chain extensibility at the molecular level and involve constraints on the deformation. The amount of shearing that tubes composed of such materials can sustain is limited by the constraint. Numerical results are also obtained for an exponential strain-energy that exhibits a less abrupt strain-stiffening effect. Potential applications of the results to the biomechanics of soft tissues are indicated.     

Optimization of the strain energy density in linear anisotropic elasticity

The problem considered here is that of extremizing the strain energy density of a linear anisotropic material by varying the relative orientation between a fixed stress state and a fixed material symmetry. It is shown that the principal axes of stress must coincide with the principal axes of strain in order to minimize or maximize the strain energy density in this situation. Specific conditions for maxima and minima are obtained. These conditions involve the stress state and the elastic constants. It is shown that the symmetry coordinate system of cubic symmetry is the only situation in linear anisotropic elasticity for which a strain energy density extremum can exist for all stress states. The conditions for the extrema of the strain energy density for transversely isotropic and orthotropic materials with respect to uniaxial normal stress states are obtained and illustrated with data on the elastic constants of some composite materials. Not surprisingly, the results show that a uniaxial normal stress in the grain direction in wood minimizes the strain energy in the set of all uniaxial stress states. These extrema are of interest in structural and material optimization.     

An alternative form of energy density demonstrating the severe strain-stiffening in thin spherical and cylindrical shells

The present article investigates an elastic instability phenomenon for internally pressurized spherical thin balloons and thin cylindrical tubes composed of incompressible hyperelastic material. A mathematical model is formulated by proposing a new strain energy density function. In the family of limited elastic materials, many material models exhibit strain-stiffening. However, they fail to predict severe strain-stiffening in a moderate range of deformations in the stress-strain relations. The proposed energy function contains three material parameters and shows substantially improved stain stiffening properties than the limited elastic material models. The model is further applied to explore the elastic instability phenomenon in spherical and cylindrical shells. The findings are compared with other existing models and validated with experimental results. The model shows better agreement with experimental results and exhibits a substantial strain-stiffening effect than the current models.     

Plane strain bending of strain-stiffening rubber-like rectangular beams

This paper is concerned with investigation of the effects of strain-stiffening for the classical problem of plane strain bending by an end moment of a rectangular beam composed of an incompressible isotropic nonlinearly elastic material. For a variety of specific strain-energy densities that give rise to strain-stiffening in the stress–stretch response, the stresses and resultant moments are obtained explicitly. While such results are well known for classical constitutive models such as the Mooney-Rivlin and neo-Hookean models, our primary focus is on materials that undergo severe strain-stiffening in the stress–stretch response. In particular, we consider in detail two phenomenological constitutive models that reflect limiting chain extensibility at the molecular level and involve constraints on the deformation. The amount of bending that beams composed of such materials can sustain is limited by the constraint. Potential applications of the results to the biomechanics of soft tissues are indicated.     

Further study of rubber-like elasticity： elastic potentials matching biaxial data＊

By virtue of the rational interpolation procedure and logarithmic strain, a direct approach is proposed to obtain elastic potentials that exactly match uniaxial data and shear data for elastomers. This approach reduces the determination of multiaxial elastic potentials to that of two one-dimensional potentials, thus bypassing usual cumbersome procedures of identifying a number of unknown parameters. Predictions of the suggested potential are derived for a general biaxial stretch test and compared with the classical data given by Rivlin and Saunders(Rivlin, R. S. and Saunders, D. W. Large elastic deformation of isotropic materials. VII: experiments on the deformation of rubber. Phill. Trans. Royal Soc. London A, 243, 251–288(1951)). Good agreement is achieved with these extensive data.     

超弹性材料本构关系的最新研究进展

超弹性材料是工程实际中的常用材料, 具有在外力作用下经历非常大变形、在外力撤去后完全恢复至初始状态的特征. 超弹性材料是典型的非线性弹性材料, 其性能可通过材料的应变能函数予以表征. 近几十年来, 围绕应变能函数形式的构造, 已提出许多超弹性材料本构关系研究的数学模型和物理模型, 但适用于多种变形模式和全变形范围的完全本构关系仍是该领域期待解决的重要问题. 本文从3个不同角度, 对超弹性材料本构关系研究的最新进展进行了总结和分析: (1)不同体积变化模式, 包含不可压与可压两种; (2)多变形模式, 包含单轴拉伸、剪切、等双轴以及复合拉剪等多个种类; (3)全范围变形程度, 包含小变形、中等变形到较大变形范围. 超弹性材料本构关系研究的最新进展表明, 为了全面描述具体材料的实验数据并在实际问题中应用超弹性材料, 需要建立适合于多种变形模式和全变形范围的可压超弹性材料的完全本构关系. 对实际超弹性材料完全本构关系的建立及可压超弹性材料应变能函数的构造, 笔者还提出了相应的实施步骤和研究方法.      

STATE OF THE ART OF CONSTITUTIVE RELATIONS OF HYPERELASTIC MATERIALS 1)

超弹性材料是工程实际中的常用材料, 具有在外力作用下经历非常大变形、在外力撤去后完全恢复至初始状态的特征. 超弹性材料是典型的非线性弹性材料, 其性能可通过材料的应变能函数予以表征. 近几十年来, 围绕应变能函数形式的构造, 已提出许多超弹性材料本构关系研究的数学模型和物理模型, 但适用于多种变形模式和全变形范围的完全本构关系仍是该领域期待解决的重要问题. 本文从3个不同角度, 对超弹性材料本构关系研究的最新进展进行了总结和分析: (1)不同体积变化模式, 包含不可压与可压两种; (2)多变形模式, 包含单轴拉伸、剪切、等双轴以及复合拉剪等多个种类; (3)全范围变形程度, 包含小变形、中等变形到较大变形范围. 超弹性材料本构关系研究的最新进展表明, 为了全面描述具体材料的实验数据并在实际问题中应用超弹性材料, 需要建立适合于多种变形模式和全变形范围的可压超弹性材料的完全本构关系. 对实际超弹性材料完全本构关系的建立及可压超弹性材料应变能函数的构造, 笔者还提出了相应的实施步骤和研究方法.