A constitutive model for plastically anisotropic solids with non-spherical voids

Plastic constitutive relations are derived for a class of anisotropic porous materials consisting of coaxial spheroidal voids, arbitrarily oriented relative to the embedding orthotropic matrix. The derivations are based on nonlinear homogenization, limit analysis and micromechanics. A variational principle is formulated for the yield criterion of the effective medium and specialized to a spheroidal representative volume element containing a confocal spheroidal void and subjected to uniform boundary deformation. To obtain closed form equations for the effective yield locus, approximations are introduced in the limit-analysis based on a restricted set of admissible microscopic velocity fields. Evolution laws are also derived for the microstructure, defined in terms of void volume fraction, aspect ratio and orientation, using material incompressibility and Eshelby-like concentration tensors. The new yield criterion is an extension of the well known isotropic Gurson model. It also extends previous analyses of uncoupled effects of void shape and material anisotropy on the effective plastic behavior of solids containing voids. Preliminary comparisons with finite element calculations of voided cells show that the model captures non-trivial effects of anisotropy heretofore not picked up by void growth models.     

A plasticity model for pressure-dependent anisotropic cellular solids

The initial and subsequent yield surfaces for an anisotropic and pressure-dependent 2D stochastic cellular material, which represents solid foams, are investigated under biaxial loading using finite element analysis. Scalar measures of stress and strain, namely characteristic stress and characteristic strain, are used to describe the constitutive response of cellular material along various stress paths. The coupling between loading path and strain hardening is then investigated in characteristic stress–strain domain. The nature of the flow rule that best describes the plastic flow of cellular solid is also investigated. An incremental plasticity framework is proposed to describe the pressure-dependent plastic flow of 2D stochastic cellular solids. The proposed plasticity framework adopts the anisotropic and pressure-dependent yield function recently introduced by Alkhader and Vural [Alkhader M., Vural M., 2009a. An energy-based anisotropic yield criterion for cellular solids and validation by biaxial FE simulations. J. Mech. Phys. Solids 57(5), 871–890]. It has been shown that the proposed yield function can be simply calibrated using elastic constants and flow stresses under uniaixal loading. Comparison of stress fields predicted by continuum plasticity model to the ones obtained from FE analysis shows good agreement for the range of loading paths and strains investigated.     

A unified criterion for the growth and coalescence of microvoids

A yield criterion is developed which unifies void growth and void coalescence theories. Standard void growth theory assumes that plastic flow is diffuse, if not prevalent everywhere within the matrix of the elementary cell considered. On the other hand, void coalescence theory assumes states of post-localized plasticity whereby plastic flow is restricted to intervoid ligaments. The new theory accommodates both scenarios through some appropriate choice of microscopic velocity fields. An important implication for actual evolution problems is a seamless transition from void growth to void coalescence. This is in contrast with previous hybrid approaches whereby abrupt transitions are associated with the presence of unavoidable corners in the effective yield surface. More generally, the new criterion is applicable to describe yielding in porous metal plasticity for both low and high void volume fractions.     

Impact of material processing and deformation on cell morphology and mechanical behavior of polyurethane and nickel foams

The cell morphology and mechanical behavior of open-cell polyurethane and nickel foams are investigated by means of combined 3D X-ray micro-tomography and large scale finite element simulations. Our quantitative 3D image analysis and finite element simulations demonstrate that the strongly anisotropic tensile behavior of nickel foams is due to the cell anisotropy induced by the deformation of PU precursor during the electroplating and heat treatment stages of nickel foam processing. In situ tensile tests on PU foams reveal that the initial main elongation axis of the cells evolves from the foam sheet normal direction to the rolling direction of the coils. Finite element simulations of the hyperelastic behavior of PU foams based on real cell morphology confirm the observation that cell struts do not experience significant elongation after 0.15 tensile straining, thus pointing out alternative deformation mechanisms like complex strut junctions deformation. The plastic behavior and the anisotropy of nickel foams are then satisfactorily retrieved from finite element simulations on a volume element containing eight cells with a detailed mesh of all the hollow struts and junctions. The experimental and computational strategy is considered as a first step toward optimization of process parameters to tailor anisotropy of cell shape and mechanical behavior for applications in batteries or Diesel particulate filtering.     

泡沫金属压痕试验的数值模拟及其反演

