Shape memory behaviour: modelling within continuum thermomechanics

A phenomenological material model to represent the multiaxial material behaviour of shape memory alloys is proposed. The material model is able to represent the main effects of shape memory alloys: the one-way shape memory effect, the two-way shape memory effect due to external loads, the pseudoelastic and pseudoplastic behaviour as well as the transition range between pseudoelasticity and pseudoplasticity.The material model is based on a free energy function and evolution equations for internal variables. By means of the free energy function, the energy storage during thermomechanical processes is described. Evolution equations for internal variables, e.g. the inelastic strain tensor or the fraction of martensite are formulated to represent the dissipative material behaviour. In order to distinguish between different deformation mechanisms, case distinctions are introduced into the evolution equations. Thermomechanical consistency is ensured in the sense that the constitutive model satisfies the Clausius–Duhem inequality.Finally, some numerical solutions of the constitutive equations for isothermal and non-isothermal strain and stress processes demonstrate that the various phenomena of the material behaviour are well represented. This applies for uniaxial processes and for non-proportional loadings as well.     

Growth,mass transfer,and remodeling in fiber-reinforced,multi-constituent materials

We represent a biological tissue by a multi-constituent, fiber-reinforced material, in which we consider two phases: fluid, and a fiber-reinforced solid. Among the various processes that may occur in such a system, we study growth, mass transfer, and remodeling. To us, mass transfer is the reciprocal exchange of constituents between the phases, growth is the variation of mass of the system in response to interactions with the surrounding environment, and remodeling is the evolution of its internal structure. We embrace the theory according to which these events, which lead to a structural reorganization of the system and anelastic deformations, require the introduction of balance laws, which complete the physical picture offered by the standard ones. The former are said to be non-standard. Our purposes are to determine the rates of anelastic deformation related to mass transfer and growth, and to study fiber reorientation in the case of a statistical distribution of fibers. In particular, we discuss the use of the non-standard balance laws in modeling transfer of mass, and compare our results with a formulation in which such balance laws are not introduced.     

Global view of microstructural evolution: Energetics,kinetics and dynamical systems

An evolving material structure is in a non-equilibrium state, with free energy expressed by the generalized coordinates. A global approach leads to robust computations for the generalized thermodynamic forces. Those forces drive various kinetic processes, causing dissipation at spots, along curves, surfaces and interfaces, and within volumetric regions. The actual evolution path, and therefore the final equilibrium state, is determined by the energetics and kinetics. A virtual work principle links the free energy landscape and the kinetic processes, and assigns a viscous environment to every point on the landscape. The approach leads to a dynamical system that governs the evolution of generalized coordinates. The microstructural evolution is globally characterized by a basin map in the coordinate space; and by a diversity map and a variety map in the parameter space. The control of basin boundaries raises the issue of energetic and kinetic bifurcations. The variation of basin boundaries under different sets of controlling parameters provides an analytical way to plot the diversity maps of structural evolution. The project supported by the National Science Foundation (USA) through grant MSS-9258115, and by the National Natural Science Foundation of China     

Infiltration of initially dry, deformable porous media

The present paper studies the infiltration of an incompressible liquid in an initially dry (or partially dry), deformable sponge-like material made of an incompressible constituent in the slug-flow approximation having in mind the application to some industrial processes involving flow through sponge-like materials and, in particular, some composite materials manufacturing processes. The resulting initial-boundary value problem is of Stefan type, with suitable interface conditions and evolution equations describing the position of the interfaces delimiting the saturated region within the porous material. Different models are then suggested in the saturated region, depending on the importance of the inertial terms and on the constitutive equation for the stress. Comparison of the simulation with known experimental results is satisfactory.     

The onset and progression of damage in isotropic paper sheets

