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.

A field model interpretation of crack initiation and growth behavior in ferroelectric ceramics: change of poling direction and boundary condition

Strain energy density expressions are obtained from a field model that can qualitatively exhibit how the electrical and mechanical disturbances would affect the crack growth behavior in ferroelectric ceramics. Simplification is achieved by considering only three material constants to account for elastic, piezoelectric and dielectric effects. Cross interaction of electric field (or displacement) with mechanical stress (or strain) is identified with the piezoelectric effect; it occurs only when the pole is aligned normal to the crack. Switching of the pole axis by 90degrees and 180degrees is examined for possible connection with domain switching. Opposing crack growth behavior can be obtained when the specification of mechanical stress sigma(infinity) and electric field E-infinity or (sigma(infinity), E-infinity) is replaced by strain e and electric displacement D-infinity or (epsilon(infinity), D-infinity). Mixed conditions (sigma(infinity),D-infinity) and (epsilon(infinity),E-infinity) are also considered. In general, crack growth is found to be larger when compared to that without the application of electric disturbances. This includes both the electric field and displacement. For the eight possible boundary conditions, crack growth retardation is identified only with (E-y(infinity),sigma(y)(infinity)) for negative E-y(infinity) and (D-y(infinity), epsilon(y)(infinity)) for positive D-y(infinity) while the mechanical conditions sigma(y)(infinity) or epsilon(y)infinity are not changed. Suitable combinations of the elastic, piezoelectric and dielectric material constants could also be made to suppress crack growth. (C) 2002 Published by Elsevier Science Ltd.

Nonlinear fatigue damage model based on the residual strength degradation law

In this paper, a logarithmic expression to describe the residual strength degradation process is developed in order to fatigue test results for normalized carbon steel. The definition and expression of fatigue damage due to symmetrical stress with a constant amplitude are also given. The expression of fatigue damage can also explain the nonlinear properties of fatigue damage. Furthermore, the fatigue damage of structures under random stress is analyzed, and an iterative formula to describe the fatigue damage process is deduced. Finally, an approximate method for evaluating the fatigue life of structures under repeated random stress blocking is presented through various calculation examples.

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.

Micromechanics analysis on evolution of crack in viscoelastic materials

A preliminary analysis on crack evolution in viscoelastic materials was presented. Based on the equivalent inclusion concept of micro-mechanics theory, the explicit expressions of crack opening displacement delta and energy release rate G were derived, indicating that both delta and G are increasing with time. The equivalent modulus of the viscoelastic solid comprising cracks was evaluated. It is proved that the decrease of the modulus comes from two mechanisms: one is the viscoelasticity of the material; the other is the crack opening which is getting larger with time.

Experimental study of gaseous detonation propagation over acoustically absorbing walls

The paper presents results front an experimental investigation of the propagation of gaseous detonation waves over tube sections lined with acoustically absorbent materials. The measurements were compared with results from control tests in a smooth wall section. The results show the increasing effectiveness of a perforated steel plate, wire mesh and steel wool in attenuating detonation.

Singularity of dynamic stress fields around an interface crack between viscoelastic bodies

The singular nature of the dynamic stress fields around an interface crack located between two dissimilar isotropic linearly viscoelastic bodies is studied. A harmonic load is imposed on the surfaces of the interface crack. The dynamic stress fields around the crack are obtained by solving a set of simultaneous singular integral equations in terms of the normal and tangent crack dislocation densities. The singularity of the dynamic stress fields near the crack tips is embodied in the fundamental solutions of the singular integral equations. The investigation of the fundamental solutions indicates that the singularity and oscillation indices of the stress fields are both dependent upon the material constants and the frequency of the harmonic load. This observation is different from the well-known -1/2 oscillating singularity for elastic bi-materials. The explanation for the differences between viscoelastic and elastic bi-materials can be given by the additional viscosity mismatch in the case of viscoelastic bi-materials. As an example, the standard linear solid model of a viscoelastic material is used. The effects of the frequency and the material constants (short-term modulus, long-term modulus and relaxation time) on the singularity and the oscillation indices are studied numerically.

