Prediction of failure in bending of an aluminium sheet alloy

This work is dedicated to numerical prediction of the bending of thin aluminium alloy sheets, with a focus on the material parameter identification and the prediction of rupture with or without pre-strains in tension prior to bending. The experimental database consists of i) mechanical tests at room temperature, such as tension and simple shear, performed at several orientations to the rolling direction and biaxial tension ii) air bending tests of rectangular samples after (or not) pre-straining in tension. The mechanical model is composed of the Yld2004-18p anisotropic yield criterion (Barlat et al. [3]) associated with a mixed hardening rule. The material parameters (altogether 21) are optimized with an inverse approach, in order to minimize the gap between experimental data and model predictions. Then, the Hosford-Coulomb rupture criterion is used in an uncoupled way, and the parameters are determined from tensile tests, both uniaxial and biaxial, with data up to rupture. In a second step, numerical simulations of the bending tests are performed, either on material in its original state or after pre-straining in tension, with pre-strain magnitudes increasing from 0.19 up to 0.3. The comparisons are performed on different outputs: load evolution, strain field and prediction of the rupture. A very good correlation is obtained over all the tests, in the identification step as well as in the validation one. Moreover, the fracture criterion proves to be successful whatever the amount of pre-strain may be. A convincing representation of the mechanical behavior at room temperature for an aluminium alloy is thus obtained.

Families of asymptotic tensile crack tip stress fields in elastic-perfectly plastic single crystals

In this work, two families of asymptotic near-tip stress fields are constructed in an elastic-ideally plastic FCC single crystal under mode I plane strain conditions. A crack is taken to lie on the (010) plane and its front is aligned along the [(1) over bar 01] direction. Finite element analysis is first used to systematically examine the stress distributions corresponding to different constraint levels. The general framework developed by Rice (Mech Mater 6:317-335, 1987) and Drugan (J Mech Phys Solids 49:2155-2176, 2001) is then adopted to generate low triaxiality solutions by introducing an elastic sector near the crack tip. The two families of stress fields are parameterized by the normalized opening stress (tau(A)(22)/tau(o)) prevailing in the plastic sector in front of the tip and by the coordinates of a point where elastic unloading commences in stress space. It is found that the angular stress variations obtained from the analytical solutions show good agreement with finite element analysis.

A comparison of stability and bifurcation criteria for inflated spherical elastic shells

The nonlinear stability analysis introduced by Chen and Haughton [1] is employed to study the full nonlinear stability of the non-homogeneous spherically symmetric deformation of an elastic thick-walled sphere. The shell is composed of an arbitrary homogeneous, incompressible elastic material. The stability criterion ultimately requires the solution of a third-order nonlinear ordinary differential equation. Numerical calculations performed for a wide variety of well-known incompressible materials are then compared with existing bifurcation results and are found to be identical. Further analysis and comparison between stability and bifurcation are conducted for the case of thin shells and we prove by direct calculation that the two criteria are identical for all modes and all materials.

Wave propagation in elastic solids undergoing damage and growth process

Modeling and analysis of wave propagation in elastic solids undergoing damage and growth process are reported in this paper. Two types of diagnostic problems, (1) the propagation of waves in the presence of a slow growth process and (2) the propagation of waves in the presence of a fast growth process, are considered. The proposed model employs a slow and a fast time scale and a homogenization technique in the wavelength scale. A detailed analysis of wave dispersion is carried out. A spectral analysis reveals certain low-frequency bands, where the interaction between the wave and the growth process produces acoustic metamaterial-like behavior. Various practical issues in designing an efficient method of acousto-ultrasonic wave based diagnostics of the growth process are discussed. Diagnostics of isotropic damage in a ductile or quasi-brittle solid by using a micro-second pulsating signal is considered for computer simulations, which is to illustrate the practical application of the proposed modeling and analysis. The simulated results explain how an estimate of signal spreading can be effectively employed to detect the presence of a steady-state damage or the saturation of a process.

Stress-induced phase transformation and pseudo-elastic/pseudo-plastic recovery in intermetallic Ni-Al nanowires