An experimental investigation is performed and analyzed in order to examine the onset and evolution of damage processes in thin isotropic paper sheets made of mechanical pulp. A microscopy technique has been used to estimate the relative fraction of bond and fibre breaks. It has been found that the active damage mechanism is bond failure, hence supporting the assumption of an isotropic scalar valued damage variable.All experiments have been performed by simultaneous with the mechanical loading monitoring the acoustic emission activity. Three different experimental setups have been designed offering the possibility to analyze the influence of stress gradients, as well as different levels of the ratios between the in-plane normal stresses, on the onset of damage. It is concluded that stress gradients in the paper specimens have a large influence on the onset of damage. When stress gradients are present a non-local theory has to be used in the analysis. In this way compliance with an isotropic damage criterion is achieved. The characteristic length, determining the gradient sensitivity, has been found to be of the same order of magnitude as some average fibre length.To study the evolution of the damage processes, wide and short specimens have been loaded in tension resulting in stable damage processes. With the assumptions made regarding the mechanical behavior of the paper material after onset of damage, the damage and the cumulative number of acoustic events curve correlates very well. The experimentally obtained data is used to determine material parameters in a proposed damage evolution law. It is found that the assumed damage evolution law can, for isotropic paper materials with bond rupture as the prevalent failure mechanism, be further simplified as only one specific material dependent damage evolution parameter has to be determined in experiments.     

The fundamental role of nonlocal and local balance laws of material forces in finite elastoplasticity and damage mechanics

In this paper the fundamental role of independent balance laws of material forces acting on dislocations and microdefects is shown. They enable a thermodynamically consistent formulation of dissipative deformation processes of continua with dislocation motion and defect evolution in the material space on meso- and microlevel.The balance laws of material forces together with the classical balance laws of physical forces and couples, first and second laws of thermodynamics for physical and material space and general constitutive equations are the basis to develop a thermodynamically consistent framework of nonlocal finite elastoplasticity and brittle and ductile damage.It is shown that a weakly-nonlocal formulation of the balance laws of material forces leads to gradient theories, where local theories are obtained, if all gradient contributions are assumed to be small. In this case the local balance laws of material forces together with the constitutive equations represent evolution laws of the material forces. In the classical approach of internal variables they are assumed from the outset with the result that there is a large number of different propositions in the literature.The well-known splitting test of a circular cylinder of concrete is simulated numerically, where the process of deformation in the physical space and defect and plastic evolution in the material space is represented.     

A rate-dependent two-dimensional free energy model for ferroelectric single crystals

The one-dimensional free energy model for ferroelectric materials developed by Smith et al. [29–31] is generalized to two dimensions. The two-dimensional free energy potential proposed in this paper consists of four energy wells that correspond to four variants of the material. The wells are separated by four saddle points, representing the barriers for 90°-switching processes, and a local maximum, across which 180°-switching processes take place. The free energy potential is combined with evolution equations for the variant fractions based on the theory of thermally activated processes. The model is compared to recent measurements on BaTiO₃ single crystals by Burcsu et al. [8], and predicitions are made concerning the response to the application of in-plane multi-axial electric fields at various frequencies and loading directions. The kinetics of the 90°- and 180°-switching processes are discussed in detail.     

Dynamic effects in materials of complex structure

We suggest a two-component model of a material experiencing structural transformations and use it to study how dynamic and evolution processes affect these transformations. We note that an essential role is played by internal interactions caused by internal forces arising between the components as well as by exchange processes describing variations in the composition of both components. We illustrate the model by specific examples.     

Creep flow, diffusion, and electromigration in small scale interconnects

This paper proposes a three-dimensional electromigration model for void evolution in small scale interconnects. Concurrent kinetics of creep flow and surface diffusion as well as the effect of the surrounding material is considered to provide better understanding of the evolution process. The multiple kinetics and energetics are incorporated into a diffusive interface model. A semi-implicit Fourier spectral method and the preconditioned biconjugate-gradient method are proposed for the computations to achieve high efficiency and numerical stability. We systematically studied kinetic processes in diffusion dominated to creep dominated regime. Which process dominates, as revealed by the analysis, is determined by a combination of viscosity, mobility, interconnect thickness, and void radius. Previous studies on electromigration suggest that a circular void subjected to an electron wind force and surface diffusion is always stable against any small shape perturbation. Our simulations show that a shape that is stable in surface diffusion can become unstable in a creep dominated process, which leads to a quite different void morphology. A spherical void can evolve into a bowl shape and further break into smaller voids. It is also shown that the interconnect geometry has an important effect.     

地圈动力学的THMCB多元关联过程及其分析

