Fracture mechanics of piezoelectric materials

This paper presents an analysis of crack problems in homogeneous piezoelectrics or on the interfaces between two dissimilar piezoelectric materials based on the continuity of normal electric displacement and electric potential across the crack faces. The explicit analytic solutions are obtained for a single crack in an infinite piezoelectric or on the interface of piezoelectric bimaterials. For homogeneous materials it is found that the normal electric displacement D-2, induced by the crack, is constant along the crack faces which depends only on the remote applied stress fields. Within the crack slit, the perturbed electric fields induced by the crack are also constant and not affected by the applied electric displacement fields. For bimaterials, generally speaking, an interface crack exhibits oscillatory behavior and the normal electric displacement D-2 is a complex function along the crack faces. However, for bimaterials, having certain symmetry, in which an interface crack displays no oscillatory behavior, it is observed that the normal electric displacement D-2 is also constant along the crack faces and the electric field E-2 has the singularity ahead of the crack tip and has a jump across the interface. Energy release rates are established for homogeneous materials and bimaterials having certain symmetry. Both the crack front parallel to the poling axis and perpendicular to the poling axis are discussed. It is revealed that the energy release rates are always positive for stable materials and the applied electric displacements have no contribution to the energy release rates.

Dislocation theory of the fracture criterion for anisotropic solids

A dislocation theory of fracture criterion for the mixed dislocation emission and cleavage process in an anisotropic solid is developed in this paper. The complicated cases involving mixed-mode loading are considered here. The explicit formula for dislocations interaction with a semi-infinite crack is obtained. The governing equation for the critical condition of crack cleavage in an anisotropic solid after a number dislocation emissions is established. The effects of elastic anisotropy, crack geometry and load phase angle on the critical energy release rate and the total number of the emitted dislocations at the onset of cleavage are analysed in detail. The analyses revealed that the critical energy release rates can increase to one or two magnitudes larger than the surface energy because of the dislocation emission. It is also found elastic anisotropy and crystal orientation have significant effects on the critical energy release rates. The anisotropic values can be several times the isotropic value in one crack orientation. The values may be as much as 40% less than the isotropic value in another crack orientation and another anisotropy parameter. Then the theory is applied to a fee single crystal. An edge dislocation can emit from the crack tip along the most highly shear stressed slip plane. Crack cleavage can occur along the most highly stressed slip plane after a number of dislocation emissions. Calculation is carried out step by step. Each step we should judge by which slip system is the most highly shear stressed slip system and which slip system has the largest energy release rate. The calculation clearly shows that the crack orientation and the load phase angle have significant effects on the crystal brittle-ductile behaviours.

Analysis of strip electric saturation model of crack problem in piezoelectric materials

This paper presents a fully anisotropic analysis of strip electric saturation model proposed by Gao et al. (1997) (Gao, H.J., Zhang, T.Y., Tong, P., 1997. Local and global energy release rates for an electrically yielded crack in a piezoelectric ceramic. J. Mech. Phys. Solids, 45, 491-510) for piezoelectric materials. The relationship between the size of the strip saturation zone ahead of a crack tip and the applied electric displacement field is established. It is revealed that the critical fracture stresses for a crack perpendicular to the poling axis is linearly decreased with the increase of the positive applied electric field and increases linearly with the increase of the negative applied electric field. For a crack parallel to the poring axis, the failure stress is not effected by the parallel applied electric field. In order to analyse the existed experimental results, the stress fields ahead of the tip of an elliptic notch in an infinite piezoelectric solid are calculated. The critical maximum stress criterion is adopted for determining the fracture stresses under different remote electric displacement fields. The present analysis indicates that the crack initiation and propagation from the tip of a sharp elliptic notch could be aided or impeded by an electric displacement field depending on the field direction. The fracture stress predicted by the present analysis is consistent with the experimental data given by Park and Sun (1995) (Park, S., Sun, C.T., 1995. Fracture criteria for piezoelectric materials. J. Am. Ceram. Soc 78, 1475-1480).

The Displacement and Strain Field of Three-Dimensional Rheologic Model of Earthquake Preparation