Extensive molecular dynamics (MD) simulations have been performed in a B2-NiAl nanowire using an embedded atom method (EAM) potential. We show a stress induced B2 -> body-centered-tetragonal (BCT) phase transformation and a novel temperature and cross-section dependent pseudo-elastic/pseudo-plastic recovery from such an unstable BCT phase with a recoverable strain of similar to 30% as compared to 5-8% in polycrystalline materials. Such a temperature and cross-section dependent pseudo-elastic/pseudo-plastic strain recovery can be useful in various interesting applications of shape memory and strain sensing in nanoscale devices. Effects of size, temperature, and strain rate on the structural and mechanical properties have also been analyzed in detail. For a given size of the nanowire the yield stress of both the B2 and the BCT phases is found to decrease with increasing temperature, whereas for a given temperature and strain rate the yield stress of both the B2 and the BCT phase is found to increase with increase in the cross-sectional dimensions of the nanowire. A constant elastic modulus of similar to 80 GPa of the B2 phase is observed in the temperature range of 200-500 K for nanowires of cross-sectional dimensions in the range of 17.22-28.712 angstrom, whereas the elastic modulus of the BCT phase shows a decreasing trend with an increase in the temperature.

Strain gradient theory with couple stress for crystalline solids

A new phenomenological strain gradient theory for crystalline solid is proposed. It fits within the framework of general couple stress theory and involves a single material length scale Ics. In the present theory three rotational degrees of freedom omega (i) are introduced, which denote part of the material angular displacement theta (i) and are induced accompanying the plastic deformation. omega (i) has no direct dependence upon u(i) while theta = (1 /2) curl u. The strain energy density omega is assumed to consist of two parts: one is a function of the strain tensor epsilon (ij) and the curvature tensor chi (ij), where chi (ij) = omega (i,j); the other is a function of the relative rotation tensor alpha (ij). alpha (ij) = e(ijk) (omega (k) - theta (k)) plays the role of elastic rotation reason The anti-symmetric part of Cauchy stress tau (ij) is only the function of alpha (ij) and alpha (ij) has no effect on the symmetric part of Cauchy stress sigma (ij) and the couple stress m(ij). A minimum potential principle is developed for the strain gradient deformation theory. In the limit of vanishing l(cs), it reduces to the conventional counterparts: J(2) deformation theory. Equilibrium equations, constitutive relations and boundary conditions are given in detail. For simplicity, the elastic relation between the anti-symmetric part of Cauchy stress tau (ij), and alpha (ij) is established and only one elastic constant exists between the two tensors. Combining the same hardening law as that used in previously by other groups, the present theory is used to investigate two typical examples, i.e., thin metallic wire torsion and ultra-thin metallic beam bend, the analytical results agree well with the experiment results. While considering the, stretching gradient, a new hardening law is presented and used to analyze the two typical problems. The flow theory version of the present theory is also given.

Effective elastic moduli of two-dimensional solids with multiple cracks

An accurate method which directly accounts for the interactions between different microcracks is used for analyzing the elastic problem of multiple cracks solids. The effective elastic moduli for randomly oriented cracks and parallel cracks are evaluated for the representative volume element (RVE) with microcracks in infinite media. The numerical results are compared with those from various micromechanics models and experimental data. These results show that the present method is simple and provides a direct and efficient approach to dealing with elastic solids containing multiple cracks.

Nonlinear Analysis of Oscillatory Indentation in Elastic and Viscoelastic Solids

Determining the mechanical properties at micro- and nanometer length scales using nanoindentation or atomic force microscopy is important to many areas of science and engineering. Here we establish equations for obtaining storage and loss modulus from oscillatory indentations by performing a nonlinear analysis of conical and spherical indentation in elastic and viscoelastic solids. We show that, when the conical indenter is driven by a sinusoidal force, the square of displacement is a sinusoidal function of time, not the displacement itself, which is commonly assumed. Similar conclusions hold for spherical indentations. Well-known difficulties associated with measuring contact area and correcting thermal drift may be circumvented using the newly derived equations. These results may help improve methods of using oscillatory indentation for determining elastic and viscoelastic properties of solids.

Effective elastic moduli of inhomogeneous solids by embedded cell model

