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.

Damage analysis of tetragonal perovskite structure ceramics implicated by asymptotic field solutions and boundary conditions

Cracking of ceramics with tetragonal perovskite grain structure is known to appear at different sites and scale level. The multiscale character of damage depends on the combined effects of electromechanical coupling, prevailing physical parameters and boundary conditions. These detail features are exhibited by application of the energy density criterion with judicious use of the mode I asymptotic and full field solution in the range of r/a = 10(-4) to 10(-2) where r and a are, respectively, the distance to the crack tip and half crack length. Very close to the stationary crack tip, bifurcation is predicted resembling the dislocation emission behavior invoked in the molecular dynamics model. At the macroscopic scale, crack growth is predicted to occur straight ahead with two yield zones to the sides. A multiscale feature of crack tip damage is provided for the first time. Numerical values of the relative distances and bifurcation angles are reported for the PZT-4 ceramic subjected to different electric field to applied stress ratio and boundary conditions that consist of the specification of electric field/mechanical stress, electric displacement/mechanical strain, and mixed conditions. To be emphasized is that the multiscale character of damage in piezoceramics does not appear in general. It occurs only for specific combinations of the external and internal field parameters, elastic/piezoelectric/dielectric constants and specified boundary conditions. (C) 2002 Published by Elsevier Science Ltd.

First-passage failure of quasi-integrable Hamiltonian systems

The first-passage failure of quasi-integrable Hamiltonian si-stems (multidegree-of-freedom integrable Hamiltonian systems subject to light dampings and weakly random excitations) is investigated. The motion equations of such a system are first reduced to a set of averaged Ito stochastic differential equations by using the stochastic averaging method for quasi-integrable Hamiltonian systems. Then, a backward Kolmogorov equation governing the conditional reliability function and a set of generalized Pontryagin equations governing the conditional moments of first-passage time are established. Finally, the conditional reliability function, and the conditional probability density and moments of first-passage time are obtained by solving these equations with suitable initial and boundary conditions. Two examples are given to illustrate the proposed procedure and the results from digital simulation are obtained to verify the effectiveness of the procedure.

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.

Electrothermal stress in conductive body with collinear cracks

The temperature and stress field in a thin plate with collinear cracks interrupting an electric current field are determined. This is accomplished by using a complex function method that allows a direct means of finding the distribution of the electric current, the temperature and stress field. Temperature dependency for the heat-transfer coefficient, coefficient of linear expansion and the elastic modulus are considered. As an example, temperature distribution is calculated for an alloy (No. GH2132) plate with two collinear cracks under high temperature. Relationships between the stress, temperature, electric density and crack length are obtained. Crack trajectories emanating from existing crack are predicted by application of the strain energy density criterion which can also be used for finding the load carrying capacity of the cracked plate. (C) 2003 Elsevier Ltd. All rights reserved.

Semi-weight function method on computation of stress intensity factors in dissimilar materials

Semi-weight function method is developed to solve the plane problem of two bonded dissimilar materials containing a crack along the bond. From equilibrium equation, stress and strain relationship, conditions of continuity across interface and free crack surface, the stress and displacement fields were obtained. The eigenvalue of these fields is lambda. Semi-weight functions were obtained as virtual displacement and stress fields with eigenvalue-lambda. Integral expression of fracture parameters, K-I and K-II, were obtained from reciprocal work theorem with semi-weight functions and approximate displacement and stress values on any integral path around crack tip. The calculation results of applications show that the semi-weight function method is a simple, convenient and high precision calculation method.

Damage evaluation and damage localization of rock

Knowledge of damage accumulation and corresponding failure evolution are prerequisite for effective maintenance of civil engineering so as to avoid disaster. Based on statistical mesoscopic damage mechanics, it was revealed that there are three stages in the process of deformation, damage and failure of multiscale heterogeneous elastic-brittle medium. These are uniformly distributed damage, localized damage and catastrophic failure. In order to identify the transitions from scattering damage to macroscopically localized one, a condition for damage localization was given. The experiments of rock under uniaxial compression with the aid of observations of acoustic emission and speckle correlation do support the concept of localization. This provides a potential approach to properly evaluate damage accumulation in practice. In addition, it is found in the experiments that catastrophic failure displays critical sensitivity. This gives a helpful clue to the prediction of catastrophic failure. (C) 2004 Elsevier Ltd. All rights reserved.

Short-and resonant-range interactions between scales in turbulence and their applications

