Electromechanical influence of crack velocity at bifurcation for poled ferroelectric materials

The piezoelastodynamic field equations are solved to determine the crack velocity at bifurcation for poled ferroelectric materials where the applied electrical field and mechanical stress can be varied. The underlying physical mechanism, however, may not correspond to that assumed in the analytical model. Bifurcation has been related to the occurrence of a pair of maximum circumferential stress oriented symmetrically about the moving crack path. The velocity at which this behavior prevails has been referred to as the limiting crack speed. Unlike the classical approach, bifurcation will be identified with finite distances ahead of a moving crack. Nucleation of microcracks can thus be modelled in a single formulation. This can be accomplished by using the energy density function where fracture initiation is identified with dominance of dilatation in relation to distortion. Poled ferroelectric materials are selected for this study because the microstructure effects for this class of materials can be readily reflected by the elastic, piezoelectic and dielectric permittivity constants at the macroscopic scale. Existing test data could also shed light on the trend of the analytical predictions. Numerical results are thus computed for PZT-4 and compared with those for PZT-6B in an effort to show whether the branching behavior would be affected by the difference in the material microstructures. A range of crack bifurcation speed upsilon(b) is found for different r/a and E/sigma ratios. Here, r and a stand for the radial distance and half crack length, respectively, while E and a for the electric field and mechanical stress. For PZT-6B with upsilon(b) in the range 100-1700 m/s, the bifurcation angles varied from +/-6degrees to +/-39degrees. This corresponds to E/sigma of -0.072 to 0.024 V m/N. At the same distance r/a = 0.1, PZT-4 gives upsilon(b) values of 1100-2100 m/s; bifurcation angles of +/-15degrees to +/-49degrees; and E/sigma of -0.056 to 0.059 V m/N. In general, the bifurcation angles +/-theta(0) are found to decrease with decreasing crack velocity as the distance r/a is increased. Relatively speaking, the speed upsilon(b) and angles +/-theta(0) for PZT-4 are much greater than those for PZT-6B. This may be attributed to the high electromechanical coupling effect of PZT-4. Using upsilon(b)(0) as a base reference, an equality relation upsilon(b)(-) < upsilon(b)(0) < upsilon(b)(+) can be established. The superscripts -, 0 and + refer, respectively, to negative, zero and positive electric field. This is reminiscent of the enhancement and retardation of crack growth behavior due to change in poling direction. Bifurcation characteristics are found to be somewhat erratic when r/a approaches the range 10(-2)-10(-1) where the kinetic energy densities would fluctuate and then rise as the distance from the moving crack is increased. This is an artifact introduced by the far away condition of non-vanishing particle velocity. A finite kinetic energy density prevails at infinity unless it is made to vanish in the boundary value problem. Future works are recommended to further clarify the physical mechanism(s) associated with bifurcation by means of analysis and experiment. Damage at the microscopic level needs to be addressed since it has been known to affect the macrocrack speeds and bifurcation characteristics. (C) 2002 Published by Elsevier Science Ltd.

Torsional impact response of a penny-shaped interface crack in bonded materials with a graded material interlayer

In this paper, the dynamic response of a penny-shaped interface crack in bonded dissimilar homogeneous half-spaces is studied. It is assumed that the two materials are bonded together with such a inhomogeneous interlayer that makes the elastic modulus in the direction perpendicular to the crack surface is continuous throughout the space. The crack surfaces art assumed to be subjected to torsional impact loading. Laplace and Hankel integral transforms are applied combining with a dislocation density,function to reduce the mixed boundary value problem into a singular integral equation with a generalized Cauchy kernel in Laplace domain. By solving the singular integral equation numerically, and using a numerical Laplace inversion technique, the dynamic stress intensity factors art obtained. The influences of material properties and interlayer thickness on the dynamic stress intensity factor are investigated.

A constitutive model for cavitation and cavity growth in rubber-like materials under arbitrary tri-axial loading

