On the wrinkling of a pre-stressed annular thin film in tension

Asymptotic properties of the neutral stability curves for a linear boundary eigenvalue problem which models the wrinkling instability of an annular thin film in tension are considered. The film is subjected to imposed radial displacement fields on its inner and outer boundaries and, when these loads are sufficiently large, the film is susceptible to wrinkling. The critical values at which this onset occurs are dictated by the solution of a fourth-order ordinary differential eigensystem whose eigenvalue λ is a function of μ(?1), a quantity inversely proportional to the non-dimensional bending stiffness of the film, and n, the number of half-waves of the wrinkling pattern that sets in around the annular domain. Previously, Coman and Haughton [2006. Localised wrinkling instabilities in radially stretched annular thin films. Acta Mech. 185, 179-200] employed the compound matrix method together with a WKB technique to characterise the form of λ(μ,n) which essentially is related to a turning point in a reduced eigenproblem. The asymptotic analysis conducted therein pertained to the case when this turning point was not too close to the inner edge of the annulus. However, in the thin film limit μ→∞, the wrinkling load and the preferred instability mode are given by a modified eigenvalue problem that involves a turning point asymptotically close to the inner rim. Here WKB and boundary-layer asymptotic methods are used to examine these issues and comparisons with direct numerical simulations made.     

Path-dependence of the forming limit stresses in a sheet metal

The effect of changing strain paths on the forming limit stresses of sheet metals is investigated using the Marciniak–Kuczyński model and a phenomenological plasticity model with non-normality effects [Kuroda, M., Tvergaard, V., 2001. A phenomenological plasticity model with non-normality effects representing observations in crystal plasticity. J. Mech. Phys. Solids 49, 1239–1263]. Forming limits are simulated for linear stress paths and two types of combined loading: a combined loading consisting of two linear stress paths in which unloading is included between the first and second loadings (combined loading A), and combined loading in which the strain path is abruptly changed without unloading (combined loading B). The forming limit stresses calculated for combined loading A agree well with those calculated for the linear stress paths, while the forming limit curves in strain space depend strongly on the strain paths. The forming limit stresses calculated for the combined loading B do not, however, coincide with those calculated for the linear stress paths. The strain-path dependence of the forming limit stress is discussed in detail by observing the strain localization process.     

Anisotropic finite elastoplasticity with nonlinear kinematic and isotropic hardening and application to sheet metal forming

The paper discusses the derivation and the numerical implementation of a finite strain material model for plastic anisotropy and nonlinear kinematic and isotropic hardening. The model is derived from a thermodynamic framework and is based on the multiplicative split of the deformation gradient in the context of hyperelasticity. The kinematic hardening component represents a continuum extension of the classical rheological model of Armstrong–Frederick kinematic hardening. Introducing the so-called structure tensors as additional tensor-valued arguments, plastic anisotropy can be modelled by representing the yield surface and the plastic flow rule as functions of the structure tensors. The evolution equations are integrated by a new form of the exponential map that preserves plastic incompressibility and uses the spectral decomposition to evaluate the exponential tensor functions in closed form. Finally, the applicability of the model is demonstrated by means of simulations of several deep drawing processes and comparisons with experiments.     

Surface wrinkling of mucosa induced by volumetric growth: Theory, simulation and experiment

Mechanics of living tissues focusing on the relationships between growth, morphology and function is not only of theoretical interest but can also be useful for diagnosis of certain diseases. In this paper, we model the surface wrinkling morphology of mucosa, the moist tissue that commonly lines organs and cavities throughout the body, induced by either physiological or pathological volumetric growth. A theoretical framework of finite deformation is adopted to analyze the deformation of a cylindrical cavity covered by mucosal and submucosal layers. It is shown that compressive residual stresses induced by the confined growth of mucosa can destabilize the tissue into various surface wrinkling patterns. A linear stability analysis of the critical condition and characteristic buckling patterns indicates that the wrinkling mode is sensitive to the thicknesses of the mucosal and submucosal layers, as well as the properties of the tissues. The thinner the mucosal layer and the lower its elastic modulus, the shorter the buckling wavelength. A series of finite element simulations are performed to validate the theoretical predictions and to study local wrinkling or non-uniform patterns associated with inhomogeneous growth. Our postbuckling analysis shows that the surface pattern may evolve towards a period-doubling morphology due to continuous growth of mucosa or submucosa beyond the critical state. Finally, the theoretical predictions and numerical simulations are compared to experimental observations.     