在理论研究的基础上，将泡 沫金属压痕试验的有限元数值模拟结果与用无量纲分析法构造出的一系列无量纲函数相结 合，建立了泡沫金属压痕试验中载荷-压痕深度关系曲线与泡沫金属的弹塑性材料参数之 间的联系. 利用这种联系，就可以实现用压痕试验通过反演分析来确定泡沫金属的材料参数. 研究结果表明，泡沫金属材料的杨氏模量，屈服强度及塑性可压缩因子等参数均可由其压痕 试验唯一的确定，但其塑性平台区终点应变的确定还需进一步的研究.     

Characterization of yielding behavior of polycrystalline metals with single crystal plasticity based on representative characteristic length

Single crystal plasticity based on a representative characteristic length is proposed and introduced into a homogenization approach based on finite element analyses, which are applied to characterization of distinctive yielding behaviors of polycrystalline metals, yield-point elongation, and grain size strengthening. The computational manner for an implicit stress update is derived with the framework of a standard multi-surface plasticity at finite strain, where the evolution of the characteristic lengths are numerically converted from the accumulated slips of all of slip systems by exploiting the mathematical feature of the characteristic length as the intermediate function of the plastic internal variables. Furthermore, a constitutive model for a single crystal reproduces the stress–strain curve divided into three parts. Using two-scale finite element analysis, the macroscopic stress–strain response with yield-point elongation under a situation of low dislocation density is reproduced. Finally, the grain size effect on the yield strength is analyzed with modeling of the grain boundary in the context of the proposed constitutive model and is discussed from both macroscopic and microscopic views.     

Mechanism-based multi-surface plasticity model for ideal truss lattice materials

The mechanical behavior of ideal truss lattice materials is controlled by the so-called direct action mechanism at the microscale which involves the uniform stretching and compressing of individual truss members. Standard homogenization techniques have been employed to develop a general micromechanics-based finite-strain constitutive model for truss lattice materials. Furthermore, a specialized small-strain plasticity model has been derived. Both models have been implemented in a finite-element program and used to simulate the anisotropic plastic behavior of the octet-truss lattice material in various applications including cyclic uniaxial loading, pure shear, and three-point bending. The constitutive model predictions agree well with the results obtained from discrete finite element models. Regarding the plasticity of the octet-truss lattice material, it has been found that the elastic domain is constrained by twelve pairwise parallel hyperplanes in the six-dimensional stress space. Moreover, the mechanism-based small-strain formulation reveals that the direction of plastic flow is normal to the pressure-dependent macroscopic yield surfaces.     

A stochastic model for the temporal aspects of flow intermittency in micropillar compression

Systematic experimental investigations have demonstrated that the plastic deformation of micropillar proceeds through a sequence of intermittent bursts, the sizes of which follow power-law statistics. In this study, a stochastic model based on the power-law distribution of burst size is formulated in the framework of crystal plasticity in order to investigate the temporal aspects of flow intermittency in micropillar compression. A Monte Carlo simulation scheme is developed to determine the burst size when a burst activity is captured. This burst size is considered as the displacement boundary condition of burst deformation. Three-dimensional finite element analysis of the model is performed and its predictions are validated by comparison with results from both micro-compression experiments and simulation tests of bulk crystals using the classic crystal plasticity finite element method (CPFEM). The model provides a reasonable prediction of stress–strain responses both at the macroscopic and microscopic scales. Finally, the capability of this model is shown with applications to the intermittent plastic deformation in micropillar compressions, in particular for their burst time durations and burst velocities. The results from such stochastic finite element analysis are shown to be consistent with earlier experimental findings and results of mean-field theory.     

具有恒定冲击载荷的梯度泡沫金属材料设计

多胞材料可通过大变形大量地吸收冲击能量,引入密度梯度可进一步提高其耐撞性。梯度多胞材料的宏观力学响应对材料密度分布极为敏感,不同类型的细观构型的影响也极为不同。已有的研究工作主要局限在对给定的密度梯度分析其动态响应,较少对耐撞性设计方法进行研究。本文针对梯度闭孔泡沫金属材料,基于非线性塑性冲击波模型发展了耐撞性反向设计方法,以维持冲击物受载恒定为目标,运用级数法获得了简化模型和渐近解。利用变胞元尺寸法构建了连续梯度变化的三维Voronoi细观有限元模型,并利用ABAQUS/Explicit有限元软件对理论设计进行数值验证。结果表明,反向设计理论简化模型的渐近解对于梯度闭孔泡沫金属材料的耐撞性设计是有效的,所提出的耐撞性设计方法在控制冲击吸能过程和冲击物受载方面具有指导意义。     

