The mixed-type integral equations for transient elastodynamic plane crack problems

In the present paper, by use of the boundary integral equation method and the techniques of Green fundamental solution and singularity analysis, the dynamic infinite plane crack problem is investigated. For the first time, the problem is reduced to solving a system of mixed-typed integral equations in Laplace transform domain. The equations consist of ordinary boundary integral equations along the outer boundary and Cauchy singular integral equations along the crack line. The equations obtained are strictly proved to be equivalent with the dual integral equations obtained by Sih in the special case of dynamic Griffith crack problem. The mixed-type integral equations can be solved by combining the numerical method of singular integral equation with the ordinary boundary element method. Further use the numerical method for Laplace transform, several typical examples are calculated and their dynamic stress intensity factors are obtained. The results show that the method proposed is successful and can be used to solve more complicated problems.

Lattice instability at a fast moving crack tip

A molecular dynamics method is used to analyze the dynamic propagation of an atomistic crack tip. The simulation shows that the crack propagates at a relatively constant global velocity which is well below the Rayleigh wave velocity. However the local propagation velocity oscillates violently, and it is limited by the longitudinal wave velocity. The crack velocity oscillation is caused by a repeated process of crack tip blunting and sharpening. When the crack tip opening displacement exceeds a certain critical value, a lattice instability takes place and results in dislocation emissions from the crack tip. Based on this concept, a criterion for dislocation emission from a moving crack tip is proposed. The simulation also identifies the emitted dislocation as a source for microcrack nucleation. A simple method is used to examine this nucleation process. (C) 1996 American Institute of Physics.

Molecular dynamics simulation of interaction of a dislocation array from a crack tip with grain boundaries

The interaction of a dislocation array emitted from a crack tip under mode II loading with asymmetric tilt grain boundaries (GBs) is analysed by the molecular dynamics method. The GBs can generally be described by planar and linear matching zones and unmatching zones. All GBs are observed to emit dislocations. The GBs migrated easily due to their planar and linear matching structure and asymmetrical type. The diffusion induced by stress concentration is found to promote the GB migration. The transmissions of dislocations are either along the matched plane or along another plane depending on tilt angle theta. Alternate processes of stress concentration and stress relaxation take place ahead of the pileup. The stress concentration can be released either by transmission of dislocations, by atom diffusion along GBs, or by migration of GBs by formation of twinning bands. The simulated results also unequivocally demonstrate two processes, i.e. asymmetrical GBs evolving into symmetrical ones and unmatching zones evolving into matching ones during the loading process.

Micromechanics analyses of interface crack

A new mechanics model based on Peierls concept is presented in this paper, which can clearly characterize the intrinsic features near a tip of an interfacial crack. The stress and displacement fields are calculated under general combined tensile and shear loadings. The near tip stress fields show some oscillatory behaviors but without any singularity and the crack faces open completely without any overlapping when remote tensile loading is comparable with remote shear loading. A fracture criterion for predicting interface toughness has been also proposed, which takes into account for the shielding effects of emitted dislocations. The theoretical toughness curve gives excellent prediction, as compared with the existing experiment data.

The annular crack in a nonhomogeneous matrix surrounding a fiber under torsional loading .2. The crack going through the bimaterial interface into the fiber

The axisymmetric problem of an elastic fiber perfectly bonded to a nonhomogeneous elastic matrix which contains an annular crack going through the interface into the fiber under axially symmetric shear stress is considered. The nature of the stress singularity is studied. It is shown that at the irregular point on the interface, whether the shear modulus is continuous or discontinuous the stresses are bounded. The problem is formulated in terms of a singular integral equation and can be solved by a regular method. The stress intensity factors and crack surface displacement are given.

Experimental investigation on low-temperature fatigue crack propagation in offshore structural steel A131 under random sea ice loading

An apparatus of low-temperature controlling for fatigue experiments and its crack measuring system were developed and used for offshore structural steel A131 under conditions of both low temperature and random sea ice. The experimental procedures and data processing were described, and a universal random data processing software for FCP under spectrum loading was written. Many specific features of random ice-induced FCP which differed with constant amplitude FCP behaviours were proposed and temperature effect on ice-induced FCP was pointed out with an easily neglected aspect in designing for platforms in sea ice emphasized. In the end, differences of FCP behaviours between sea ice and ocean wave were presented.

Transient effect on ideally plastic stationary crack-tip fields

This paper analyses the transient effect on ideally plastic stationary crack-tip fields under mode I plane strain conditions, when the inertial forces are not negligible. It is shown that the governing equation for such a problem can be expressed in formal simplicity when referred to a system of moving curvilinear coordinates, which is a generalization of the system defined by the slip-line field in quasi-static plasticity. A perturbation method of solving the equations is described and illustrated by application to problems of ideally plastic stationary crack-tip fields when the inertia forces are not negligible.

