Some advances in study of crack identification problems

In the present paper, the crack identification problems are investigated. This kind of problems belong to the scope of inverse problems and are usually ill-posed on their solutions. The paper includes two parts: (1) Based on the dynamic BIEM and the optimization method and using the measured dynamic information on outer boundary, the identification of crack in a finite domain is investigated and a method for choosing the high sensitive frequency region is proposed successfully to improve the precision. (2) Based on 3-D static BIEM and hypersingular integral equation theory, the penny crack identification in a finite body is reduced to an optimization problem. The investigation gives us some initial understanding on the 3-D inverse problems.

Finite Boundary Effects in Problem of a Crack Perpendicular to and Terminating at a Bimaterial Interface

In this paper, the problem of a crack perpendicular to and terminating at an interface in bimaterial structure with finite boundaries is investigated. The dislocation simulation method and boundary collocation approach are used to derive and solve the basic equations. Two kinds of loading form are considered when the crack lies in a softer or a stiffer material, one is an ideal loading and the other one fits to the practical experiment loading. Complete solutions of the stress field including the T stress are obtained as well as the stress intensity factors. Influences of T stress on the stress field ahead of the crack tip are studied. Finite boundary effects on the stress intensity factors are emphasized. Comparisons with the problem presented by Chen et al. (Int. J. Solids and Structure, 2003, 40, 2731-2755) are discussed also.

Stress state in front of a crack perpendicular to bimaterial interface

Using a dislocation simulation approach, the basic equation for a crack perpendicular to a bimaterial interface is formulated in this paper. A novel expansion method is proposed for solving the problem. The complete solution for the problem, including the T stress ahead of the crack tip and the stress intensity factors are presented. The stress field characteristics are analyzed in detail. It is found that ahead of the crack tip and near the interface the normal stress, perpendicular to the crack plane, sigma(x), is characterized by the K fields and the normal stress sigma(y) is dominated by the K field plus T stress in the region of 0 < r/b < 0.4 for b/a(0) less than or equal to 0.1, where b is the distance from the crack tip to the interface.

Steady-state crack growth and fracture work based on the theory of mechanism-based strain gradient plasticity

Mode I steady-state crack growth is analyzed under plane strain conditions in small scale yielding. The elastic-plastic solid is characterized by the mechanism-based strain gradient (MSG) plasticity theory [J. Mech. Phys. Solids 47 (1999) 1239, J. Mech. Phys. Solids 48 (2000) 99]. The distributions of the normal separation stress and the effective stress along the plane ahead of the crack tip are computed using a special finite element method based on the steady-state fundamental relations and the MSG flow theory. The results show that during the steady-state crack growth, the normal separation stress on the plane ahead of the crack tip can achieve considerably high value within the MSG strain gradient sensitive zone. The results also show that the crack tip fields are insensitive to the cell size parameter in the MSG theory. Moreover, in the present research, the steady-state fracture toughness is computed by adopting the embedded process zone (EPZ) model. The results display that the steady-state fracture toughness strongly depends on the separation strength parameter of the EPZ model and the length scale parameter in the MSG theory. Furthermore, in order for the results of steady crack growth to be comparable, an approximate relation between the length scale parameters in the MSG theory and in the Fleck-Hutchinson strain gradient plasticity theory is obtained.

Instability Of Crack Propagation In Brittle Bulk Metallic Glass

The instability of the crack tip in brittle Mg-based bulk metallic glass (BMG) is studied. The formation of various fractographic surfaces of the BMG is associated with the instability of the fluid meniscus, which is due to viscous fluid matter being present on the fracture process zone. Depending on the values of the wavelength of the initial perturbation of the fluid meniscus and the local stress intensity factor, different fracture surface profiles, i.e. a dimple-like structure, a periodic corrugation pattern and a pure mirror zone are formed. The fractographic evolution is significantly affected by the applied stress. A decreased fracture Surface roughness is observed under a low applied stress. An increased fracture surface roughness, which has frequently been reported by other researchers, is also observed in the present studies under a high applied stress. Unique fractographic features are attributed to the non-linear hyperelastic stiffening for less softening) mechanism. (C) 2008 Acta Materialia Inc. Published by Elsevier Ltd. All rights reserved.

Inverse Grain-Size Effect On Twinning In Nanocrystalline Ni

A long-standing controversy exists between molecular dynamics simulations and experiments on the twinning propensity of nanocrystalline (NC) face-centered-cubic metals. For example, three-dimensional molecular dynamics simulations rarely observed twins in NC Ni, whereas experiments readily observed them. Here this discrepancy is resolved by experimental observation of an inverse grain-size effect on twinning. Specifically, decreasing the grain size first promotes twinning in NC Ni and then hinders twinning due to the inverse grain-size effect. Interestingly, no inverse grain-size effect exists on stacking fault formation. These observations are explained by generalized planar fault energies and grain-size effect on partial emissions.

Thermal fatigue on pistons induced by shaped high power laser. Part I: Experimental study of transient temperature field and temperature oscillation

