Effects of pre-stress on crack-tip fields in elastic,incompressible solids

A closed-form asymptotic solution is provided for velocity fields and the nominal stress rates near the tip of a stationary crack in a homogeneously pre-stressed configuration of a nonlinear elastic, incompressible material. In particular, a biaxial pre-stress is assumed with stress axes parallel and orthogonal to the crack faces. Two boundary conditions are considered on the crack faces, namely a constant pressure or a constant dead loading, both preserving an homogeneous ground state. Starting from this configuration, small superimposed Mode I or Mode II deformations are solved, in the framework of Biot's incremental theory of elasticity. In this way a definition of an incremental stress intensity factor is introduced, slightly different for pressure or dead loading conditions on crack faces. Specific examples are finally developed for various hyperelastic materials, including the J₂-deformation theory of plasticity. The presence of pre-stress is shown to strongly influence the angular variation of the asymptotic crack-tip fields, even if the nominal stress rate displays a square root singularity as in the infinitesimal theory. Relationships between the solution with shear band formation at the crack tip and instability of the crack surfaces are given in evidence.     

Evolution of subsequent yield surfaces and elastic constants with finite plastic deformation. Part III: Yield surface in tension–tension stress space (Al 6061–T 6511 and annealed 1100 Al)

Subsequent yield surfaces for aluminum alloys are determined for three proportional loading paths (i.e., axial, hoop, and combined hoop and axial stress) using 10 με deviation from linearity as the definition of yield. This paper is in continuation with Parts I and II of the author’s previous papers on subsequent yield surfaces under tension–torsion (σ11–√3σ₁₂) stress space using similar small offset definition of yield. In the current paper comprehensive experimental results on subsequent yield surfaces under tension–tension (σ₁₁–σ₂₂) stress space are presented. Comparison of subsequent yield surfaces under (σ₁₁–√3σ₁₂) stress space, investigated in the earlier papers, clearly indicated distinctive hardening behavior under various loading paths. However, subsequent yield surfaces for Al 6061–T 6511 (a low work hardening alloy) showed contraction and negative cross-effect with finite deformation as compared to the annealed 1100 Al (a high work hardening alloy) where expansion and positive cross-effect was observed.     

The torsional buckling of a cruciform column under compressive load with a vertex plasticity model

The torsional buckling of a plastically deforming cruciform column under compressive load is investigated. The problem is solved analytically based on the von Kármán shallow shell theory and the virtual work principle. Solutions found in the literature are extended for path-dependent incremental behaviour as typically found in the presence of the vertex effect that is present in metallic polycrystals.At the critical load for buckling the direction of straining changes by an additional shear component. It is shown that the incremental elastic–plastic moduli are spatially nonuniform for such situations, contrary to the classical J₂ flow and deformation theories. The critical shear modulus that governs the buckling equation is obtained as a weighted average of the incremental elastic–plastic moduli over the cross-section of the cruciform.Using a plasticity model proposed by the authors, that includes the vertex effect, the buckling-critical load is computed for a aluminium column both with the analytical model and a FEM-based eigenvalue buckling analysis. The stable post-buckling path is determined by the energy criterion of path-stability. A comparison with the experimentally obtained classical results by Gerard and Becker (1957) shows good agreement without relying on artificial imperfections as necessary in the classical J₂ flow theory.     

Plasticity of a high strength steel alloy

Experimentally determined plastic constitutive equations of the parabolic form (σ − σ_y = β(ε − ε_y)1/2) are presented for a high strength alloy steel. Two deformation moduli β were required to describe the quasistatic compression data, both of which, as well as the point of change, were predicted by a mode index and transition strain structure of a general theory of plasticity. Dynamic strain, duration of impact, and final strain distribution were measured on specimen rods subjected to axial, symmetric, and constant velocity impact. The dynamic yield stress was 16% higher than the quasistatic value. The dynamic response function, deduced from a simple wave propagation theory, was also parabolic and a single deformation modulus, equal to the initial quasistatic value, applied. Thus, it was shown that the form of the quasistatic response function was preserved in dynamic loading, and that the increase of the dynamic stress was due to the increase of the yield stress.     

Analyses of steady quasistatic crack growth in plane strain tension in elastic-plastic materials with non-isotropic hardening

