Buckling of nano-fibre reinforced composites: a re-examination of elastic buckling

Elastic buckling of layered/fibre reinforced composites is investigated. Assuming the existence of both shear and transverse modes of failure, the fibre is analysed as a layer embedded in a matrix. Interacting stresses, acting at the interfaces are determined from an exact derived stress field in the matrix. It is shown that buckling can occur only in the shear buckling mode and that the transverse buckling mode is spurious. As opposed to the well known Rosen shear buckling mode solution (predicated on an infinite buckling wavelength), shear buckling is shown to exist under two régimes: buckling of dilute composites with finite wavelengths and buckling of non-dilute composites with infinite wavelengths. Based on the analysis, a model is constructed which defines the fibre concentration at which the transition between the two régimes occurs. The buckling strains are shown to be (approximately) constant for dilute composites and, in the case of very stiff fibres, to have realistic values compatible with elastic behaviour. For the case of non-dilute composites, the strains are found to be in agreement with those given by the Rosen shear buckling solution. Numerical results for the buckling strains and stresses are presented and compared with the Rosen solution. These reveal that the Rosen solution is valid only for the case of non-dilute composites. The investigation demonstrates that elastic buckling may be a dominant failure mechanism of composites consisting of very stiff fibres fabricated in the framework of nano-technology.     

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.     

Microscopic and macroscopic instabilities in finitely strained fiber-reinforced elastomers

The present work is a detailed study of the connections between microstructural instabilities and their macroscopic manifestations — as captured through the effective properties — in finitely strained fiber-reinforced elastomers, subjected to finite, plane-strain deformations normal to the fiber direction. The work, which is a complement to a previous and analogous investigation by the same authors on porous elastomers, (Michel et al., 2007), uses the linear comparison, second-order homogenization (S.O.H.) technique, initially developed for random media, to study the onset of failure in periodic fiber-reinforced elastomers and to compare the results to more accurate finite element method (F.E.M.) calculations. The influence of different fiber distributions (random and periodic), initial fiber volume fraction, matrix constitutive law and fiber cross-section on the microscopic buckling (for periodic microgeometries) and macroscopic loss of ellipticity (for all microgeometries) is investigated in detail. In addition, constraints to the principal solution due to fiber/matrix interface decohesion, matrix cavitation and fiber contact are also addressed. It is found that both microscopic and macroscopic instabilities can occur for periodic microstructures, due to a symmetry breaking in the periodic arrangement of the fibers. On the other hand, no instabilities are found for the case of random microstructures with circular section fibers, while only macroscopic instabilities are found for the case of elliptical section fibers, due to a symmetry breaking in their orientation.     

Microscopic and macroscopic instabilities in finitely strained porous elastomers

The present work is an in-depth study of the connections between microstructural instabilities and their macroscopic manifestations—as captured through the effective properties—in finitely strained porous elastomers. The powerful second-order homogenization (SOH) technique initially developed for random media, is used for the first time here to study the onset of failure in periodic porous elastomers and the results are compared to more accurate finite element method (FEM) calculations. The influence of different microgeometries (random and periodic), initial porosity, matrix constitutive law and macroscopic load orientation on the microscopic buckling (for periodic microgeometries) and macroscopic loss of ellipticity (for all microgeometries) is investigated in detail. In addition to the above-described stability-based onset-of-failure mechanisms, constraints on the principal solution are also addressed, thus giving a complete picture of the different possible failure mechanisms present in finitely strained porous elastomers.     

Microscopic symmetric bifurcation condition of cellular solids based on a homogenization theory of finite deformation

In this paper, we establish a homogenization framework to analyze the microscopic symmetric bifurcation buckling of cellular solids subjected to macroscopically uniform compression. To this end, describing the principle of virtual work for infinite periodic materials in the updated Lagrangian form, we build a homogenization theory of finite deformation, which satisfies the principle of material objectivity. Then, we state a postulate that at the onset of microscopic symmetric bifurcation, microscopic velocity becomes spontaneous, yet changing the sign of such spontaneous velocity has no influence on the variation in macroscopic states. By applying this postulate to the homogenization theory, we derive the conditions to be satisfied at the onset of microscopic symmetric bifurcation. The resulting conditions are verified by analyzing numerically the in-plane biaxial buckling of an elastic hexagonal honeycomb. It is thus shown that three kinds of experimentally observed buckling modes of honeycombs i.e., uniaxial, biaxial and flower-like modes, are attained and classified as microscopic symmetric bifurcation. It is also shown that the multiplicity of bifurcation gives rise to the complex cell-patterns in the biaxial and flower-like modes.     

Onset of necking in electro-magnetically formed rings

