On the overall behavior, microstructure evolution, and macroscopic stability in reinforced rubbers at large deformations: II—Application to cylindrical fibers

In Part I of this paper, we presented a general homogenization framework for determining the overall behavior, the evolution of the underlying microstructure, and the possible onset of macroscopic instabilities in fiber-reinforced elastomers subjected to finite deformations. In this work, we make use of this framework to generate specific results for general plane-strain loading of elastomers reinforced with aligned, cylindrical fibers. For the special case of rigid fibers and incompressible behavior for the matrix phase, closed-form, analytical results are obtained. The results suggest that the evolution of the microstructure has a dramatic effect on the effective response of the composite. Furthermore, in spite of the fact that both the matrix and the fibers are assumed to be strongly elliptic, the homogenized behavior is found to lose strong ellipticity at sufficiently large deformations, corresponding to the possible development of macroscopic instabilities [Geymonat, G., Müller, S., Triantafyllidis, N., 1993. Homogenization of nonlinearly elastic materials, macroscopic bifurcation and macroscopic loss of rank-one convexity. Arch. Rat. Mech. Anal. 122, 231-290]. The connection between the evolution of the microstructure and these macroscopic instabilities is put into evidence. In particular, when the reinforced elastomers are loaded in compression along the long, in-plane axis of the fibers, a certain type of “flopping” instability is detected, corresponding to the composite becoming infinitesimally soft to rotation of the fibers.     

Homogenization-based constitutive models for porous elastomers and implications for macroscopic instabilities: II—Results

In Part I of this paper, we developed a homogenization-based constitutive model for the effective behavior of isotropic porous elastomers subjected to finite deformations. In this part, we make use of the proposed model to predict the overall response of porous elastomers with (compressible and incompressible) Gent matrix phases under a wide variety of loading conditions and initial values of porosity. The results indicate that the evolution of the underlying microstructure—which results from the finite changes in geometry that are induced by the applied loading—has a significant effect on the overall behavior of porous elastomers. Further, the model is in very good agreement with the exact and numerical results available from the literature for special loading conditions and generally improves on existing models for more general conditions. More specifically, we find that, in spite of the fact that Gent elastomers are strongly elliptic materials, the constitutive models for the porous elastomers are found to lose strong ellipticity at sufficiently large compressive deformations, corresponding to the possible onset of “macroscopic” (shear band-type) instabilities. This capability of the proposed model appears to be unique among theoretical models to date and is in agreement with numerical simulations and physical experience. The resulting elliptic and non-elliptic domains, which serve to define the macroscopic “failure surfaces” of these materials, are presented and discussed in both strain and stress space.     

Homogenization-based constitutive models for porous elastomers and implications for macroscopic instabilities: I—Analysis

The purpose of this paper is to provide homogenization-based constitutive models for the overall, finite-deformation response of isotropic porous rubbers with random microstructures. The proposed model is generated by means of the “second-order” homogenization method, which makes use of suitably designed variational principles utilizing the idea of a “linear comparison composite.” The constitutive model takes into account the evolution of the size, shape, orientation, and distribution of the underlying pores in the material, resulting from the finite changes in geometry that are induced by the applied loading. This point is key, as the evolution of the microstructure provides geometric softening/stiffening mechanisms that can have a very significant effect on the overall behavior and stability of porous rubbers. In this work, explicit results are generated for porous elastomers with isotropic, (in)compressible, strongly elliptic matrix phases. In spite of the strong ellipticity of the matrix phases, the derived constitutive model may lose strong ellipticity, indicating the possible development of shear/compaction band-type instabilities. The general model developed in this paper will be applied in Part II of this work to a special, but representative, class of isotropic porous elastomers with the objective of exploring the complex interplay between geometric and constitutive softening/stiffening in these materials.     

Multiscale modeling of oriented thermoplastic elastomers with lamellar morphology

Thermoplastic elastomers (TPEs) are block copolymers made up of “hard” (glassy or crystalline) and “soft” (rubbery) blocks that self-organize into “domain” structures at a length scale of a few tens of nanometers. Under typical processing conditions, TPEs also develop a “polydomain” structure at the micron level that is similar to that of metal polycrystals. Therefore, from a continuum point of view, TPEs may be regarded as materials with heterogeneities at two different length scales. In this work, we propose a constitutive model for highly oriented, near-single-crystal TPEs with lamellar domain morphology. Based on small-angle X-ray scattering (SAXS) and transmission electron microscopy (TEM) observations, we consider such materials to have a granular microstructure where the grains are made up of the same, perfect, lamellar structure (single crystal) with slightly different lamination directions (crystal orientations). Having identified the underlying morphology, the overall finite-deformation response of these materials is determined by means of a two-scale homogenization procedure. Interestingly, the model predictions indicate that the evolution of microstructure—especially the rotation of the layers—has a very significant, but subtle effect on the overall properties of near-single-crystal TPEs. In particular, for certain loading conditions—namely, for those with sufficiently large compressive deformations applied in the direction of the lamellae within the individual grains—the model becomes macroscopically unstable (i.e., it loses strong ellipticity). By keeping track of the evolution of the underlying microstructure, we find that such instabilities can be related to the development of “chevron” patterns.     