Steady state crack propagation problems of elastic-plastic materials in Mode I, plane strain under small scale yielding conditions were investigated with the aid of the finite element method. The elastic-perfectly plastic solution shows that elastic unloading wedges subtended by the crack tip in the plastic wake region do exist and that the stress state around the crack tip is similar to the modified Prandtl fan solution. To demonstrate the effects of a vertex on the yield surface, the small strain version of a phenomenological J₂, corner theory of plasticity (Christoffersen, J. and Hutchinson, J. W. J. Mech. Phys. Solids,27, 465 C 1979) with a power law stress strain relation was used to govern the strain hardening of the material. The results are compared with the conventional J₂ incremental plasticity solution. To take account of Bauschinger like effects caused by the stress history near the crack tip, a simple kinematic hardening rule with a bilinear stress strain relation was also studied. The results are again compared with the smooth yield surface isotropic hardening solution for the same stress strain curve. There appears to be more potential for steady state crack growth in the conventional J₂ incremental plasticity material than in the other two plasticity laws considered here if a crack opening displacement fracture criterion is used. However, a fracture criterion dependent on both stress and strain could lead to a contrary prediction.     

An experimental investigation on the tensile moduli and strengths of graphite/epoxy laminates

The results of a series of tensile tests on graphite/epoxy [0°]_8S and [0°/±30°/0°]_2S laminates at rates varying from 0.002 in/min to 2 in./min are reported. The loads are applied at various angles to the fiber directions in each case. The rate-dependent behavior of the stress-strain response is assessed. Evidence is presented to indicate that failure first occurs on inner plies. Also, evidence is presented to indicate that, in some cases, moduli increase with increased stress or strain level. Lamination theory is used to predict moduli and comparisons with experiment are given. This theory is also used in conjunction with three failure theories to predict ultimate strengths with varying degrees of success. Further, two approaches to ply unloading after first-ply failure are used and discussed. One is a standard method found in the literature while the other is a proposed ‘strength-of-materials’ type of technique which is computationally much simpler.     

The temperature dependence of yield

A theoretically rigorous expression is derived from the fundamental equation of rate processes to describe the temperature dependence of the upper and lower yield stress in terms of physical quantities. The effect of the activation volume and activation energy on the yield behavior was investigated over a wide temperature range. The analysis shows that the yield stress is a sensitive function of the temperature when the activation volume is small and that the temperature sensitivity is relatively weak at activation volumes V _f>10³ b ³. The theory predicts that the yield drop has a maximum at an intermediate temperature range and vanishes both at very low and very high temperatures. The calculated behavior was found to be in good qualitative agreement with the results of the measurements obtained by several investigators. It is concluded that the proposed theory can explain fully the temperature dependence of the yield behavior.     

Plane strain slip line theory for anisotropic rigid/plastic materials

A slip line theory governing states of incipient plane flow at the yield point load is developed for anisotropic rigid/plastic materials which exhibit a reduced yield criterion, governing states of plane flow, that depends only on the deviatoric parts of the in-plane stress tensor. It is shown that every homogeneous, incompressible material which complies with the principle of maximum plastic work, but which is of otherwise arbitrary anisotropy, is of this class. The stress equilibrium requirements are seen to take a remarkably simple form expressing the constancy of the quantities mean in-plane normal stress plus or minus arc length around the governing yield contour in a Mohr stress plane along members of the two slip line families. Further, this generalization of the Hencky equations is valid for every material of the considered class. Some special features of yield contours containing corners and flat segments are discussed, and velocity equations are given for materials complying with the maximum work inequality. The theory is applied to obtain the solution for indentation of an arbitrarily anisotropic half-space with a flat-ended, rigid, frictionless punch. A simple, universal formula, involving only geometrical dimensions of the governing yield contour, is derived for the yield point indentation pressure.     

Compressive response and failure of braided textile composites: Part 1—experiments

