An elastoplastic damage model for metal matrix composites considering progressive imperfect interface under transverse loading

An elastoplastic damage model considering progressive imperfect interface is proposed to predict the effective elastoplastic behavior and multi-level damage progression in fiber-reinforced metal matrix composites (FRMMCs) under transverse loading. The modified Eshelby’s tensor for a cylindrical inclusion with slightly weakened interface is adopted to model fibers having mild or severe imperfect interfaces [Lee, H.K., Pyo, S.H., 2009. A 3D-damage model for fiber-reinforced brittle composites with microcracks and imperfect interfaces. J. Eng. Mech. ASCE. doi:10.1061/(ASCE)EM.1943-7889.0000039]. An elastoplastic model is derived micromechanically on the basis of the ensemble-volume averaging procedure and the first-order effects of eigenstrains. A multi-level damage model [Lee, H.K., Pyo, S.H., 2008a. Multi-level modeling of effective elastic behavior and progressive weakened interface in particulate composites. Compos. Sci. Technol. 68, 387–397] in accordance with the Weibull’s probabilistic function is then incorporated into the elastoplastic multi-level damage model to describe the sequential, progressive imperfect interface in the composites. Numerical examples corresponding to uniaxial and biaxial transverse tensile loadings are solved to illustrate the potential of the proposed micromechanical framework. A series of parametric analysis are carried out to investigate the influence of model parameters on the progression of imperfect interface in the composites. Furthermore, a comparison between the present prediction and experimental data in the literature is made to assess the capability of the proposed micromechanical framework.     

Modeling of damage in unidirectional ceramic matrix composites and multi-scale experimental validation on third generation SiC/SiC minicomposites

The purpose of this paper is to experimentally validate a 1D probabilistic model of damage evolution in unidirectional SiC/SiC composites. The key point of this approach lies in the identification and validation at both local and macroscopic scales. Thus, in addition to macroscopic tensile tests, the evolution of microscopic damage mechanisms – in the form of matrix cracks and fiber breaks – is experimentally analyzed and quantified through in-situ scanning electron microscope and computed tomography tensile tests. A complete model, including both matrix cracking and fiber breaking, is proposed on the basis of existing modeling tools separately addressing these mechanisms. It is based on matrix and fiber failure probability laws and a stress redistribution assumption in the vicinity of matrix cracks or fiber breaks. The identification of interfacial parameters is conducted to fit the experimental characterization, and shows that conventional assumptions of 1D probabilistic models can adequately describe matrix cracking at both macro- and microscopic scales. However, it is necessary to enrich them to get a proper prediction of ultimate failure and fiber break density for Hi-Nicalon type S fiber-reinforced SiC/SiC minicomposites.     

Global load-sharing model for unidirectional hybrid fibre-reinforced composites

A promising strategy to increase the tensile failure strain of carbon fibre-reinforced composites is to hybridise carbon fibres with other, higher-elongation fibres. The resulting increase in failure strain is known as the hybrid effect. In the present article, a global load-sharing model for hybrid composites is developed and used to carry out a parametric study for carbon/glass hybrids. Hybrid effects of up to 15% increase in failure strain are predicted, corresponding reasonably well to literature data. Scatter in the carbon fibre strength is shown to be crucial for the hybrid effect, while the scatter in glass fibre strength is much less important. In contrast to reports in earlier literature, the ratio of failure strains of the two fibres has only a small influence on the hybrid effect. The results provide guidelines for designing optimal hybrid composites.     

Compressive strength of fibre composites with random fibre waviness

The compressive strength of unidirectional long fibre composites is predicted for plastic microbuckling from a random two-dimensional distribution of fibre waviness. The effect of the physical size of waviness is addressed by using couple stress theory, with the fibre bending resistance scaling with the fibre diameter d. The predicted statistical distribution of compressive strength is found using a Monte Carlo method. An ensemble of fibre waviness profiles is generated from an assumed spectral density of waviness and the compressive strength for each such realisation is calculated directly by the finite element method. The average predicted strength agrees reasonably with practical values, confirming the hypothesis that microbuckles can be initiated by fibre misalignment. It is found that the probability distribution of strength is well matched by a Weibull fit, and the dependence of the Weibull parameters upon the spectral density of waviness is determined. For the practical range of fibre distributions considered, it is concluded that the strength depends mainly upon the root mean square amplitude of fibre misalignment, with the shape of the power spectral density function playing only a minor role. An engineering model for predicting the compressive strength is proposed, akin to weakest link theory for materials containing flaws. A specimen containing randomly distributed waviness is examined to locate regions of high-fibre misalignment. The strength of each of these weak regions is estimated from a look-up table derived from calculations with idealised circular or elliptical patches of waviness. The strength of the composite is given by the failure stress associated with the weakest such patch. For random distributions of waviness, the predictions using this engineering approach are in good agreement with the direct calculations of strength using the finite element method.     