Homogenization-based constitutive models for magnetorheological elastomers at finite strain

This paper proposes a new homogenization framework for magnetoelastic composites accounting for the effect of magnetic dipole interactions, as well as finite strains. In addition, it provides an application for magnetorheological elastomers via a “partial decoupling” approximation splitting the magnetoelastic energy into a purely mechanical component, together with a magnetostatic component evaluated in the deformed configuration of the composite, as estimated by means of the purely mechanical solution of the problem. It is argued that the resulting constitutive model for the material, which can account for the initial volume fraction, average shape, orientation and distribution of the magnetically anisotropic, non-spherical particles, should be quite accurate at least for perfectly aligned magnetic and mechanical loadings. The theory predicts the existence of certain “extra” stresses—arising in the composite beyond the purely mechanical and magnetic (Maxwell) stresses—which can be directly linked to deformation-induced changes in the microstructure. For the special case of isotropic distributions of magnetically isotropic, spherical particles, the extra stresses are due to changes in the particle two-point distribution function with the deformation, and are of order volume fraction squared, while the corresponding extra stresses for the case of aligned, ellipsoidal particles can be of order volume fraction, when changes are induced by the deformation in the orientation of the particles. The theory is capable of handling the strongly nonlinear effects associated with finite strains and magnetic saturation of the particles at sufficiently high deformations and magnetic fields, respectively.     

A general hyperelastic model for incompressible fiber-reinforced elastomers

This work presents a new constitutive model for the effective response of fiber-reinforced elastomers at finite strains. The matrix and fiber phases are assumed to be incompressible, isotropic, hyperelastic solids. Furthermore, the fibers are taken to be perfectly aligned and distributed randomly and isotropically in the transverse plane, leading to overall transversely isotropic behavior for the composite. The model is derived by means of the “second-order” homogenization theory, which makes use of suitably designed variational principles utilizing the idea of a “linear comparison composite.” Compared to other constitutive models that have been proposed thus far for this class of materials, the present model has the distinguishing feature that it allows consideration of behaviors for the constituent phases that are more general than Neo-Hookean, while still being able to account directly for the shape, orientation, and distribution of the fibers. In addition, the proposed model has the merit that it recovers a known exact solution for the special case of incompressible Neo-Hookean phases, as well as some other known exact solutions for more general constituents under special loading conditions.     

Homogenization estimates for fiber-reinforced elastomers with periodic microstructures

This work presents a homogenization-based constitutive model for the mechanical behavior of elastomers reinforced with aligned cylindrical fibers subjected to finite deformations. The proposed model is derived by making use of the second-order homogenization method [Lopez-Pamies, O., Ponte Castañeda, P., 2006a. On the overall behavior, microstructure evolution, and macroscopic stability in reinforced rubbers at large deformations: I—theory. J. Mech. Phys. Solids 54, 807–830], which is based on suitably designed variational principles utilizing the idea of a “linear comparison composite.” Specific results are generated for the case when the matrix and fiber materials are characterized by generalized Neo-Hookean solids, and the distribution of fibers is periodic. In particular, model predictions are provided and analyzed for fiber-reinforced elastomers with Gent phases and square and hexagonal fiber distributions, subjected to a wide variety of three-dimensional loading conditions. It is found that for compressive loadings in the fiber direction, the derived constitutive model may lose strong ellipticity, indicating the possible development of macroscopic instabilities that may lead to kink band formation. The onset of shear band-type instabilities is also detected for certain in-plane modes of deformation. Furthermore, the subtle influence of the distribution, volume fraction, and stiffness of the fibers on the effective behavior and onset of macroscopic instabilities in these materials is investigated thoroughly.     

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 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.     

Finite deformation of incompressible fiber-reinforced elastomers: A computational micromechanics approach