Deformation localization and shear band fracture in strong anisotropy sheet tension

The tensile deformation localization and the shear band fracture behaviors of sheet metals with strong anisotropy are numerically simulated by using Updating Lagrange finite element method, Quasi-flow plastic constitutive theory^[1] and B-L planar anisotropy yield criterion^[2]. Simulated results are compared with experimental ones. Very good consistence is obtained between numerical and experimental results. The relationship between the anisotropy coefficientR and the shear band angle θ is found. The project supported by the National Natural Science Foundation of China and the Excellent Youth Teacher Foundation of the State Education Commission of China     

Marciniak–Kuczynski type modelling of the effect of Through-Thickness Shear on the forming limits of sheet metal

The Marciniak–Kuczynski (MK) forming limit model is extended in order to predict localized necking in sheet metal forming operations in which Through-Thickness Shear (TTS), also known as out-of-plane shear, occurs. An example of such a forming operation is Single Point Incremental Forming. The Forming Limit Diagram (FLD) of a purely plastic, isotropic hardening material with von Mises yield locus is discussed, for monotonic deformation paths that include TTS. If TTS is present in the plane containing the critical groove direction in the MK model, it is seen that formability is increased for all in-plane strain modes, except equibiaxial stretching. The increase in formability due to TTS is explained through a detailed study of some selected deformation modes. The underlying mechanism is a change of the stress mode in the groove that results in a delay of the onset of localized necking.     

Failure in notched tension bars due to high-temperature creep: Interaction between nucleation controlled cavity growth and continuum cavity growth

For axi-symmetrically notched tension bars [Dyson, B.F., Loveday, M.S., 1981, Creep Fracture in Nimonic 80A under Tri-axial Tensile Stressing, In: Ponter A.R.S., Hayhurst, D.R. (Eds.), Creep in Structures, Springer-Verlag, Berlin, pp. 406–421] show two types of damage propagation are shown: for low stress, failure propagates from the outside notch surface to the centre-line; and for high stress, failure propagates from the centre-line to the outside notch surface. The objectives of the paper are to: identify the physics of the processes controlling global failure modes; and, describe the global behaviour using physics-based constitutive equations.Two sets of constitutive equations are used to model the softening which takes place in tertiary creep of Nimonic 80A at 750 °C. Softening by multiplication of mobile dislocations is firstly combined, for low stress, with softening due to nucleation controlled creep constrained cavity growth; and secondly combined, for high stress, with softening due to continuum void growth. The Continuum Damage Mechanics, CDM, Finite Element Solver DAMAGE XX has been used to study notch creep fracture. Low stress notch behaviour is accurately predicted provided that the constitutive equations take account of the effect of stress level on creep ductility. High stress notch behaviour is accurately predicted from a normalized inverse cavity spacing d/2? = 6, and an initial normalized cavity radius r_hi/? = 3.16 × 10^?3, where 2? is the cavity spacing, and d is the grain size; however, the constants in the strain rate equation required recalibration against high stress notch data. A void nucleation mechanism is postulated for high stress behaviour which involves decohesion where slip bands intersect second phase grain boundary particles. Both equation sets accurately predict experimentally observed global failure modes.     

Initiation and propagation of localized deformation in elasto-plastic strips under uniaxial tension

