Historical review of internal state variable theory for inelasticity

A review of the development and the usages of internal state variable (ISV) theory are presented in this paper. The history of different developments leading up the formulation of the watershed paper by Coleman and Gurtin is discussed. Following the Coleman and Gurtin thermodynamics, different researchers have employed the ISV theory for dislocations, creep, continuum damage mechanics (CDM), unified-creep-plasticity (UCP), polymers, composites, biomaterials, particulate materials, multiphase and multiphysics materials, materials processing, multiscale modeling, integrating materials science (structure–property relations) into applied mechanics formulations, and design optimization under uncertainty for use in practical engineering applications.     

Constitutive modeling of orthotropic sheet metals by presenting hardening-induced anisotropy

An essential work on the constitutive modeling of rolled sheet metals is the consideration of hardening-induced anisotropy. In engineering applications, we often use testing results of four specified experiments, three uniaxial-tensions in rolling, transverse and diagonal directions and one equibiaxial-tension, to describe the anisotropic features of rolled sheet metals. In order to completely take all these experimental results, including stress-components and strain-ratios, into account in the constitutive modeling for presenting hardening-induced anisotropy, an appropriate yield model is developed. This yield model can be characterized experimentally from the offset of material yield to the end of material hardening. Since this adaptive yield model can directly represent any subsequent yielding state of rolled sheet metals without the need of an artificially defined “effective stress”, it makes the constitutive modeling simpler, clearer and more physics-based. This proposed yield model is convex from the initial yield state till the end of strain-hardening and is well-suited in implementation of finite element programs.     

用ABAQUS软件处理管土相互作用中的接触面问题

采用ABAQUS软件处理管土相互作用中的接触面问题,利用ABAQUS软件中的主控-从属接触算法,使管道和海床形成一个接触对;并且建立了管土系统有限单元模型．海床土体分别采用非线性弹性模型、多孔弹性模型、Ramberg-Osgood塑性模型对管土系统进行计算通过分析计算,得到了管道沉降量与管重间的关系,以及由于管道沉陷而形成的土体楔形,土体楔形的存在,增加了管道的稳定性．计算结果和有关试验结果相符合,说明采用该软件进行管土相互作用问题分析是可行的．     

Modelling of the interface between a thin film and a substrate within a strain gradient plasticity framework

Interfaces play an important role for the plastic deformation at the micron scale. In this paper, two types of interface models for isotropic materials are developed and applied in a thin film analysis. The first type, which can also be motivated from dislocation theory, assumes that the plastic work at the interface is stored as a surface energy that is linear in plastic strain. In the second model, the plastic work is completely dissipated and there is no build-up of a surface energy. Both formulations introduce one length scale parameter for the bulk material and one for the interface, which together control the film behaviour. It is demonstrated that the two interface models give equivalent results for a monotonous, increasing load. The combined influence of bulk and interface is numerically studied and it is shown that size effects are obtained, which are controlled by the length scale parameters of bulk and interface.     

Constitutive modeling and the trousers test for fracture of rubber-like materials

The classical constitutive modeling of incompressible hyperelastic materials such as vulcanized rubber involves strain-energy densities that depend on the first two invariants of the strain tensor. The most well-known of these is the Mooney-Rivlin model and its specialization to the neo-Hookean form. While each of these models accurately predicts the mechanical behavior of rubber at moderate stretches, they fail to reflect the severe strain-stiffening and effects of limiting chain extensibility observed in experiments at large stretch. In recent years, several constitutive models that capture the effects of limiting chain extensibility have been proposed. Here we confine attention to two such phenomenological models. The first, proposed by Gent in 1996, depends only on the first invariant and involves just two material parameters. Its mathematical simplicity has facilitated the analytic solution of a wide variety of basic boundary-value problems. A modification of this model that reflects dependence on the second invariant has been proposed recently by Horgan and Saccomandi. Here we discuss the stress response of the Gent and HS models for some homogeneous deformations and apply the results to the fracture of rubber-like materials. Attention is focused on a particular fracture test, namely the trousers test where two legs of a cut specimen are pulled horizontally apart. It is shown that the cut position plays a key role in the fracture analysis, and that the effect of the cut position depends crucially on the constitutive model employed. For stiff rubber-like or biological materials, it is shown that the influence of the cut position is diminished. In fact, for linearly elastic materials, the critical driving force for fracture is independent of the cut position. It is also shown that the limiting chain extensibility models predict finite fracture toughness as the cut position approaches the edge of the specimen whereas classical hyperelastic models predict unbounded toughness in this limit. The results are relevant to the structural integrity of rubber components such as vibration isolators, vehicle tires, earthquake bearings, seals and flexible joints.     