Plastic spin associated with a non-normality theory of plasticity

A plasticity model using a vertex-type plastic flow rule on a smooth yield surface for an anisotropic solid has been proposed recently. This model is here completed by incorporating the effect of plastic spin. Simple shear with a large shear strain is one of the hardest tests on finite strain anisotropic plasticity models, and it is here shown which plastic spin expression is needed to avoid unrealistic oscillatory behavior of the shear stress under large shear strains. The idea of using non-normality with a smooth yield surface originates from a recent proposal of using an abrupt strain path change to determine the subsequent yield surface shape. For this method both polycrystal plasticity calculations and experiments have shown a vertex-type response on the apparently smooth yield surface.     

Continuum modeling of a porous solid with pressure-sensitive dilatant matrix

The pressure-sensitive plastic response of a material has been studied in terms of the intrinsic sensitivity of its yield stress to pressure and the presence and growth of cavities. This work focuses on the interplay between these two distinctly different mechanisms and the attendant material behavior. To this end, a constitutive model is proposed taking both mechanisms into account. Using Gurson's homogenization, an upper bound model is developed for a voided solid with a plastically dilatant matrix material. This model is built around a three-parameter axisymmetric velocity field for a unit sphere containing a spherical void. The void is also subjected to internal pressure; this can be relevant for polymeric adhesives permeated by moisture that vaporizes at elevated temperatures. The plastic response of the matrix material is described by Drucker–Prager's yield criterion and an associated flow rule. The resulting yield surface and porosity evolution law of the homogenized constitutive model are presented in parametric form. Using the solutions to special cases as building blocks, approximate models with explicit forms are proposed. The parametric form and an approximate explicit form are compared against full-field solutions obtained from finite element analysis. They are also studied for loading under generalized tension conditions. These computational simulations shed light on the interplay between the two mechanisms and its enhanced effect on yield strength and plastic flow. Among other things, the tensile yield strength of the porous solid is greatly reduced by the internal void pressure, particularly when a liquid/vapor phase is the source of the internal pressure.     

考虑材料循环塑性的疲劳裂纹扩展模拟

提出了一种考虑材料循环塑性性能的研究疲劳裂纹扩展与闭合行为的有限元模拟方法．对所选用的循环塑性本构关系进行了基本实验检验．探讨了在疲劳裂纹扩展有限元分析中网格尺寸的影响，给出了网格优化准则．研究了在循环硬化条件下考虑裂纹合效应时裂纹面张开廓形、裂纹尖端应力、应变场和正反向塑性区的演变规律．对于循环硬化和不同循环应力比Ｒ等因素对裂纹张开应力水平的影响也作了考察     

Determination of the effective elastic–plastic response of metal–ceramic composites

The mechanical response of metal–ceramic composites is analysed through a homogenization model accounting for the mechanical behaviour of the constituent materials. In order to achieve this purpose a nonlinear homogenization method based on the phase field approach has been suitably implemented into a numerical code. A prescribed homogenized strain state is applied to a unit volume element of a metal–ceramic composite with proportional loading in which all components of the strain tensor are proportional to one scalar parameter. The mechanical response of the material has been modeled by considering a von Mises plasticity model for the metal phase and a Drucker–Prager associative elastic–plastic material model for the ceramic phase. A two stages plasticity has been obtained in which inelastic strain develops in the metal phase followed by a fully plastic response. A comparison with a finite element model of the stress–strain response of an axisymmetric unit cell has been carried out with the purpose to validate the homogenization based modeling presented in the paper. Plastic parameters of a Drucker–Prager yield surface for the homogenized composite have been calculated at different materials compositions. Associative Drucker–Prager plasticity has been found to be accurate for high ceramic content.     

Three-dimensional modeling of the mechanical property of linearly elastic open cell foams