In this paper, we attempted to construct a constitutive model to deal with the phenomenon of cavitation and cavity growth in a rubber-like material subjected to an arbitrary tri-axial loading. To this end, we considered a spherical elementary representative volume in a general Rivlin's incompressible material containing a central spherical cavity. The kinematics proposed by [Hou, H.S., Abeyaratne, R., 1992. Cavitation in elastic and elastic-plastic solids. J. Mech. Phys. Solids 40, 571-722] was adopted in order to construct an approximate but optimal field. In order to establish a suitable constitutive law for this class of materials, we utilized the homogenisation technique that permits us to calculate the average strain energy density of the volume. The cavity growth was considered through a physically realistic failure criterion. Combination of the constitutive law and the failure criterion enables us to describe correctly the global behaviour and the damage evolution of the material under tri-axial loading. It was shown that the present models can efficiently reproduce different stress states, varying from uniaxial to tri-axial tensions, observed in experimentations. Comparison between predicted results and experimental data proves that the proposed model is accurate and physically reasonable. Another advantage is that the proposed model does not need special identification work, the initial Rivlin's law for the corresponding incompressible material is sufficient to form the new law for the compressible material resulted from cavitation procedure. (C) 2007 Elsevier Ltd. All rights reserved.

Unloading characteristics of anti-plane shear crack in softening materials

The local characteristics of the anti-plane shear stress and strain field are determined for a material where the stress increases linearly with strain up to a limit and then softens nonlinearly. Two unloading models are considered such that the unloading path always returns to the origin while the other assumes the unloading modulus to be that of the initial shear modulus. As the applied shear increases, an unloading zone is found to prevail between a zone in which the material softens and another zone in which the material is linear-elastic even though the crack does not propagate. The divisions of these zones are displayed graphically.

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

In this paper, we present an exact higher-order asymptotic analysis on the near-crack-tip fields in elastic-plastic materials under plane strain, Mode I. A four- or five-term asymptotic series of the solutions is derived. It is found that when 1.6 < n less-than-or-equal-to 2.8 (here, n is the hardening exponent), the elastic effect enters the third-order stress field; but when 2.8< n less-than-or-equal-to 3.7 this effect turns to enter the fourth-order field, with the fifth-order field independent. Moreover, if n>3.7, the elasticity only affects the fields whose order is higher than 4. In this case, the fourth-order field remains independent. Our investigation also shows that as long as n is larger than 1.6, the third-order field is always not independent, whose amplitude coefficient K3 depends either on K1 or on both K1 and K2 (K1 and K2 arc the amplitude coefficients of the first- and second-order fields, respectively). Firmly, good agreement is found between our results and O'Dowd and Shih's numerical ones[8] by comparison.

Kinking of an interface crack between 2 dissimilar anisotropic elastic solids

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

Elastic-plastic finite element analyses of short cracks

Stress and strain distributions and crack opening displacement characteristics of short cracks have been studied in single edge notch bend and centre cracked panel specimens using elastic–plastic finite element analyses incorporating both a non strain hardening and a power law hardening behaviour. J contour integral solutions to describe stress strain conditions at crack tips for short cracks differ from those for long cracks. The analyses show that (i) short cracks can propagate at stress levels lower than those required for long cracks and (ii) a two-parameter description of crack tip fields is necessary for crack propagation.

Elastic and statistically brittle (ESB) model, damage localization and catastrophic rupture

In this paper, an elastic and statistically brittle (ESB) model is applied to the process of damage evolution induced catastrophic rupture and the influence of localization and softening on catastrophic rupture is discussed. According to the analysis, the uncertainty of catastrophic rupture should be attributed to the unknown scale of localized zone. Based on the elastic and statistically brittle model but local mean field approximation, the relation between the scale of localized zone and catastrophic rupture is obtained and then justified with experiments. These results can not only give a deeper understanding of the mechanism governing catastrophic rupture, but also provide a possible tool to foresee the occurrence of catastrophic rupture.

Adhesive behavior of two-dimensional power-law graded materials

In this paper, we investigate the adhesive contact between a rigid cylinder of radius R and a graded elastic half-space with a Young's modulus varying with depth according to a power-law, E = E-0(y/c(0))(k) (0 < k < 1), while the Poisson's ratio v remains constant. The results show that, for a given value of ratio R/C-0, a critical value of k exists at which the pull-off force attains a maximum; for a fixed value of k, the larger the ratio R/c(0), the larger the pull-off force is. For Gibson materials (i.e., k = 1 and v = 0.5), closed-form analytical solutions can be obtained for the critical contact half-width at pull-off and pull-off force. We further discuss the perfect stick case with both externally normal and tangential loads.

