A constitutive model for cavitation and cavity growth in rubber-like materials under arbitrary tri-axial loading

In this paper, we attempted to construct a constitutive model to deal with the phenomenon of cavitation and cavity growth in a rubber-like material subjected to an arbitrary tri-axial loading. To this end, we considered a spherical elementary representative volume in a general Rivlin’s incompressible material containing a central spherical cavity. The kinematics proposed by [Hou, H.S., Abeyaratne, R., 1992. Cavitation in elastic and elastic-plastic solids. J. Mech. Phys. Solids 40, 571–722] was adopted in order to construct an approximate but optimal field. In order to establish a suitable constitutive law for this class of materials, we utilized the homogenisation technique that permits us to calculate the average strain energy density of the volume. The cavity growth was considered through a physically realistic failure criterion. Combination of the constitutive law and the failure criterion enables us to describe correctly the global behaviour and the damage evolution of the material under tri-axial loading. It was shown that the present models can efficiently reproduce different stress states, varying from uniaxial to tri-axial tensions, observed in experimentations. Comparison between predicted results and experimental data proves that the proposed model is accurate and physically reasonable. Another advantage is that the proposed model does not need special identification work, the initial Rivlin’s law for the corresponding incompressible material is sufficient to form the new law for the compressible material resulted from cavitation procedure.     

Simulation of damage evolution in ductile metals undergoing dynamic loading conditions

A set of constitutive equations for large rate-dependent elastic-plastic-damage materials at elevated temperatures is presented to be able to analyze adiabatic high strain rate deformation processes for a wide range of stress triaxialities. The model is based on the concepts of continuum damage mechanics. Since the material macroscopic thermo-mechanical response under large strain and high strain rate deformation loading is governed by different physical mechanisms, a multi-dissipative approach is proposed. It incorporates thermo-mechanical coupling effects as well as internal dissipative mechanisms through rate-dependent constitutive relations with a set of internal variables. In addition, the effect of stress triaxiality on the onset and evolution of plastic flow, damage and failure is discussed.Furthermore, the algorithm for numerical integration of the coupled constitutive rate equations is presented. It relies on operator split methodology resulting in an inelastic predictor-elastic corrector technique. The explicit finite element program LS-DYNA augmented by an user-defined material subroutine is used to approximate boundary-value problems under dynamic loading conditions. Numerical simulations of dynamic experiments with different specimens are performed and good correlation of numerical results and published experimental data is achieved. Based on numerical studies modified specimens geometries are proposed to be able to detect complex damage and failure mechanisms in Hopkinson-Bar experiments.     

A constitutive model dealing with damage due to cavity growth and the Mullins effect in rubber-like materials under triaxial loading

In this work, we attempted to describe the evolution of damage in rubber-like materials due to the Mullins effect and the cavity growth process. To this end we introduced two distinct internal variables into the constitutive laws; the first one essentially describes the Mullins damage and the second describes the cavity growth. The Mullins effect was considered as a continuous type of damage that can be modelled within the continuum damage theory. The cavity growth, being discontinuous at the microscopic scale, was also modelled by a continuous variable after a homogenization procedure. These analyses allow the establishment of a compressible constitutive law describing the strain-softening phenomena for rubber-like materials. In order to identify the material parameters and to verify the efficiency of the model, we carried out experimental studies involving uniaxial, biaxial, and hydrostatic tensions under monotonic and cyclic loading. Comparison between the model-predicted results and the experimental data shows that the present model can efficiently describe both the Mullins damage and the porosity evolution of rubber-like materials under triaxial monotonic or cyclic loading with a satisfactory accuracy. The proposed concept is simple and easy to apply to engineering calculations.     

Upper and lower bounds for constitutive relation of crack-weakened rock masses under dynamic compressive loads

