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.

Multiple Axial Cracks in A Coated Hollow Cylinder due to Thermal Shock

The transient thermal stress problem of an inner-surface-coated hollow cylinder with multiple pre-existing surface cracks contained in the coating is considered. The transient temperature, induced thermal stress, and the crack tip stress intensity factor (SIF) are calculated for the cylinder via finite element method (FEM), which is exposed to convective cooling from the inner surface. As an example, the material pair of a chromium coating and an underlying steel substrate 30CrNi2MoVA is particularly evaluated. Numerical results are obtained for the stress intensity factors as a function of normalized quantities such as time, crack length, convection severity, material constants and crack spacing. (c) 2005 Elsevier Ltd. All rights reserved.

The Flow Theory of Mechanism-Based Strain Gradient Plasticity

The flow theory of mechanism-based strain gradient (MSG) plasticity is established in this paper following the same multiscale, hierarchical framework for the deformation theory of MSG plasticity in order to connect with the Taylor model in dislocation mechanics. We have used the flow theory of MSG plasticity to study micro-indentation hardness experiments. The difference between deformation and flow theories is vanishingly small, and both agree well with experimental hardness data. We have also used the flow theory of MSG plasticity to investigate stress fields around a stationary mode-I crack tip as well as around a steady state, quasi-statically growing crack tip. At a distance to crack tip much larger than dislocation spacings such that continuum plasticity still applies, the stress level around a stationary crack tip in MSG plasticity is significantly higher than that in classical plasticity. The same conclusion is also established for a steady state, quasi-statically growing crack tip, though only the flow theory can be used because of unloading during crack propagation. This significant stress increase due to strain gradient effect provides a means to explain the experimentally observed cleavage fracture in ductile materials [J. Mater. Res. 9 (1994) 1734, Scripta Metall. Mater. 31 (1994) 1037; Interface Sci. 3(1996) 169].

Mechanism on Progressive Failure of a Faulted Rock Slope Due to Slip-Weakening

Slip-weakening is one of the characteristics of geological materials under certain loadings. Non-uniform rock structure may exist in the vicinity of the slip surface for a rock slope. Some portion of the slip surface may be penetrated but the other not. For the latter case, the crack or the fault surface will undergo shear deformation before it becomes a successive surface under a certain loading. As the slipped portion advances,slip-weakening occurs over a distance behind the crack tip. In the weakening zone, the shear strength will decrease from its peak value to residual friction level. The stress will redistribute along the surface of crack and in the weakening zone. Thus the changed local stress concentration leads the crack to extend and the ratio of penetration of the slip surface to increase. From the view of large-scale for the whole slip surface, the shear strength will decrease due to the damage of interior rock structure, and the faulted rock behaves as a softening material. Such a kind of mechanism performs in a large number of practical landslides in the zones experienced strong earthquakes. It should be noted that the mechanism mentioned above is different from that of the breakage of structural clay,in which the geological material is regarded as a medium containing structural lumps and structural bands. In this paper, the softening behavior of a faulted rock should be regarded as a comprehensive result of the whole complicated process including slip-weakening, redistribution of stress, extension of crack tip, and the penetration of the slip surface. This process is accompanied by progressive failure and abrupt structural damage. The size of slip-weakening zone is related to the undergoing strain. Once the relative slide is initiated (local or integrated), the effect of slip-weakening will behave in a certain length behind the crack tip until the formation of the whole slip surface.

Correlative reference model and molecular dynamics simulation of dislocation emission process

A correlative reference model for computer molecular dynamics simulations is proposed. Based on this model, a flexible displacement boundary scheme is introduced and the dislocations emitted from a crack tip can continuously pass through the border of the inner discrete atomic region and pile up at the outer continuum region. The effect of the emitted dislocations within the plastic zone on the inner atomistic region can be clearly demonstrated. The simulations for a molybdinum crystal show that a full dislocation in a bcc crystal is dissociated into three partial dislocations and interaction between the crack and the emitted dislocations results in gradual decrease of the local stress intensity factor.