The in-plane finite deformation of incompressible fiber-reinforced elastomers was studied using computational micromechanics. Composite microstructure was made up of a random and homogeneous dispersion of aligned rigid fibers within a hyperelastic matrix. Different matrices (Neo-Hookean and Gent), fibers (monodisperse or polydisperse, circular or elliptical section) and reinforcement volume fractions (10–40%) were analyzed through the finite element simulation of a representative volume element of the microstructure. A successive remeshing strategy was employed when necessary to reach the large deformation regime in which the evolution of the microstructure influences the effective properties. The simulations provided for the first time “quasi-exact” results of the in-plane finite deformation for this class of composites, which were used to assess the accuracy of the available homogenization estimates for incompressible hyperelastic composites.     

Brain tissue deforms similarly to filled elastomers and follows consolidation theory

Slow, large deformations of human brain tissue—accompanying cranial vault deformation induced by positional plagiocephaly, occurring during hydrocephalus, and in the convolutional development—has surprisingly received scarce mechanical investigation. Since the effects of these deformations may be important, we performed a systematic series of in vitro experiments on human brain tissue, revealing the following features. (i) Under uniaxial (quasi-static), cyclic loading, brain tissue exhibits a peculiar nonlinear mechanical behaviour, exhibiting hysteresis, Mullins effect and residual strain, qualitatively similar to that observed in filled elastomers. As a consequence, the loading and unloading uniaxial curves have been found to follow the Ogden nonlinear elastic theory of rubber (and its variants to include Mullins effect and permanent strain). (ii) Loaded up to failure, the “shape” of the stress/strain curve qualitatively changes, evidencing softening related to local failure. (iii) Uniaxial (quasi-static) strain experiments under controlled drainage conditions provide the first direct evidence that the tissue obeys consolidation theory involving fluid migration, with properties similar to fine soils, but having much smaller volumetric compressibility. (iv) Our experimental findings also support the existence of a viscous component of the solid phase deformation.Brain tissue should, therefore, be modelled as a porous, fluid-saturated, nonlinear solid with very small volumetric (drained) compressibility.     

Analytical study of spherical cloak/anti-cloak interactions

The intriguing concept of “anti-cloaking” has been recently introduced within the framework of transformation optics (TO), first as a “countermeasure” to invisibility-cloaking (i.e., to restore the scattering response of a cloaked target), and more recently in connection with “sensor invisibility” (i.e., to strongly reduce the scattering response while maintaining the field-sensing capabilities). In this paper, we extend our previous studies, which were limited to a two-dimensional cylindrical scenario, to the three-dimensional spherical case. More specifically, via a generalized (coordinate-mapped) Mie-series approach, we derive a general analytical full-wave solution pertaining to plane-wave-excited configurations featuring a spherical object surrounded by a TO-based invisibility cloak coupled via a vacuum layer to an anti-cloak, and explore the various interactions of interest. With a number of selected examples, we illustrate the cloaking and field-restoring capabilities of various configurations, highlighting similarities and differences with respect to the cylindrical case, with special emphasis on sensor-cloaking scenarios and ideas for approximate implementations that require the use of double-positive media only.     

An atomistic investigation of the kinetics of detwinning

This paper presents a “first principles” atomistic study of the dynamics of detwinning in a shape-memory alloy. In order to describe the macroscopic motion of twin boundaries, the continuum theory of twinning must be provided with a “kinetic relation”, i.e. a relation between the driving force and the propagation speed. This kinetic relation is a macroscopic characterization of the underlying atomistic processes. The goal of the present atomistic study is to provide the continuum theory with this kinetic relation by extracting the essential macroscopic features of the dynamics of the atoms. It also aims to elucidate the mechanism underlying the process of detwinning.The material studied is stoichiometric nickel-manganese, and interatomic interactions are described using three physically motivated Lennard-Jones potentials. The effect of temperature and shear stress on detwinning — specifically on the rate of transformation from one variant of martensite to the other — is examined using molecular dynamics. An explicit formula for this (kinetic) relation is obtained by fitting an analytic expression to the simulation results. The numerical experiments also verify that transverse ledge propagation is the mechanism underlying twin-boundary motion. All calculations are carried out in a two-dimensional setting.     

On the possibility of piezoelectric nanocomposites without using piezoelectric materials

