On the evolution of lattice deformation in austenitic stainless steels—The role of work hardening at finite strains

In this work, a three dimensional crystal plasticity-based finite element model is presented to examine the micromechanical behaviour of austenitic stainless steels. The model accounts for realistic polycrystal micromorphology, the kinematics of crystallographic slip, lattice rotation, slip interaction (latent hardening) and geometric distortion at finite deformation. We utilise the model to predict the microscopic lattice strain evolution of austenitic stainless steels during uniaxial tension at ambient temperature with validation through in situ neutron diffraction measurements. Overall, the predicted lattice strains are in very good agreement with those measured in both longitudinal and transverse directions (parallel and perpendicular to the tensile loading axis, respectively). The information provided by the model suggests that the observed nonlinear response in the transverse {200} grain family is associated with a competitive bimodal evolution of strain during inelastic deformation. The results associated with latent hardening effects at the microscale also indicate that in situ neutron diffraction measurements in conjunction with macroscopic uniaxial tensile data may be used to calibrate crystal plasticity models for the prediction of the inelastic material deformation response.     

Study of lattice strains in magnesium alloy AZ31 based on a large strain elastic-viscoplastic self-consistent polycrystal model

The recently developed large strain elastic visco-plastic self-consistent (EVPSC) model, which incorporates both slip and twinning deformation mechanisms, is used to study the lattice strain evolution in extruded magnesium alloy AZ31 under uniaxial tension and compression. The results are compared against in-situ neutron diffraction measurements done on the same alloy. For the first time, the effects of stress relaxation and strain creep on lattice strain measurements in respectively displacement controlled and load controlled in-situ tests are numerically assessed. It is found that the stress relaxation has a significant effect on the lattice strain measurements. It is also observed that although the creep does not significantly affect the trend of the lattice strain evolution, a better agreement with the experiments is found if creep is included in the simulations.     

Mechanical Testing of Thin Sheet Magnesium Alloys in Biaxial Tension and Uniaxial Compression

Tension and compression experiments on magnesium rolled sheets and extruded products of AZ31 (Mg?+?3%Al?+?1%Zn) and ZE10 (Mg?+?1%Zn?+?0.3%Ce based misc metal) were performed at room temperature. The tests were conducted along the longitudinal and the transverse direction to quantify the in-plane anisotropy. Samples built from adhesively-bonded layers of sheets were used for in-plane as well as through-thickness compression testing. It was verified that this simple testing method leads to identical results as using comb-like dies and equi-biaxial bulge testing, respectively. In the case of uniaxial loading, the longitudinal and transverse strain components were measured using independent extensometers. R-values were calculated from these signals. The mechanical responses were correlated to the microstructure and the texture. The recorded differences between tensile and compressive response reveal the strength differential effect of the materials. The distortional character of the plastic behaviour is evidenced through their responses to equi-biaxial tensile loading. Significant differences in the compressive responses of the two alloys were identified by comparing the respective hardening rates.     

Modeling the effects of material non-linearity using moving window micromechanics

In the analysis of materials with random heterogeneous microstructure the assumption is often made that material behavior can be represented by homogenized or effective properties. While this assumption yields accurate results for the bulk behavior of composite materials, it ignores the effects of the random microstructure. The spatial variations in these microstructures can focus, initiate and propagate localized non-linear behavior, subsequent damage and failure. In previous work a computational method, moving window micromechanics (MW), was used to capture microstructural detail and characterize the variability of the local and global elastic response. Digital images of material microstructure described the microstructure and a local micromechanical analysis was used to generate spatially varying material property fields. The strengths of this approach are that the material property fields can be consistently developed from digital images of real microstructures, they are easy to import into finite element models (FE) using regular grids, and their statistical characterizations can provide the basis for simulations further characterizing stochastic response. In this work, the moving window micromechanics technique was used to generate material property fields characterizing the non-linear behavior of random materials under plastic yielding; specifically yield stress and hardening slope, post yield. The complete set of material property fields were input into FE models of uniaxial loading. Global stress strain curves from the FE–MW model were compared to a more traditional micromechanics model, the generalized method of cells. Local plastic strain and local stress fields were produced which correlate well to the microstructure. The FE–MW method qualitatively captures the inelastic behavior, based on a non-linear flow rule, of the sample continuous fiber composites in transverse uniaxial loading.     

Incorporation of twinning into a crystal plasticity finite element model: Evolution of lattice strains and texture in Zircaloy-2

