Coupling effects of void size and void shape on the growth of prolate ellipsoidal microvoid

The combined effects of void size and void shape on the void growth are studied by using the classical spectrum method. An infinite solid containing an isolated prolate spheroidal void is considered to depict the void shape effect and the Fleck-Hutchinson phenomenological strain gradient plasticity theory is employed to capture the size effects. It is found that the combined effects of void size and void shape are mainly controlled by the remote stress triaxiality. Based on this, a new size-dependent void growth model similar to the Rice-Tracey model is proposed and an important conclusion about the size-dependent void growth is drawn: the growth rate of the void with radius smaller than a critical radius r_c may be ignored. It is interesting that r_c is a material constant independent of the initial void shape and the remote stress triaxiality.The project supported by the National Natural Science Foundation of China (A10102006) and the New Century Excellent Talents in Universities of China. The English text was polished by Keren Wang.     

A note on the effect of surface energy and void size to void growth

Recent studies revealed that rapid void growth is the dominant failure mechanism in an elasto-plastic solid under high mean tensile stress. This paper studies the effect of the surface energy and void size to the void growth. The models of a thick spherical shell and a thick spherical column in void growth are analyzed and numerically estimated. The main conclusion from this study is that, for typical metals, the surface energy effect is negligible for voids larger than 100 nm in size, but it may become significant when the void size is on the order of 10 nm.     

The effects of void cluster size on ductile fracture

The effects of void clustering on ductile fracture are studied by modeling a discrete set of randomly distributed clusters. Each cluster consists of four, equally-spaced, cylindrical voids. The spacing between the clusters is held constant while the spacing between the voids is varied. A Eulerian finite element program is used to numerically solve the boundary value problems. A salient feature of the previous investigations is that both the ultimate stress and the fracture strain are functions of the void distribution. In contrast, the ultimate stress remains constant while the fracture strain changes with the void cluster diameter in the current investigation.     

Micromechanical finite element calculations of temperature and void configuration effects on void growth and coalescence

We present micromechanical finite element results that quantify coalescence effects based upon temperature and different spatial arrangements of voids. We propose a critical intervoid ligament distance (ILD) to define void coalescence that is derived from micromechanical simulations in which void volume fraction evolves as a function of strain. Several parameters were varied using the temperature and strain rate internal variable plasticity model of Bammann–Chiesa–Johnson to determine the coalescence effects. The parameters include two types of materials with different work hardening rates (304L stainless steel and 6061T6 aluminum), three different temperatures (298, 400, and 600 K), several boundary conditions (force and displacement: uniaxial, plane strain, and biaxial), type of element used (plane strain and axisymmetric), different ILDs, and the number of voids (one and two void configurations). The present study provides a basis for macroscale modeling of coalescence which is briefly discussed.     

Vapor pressure and void size effects on failure of a constrained ductile film

To achieve certain properties, semiconductor adhesives and molding compounds are made by blending filler particles with polymer matrix. Moisture collects at filler particle/polymer matrix interfaces and within voids of the composite. At reflow temperatures, the moisture vaporizes. The rapidly expanding vapor creates high internal pressure on pre-existing voids and particle/matrix interfaces. The simultaneous action of thermal stresses and internal vapor pressure drives both pre-existing and newly nucleated voids to grow and coalesce causing material failure. Particularly susceptible are polymeric films and adhesives joining elastic substrates, e.g. Ag filled epoxy. Several competing failure mechanisms are studied including: near-tip void growth and coalescence with the crack; extensive void growth and formation of an extended damaged zone emanating from the crack; and rapid void growth at highly stressed sites at large distances ahead of the crack, leading to multiple damaged zones. This competition is driven by the interplay between stress elevation induced by constrained plastic flow and stress relaxation due to vapor pressure assisted void growth.A model problem of a ductile film bonded between two elastic substrates, with a centerline crack, is studied. The computational study employs a Gurson porous material model incorporating vapor pressure effects. The formation of multiple damaged zones is favored when the film contains small voids or dilute second-phase particle distribution. The presence of large voids or high vapor pressure favor the growth of a self-similar damage zone emanating from the crack. High vapor pressure accelerates film cracking that can cause device failures.     

FCC晶体中不均匀变形对空洞长大的影响

通过编制率相关有限元用户子程序,采用一个单胞模型研究了FCC晶体中孔洞在单晶及晶界的长大行为,分析了由于晶体取向及变形失配对孔洞长大和聚合的影响。研究结果表明：孔洞的形状和长大方向与晶体取向密切相关;晶界上孔洞的长大速度大于单晶中孔洞的长大速度;晶粒间的变形失配加速了晶界上孔洞的长大趋势,因而使材料易发生沿晶断裂,随着晶粒间取向因子差异的增加,孔洞越易沿着晶界长大。     