Three-dimensional Voronoi models are developed to investigate the mechanical behavior of linearly elastic open cell foams. Dependence of the Young’s modulus, Poisson’s ratio and bulk modulus of the foams on the relative density is evaluated through finite element analysis. Obtained results show that in the low density regime the Young’s modulus and bulk modulus of random Voronoi foams can be well represented by those of Kelvin foams, and are sensitive to the geometric imperfections inherent in the microstructure of foams. In contrast, the compressive plateau stress of the foams is less sensitive to the imperfections. Failure surface of the foams subject to multi-axial compression is determined and is found to comply with the maximum compressive principal stress criterion, consistent with available experimental observations on polymer foams. Numerical results also show that elastic buckling of cell edges at microscopic level is the dominant mechanism responsible for the compressive failure of elastic open cell foams.     

空心颗粒填充复合材料的建模与弹塑性仿真计算

空心玻璃微珠（hollow glass microballon,简称HGM）填充复合材料称为复合泡沫材料,近来已被广泛应用在工业中。利用RSA(Random Sequential Adsorption)方法生成了含不同体积比的HGMs填充代表体元模型,然后用有限元方法计算得出了材料的应力-应变关系,将材料属性简化为双线性随动强化模型,分析了HGM填充比、壁厚对材料的有效弹性常数、屈服极限的影响,并分析了材料内部细观应力场和塑性应变分布情况。结果发现,HGM壁厚比对材料有效弹性模量和屈服极限的增减起着决定性作用,而对于任意填充模式来说,其比模量和比强度总是大于纯树脂,这一点体现了该材料轻质的优良特性。材料基体的应力集中部位分布以及塑性应变区域的分布也取决于HGM的壁厚比。     

STOCHASTIC ELASTO-PLASTIC FRACTURE ANALYSIS OF ALUMINUM FOAMS

Model I quasi-static nonlinear fracture of aluminum foams is analyzed by considering the effect of microscopic heterogeneity. Firstly, a continuum constitutive model is adopted to account for the plastic compressibility of the metallic foams. The yield strain modeled by a two- parameter Weibull-type function is adopted in the constitutive model. Then, a modified cohesive zone model is established to characterize the fracture behavior of aluminum foams with a cohesive zone ahead of the initial crack. The tensile traction versus local crack opening displacement relation is employed to describe the softening characteristics of the material. And a Weibull statistical model for peak bridging stress within the fracture process zone is used for considering microscopic heterogeneity of aluminum foams. Lastly, the influence of stochastic parameters on the curve of stress-strain is given. Numerical examples are given to illustrate the numerical model presented in this paper and the effects of Weibull parameters and material properties on J-integral are discussed.     

A yield function for single crystals containing voids

A yield function for single crystals containing voids has been developed based on a variational approach. This first yield function is phenomenologically extended by modifying the dependence on the mean stress and introducing three adjustable parameters. Unit cell finite element calculations are performed for various stress triaxiality ratios, main loading directions and porosity levels in the case of a perfectly plastic FCC single crystal. The three model parameters are adjusted on the unit cell calculations so that a very good agreement between simulation results and the proposed model is obtained.     

Nonlocal versus local elastoplastic behaviour of heterogeneous materials

The scale transition methods have been developed for many years in order to obtain the overall behavior of polycrystalline materials from their microscopic behavior and their microstructure. Nevertheless, some basic aspects are absent from such formalisms. The most significant one seems to be the heterogeneization by plastic straining which involves nonlocality of hardening. In this article, a nonlocal theory based upon crystalline plasticity is developed from which a nonlocal constitutive equation at the grain level is derived. With regard to the polycrystal, in order to deduce the behavior of a local equivalent homogeneous medium, an integral equation is proposed and solved for nonlocal inhomogeneous materials by the self-consistent approximation. This scheme is developed in case of a two-phase nonlocal material representing the dislocation cell structure induced during plastic straining. Numerical simulations based on a simplified model show significant effects on the intragranular heterogeneization.     

A novel yield function for single crystals based on combined constraints optimization

A novel yield function representing the overall plastic deformation in a single crystal is developed using the concept of optimization. Based on the principle of maximum dissipation during a plastic deformation, the problem of single crystal plasticity is first considered as a constrained optimization problem in which constraints are yield functions for slip systems. To overcome the singularity that usually arises in solving the above problem, a mathematical technique is used to replace the above constrained optimization problem with an equivalent problem which has only one constraint. This single constraint optimization problem, the so-called combined constraints crystal plasticity (CCCP) model, is implemented into a finite element code and the results of modeling the uniaxial tensions of the single crystal copper along different crystallographic directions and also hydroforming of aluminum tubes proved the capability of the proposed CCCP model in accurately predicting the deformation in polycrystalline materials.