Transient effect on ideally plastic stationary crack-tip fields

This paper analyses the transient effect on ideally plastic stationary crack tip fields under mode I plane strain conditions, when the inertial forces are not negligible. It is shown that the governing equation for such a problem can be expressed in formal simplicity when referred to a system of moving curvilinear coordinates, which is a generalization of the system defined by the slip-line field in quasi-static plasticity. A perturbation method of solving the equations is described and illustrated by application to problems of ideally plastic stationary crack tip fields when the inertial forces are not negligible.

Molecular dynamics simulation of crack-tip processes in copper

The crack tip processes in copper under mode II loading have been simulated by a molecular dynamics method. The nucleation, emission, dislocation free zone (DFZ) and pile-up of the dislocations are analyzed by using a suitable atom lattice configuration and Finnis & Sinclair potential. The simulated results show that the dislocation emitted always exhibits a dissociated fashion. The stress intensity factor for dislocation nucleation, DFZ and dissociated width of partial dislocations are strongly dependent on the loading rate. The stress distributions are in agreement with the elasticity solution before the dislocation emission, but are not in agreement after the emission. The dislocation can move at subsonic wave speed (less than the shear wave speed) or at transonic speed (greater than the shear wave speed but less than the longitudinal wave speed), but at the longitudinal wave speed the atom lattice breaks down.

Dislocation nucleation and emission from crack-tip

The problems of dislocation nucleation and emission from a crack tip are analysed based on Peierls model. The concept adopted here is essentially the same as that proposed by Rice. A slight modification is introduced here to identify the pure linear elastic response of material. A set of new governing equations is developed, which is different from that used by Beltz and Rice. The stress field and the dislocation density field can be expressed as the first and second Chebyshev polynomial series respectively. Then the opening and slip displacements can be expanded as the trigonometric series. The Newton-Raphson Method is used to solve a set of nonlinear algebraic equations. The new governing equations allow us to extend the analyses to the case of dislocation emission. The calculation results for pure shearing, pure tension and combined tension and shear loading are given in detail.

Elastodynamic stress-intensity factors for a semiinfinite crack under 3-d combined mode loading

The dynamic stress intensity factor histories for a half plane crack in an otherwise unbounded elastic body are analyzed. The crack is subjected to a traction distribution consisting of two pairs of suddenly-applied shear point loads, at a distance L away from the crack tip. The exact expression for the combined mode stress intensity factors as the function of time and position along the crack edge is obtained. The method of solution is based on the direct application of integral transforms together with the Wiener-Hopf technique and the Cagniard-de Hoop method, which were previously believed to be inappropriate. Some features of solutions are discussed and the results are displayed in several figures.

The effect of thermal-activation on dislocation processes at an atomistic crack-tip

The effects of thermal activation on the dislocation emission from an atomistic crack tip are discussed, Molecular dynamics simulations at different constant temperatures are carried out to investigate the thermal effects. The simulated results show that the processes of the partial dislocation generation and emission are temperature dependent. As the temperature increases, the incipient duration of the partial dislocation nucleation becomes longer, the critical stress intensity factor for partial dislocation emission is reduced and, at the same loading level, more dislocations are emitted. The dislocation velocity moving away from the crack tip and the separations of partial dislocations are apparently not temperature dependent. The simulated results also show that, as the temperature increases, the stress distribution along the crack increases slightly. Therefore stress softening at the crack tip induced by thermal activation does not exist in the present simulation. A simple model is proposed to evaluate the relation of the critical stress intensity factor versus temperature. The obtained relation is in good agreement with our molecular dynamics results.

Unloading characteristics of anti-plane shear crack in softening materials

The local characteristics of the anti-plane shear stress and strain field are determined for a material where the stress increases linearly with strain up to a limit and then softens nonlinearly. Two unloading models are considered such that the unloading path always returns to the origin while the other assumes the unloading modulus to be that of the initial shear modulus. As the applied shear increases, an unloading zone is found to prevail between a zone in which the material softens and another zone in which the material is linear-elastic even though the crack does not propagate. The divisions of these zones are displayed graphically.

Elastodynamic stress intensity factor history for a semi-infinite crack under three-dimensional transient loading

The dynamic stress intensity factor history for a semi-infinite crack in an otherwise unbounded elastic body is analyzed. The crack is subjected to a pair of suddenly-applied point loadings on its faces at a distance L away from the crack tip. The exact expression for the mode I stress intensity factor as a function of time is obtained. The method of solution is based on the direct application of integral transforms, the Wiener-Hopf technique and the Cagniard-de Hoop method. Due to the existence of the characteristic length in loading this problem was long believed a knotty problem. Some features of the solutions are discussed and graphical result for numerical computation is presented.