In this work, predicated on nanoscale size-effects, we explore the tantalizing possibility of creating apparently piezoelectric composites without using piezoelectric constituent materials. In a piezoelectric material an applied uniform strain can induce an electric polarization (or vice-versa). Crystallographic considerations restrict this technologically important property to non-centrosymmetric systems. Non-uniform strain can break the inversion symmetry and induce polarization even in non-piezoelectric dielectrics. The key concept is that all dielectrics (including non-piezoelectric ones) exhibit the aforementioned coupling between strain gradient and polarization—an experimentally verified phenomenon known in some circles as the flexoelectric effect. This flexoelectric coupling, however, is generally very small and evades experimental detection unless very large strain gradients (or conversely polarization gradients) are present. Based on a field theoretic framework and the associated Greens function solutions developed in prior work, we quantitatively demonstrate the possibility of “designing piezoelectricity,” i.e. we exploit the large strain gradients present in the interior of composites containing nanoscale inhomogeneities to achieve an overall non-zero polarization even under an uniformly applied stress. We prove that the aforementioned effect may be realized only if both the shapes and distributions of the inhomogeneities are non-centrosymmetric. Our un-optimized quantitative results, based on limited material data and restrictive assumptions on inhomogeneity shape and distribution, indicate that apparent piezoelectric behavior close to 10% of Quartz may be achievable for inhomogeneity sizes in the 4 nm range. In future works, it is not unreasonable to expect enhanced performance based on optimization of shape, topology and appropriate material selection.     

Onset of macroscopic instabilities in fiber-reinforced elastomers at finite strain

In this paper, we investigate theoretically the possible development of instabilities in fiber-reinforced elastomers (and other soft materials) when they are subjected to finite-strain loading conditions. We focus on the physically relevant class of “macroscopic” instabilities, i.e., instabilities with wavelengths that are much larger than the characteristic size of the underlying microstructure. To this end, we make use of recently developed homogenization estimates, together with a fundamental result of Geymonat, Müller and Triantafyllidis linking the development of these instabilities to the loss of strong ellipticity of the homogenized constitutive relations. For the important class of material systems with very stiff fibers and random microstructures, we derive a closed-form formula for the critical macroscopic deformation at which instabilities may develop under general loading conditions, and we show that this critical deformation is quite sensitive to the loading orientation relative to the fiber direction. The result is also confronted with classical estimates (including those of Rosen) for laminates, which have commonly been used as two-dimensional (2-D) approximations for actual fiber-reinforced composites. We find that while predictions based on laminate models are qualitatively correct for certain loadings, they can be significantly off for other more general 3-D loadings. Finally, we provide a parametric analysis of the effects of the matrix and fiber properties and of the fiber volume fraction on the onset of instabilities for various loading conditions.     

Size effects in generalised continuum crystal plasticity for two-phase laminates

The solutions of a boundary value problem are explored for various classes of generalised crystal plasticity models including Cosserat, strain gradient and micromorphic crystal plasticity. The considered microstructure consists of a two-phase laminate containing a purely elastic and an elasto-plastic phase undergoing single or double slip. The local distributions of plastic slip, lattice rotation and stresses are derived when the microstructure is subjected to simple shear. The arising size effects are characterised by the overall extra back stress component resulting from the action of higher order stresses, a characteristic length l_c describing the size-dependent domain of material response, and by the corresponding scaling law lⁿ as a function of microstructural length scale, l. Explicit relations for these quantities are derived and compared for the different models. The conditions at the interface between the elastic and elasto-plastic phases are shown to play a major role in the solution. A range of material parameters is shown to exist for which the Cosserat and micromorphic approaches exhibit the same behaviour. The models display in general significantly different asymptotic regimes for small microstructural length scales. Scaling power laws with the exponent continuously ranging from 0 to −2 are obtained depending on the values of the material parameters. The unusual exponent value −2 is obtained for the strain gradient plasticity model, denoted “curl H^p” in this work. These results provide guidelines for the identification of higher order material parameters of crystal plasticity models from experimental data, such as precipitate size effects in precipitate strengthened alloys.     

One, no one, and one hundred thousand crack propagation laws: A generalized Barenblatt and Botvina dimensional analysis approach to fatigue crack growth

Barenblatt and Botvina with elegant dimensional analysis arguments have elucidated that Paris’ power-law is a weak form of scaling, so that the Paris’ parameters C and m should not be taken as material constants. On the contrary, they are expected to depend on all the dimensionless parameters of the problem, and are really “constants” only within some specific ranges of all these. In the present paper, the dimensional analysis approach by Barenblatt and Botvina is generalized to explore the functional dependencies of m and C on more dimensionless parameters than the original Barenblatt and Botvina, and experimental results are interpreted for a wider range of materials including both metals and concrete. In particular, we find that the size-scale dependencies of m and C and the resulting correlation between C and m are quite different for metals and for quasi-brittle materials, as it is already suggested from the fact the fatigue crack propagation processes lead to m=2-5 in metals and m=10-50 in quasi-brittle materials. Therefore, according to the concepts of complete and incomplete self-similarities, the experimentally observed breakdowns of the classical Paris’ law are discussed and interpreted within a unified theoretical framework. Finally, we show that most attempts to address the deviations from the Paris’ law or the empirical correlations between the constants can be explained with this approach. We also suggest that “incomplete similarity” corresponds to the difficulties encountered so far by the “damage tolerant” approach which, after nearly 50 years since the introduction of Paris’ law, is still not a reliable calculation of damage, as Paris himself admits in a recent review.     

