Finite telescopic shear of a compressible hyperelastic tube

Finite telescopic shear of a compressible hyperelastic tube is considered. It is shown that solutions with isochoric deformation fields exist for a class of strain energy functions. A numerical method is proposed for the analysis of the problem when a solution with an isochoric deformation field does not exist. Numerical results obtained by using a programmable desk calculator are presented graphically for two strain energy functions.     

Second-order estimate of the macroscopic behavior of periodic hyperelastic composites: theory and experimental validation

This paper deals with some theoretical and experimental aspects of the behavior of periodic hyperelastic composites. We focus here on composites consisting of an elastomeric matrix periodically reinforced by long fibers. The paper is composed of three parts. The first part deals with the theoretical aspects of compressible behavior. The second-order theory of Ponte Castañeda (J. Mech. Phys. Solids 44 (1996) 827) is considered and extended to periodic microstructures. Comparisons with results obtained by the finite element method show that the composite behavior predicted by the present model is much more accurate for compressible than for incompressible materials. The second part deals with the extension of the method to incompressible behavior. A mixed formulation (displacement-pressure) is used which improves the accuracy of the estimate given by the model. The third part presents experimental results. The composite tested is made of a rubber matrix reinforced by steel wires. Firstly, the matrix behavior is identified with a tensile test and a shear test carried out on homogeneous samples. Secondly, the composite is tested under shearing. The experimentally measured homogenized stress is then compared with the predictions of the model.     

A composites-based hyperelastic constitutive model for soft tissue with application to the human annulus fibrosus

This paper presents a composites-based hyperelastic constitutive model for soft tissue. Well organized soft tissue is treated as a composite in which the matrix material is embedded with a single family of aligned fibers. The fiber is modeled as a generalized neo-Hookean material in which the stiffness depends on fiber stretch. The deformation gradient is decomposed multiplicatively into two parts: a uniaxial deformation along the fiber direction and a subsequent shear deformation. This permits the fiber-matrix interaction caused by inhomogeneous deformation to be estimated by using effective properties from conventional composites theory based on small strain linear elasticity and suitably generalized to the present large deformation case. A transversely isotropic hyperelastic model is proposed to describe the mechanical behavior of fiber-reinforced soft tissue. This model is then applied to the human annulus fibrosus. Because of the layered anatomical structure of the annulus fibrosus, an orthotropic hyperelastic model of the annulus fibrosus is developed. Simulations show that the model reproduces the stress-strain response of the human annulus fibrosus accurately. We also show that the expression for the fiber-matrix shear interaction energy used in a previous phenomenological model is compatible with that derived in the present paper.     

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

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

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

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

Distortion of Anisotropic Hyperelastic Solids Under Pure Pressure Loading: Compressibility, Incompressibility and Near-Incompressibility

The purpose of this note is to examine distortion during pure pressure loading for anisotropic hyperelastic solids. We contrast the corresponding issues in compressible and incompressible hyperelasticity, and then use these results to examine nearly incompressible materials. An anisotropic compressible hyperelastic solid will generally exhibit both volume change and distortion under hydrostatic pressure loading. In contrast, an incompressible hyperelastic solid—both isotropic and anisotropic—exhibits no change to its current state of deformation as the hydrostatic pressure is varied. Nearly incompressible hyperelastic materials are compressible, but approach an incompressible response in an appropriate limit. We examine this limiting process in the context of transverse isotropy. The issue arises as to how to implement a nearly incompressible version of a given truly incompressible material model. Here we examine how certain implementations eliminate distortion under pure pressure loading and why alternative implementations do not eliminate the distortion.     

Finite deformation of incompressible fiber-reinforced elastomers: A computational micromechanics approach

The in-plane finite deformation of incompressible fiber-reinforced elastomers was studied using computational micromechanics. Composite microstructure was made up of a random and homogeneous dispersion of aligned rigid fibers within a hyperelastic matrix. Different matrices (Neo-Hookean and Gent), fibers (monodisperse or polydisperse, circular or elliptical section) and reinforcement volume fractions (10–40%) were analyzed through the finite element simulation of a representative volume element of the microstructure. A successive remeshing strategy was employed when necessary to reach the large deformation regime in which the evolution of the microstructure influences the effective properties. The simulations provided for the first time “quasi-exact” results of the in-plane finite deformation for this class of composites, which were used to assess the accuracy of the available homogenization estimates for incompressible hyperelastic composites.     

多孔硅橡胶有限变形的粘弹性行为

针对孔隙度较大(孔隙度大于50%)的硅橡胶材料在有限变形时的粘弹性行为,从建立描述材料粘弹性特征的松弛函数和变形特征的应变能函数出发,提出了适合多孔隙、可压硅橡胶材料的非线性粘弹性力学行为的本构关系,松弛函数和应变能函数可解耦为等容和体积变形两部分,并引入了拟时间的概念来反映变形对材料特征时间的影响.利用硅橡胶材料的单轴压缩松弛实验与材料模型进行了对比,讨论了多孔硅橡胶的等容变形和体积变形对应力松弛的影响.     