Interactions between different scales in turbulence were studied starting from the incompressible Navier-Stokes equations. The integral and differential formulae of the short-range viscous stresses, which express the short-range interactions between contiguous scales in turbulence, were given. A concept of the resonant-range interactions between extreme contiguous scales was introduced and the differential formula of the resonant-range viscous stresses was obtained. The short- and resonant-range viscous stresses were applied to deduce the large-eddy simulation (LES) equations as well as the multiscale equations, which are approximately closed and do not contain any empirical constants or relations. The properties and advantages of using the multiscale equations to compute turbulent flows were discussed. The short-range character of the interactions between the scales in turbulence means that the multiscale simulation is a very valuable technique for the calculation of turbulent flows. A few numerical examples were also given.

Observation and modeling for terrestrial processes in alpine meadow

The water-heat transfer process between land and atmosphere in Haibei alpine meadow area has been systematically observed. A multi-layer coupling model for land-atmosphere interaction was presented with special attention paid to the moisture transfer in leaf stomata under unsaturated condition. A profound investigation on the physical process of turbulent transfer inside the vegetation has been performed with a revised formula of water absorption for root system. The present model facilitates the study of vertically distributed physical variables in detail. Numerical simulation was conducted according to the transfer process of Kinesia humility meadow in the area of Haibei Alpine Meadow Ecosystem Station, CAS. The calculated results agree well with observation.

Coupling effects of heterogeneity and stress fluctuation on rupture

The coupling of mesoscopic strength distribution and stress fluctuation on damage evolution and rupture are examined. The numerical simulations show that there is only weak stress fluctuation at the initial damage stage when the mean field approximation is in effect. As the damage fraction becomes larger than the threshold value, the fluctuation is amplified significantly, and damage localization appears. The coupling between stress fluctuation, disordered heterogeneity and the damage localization may play an essential role in catastrophic rupture. (C) 2003 Elsevier Ltd. All rights reserved.

Lattice type of fracture model for concrete

Concrete is usually described as a three-phase material, where matrix, aggregate and interface zones are distinguished. The beam lattice model has been applied widely by many investigators to simulate fracture processes in concrete. Due to the extremely large computational effort, however, the beam lattice model faces practical difficulties. In our investigation, a new lattice called generalized beam (GB) lattice is developed to reduce computational effort. Numerical experiments conducted on a panel subjected to uniaxial tension show that the GB lattice model can reproduce the load-displacement curves and crack patterns in agreement to what are observed in tests. Moreover, the effects of the particle overlay on the fracture process are discussed in detail. (C) 2007 Elsevier Ltd. All rights reserved.

A moving-coil designed micro-mechanics tester with application on MEMS

A moving-coil designed micro-mechanics tester, named as MicroUTM (universal testing machine), is in-house developed in this paper for micro-mechanics tests. The main component is a moving coil suspended in a uniform magnetic field through a set of springs. When a current passes through the coil, the electromagnetic force is proportional to the magnitude of the current, so the load can easily be measured by the current. The displacement is measured using a capacitive sensor. The load is calibrated using a Sartorius BP211D analytical balance, with a resolution/range of 0.01 mg/80 g or 0.1 mg/210 g. The displacement is calibrated using a HEIDENHAIN CT-6002 length gauge with an accuracy of +/- 0.1 mu m. The calibration results show that the load range is +/- 1 N and the displacement range is +/- 300 mu m. The noise levels of the load and displacement are 50 mu N and 150 nm, respectively. The nonlinearity of the load is only 0.2%. Several in-plane load tests of the MEMS micro-cantilever are performed using this tester. Experimental results, with excellent repeatability, demonstrate the reliability of the load measurement as well as the flexible function of this tester.

Dynamic equations for curved submerged floating tunnel

In virtue of reference Cartesian coordinates, geometrical relations of spatial curved structure are presented in orthogonal curvilinear coordinates. Dynamic equations for helical girder are derived by Hamilton principle. These equations indicate that four generalized displacements are coupled with each other. When spatial structure degenerates into planar curvilinear structure, two generalized displacements in two perpendicular planes are coupled with each other. Dynamic equations for arbitrary curvilinear structure may be obtained by the method used in this paper.

Relationship between strain localization and catastrophic rupture

In order to explore a prior warning to catastrophic rupture of heterogeneous media, like rocks, the present study investigates the relationship between surface strain localization and catastrophic rupture. Instrumented observations on the evolution of surface strain field and the catastrophic rupture of a rock under uniaxial compression were carried out. It is found that the evolution of surface strain field displays two phases: at the early stage, the strain field keeps nearly uniform with weak fluctuations increasing slowly; but at the stage prior to catastrophic rupture, a certain accelerating localization develops and a localized zone emerges. Based on the measurements, an analysis was performed with local mean-field approximation. More importantly, it is found that the scale of localized zone is closely related to the catastrophic rupture strain and the rupture strain can be calculated in accord with the local-mean-field model satisfactorily. This provides a possible clue to the forecast of catastrophic rupture. (c) 2007 Elsevier Ltd. All rights reserved.