Integration of self-consistent polycrystal plasticity with dislocation density based hardening laws within an implicit finite element framework: Application to low-symmetry metals

We present an implementation of the viscoplastic self-consistent (VPSC) polycrystalline model in an implicit finite element (FE) framework, which accounts for a dislocation-based hardening law for multiple slip and twinning modes at the micro-scale grain level. The model is applied to simulate the macro-scale mechanical response of a highly anisotropic low-symmetry (orthorhombic) crystal structure. In this approach, a finite element integration point represents a polycrystalline material point and the meso-scale mechanical response is obtained by the mean-field VPSC homogenization scheme. We demonstrate the accuracy of the FE-VPSC model by analyzing the mechanical response and microstructure evolution of α-uranium samples under simple compression/tension and four-point bending tests. Predictions of the FE-VPSC simulations compare favorably with experimental measurements of geometrical changes and microstructure evolution. Specifically, the model captures accurately the tension–compression asymmetry of the material associated with twinning, as well as the rigidity of the material response along the hard-to-deform crystallographic orientations.     

Cavitation in elastomeric solids: I—A defect-growth theory

It is by now well established that loading conditions with sufficiently large triaxialities can induce the sudden appearance of internal cavities within elastomeric (and other soft) solids. The occurrence of such instabilities, commonly referred to as cavitation, can be attributed to the growth of pre-existing defects into finite sizes. This paper introduces a new theory to study the phenomenon of cavitation in soft solids that: (i) allows to consider general 3D loading conditions with arbitrary triaxiality, (ii) applies to large (including compressible and anisotropic) classes of nonlinear elastic solids, and (iii) incorporates direct information on the initial shape, spatial distribution, and mechanical properties of the underlying defects at which cavitation can initiate. The basic idea is to first cast cavitation in elastomeric solids as a homogenization problem of nonlinear elastic materials containing random distributions of zero-volume cavities, or defects. This problem is then addressed by means of a novel iterated homogenization procedure, which allows to construct solutions for a specific, yet fairly general, class of defects. These include solutions for the change in size of the defects as a function of the applied loading conditions, from which the onset of cavitation — corresponding to the event when the initially infinitesimal defects suddenly grown into finite sizes — can be readily determined. In spite of the generality of the proposed approach, the relevant calculations amount to solving tractable Hamilton-Jacobi equations, in which the initial size of the defects plays the role of “time” and the applied load plays the role of “space”. When specialized to the case of hydrostatic loading conditions, isotropic solids, and defects that are vacuous and isotropically distributed, the proposed theory recovers the classical result of Ball (1982) for radially symmetric cavitation. The nature and implications of this remarkable connection are discussed in detail.     

A mathematical model for cylindrical,fiber reinforced electro-pneumatic actuators

In recent years, dielectric elastomers have received increasing attention due to their unparalleled large strain actuation response (>100%). The force output, however, has remained a major limiting factor for many applications. To address this limitation, a model for a fiber reinforced dielectric elastomer actuator based on the deformation mechanism of McKibben actuators is presented. In this novel configuration, the outer cylindrical surface of a dielectric elastomer is enclosed by a network of helical fibers that are thin, flexible and inextensible. This configuration yields an axially contractile actuator, in contrast to unreinforced actuators which extend. The role of the fiber network is twofold: (i) to serve as reinforcement to improve the load-bearing capability of dielectric elastomers, and (ii) to render the actuator inextensible in the axial direction such that the only free deformation path is simultaneous radial expansion and axial contraction. In this paper, a mathematical model of the electromechanical response of fiber reinforced dielectric elastomers is derived. The model is developed within a continuum mechanics framework for large deformations. The cylindrical electro-pneumatic actuator is modeled by adapting Green and Adkins’ theory of reinforced cylinders to account for the applied electric field. Using this approach, numerical solutions are obtained assuming a Mooney–Rivlin material model. The results indicate that the relationship between the contractile force and axial shortening is bilinear within the voltage range considered. The characteristic response as a function of various system parameters such as the fiber angle, inflation pressure, and the applied voltage are reported. In this paper, the elastic portion of the modeling approach is validated using experimental data for McKibben actuators.