A New Finite Element Method for Strain Gradient Theories and Applications to Fracture Analyses

A new compatible finite element method for strain gradient theories is presented. In the new finite element method, pure displacement derivatives are taken as the fundamental variables. The new numerical method is successfully used to analyze the simple strain gradient problems – the fundamental fracture problems. Through comparing the numerical solutions with the existed exact solutions, the effectiveness of the new finite element method is tested and confirmed. Additionally, an application of the Zienkiewicz–Taylor C1 finite element method to the strain gradient problem is discussed. By using the new finite element method, plane-strain mode I and mode II crack tip fields are calculated based on a constitutive law which is a simple generalization of the conventional J2 deformation plasticity theory to include strain gradient effects. Three new constitutive parameters enter to characterize the scale over which strain gradient effects become important. During the analysis the general compressible version of Fleck–Hutchinson strain gradient plasticity is adopted. Crack tip solutions, the traction distributions along the plane ahead of the crack tip are calculated. The solutions display the considerable elevation of traction within the zone near the crack tip.

Molecular dynamics simulation of hydrogen enhancing dislocation emission

The interactive pair potential between Al and H is obtained based on the ab initio calculation and the Chen-Mobius 3D lattice inversion formula. By utilizing the pair potentials calculated, the effects of hydrogen on the dislocation emission from crack tip have been studied. The simulated result shows that hydrogen can reduce the cohesive strength for Al single crystal, and then the critical stress intensity factor for partial dislocation emission decreases from 0.11 MPa root m (C-H = 0) to 0.075 MPa root m (C-H=0.72%) and 0.06 MPa root m (C-H = 1.44%). This indicates thar hydrogen can enhance the dislocation emission. The simulation also shows that atoms of hydrogen can gather and turn into small bubbles, resulting in enhancement of the equilibrium vacancy concentration.

压电材料平面问题的尖端场和应力强度因子的求解

利用Lekhnitskii理论和Stroh理论的相互联系,把已知的基于Lekhnitskii理论平面应变结果转化为Stroh理论形式的结果,直接获得Stroh公式中A,B的显式表达式,此方法可扩展到平面应力情况,然后导出压电材料平面应变问题的尖端场Williams形式的展开式,采用半权函数法计算有限大压电体平面问题应力和电位移强度因子.对无穷大板含中心裂纹的情况下本文结果和已有结果进行了比较,表明本文方法得到的结果精度可靠.本文方法的最大优点是可以求解有限压电体的应力强度因子,并且需要的单元少,精度高,实用性好.

Modeling Nonlinear Peeling of Ductile Thin Films-Critical Assessment of Analytical Bending Models Using FE Simulations

Three analytical double-parameter criteria based on a bending model and a two-dimensional finite element analysis model are presented for the modeling of ductile thin film undergoing a nonlinear peeling process. The bending model is based on different governing parameters: (1) the interfacial fracture toughness and the separation strength, (2) the interfacial fracture toughness and the crack tip slope angle, and (3) the interfacial fracture toughness and the critical Mises effective strain of the delaminated thin film at the crack tip. Thin film nonlinear peeling under steady-state condition is solved with the different governing parameters. In addition, the peeling test problem is simulated by using the elastic-plastic finite element analysis model. A critical assessment of the three analytical bending models is made by comparison of the bending model solutions with the finite element analysis model solutions. Furthermore, through analyses and comparisons for solutions based on both the bending model and the finite element analysis model, some connections between the bending model and the finite element analysis model are developed. Moreover, in the present research, the effect of different selections for cohesive zone shape on the ductile film peeling solutions is discussed.

Dynamic stress intensity factor of a functionally graded material under antiplane shear loading