A crystal plasticity finite element code is developed to model lattice strains and texture evolution of HCP crystals. The code is implemented to model elastic and plastic deformation considering slip and twinning based plastic deformation. The model accounts for twinning reorientation and growth. Twinning, as well as slip, is considered to follow a rate dependent formulation. The results of the simulations are compared to previously published in situ neutron diffraction data. Experimental results of the evolution of the texture and lattice strains under uniaxial tension/compression loading along the rolling, transverse, and normal direction of a piece of rolled Zircaloy-2 are compared with model predictions. The rate dependent formulation introduced is capable of correctly capturing the influence of slip and twinning deformation on lattice strains as well as texture evolution.     

On the transverse compression response of Kevlar KM2 using fiber-level finite element model

Flexible textile composites like woven Kevlar fabrics are widely used in high velocity impact (HVI) applications. Upon HVI they are subjected to both longitudinal tensile and transverse compressive loads. To understand the role of transverse properties, the single fiber and tow transverse compression response (SFTCR and TTCR) of Kevlar KM2 fibers are numerically analyzed using plane strain finite element (FE) models. A finite strain formulation with a minimum number of 84 finite elements is determined to be required for the fiber cross section to capture the finite strain SFTCR through a mesh convergence study. Comparison of converged numerical solution to the experimental results indicates the dominant role of geometric stiffening at finite strains due to growth in contact width. The TTCR is studied using a fiber length scale FE model of a single tow comprised of 400 fibers transversely loaded between rigid platens. This study along with micrographs of yarn after mechanical compaction illustrates fiber spreading and fiber–fiber contact friction interactions are important deformation mechanisms at finite strains. The TTCR is also studied using homogenized yarn level models with properties from the literature. Comparison of TTCR between fiber length scale and homogenized yarn length scale models indicate the need for a nonlinear material model for homogenized approaches to accurately predict the transverse compression response of the fabrics.     

Compressive response and failure of braided textile composites: Part 2—computations

In this, the second part of a two part paper, results obtained by using the finite element (FE) method in conjunction with micromechanics to predict the effective elastic stiffness and strength of a carbon 2D triaxially braided composite (2DTBC), are presented. The 3D FE based micromechanics study was carried out on one representative unit cell (RUC) of the carbon 2DTBC (the “micromodel”). The FE models were first used to determine the macroscopic elastic orthotropic stiffnesses of the 2DTBC. The micromodel was deemed acceptable (in terms of the number of elements used in the mesh of the micromodel) if the elastic stiffnesses it displayed were within 5% of the elastic properties found experimentally. Subsequently, buckling eigenmodes were determined for the FE RUC under uniaxial and biaxial loading states, corresponding to the experimental investigation reported in part I of this two part paper. The lowest symmetric modes were identified and these mode shapes were used as imperfections to the FE model for a subsequent nonlinear response analysis using an arc-length method in conjunction with the ABAQUS commercial FE code. The magnitude of the imperfections was left as a parameter and its effect on the predicted response was quantified. The present micromechanics computational model provides a means to assess the compressive and compressive/tensile biaxial strength of the braided composites and its dependence on various microstructural parameters. It also serves as a tool to assess the most significant parameter that affects compressive strength.     

Dynamic tensile characterization of pig skin

The strain-rate dependent response of porcine skin oriented in the fiber direction is explored under tensile loading. Quasi-static response was obtained at strain rates in the range of 10-3s-1to 25 s-1. Characterization of the response at even greater strain rates is accomplished by measuring the spatio-temporal evolution of the particle velocity and strain in a thin strip subjected to high speed impact loading that generates uniaxial stress conditions. These experiments indicate the formation of shock waves; the shock Hugoniot that relates particle velocity to the shock velocity and the dynamic stress to dynamic strain is obtained directly through experimental measurements, without any assumptions regarding the constitutive properties of the material.     

半结晶聚合物损伤演化的实验表征与数值模拟

损伤本构模型对研究材料的断裂失效行为有重要意义, 但聚合物材料损伤演化的定量表征实验研究相对匮乏. 通过4种高密度聚乙烯(high density polythylene, HDPE)缺口圆棒试样的单轴拉伸实验获得了各类试样的载荷-位移曲线和真应力-应变曲线, 采用实验和有限元模拟相结合的方法确定了HDPE材料不同应力状态下的本构关系, 并建立了缺口半径与应力三轴度之间的关系;采用两阶段实验法定量描述了4种HDPE试样单轴拉伸过程中的弹性模量变化, 并建立了基于弹性模量衰减的损伤演化方程, 结合中断实验和扫描电子显微镜分析了应力状态对HDPE材料微观结构演化的影响. 结果表明缺口半径越小, 应力三轴度越大, 损伤起始越早、演化越快; 微观表现为: 高应力三轴度促进孔洞的萌生和发展, 但抑制纤维状结构的产生;基于实验和有限元模拟获得的断裂应变、应力三轴度、损伤演化方程等信息提出了一种适用于聚合物的损伤模型参数确定方法, 最后将本文获得的本构关系和损伤模型用于HDPE平板的冲压成形模拟, 模拟结果与实验结果吻合良好.      