Experimental results obtained by examining the planar biaxial compression/tension response of carbon 2D triaxial braided composites (2DTBC) are reported in this paper. These experiments were motivated by a need to examine the failure of 2DTBC in a state of stress that would be similar to what is experienced by the walls of a tubular member under compressive crush loads. Results obtained from a series of biaxial tests that were conducted with different proportional displacement loading ratio combinations of compression and tension are reported. In all cases, the dominant failure mechanism under such a stress state is the buckling of the bias and axial tows within the composite. Full field surface displacement data is acquired concurrently during all biaxial and some uniaxial tests using the technique of digital speckle photography. Digital images of the specimen surface that is illuminated with a He-Ne laser are acquired at discrete time intervals during the loading history using a high-resolution digital camera. These images are stored and analyzed to obtain the incremental inplane surface displacement field, Δu(x,y) and Δv(x,y). From these, the incremental inplane surface strains Δε_x, Δε_y and Δγ_xy are obtained by numerical differentiation. The present paper, which is the first in a two part series, is devoted to the biaxial experimental results pertaining to 2DTBC failure.     

Thermodynamic driving force in ferroelectric crystals with a rank-2 laminated domain pattern, and a study of enhanced electrostriction

The thermodynamic driving force for domain growth in a rank-2 laminated ferroelectric crystal is derived in this article, and we used it, together with a homogenization theory, to study the issue of enhanced electrostrictive actuation recently reported by Burcsu et al. [2004. Large electrostrictive actuation of barium titanate single crystals. J. Mech. Phys. Solids 52, 823-846]. We derived this force from the reduction of Gibbs free energy with respect to the increase of domain concentration. It is shown that both the free energy and the thermodynamic force consist of three parts: the first arises from the difference in M₀ and M₁, the linear electromechanical compliances of the parent and product domains, respectively, at a given level of applied stress and electric field, the second stems from the electromechanical work associated with the change of spontaneous strain and spontaneous polarization during domain switch, and the third from the internal energy due to the distribution of polarizations strain and electric polarization inside the crystal. We prove that the first term is substantially lower than the second one, and the third one is identically zero with compatible domain pattern. The second one is, however, not exactly equal to the commonly written sum of the products of stress with strain, and electric field with polarization during switch, unless both domains have identical moduli in the common global axes. We also show that, with compatible domain patterns and when M₁=M₀, this driving force is identical to Eshelby's driving force acting on a flat interface due to the jump of energy-momentum tensor. Applications of the theory to a BaTiO₃ crystal subjected to a fixed axial compression and decreasing electric field from the [0 0 1] state reveal that the crystal undergoes a three-stage switching process: (i) the 0→90° switch to form a rank-1 laminate, (ii) the 0→180° switch inside the 0° domain to form a laminate I with a concurrent 90°→−90° switch inside the 90° domain to form laminate II, creating a rank-2-laminated domain pattern, and (iii) finally the 90→180° switch. It is the exchange of stability between the 0, 90°, and 180° domains under compression and electric field that is the origin of the enhanced actuation. We illustrate these intrinsic features by showing the evolution of these domains, and demonstrate how the reported large actuation strain can be attained with a rank-2 laminate.     

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

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

Measurement of stress distribution in sandwich beams under four-point bending

The bending-stress distributions through thickness in sandwich-composite beams are different from those obtained by conventional composite-beam theory because of the shear effect of the core, especially when the ratio of elastic moduli of face to core (k=E _f /E _c ) is large. Accordingly, the stress distributions in sandwich beams of composite materials with various combinations of face and core materials subjected to four-point bending are analyzed by introducing the multilayer-buitup theory. The bending stiffnesses of face and core, and the relative displacement between both faces are taken into consideration in the analysis. Photoelastic measurements were carried out on model specimens having four differentk-values and the applicable ranges of the two theories are discussed on the basis of the experimental results. It is shown that the experimental-stress distributions in sandwich-composite beams havingk-values larger than 120 can be well explained by the multilayer-builtup theory. The ratio of the coupling moment due to the axial forces in the two faces to the applied total moment, which denotes the sandwich structural efficiency, can be well estimated by the multilayer-builtup theory. The availability of this simple onedimensional theory should be useful in the structural design of sandwich beams with a small-core rigidity.     

Prediction of the initial thickness of shear band localization based on a reduced strain gradient theory

