确定复合材料宏观屈服准则的细观力学方法

运用细观力学中的均匀化方法，分析了含周期性微结构复合材料的宏观屈服准则，并对Hill-Tsai准则进行了修正。从基于复合材料细观结构的代表性胞元入手，运用塑性极限理论中的机动分析以及有限元方法，计算了细观结构的极限载荷域。通过宏细观尺度对应关系，得到复合材料的宏观屈服准则。     

Plastic flow for non-monotonic loading conditions of an aluminum alloy sheet sample

Non-linear deformation paths obtained using uniaxial tension followed by simple shear tests were performed for a 1050-O aluminum alloy sheet sample in different specimen orientations with respect to the material symmetry axes. In order to eliminate the time influence, the time interval between the first and second loading steps was kept constant for all the tests. Monotonic uniaxial tension tests interrupted during loading were used to assess the recovery that takes place during this time. In order to eliminate the influence of the initial plastic anisotropy and to compare the results as if the material hardening was isotropic, the flow stress was represented as a function of the plastic work. The behavior of the material after reloading was analyzed in terms of dislocation microstructure and crystallographic texture evolutions. For more quantitative assessment, the full constraints [Int. J. Plasticity 13 (1997) 75] and visco-plastic self-consistent [Acta Metall. Mater. 41 (1993) 2611] polycrystal models were used to simulate the material behavior in the non-linear deformation paths. Based on experimental and simulation results, the relative contributions of the crystallographic texture and dislocation microstructure evolution to the anisotropic hardening behavior of the material were discussed.     

金属材料在复杂应力状态下的塑性流动特性及本构模型

工程应用中,金属材料和结构往往处于复杂应力状态。材料的塑性行为会受到应力状态的影响,要精确描述材料在复杂应力状态下的塑性流动行为,必须在本构模型中考虑应力状态效应的影响。然而,由于在动态加载下材料的应变率效应和应力状态效应相互耦合、难以分离,给应力状态效应的研究和模型的建立造成很大困难。通过对Ti-6Al-4V钛合金材料开展不同加载条件下的力学性能测试,提出了一个包含应力三轴度和罗德角参数影响的新型本构模型,并通过VUMAT用户子程序嵌入ABAQUS/Explicit软件。分别采用新提出的塑性模型和Johnson-Cook模型对压剪复合试样的动态实验进行了数值模拟。结果表明,新模型不仅在对材料本构曲线的拟合方面具有较强的优势,而且由该模型所得到的透射脉冲和载荷-位移曲线均更加准确。因此,该模型能够更精确地描述和预测金属材料在复杂应力状态下的塑性流变行为。     

Folding of fiber composites with a hyperelastic matrix

This paper presents an experimental and numerical study of the folding behavior of thin composite materials consisting of carbon fibers embedded in a silicone matrix. The soft matrix allows the fibers to microbuckle without breaking and this acts as a stress relief mechanism during folding, which allows the material to reach very high curvatures. The experiments show a highly non-linear moment vs. curvature relationship, as well as strain softening under cyclic loading. A finite element model has been created to study the micromechanics of the problem. The fibers are modeled as linear-elastic solid elements distributed in a hyperelastic matrix according to a random arrangement based on experimental observations. The simulations obtained from this model capture the detailed micromechanics of the problem and the experimentally observed non-linear response. The proposed model is in good quantitative agreement with the experimental results for the case of lower fiber volume fractions but in the case of higher volume fractions the predicted response is overly stiff.     

Homogenization in probabilistic terms: The variational principle and some approximate solutions

Studying of materials with evolving random microstructures requires the knowledge of probabilistic characteristics of local fields because the path of the microstructure evolution is controlled by the local fields. The probabilistic characteristics of local fields are determined by the probabilistic characteristics of material properties. In this paper it is considered the problem of finding the probabilistic characteristic of local fields, if the probabilistic characteristics of material properties are given. The probabilistic characteristics of local fields are sought from the variational principle for probabilistic measure. Minimizers of this variational problem provide all statistical information of local fields as well as the effective coefficients. Approximate solutions are obtained for electric current in composites for two cases: multi-phase isotropic composites with lognormal distribution of conductivities and two-phase isotropic composites. The solutions contain a lot of statistical information that has not been available previously by analytical treatments.     