An embedded cell model is presented to obtain the effective elastic moduli for three-dimensional two-phase composites which is an exact analytic formula without any simplified approximation and can be expressed in an explicit form. For the different cells such as spherical inclusions and cracks surrounded by sphere and oblate ellipsoidal matrix, the effective elastic moduli are evaluated and the results are compared with those from various micromechanics models. These results show that the present model is direct, simple and efficient to deal with three-dimensional tyro-phase composites.

The inhomogeneity of plastic-deformation in ductile single-crystals

A discrete slip model which characterizes the inhomogeneity of material properties in ductile single crystals is proposed in this paper. Based on this model rate-dependent finite element investigations are carried out which consider the finite deformation, finite rotation, latent hardening effect and elastic anisotropy. The calculation clearly exhibits the process from microscopic inhomogeneous and localized deformation to necking and the formation of LSBS and reveals several important features of shear localization. For example, the inhomogeneous deformation is influenced by the imperfections and initial non-uniformities of material properties. The inhomogeneous deformation may either induce necking which results in the lattice rotation and leads to geometrical softening, which in turn promotes the formation of CSBS, or induces heavily localized deformation. The microscopic localized deformation eventually develops into the LSBS and results in a failure. These results are in close agreement with experiment. Our calculations also find that the slip lines on the specimen's surface at necking become curved and also find that if the necking occurs before the formation of LSBS, this band must be misoriented from the operative slip systems. In this case, the formation of LSBS must involve non-crystallographic effects. These can also be indirectly confirmed by experiment. All these suggest that our present discrete slip model offers a correct description of the inhomogeneous deformation characterization in ductile crystals.

Kinking of an interface crack between 2 dissimilar anisotropic elastic solids

A detailed analysis of kinking of an interface crack between two dissimilar anisotropic elastic solids is presented in this paper. The branched crack is considered as a distributed dislocation. A set of the singular integral equations for the distribution function of the dislocation density is developed. Explicit formulas of the stress intensity factors and the energy release rates for the branched crack are given for orthotropic bimaterials and misoriented orthotropic bicrystals. The role of the stress parallel to the interface, sigma0 is taken into account in these formulas. The interface crack can advance either by continued extension along the interface or by kinking out of the interface into one of the adjoining materials. This competition depends on the ratio of the energy release rates for interface cracking and for kinking out of the interface and the ratio of interface toughness to substrate toughness. Throughout the paper, the influences of the inplane stress sigma0 on the stress intensity factors and the energy release rates for the branched crack, which can significantly alter the conditions for interface cracking, are emphasized.

Higher-order analysis of crack tip fields in elastic power-law hardening materials

A HIGHER-ORDER asymptotic analysis of a stationary crack in an elastic power-law hardening material has been carried out for plane strain, Mode 1. The extent to which elasticity affects the near-tip fields is determined by the strain hardening exponent n. Five terms in the asymptotic series for the stresses have been derived for n = 3. However, only three amplitudes can be independently prescribed. These are K1, K2 and K5 corresponding to amplitudes of the first-, second- and fifth-order terms. Four terms in the asymptotic series have been obtained for n = 5, 7 and 10; in these cases, the independent amplitudes are K1, K2 and K4. It is found that appropriate choices of K2 and K4 can reproduce near-tip fields representative of a broad range of crack tip constraints in moderate and low hardening materials. Indeed, fields characterized by distinctly different stress triaxiality levels (established by finite element analysis) have been matched by the asymptotic series. The zone of dominance of the asymptotic series extends over distances of about 10 crack openings ahead of the crack tip encompassing length scales that are microstructurally significant. Furthermore, the higher-order terms collectively describe a spatially uniform hydrostatic stress field (of adjustable magnitude) ahead of the crack. Our results lend support to a suggestion that J and a measure of near-tip stress triaxiality can describe the full range of near-tip states.

Micromechanics of elastic solids weakened by a large number of cracks

This paper presents a micromechanics analysis of the elastic solids weakened by a large number of microcracks in a plane problem. A new cell model is proposed. Each cell is an ellipse subregion and contains a microcrack. The effective moduli and the stress intensity factors for an ellipse cell are obtained. The analytic closed formulas of concentration factor tensor for an isotropic matrix containing an anisotropic inclusion are derived. Based on a self-consistent method, the effective elastic moduli of the solids weakened by randomly oriented microcracks are obtained.