Debonding failure and size effects in micro-reinforced composites

Failure in micro-reinforced composites is investigated numerically using the strain-gradient plasticity theory of Gudmundson [Gudmundson, P., 2004. A unified treatment of strain gradient plasticity. Journal of the Mechanics and Physics of Solids 52 (6) 1379–1406] in a plane strain visco-plastic formulation. Bi-axially loaded unit cells are used and failure is modeled using a cohesive zone at the reinforcement interface. During debonding a sudden stress drop in the overall average stress–strain response is observed. Adaptive higher-order boundary conditions are imposed at the reinforcement interface for realistically modeling the restrictions on moving dislocations as debonding occurs. It is found that the influence of the imposed higher-order boundary conditions at the interface is minor. If strain-gradient effects are accounted for a void with a smooth shape develops at the reinforcement interface while a smaller void having a sharp tip nucleates if strain-gradient effects are excluded. Using orthogonalization of the plastic strain gradient with three corresponding material length scales it is found that, the first length scale dominates the evaluated overall average stress–strain response, the second one only has a small effect and the third one has an intermediate effect. Finally, studies of reinforcement having elliptical cross-sections show rather significant gradients of stress which is not seen for the corresponding circular cross-sections. Also, an increased drop in the overall load carrying capacity is observed for cross-sections elongated perpendicular to the principal tensile direction compared to the corresponding circular cross-sections.     

Strain gradient plasticity effects in whisker-reinforced metals

A metal reinforced by fibers in the micron range is studied using the strain gradient plasticity theory of Fleck and Hutchinson (J. Mech. Phys. Solids 49 (2001) 2245). Cell-model analyses are used to study the influence of the material length parameters numerically, for both a single parameter version and the multiparameter theory, and significant differences between the predictions of the two models are reported. It is shown that modeling fiber elasticity is important when using the present theories. A significant stiffening effect when compared to conventional models is predicted, which is a result of a significant decrease in the level of plastic strain. Moreover, it is shown that the relative stiffening effect increases with fiber volume fraction. The higher-order nature of the theories allows for different higher-order boundary conditions at the fiber-matrix interface, and these boundary conditions are found to be of importance. Furthermore, the influence of the material length parameters on the stresses along the interface between the fiber and the matrix material is discussed, as well as the stresses within the elastic fiber which are of importance for fiber breakage.     

Experimental investigation of initial and subsequent yield surfaces for laminated metal matrix composites

The aim of this work is to construct yield surfaces to describe initial yielding and characterize hardening behavior of a highly anisotropic material. A methodology for constructing yield surfaces for isotropic materials using axial–torsion loading is extended to highly anisotropic materials. The technique uses a sensitive definition of yielding based on permanent strain rather than offset strain, and enables multiple yield points and multiple yield surfaces to be conducted on a single specimen. A target value of 20 × 10⁻⁶ is used for Al₂O₃ fiber reinforced aluminum laminates having a fiber volume fraction of 0.55. Sixteen radial probes are used to define the yield locus in the axial–shear stress plane. Initial yield surfaces for [0₄], [90₄], and [0/90]₂ fibrous aluminum laminates are well described by ellipses in the axial–shear stress plane having aspect ratios of 10, 2.5, and 3.3, respectively. For reference, the aspect ratio of the Mises ellipse for an isotropic material is 1.73. Initial yield surfaces do not have a tension–compression asymmetry. Four overload profiles (plus, ex, hourglass, and zee) are applied to characterize hardening of a [0/90]₂ laminate by constructing 30 subsequent yield surfaces. Parameters to describe the center and axes of an ellipse are regressed to the yield points. The results clearly indicate that kinematic hardening dominates so that material state evolution can be described by tracking the center of the yield locus. For a nonproportional overload of (σ, τ) = (500, 70) MPa, the center of the yield locus translated to (σ, τ) = (430, 37) MPa and the ellipse major axis was only 110 MPa.     

A two-dimensional finite element method for simulating the constitutive response and microstructure of polycrystals during high temperature plastic deformation

