伴随微孔洞生长的裂纹弹塑性定常扩展

裂纹在韧性材料中扩展时,将们随着微孔洞的萌生和生长,孔洞的萌生和深化将直接影响着材料的总体断裂韧性和强度,以往的研究主要集中在将裂纹的扩展刻划为微孔洞的萌生、生长和汇合这样一个过程。从传统的断裂过程区模型出发研究微孔洞的萌生和生长对材料总体断裂韧性的影响,通过采用Ｇｕｒｓｏｎ模型,建立塑性增量本构关系,然后针对定常扩展情况直接进行分析,孔洞对材料断裂韧性的影响由本构关系刻划,而在孔洞汇合模型中,上     

Dynamic void nucleation and growth in solids: A self-consistent statistical theory

We present a framework for a self-consistent theory of spall fracture in ductile materials, based on the dynamics of void nucleation and growth. The constitutive model for the material is divided into elastic and “plastic” parts, where the elastic part represents the volumetric response of a porous elastic material, and the “plastic” part is generated by a collection of representative volume elements (RVEs) of incompressible material. Each RVE is a thick-walled spherical shell, whose average porosity is the same as that of the surrounding porous continuum, thus simulating void interaction through the resulting lowered resistance to further void growth. All voids nucleate and grow according to the appropriate dynamics for a thick-walled sphere made of incompressible material. The macroscopic spherical stress in the material drives the response in all volume elements, which have a distribution of critical stresses for void nucleation, and the statistically weighted sum of the void volumes of all RVEs generates the global porosity. Thus, macroscopic pressure, porosity, and a distribution of growing microscopic voids are fully coupled dynamically. An example is given for a rate-independent, perfectly plastic material. The dynamics of void growth gives rise to a rate effect in the macroscopic material even though the parent material is rate independent.     

Two mechanisms of ductile fracture: void by void growth versus multiple void interaction

Two distinct mechanisms of crack initiation and advance by void growth have been identified in the literature on the mechanics of ductile fracture. One is the interaction a single void with the crack tip characterizing initiation and the subsequent void by void advance of the tip. This mechanism is represented by the early model of Rice and Johnson and the subsequent more detailed numerical computations of McMeeking and coworkers on a single void interacting with a crack tip. The second mechanism involves the simultaneous interaction of multiple voids on the plane ahead of the crack tip both during initiation and in subsequent crack growth. This mechanism is revealed by models with an embedded fracture process zone, such as those developed by Tvergaard and Hutchinson. While both mechanisms are based on void nucleation, growth and coalescence, the inferences from them with regard to crack growth initiation and growth are quantitatively different. The present paper provides a formulation and numerical analysis of a two-dimensional plane strain model with multiple discrete voids located ahead of a pre-existing crack tip. At initial void volume fractions that are sufficiently low, initiation and growth is approximately represented by the void by void mechanism. At somewhat higher initial void volume fractions, a transition in behavior occurs whereby many voids ahead of the tip grow at comparable rates and their interaction determines initiation toughness and crack growth resistance. The study demonstrates that improvements to be expected in fracture toughness by reducing the population of second phase particles responsible for nucleating voids cannot be understood in terms of trends of one mechanism alone. The transition from one mechanism to the other must be taken into account.     

Applications of planar anisotropic yield criteria to porous sheet metal forming simulations

Planar anisotropic yield functions, with rounded vertexes especially near the equal-biaxial direction of the corresponding yield loci, appropriate for some structure metals are employed for the matrix surrounding voids in the present study. The widely adopted Hill anisotropic yield functions are also implemented into the matrix for comparisons. Mechanisms of the void growth, void nucleation, and void coalescence are simultaneously considered here. Effects of the yield function of the corresponding matrix on the sheet metal under two typical sheet forming operations, a hemispherical punch stretching operation and a cup drawing operation, are investigated via a finite element analysis. Thickness strains in various orientations of the sheet are then evaluated. Simulation results show that the yield function of the corresponding matrix plays important roles on the strain distribution and the strain localization as well. Early localization would be found for the sheet with relatively small initial void volume fraction in two operations. Yield functions of the matrix rather influence the earing phenomenon under the cup drawing procedure even similar displacement profiles of the outer boundary could be observed.     