An Improved Ce/Se Scheme For Multi-Material Elastic-Plastic Flows And Its Applications

In this paper, a new definition of SE and CE, which is based on the hexahedron mesh and simpler than Chang's original CE/SE method (the space-time Conservation Element and Solution Element method), is proposed and an improved CE/SE scheme is constructed. Furthermore, the improved CE/SE scheme is extended in order to solve the elastic-plastic flow problems. The hybrid particle level set method is used for tracing the interfaces of materials. Proper boundary conditions are presented in interface tracking. Two high-velocity impact problems are simulated numerically and the computational results are carefully compared with the experimental data, as well as the results from other literature and LS-DYNA software. The comparisons show that the computational scheme developed currently is clear in physical concept, easy to be implemented and high accurate and efficient for the problems considered. (C) 2008 Elsevier Ltd. All rights reserved.

I. Rigid body penetration into brittle materials. II. Phase change effect on shock wave propagation

Part I.

We have developed a technique for measuring the depth time history of rigid body penetration into brittle materials (hard rocks and concretes) under a deceleration of ~ 10⁵ g. The technique includes bar-coded projectile, sabot-projectile separation, detection and recording systems. Because the technique can give very dense data on penetration depth time history, penetration velocity can be deduced. Error analysis shows that the technique has a small intrinsic error of ~ 3-4 % in time during penetration, and 0.3 to 0.7 mm in penetration depth. A series of 4140 steel projectile penetration into G-mixture mortar targets have been conducted using the Caltech 40 mm gas/ powder gun in the velocity range of 100 to 500 m/s.

We report, for the first time, the whole depth-time history of rigid body penetration into brittle materials (the G-mixture mortar) under 10⁵ g deceleration. Based on the experimental results, including penetration depth time history, damage of recovered target and projectile materials and theoretical analysis, we find:

1. Target materials are damaged via compacting in the region in front of a projectile and via brittle radial and lateral crack propagation in the region surrounding the penetration path. The results suggest that expected cracks in front of penetrators may be stopped by a comminuted region that is induced by wave propagation. Aggregate erosion on the projectile lateral surface is < 20% of the final penetration depth. This result suggests that the effect of lateral friction on the penetration process can be ignored.

2. Final penetration depth, P_max, is linearly scaled with initial projectile energy per unit cross-section area, e_s , when targets are intact after impact. Based on the experimental data on the mortar targets, the relation is P_max(mm) 1.15e_s (J/mm² ) + 16.39.

3. Estimation of the energy needed to create an unit penetration volume suggests that the average pressure acting on the target material during penetration is ~ 10 to 20 times higher than the unconfined strength of target materials under quasi-static loading, and 3 to 4 times higher than the possible highest pressure due to friction and material strength and its rate dependence. In addition, the experimental data show that the interaction between cracks and the target free surface significantly affects the penetration process.