Modeling lattice strain evolution at finite strains and experimental verification for copper and stainless steel using in situ neutron diffraction

An elasto-plastic self-consistent (EPSC) polycrystal model is extended to account, in an approximate fashion, for the kinematics of large strains, rigid body rotations, texture evolution and grain shape evolution. In situ neutron diffraction measurements of the flow stress, internal strain, texture and diffraction peak intensity evolutions were performed on polycrystalline copper and stainless steel, up to true tensile strains of ε = 0.3. Suitably adjusted slip system hardening model parameters enable the model to quantitatively describe the flow stress of the polycrystalline aggregate. Quantitative predictions of the texture evolution and the internal strain evolution along the stress axis are good, while predictions of transverse internal strains (perpendicular to the tensile loading direction) are less satisfactory. The latter exhibit a large dispersion from grain to grain around a macroscopic average, and the implications of this finding for the interpretation of in situ neutron diffraction method are explored. Finally, as a demonstration of the applicability of the model to problems involving finite rotation, as well as deformation, simulations of simple shear were conducted which predict a texture evolution in agreement with published experimental data, and other modeling approaches as well.     

Mechanical response and texture evolution of AZ31 alloy at large strains for different strain rates and temperatures

In order to study the behavior of material under finite deformation at various strain rates, the responses of AZ31 Mg sheet are measured under uniaxial (tension and compression) and multiaxial (simple shear) loadings along rolling direction (RD), 45° to rolling direction (DD), 90° to rolling direction (TD), and normal to the sheet (ND) to large strains. The material exhibits positive strain rate sensitivity (SRS) at room and elevated temperatures; the SRS is more pronounced at high temperatures and lower strain rates. The r-value of the material under tensile loading at room temperatures is higher in TD at lower strain rate. Texture measurements on several failed specimens are reported under tension and simple shear after finite plastic deformation of about 20% equivalent strain. The as-received material exhibits a strong fiber with equal fractions of grains having the c-axis slightly tilted away from the sheet normal towards both +RD and −RD. Pole figures obtained after tensile loading along the rolling direction (RD) show that the texture of the material strengthens even at low strains, with c-axis perpendicular to the sheet plane and prism planes lining up in a majority of grains. However, the tensile loading axis along TD does not lead to similar texture strengthening; the c-axis distribution appears to be virtually unchanged from the virgin state. The pole figures obtained after in-plane compression along RD brings the c-axes of the grains parallel to the loading direction. The pole figures after simple shear loading show that the c-axis rotates to lie on the sheet plane consistent with a compression axis 45° away on the sheet plane.     

Micromechanical modeling of low cycle fatigue under complex loadings — Part II. Applications

Motivated by the effect of the different loading paths (uniaxial and multiaxial) and the mean stress or strain effects, fatigue life is qualitatively evaluated by employing a micromechanical approach developed in the first part of this paper. Taking into account these loading conditions, a number of numerical tests are carried out to demonstrate the ability of this model to describe the inelastic-damaged behavior of low cyclic fatigue. Theoretically, the internal variables evolution on the Crystallographic Slip Systems (CSS) are recorded under these loadings up to the final fracture of the Representative Volume Element (RVE). The predicted macro-damage D^T as well as the applied macro-stress Σ are recorded and analyzed during fatigue life. It is demonstrated that, at least qualitatively, the developed model gives successful results describing the fatigue lives under employed loading paths. Quantitatively, after the identification of the material constants of the over-aged Waspaloy, two examples are taken into consideration representing the uniaxial and biaxial cyclic behavior of this alloy. It is shown that the predicted responses are in good agreement with the experimental results.     

Anisotropic finite-strain models for porous viscoplastic materials with microstructure evolution