A micro-mechanics-based model is proposed to investigate the rate-dependent constitutive relation for crack-weakened rock masses subjected to dynamic compressive loads. The present micro-mechanical model reveals that the nucleation, growth and coalescence of sliding cracks dominate the failure and macroscopic properties of crack-weakened rock masses subjected to dynamic compressive loads. The interactions among multiple parallel sliding cracks in crack-weakened rock masses subjected to dynamic compressive loads are examined asymptotically in an explicit and quantitative manner in order to reveal fully their so-called shielding and magnification effects on the stress–strain relation. Based on the micro-mechanical framework and the asymptotic analysis, analytical upper and lower bounds are proposed for the rate-relation for rock masses containing multiple rows of echelon cracks subjected to dynamic compressive loads. The factors that affect the rate-dependent properties of crack-weakened rock masses have been analyzed. The strain energy density factor approach, which is related to crack growth velocity and dynamic fracture toughness of rock material, is employed in the analysis. The rate-dependent constitutive relation of crack-weakened rock masses is derived from micro-mechanical framework and the asymptotic analysis. The closed-form explicit expression for the rate-dependent constitutive relation of rock masses containing echelon cracks subjected to dynamic compressive loads is obtained. Finally, the present model is used to analyze the complete stress–strain relation and strength for jointed rock masses at shiplock slope of the Three Gorges Dam.     

Rate-dependent transition from tensile damage to discrete fracture in dynamic brittle failure

Based on the work for a combined damage/plasticity model of geologic materials and the bifurcation analysis of material failure, an analytical framework is established to study the rate-dependent transition from continuum damage to discrete fracture in dynamic brittle failure. Because of the simple formulation, a vectorized constitutive model solver can be designed for large-scale computer simulation. A continuum tangent stiffness tensor is invoked for the tensile damage evolution such that the bifurcation analysis can be performed to identify the initiation and orientation of tensile failure. It is shown that the orientation of tensile failure is rate-independent although the limit state is rate-dependent for the rate-dependent tensile damage model. Sample problems are considered to demonstrate the features of the proposed approach.     

冲击载荷下猪后腿肌肉的横向同性本构模型

基于纤维增强复合材料连续介质力学理论及粘弹性理论,提出了猪后腿肌肉的率相关本构模型。 通过拟合以往研究中猪后腿肌肉的SHPB和SHTB实验应力应变曲线,确定了本构模型的相关参数。结果 表明:提出的本构模型既能描述猪后腿肌肉沿纤维方向的动态压缩力学性能又能描述其动态拉伸性能,理论 模型与实验模型有较好的一致性。该结果可为安全防护数值模拟提供一定的理论依据。     

橡胶材料弹性的一种新的螺旋管模型

建立统计力学模型正确描述材料微观结构与宏观力学特性之间的关系是软物质类材料的最大挑战之一,已有的橡胶材料统计模型尚存在一些不足.文章根据橡胶类材料宏观各向同性、连续均匀和不可压缩特性,结合分子链的非高斯统计模型,提出一种橡胶材料网络结构的力学特性模型.该模型将代表体元上对应点之间的传力路径用一个类螺旋管区域约束的分子链子网络来描述,螺旋管的表面随材料的宏观变形做仿射变形,分子链子网络由方向和长度随机的分子链或链段首尾链接而成,在此基础上由分子链的熵推导出描述材料宏观力学特性的本构关系.通过大量的材料测试数据对本构模型进行拟合验证,拟合结果表明该模型具有非常好的精度,并且在采用两个参数时模型具有非常高的可靠性,仅用单轴拉伸实验数据拟合模型就能准确预测全部3类实验数据.该模型使用了仿射的弯曲管假设,能从微观结构尺度上说明材料的不可压缩特性,避免了直管模型的近似性,为微观尺度的随机性和宏观的均匀性的联系提出一个新的模型.     

On the interplay between material flaws and dynamic necking

In this paper we investigate the interplay between material defects and flow localization in elastoplastic bars subjected to dynamic tension. For that task, we have developed a 1D finite difference scheme within a large deformation framework in which the material is modelled using rate-dependent J₂ plasticity. A perturbation of the initial yield stress is introduced in each node of the finite difference mesh to model localized material flaws. Numerical computations are carried out within a wide spectrum of strain rates ranging from 500 s⁻¹ to 2500 s⁻¹. On the one hand, our calculations reveal the effect of the material defects in the necking process. On the other hand, our results show that the necking inception, instead of being a random type process, is the deterministic result of the interplay between the mechanical behaviour of the material and the boundary conditions. This conclusion agrees with the experimental evidence reported by Rittel et al. [1] and Rotbaum et al. [2].     

基于近场动力学理论的准脆性材料的本构力函数的构建