Thermal fatigue behavior is one of the foremost considerations in the design and operation of diesel engines. It is found that thermal fatigue is closely related to the temperature field and temperature fluctuation in the structure. In this paper, spatially shaped high power laser was introduced to simulate thermal loadings on the piston. The incident Gaussian beam was transformed into concentric multi-circular beam of specific intensity distribution with the help of diffractive optical element (DOE), and the transient temperature fields in the piston similar to those under working conditions could be achieved by setting up appropriate loading cycles. Simulation tests for typical thermal loading conditions, i.e., thermal high cycle fatigue (HCF) and thermal shock (or thermal low cycle fatigue, LCF) were carried out. Several important parameters that affect the transient temperature fields and/or temperature oscillations, including controlling mode, intensity distribution of shaped laser, laser power, temporal profile of laser pulse, heating time and cooling time in one thermal cycle, etc., were investigated and discussed. The results show that as a novel method, the shaped high power laser can simulate thermal loadings on pistons efficiently, and it is helpful in the study of thermal fatigue behavior in pistons. (C) 2007 Elsevier Ltd. All rights reserved.

Thermal fatigue on pistons induced by shaped high power laser. Part II: Design of spatial intensity distribution via numerical simulation

In the laser induced thermal fatigue simulation test on pistons, the high power laser was transformed from the incident Gaussian beam into a concentric multi-circular pattern with specific intensity ratio. The spatial intensity distribution of the shaped beam, which determines the temperature field in the piston, must be designed before a diffractive optical element (DOE) can be manufactured. In this paper, a reverse method based on finite element model (FEM) was proposed to design the intensity distribution in order to simulate the thermal loadings on pistons. Temperature fields were obtained by solving a transient three-dimensional heat conduction equation with convective boundary conditions at the surfaces of the piston workpiece. The numerical model then was validated by approaching the computational results to the experimental data. During the process, some important parameters including laser absorptivity, convective heat transfer coefficient, thermal conductivity and Biot number were also validated. Then, optimization procedure was processed to find favorable spatial intensity distribution for the shaped beam, with the aid of the validated FEM. The analysis shows that the reverse method incorporated with numerical simulation can reduce design cycle and design expense efficiently. This method can serve as a kind of virtual experimental vehicle as well, which makes the thermal fatigue simulation test more controllable and predictable. (C) 2007 Elsevier Ltd. All rights reserved.

Dynamic mode I perturbation solution for a moving crack unsteadily

Rice et al. (Jounal of Mechanics and Physics of Solids 42, 813-843) analyze the propagation of a planar crack with a nominally straight front in a model elastic solid with a single displacement component. Using the form of Willis er al. (Journal of the Mechanics and Physics of Solids 43, 319-341), of dynamic mode I weight functions for a moving crack, we address that problem solved by Rice ei al. in the 3D context of elastodynamic theory. Oscillatory crack tip motion results from constructive-destructive interference of stress intensity waves. Those waves, including system of the dilatational, shear and Rayleigh waves, interact on each other and with moving edge of crack, can lead to continuing fluctuations of the crack front and propagation velocity. (C) 1997 Elsevier Science Ltd.

Crack propagation in the power-law nonlinear viscoelastic material

An analysis on crack creep propagation problem of power-law nonlinear viscoelastic materials is presented. The creep incompressilility assumption is used. To simulate fracture behavior of craze region, it is assumed that in the fracture process zone near the crack tip, the cohesive stress sigma(f) acts upon the crack surfaces and resists crack opening. Through a perturbation method, i. e., by superposing the Mode-I applied force onto a referential uniform stress state, which has a trivial solution and gives no effect on the solution of the original problem, the nonlinear viscoelastic problem is reduced to linear problem. For weak nonlinear materials, for which the power-law index n similar or equal to 1, the expressions of stress and crack surface displacement are derived. Then, the fracture process zone local energy criterion is proposed and based on which the formulas of cracking incubation time t

Asymptotic analysis of plane-strain mode I steady-state crack growth in transformation toughening ceramics (II)

Based on a constitutive law which includes the shear components of transformation plasticity, the asymptotic solutions to near-tip fields of plane-strain mode I steadity propagating cracks in transformed ceramics are obtained for the case of linear isotropic hardening. The stress singularity, the distributions of stresses and velocities at the crack tip are determined for various material parameters. The factors influencing the near-tip fields are discussed in detail.

Stress intensity factor histories of a steadily propagating crack under 3-D combined mode traction

Elastodynamic stress intensity factor histories of an unbounded solid containing a semi-infinite plane crack that propagates at a constant velocity under 3-D time-independent combined mode loading are considered. The fundamental solution, which is the response of point loading, is obtained. Then, stress intensity factor histories of a general loading system are written out in terms of superposition integrals. The methods used here are the Laplace transform methods in conjunction with the Wiener-Hopf technique.

A contact crack problem in an infinite plate

The problem of an infinite plate with crack of length 2a loaded by the remote tensile stress P and a pair of concentrated forces Q is discussed. The value of the force Q for the initial contact of crack face is investigated and the contact length elevated, while the Q force increases. The problem is solved assuming that the stress intensity factor vanishes at the end point of the contact portion. By the Fredholm integral equation for the multiple cracks, the reduction of stress intensity factor due to Q is found. (C) 1999 Elsevier Science Ltd. All rights reserved.

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.