The electromagnetic forming (EMF) process is an attractive manufacturing technique, which uses electromagnetic (Lorentz) body forces to fabricate metallic parts. One of the many advantages of EMF is the considerable ductility increase observed in several metals, with aluminum featuring prominently among them. The quantitative explanation of this phenomenon is the primary motivation of this work.To study the ductility increase due to EMF, an aluminum ring is placed outside a fixed coil, which is connected to a capacitor. Upon the capacitor's discharge, the time varying current in the coil induces a larger current in the ring specimen and the resulting Lorentz forces make it expand. The coupled coil-ring electromagnetic and thermomechanical problem is solved, using an experimentally obtained constitutive model for a particular aluminum alloy. Our results show that for realistic imperfections, the EMF ring starts necking at strains about six times larger than its static counterpart, as observed experimentally. This study establishes the importance of solving the fully coupled electromagnetic and thermomechanical problem and provides insight on how different constitutive parameters influence ductility in an EMF process.     

Effects of texture on shear band formation in plane strain tension/compression and bending

In this study, effects of typical texture components observed in rolled aluminum alloy sheets on shear band formation in plane strain tension/compression and bending are systematically studied. The material response is described by a generalized Taylor-type polycrystal model, in which each grain is characterized in terms of an elastic–viscoplastic continuum slip constitutive relation. First, a simple model analysis in which the shear band is assumed to occur in a weaker thin slice of material is performed. From this simple model analysis, two important quantities regarding shear band formation are obtained: i.e. the critical strain at the onset of shear banding and the corresponding orientation of shear band. Second, the shear band development in plane strain tension/compression is analyzed by the finite element method. Predictability of the finite element analysis is compared to that of the simple model analysis. Third, shear band developments in plane strain pure bending of a sheet specimen with the typical textures are studied. Regions near the surfaces in a bent sheet specimen are approximately subjected to plane strain tension or compression. From this viewpoint, the bendability of a sheet specimen may be evaluated, using the knowledge regarding shear band formation in plane strain tension/compression. To confirm this and to encompass overall deformation of a bent sheet specimen, including shear bands, finite element analyses of plane strain pure bending are carried out, and the predicted shear band formation in bent specimens is compared to that in the tension/compression problem. Finally, the present results are compared to previous related studies, and the efficiency of the present method for materials design in future is discussed.     

Experimental observations of evolving yield loci in biaxially strained AA5754-O

Experimental measurement of the plastic biaxial mechanical response for an aluminum alloy (AA5754-O) sheet metal is presented. Traditional methods of multiaxial sheet metal testing require the use of finite element analysis (FEA) or other assumptions of material response to determine the multiaxial true stress versus true strain behavior of the as-received sheet material. The method used here strives to produce less ambiguous measurements of data for a larger strain range than previously possible, through a combination of the Marciniak flat bottom ram test and an X-ray diffraction technique for stress measurement. The study is performed in conjunction with a study of the microstructural changes that occur during deformation, and these microstructural results are briefly mentioned in this work. Issues of calibration and applicability are discussed, and results are presented for uniaxial (U), plane strain (PS), and balanced biaxial (BB) extension. The results show repeatable behavior (within quantified uncertainties) for U to 20%, PS to almost 15%, and BB to above 20% in-plane strains. The results are first compared with three common yield locus models (von Mises’, Hill’48, and Hosford’79), and show some unexpected results in the shape change of the yield locus at high strain levels (>5% strain). These changes include the rotation of the locus toward the von Mises surface and elongation in the balanced biaxial direction. Comparison with a more complex yield locus model (Yld2000-2d with eight adjustable parameters) showed that the locus elongation in the biaxial direction could be fit well (for a specific level of work), but at the detriment of fit to the plane strain data. Artificially large plastic strain ratios would be needed to match both the biaxial and plane strain behavior even with this more complex model.     

Effective pressure-sensitive elastoplastic behavior of particle-reinforced composites and porous media under isotropic loading

The composite under investigation consists of an elastoplastic matrix reinforced by elastic particles or weakened by pores. The material forming the matrix is pressure-sensitive. The Drucker–Prager yield criterion and a one-parameter non-associated flow rule are employed to formulate the yield behavior of the matrix. The objective of this work is to estimate the effective elastoplastic behavior of the composite under isotropic tensile and compressive loadings. To achieve this objective, the composite sphere assemblage model of Hashin [Z. Hashin, The elastic moduli of heterogeneous materials, ASME J. Appl. Mech. 29 (1962) 143–150] is used. Exact solutions are thus derived as estimations for the effective secant and tangent bulk moduli of the composite. The effects of the loading modes and phase properties on the effective elastoplastic behavior of the composite are analytically and numerically evaluated.     