Based on the three-dimensional elastic inclusion model proposed by Dobrovolskii, we developed a rheological inclusion model to study earthquake preparation processes. By using the Corresponding Principle in the theory of rheologic mechanics, we derived the analytic expressions of viscoelastic displacement U(r, t) , V(r, t) and W(r, t), normal strains epsilon(xx) (r, t), epsilon(yy) (r, t) and epsilon(zz) (r, t) and the bulk strain theta (r, t) at an arbitrary point (x, y, z) in three directions of X axis, Y axis and Z axis produced by a three-dimensional inclusion in the semi-infinite rheologic medium defined by the standard linear rheologic model. Subsequent to the spatial-temporal variation of bulk strain being computed on the ground produced by such a spherical rheologic inclusion, interesting results are obtained, suggesting that the bulk strain produced by a hard inclusion change with time according to three stages (alpha, beta, gamma) with different characteristics, similar to that of geodetic deformation observations, but different with the results of a soft inclusion. These theoretical results can be used to explain the characteristics of spatial-temporal evolution, patterns, quadrant-distribution of earthquake precursors, the changeability, spontaneity and complexity of short-term and imminent-term precursors. It offers a theoretical base to build physical models for earthquake precursors and to predict the earthquakes.

Formation, Thermal Stability and Deformation Behavior of Graphite-Flakes Reinforced Cu-based Bulk Metallic Glass Matrix Composites

Graphite-flake reinforced Cu47Ti34Zr11 Ni-8 bulk metallic glass matrix composite was fabricated by water-cooled copper mould cast. Most of the graphite flakes still keep unreacted and distribute uniformly in the amorphous matrix except that some reactive wetting occurs by the formation of TiC particles around the flakes. It reveals that the presence of graphite flakes does not affect the onset of the glass transition temperature, crystallization reaction and liquidus of the metallic glass. The resulting material shows obvious serrated flow and higher fracture strength under room temperature compressive load, comparing with the monolithic bulk metallic glass (BMG). Three types of interaction between the shear bands and graphite flakes, namely, shear band termination, shear bands branching and new shear bands formation near the graphite flakes can be observed by quasi-static uniaxial compression test and bonded interface technique through Vickers indentation.

Transonic Aeroelastic Numerical Simulation in Aeronautical Engineering

A lower-upper symmetric Gauss-Seidel (LU-SGS) subiteration scheme is constructed for time-marching of the fluid equations. The Harten-Lax-van Leer-Einfeldt-Wada (HLLEW) scheme is used for the spatial discretization. The same subiteration formulation is applied directly to the structural equations of motion in generalized coordinates. Through subiteration between the fluid and structural equations, a fully implicit aeroelastic solver is obtained for the numerical simulation of fluid/structure interaction. To improve the ability for application to complex configurations, a multiblock grid is used for the flow field calculation and transfinite interpolation (TFI) is employed for the adaptive moving grid deformation. The infinite plate spline (IPS) and the principal of virtual work are utilized for the data transformation between the fluid and structure. The developed code was first validated through the comparison of experimental and computational results for the AGARD 445.6 standard aeroelastic wing. Then, the flutter character of a tail wing with control surface was analyzed. Finally, flutter boundaries of a complex aircraft configuration were predicted.

Stress Wave Propagation in a Gradient Elastic Medium

The gradient elastic constitutive equation incorporating the second gradient of the strains is used to determine the monochromatic elastic plane wave propagation in a gradient infinite medium and thin rod. The equation of motion, together with the internal material length, has been derived. Various dispersion relations have been determined. We present explicit expressions for the relationship between various wave speeds, wavenumber and internal material length.

A Crack Perpendicular to the Bimaterial Interface in Finite Solid

The dislocation simulation method is used in this paper to derive the basic equations for a crack perpendicular to the bimaterial interface in a finite solid. The complete solutions to the problem, including the T stress and the stress intensity factors are obtained. The stress field characteristics are investigated in detail. It is found that when the crack is within a weaker material, the stress intensity factor is smaller than that in a homogeneous material and it decreases when the distance between the crack tip and interface decreases. When the crack is within a stiffer material, the stress intensity factor is larger than that in a homogeneous material and it increases when the distance between the crack tip and interface decreases. In both cases, the stress intensity factor will increase when the ratio of the size of a sample to the crack length decreases. A comparison of stress intensity factors between a finite problem and an infinite problem has been given also. The stress distribution ahead of the crack tip, which is near the interface, is shown in details and the T stress effect is considered.