In this work, we propose a constitutive model for the finite-strain, macroscopic response of porous viscoplastic solids, accounting for deformation-induced changes in the size, shape and distribution of the voids. The model makes use of consistent homogenization estimates obtained by the “iterated variational linear comparison” procedure of Agoras and Ponte Castañeda (2013) to characterize both the instantaneous effective response of the porous material and the evolution of the underlying microstructure. The proposed model applies for general, three-dimensional loading conditions and can be implemented numerically for use in standard FEM codes. We also investigate the interplay between the evolution of the microstructure and the macroscopic stress–strain response, in the context of displacement-controlled, plane strain loading (bi-axial straining) of initially isotropic, porous, rigid-plastic materials with power-law hardening. We focus on the effect of strain triaxiality, and consider both extensional and contractile loading conditions leading to porosity growth and collapse, respectively. For both types of loadings, it is found that the macroscopic behavior of the material exhibits an initial hardening regime followed by a softening regime at sufficiently large strains. Consistent with earlier models and experimental results, the softening regime for extensional loadings is a consequence of porosity growth. On the other hand, the softening behavior predicted for contractile loadings is found to be a consequence of void collapse, i.e., of rapid changes in the average shape of the pores leading to crack-like shapes. For both types of loadings, the transition from hardening to softening in the macroscopic response can be identified with the onset of macroscopic strain localization. The associated critical conditions at the onset of localization are determined as a function of the strain triaxiality. The type of localization band ranges from dilatational to compaction bands as the bi-axial straining varies from uniaxial extension to uniaxial contraction.     

Size-effects on yield surfaces for micro reinforced composites

Size effects in heterogeneous materials are studied using a rate independent higher order strain gradient plasticity theory, where strain gradient effects are incorporated in the free energy of the material. Numerical studies are carried out using a finite element method, where the components of the plastic strain tensor appear as free variables in addition to the displacement variables. Non-conventional boundary conditions are applied at material interfaces to model a constraint on plastic flow due to dislocation blocking. Unit cell calculations are carried out under generalized plane strain conditions. The homogenized response of a material with cylindrical reinforcing fibers is analyzed for different values of the internal material length scale and homogenized yield surfaces are presented. While the main focus is on initial yield surfaces, subsequent yield surfaces are also presented. The center of the yield surface is tracked under uniaxial loading both in the transverse and longitudinal directions and an anisotropic Bauschinger effect is shown to depend on the size of the fibers. Results are compared to conventional predictions, and size-effects on the kinematic hardening are accentuated.     

Coupled experiments and simulations of microstructural damage in wood

In this paper, we explore ways to couple experimental measurements with the numerical simulations of the mechanical properties of wood. For our numerical simulations, we have adopted a lattice approach, where wood fibers or bundles of wood fibers are modeled as discrete structural elements connected by a lattice of spring elements. Element strength and stiffness properties are determined from bulk material properties. Damage is represented by broken lattice elements, which cause both stiffness and strength degradation. The modeling approach was applied to small specimens of spruce subjected to transverse uniaxial tension, and mode I transverse splitting. The model was found to be good at predicting the load-deformation response of both notched and unnotched specimens, including the post-peak softening response. In addition, the damage patterns predicted by the model are consistent with those observed in the experiments.     

An analytical delamination model for laminated plates including bridging effects

A penalised interface model, whose strain energy is the penalty functional related to interface adhesion constraint, is introduced in conjunction with a damageable interface whose local constitutive law, in turn, represents bridging stress effects, in order to analyse delamination and bridging phenomena in laminated plates. The laminate is modelled by means of first-order shear deformable layer-wise kinematics and the governing equations are formulated in the form of a non-linear differential system with moving intermediate boundary conditions related to opportune delamination and bridging growth conditions. The problem is solved through an analytical approach. The model leads to an accurate and self-consistent evaluation of the energy release rate and its mode components due to the inclusion of significant contributions arising from coupling between in-plane and transverse shear stresses, and to an asymptotic estimate of interlaminar stresses. The salient features of the proposed model are investigated in the context of an energy balance approach and of a J-integral formulation, thus providing simple results useful to model delamination growth and bridging behaviour when mixed mode loading is involved. The accuracy of the proposed model is substantiated through comparisons with results from continuum analysis obtained by a finite element (FE) procedure. The effectiveness of the proposed model is highlighted by showing the solution of a two-layered plate scheme subjected to pure and mixed mode loading conditions and to fibre bridging stresses. The results point out that the present model, despite its low computational cost in comparison with more complex FE analyses, is an efficient tool to predict delamination and bridging evolution.     

Stochastic damage model for concrete based on energy equivalent strain