This paper deals with the evolution of inhomogeneous deformation in shape memory alloy strips and mild steel strips under uniaxial tension. New experiments on NiTi strips, which initially are in an austenitic phase, show that at a critical stress level martensite nucleates in sharp bands inclined at 55° to the axis of loading. Under prescribed end displacement martensite subsequently spreads either by steady-state propagation of inclined transition fronts or via a criss-cross pattern of finger-like features. Similar events have been reported in the literature regarding the evolution of Lüders bands in fine grained steel strips and wires. The similarity of macroscopic events, despite the different mechanisms of instability at the micro-level, prompted us to approximate the material behavior as a finitely deforming elasto-plastic solid with a trilinear up-down-up nominal stress-strain response. Two such stress-strain responses were used in finite element simulations of strip tension tests. In the first the true stress-strain response maintains its stability and in the second the intermediate branch has a negative slope. While both material models produced inhomogeneous deformations with features similar to those of the experiments, the larger initiation peak associated with the second gave results which closely resembled specific experiments. The numerical simulations confirmed that the evolution of events seen in experiments on SMAs and mild steels is strongly influenced by overall geometric (structural) effects. Furthermore, the success of this simple continuum constitutive model strongly suggests that continuum level events remain dominant players in such fine grained materials.     

Identification and characterization of delamination in polymer coated metal sheet

In this paper, a mixed numerical-experimental approach is adopted to quantitatively investigate and characterize delamination in polymer coated steel. The integral bi-material system is analysed and special attention is given to the constitutive modelling of the polymer coating, the interface and the determination of all involved parameters. An extended cohesive zone model for large displacements is proposed, allowing for a mode-dependent behaviour in large deformations. Finally, peel tests are used to characterize the interface, whereby the interfacial properties are determined through an inverse parameter identification procedure.     

Micro/macro-crack growth due to creep–fatigue dependency on time–temperature material behavior

The implicit character of micro-structural degradation is determined by specifying the time history of crack growth caused by creep–fatigue interaction at high temperature. A dual scale micro/macro-equivalent crack growth model is used to illustrate the underlying principle of multiscaling which can be applied equally well to nano/micro. A series of dual scale models can be connected to formulate triple or quadruple scale models. Temperature and time-dependent thermo-mechanical material properties are developed to dictate the design time history of creep–fatigue cracking that can serve as the master curve for health monitoring.In contrast to the conventional procedure of problem/solution approach by specifying the time- and temperature-dependent material properties as a priori, the desired solution is then defined for a class of anticipated loadings. A scheme for matching the loading history with the damage evolution is then obtained. The results depend on the initial crack size and the extent of creep in proportion to fatigue damage. The path dependent nature of damage is demonstrated by showing the range of the pertinent parameters that control the final destruction of the material. A possible scenario of 20 yr of life span for the 38Cr2Mo2VA ultra-high strength steel is used to develop the evolution of the micro-structural degradation. Three micro/macro-parameters μ^*, d^* and σ^* are used to exhibit the time-dependent variation of the material, geometry and load effects. They are necessary to reflect the scale transitory behavior of creep–fatigue damage. Once the algorithm is developed, the material can be tailor made to match the behavior. That is a different life span of the same material would alter the time behavior of μ^*, d^* and σ^* and hence the micro-structural degradation history. The one-to-one correspondence of the material micro-structure degradation history with that of damage by cracking is the essence of path dependency. Numerical results and graphs are obtained to demonstrate how the inherently implicit material micro-structure parameters can be evaluated from the uniaxial bulk material properties at the macroscopic scale.The combined behavior of creep and fatigue can be exhibited by specifying the parameter ξ with reference to the initial defect size a₀. Large ξ (0.90 and 0.85) gives critical crack size a_cr = 11–14 mm (at t < 20 yr) for a₀ about 1.3 mm. For small ξ (0.05 and 0.15), there results critical a_cr = 6–7 mm (at t < 20 yr) for a₀ about 0.7–0.8 mm. The initial crack is estimated to increase its length by an order of magnitude before triggering global to the instability. This also applies ξ ≈ 0.5 where creep interacts severely with fatigue. Fine tuning of a_cr and a₀ can be made to meet the condition oft = 20 yr.Trade off among load, material and geometric parameters are quantified such that the optimum conditions can be determined for the desired life qualified by the initial–final defect sizes. The scenario assumed in this work is indicative of the capability of the methodology. The initial–final defect sizes can be varied by re-designing the time–temperature material specifications. To reiterate, the uniqueness of solution requires the end result to match with the initial conditions for a given problem. This basic requirement has been accomplished by the dual scale micro/macro-crack growth model for creep and fatigue.     