近场动力学理论(PD)是基于非局部思想的连续介质力学新理论,用于研究材料破坏问题。根据准脆性材料破坏的线性和非线性的力学行为,在初始微观弹脆性材料(PMB)的本构力函数中引入了键的损伤模型,将键的断裂过程分成了线性的弹性变形阶段和非线性的损伤变形阶段,以此构建了准脆性材料的本构力函数的基本形式。以典型的准脆性材料为例构建了其本构力函数,通过在压缩载荷下对含预制不同角度单裂纹缺陷的类岩材料的裂纹扩展进行PD数值模拟仿真,裂纹起裂位置和扩展方向与试样试验结果在一定程度上保持了一致,证明了该基于近场动力学理论的典型准脆性材料的本构力函数可用于该类材料的破坏分析。     

Analysis of cavitation problem of heated elastic composite ball

The cavitation problem of a composite ball under a uniform temperature is investigated, and the ball is composed of two elastic solid materials. The nonlinear mathematical model of the problem is established with the finite logarithmic strain measure for a large geometric deformation and by the Hooke law for elastic materials. The analytic solutions in a parametric form are derived for the thermal dilatation of the composite ball with a large elastic deformation. Solution curves are given to describe the variations of the critical temperature in the cavitation with the geometric and material parameters. The bifurcation curve is also given to reveal the cavity growth after void nucleation. The numerical results for a computational example indicate that the radius of the cavity will rapidly grow above the critical temperature, and the loop stress will become infinite when void nucleation. This means that the materials near the cavity will produce a plastic deformation leading to local failure and fracture if the material of the internal ball is elastoplastic. In addition, the cavitation of the composite ball appears at a low temperature if the elastic property of the material of the internal ball is nearly uncompressible.     

An endochronic plasticity formulation for filled rubber

The nature of elastomeric material demands the consideration of finite deformations, nonlinear elasticity including damage as well as rate-dependent and rate-independent dissipative properties. While many models accounting for these effects have been refined over time to do better justice to the real behavior of rubber-like materials, the realistic simulation of the elastoplastic characteristics for filled rubber remains challenging.The classical elastic-ideal-plastic formulation exhibits a distinct yield-surface, whereas the elastoplastic material behavior of filled rubber components shows a yield-surface free plasticity. In order to describe this elastoplastic deformation of a material point adequately, a physically based endochronic plasticity model was developed and implemented into a Finite Element code. The formulation of the ground state elastic characteristics is based on Arruda and Boyce (1993) eight-chain model. The evolution of the constitutive equations for the nonlinear endochronic elastoplastic response are derived in analogy to the Bergström–Boyce finite viscoelasticity model discussed by Dal and Kaliske (2009).     

一种新的橡胶材料弹性本构模型

橡胶类材料本构关系对于科学研究和工程应用具有重要意义,但已有的橡胶模型的拟合能力和可靠性需要进一步提高.为解决此问题,本文提出了一种新的橡胶材料的各向同性、不可压缩柯西弹性模型.研究了橡胶材料本构关系的模型形式,基于平面应力变形状态,提出了一种以较大的两个伸长率为自变量、适用于一般变形状态的橡胶材料弹性本构模型形式;研究了橡胶材料在侧面受约束条件下的变形规律,分析了橡胶材料本构关系需要满足的约束条件;在此基础上,结合一个可以通过实验确定的描述平面拉伸变形状态下的橡胶材料力学特性函数,提出一种将该函数拓展为平面应力状态一般模型的方法,并给出了一个具体的函数形式,形成了一个新的不可压缩、各向同性的橡胶材料弹性本构模型.使用5组包含3种类型实验的数据和一组较全面的双轴测试数据对该模型进行了参数拟合,结果表明:该模型具有很好的拟合精度和更高的可靠性,仅用一种类型实验数据,如单轴拉伸或者平面拉伸等,也能获得较好的拟合结果.      

基于石膏材料的小结构重力坝损伤本构模型研究

大坝模型试验的模型材料力学性能是试验成功与否的关键影响因素。基于石膏模型材料基本力学性质和混凝土材料基本相似的实际情况,考虑石膏模型材料的孔隙率,在能量损伤理论模型的基础上,结合工程规范,建立了石膏模型材料的受拉损伤本构模型。同时,分析了模型材料应变率对材料峰值应力的影响,引入动应力提高系数,并给出了与应变率相关的模型材料受拉损伤本构模型。利用该本构模型模拟了动荷载下石膏模型重力坝的破坏过程,并与模型试验结果进行对比。结果表明：数值模拟得到的重力坝开裂破坏的位置、形式与模型试验结果较为接近。     

基于Khan-Huang模型高强度钢DP600率相关特性实验与本构模型研究