Separation at the interface of non-linearly elastic spinning tube-rigid shaft subjected to circumferential shear

Separation at the interface of homogeneous, isotropic, compressible, hyperelastic, spinning cylindrical tube-rigid shaft subjected to circumferential shear is investigated within the context of the finite elasticity theory. The compressible, hyperelastic spinning tube with a uniform wall thickness is assumed to be tautly fitted to a rigid shaft along its inner curved surface. The outer surface of the tube is subjected to a constant uniformly distributed circumferential shearing stress while the rigid shaft is assumed to spin with an angular speed. The state when a separation occurs at the interface of the shaft and the tube is investigated. The critical values are given for slightly compressible rubbers and nearly incompressible rubbers.     

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 Deformations with the Same Principal Stretches

In the context of the finite strain theory, plane isochoric homogeneous deformations are considered. Inspired by two examples (plane elliptical and plane hyperbolic deformations), it is seen that for any such isochoric deformation the corresponding principal stretches are equal to those of simple shear provided there is a certain relation between the amount of shear of the simple shear and the parameters of the general plane deformation. Then, the link is established between any two homogeneous deformations which have identical principal stretches. It involves two rotations, one in the undeformed state and the other in the deformed state. These rotations are determined explicitly for an arbitrary isochoric homogeneous deformation and the simple shear with the same principal stretches.     

A unified exact analysis for the Poynting effects of cylindrical tubes made of Hill’s class of Hookean compressible elastic materials at finite strain

A direct, natural extension of Hooke’s law to finite strain was achieved by R. Hill in 1978, employing the notion of work-conjugate measures of stress and strain. With Seth-Hill (Doyle-Ericksen) class of finite strain measures, this extension actually defines a broad class of compressible hyperelastic materials at finite strain, each of which retains the simple linear structure of Hooke’s law as stress–strain relationship. Several known simple elasticity models at finite strain are included as its particular examples. With a novel idea of utilizing a suitable parametric variable, here we present a unified study of the free-end torsion problem (Poynting effects) of thin-walled cylindrical tubes made of the foregoing Hill’s class of Hookean type hyperelastic materials. We show that it is possible to derive a unified exact solution to the nonlinear coupling equations relating the torque (the shear stress) and the controlling deformation quantities including, in particular, the axial length change. Discussions and comparisons concerning various Hookean type elasticity models are made based on the exact solution obtained.     

Radial deflections of thin precompressed cylindrical rubber bush mountings

For an isotropic incompressible hyperelastic Varga material the plane stress (membrane) theory of thin sheets is employed to formulate the load-deflection relation for small superimposed radial deflections of a cylindrical rubber bush which is precompressed by a large uniform radial inflation. The Varga material is a prototype for rubber over a limited range of deformation and the load-deflection relation obtained provides an extreme lower bound to the practical situation. Moreover this relation complements existing results for cylindrical rubber bushes so that now at least some assessment can be made of the effect of precompression on the radial mode of deflection for bushes of finite length. Typical numerical values are given and are contrasted with corresponding values obtained from existing plane strain radial load-deflection relations for long precompressed cylindrical rubber bush mountings.     

Finite strain––isotropic hyperelasticity

This paper presents a strain energy density for isotropic hyperelastic materials. The strain energy density is decomposed into a compressible and incompressible component. The incompressible component is the same as the generalized Mooney expression while the compressible component is shown to be a function of the volume invariant J only. The strain energy density proposed is used to investigate problems involving incompressible isotropic materials such as rubber under homogeneous strain, compressible isotropic materials under high hydrostatic pressure and volume change under uniaxial tension. Comparison with experimental data is good. The formulation is also used to derive a strain energy density expression for compressible isotropic neo-Hookean materials. The constitutive relationship for the second Piola–Kirchhoff stress tensor and its physical counterpart, involves the contravariant Almansi strain tensor. The stress stretch relationship comprises of a component associated with volume constrained distortion and a hydrostatic pressure which results in volumetric dilation. An important property of this constitutive relationship is that the hydrostatic pressure component of the stress vector which is associated with volumetric dilation will have no shear component on any surface in any configuration. This same property is not true for a neo-Hookean Green’s strain–second Piola–Kirchhoff stress tensor formulation.     

On Compressible Materials Capable of Sustaining Axisymmetric Shear Deformations. Part 4: Helical Shear of Anisotropic Hyperelastic Materials