We describe a finite element method designed to model the mechanisms that cause superplastic deformation. Our computations account for grain boundary sliding, grain boundary diffusion, grain boundary migration, and surface diffusion, as well as thermally activated dislocation creep within the grains themselves. Front tracking and adaptive mesh generation are used to follow changes in the grain structure. The method is used to solve representative boundary value problems to illustrate its capabilities.     

Electroelastic field concentrations ahead of electrodes in multilayer piezoelectric actuators: experiment and finite element simulation

The effects of applied voltage on the electroelastic field concentrations ahead of electrodes in multilayer piezoelectric actuators were examined in a combined experimental and numerical investigation. Experiments were performed to measure the strain near internal and surface electrodes at various electrical loading conditions. The finite element method was also used to solve the coupled electro-elastic boundary value problem. The strain, stress and electric displacement concentrations were calculated and a non-linear behavior induced by localized polarization switching was discussed. A comparison of strain concentration was made between experiment and simulation.     

Effects of spatial grain orientation distribution and initial surface topography on sheet metal necking

The finite element method is used to numerically simulate localized necking in AA6111-T4 under stretching. The measured EBSD data (grain orientations and their spatial distributions) are directly incorporated into the finite element model and the constitutive response at an integration point is described by the single crystal plasticity theory. We assume that localized necking is associated with surface instability, the onset of unstable growth in surface roughening. It is demonstrated that such a surface instability/necking is the natural outcome of the present approach, and the artificial initial imperfection necessitated by the macroscopic M–K approach [Marciniak and Kuczynski (1967). Int. J. Mech. Sci. 9, 609–620] is not relevant in the present analysis. The effects of spatial orientation distribution, material strain rate sensitivity, texture evolution, and initial surface topography on necking are discussed. It is found that localized necking depends strongly on both the initial texture and its spatial orientation distribution. It is also demonstrated that the initial surface topography has only a small influence on necking.     

On the fracture toughness of ferroelastic materials

The toughness enhancement due to domain switching near a steadily growing crack in a ferroelastic material is analyzed. The constitutive response of the material is taken to be characteristic of a polycrystalline sample assembled from randomly oriented tetragonal single crystal grains. The constitutive law accounts for the strain saturation, asymmetry in tension versus compression, Bauschinger effects, reverse switching, and strain reorientation that can occur in these materials due to the non-proportional loading that arises near a propagating crack. Crack growth is assumed to proceed at a critical level of the crack tip energy release rate. Detailed finite element calculations are carried out to determine the stress and strain fields near the growing tip, and the ratio of the far field applied energy release rate to the crack tip energy release rate. The results of the finite element calculations are then compared to analytical models that assume the linear isotropic K-field solution holds for either the near tip stress or strain field. Ultimately, the model is able to account for the experimentally observed toughness enhancement in ferroelastic ceramics.     

FE analysis of stress and strain fields in finite random composite bodies

In this paper the mechanical behaviour of finite random heterogeneous bodies is considered. The analysis of non-local interactions between heterogeneities in microscopically heterogeneous materials is necessary when the spatial variation of the load or the dimensions of the body, relative to the scale of the microstructure, cannot be ignored. Microstructures can be periodic but generically they are random. In the first case, an exact calculation can be performed but in the second case recourse has to be made either to simulation or to some scheme of approximation. One such scheme is based on a stochastic variational principle. The novelty of the present work is that a stochastic variational principle is projected directly onto a finite-element basis so that all subsequent analysis is performed within a finite-element framework. The proposed formulation provides expressions for the local stress and strain fields in any realization of the medium, from which expressions for statistically-averaged quantities can be derived. Then an approximation of Hashin-Shtrikman type is developed, which generates a FE-based numerical procedure able to take account of interactions between random inclusions and boundary layer effects in finite composite structures. Finally, two examples are presented, namely a cylinder with square cross-section subjected to mixed boundary conditions of different types on different faces and a rectangular body containing a centre crack. The results show that in the vicinity of the boundary or close to the crack tip, the strain and the stress in the matrix and in the inclusions differ considerably from those obtained by the formal application of conventional homogenization.     

A non-associated constitutive model with mixed iso-kinematic hardening for finite element simulation of sheet metal forming