Plastic flow localization in ductile materials subjected to pure shear loading and uniaxial tension is investigated respectively in this paper using a reduced strain gradient theory, which consists of the couple-stress (CS) strain gradient theory proposed by Fleck and Hutchinson (1993) and the strain gradient hardening (softening) law (C–W) proposed by Chen and Wang (2000). Unlike the classical plasticity framework, the initial thickness of the shear band and the strain rate distribution in both cases are predicted analytically using a bifurcation analysis. It shows that the strain rate is obviously non-uniform inside the shear band and reaches a maximum at the center of the shear band. The initial thickness of the shear band depends on not only the material intrinsic length l_cs but also the material constants, such as the yield strength, ultimate tension strength, the linear hardening and softening shear moduli. Specially, in the uniaxial tension case, the most possible tilt angle of shear band localization is consistent qualitatively with the existing experimental observations. The results in this paper should be useful for engineers to predict the details of material failures due to plastic flow localization.     

Dispersion and Dispersivity Tensors in Saturated Porous Media with Uniaxial Symmetry

The coefficient of dispersion, D _ij , and the dispersivity, a _ijkl , appear in the expression for the flux of a solute in saturated flow through porous media. We present a detailed analysis of these tensors in an axially symmetric porous medium, e.g., a stratified porous medium, with alternating layers, and show that in such a medium, the dispersivity is governed by six independent moduli. We present also the constraints that have to be satisfied by these moduli. We also show that at least two independent experiments are required in order to obtain the values of these coefficients for any three-dimensional porous medium domain.     

Shakedown of unidirectional composites

The shakedown problem for a composite lamina made of an elastic-plastic matrix and elastic cylindrical fibers is studied. The plastic deformation modes of the lamina are reviewed, and it is concluded that significant shakedown effects can be caused only by the I₁ = 1/2(T₁₁ + T₂₂) and I₂ = T₃₃ components of the remotely applied stress field which are symmetric about the axis x₃ of the fiber; T₁₁ and T₂₂ are the normal composite stresses in the transverse plane. It is shown that the I₁I₂ stress system is needed also to represent thermal loads caused by a uniform change of temperature in the composite.Two methods for evaluation of shakedown limits in the I₁I₂-plane are described. First, the classical approach involving the determination of parametric families of self-stress fields and the solution of mathematical programming problems is used. Results are presented for selected B-Al, Be-Al, B-Ti and B-Mg composites.In the second method, the shakedown problem is related to the recently developed kinematic hardening rules for fibrous composites. It is shown that the composite will shake down for any loading program within a prescribed domain in the I₁I₂-plane, providing that the domain can be contained within a translated initial yield surface. This approach leads to a closed-form evaluation of shakedown limits for any arbitrary combination of mechanical and thermal variable cyclic loads in fibrous composites with temperaturedependent matrix yield strengths.The relationship between shakedown and fatigue in metal matrix composites is discussed.     

The shift of Curie temperature and evolution of ferroelectric domain in ferroelectric crystals

A micromechanics-based thermodynamic model for the phase transition of ferroelectric crystals is developed and, with it, the shift of Curie temperature and evolution of ferroelectric phase upon cooling are examined. This approach differs from the classical phenomenological one in that the evolution of new domain concentration can be predicted. We start out by formulating the Gibbs free energy of a generic, two-phase crystal consisting of the parent paraelectric phase and the transformed ferroelectric phase, at a given level of temperature, stress, and electric field. The thermodynamic driving force for domain growth is then derived and, together with the resistance force, a kinetic equation is established. The derived driving force is found to arise from three different sources of Gibbs free energy: (i) the interaction energy due to the heterogeneity of electromechanical moduli of the parent and product phases, (ii) the energy dissipation due to spontaneous polarization, and (iii) the self-energy of the dual-phase system due to the existence of polarization strain and electric polarization. For a BaTiO₃ crystal the electromechanical heterogeneity is found to play a rather significant role that seems not to have been recognized before. The derived shift recovers to the Clausius-Clapeyron relation if such heterogeneity disappears. We have examined in detail several factors that affect the shift of Curie temperature, and calculated the evolution of overall polarization and dielectric constant of a BaTiO₃ crystal. The results are found to be consistent with available test data.     

On the overall behavior, microstructure evolution, and macroscopic stability in reinforced rubbers at large deformations: I—Theory