The rational mechanics and thermodynamics of polymeric fluids based upon the concept of a variable relaxed state

The conceptual framework of polymer continuum mechanics based upon Eckart's idea of a variable relaxed state is developed. No constitutive models are explicitly used. The theory admits four constitutive functions only, the scalar specific internal energy, the vectorial heat flux, and two tensorial fluxes representing non-elastic stress and flow (slippage). The non-linearity of the constitutive relations includes self-induced anisotropy (Leonov) with Reiner-Rivlin's equation representing a special example for this. — The effectiveness of this non-linear theory is demonstrated by treating elongational flows of polymer melts.     

Micromechanical modeling of the elasto-viscoplastic behavior of semi-crystalline polymers

A micromechanically based constitutive model for the elasto-viscoplastic deformation and texture evolution of semi-crystalline polymers is developed. The model idealizes the microstructure to consist of an aggregate of two-phase layered composite inclusions. A new framework for the composite inclusion model is formulated to facilitate the use of finite deformation elasto-viscoplastic constitutive models for each constituent phase. The crystalline lamellae are modeled as anisotropic elastic with plastic flow occurring via crystallographic slip. The amorphous phase is modeled as isotropic elastic with plastic flow being a rate-dependent process with strain hardening resulting from molecular orientation. The volume-averaged deformation and stress within the inclusions are related to the macroscopic fields by a hybrid interaction model. The uniaxial compression of initially isotropic high density polyethylene (HDPE) is taken as a case study. The ability of the model to capture the elasto-plastic stress-strain behavior of HDPE during monotonic and cyclic loading, the evolution of anisotropy, and the effect of crystallinity on initial modulus, yield stress, post-yield behavior and unloading-reloading cycles are presented.     

A constitutive model for the formation of martensite in austenitic steels under large strain plasticity

A constitutive model for diffusionless phase transitions in elastoplastic materials undergoing large deformations is developed. The model takes basic thermodynamic relations as its starting point and the phase transition is treated through an internal variable (the phase fractions) approach. The usual yield potential is used together with a transformation potential to describe the evolution of the new phase. A numerical implementation of the model is presented, along with the derivation of a consistent algorithmic tangent modulus. Simulations based on the presented model are shown to agree well with experimental findings. The proposed model provides a robust tool suitable for large-scale simulations of phase transformations in austenitic steels undergoing extensive deformations, as is demonstrated in simulations of the necking of a bar under tensile loading and also in simulations of a cup deep-drawing process.     

On the definitions of failure and critical states in hypoplasticity

A granular body is said to be at failure or in a critical state if the stress state does not change while the body is continuously deformed. Within the framework of hypoplasticity, such states, generally called stationary states,are conventionally defined by the condition that an objective (the Jaumann) stress rate vanishes. However, not all stationary states attained under monotonic deformation lie within the scope of this definition. Simple shear is an example. In fact, stationary states are characterized by zero material time derivative of the stress tensor rather than zero Jaumann rate. In the present paper, we give a generalized definition of stationarity by the condition of zero material time derivative of the stress tensor. The new definition extends the set of possible stationary states and includes those which are not covered by the previous definition. Stationary states are analysed quantitatively using calibrated hypoplastic equations. It is shown numerically that, if the norm of the spin tensor is of the same order as, or smaller than, the norm of the stretching tensor, the old definition approximates all possible sationary states with sufficient accuracy.      