A non-linear viscoelastic model for ceramics at high temperatures

High-temperature mechanical behavior of ceramics is characterized by non-linear rate dependent responses, asymmetric behavior in tension and compression, and nucleation and coalescence of voids leading to rupture. Moreover, rupture experiments show considerable scatter or randomness in fatigue lives of nominally equal specimens. To capture the non-linear, asymmetric time-dependent behavior, a new non-linear viscoelastic model is proposed. Non-linearity and asymmetry are introduced in the volumetric component. To model the random formation and coalescence of voids, each element is assigned a failure strain sampled from a lognormal distribution. An element is deleted when its volumetric strain exceeds its failure strain. Temporal increases in strains produce a sequential loss of elements (a model for void nucleation and growth), which in turn leads to failure. Non-linear viscoelastic model parameters are determined from uniaxial tensile and compressive creep experiments on silicon nitride. The model is then used to predict the deformation of four-point bending and ball-on-ring specimens. Simulation is used to predict statistical moments of rupture lives. Numerical simulation results compare well with results of four-point bending experiments.     

A physical model for nucleation and early growth of voids in ductile materials under dynamic loading

Spall fracture and other rapid tensile failures in ductile materials are often dominated by the rapid growth of voids. Recent research on the mechanics of void growth clearly shows that void nucleation may be represented as a bifurcation phenomenon, wherein a void forms spontaneously followed by highly localized plastic flow around the new void. Although thermal, viscoplastic, and work hardening effects all play an essential role in the earliest stages of nucleation and growth, the flow becomes dominated by spherical radial inertia, which soon causes all voids to grow asymptotically at the same rate, regardless of differences in initial conditions or constitutive details, provided only that there is the same density of matrix material and the same excess loading history beyond the cavitation stress.These two facts, initiation by bifurcation at a cavitation stress, at which a void first appears, and rapid domination by inertia, are used to postulate a simple, but physically realistic, model for nucleation and early growth of voids in a ductile material under rapid tensile loading. A reasonable statistical distribution for the cavitation stress at various nucleation sites and a simple similarity solution for inertially dominated void growth permit a simple calculation of the initiation and early growth of porosity in the material.Parametric analyses are presented to show the effect that loading rate, peak loading stress, density of nucleation sites, physical properties of the material, etc. have on the applied pressure and distribution of void sizes when a critical porosity is reached.     

Creep cavitation in metals

This review concisely describes the state-of-the-art of the understanding of cavity, or r-type void, formation during stages I and II (primary and secondary) creep in polycrystalline metals and alloys, particularly at elevated temperatures. These cavities can directly lead to Stage III, or tertiary, creep and the eventual failure of metals. There have been, in the past, a variety of creep fracture reviews that omitted important developments relevant to creep cavitation or are less than balanced in their discussions of conflicting ideas or theories regarding various aspects of cavity nucleation and growth. This concise, comprehensive, review discusses all of the important developments over the past several decades relating to both the nucleation and growth of cavities. The nucleation section discusses the details and limitations of the approaches based on “classic” nucleation theory, slip-induced nucleation as well as grain boundary sliding effects. Growth is discussed starting from the Hull–Rimmer diffusion controlled cavity growth (DCCG) model. This will be followed by refinements to DCCG by others. Next, there will be a discussion of plastic cavity growth and diffusion-plasticity coupling theories. This will be followed by the particularly important development of constrained cavity growth, initially proposed by Dyson, and probably under-appreciated. Other growth effects by grain boundary sliding will also be discussed. All of these mechanisms will be compared with their predictions in terms of creep fracture phenomenology such as the Monkman–Grant relationship. Finally, there will be a discussion of creep crack propagation by cavitation ahead of the crack tip.     

Modeling stress state dependent damage evolution in a cast Al–Si–Mg aluminum alloy