地圈中广泛发育的动力学过程属于多元关联过程,它们是在岩土圈与水圈、大气圈、生物圈相互作用中进行的。它们的活动可以概括为热力学、流体力学、固体力学、化学和生物学等基本物质运动的作用及系统的演化。由于它们之间有着互动影响和制约,以及由此引发的能量的转换和参数的变异,因而密切相关。本文建议确立THMCB多元关联过程的概念,开展多元关联过程的研究和分析,以期对地圈动力学过程进行系统有效的综合评价和预测。     

A finite strain constitutive model for metal powder compaction using a unique and convex single surface yield function

Constitutive modelling of metal powder compaction processes is a challenge in view of realistic simulations. To this end, the article under consideration has two objectives: the first goal is to present a new unique and convex single surface yield function for pressure dependent materials, which is also applicable to other areas of granular materials such as soils or concrete. The flexibility is shown at various materials. The yield function is based on a log-interpolation of two known simple yield functions. A convexity proof of the new yield function is provided. The second objective is to propose a new rate-independent finite strain plasticity model for metal powder compaction, which is based on the multiplicative decomposition of the deformation gradient into an elastic and a plastic part with evolution equations for internal variables representing the basic behaviour of powder materials under compaction conditions. These variables are used for the evolution of the yield function in order to represent the compressible hardening behaviour of powder materials. On the basis of the constitutive model, the material parameters are identified at experimental data of copper powder.     

Modeling multi-axial deformation of planar anisotropic elasto-plastic materials,part I: Theory

In sheet metal forming processes local material points can experience multi-axial and multi-path loadings. Under such loading conditions, conventional phenomenological material formulations are not capable to predict the deformation behavior within satisfying accuracy. While micro-mechanical models have significantly improved the understanding of the deformation processes under such conditions, these models require large sets of material data to describe the micromechanical evolution and quite enormous computation expenses for industrial applications. To reduce the drawbacks of phenomenological material models under the multi-path loadings a new anisotropic elasto-plastic material formulation is suggested. The model enables the anisotropic yield surface to grow (isotropic hardening), translate (kinematic hardening) and rotate (rotation of the anisotropy axes) with respect to the deformation, while the shape of the yield surface remains essentially unchanged.Essentially, the model is formulated on the basis of an Armstrong–Frederick type kinematic hardening, the plastic spin theory for the reorientation of the symmetry axes of the anisotropic yield function, and additional terms coupling these expressions. The capability of the model is illustrated with multi-path loading simulations in ‘tension-shear’ and ‘reverse-shear’ to assess its performance with ‘cross’ hardening and ‘Bauschinger’ effects.     

爆轰加载下金属锡层裂破碎数值模拟

对爆轰加载下低熔点金属锡的层裂破碎问题开展了数值模拟。在利用实验数据对所采用数值方法和材料模型开展对比验证的基础上,通过对样品内部物理量时间及空间分布演化对比分析,剖析了冲击加-卸载中样品内部应力波与材料相互作用过程。此外,通过对比分析不同厚度锡样品在爆轰加载下的动态行为特征,进一步认识了自由面反射稀疏波、边侧稀疏波和入射稀疏波共同作用下层裂破碎演化机制。结果表明,当样品较薄时,层裂破碎行为由反射稀疏波主导;随着样品厚度的增大,反射稀疏波主导区缩小,入射稀疏波和边侧稀疏波主导区逐渐增大。     

Irreversible deformation with subsequent unloading of a spherical viscoelastoplastic layer

Analytic solutions are obtained for a sequence of one-dimensional quasistatic problems describing viscoelastic deformation processes in the material of a hollow ball and the plastic flow nucleation and evolution processes occurring in the ball as the pressure on the outer boundary increases. The unloading process under slow removal of the loading pressure is considered as well. The stress fields and the elastic and plastic strain fields in the spherical layer material, the law of motion of the elastoplastic boundary, and the residual stress level and distribution are computed. It is assumed that at the stage preceding the plastic flow the material obeys the viscoelastic Voigt model and the loading surface is determined by the von Mises plastic flow condition.     

Overview on the Prediction Models for Sheet Metal Forming Failure:Necking and Ductile Fracture