On electromagnetic forming processes in finitely strained solids: Theory and examples

The process of electromagnetic forming (EMF) is a high velocity manufacturing technique that uses electromagnetic (Lorentz) body forces to shape sheet metal parts. EMF holds several advantages over conventional forming techniques: speed, repeatability, one-sided tooling, and most importantly considerable ductility increase in several metals. Current modeling techniques for EMF processes are not based on coupled variational principles to simultaneously account for electromagnetic and mechanical effects. Typically, separate solutions to the electromagnetic (Maxwell) and motion (Newton) equations are combined in staggered or lock-step methods, sequentially solving the mechanical and electromagnetic problems.The present work addresses these issues by introducing a fully coupled Lagrangian (reference configuration) least-action variational principle, involving magnetic flux and electric potentials and the displacement field as independent variables. The corresponding Euler-Lagrange equations are Maxwell's and Newton's equations in the reference configuration, which are shown to coincide with their current configuration counterparts obtained independently by a direct approach. The general theory is subsequently simplified for EMF processes by considering the eddy current approximation.Next, an application is presented for axisymmetric EMF problems. It is shown that the proposed variational principle forms the basis of a variational integration numerical scheme that provides an efficient staggered solution algorithm. As an illustration a number of such processes are simulated, inspired by recent experiments of freely expanding uncoated and polyurea-coated aluminum tubes.     

Elastoplastic microscopic bifurcation and post-bifurcation behavior of periodic cellular solids

In this study, a general framework is developed to analyze microscopic bifurcation and post-bifurcation behavior of elastoplastic, periodic cellular solids. The framework is built on the basis of a two-scale theory, called a homogenization theory, of the updated Lagrangian type. We thus derive the eigenmode problem of microscopic bifurcation and the orthogonality to be satisfied by the eigenmodes. The orthogonality allows the macroscopic increments to be independent of the eigenmodes, resulting in a simple procedure of the elastoplastic post-bifurcation analysis based on the notion of comparison solids. By use of this framework, then, bifurcation and post-bifurcation analysis are performed for cell aggregates of an elastoplastic honeycomb subject to in-plane compression. Thus, demonstrating a basic, long-wave eigenmode of microscopic bifurcation under uniaxial compression, it is shown that the eigenmode has the longitudinal component dominant to the transverse component and consequently causes microscopic buckling to localize in a cell row perpendicular to the loading axis. It is further shown that under equi-biaxial compression, the flower-like buckling mode having occurred in a macroscopically stable state changes into an asymmetric, long-wave mode due to the sextuple bifurcation in a macroscopically unstable state, leading to the localization of microscopic buckling in deltaic areas.     

Kinetic wrinkling of an elastic film on a viscoelastic substrate

A compressed elastic film on a compliant substrate can form wrinkles. On an elastic substrate, equilibrium and energetics set the critical condition and select the wrinkle wavelength and amplitude. On a viscous substrate, wrinkle grows over time and the kinetics selects the fastest growing wavelength. More generally, on a viscoelastic substrate, both energetics and kinetics play important roles in determining the critical condition, the growth rate, and the wavelength. This paper studies the wrinkling process of an elastic film on a viscoelastic layer, which in turn lies on a rigid substrate. The film is elastic and modeled by the nonlinear von Karman plate theory. The substrate is linear viscoelastic with a relaxation modulus typical of a cross-linked polymer. Beyond a critical stress, the film wrinkles by the out-of-plane displacement but remains bonded to the substrate. This study considers plane strain wrinkling and neglects the in-plane displacement. A classification of the wrinkling behavior is made based on the critical conditions at the elastic limits, the glassy and rubbery states of the viscoelastic substrate. Linear perturbation analyses are conducted to reveal the kinetics of wrinkling in films subjected to intermediate and large compressive stresses. It is shown that, depending on the stress level, the growth of wrinkles at the initial stage can be exponential, accelerating, linear, or decelerating. In all cases, the wrinkle amplitude saturates at an equilibrium state after a long time. Subsequently, both amplitude and wavelength of the wrinkle evolve, but the process is kinetically constrained and slow compared to the initial growth.     

An experimental investigation of shock wave propagation in periodically layered composites