Internal state variable rate equations are cast in a continuum framework to model void nucleation, growth, and coalescence in a cast Al–Si–Mg aluminum alloy. The kinematics and constitutive relations for damage resulting from void nucleation, growth, and coalescence are discussed. Because damage evolution is intimately coupled with the stress state, internal state variable hardening rate equations are developed to distinguish between compression, tension, and torsion straining conditions. The scalar isotropic hardening equation and second rank tensorial kinematic hardening equation from the Bammann–Chiesa–Johnson (BCJ) Plasticity model are modified to account for hardening rate differences under tension, compression, and torsion. A method for determining the material constants for the plasticity and damage equations is presented. Parameter determination for the proposed phenomenological nucleation rate equation, motivated from fracture mechanics and microscale physical observations, involves counting nucleation sites as a function of strain from optical micrographs. Although different void growth models can be included, the McClintock void growth model is used in this study. A coalescence model is also introduced. The damage framework is then evaluated with respect to experimental tensile data of notched Al–Si–Mg cast aluminum alloy specimens. Finite element results employing the damage framework are shown to illustrate its usefulness.     

A constitutive model for elastoplastic solids containing primary and secondary voids

In many ductile metallic alloys, the damage process controlled by the growth and coalescence of primary voids nucleated on particles with a size varying typically between 1 and 100 μm, is affected by the growth of much smaller secondary voids nucleated on inclusions with a size varying typically between 0.1 and 3 μm. The goal of this work is first to quantify the potential effect of the growth of these secondary voids on the coalescence of primary voids using finite element (FE) unit cell calculations and second to formulate a new constitutive model incorporating this effect. The nucleation and growth of secondary voids do essentially not affect the growth of the primary voids but mainly accelerate the void coalescence process. The drop of the ductility caused by the presence of secondary voids increases if the nucleation strain decreases and/or if their volume fraction increases and/or if the primary voids are flat. A strong coupling is indeed observed between the shape of the primary voids and the growth of the second population enhancing the anisotropy of the ductility induced by void shape effects. The new micromechanics-based coalescence condition for internal necking introduces the softening induced by secondary voids growing in the ligament between two primary voids. The FE cell calculations were used to guide and assess the development of this model. The use of the coalescence condition relies on a closed-form model for estimating the evolution of the secondary voids in the vicinity of a primary cavity. This coalescence criterion is connected to an extended Gurson model for the first population including the effect of the void aspect ratio. With respect to classical models for single void population, this new constitutive model improves the predictive potential of damage constitutive models devoted to ductile metal while requiring only two new parameters, i.e. the initial porosity of second population and a void nucleation stress, without any additional adjustment.     

Modelling of dynamic ductile fracture and application to the simulation of plate impact tests on tantalum

A model of dynamic damage by void nucleation and growth is proposed for elastic-viscoplastic materials sustaining intense loading. The model is dedicated to ductile materials for which fracture is caused by microvoiding. The material contains potential nucleation sites where microvoids are generated when the local pressure overcomes the nucleation pressure. A probability density function is adopted to describe the fluctuation of the nucleation pressure within the material. The void growth is described by using a hollow sphere model where micro-inertia effects are accounted for. The matrix weakening due to void growth is also included.The model has been first tested under uniaxial deformation. When the strain rate is assumed constant, the pressure inside the material has nearly a linear response up to a maximum. An analytical expression for the maximum pressure is proposed.Finite element simulations of plate impact tests have been carried out and compared to experiments on tantalum. From simulations based on the proposed model, an increase of the spall strength is observed with higher shock intensities. Therefore, the relationship between the velocity pullback and spall strength usually assumed in the literature (based on the acoustic approach) seems to be inadequate. Velocity profiles are simulated for different flyer thicknesses and different impact velocities with close agreement with experiments.     

Prediction of central bursting during axisymmetric cold extrusion of a metal alloy containing particles

