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

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

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

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

横观各向同性材料三维裂纹问题的数值分析

严格从三维横观各向同性材料弹性空间问题的Green函数出发,采用Hadamard有限部积分概念,导出了三维状态下单位位移间断(位错)集度的基本解.在此基础上,将三维任意形状的片状裂纹问题归结为求解-组以未知位移间断表示的超奇异积分方程;并给出了边界元离散形式.对方程中出现的超奇异积分,采用了Had-alnard定义的有限部积分来处理.论文最后给出了若干典型片状裂纹问题的数值算例,数值结果表明了本文方法是非常有效的.     

立方晶粒任意集合下的格林函数和有效弹性刚度张量

为了推导多晶体材料的有效弹性刚度张量,给出立方晶粒任意集合的格林函数封闭但近似的表达式,该格林函数表达式包含三个单晶弹性常数和多晶体材料五个织构系数,它考虑取向分布函数的影响直至织构系数的线性项,它适用于弱织构多晶体材料或具有弱各向异性晶粒的多晶体材料(如金属铝),它与Nishioka格林函数近似式的比较通过三个算例给出;Synge的格林函数积分式则直接通过数值计算完成,它可作为问题的精确解供参考．该文还简单介绍了多晶体材料有效弹性刚度张量的推导过程,并把所得结果和有限元计算结果进行比较。     

高阶边界元奇异积分的一种通用高效计算方法

针对三维边界元法中曲面单元上的(弱、强、超)奇异积分提出了一种通用高效的计算方法。经极坐标变换,将奇异积分转化为常规积分;采用数值方法计算Cauchy主值积分和Hadamard有限项积分系数;引入保角变换和反曲变换消除因单元畸形或因积分点靠近单元边界而引起的周向积分奇异性。该方法可以统一处理(弱、强、超)奇异积分,并且只需要知道核函数的奇异阶数和少数几个点上的被积函数值,不依赖于积分和函数的具体选取;所需的积分点少,精度高,并且受单元畸形程度影响较小,稳定性好。采用该方法计算了声学和弹性力学中的典型奇异积分,并结合二阶Nystrm方法求解了弹性力学的边界积分方程,验证了方法的高精度和高效性。本文数值积分程序可向作者索取。     

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

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

横观各向同性材料的三维断裂力学问题

从三维横观各向同性材料弹性力学理论出发, 使用Hadamard有限部积分概念, 导出了三维状态下单位位移间断(位错)集度的基 本解. 在此基础上, 进一步运用极限理论, 将任意载荷作用下, 三维无限大横观各向 同性材料弹性体中, 含有一个位于弹性对称面内的任意形状的片状裂纹问题, 归结为求 解一组超奇异积分方程的问题. 通过二维超奇异积分的主部分析方法, 精确地求得了裂纹前沿光滑点附近的应力奇异指数和奇异应力场, 从而找到了以裂纹表面位移间断表示的应力强度因子表达式及裂纹局部扩展所提供 的能量释放率. 作为以上理论的实际应用,最后给出了一个圆形片状裂纹问题 的精确解例和一个正方形片状裂纹问题的数值解例. 对受轴对称法向均布载荷作用下圆形片状裂纹问题, 讨论了超奇异积分方程的精确求解方法, 并获得了位移间断和应力强度因子的封闭解, 此结果与现有理论解完全一致.     

不可压缩超弹性止水材料的粘弹性计算方法研究

发展改进了研究不可压缩超弹性材料的Mooney-Rivlin公式,把材料参数假设成是随时间变化的粘弹性函数,并讨论了粘弹性函数的推导方法,通过工程实例进行了实验计算,并将计算结果与应力松弛试验进行了比较.结果表明,用改进Mooney-Rivlin公式可以简便有效地计算该类材料的粘弹性问题,为解决工程实际问题提供了新的途径.     

三维有限变形弹性问题的分析方法

至今还未见到用通常的应力函数和位移函数分析三维有限变形弹性问题的报导。利用Hasegawa的工作和Adkins的摄动法,本文将位移函数用于求解表面力或体积力作用下的有限变形轴对称弹性问题,提出一个分析可压缩和不可压缩材料三维弹性问题的新的解析法,并用两个简单例子验证了这种分析方法。     