In heterogeneous media, scattering due to interfaces/microstructure between dissimilar materials could play an important role in shock wave dissipation and dispersion. In this work, the influence of interface scattering on finite-amplitude shock waves was experimentally investigated by impacting flyer plates onto periodically layered polycarbonate/6061 aluminum, polycarbonate/304 stainless steel and polycarbonate/glass composites. Experimental results (obtained using velocity interferometer and stress gage) show that these periodically layered composites can support steady structured shock waves. Due to interface scattering, the effective shock viscosity increases with the increase of interface impedance mismatch, and decreases with the increase of interface density (interface area per unit volume) and loading amplitude. For the composites studied here, the strain rate within the shock front is roughly proportional to the square of the shock stress. This indicates that layered composites have much larger shock viscosity due to the interface/microstructure scattering in comparison with the increase of shock strain rate by the fourth power of the shock stress for homogeneous metals. Experimental results also show that due to the scattering effects, shock propagation in the layered composites is dramatically slowed down and the shock speed in composites can be lower than that either of its components.     

Modelling of laminated microstructures in stress-induced martensitic transformations

This paper is concerned with micromechanical modelling of stress-induced martensitic transformations in crystalline solids, with the focus on distinct elastic anisotropy of the phases and the associated redistribution of internal stresses. Micro-macro transition in stresses and strains is analysed for a laminated microstructure of austenite and martensite phases. Propagation of a phase transformation front is governed by a time-independent thermodynamic criterion. Plasticity-like macroscopic constitutive rate equations are derived in which the transformed volume fraction is incrementally related to the overall strain or stress. As an application, numerical simulations are performed for cubic β₁ (austenite) to orthorhombic γ₁′ (martensite) phase transformation in a single crystal of Cu-Al-Ni shape memory alloy. The pseudoelasticity effect in tension and compression is investigated along with the corresponding evolution of internal stresses and microstructure.     

Comparison of the work-hardening of metallic sheets using tensile and shear strain paths

This work deals with the characterization of the kinematic work-hardening of a bake-hardening steel. A shear test device has been designed and its use for the characterization of the work-hardening of sheet metals is described. Two main results are presented. Firstly, a local strain measurement, based on the following of three dots drawn on the gauge area, gives the evolution of the strain tensor eigenvalues during the test. It is shown, by comparing the theoretical kinematics of simple shear with a slightly perturbated one, that the strain state is close to the ideal one in the center of the gauge area. Secondly, reversal of the shear direction is performed after several prestrain and the evolution of the kinematic work-hardening with the equivalent plastic strain has been identified using an anisotropic elasto-viscoplastic model of Hill 1948 type. Isotropic and kinematic contributions of the work-hardening are also calculated from loading–unloading tensile tests and are compared to those obtained from the simple shear tests. The results show a discrepancy between both identification for the isotropic and the kinematic hardening. However, they are in agreement concerning the evolution of the global work-hardening.     

Kink band instability in layered structures

A recent two-dimensional prototype model for the initiation of kink banding in compressed layered structures is extended to embrace the two propagation mechanisms of band broadening and band progression. As well as interlayer friction, overburden pressure and layer bending energy, the characteristics of transverse layer compressibility and foundation stiffness are now included. Experiments on constrained layers of paper show good agreement with the predictions of angle of orientation, kink band width and post-kink load-deflection response obtained from the model.     

Multiscale modeling of plasticity based on embedding the viscoplastic self-consistent formulation in implicit finite elements

This paper is concerned with the multiscale simulation of plastic deformation of metallic specimens using physically-based models that take into account their polycrystalline microstructure and the directionality of deformation mechanisms acting at single-crystal level. A polycrystal model based on self-consistent homogenization of single-crystal viscoplastic behavior is used to provide a texture-sensitive constitutive response of each material point, within a boundary problem solved with finite elements (FE) at the macroscale. The resulting constitutive behavior is that of an elasto-viscoplastic material, implemented in the implicit FE code ABAQUS. The widely-used viscoplastic selfconsistent (VPSC) formulation for polycrystal deformation has been implemented inside a user-defined material (UMAT) subroutine, providing the relationship between stress and plastic strain-rate response. Each integration point of the FE model is considered as a polycrystal with a given initial texture that evolves with deformation. The viscoplastic compliance tensor computed internally in the polycrystal model is in turn used for the minimization of a suitable-designed residual, as well as in the construction of the elasto-viscoplastic tangent stiffness matrix required by the implicit FE scheme.Uniaxial tension and simple shear of an FCC polycrystal have been used to benchmark the accuracy of the proposed implicit scheme and the correct treatment of rotations for prediction of texture evolution. In addition, two applications are presented to illustrate the potential of the multiscale strategy: a simulation of rolling of an FCC plate, in which the model predicts the development of different textures through the thickness of the plate; and the deformation under 4-point bending of textured HCP bars, in which the model captures the dimensional changes associated with different orientations of the dominant texture component with respect to the bending plane.     

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.