Conditions on the form of the strain energy function in order that homogeneous, compressible and isotropic hyperelastic materials may sustain controllable static, axisymmetric anti-plane shear, azimuthal shear, and helical shear deformations of a hollow, circular cylinder have been explored in several recent papers. Here we study conditions on the strain energy function for homogeneous and compressible, anisotropic hyperelastic materials necessary and sufficient to sustain controllable, axisymmetric helical shear deformations of the tube. Similar results for separate axisymmetric anti-plane shear deformations and rotational shear deformations are then obtained from the principal theorem for helical shear deformations. The three theorems are illustrated for general compressible transversely isotropic materials for which the isotropy axis coincides with the cylinder axis. Previously known necessary and sufficient conditions on the strain energy for compressible and isotropic hyperelastic materials in order that the three classes of axisymmetric shear deformations may be possible follow by specialization of the anisotropic case. It is shown that the required monotonicity condition for the isotropic case is much simpler and less restrictive. Restrictions necessary and sufficient for anti-plane and rotational shear deformations to be possible in compressible hyperelastic materials having a helical axis of transverse isotropy that winds at a constant angle around the tube axis are derived. Results for the previous case and for a circular axis of transverse isotropy are included as degenerate helices. All of the conditions derived here have essentially algebraic structure and are easy to apply. The general rules are applied in several examples for specific strain energy functions of compressible and homogeneous transversely isotropic materials having straight, circular, and helical axes of material symmetry.     

Shear-band bifurcation of an isotropic compressible hyperelastic solid

This paper concerns shear-band bifurcations from the homogeneous finite plane deformation of an isotropic compressible hyperelastic solid. The governing equations for the incremental plane deformation superposed on the initial finite deformation are derived and then the equilibrium equations in terms of incremental displacements are classified into the elliptic type, parabolic type, etc. From this classification follows a restriction which should be placed on the strain-energy function in order that the equilibrium equations may be either elliptic or parabolic for all principal stretches. For the hyperelastic solid complying with this restriction, the condition for the shear-band bifurcation is obtained. Finally, the incremental displacement field of an infinite series of shear bands in a slab sandwiched between slippery rigid layers is established.     

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

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

Three-dimensional micromechanical modeling of voided polymeric materials

A three-dimensional micromechanical unit cell model for particle-filled materials is presented. The cell model is based on a Voronoi tessellation of particles arranged on a body-centered cubic (BCC) array. The three-dimensionality of the present cell model enables the study of several deformation modes, including uniaxial, plane strain and simple shear deformations, as well as arbitrary principal stress states.The unit cell model is applied to studies on the micromechanical and macromechanical behavior of rubber-toughened polycarbonate. Different load cases are examined, including plane strain deformation, simple shear deformation and principal stress states. For a constant macroscopic strain rate, the different load cases show that the macroscopic flow strength of the blend decreases with an increase in void volume fraction, as expected. The main mechanism for plastic deformation is broad shear banding across inter-particle ligaments. The distributed nature of plastic straining acts to reduce the amount of macroscopic strain softening in the blend as the initial void volume fraction is increased. In the case of plane strain deformation, the plastic flow is observed to initiate across inter-particle ligaments in the direction of constraint. This particular mode of deformation could not have been captured using a two-dimensional, plane strain idealization of cylindrical voids in a matrix.The potential for localized crazing and/or cavitation in the matrix is addressed. It is observed that the introduction of voids acts to relieve hydrostatic stress in the matrix material, compared to the homopolymer. It is also seen that the predicted peak hydrostatic stress in the matrix is higher under plane strain deformation than under triaxial tension (with equal lateral stresses), for the same macroscopic stress triaxiality.The effect of void volume fraction on the macroscopic uniaxial tension behavior of the different blends is examined using a Considère construction for dilatant materials. The natural draw ratio was predicted to decrease with an increase in void volume fraction.     

A shearable and thickness stretchable finite strain beam model for soft structures

Soft materials and structures have recently attracted lots of research interests as they provide paramount potential applications in diverse fields including soft robotics, wearable devices, stretchable electronics and biomedical engineering. In a previous work, an Euler–Bernoulli finite strain beam model with thickness stretching effect was proposed for soft thin structures subject to stiff constraint in the width direction. By extending that model to account for the transverse shear effect, a Timoshenko-type finite strain beam model within the plane-strain context is developed in the present work. With some kinematic hypotheses, the finite deformation of the beam is analyzed, constitutive equations are deduced from the theory of finite elasticity, and by employing the standard variational method, the equilibrium equations and associated boundary conditions are derived. In the limit of infinitesimal strain, the new model degenerates to the classical extensible and shearable elastica model. The corresponding incremental equilibrium equations and associated boundary conditions are also obtained. Based on the new beam model, analytical solutions are given for simple deformation modes, including uniaxial tension, simple shear, pure bending, and buckling under an axial load. Furthermore, numerical solution procedures and results are presented for cantilevered beams and simply supported beams with immovable ends. The results are also compared with the previously developed finite strain Euler–Bernoulli beam model to demonstrate the significance of transverse shear effect for soft beams with a small length-to-thickness ratio. The developed beam model will contribute to the design and analysis of soft robots and soft devices.