Lattice orientation effects on void growth and coalescence in fcc single crystals

Void growth and coalescence in fcc single crystals were studied using crystal plasticity under uniaxial and biaxial loading conditions and various orientations of the crystalline lattice. A 2D plane strain unit cell with one and two cylindrical voids was employed using three-dimensional 12 potentially active slip systems. The results were compared to five representative orientations of the tensile axis on the stereographic triangle. For uniaxial tension conditions, the void volume fraction increase under the applied load is strongly dependent on the crystallographic orientation with respect to the tensile axis. For some orientations of the tensile axis, such as [1 0 0] or [1 1 0], the voids exhibited a growth rate twice as fast compared with other orientations ([1 0 0], [2 1 1]). Void growth and coalescence simulations under uniaxial loading indicated that during deformation along some orientations with asymmetry of the slip systems, the voids experienced rotation and shape distortion, due mainly to lattice reorientation. Coalescence effects are shown to diminish the influence of lattice orientation on the void volume fraction increase, but noteworthy differences are still present. Under biaxial loading conditions, practically all differences in the void volume fraction for different orientations of the tensile axes during void growth vanish. These results lead to the conclusion that at microstructural length scales in regions under intense biaxiality/triaxiality conditions, such as crack tip or notched regions, the plastic anisotropy due to the initial lattice orientation has only a minor role in influencing the void growth rate. In such situations, void growth and coalescence are mainly determined by the stress triaxiality, the magnitude of accumulated strain, and the spatial localization of such plastic strains.     

Experimental studies on the yield behavior of ductile and brittle aluminum foams

Specimen size effects are a major cause of the unreliability of foam models in finite element codes. Here, the modified Arcan apparatus is used to investigate the biaxial yielding of ductile and brittle Al foams. This apparatus subjects a central section of a “butterfly-shaped” specimen to a uniform state of plane stress. The stresses have local maxima at the central section, thus ensuring that yielding occurs there. A yield envelope, which directly relates to the crushing process, can then be determined. Size effects are introduced when using conventional methods such as tri-axial or plate-shear tests. In such tests, averages of stress and strain are measured. These measures do not represent the actual yield event, because foam's internal structure is inhomogeneous and so is the deformation field. Strain localization and failure can occur at any weak layer of cells in the bulk. In this study, we have performed a series of biaxial tests on isotropic Alporas and anisotropic Hydro closed-cell Al foams of approximately equal densities. Alporas failed locally by a ductile phenomenon of progressive crushing of cells. It also possessed uniaxial strength asymmetry. Hydro specimens parallel and perpendicular to ‘foam rise’ were investigated. The Hydro foam developed a local, characteristic brittle crack at loads in the vicinity of the yield point. Phenomenological yield surfaces, which incorporate these features are obtained for the foams, and show dependence on both the deviatoric and hydrostatic stresses. We also provide expressions for the shear and hydrostatic strengths in terms of the uniaxial strengths. Finally, the size-independence of the yield surface is verified using the uniaxial compression of tapered specimens.     

A microscopic damage model considering the change of void shape and application in the void closing

A microscopic damage model of ellipsoidal body containing ellipsoidal void for nonlinear matrix materials is developed under a particular coordinate. The change of void shape is considered in this model. The viscous restrained equation obtained from the model is affected by stress ∑_(ij), void volume fraction f, material strain rate exponent m as well as the void shape. Gurson's equation is modified from the numerical solution. The modified equation is suitable for the case of nonlinear matrix materials and changeable voids. Lastly, the model is used to analyze the closing process of voids.     

Influences of microscope size effect of matrix on plastic potential and void growth of porous materials

Based on approximate theoretical analyses on a typical spherical cell containing a spherical microvoid, the influences of matrix materials' microscopic scale on the macroscopic constitutives potential theory of porous material and microvoid growth have been investigated in detail. By assuming that the plastic deformation behavior of matrix materials follows the strain gradient (SG) plastic theory involving the stretch and rotation gradients, the ratio (λ=l/a) of the matrix materials' intrinsic characteristic lengthl to the microvoid radiusa is introduced into the plastic constitutives potential and the void growth law. The present results indicate that, when the radiusa of microvoids is comparable with the intrinsic characteristic lengthl of the matrix materials, the influence of microscopic size effect on neither the constitutive potential nor the microvoid evolution predicted can be ignored. And when the void radiusa is much lager than the intrinsic characteristic lengthl of the matrix materials, the present model can retrogress automatically to the improved Gurson model that takes into account the strain hardening effect of matrix materials. Project supported by the National Natural Science Foundation of China (No. A10102006).     