The formability of the material determines the amount of available deformation before failure and thus is important for the production of various structural components in industries. The workability of materials is commonly evaluated by different forms of failure models during sheet metal forming(SMF) processes. In order to provide a whole picture about the prediction models for SMF failure, necking-related formability and ductile fracture-related formability studies in SMF processes are systematically summarized, the applicability and limitation of each model are highlighted, and the link between forming limit diagram and ductile fracture criterion is pointed out. Conclusions about some critical issues on failure in SMF are made.     

On the coupling of anisotropic damage and plasticity models for ductile materials

In this contribution various aspects of an anisotropic damage model coupled to plasticity are considered. The model is formulated within the thermodynamic framework and implements a strong coupling between plasticity and damage. The constitutive equations for the damaged material are written according to the principle of strain energy equivalence between the virgin material and the damaged material. The damaged material is modeled using the constitutive laws of the effective undamaged material in which the nominal stresses are replaced by the effective stresses. The model considers different interaction mechanisms between damage and plasticity defects in such a way that two-isotropic and two-kinematic hardening evolution equations are derived, one of each for the plasticity and the other for the damage. An additive decomposition of the total strain into elastic and inelastic parts is adopted in this work. The elastic part is further decomposed into two portions, one is due to the elastic distortion of the material grains and the other is due to the crack closure and void contraction. The inelastic part is also decomposed into two portions, one is due to nucleation and propagation of dislocations and the other is due to the lack of crack closure and void contraction. Uniaxial tension tests with unloadings have been used to investigate the damage growth in high strength steel. A good agreement between the experimental results and the model is obtained.     

A brick model for asperity sintering and creep of APS TBCs

A micromechanical model is developed for the microstructural evolution of an air plasma sprayed (APS), thermal barrier coating: discrete, brick-like splats progressively sinter together at contacting asperities and also undergo Coble creep within each splat. The main microstructural features are captured: the shape, orientation and distribution of asperities between disc-shaped splats, and the presence of columnar grains within each splat. Elasticity is accounted for at the asperity contacts and within each splat, and the high contact compliance explains the fact that APS coatings have a much lower modulus (and thermal conductivity) than that of the parent, fully dense solid. The macroscopic elastic, sintering and creep responses are taken to be transversely isotropic, and remain so with microstructural evolution. Despite the large number of geometric and kinetic parameters, the main features of the behaviour are captured by a small number of characteristic material timescales: these reveal the competition between the deformation mechanisms and identify the rate controlling processes for both free and constrained sintering. The evolution of macroscopic strain, moduli and asperity size is compared for free and constrained sintering, and the level of in-plane stress within a constrained coating is predicted.     

Numerical study of strength properties for a composite material with short reinforcing fibers

A numerical approach to the determination of strength properties for concrete with short reinforcing fibers on the basis of the finite element method is proposed. The mathematical model takes into account various processes of nonlinear deformation of the concrete matrix under compressive and tensile loads, the possibility of developing the inelastic strains in the concrete matrix and reinforcing fibers and the nonlinear interaction between them. The effect of fiber concentration, various loading surfaces for the material matrix, and the bonding type on the deformation of a composite material is analyzed. Numerical examples of strength analysis are given.     

Statistical aspects of the identification of material parameters for elasto-plastic models

Summary The presented method to identify material parameters for inelastic deformation laws is based on the numerical analysis of inhomogeneous stress and strain fields received from suitable experiments. Tensile and bending tests were carried out to obtain elastic and hardening parameters. The deformation law for small elasto-plastic strains is presented as a system of nonlinear differential and algebraic equations (DAE) consisting of the stress–strain relation, evolution equations for the internal variables and the yield condition. Different rules for the evolution equations of isotropic, kinematic and distorsional hardening are proposed. The DAE are discretized using an implicit Euler method, and the resulting system of nonlinear algebraic equations is solved using the Newton method. Deterministic optimization procedures are preferred to identify material parameters from a least-squares functional of numerical and measured comparative quantities. The gradient of the objective function was calculated using a semianalytical sensitivity analysis. Due to measurement errors, the optimal sets of material parameters are non unique. The approximate estimation of confidence regions and the calculation of correlation coefficients is presented. The results of several optimization processes for material parameters of elasto-plastic deformation laws show a good agreement between measured and calculated values, but they show also problems which may occur if systematic errors will not be recognized and deleted. Received 30 September 1999; accepted for publication 8 March 2000