Orthotropic yield criteria for description of the anisotropy in tension and compression of sheet metals

In this paper, yield functions describing the anisotropic behavior of textured metals are proposed. These yield functions are extensions to orthotropy of the isotropic yield function proposed by Cazacu et al. (Cazacu, O., Plunkett, B., Barlat, F., 2006. Orthotropic yield criterion for hexagonal close packed metals. Int. J. Plasticity 22, 1171–1194). Anisotropy is introduced using linear transformations of the stress deviator. It is shown that the proposed anisotropic yield functions represent with great accuracy both the tensile and compressive anisotropy in yield stresses and r-values of materials with hcp crystal structure and of metal sheets with cubic crystal structure. Furthermore, it is demonstrated that the proposed formulations can describe very accurately the anisotropic behavior of metal sheets whose tensile and compressive stresses are equal.     

A constitutive model dealing with damage due to cavity growth and the Mullins effect in rubber-like materials under triaxial loading

In this work, we attempted to describe the evolution of damage in rubber-like materials due to the Mullins effect and the cavity growth process. To this end we introduced two distinct internal variables into the constitutive laws; the first one essentially describes the Mullins damage and the second describes the cavity growth. The Mullins effect was considered as a continuous type of damage that can be modelled within the continuum damage theory. The cavity growth, being discontinuous at the microscopic scale, was also modelled by a continuous variable after a homogenization procedure. These analyses allow the establishment of a compressible constitutive law describing the strain-softening phenomena for rubber-like materials. In order to identify the material parameters and to verify the efficiency of the model, we carried out experimental studies involving uniaxial, biaxial, and hydrostatic tensions under monotonic and cyclic loading. Comparison between the model-predicted results and the experimental data shows that the present model can efficiently describe both the Mullins damage and the porosity evolution of rubber-like materials under triaxial monotonic or cyclic loading with a satisfactory accuracy. The proposed concept is simple and easy to apply to engineering calculations.     

A three-dimensional model of anisotropic hardening in metals and its application to the analysis of sheet metal formability

The hardening model proposed by Z. Mróz based on the uniaxial fatigue behavior of many metals is adopted to derive an incremental constitutive equation for general three-dimensional problems. This constitutive law is then employed in the analysis of metal forming problems to assess the influence of loading cycles, of the types involved in standard forming processes, on the ultimate formability of sheet metals. The predicted forming limit curves differ quantitatively from results obtained via an isotropie hardening model and differ qualitatively from those obtained via a kinematic model. Also investigated are the effects of such loading cycles on material response to simple tensile loading, which is often used to characterize a material. Significant differences between the present model and the other two models considered are observed in such characterizers of simple tensile behavior as the stress-strain curve, the anisotropy parameter and the uniform elongation. These differences suggest a rather simple experiment to identify the proper material model to be used in analyses of problems which involve loading cycles. Comparisons with some experimental results reveal that the employment of an anisotropic hardening model, such as the generalized Mróz model derived herein, is indeed crucial in accurately predicting material response to complicated loading histories.     

Model identification and FE simulations: Effect of different yield loci and hardening laws in sheet forming

The bi-axial experimental equipment [Flores, P., Rondia, E., Habraken, A.M., 2005a. Development of an experimental equipment for the identification of constitutive laws (Special Issue). International Journal of Forming Processes] developed by Flores enables to perform Bauschinger shear tests and successive or simultaneous simple shear tests and plane strain tests. Flores investigates the material behavior with the help of classical tensile tests and the ones performed in his bi-axial machine in order to identify the yield locus and the hardening model. With tests performed on one steel grade, the methods applied to identify classical yield surfaces such as [Hill, R., 1948. A theory of the yielding and plastic flow of anisotropic materials. Proceedings of the Royal Society of London A 193, 281–297; Hosford, W.F., 1979. On yield loci of anisotropic cubic metals. In: Proceedings of the 7th North American Metalworking Conf. (NMRC), SME, Dearborn, MI, pp. 191–197] ones as well as isotropic Swift type hardening, kinematic Armstrong–Frederick or Teodosiu and Hu hardening models are explained. Comparison with the Taylor–Bishop–Hill yield locus is also provided. The effect of both yield locus and hardening model choices is presented for two applications: plane strain tensile test and Single Point Incremental Forming (SPIF).     