Modelling long-term deformation of granular soils incorporating the concept of fractional calculus

Many constitutive models exist to characterise the cyclic behaviour of granular soils but can only simulate deformations for very limited cycles. Fractional derivatives have been regarded as one potential instrument for modelling memory-dependent phenomena. In this paper, the physical connection between the fractional derivative order and the fractal dimension of granular soils is investigated in detail.Then a modified elasto-plastic constitutive model is proposed for evaluating the long-term deformation of granular soils under cyclic loading by incorporating the concept of factional calculus. To describe the flow direction of granular soils under cyclic loading, a cyclic flow potential considering particle breakage is used. Test results of several types of granular soils are used to validate the model performance.     

基于扰动状态概念与E-B模型的粗粒料力学行为模拟

考虑到基于复杂弹塑性模型与扰动状态概念的本构模型在数值积分时的困难,以Duncan-Chang(E-B)模型描述材料在相对完整状态(RI)下的响应,用临界状态模型描述材料在完全调整状态(FA)下的力学行为。依据常规三轴试验曲线判断材料状态的转变并建议了相应的扰动因子计算方法。利用Abaqus软件实现了相应的程序代码,数值算例模拟了堆石料的三轴试验,结果表明邓肯-张模型与扰动状态概念的结合在模拟应变软化、剪缩与剪胀方面具有良好能力。在此基础上进一步模拟了某混凝土面板堆石坝的施工、蓄水、退水的过程,数值结果表明变形分布符合一般规律,且在定量上与原型观测吻合较好。     

Multiscale modeling of the plasticity in an aluminum single crystal

This paper describes a numerical, hierarchical multiscale modeling methodology involving two distinct bridges over three different length scales that predicts the work hardening of face centered cubic crystals in the absence of physical experiments. This methodology builds a clear bridging approach connecting nano-, micro- and meso-scales. In this methodology, molecular dynamics simulations (nanoscale) are performed to generate mobilities for dislocations. A discrete dislocations numerical tool (microscale) then uses the mobility data obtained from the molecular dynamics simulations to determine the work hardening. The second bridge occurs as the material parameters in a slip system hardening law employed in crystal plasticity models (mesoscale) are determined by the dislocation dynamics simulation results. The material parameters are computed using a correlation procedure based on both the functional form of the hardening law and the internal elastic stress/plastic shear strain fields computed from discrete dislocations. This multiscale bridging methodology was validated by using a crystal plasticity model to predict the mechanical response of an aluminum single crystal deformed under uniaxial compressive loading along the [4 2 1] direction. The computed strain-stress response agrees well with the experimental data.     

A macroscopic multi-mechanism based constitutive model for the thermo-mechanical cyclic degeneration of shape memory effect of NiTi shape memory alloy

A macroscopic based multi-mechanism constitutive model is constructed in the framework of irreversible thermodynamics to describe the degeneration of shape memory effect occurring in the thermo-mechanical cyclic deformation of NiTi shape memory alloys(SMAs). Three phases,austenite A, twinned martensite Mtand detwinned martensite M~d, as well as the phase transitions occurring between each pair of phases( A → M~t, M~t→ A, A → M~d,M~d→ A, and M~t→ M~d) are considered in the proposed model. Meanwhile, two kinds of inelastic deformation mechanisms, martensite transformation-induced plasticity and reorientation-induced plasticity, are used to explain the degeneration of shape memory effects of NiTi SMAs. The evolution equations of internal variables are proposed by attributing the degeneration of shape memory effect to the interaction between the three phases(A, M~t, and M~d) and plastic deformation. Finally, the capability of the proposed model is verified by comparing the predictions with the experimental results of NiTi SMAs. It is shown that the degeneration of shape memory effect and its dependence on the loading level can be reasonably described by the proposed model.     