This work presents an analytical framework for determining the overall constitutive response of elastomers that are reinforced by rigid or compliant fibers, and are subjected to finite deformations. The framework accounts for the evolution of the underlying microstructure, including particle rotation, which results from the finite changes in geometry that are induced by the applied loading. In turn, the evolution of the microstructure can have a significant geometric softening (or hardening) effect on the overall response, leading to the possible development of macroscopic instabilities through loss of strong ellipticity of the homogenized incremental moduli. The theory is based on a recently developed “second-order” homogenization method, which makes use of information on both the first and second moments of the fields in a suitably chosen “linear comparison composite,” and generates fairly explicit estimates—linearizing properly—for the large-deformation effective response of the reinforced elastomers. More specific applications of the results developed in this paper will be presented in Part II.     

Creep behavior of a poled PZT wafer under longitudinal tensile stress and through thickness electric field

A commercially available soft PZT wafer that is poled in thickness direction is subjected to three different sets of loading environment, and variations of electric displacement in thickness direction and longitudinal/transverse strains are measured over time. Pure tensile stress creep experiments are made in short and open-circuit conditions. Different material responses in the two electrical boundary conditions are explained by the effects of piezoelectrically produced internal electric field on linear material moduli and domain switching mechanisms. The material responses under pure antiparallel electric field load and combined load are compared with each other and the differences are explained by the opposite effects of longitudinal tensile stress on 90° domain switching.     

Theoretical and experimental studies in fracture mechanics of crack-like defects under biaxial loading

It is a common point of view in fracture mechanics that, for any geometry of the body with a crack and any boundary conditions for the loading acting in the body plane, the stress and displacement components near the crack tip can be approximated in the framework of the theory of elasticity by a one-parameter or one-term representation, i.e., strictly in terms of the stress intensity coefficients K _I and K _II for an arbitrary failure crack [1, 2]. The authors of [2] specified the Westergaard function of the singular solution for a central crack under the biaxial loading of a plate. This approximate two-component solution has satisfactory accuracy. It is clear from [2] that this method cannot be admitted as a general statement [1], although it has long been assumed to be correct. The cause is that one cannot reasonably justify neglecting the second term in the Williams representation of the stress components in the plane case in the form of eigenfunction series; the contribution of this term in the rectangular coordinate system x, y is independent of the distance from the crack tip. This method may result in a serious mistake, from both the qualitative and quantitative viewpoints, in the prediction of local stresses, displacements, and related variables that are of interest. Apparently, this can best be demonstrated by an example of biaxial loading of a plate with a crack [1]. The unfounded neglect of the second term (whose contribution is independent of the distance from the crack tip) in the series representing the stress components is the source of the above-mentioned difficulties. In this problem, the influence of the load applied in the direction parallel to the crack plane manifests itself only in the second term of the series [3]. Therefore, this term should be clearly determined and studied in detail in the case of technological welding defects (faulty fusions, incomplete fusions, undercuts, and slag inclusions) and crack-like defects (scratches and cuts) in the base metal. The influence of the stress σ _OX along the crack axis on the stress tensor σ _x , σ _y , τ _xy and on the displacements u _x and u _y is confirmed by experimental studies of cracks by the photoelasticity method [4].     

Exact solutions for radial deformations of a functionally graded isotropic and incompressible second-order elastic cylinder

We find closed-form solutions for axisymmetric plane strain deformations of a functionally graded circular cylinder comprised of an isotropic and incompressible second-order elastic material with moduli varying only in the radial direction. Cylinder's inner and outer surfaces are loaded by hydrostatic pressures. These solutions are specialized to cases where only one of the two surfaces is loaded. It is found that for a linear through-the-thickness variation of the elastic moduli, the hoop stress for the first-order solution (or in a cylinder comprised of a linear elastic material) is a constant but that for the second-order solution varies through the thickness. The radial displacement, the radial stress and the hoop stress do not depend upon the second-order elastic constant but the hydrostatic pressure and hence the axial stress depends upon it. When the two elastic moduli vary as the radius raised to the power two or four, the radial and the hoop stresses in an infinite space with a pressurized cylindrical cavity equal the pressure in the cavity. For an affine variation of the elastic moduli, the hoop stress in an internally loaded cylinder made of a linear elastic isotropic and incompressible material at the point is the same as that in a homogeneous cylinder. Here R_in and R_ou equal, respectively, the inner and the outer radius of the undeformed cylinder and R the radial coordinate of a point in the unstressed reference configuration.