Experimental study of microhole growth and coalescence in ductile materials

The growth and coalescence of two microholes in copper foil were studied experimentally byin situ tensile tests under a scanning electronic microscope. Two microholes of 15–35 μm in diameter were arranged in different distances and orientations. It was found that the mechanisms of microhole evolution were represented by slipping band creation, and then crack initiation and propagation along the slipping bands in ligament. The process of microhole growth and coalescence was influenced by the inter-center distance and orientation of microholes. The critical surface of microholes at coalescence is about 2–2.5 times that of the initial one. The variation of both the inter-center distance and orientation depends on the initial angle.     

Damage growth by debonding in a single fibre metal matrix composite: Elastoplasticity and strain energy density criterion

A displacement-based finite element-based numerical approach has been employed to study the damage growth in a unidirectional SiC/Al composite containing a pre-existing crack along the fibre/matrix interface. The composite is modeled as a two-material cylinder subjected to uniform displacement. A detailed analysis is made for the stress field in the vicinity of the debond crack tip. This approach incorporates an elastic-plastic analysis combined with a strain energy density criterion to predict debonded crack growth direction, extended stable growth and final termination. The influence of contact taking place between the debonded surfaces is also considered. It is shown that such surface contact leads to reduced stress and strain fields around the crack tip, while the extent of reduction is increased with debonding length. By combining the reduced stress field with the strain energy density criterion, a limiting value for the debonding extension can be calculated for the critical applied displacement that led to fibre fracture.     

Advances in experiments on metal sheets and tubes in support of constitutive modeling and forming simulations

This work is a review of experimental methods for observing and modeling the anisotropic plastic behavior of metal sheets and tubes under a variety of loading paths, such as biaxial compression tests; biaxial tension tests on metal sheets and tubes using closed-loop electrohydraulic testing machines; the abrupt strain path change method for detecting a yield vertex and subsequent yield loci without unloading; in-plane stress reversal tests on metal sheets; and multistage tension tests. Observed material responses are compared with the predictions of phenomenological plasticity models. Special attention is paid to the plastic deformation behavior of materials commonly used in industry, and to verifying the validity of conventional anisotropic yield criteria for those materials and associated flow rules at large plastic strains. The effects of using appropriate anisotropic yield criteria on the accuracy of simulations of forming defects, such as large springback and fracture, are also presented to highlight the importance of accurate material testing and modeling.     

Development of an Internal Bone Remodelling Theory and Applications to Some Problems in Orthopaedic Biomechanics*

Bone tissue is a porous, heterogeneous and anisotropic material, which adapts its mechanical properties depending of the local stress level. This evolutive behaviour of the bone is normally known as bone remodelling. In this work, a bone-remodelling theory, based on the principles of continuum damage mechanics, is presented. The corresponding mathematical formulation has been implemented in a finite element code in order to predict the bone response after implantation of a prosthesis or fixation. Although the present model is not based on experimental verification, the model predicts important qualitative experimental results, being still necessary to test against experimental/clinical work. The main aim of this paper has been, therefore, the qualitative study of the long-term bone evolution, especially of the human femur when different types of implants are employed. A comparative analysis between two widely used hip prostheses (the Exeter and the SHP), has been performed. We have also studied the treatment of proximal femoral fractures by means of extramedullary and intramedullary implants.     

Cavity nucleation, growth and stress field in elastic-plastic medium

A finite element method is used to study the cavity nucleation and growth in an elastic-plastic medium. The critical cohesion strength, _c at IPM (interface between second phase particle and matrix) is employed as the criterion of cavity formation. Three different values of _c are taken to examine their influence on the overall mechanical behaviour and process of cavity formation.