对汽车用高强钢DP600在10-4s-1~103s-1应变率范围内进行力学拉伸实验,在此基础上,采用唯象的方法对Khan-Huang本构模型进行修正建立DP600考虑率敏感效应的本构模型,对高应变率下试样的拉伸过程进行数值模拟验证模型有效性.结果表明,高强钢常温下具有明显的应变率敏感特性,高应变率下材料的屈服应力接近低应变率下的两倍.对Khan-Huang本构模型进行修正后可以比较准确地描述材料在不同应变率下力学行为,且参数容易确定,便于在有限元软件中实现接口,因此该模型可进一步应用于汽车的碰撞安全数值仿真中.     

Shot peening and peen forming finite element modelling – Towards a quantitative method

Peen forming is commonly used on aluminium alloys in the aerospace industry for wing skin shaping. Numerous analytical, numerical, and experimental studies have been made to better understand the effects of various peening parameters on the final material state and to predict deformed shapes, but conclusions were often limited to trends. The purpose of this study is therefore to develop and verify experimentally quantitative numerical tools for peen forming applications by studying the simple case of peening an Almen-sized AA-2024 aluminium strip in an Almen holder. The first step consisted in improving an existing random dynamic model by determining optimal dimensions. The AA-2024 target mechanical behaviour was characterized experimentally and a combined isotropic-kinematic hardening law was selected to model the material behaviour. The dynamic impact model and material constitutive law provided good prediction of peening-induced stresses in thick AA-2024 for two shot velocities. The sequence-sensitive aspect of the forming process was also investigated and a new shell-based finite element model was proposed. Numerical and experimental results for three shot velocities were compared to evaluate the validity of this numerical simulation method and promising agreement was observed.     

On criteria for dynamic adiabatic shear band propagation

In this work, we postulate the physical criterion for dynamic shear band propagation, and based on this assumption, we implement a numerical algorithm and a computation criterion to simulate initiation and propagation of dynamic adiabatic shear bands (ASBs). The physical criterion is based on the hypothesis that material inside the shear band region undergoes a dynamic recrystallization process during deformation under high temperature and high strain-rate conditions. In addition to providing a new perspective to the physics of the adiabatic shearbanding process and identifying material properties that play a crucial role in defining the material's susceptibility to ASBs, the proposed criterion is instrumental in numerical simulations of the propagation of ASBs when multi-physics models are adopted to describe and predict the complex constitutive behavior of ASBs in ductile materials. Systematic and large scale meshfree simulations have been conducted to test and validate the proposed criterion by examining the formation, propagation, and post-bifurcation behaviors of ASBs in two materials, 4340 steel and OFHC copper. The effects of heat conduction, in particular the length scale introduced by heat conduction, are also studied. The results of the numerical simulations are compared with experimental observations and a close agreement is found for various characteristic features of ASBs, such as the shear band width, speed of propagation, and maximum temperature.     

Mechanical modelling of ductile fracture

Mechanical models of material failure by void growth to coalescence are described to give a brief overview of methods applied in the analysis of ductile fracture. Approximate constitutive relations for porous ductile materials are discussed, modelling both the nucleation and growth of voids. The application of the material models is illustrated by numerical analyses for a tensile test specimen and for dynamic, ductile crack growth. Unstable void growth is a relevant mechanism in ductile materials subject to a high level of triaxial tension. The analysis of such cavitation instabilities in elastic-perfectly plastic materials is discussed for axisymmetric stress states, and the relevance to metal/ceramic components is emphasized.General Lecture presented at the 10th Italian National Congress of Theoretical and Applied Mechanics; AIMETA, Pisa, October 1990.     

Spherical indentation of incompressible rubber-like materials

In recent years, indentation tests have been proven very useful in probing mechanical properties of small volumes of materials. However, a class of materials that very little has been done in this direction is rubber-like materials (elastomers). The present work investigates the spherical indentation of incompressible rubber-like materials. The analysis is performed in the context of second-order hyperelasticity and is accompanied by finite element computations and an extensive experimental program with spherical indentors of different radii. Uniaxial tensile tests were also performed and it was found that the initial elastic modulus correlates well with the indentation response. The experiments suggest stiffer indentation response than that predicted by linear elasticity, which is somehow counter-intuitive, if the uniaxial material response is to be considered. Regarding the uniqueness of the inverse problem, that is to establish material properties from spherical indentation tests, the answer is disappointing. We prove that the inverse problem does not give unique answer regarding the constitutive relation, except for the initial stiffness.