Singularity characteristics for a lip-shaped crack subjected to remote biaxial loading

In this paper, the conformal mapping method was adopted to solve the problem of an infinite plate containing a central lip-shaped crack subjected to remote biaxial loading. A kind of leaf-shaped configuration was also constructed in order to solve the problem. The analytical result showed that the singularity order of the stress field at the tip of a lip-shaped crack remains -1/2, despite the difference in notch-crack width.

Saltation in wind blown sand

The Boltzmann equation of the sand particle velocity distribution function in wind-blown sand two-phase flow is established based on the motion equation of single particle in air. And then, the generalized balance law of particle property in single phase granular flow is extended to gas-particle two-phase flow. The velocity distribution function of particle phase is expanded into an infinite series by means of Grad's method and the Gauss distribution is used to replace Maxwell distribution. In the case of truncation at the third-order terms, a closed third-order moment dynamical equation system is constructed. The theory is further simplified according to the measurement results obtained by stroboscopic photography in wind tunnel tests.

Dynamic SIF of interface crack between two dissimilar viscoelastic bodies under impact loading

In this paper, the transient dynamic stress intensity factor (SIF) is determined for an interface crack between two dissimilar half-infinite isotropic viscoelastic bodies under impact loading. An anti-plane step loading is assumed to act suddenly on the surface of interface crack of finite length. The stress field incurred near the crack tip is analyzed. The integral transformation method and singular integral equation approach are used to get the solution. By virtue of the integral transformation method, the viscoelastic mixed boundary problem is reduced to a set of dual integral equations of crack open displacement function in the transformation domain. The dual integral equations can be further transformed into the first kind of Cauchy-type singular integral equation (SIE) by introduction of crack dislocation density function. A piecewise continuous function approach is adopted to get the numerical solution of SIE. Finally, numerical inverse integral transformation is performed and the dynamic SIF in transformation domain is recovered to that in time domain. The dynamic SIF during a small time-interval is evaluated, and the effects of the viscoelastic material parameters on dynamic SIF are analyzed.

Effective elastic moduli of two-dimensional solids with multiple cracks

An accurate method which directly accounts for the interactions between different microcracks is used for analyzing the elastic problem of multiple cracks solids. The effective elastic moduli for randomly oriented cracks and parallel cracks are evaluated for the representative volume element (RVE) with microcracks in infinite media. The numerical results are compared with those from various micromechanics models and experimental data. These results show that the present method is simple and provides a direct and efficient approach to dealing with elastic solids containing multiple cracks.

Diffusive Wave Solutions for Open Channel Flows with Uniform and Concentrated Lateral Inflow

The diffusive wave equation with inhomogeneous terms representing hydraulics with uniform or concentrated lateral inflow intoa river is theoretically investigated in the current paper. All the solutions have been systematically expressed in a unified form interms of response function or so called K-function. The integration of K-function obtained by using Laplace transform becomesS-function, which is examined in detail to improve the understanding of flood routing characters. The backwater effects usuallyresulting in the discharge reductions and water surface elevations upstream due to both the downstream boundary and lateral infloware analyzed. With a pulse discharge in upstream boundary inflow, downstream boundary outflow and lateral inflow respectively,hydrographs of a channel are routed by using the S-functions. Moreover, the comparisons of hydrographs in infinite, semi-infiniteand finite channels are pursued to exhibit the different backwater effects due to a concentrated lateral inflow for various channeltypes.

Behavior of Multiple Shear Bands in Zr-Based Bulk Metallic Glass

The deformation behavior of Zr41.2Ti13.8Cu12.5Ni10Be22.5 bulk metallic glass was studied by in situ scanning electron microscopy (SEM) quasi-static uniaxial compression tests at room temperature. Multiple shear bands were observed with a large plasticity. Microscopic examination demonstrates that slipping, branching and intersecting of multiple shear bands are the main mechanisms for enhancing the plasticity of this metallic glass. Additionally, nano/micro-scale voids and cracks at the intersecting sites of shear bands and preferential etching of shear bands were observed as well. These observations demonstrated that the formation of shear bands in bulk metallic glasses is resulted mainly from local free volume coalescence.