A simple cell model based damage dependent yield surface is used to model the effect of void nucleation and growth in an aluminum alloy during an axisymmetric cold extrusion process. Material parameters for characterization of the yield surface are determined through a physically consistent micromechanical cell modeling technique. The model can account for the behavior of a void containing a particle under severe compressive processing conditions. The formation of distinct, equally spaced, arrowhead shaped ‘central burst’ defects is observed during simulation of the extrusion process. Application of the model to a two-stage rolling process is also briefly illustrated. Formation of central bursts during extrusion and edge cracking during rolling is explained in terms of the hydrostatic stress distribution and the related void growth. The affects of material hardening, surface friction and die geometry are examined in the case of extrusion. Correlation is found between the simulations and analytical and experimental results, confirming the suitability of the constitutive model.     

A void nucleation and growth based damage framework to model failure initiation ahead of a sharp notch in irradiated bcc materials

A continuum damage framework is developed and coupled with an existing crystal plasticity framework, to model failure initiation in irradiated bcc polycrystalline materials at intermediate temperatures. Constitutive equations for vacancy generation due to inelastic deformation, void nucleation due to vacancy condensation, and diffusion-assisted void growth are developed. The framework is used to simulate failure initiation at dislocation channel interfaces and grain boundaries ahead of a sharp notch. Evolution of the microstructure is considered in terms of the evolution of inelastic deformation, vacancy concentration, and void number density and radius. Evolution of the damage, i.e., volume fraction of the voids, is studied as a function of applied deformation. Effects of strain rate and temperature on failure initiation are also studied. The framework is used to compute the fracture toughness of irradiated specimens for various loading histories and notch geometries. Crack growth resistance of the irradiated specimens are computed and compared to that of virgin specimens. Results are compared to available experimental data.     

Grain–boundary interactions and orientation effects on crack behavior in polycrystalline aggregates

A dislocation-density grain–boundary interaction scheme has been developed to account for the interrelated dislocation-density interactions of emission, absorption and transmission in GB regions. The GB scheme is based on slip-system compatibility, local resolved shear stresses, and immobile and mobile dislocation-density accumulation at critical GB locations. To accurately represent dislocation-density evolution, a conservation law for dislocation-densities is used to balance dislocation-density absorption, transmission and emission from the GB. The behavior of f.c.c. polycrystalline copper, with different random low and high angle GBs, are investigated for different crack lengths. For aggregates with random low angle GBs, dislocation-density transmission dominates at the GBs, which can indicate that the low angle GB will not significantly change crack growth directions. For aggregates with random high angle GBs, extensive dislocation-density absorption and pile-ups occur. The high stresses associated with this behavior, along the GBs, can result in intergranular crack growth due to potential crack nucleation sites in the GB.     

Numerical investigation of condensing steam flow in boundary layers

The paper describes a numerical method for the prediction of condensing steam flow within compressible boundary layers. The method is based on a simple stream function technique, which enables straightforward integration of the nucleation and droplet growth equations in a Lagrangian frame of reference. Calculations show how viscous dissipation and reduced expansion rate within a typical boundary layer influence nucleation and growth, leading to droplet radii and size distributions that differ substantially from those predicted in inviscid flow. The impact of condensation on temperature and velocity profiles, and the implications for thermodynamic loss are also considered.     

Site of ductile fracture initiation in cold forging: A finite element model

Ductile fracture occurs due to micro-void nucleation, growth and finally coalescence into micro-crack. In this study a new ductile fracture condition that based on the microscopic phenomena of void nucleation, growth and coalescence was proposed. Using this condition and combining with finite element model to predict the fracture locations in bulk metal forming, the results show that it is in close accordance with observations of some experimental specimens. However, in order to obtaining the high trustiness many experiments have to be carried out.     

Determination of Ductile Damage Parameters by Local Deformation Fields: Measurement and Simulation