湿热作用下热超弹性材料在电子封装中的分层失效问题

研究了具有Gent-Thomas特征的热超弹性材料构成的高聚物电子封装件在回流焊过程中由于吸湿所引发的蒸汽压力以及由于材料的热失配引发的热应力共同作用下而导致的“爆米花”式的分层失效问题.利用超弹性材料空穴生成和增长的理论给出了此类封装材料在回流焊过程中孔穴的生成及增长与蒸汽压力和热应力之间的解析关系.分析结果表明,当...     

On the surface instability of a highly elastic half-space

The phenomenon of surface instability of an isotropic half-space under biaxial plane stress is studied for compressible elastic materials in finite strain. Euler's method is used to derive the general form of the stability criterion, and analytical details are exhibited by special application to the class of hyperelastic Hadamard materials in two complementary cases: (i) the full solution is derived for the compressible, neo-Hookean members, and (ii) the plane deformation solution is provided for every isotropic, elastic material and specific results are presented for the full Hadamard class. Results appropriate to incompressible Mooney-Rivlin materials are herein obtained as special limit cases. Several theorems are established and some of the conclusions are illustrated graphically.     

On instabilities of single-integral constitutive equations for viscoelastic liquids

Various types of instabilities are exposed in this paper for time-strain separable single-integral viscoelastic constitutive equations (CE's). They were distinguished into two groups and defined as Hadamard and dissipative type of instabilities. As for the Hadamard-type, previously obtained criteria are found to be necessary only. They are necessary and sufficient only for thermodynamic stability. Improved, stricter Hadamard stability criteria are described briefly in this paper, and then applied to study of stability of several CE's. It is shown that the Currie potential with the K-BKZ equation and the model proposed by Papanastasiou et al. are Hadamard unstable. In the case of dissipative stability, the necessary and sufficient condition for stress boundedness in any regular flow with a given history, is proved. Then, this criterion was applied to the neoHookean, Mooney, and Yen and McIntire specifications of the general K-BKZ model, to exhibit unbounded solutions. In addition, Larson-Monroe potential which is later proved to be Hadamard unstable but satisfies the above criterion of boundedness, is shown to have unstable decreasing branch in steady simple shear flow. At present, to the authors' knowledge, there is no viscoelastic single-integral CE of factorable type proposed in the literature which can satisfy all the Hadamard and dissipative stability criteria.     

On non-linear stability in unconstrained non-linear elasticity

A method of obtaining a full three-dimensional non-linear Hadamard stability analysis of inhomogeneous deformations of arbitrary, unconstrained, hyperelastic materials is presented. The analysis is an extension of that given by Chen and Haughton (Proc. Roy. Soc. London A 459 (2003) 137) for two-dimensional incompressible problems. The process that we present replaces the second variation condition expressed as an integral involving a quadratic in three arbitrary perturbations, with an equivalent sixth-order system of ordinary differential equations. The positive definiteness condition is thereby reduced to the simple numerical evaluation of zeros of a well-behaved function. The general theory is illustrated by applying it to the problem of the inflation of a thick-walled spherical shell. The present analysis provides a simpler alternative approach to bifurcation problems approached by using the incremental equations of non-linear elasticity.     

Finite-time stability of Hadamard fractional differential equations in weighted Banach spaces

Nonlinear Dynamics - The main purpose of this paper is to investigate the finite-time stability of Hadamard fractional differential equations (HFDEs). Firstly, the standard definitions of...     

Estimation of ultimate safe loads in elastic stability theory

Ultimate safe load estimates are derived for general isotropic materials. Conditions sufficient to assure stability of equilibrium in the dead load traction boundary value problem are obtained from an energy type criterion by application of certain inequalities and by subsequent separation of the incremental shear strain components from the others. Specific estimates are provided for stability of the undistorted states of every isotropic elastic material, and it is shown that the inequalities based upon the estimates derived here imply those given elsewhere; in the natural state of vanishing stress and strain the conclusions are coincident with the classical inequalities for the shear and bulk moduli. General sufficient conditions for stability of equilibrium of an incompressible Mooney-Rivlin body are derived. Simple shear and simple extension analyses for isotropic Hadamard elastic materials reveal that several classical relations can be obtained for this class of materials in the theory of finite strain, and explicit formulae for determination of certain elastic constants and familiar moduli are provided. Sufficient conditions for stability for every such Hadamard elastic material are obtained by special application of the general formulae developed earlier in the work, and particular results for a compressible Mooney-Rivlin type material are given. It is found that these include as a special case our earlier results for the familiar incompressible Mooney-Rivlin body. It happens also that several inequalities usually obtained from certain requirements that assure resonable material response to the loading, but which otherwise have been unrelated to stability, appear here as natural sufficient conditions for stability in the dead load problem. Finally, several of the results are applied to study the problem of stability of an Euler column, and physical implications of the analyses are discussed.

---

Zusammenfassung Für allgemeine isotrope Materialien werden völlig sichere Lastabschätzungen hergeleitet. Bei Randwertproblemen nicht mitgehender Lasten werden hinreichende Bedingungen zur Gewährleistung der Stabilität des Gleichgewichts von einem Energie-Kriterium durch Anwendung verschiedener Ungleichungen und durch nachfolgende Trennung des Zuwachses der Schubspannungskomponenten von den übrigen erhalten. Für die Stabilität unverformter Zustände beliebiger isotroper elastischer Materialien werden spezielle Abschätzungen geliefert. Dabei wird gezeigt, daß die Ungleichungen, die sich aus den hergeleiteten Abschätzungen ergeben, diejenigen anderer Autoren enthalten. Die Ergebnisse stimmen für den spannungs- und verformungsfreien Zustand mit den klassischen Ungleichungen für Schub- und Kompressionsmodul überein. Es werden allgemeine hinreichende Bedingungen für die Stabilität des Gleichgewichts eines inkompressiblen Körpers vom Mooney-Rivlin-Typ abgeleitet. Einfache Schub- und Dehnungsberechnungen für isotrope Hadamard-elastische Materialien zeigen, daß einige klassische Beziehungen in der Theorie endlicher Dehnungen für diese Materialien erhalten werden können, und es werden explizite Formeln für die Bestimmung gewisser elastischer Konstanten und üblicher Moduln angegeben. Hinreichende Stabilitätskriterien werden durch eine spezielle Anwendung der allgemeinen Gleichungen gewonnen und Ergebnisse für ein kompressibles Mooney-Rivlin-Material dargestellt. Es zeigt sich, daß diese unsere früheren Ergebnisse für den inkompressiblen Mooney-Rivlin-Körper als Spezialfall enthalten. Es tritt auch der Fall auf, daß einzelne Ungleichungen, die gewöhnlich aus gewissen Forderungen zur Gewährleistung einer sinnvollen Reaktion des Materials auf die Belastung folgen, die aber auf der anderen Seite keinen Bezug zur Stabilität besitzen, hier als natürliche hinreichende Stabilitätskriterien bei Problemen nicht mitgehender Lasten auftreten. Schließlich werden einige Ergebnisse zur Untersuchung der Stabilität eines Euler-Stabes herangezogen und die physikalische Bedeutung der Rechnungen diskutiert.

     

On the mechanical modeling of anisotropic biological soft tissue and iterative parallel solution strategies

Biological soft tissues appearing in arterial walls are characterized by a nearly incompressible, anisotropic, hyperelastic material behavior in the physiological range of deformations. For the representation of such materials we apply a polyconvex strain energy function in order to ensure the existence of minimizers and in order to satisfy the Legendre–Hadamard condition automatically. The 3D discretization results in a large system of equations; therefore, a parallel algorithm is applied to solve the equilibrium problem. Domain decomposition methods like the Dual-Primal Finite Element Tearing and Interconnecting (FETI-DP) method are designed to solve large linear systems of equations, that arise from the discretization of partial differential equations, on parallel computers. Their numerical and parallel scalability, as well as their robustness, also in the incompressible limit, has been shown theoretically and in numerical simulations. We are using a dual-primal FETI method to solve nonlinear, anisotropic elasticity problems for 3D models of arterial walls and present some preliminary numerical results.     

Analysis of material instabilities in inelastic solids by incremental energy minimization and relaxation methods: evolving deformation microstructures in finite plasticity

We propose an approach to the definition and analysis of material instabilities in rate-independent standard dissipative solids at finite strains based on finite-step-sized incremental energy minimization principles. The point of departure is a recently developed constitutive minimization principle for standard dissipative materials that optimizes a generalized incremental work function with respect to the internal variables. In an incremental setting at finite time steps this variational problem defines a quasi-hyperelastic stress potential. The existence of this potential allows to be recast a typical incremental boundary-value problem of quasi-static inelasticity into a principle of minimum incremental energy for standard dissipative solids. Mathematical existence theorems for sufficiently regular minimizers then induce a definition of the material stability of the inelastic material response in terms of the sequentially weakly lower semicontinuity of the incremental variational functional. As a consequence, the incremental material stability of standard dissipative solids may be defined in terms of the quasi-convexity or the rank-one convexity of the incremental stress potential. This global definition includes the classical local Hadamard condition but is more general. Furthermore, the variational setting opens up the possibility to analyze the post-critical development of deformation microstructures in non-stable inelastic materials based on energy relaxation methods. We outline minimization principles of quasi- and rank-one convexifications of incremental non-convex stress potentials for standard dissipative solids. The general concepts are applied to the analysis of evolving deformation microstructures in single-slip plasticity. For this canonical model problem, we outline details of the constitutive variational formulation and develop numerical and semi-analytical solution methods for a first-level rank-one convexification. A set of representative numerical investigations analyze the development of deformation microstructures in the form of rank-one laminates in single slip plasticity for homogeneous macro-deformation modes as well as inhomogeneous macroscopic boundary-value problems. The well-posedness of the relaxed variational formulation is indicated by an independence of typical finite element solutions on the mesh-size.     

Ill-posed problems in mechanics

The notion of ill-posed initial and boundary value problems for partial differential equations was introduced by Hadamard, who also presented the first example of an ill-posed problem for a specific partial differential equation. At the same time, there are numerous examples of ill-posed problems in any field of mechanics.Hadamard and some of his successors believed that any ill-posed problem has no physical meaning and such problems should not be posed.The present paper contains several examples of ill-posed problems. We show that if one deals with an applied problem, then overcoming the ill-posedness mathematically can help one to improve the structure in practice, which justifies the study of ill-posed problems.     

Constraints, Constitutive Limits, and Instability in Finite Thermoelasticity

This paper examines all the possible types of thermomechanical constraints in finite-deformational elasticity. By a thermomechanical constraint we mean a functional relationship between a mechanical variable, either the deformation gradient or the stress, and a thermal variable, temperature, entropy or one of the energy potentials; internal energy, Helmholtz free energy, Gibbs free energy or enthalpy. It is shown that for the temperature-deformation, entropy-stress, enthalpy-deformation, and Helmholtz free energy-stress constraints equilibrium states are unstable, in the sense that certain perturbations of the equilibrium state grow exponentially. By considering the constrained materials as constitutive limits of unconstrained materials, it is shown that the instability is associated with the violation of the Legendre–Hadamard condition on the internal energy. The entropy-deformation, temperature-stress, internal energy-stress, and Gibbs free energy-deformation constraints do not exhibit this instability. It is proposed that stability of the rest state (or equivalently convexity of internal energy) is a necessary requirement for a model to be physically valid, and hence entropy-deformation, temperature-stress, internal energy-stress, and Gibbs free energy-deformation constraints are physical, whereas temperature-deformation constraints (including the customary notion of thermal expansion that density is a function of temperature only), entropy-stress constraints, enthalpy-deformation constraints, and Helmholtz free energy-stress constraints are not.     

Plate waves in Hadamard materials

The dispersion equation for elastic waves of small amplitude in a pre-stressed plate of restricted Hadamard material is derived and its solutions are investigated in detail. The implications for stability are discussed.