Statistical characterization of meso-scale uniaxial compressive strength in brittle materials with randomly occurring flaws

Increasingly fine spatial resolution in numerical models of brittle materials promises to improve prediction and characterization of dynamic failure in these materials. However, as the resolution of these numerical models begins to approach the material micro-scale, the associated discretization requires a definitive connection to the microstructure. In many cases a numerical model (e.g., a finite element mesh) that explicitly resolves each flaw within the material is not feasible for macro-scale analyses. As an alternative, each element can be treated as a meso-scale continuum with constitutive properties that reflect the characteristics of the underlying microstructure. Small scale elements will exhibit random variations in the constitutive properties as a result of the random variations in the number and types of flaws and the flaw sizes contained within each element. The present paper proposes a technique for assigning probability distributions to these element properties, which can be thought of as the meso-scale constitutive properties. In particular, the strain-rate dependent compressive uniaxial strength of a ceramic is modeled using a two-dimensional analytical model developed by Paliwal and Ramesh (2008). The effect on the probability distribution of meso-scale (or element-level) strength from flaw density, flaw size distribution, flaw clustering, and strain rate are studied. Higher strain rates, more flaw clustering, and decreasing element size all contribute to greater scatter in uniaxial compressive strength. Variations in flaw size increase the scatter in the strength more for low strain rate loadings and less clustered microstructures. The results provide interesting comparisons to the classical assumption of a two-parameter Weibull-distributed strength, showing that a three-parameter Weibull distribution and even a lognormal distribution fit better with the simulated strength data.     

A computational micromechanics constitutive model for the unloading behavior of paper

In this study, a computational micromechanics material model for the unloading behavior of paper and other nonwoven materials is presented. The asymptotic fiber and bond (AFB) model for paper elastic–plastic behavior [Sinha, S.K., Perkins, R.W., 1995. Micromechanics constitutive model for use in finite element analysis, In: Proceedings of the 1995, Joint ASME Applied Mechanics and Materials Summer Meeting, Los Angeles, CA, USA, Jun 28–30, 1995] has been extended to model the unloading process through a computational algorithm and implemented using the UMAT subroutine in ABAQUS finite element code. For every unloading increment, the material model assumes elastic unloading with a slope equal to the initial elastic modulus. The Jacobian matrix of the constitutive model is updated at every unloading increment by applying the incremental form of AFB model for a planar element with an elastic fiber and bond condition. A uniaxial tensile and a biaxial Mullen burst loading–unloading experiments were carried out for a paperboard sample and simulated using the model. The stress–strain curve and residual strain for the uniaxial loading were in good agreement with experimental results. The finite element model of the burst test with the AFB unloading material model predicted the general shape of the pressure versus deflection curve. However, the model over predicted the residual deflection by more than 50%. The loading portion of the pressure–deflection curve had a significant offset from experimental curves, and the nonlinearity in the unloading curve towards the end was not predicted. The discrepancies with experimental results are attributed to the burst test itself, model parameter estimation inadequacies, boundary conditions used in the FEA, and neglecting time-dependant effects. Nevertheless, the model can be useful in parametric studies relating microstructure to unloading behavior in structural problems.     

Effect of a nonuniform distribution of voids on the plastic response of voided materials: a computational and statistical analysis

This study investigates the overall and local response of porous media composed of a perfectly plastic matrix weakened by stress-free voids. Attention is focused on the specific role played by porosity fluctuations inside a representative volume element. To this end, numerical simulations using the Fast Fourier Transform (FFT) are performed on different classes of microstructure corresponding to different spatial distributions of voids. Three types of microstructures are investigated: random microstructures with no void clustering, microstructures with a connected cluster of voids and microstructures with disconnected void clusters. These numerical simulations show that the porosity fluctuations can have a strong effect on the overall yield surface of porous materials. Random microstructures without clusters and microstructures with a connected cluster are the hardest and the softest configurations, respectively, whereas microstructures with disconnected clusters lead to intermediate responses. At a more local scale, the salient feature of the fields is the tendency for the strain fields to concentrate in specific bands. Finally, an image analysis tool is proposed for the statistical characterization of the porosity distribution. It relies on the distribution of the ‘distance function’, the width of which increases when clusters are present. An additional connectedness analysis allows us to discriminate between clustered microstructures.     

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.     