In this paper an anisotropic material model based on non-associated flow rule and mixed isotropic–kinematic hardening was developed and implemented into a user-defined material (UMAT) subroutine for the commercial finite element code ABAQUS. Both yield function and plastic potential were defined in the form of Hill’s [Hill, R., 1948. A theory of the yielding and plastic flow of anisotropic metals. Proc. R. Soc. Lond. A 193, 281–297] quadratic anisotropic function, where the coefficients for the yield function were determined from the yield stresses in different material orientations, and those of the plastic potential were determined from the r-values in different directions. Isotropic hardening follows a nonlinear behavior, generally in the power law form for most grades of steel and the exponential law form for aluminum alloys. Also, a kinematic hardening law was implemented to account for cyclic loading effects. The evolution of the backstress tensor was modeled based on the nonlinear kinematic hardening theory (Armstrong–Frederick formulation). Computational plasticity equations were then formulated by using a return-mapping algorithm to integrate the stress over each time increment. Either explicit or implicit time integration schemes can be used for this model. Finally, the implemented material model was utilized to simulate two sheet metal forming processes: the cup drawing of AA2090-T3, and the springback of the channel drawing of two sheet materials (DP600 and AA6022-T43). Experimental cyclic shear tests were carried out in order to determine the cyclic stress–strain behavior and the Bauschinger ratio. The in-plane anisotropy (r-value and yield stress directionalities) of these sheet materials was also compared with the results of numerical simulations using the non-associated model. These results showed that this non-associated, mixed hardening model significantly improves the prediction of earing in the cup drawing process and the prediction of springback in the sidewall of drawn channel sections, even when a simple quadratic constitutive model is used.     

On oblique contact of creeping solids

Oblique indentation of power-law creeping solids by a rigid die is analysed in three dimensions with perfectly plastic behaviour emerging as an asymptotic case. Indenter profiles are prescribed to be axisymmetric for simplicity but not by necessity. Invariance and generality is aimed at, as the problem is governed by only four essential parameters, i.e. the die profile, p, the indentation angle, γ, the power-law exponent, n, and the coefficient of friction, μ. The solution strategy is based on a self-similar transformation resulting in a reduced problem corresponding to flat die indentation of complete contact. The reduced auxiliary problem, being independent of loading, history and time, was solved by a three-dimensional finite element analysis characterized by high accuracy. Subsequently, cumulative superposition was used to resolve the original problem and global and invariant relations between force, depth and contact area were determined. Detailed results are given for the location and shape of the contact region and stick/slip contours as well as for local states of surface stresses and deformation at flat and spherical indenters. Due to the asymmetry prevailing, it was found that in the spherical case, contact contours proved to be oval and shifted, although with normal and tangential forces only weakly coupled. Finite friction as compared to full adhesion proved to have only a minor effect on global relations. The framework laid down may be applied to the contact of structural assemblies subjected especially to elevated temperatures and also to various issues such as compaction of powder aggregates, flattening of rough surfaces and plastic impact.     

Dynamic plasticity and fracture in high density polycrystals: constitutive modeling and numerical simulation

Presented is a constitutive framework for modeling the dynamic response of polycrystalline microstructures, posed in a thermodynamically consistent manner and accounting for finite deformation, strain rate dependence of flow stress, thermal softening, thermal expansion, heat conduction, and thermoelastic coupling. Assumptions of linear and square-root dependencies, respectively, of the stored energy and flow stresses upon the total dislocation density enable calculation of the time-dependent fraction of plastic work converted to heat energy. Fracture at grain boundary interfaces is represented explicitly by cohesive zone models. Dynamic finite element simulations demonstrate the influences of interfacial separation, random crystallographic orientation, and grain morphology on the high-rate tensile response of a realistic two-phase material system consisting of comparatively brittle pure tungsten (W) grains embedded in a more ductile matrix of tungsten-nickel iron (W-Ni-Fe) alloy. Aspects associated with constitutive modeling of damage and failure in the homogenized material system are discussed in light of the computational results.     

The effect of time-dependent twinning on low temperature (<0.25 ∗ Tm) creep of an alpha-titanium alloy

The low-temperature (less than one-fourth of the melting temperature) creep deformation behavior of hexagonally close-packed (HCP) α-Ti–1.6 wt.% V was investigated. Creep tests were performed at various temperatures between room temperature and 205 °C at 95% of the respective yield stress at the different temperatures. The creep strain rate was found to increase with increasing temperature. Scanning and transmission electron microscopy revealed that slip and unusually slow twin growth, or time-dependent twinning, are active deformation mechanisms for the entire temperature range of this investigation. The activation energy for creep of this alloy was calculated to identify the rate-controlling deformation mechanism, and was found to increase with increasing creep strain. At low strain, the activation energy for creep was found to be close to the previously calculated activation energy for slip. At high strain, the calculated activation energy indicates that both slip and twinning are significant deformation mechanisms. The appearance of twinning at high strains is explained by a model for twin nucleation by dislocation pileups.