The dynamic response of a finite crack in an unbounded Functionally Graded Material (FGM) subjected to an antiplane shear loading is studied in this paper. The variation of the shear modulus of the functionally graded material is modeled by a quadratic increase along the direction perpendicular to the crack surface. The dynamic stress intensity factor is extracted from the asymptotic expansion of the stresses around the crack tip in the Laplace transform plane and obtained in the time domain by a numerical Laplace inversion technique. The influence of graded material property on the dynamic intensity factor is investigated. It is observed that the magnitude of dynamic stress intensity factor for a finite crack in such a functionally graded material is less than in the homogeneous material with a property identical to that of the FGM crack plane.

Molecular dynamics simulation on dislocation emission process

A correlative reference model for a computer simulation of molecular dynamics is proposed in this paper. Based on this model, a flexible displacement boundary scheme is naturally introduced and the dislocations emitted from a crack tip are presumed to continuously pass through the border of an inner discrete atomic region to pile up at an outer continuum region. The simulations for a Mo crystal show that the interaction between a crack and emitted dislocations results in the decrease in local stress intensity factor gradually.

Assessment of Microscale Test Methods of Peeling and Splitting along Surface of Thin-Film/Substrates

Peel test methods are assessed through being applied to a peeling analysis of the ductile film/ceramic substrate system. Through computing the fracture work of the system using the either beam bend model (BB model) or the general plane analysis model (GPA model), surprisingly, a big difference between both model results is found. Although the BB model can capture the plastic dissipation phenomenon for the ductile film case as the GPA model can, it is much sensitive to the choice of the peeling criterion parameters, and it overestimates the plastic bending effect unable to capture crack tip constraint plasticity. In view of the difficulty of measuring interfacial toughness using peel test method when film is the ductile material, a new test method, split test, is recommended and analyzed using the GPA model. The prediction is applied to a wedge-loaded experiment for Al-alloy double-cantilever beam in literature.

Failure of SiC particulate-reinforced metal matrix composites induced by laser thermal shock

Thermal failure of SiC particulate-reinforced 6061 aluminum alloy composites induced by both laser thermal shock and mechanical load has been investigated. The specimens with a single-edge notch were mechanically polished to 0.25 mm in thickness. The notched-tip region of the specimen is subjected to laser beam rapid heating. In the test, a pulsed Nd:glass laser beam is used with duration 1.0 ms or 250 mu s, intensity 15 or 70 kW/cm(2), and spot size 5.0 mm in diameter. Threshold intensity was tested and fracture behavior was studied. The crack-tip process zone development and the microcrack formation were macroscopically and microscopically observed. It was found that in these materials, the initial crack occurred in the notched-tip region, wherein the initial crack was induced by either void nucleation, growth, and subsequent coalescence of the matrix materials or separation of the SiC particulate-matrix interface. It was further found that the process of the crack propagation occurred by the fracture of the SiC particulates.

The investigation on dynamic fracture behavour of materials under compressive-shear combined stress waves

In this paper, an improved plate impact experimental technique is presented for studying dynamic fracture mechanism of materials, under the conditions that the impacting loading is provided by a single pulse and the loading time is in the sub-microsecond range. The impacting tests are carried out on the pressure-shear gas gun. The loading rate achieved is dK/dt similar to 10(8) MPa m(1/2) s(-1). With the elimination of influence of the specimen boundary, the plane strain state of a semi-infinite crack in an infinite elastic plate is used to simulate the deformation fields of crack tip. The single pulses are obtained by using the "momentum trap" technique. Therefore, the one-time actions of the single pulse are achieved by eradicating the stress waves reflected from the specimen boundary or diffracted from the crack surfaces. In the current study, some important phenomena have been observed. The special loading of the single pulse can bring about material damage around crack tip, and affect the material behavior, such as kinking and branching of the crack propagation. Failure mode transitions from mode I to mode II crack are observed under asymmetrical impact conditions. The mechanisms of the dynamic crack propagation are consistent with the damage failure model.