4. Based on the fact that the penetration duration, t_max, increases slowly with e_s and does not depend on projectile radius approximately, the dependence of t_max on projectile length is suggested to be described by t_max(μs) = 2.08e_s (J/mm² + 349.0 x m/(πR²), in which m is the projectile mass in grams and R is the projectile radius in mm. The prediction from this relation is in reasonable agreement with the experimental data for different projectile lengths.

5. Deduced penetration velocity time histories suggest that whole penetration history is divided into three stages: (1) An initial stage in which the projectile velocity change is small due to very small contact area between the projectile and target materials; (2) A steady penetration stage in which projectile velocity continues to decrease smoothly; (3) A penetration stop stage in which projectile deceleration jumps up when velocities are close to a critical value of ~ 35 m/s.

6. Deduced averaged deceleration, a, in the steady penetration stage for projectiles with same dimensions is found to be a(g) = 192.4v + 1.89 x 10⁴, where v is initial projectile velocity in m/s. The average pressure acting on target materials during penetration is estimated to be very comparable to shock wave pressure.

7. A similarity of penetration process is found to be described by a relation between normalized penetration depth, P/P_max, and normalized penetration time, t/t_max, as P/P_max = f(t/t_max, where f is a function of t/t_max. After f(t/t_max is determined using experimental data for projectiles with 150 mm length, the penetration depth time history for projectiles with 100 mm length predicted by this relation is in good agreement with experimental data. This similarity also predicts that average deceleration increases with decreasing projectile length, that is verified by the experimental data.

8. Based on the penetration process analysis and the present data, a first principle model for rigid body penetration is suggested. The model incorporates the models for contact area between projectile and target materials, friction coefficient, penetration stop criterion, and normal stress on the projectile surface. The most important assumptions used in the model are: (1) The penetration process can be treated as a series of impact events, therefore, pressure normal to projectile surface is estimated using the Hugoniot relation of target material; (2) The necessary condition for penetration is that the pressure acting on target materials is not lower than the Hugoniot elastic limit; (3) The friction force on projectile lateral surface can be ignored due to cavitation during penetration. All the parameters involved in the model are determined based on independent experimental data. The penetration depth time histories predicted from the model are in good agreement with the experimental data.

9. Based on planar impact and previous quasi-static experimental data, the strain rate dependence of the mortar compressive strength is described by σ_f/σ⁰_f = exp(0.0905(log(έ/έ_0) ^1.14, in the strain rate range of 10^-7/s to 10³/s (σ⁰_f and έ are reference compressive strength and strain rate, respectively). The non-dispersive Hugoniot elastic wave in the G-mixture has an amplitude of ~ 0.14 GPa and a velocity of ~ 4.3 km/s.

Part II.

Stress wave profiles in vitreous GeO₂ were measured using piezoresistance gauges in the pressure range of 5 to 18 GPa under planar plate and spherical projectile impact. Experimental data show that the response of vitreous GeO₂ to planar shock loading can be divided into three stages: (1) A ramp elastic precursor has peak amplitude of 4 GPa and peak particle velocity of 333 m/s. Wave velocity decreases from initial longitudinal elastic wave velocity of 3.5 km/s to 2.9 km/s at 4 GPa; (2) A ramp wave with amplitude of 2.11 GPa follows the precursor when peak loading pressure is 8.4 GPa. Wave velocity drops to the value below bulk wave velocity in this stage; (3) A shock wave achieving final shock state forms when peak pressure is > 6 GPa. The Hugoniot relation is D = 0.917 + 1.711u (km/s) using present data and the data of Jackson and Ahrens [1979] when shock wave pressure is between 6 and 40 GPa for ρ₀ = 3.655 gj cm³ . Based on the present data, the phase change from 4-fold to 6-fold coordination of Ge⁺⁴ with O^-2 in vitreous GeO₂ occurs in the pressure range of 4 to 15 ± 1 GPa under planar shock loading. Comparison of the shock loading data for fused SiO₂ to that on vitreous GeO₂ demonstrates that transformation to the rutile structure in both media are similar. The Hugoniots of vitreous GeO₂ and fused SiO₂ are found to coincide approximately if pressure in fused SiO₂ is scaled by the ratio of fused SiO₂to vitreous GeO₂ density. This result, as well as the same structure, provides the basis for considering vitreous Ge0₂ as an analogous material to fused SiO₂ under shock loading. Experimental results from the spherical projectile impact demonstrate: (1) The supported elastic shock in fused SiO₂ decays less rapidly than a linear elastic wave when elastic wave stress amplitude is higher than 4 GPa. The supported elastic shock in vitreous GeO₂ decays faster than a linear elastic wave; (2) In vitreous GeO₂ , unsupported shock waves decays with peak pressure in the phase transition range (4-15 GPa) with propagation distance, x, as α 1/x^-3.35 , close to the prediction of Chen et al. [1998]. Based on a simple analysis on spherical wave propagation, we find that the different decay rates of a spherical elastic wave in fused SiO₂ and vitreous GeO₂ is predictable on the base of the compressibility variation with stress under one-dimensional strain condition in the two materials.