This work comprises the development, implementation and application of methods for the parameter identification of damage mechanical constitutive laws. Ductile damage is described on a continuum mechanical basis by extension of the von Mises yield condition with the Gurson–Tvergaard–Needleman as well as with the Rousselier model. The classical Rousselier model is complemented by accelerated void growth and void nucleation. The non-linear boundary and initial value problem is solved by the finite element system SPC–PMHP, which was developed in the frame of the special research program SFB393 for parallel computers. The material parameters are identified by locally measured displacement fields and measured force–displacement curves. For the material parameter identification a non-linear optimization algorithm is used, which renders the objective function to a minimum by means of a gradient based method. A useful strategy to identify the material parameters was found by careful numerical studies. Finally, using the object grating method the local displacement fields as well as the force–displacement curves are measured at notched flat bar tension specimens made of StE 690 and the parameters of the material are identified.     

磁电弹复合材料和刚性导电导磁圆柱压头的二维微动接触分析

本文求解平面应变状态下磁电弹复合材料半平面和刚性导电导磁圆柱压头的二维微动接触问题。假设压头具有良好的导电导磁性,且表面电势和磁势是常数。微动接触依赖载荷的加载历史,所以首先求解单独的法向加载问题,然后在法向加载问题的基础上求解循环变化的切向加载问题。整个接触区可以分为内部的中心粘着区和两个外部的滑移区,其中滑移区满足Coulomb摩擦法则。利用Fourier积分变换,磁电弹半平面的微动接触问题将简化为耦合的Cauchy奇异积分方程组,然后数值离散为线性代数方程组,利用迭代法求解未知的粘着/滑移区尺寸、电荷分布、磁感应强度、法向接触压力和切向接触力。数值算例给出了摩擦系数、总电荷和总磁感应强度对各加载阶段的表面接触应力、电位移和磁感应强度的影响。     

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.     

Modeling of ductile fracture: Application of the mechanism-based concepts

This paper summarizes our recent studies on modeling ductile fracture in structural materials using the mechanism-based concepts. We describe two numerical approaches to model the material failure process by void growth and coalescence. In the first approach, voids are considered explicitly and modeled using refined finite elements. In order to predict crack initiation and propagation, a void coalescence criterion is established by conducting a series of systematic finite element analyses of the void-containing, representative material volume (RMV) subjected to different macroscopic stress states and expressed as a function of the stress triaxiality ratio and the Lode angle. The discrete void approach provides a straightforward way for studying the effects of microstructure on fracture toughness. In the second approach, the void-containing material is considered as a homogenized continuum governed by porous plasticity models. This makes it possible to simulate large amount of crack extension because only one element is needed for a representative material volume. As an example, a numerical approach is proposed to predict ductile crack growth in thin panels of a 2024-T3 aluminum alloy, where a modified Gologanu–Leblond–Devaux model [Gologanu, M., Leblond, J.B., Devaux, J., 1993. Approximate models for ductile metals containing nonspherical voids – Case of axisymmetric prolate ellipsoidal cavities. J. Mech. Phys. Solids 41, 1723–1754; Gologanu, M., Leblond, J.B., Devaux, J., 1994. Approximate models for ductile metals containing nonspherical voids – Case of axisymmetric oblate ellipsoidal cavities. J. Eng. Mater. Tech. 116, 290–297; Gologanu, M., Leblond, J.B., Perrin, G., Devaux, J., 1995. Recent extensions of Gurson’s model for porous ductile metals. In: Suquet, P. (Ed.) Continuum Micromechanics. Springer-Verlag, pp. 61–130] is used to describe the evolution of void shape and void volume fraction and the associated material softening, and the material failure criterion is calibrated using experimental data. The calibrated computational model successfully predicts crack extension in various fracture specimens, including the compact tension specimen, middle crack tension specimens, multi-site damage specimens and the pressurized cylindrical shell specimen.     

金属基复合材料和强度与损伤分析

用观察计算力学的方法分析了金属基复合材料（ＭＭＣ）多重损伤与强度的关系,采用唯象的内聚力模型模拟纤维／基体界面的脱粘和采用Ｇ－Ｔ模型描述韧性基体的损伤。并用上述模型分析了长纤维增强ＭＭＣ在横向荷载作用下损伤演化的规律,讨论了不同界面性质与材料强度及损伤、破坏模式之间的关系。