A micromechanics approach to homogenizing elasto-viscoplastic heterogeneous materials

The variational asymptotic method for unit cell homogenization (VAMUCH) has emerged as a general-purpose micromechanics code capable of predicting the effective properties of heterogeneous materials and recovering the local fields. The objective of this paper is to propose a micromechanics approach enabling VAMUCH to homogenize elasto-viscoplastic heterogeneous materials. An affine formulation of the constitutive relations for an elasto-viscoplastic constituent, which exhibits viscoplastic anisotropy and combined isotropic–kinematic hardening, is derived. The weak form of the problem is derived using an asymptotic method, discretized using finite elements, and implemented into VAMUCH. The new features of VAMUCH are validated with examples such as homogenizing binary, fiber-reinforced, and particle-reinforced composites. VAMUCH is found to be capable of handling various microstructure, complex material models, complex loading conditions, and complex loading paths. More sophisticated material models can be implemented into it.     

A micromechanics analysis for the microstructure design of a two-phase pseudoelastic composite

A micromechanics analysis on the possibility of designing a two-phase pseudoelastic composite is made for the case where ductile transformable shape memory alloy plastic particles are imbedded coherently in an elastic matrix. It is demonstrated that a pseudoelastic stress-strain loop in a macroscopic loading-unloading cycle can be obtained by microscopically stress induced forward and reverse martensitic transformations in the SMA particles. The relation between the macroscopic stress-strain response and the material parameters of the constituents of this composite is quantified through the micromechanics calculations, which reveals that the best ductility and thus the greatest energy absorption capacity of this novel microstructure can be obtained by the optimum material design.     

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

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

Photovisco-elastoplastic behavior of polycarbonate material under creep and tension tests

Polymers are widely used as photomechanical models of a prototype material (often a metal). Photoplasticity is one of the methods used in order to show the behavior of plastic materials stressed beyond the linear elastic limit. To illustrate this process we have analyzed the photovisco-elastoplastic behavior of polycarbonate as a photoplastic material. In this paper a technique for local and simultaneous measurement of birefringence and principal strains is presented. The mechanical and optical properties, at room temperature, have been evaluated by means of uniaxial tension tests. A series of creep tests has been carried out in order to study the photovisco-elastoplastic behavior of polycarbonate. In two different experiments we analyzed nonlinear birefringence and the amplitude of the corresponding strains. We could thus evaluate the distribution of strains and the distribution of uniaxial stress for each birefringence state and vice versa.     

Domain switching and non-linear coupling of ferroelectric composites

One of the most notable characteristics of ferroelectric materials is that they could undergo spontaneous polarization and spontaneous strain changes by applied fields. Reorientation of the spontaneous polarization and spontaneous polarization strain of ferroelectric inclusions in ferroelectric composites can change microstructures and affect effective electroelastic properties of ferroelectric composites. Based on orientation distribution function and its evolution as well as switching criterion, non-linear electromechanical coupling behaviour of ferroelectric composites is studied by application of micromechanics. A constitutive model of ferroelectric composites is developed. Comparison between analytical and experimental results shows that the model presented can describe many non-linear electromechanical coupling problems of ferroelectric composites such as polarization or depolarization, etc.     

Microstructure-dependent nonlinear viscoelasticity due to extracellular flow within cellular structures

A constitutive model that highlights a special viscoelastic property of materials with cellular microstructures is developed. We model the microstructure as a regularly arranged system of the same elastic cells that are mutually interconnected by elastic linkages. The space between cells is filled by a fluid that may flow freely within this extracellular space. The macroscopic behavior of the whole structure is studied by means of continuum mechanics using a differential scheme with internal variables. Here, the internal variables are chosen as the distances that separate neighboring cells. The evolution equations are derived from the Clausius–Planck inequality, which considers the internal dissipation to be exclusively due to the extracellular fluid movement. Special attention is paid to incompressible materials in the context of uniaxial load. In this context, the importance of the fluid viscosity on material behavior is related to microstructural parameters like the cells’ dimensions and the relative stiffness between the cells and matrix elastic reinforcement.     

Velocity dispersion in an elastic plate with microstructure: effects of characteristic length in a couple stress model