Starting with the framework of conventional elastoplastic damage mechanics, a class of stochastic damage constitutive model is derived based on the concept of energy equivalent strain. The stochastic damage model derived from the parallel element model is adopted to develop the uniaxial damage evolution function. Based on the expressions of damage energy release rates (DERRs) conjugated to the damage variables thermodynamically, the concept and its tensor formulations of energy equivalent strain is proposed to bridge the gap between the uniaxial and the multiaxial constitutive models. Furthermore, a simplified coupling model is proposed to consider the evolution of plastic strain. And the analytical expressions of the constitutive model in 2-D are established from the abstract tensor expression. Several numerical simulations are presented against the biaxial loading test results of concrete, demonstrating that the proposed models can reflect the salient features for concrete under uniaxial and biaxial loading conditions.     

An evaluation of a combined isotropic-kinematic hardening model for representation of complex strain-path changes in dual-phase steel

This paper aims at evaluating an elastoplastic constitutive model which accounts for combined isotropic-kinematic hardening for complex strain-path changes in a dual-phase steel, DP800. The capability of the model to reproduce the transient hardening phenomena under two-stage non-proportional loading has been assessed through numerical simulations of sequential uniaxial tension and notched tension/shear tests. Finite element simulations with shell elements were performed using the explicit non-linear FE code LS-DYNA. Numerical predictions of the stress–strain response were compared to the corresponding experimental data. The results from the experiments demonstrated that prior plastic deformation has certainly influenced the subsequent work-hardening behaviour of the material under biaxial or shear deformation modes. Furthermore, the numerical simulations from the two-stage uniaxial tension–notched tension and uniaxial tension–shear tests predicted the general trends of the experimental results such as transitory hardening and overall work hardening. However, some discrepancies were found in accurately describing the transient hardening behaviour subsequent to strain path changes between the experiments and numerical simulations.     

Thermomechanics of shape memory polymers: Uniaxial experiments and constitutive modeling

Shape memory polymers (SMPs) can retain a temporary shape after pre-deformation at an elevated temperature and subsequent cooling to a lower temperature. When reheated, the original shape can be recovered. Relatively little work in the literature has addressed the constitutive modeling of the unique thermomechanical coupling in SMPs. Constitutive models are critical for predicting the deformation and recovery of SMPs under a range of different constraints. In this study, the thermomechanics of shape storage and recovery of an epoxy resin is systematically investigated for small strains (within ±10%) in uniaxial tension and uniaxial compression. After initial pre-deformation at a high temperature, the strain is held constant for shape storage while the stress evolution is monitored. Three cases of heated recovery are selected: unconstrained free strain recovery, stress recovery under full constraint at the pre-deformation strain level (no low temperature unloading), and stress recovery under full constraint at a strain level fixed at a low temperature (low temperature unloading). The free strain recovery results indicate that the polymer can fully recover the original shape when reheated above its glass transition temperature (T_g). Due to the high stiffness in the glassy state (T < T_g), the evolution of the stress under strain constraint is strongly influenced by thermal expansion of the polymer. The relationship between the final recoverable stress and strain is governed by the stress–strain response of the polymer above T_g. Based on the experimental results and the molecular mechanism of shape memory, a three-dimensional small-strain internal state variable constitutive model is developed. The model quantifies the storage and release of the entropic deformation during thermomechanical processes. The fraction of the material freezing a temporary entropy state is a function of temperature, which can be determined by fitting the free strain recovery response. A free energy function for the model is formulated and thermodynamic consistency is ensured. The model can predict the stress evolution of the uniaxial experimental results. The model captures differences in the tensile and compressive recovery responses caused by thermal expansion. The model is used to explore strain and stress recovery responses under various flexible external constraints that would be encountered in applications of SMPs.     

Higher-order theory for periodic multiphase materials with inelastic phases

An extension of a recently-developed linear thermoelastic theory for multiphase periodic materials is presented which admits inelastic behavior of the constituent phases. The extended theory is capable of accurately estimating both the effective inelastic response of a periodic multiphase composite and the local stress and strain fields in the individual phases. The model is presently limited to materials characterized by constituent phases that are continuous in one direction, but arbitrarily distributed within the repeating unit cell which characterizes the material's periodic microstructure. The model's analytical framework is based on the homogenization technique for periodic media, but the method of solution for the local displacement and stress fields borrows concepts previously employed by the authors in constructing the higher-order theory for functionally graded materials, in contrast with the standard finite-element solution method typically used in conjunction with the homogenization technique. The present approach produces a closed-form macroscopic constitutive equation for a periodic multiphase material valid for both uniaxial and multiaxial loading. The model's predictive accuracy in generating both the effective inelastic stress-strain response and the local stress and inelastic strain fields is demonstrated by comparison with the results of an analytical inelastic solution for the axisymmetric and axial shear response of a unidirectional composite based on the concentric cylinder model and with finite-element results for transverse loading.