Effects of void band orientation and crystallographic anisotropy on void growth and coalescence

The effects of void band orientation and crystallographic anisotropy on void growth and linkage have been investigated. 2D model materials were fabricated by laser drilling a band of holes into the gage section of sheet tensile samples using various orientation angles with respect to the tensile axis normal. Both copper and magnesium sheets have been studied in order to examine the role of crystallographic anisotropy on the void growth and linkage processes. The samples were pulled in uniaxial tension inside the chamber of an SEM, enabling a quantitative assessment of the growth and linkage processes. The void band orientation angle has a significant impact on the growth and linkage of the holes in copper. As the void band orientation angle is increased from 0° to 45°, the processes of coalescence and linkage are delayed to higher strain values. Furthermore, the mechanism of linkage changes from internal necking to one dominated by shear localization. In contrast, the void band orientation does not have a significant impact on the void growth and linkage processes in magnesium. Void growth in these materials occurs non-uniformly due to interactions between the holes and the microstructure. The heterogeneous nature of deformation in magnesium makes it difficult to apply a coalescence criterion based on the void dimensions. Furthermore, the strain at failure does not show a relationship with the void band orientation angle. Failure associated with twin and grain boundaries interrupts the plastic growth of the holes and causes rapid fracture. Therefore, the impact of the local microstructure outweighs the effects of the void band orientation angle in this material.     

Experimental studies on yield behavior of steel under rapid heating

The yield behavior of steels under different heating rates and pre-stress levels is studied experimentally. The yield temperatures are found to increase with the heating rates and decrease with the pre-stress values. Two remarkable changes in the microstructures of tested steels are the grain elongation along the loading direction owing to the plastic deformation, and the grain refineries owing to the recrystallization. The micrographic features prove the above remarks. An equation of plastic strain rates of metals can estimate the yield temperature of metals under various heating rates. An extended Johnson–Mehl equation of recrystallization is derived, which includes the influence of heating rates.     

Effect of void shape on the macroscopic response of non-linear porous solids

The macroscopic response of an incompressible power-law matrix containing aligned spheroidal voids is investigated. The voids are assumed to be arranged in a uniform array, and the response of the solid is evaluated by isolating a typical block of the material containing a single void. The requisite boundary value problem for this “unit cell” is solved using a spectral method which is an adaption of that used by Lee and Mear “Axisymmetric Deformation of Power-law Solids containing Elliptical Inhomogeneities. Part I: Rigid Inclusions”, J. Mech. Phys. Solids, (1992) 8, 1805. Attention is restricted to axisymmetric deformation, and results for the macroscopic strain-rates (or strains) are presented for a range of void shape, void volume concentrations, hardening exponents and remote stress triaxilities.     

Effect of yield surface curvature and void nucleation on plastic flow localization

The effect of void nucleation is incorporated in a recently proposed material model that accounts for a combination of kinematic hardening and isotropic hardening of a porous ductile material. Since each of plastic dilatancy, void nucleation and yield surface curvature have a strong influence on predictions of plastic flow localization, the present material model can be used to study the interaction of these effects. Nucleation controlled by the plastic strain as well as nucleation controlled by the maximum normal stress on the particle-matrix interface are modelled. The predictions of the material model, for various combinations of parameters, are illustrated by analyses of shear band formation under plane strain or axisymmetric conditions, and by analyses of necking in biaxially stretched sheets.     

Influence of stress triaxiality and temperature on void growth and ductile/brittle transition behavior of 40 Cr steel

Examined experimentally are the influence of stress triaxiality and temperature on the growth of microvoids and the ductile/brittle transition (DBT) macrobehavior of 40 Cr steel subjected to two different heat treatments. This is accomplished by testing more than 300 smooth and notched specimens over a temperature range of 20°C to −196°C. Changes in the microstructure morphology are examined by scanning electron microscopy (SEM) and identified with fracture data on a surface constructed from the uniaxial strain ε_c at fracture, the stress triaxiality R_σ and the temperature T. While stress triaxiality has a significant influence on the DBT temperature T_c, it does not affect the ratio of the average radius of voids R_o to that of inclusions R_i. The ratio R_o/R_i is found to increase with temperature and remains constant in specimens with different notch radii regardless of the temperature. Empirical relations between T_c and R_σ⁻ and R_o/R_i and T are proposed to better understand how macrofracture parameters are influenced by microstructure entities.     