Microstructural effects become important, when dimensions of the heterogeneous material are comparable to the length scale of microstructure and the state of stress needs to be defined in a non-local manner. Linear theory of elasticity, which is associated with the concept of homogeneity of material and local stresses, cannot describe the behavior of the materials with microstructures. In this study, Couple stress theory of elasticity has been employed to capture the size effects on the propagation of Lamb waves in an elastic plate with microstructure. Effects on the dispersion curves of Lamb waves are studied, when the characteristic length of the material is comparable to cell size. The governing equations of couple stress theory, involving stresses and couple stresses are solved to study the impact of different characteristic lengths, comparable with cell size. Since bone is a material with microstructure, so for numerical calculations and graphical representation of the results, the plate is considered to have mechanical properties typically used for bones.     

A phase field model incorporating strain gradient viscoplasticity: Application to rafting in Ni-base superalloys

The first formulation of a phase field model accounting for size-dependent viscoplasticity is developed to study materials in which microstructure evolution and viscoplastic behavior are strongly coupled. Plasticity is introduced using a continuum strain gradient formalism which captures the size effect of the viscoplastic behavior. First, the influence of this size effect on the mechanical behavior of the material is discussed in static microstructures. Then, the dynamic coupling between microstructure evolution and viscoplastic activity is addressed and illustrated by the rafting of the microstructure observed in Ni-base superalloys under creep conditions. It is found that the plastic size effect has only a moderate impact on the shape of the rafts but is crucial to reproduce the macroscopic mechanical behavior of that particular material.     

A Study of Large Plastic Deformations in Dual Phase Steel Using Digital Image Correlation and FE Analysis

Large plastic deformation in sheets made of dual phase steel DP800 is studied experimentally and numerically. Shear testing is applied to obtain large plastic strains in sheet metals without strain localisation. In the experiments, full-field displacement measurements are carried out by means of digital image correlation, and based on these measurements the strain field of the deformed specimen is calculated. In the numerical analyses, an elastoplastic constitutive model with isotropic hardening and the Cockcroft–Latham fracture criterion is adopted to predict the observed behaviour. The strain hardening parameters are obtained from a standard uniaxial tensile test for small and moderate strains, while the shear test is used to determine the strain hardening for large strains and to calibrate the fracture criterion. Finite Element (FE) calculations with shell and brick elements are performed using the non-linear FE code LS–DYNA. The local strains in the shear zone and the nominal shear stress-elongation characteristics obtained by experiments and FE simulations are compared, and, in general, good agreement is obtained. It is demonstrated how the strain hardening at large strains and the Cockcroft–Latham fracture criterion can be calibrated from the in-plane shear test with the aid of non-linear FE analyses. An erratum to this article can be found at     

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.     

W25Fe25Ni25Mo25高熵合金高速侵彻细观结构演化特性

为了探究W₂₅Fe₂₅Ni₂₅Mo₂₅高熵合金弹体在侵彻过程中宏观变形行为与材料微细观结构之间的联系, 基于对两相流动模型的简化, 建立了考虑软、硬相密度、流速以及浓度差异的等截面直管两相流动演化模型. 类比宏观状态下侵彻弹体头部材料的流入流出特性, 选定分析区域, 建立两相细观结构下材料在分析区域的流入流出关系, 再结合细观结构演化方程, 给出了分析区域中浓度演化结果, 提出了表征材料浓度演化速率的流动稳定系数t/l_length. 为了对比不同细观结构弹体的侵彻行为, 选取典型两相材料钨铜合金(W₇₀Cu₃₀), 基于小口径弹道枪发射平台开展两种弹体侵彻半无限钢靶试验, 对比两种合金弹体细观结构演化行为. 结果表明, 硬相浓度分布总体上体现“中心浓, 边缘稀”的特点; 硬相的浓度越高, 密度越大, 驱动速度越快, 则流动稳定系数t/l_length值越小, 侵彻过程中弹体的流动稳定性越好, 弹体头部材料越容易形成连续的塑性流动带. 等截面直管两相流动演化模型可用于描述侵彻过程中弹体头部材料的流动稳定性, 揭示了侵彻过程中弹体头部变形与细观两相结构之间的关联机制.