An investigation of electromagnetic wave propagation in plasma by shock tube

This paper presents the electromagnetic wave propagation characteristics in plasma and the attenuation coefficients of the microwave in terms of the parameters n(e), v, w, L, w(b). The phi800 mm high temperature shock tube has been used to produce a uniform plasma. In order to get the attenuation of the electromagnetic wave through the plasma behind a shock wave, the microwave transmission has been used to measure the relative change of the wave power. The working frequency is f = (2 similar to 35) GHz (w = 2pif, wave length lambda = 15 cm similar to 8 mm). The electron density in the plasma is n(e) = (3 x 10(10) similar to 1 x 10(14)) cm(-3). The collision frequency v = (1 x 10(8) similar to 6 x 10(10)) Hz. The thickness of the plasma layer L = (2 similar to 80) cm. The electron circular frequency w(b) = eB(0)/m(e), magnetic flux density B-0 = (0 similar to 0.84) T. The experimental results show that when the plasma layer is thick (such as L/lambda greater than or equal to 10), the correlation between the attenuation coefficients of the electromagnetic waves and the parameters n(e), v, w, L determined from the measurements are in good agreement with the theoretical predictions of electromagnetic wave propagations in the uniform infinite plasma. When the plasma layer is thin (such as when both L and lambda are of the same order), the theoretical results are only in a qualitative agreement with the experimental observations in the present parameter range, but the formula of the electromagnetic wave propagation theory in an uniform infinite plasma can not be used for quantitative computations of the correlation between the attenuation coefficients and the parameters n(e), v, w, L. In fact, if w < w(p), v(2) much less than w(2), the power attenuations K of the electromagnetic waves obtained from the measurements in the thin-layer plasma are much smaller than those of the theoretical predictions. On the other hand, if w > w(p), v(2) much less than w(2) (just v approximate to f), the measurements are much larger than the theoretical results. Also, we have measured the electromagnetic wave power attenuation value under the magnetic field and without a magnetic field. The result indicates that the value measured under the magnetic field shows a distinct improvement.

Atomistic simulations of crack nucleation and intergranular fracture in bulk nanocrystalline nickel

We report large scale molecular dynamics simulations of dynamic cyclic uniaxial tensile deformation of pure, fully dense nanocrystalline Ni, to reveal the crack initiation, and consequently intergranular fracture is the result of coalescence of nanovoids by breaking atomic bonds at grain boundaries and triple junctions. The results indicate that the brittle fracture behavior accounts for the transition from plastic deformation governed by dislocation to one that is grain-boundary dominant when the grain size reduces to the nanoscale. The grain-boundary mediated plasticity is also manifested by the new grain formation and growth induced by stress-assisted grain-boundary diffusion observed in this work. This work illustrates that grain-boundary decohesion is one of the fundamental deformation mechanisms in nanocrystalline Ni.

Wave propagation and the frequency domain Green's functions in viscoelastic Biot/squirt (BISQ) media

In this paper, we examine the characteristics of elastic wave propagation in viscoelastic porous media, which contain simultaneously both the Biot-flow and the squirt-flow mechanisms (BISQ). The frequency-domain Green's functions for viscoelastic BISQ media are then derived based on the classic potential function methods. Our numerical results show that S-waves are only affected by viscoelasticity, but not by squirt-flows. However, the phase velocity and attenuation of fast P-waves are seriously influenced by both viscoelasticity and squirt-flows; and there exist two peaks in the attenuation-frequency variations of fast P-waves. In the low-frequency range, the squirt-flow characteristic length, not viscoelasticity, affects the phase velocity of slow P-waves, whereas it is opposite in the high-frequency range. As to the contribution of potential functions of two types of compressional waves to the Green's function, the squirt-flow length has a small effect, and the effects of viscoelastic parameter are mainly in the higher frequency range. Crown Copyright (C) 2006 Published by Elsevier Ltd. All rights reserved.