A diffuse-interface approach to two-phase isothermal flow of a Van der Waals fluid near the critical point

A novel numerical method for simulations of isothermal, compressible two-phase flows of one fluid component near the critical point is presented on the basis of a diffuse-interface model and a Van der Waals equation of state. Because of the non-convexity of the latter, the nature of the set of governing equations is mixed hyperbolic–elliptic. This prevents the application of standard numerical methods for compressible flow. Moreover, the Korteweg capillary stress tensor, characteristic for the diffuse-interface approach, introduces third-order spatial derivatives of mass density in the Navier–Stokes equation, resulting in a dispersive behavior of the solution. Our computational method relies on a transformation of the conserved variables, which controls dispersion, stabilizes the numerical simulation and enables the use of coarser grids. A one-dimensional simulation shows that this method provides better stability and accuracy than without transformation of variables. Two- and three-dimensional simulations for isothermal liquid–vapor flows, in particular the retraction of a liquid non-spherical drop in vapor and the binary droplet collision in vapor, show the applicability of the method. The surface tension calculated from the numerical results is in good agreement with its theoretical value if the computational grid is sufficiently fine.     

On the thermomechanical modeling of shape memory alloys

In this work, a general inelastic framework for the derivation of general three-dimensional thermomechanical constitutive laws for materials undergoing phase transformations is proposed. The proposed framework is based on the generalized plasticity theory and on some basic elements from the theory of continuum damage mechanics. More specifically, a new elaborate formulation of generalized plasticity theory capable of accommodating the multiple and interacting loading mechanisms, which occur during the phase transformations, is developed. Furthermore, the stiffness variations occurring during phase transformations are taken into account by the proposed framework. For this purpose, the free energy is decomposed into elastic and inelastic parts, not in a conventional way, but in one which resembles the elastic-damage cases. Also, a rate-dependent version of the theory is provided. The concepts presented are applied for the derivation of a three-dimensional thermomechanical constitutive model for Shape Memory Alloy materials. Numerical simulations to show qualitatively the ability of the model to capture the behavior of the shape memory alloys are also presented. Furthermore, the model has been fitted to actual experimental results from the literature.     

The effects of relative density of metal foams on the stresses and deformation of beam under bending

The exact analytic solution of the pure bending beam of metallic foams is given. The effects of relative density of the material on stresses and deformation are revealed with the Triantafillou and Gibson constitutive law (TG model) taken as the analysis basis. Several examples for individual foams are discussed, showing the importance of compressibility of the cellular materials. One of the objects of this study is to generalize Hill’s solution for incompressible plasticity to the case of compressible plasticity, and a kinematics parameter is brought into the analysis so that the velocity field can be determined. The English text was polished by Yunming Chen.     

Effects of microscopic boundary conditions on the deformation behavior of small-volume metallic glasses

Although large-volume metallic glasses (MGs) are susceptible to shear localization due to their intrinsically strain-softening response, recent experiments and molecular dynamics simulations have shown that small-volume MGs samples are able to exhibit work hardening response. Here, we seek to address two issues regarding the mechanical response of small-volume MGs at low homologous temperatures from a continuum-based modeling perspective: (1) are MGs capable of exhibiting a work hardening response, and (2) what is the physical mechanism which causes its work hardening response?Along with implementing a recently-developed finite-deformation, strain gradient plasticity-based constitutive model for MGs into a self-developed finite-element code, we study the tensile response of small-volume MG samples of various sizes through finite-element simulations. Our simulations show that small-volume MG samples are capable of exhibiting a work hardening response provided the following conditions are met: (a) the sample size is small enough, and (b) the appropriate microscopic boundary conditions for the free volume are imposed on the sample.