Dynamic multi-axial behavior of shape memory alloy nanowires with coupled thermo-mechanical phase-field models

The objective of this paper is to provide new insight into the dynamic thermo-mechanical properties of shape memory alloy (SMA) nanowires subjected to multi-axial loadings. The phase-field model with Ginzburg–Landau energy, having appropriate strain based order parameter and strain gradient energy contributions, is used to study the martensitic transformations in the representative 2D square-to-rectangular phase transformations for FePd SMA nanowires. The microstructure and mechanical behavior of martensitic transformations in SMA nanostructures have been studied extensively in the literature for uniaxial loading, usually under isothermal assumptions. The developed model describes the martensitic transformations in SMAs based on the equations for momentum and energy with bi-directional coupling via strain, strain rate and temperature. These governing equations of the thermo-mechanical model are numerically solved simultaneously for different external loadings starting with the evolved twinned and austenitic phases. We observed a strong influence of multi-axial loading on dynamic thermo-mechanical properties of SMA nanowires. Notably, the multi-axial loadings are quite distinct as compared to the uniaxial loading case, and the particular axial stress level is reached at a lower strain. The SMA behaviors predicted by the model are in qualitative agreements with experimental and numerical results published in the literature. The new results reported here on the nanowire response to multi-axial loadings provide new physical insight into underlying phenomena and are important, for example, in developing better SMA-based MEMS and NEMS devices     

Extension of the gurson model accounting for the void size effect

A continuum model of solids with cylindrical microvoids is proposed based on the Taylor dislocation model.The model is an extension of Gurson model in the sense that the void size effect is accounted for. Beside the void volume fraction f, the intrinsic material length I becomes a parameter representing voids since the void size comes into play in the Gurson model. Approximate yield functions in analytic forms are suggested for both solids with cylindrical microvoids and with spherical microvoids. The application to uniaxial tension curves shows a precise agreement between the approximate analytic yield function and the “exact” parametric form of integrals.     

Effects of pressure-sensitivity and plastic dilatancy on void growth and interaction

Hydrostatic stress can affect the non-elastic deformation and flow stress of polymeric materials and certain metallic alloys. This sensitivity to hydrostatic stress can also influence the fracture toughness of ductile materials, which fail by void growth and coalescence. These materials typically contain a non-uniform distribution of voids of varying size-scales and void shapes. In this work, the effects of void shape and microvoid interaction in pressure-sensitive materials are examined via a two-prong approach: (i) an axisymmetric unit-cell containing a single ellipsoidal void and (ii) a plane-strain unit-cell consisting of a single large void and a population of discrete microvoids. The representative material volume in both cases is subjected to physical stress states similar to highly stressed regions ahead of a crack. Results show that oblate voids and microvoid cavitation can severely reduce the critical stress of the material. These effects can be compounded under high levels of pressure-sensitivity. In some cases, the critical stress responsible for rapid void growth is reduced to levels comparable to the yield strength of the material. The contribution of void shape and pressure-sensitivity to the thermal- and moisture-induced voiding phenomenon in IC packages is also discussed.     

Simulation of convectively coupled waves using WRF: a framework for assessing the effects of mesoscales on synoptic scales

The atmospheric variability in the tropics is primarily driven by convective heating. Observations revealed that convection in the tropics is organized into a hierarchy of multiscale convective systems ranging from the individual cloud cells to planetary scale disturbances that are nested within each other like Russian dolls. Current global climate models simulate very poorly these convectively coupled waves due in part to inadequate treatment of organized convection by the underlying cumulus parameterizations. Here, we present idealized simulations of convectively coupled equatorial waves (CCWs) using the weather research and forecast model in a horizontally limited domain consisting of a 4,500 km-wide square centered at the equator at moderate horizontal resolution of 10 km. We attempted and compared various configuration options, including switching on and off the cumulus parameterization (CP) and nesting a fine resolution 3.33 km domain, a 2,000 km-wide square, in the middle of the domain. It turns out that the results without a CP are much superior than those using a CP. While the cases without a CP resulted in a coherent eastward propagating CCW, which has many common features with observed convectively coupled Kelvin waves, the cumulus parameterization tends to destroy both the coherence of the propagating waves, even in the case with a nested domain, and reduces dramatically the variability. A primary demonstration on how such results could be used to show evidence of energy exchange, through momentum transport, between small-scale circulation due to mesoscale convection and the propagating synoptic scale wave will be reported is also presented.