A new method for evaluation of stress intensities for interface cracks

A new method is presented for calculating the values of K-I and K-II in the elasticity solution at the tip of an interface crack. The method is based on an evaluation of the J-integral by the virtual crack extension method. Expressions for calculating K-I and K-II by using the displacements and the stiffness derivative of the finite element solution and asymptotic crack tip displacements are derived. The method is shown to produce very accurate solutions even with coarse element mesh.

Stress-distribution at the rigid circular-arc inclusion end

The elastic plane problem of a rigid co-circular arc inclusion under arbitrary loads is dealt with. Applying Schwarz's reflection principle integrated with the analysis of the singularity of complex stress functions, the general solution of the problem is found and several closed-form solutions to some problems of practical importance are given. Finally, the stress distribution at the arc inclusion end is examined and a comparison is made with that of the rigid line inclusion end to show the effect of curvature.

High order asymptotic field for plane stress crack problem

This paper presents an exact analysis for high order asymptotic field of the plane stress crack problem. It has been shown that the second order asymptotic field is not an independent eigen field and should be matched with the elastic strain term of the first order asymptotic field. The second order stress field ahead of the crack tip is quite small compared with the first order stress field. The stress field ahead of crack tip is characterized by the HRR field. Hence the J integral can be used as a criterion for crack initiation.

Fatigue crack-growth from a circular notch under high-levels of biaxial stress

A series of experiments have been conducted on cruciform specimens to investigate fatigue crack growth from circular notches under high levels of biaxial stress. Two stress levels (Δσ₁= 380 and 560 MPa) and five stress biaxialities (λ=+1.0, +0.5, 0, −0.5 and −1.0; where λ=σ₂/σ₁ were adopted in the fatigue tests in type 316 stainless steel having a monotonic yield strength of 243 MPa. The results reveal that fatigue crack growth rates are markedly influenced by both the stress amplitude and the stress biaxiality. A modified model has been developed to describe fatigue crack growth under high levels of biaxial stress.

Fracture behavior of axisymmetrical bars under high triaxial stress and large strain cyclic loading

Axisymmetric notched bars with notch roots of large and small radii were tested under large strain cyclic loading. The main attention is focused on the fracture behaviour of steels having cycles to failure within the range 1-100. Our study shows that a gradual transition from a static ductile nature to one of fatigue cleavage can be observed and characterized by the Coffin-Manson formula in a generalized form. Both the triaxial tensile stress within the central region of specimens and static damage caused by the first increasing load have effects on the final failure event. A generalized cyclic strain range parameter DELTAepsilon is proposed as a measure of the numerous factors affecting behaviour. Fractographs are presented to illustrate the behaviour reported in the paper.

Plastic stress-strain fields and blunting effects during large deformations

Plastic stress-strain fields of two types of steel specimens loaded to large deformations are studied. Computational results demonstrate that, owing to the fact that the hardening exponent of the material varies as strain enlarges and the blunting of the crack tip, the well known HRR stress field in the plane strain model can only hold for the stage of a small plastic strain. Plastic dilatancy is shown to have substantial effects on strain distributions and blunting. To justify the constitutive equations used for analysis and to check the precision of computations, the load-deflection of a three-point bend beam and the load-elongation of an axisymmetric bar notched by a V-shaped cut were tested and recorded. The computed curves are in good accordance with experimental data.

Stress-distribution around a rigid line in dissimilar media

The elastic plane problem of a rigid line inclusion between two dissimilar media was considered. By solving the Riemann-Hilbert problem, the closed-form solution was obtained and the stress distribution around the rigid line was investigated. It was found that the modulus of the singular behavior of the stress remains proportional to the inverse square root of the distance from the rigid line end, but the stresses possess a pronounced oscillatory character as in the case of an interfacial crack tip.

An asymptotic analysis for one dimensional stress wave propagation in damaged viscoelastic media

A regular perturbation technique is suggested to deal with the problem of one dimensional stress wave propagation in viscoelastic media with damage. Based upon the first order asymptotic solution obtained, the characteristics of wave attenuation are studied. In fact, there exist three different time-dependent phenomena featuring the dynamic response of the materials, the first expressing the characteristics of wave propagation, the second indicating the innate effect of visco-elastic matrix and the third coming from the time dependent damage. The comparision of first order asymptotic solution with the numerical results calculated by a finite difference procedure shows that the perturbation expansion technique may offer a useful approach to the problem concerned.

A discussion about a class of stress intensity factors and its verification

In this paper, new formulae of a class of stress intensity factors for an infinite plane with two collinear semi-infinite cracks are presented. The formulae differ from those gathered in several handbooks used all over the world. Some experiments and finite element calculations have been developed to verify the new formulae and the results have shown its reliability. Finally, the new formulae and the old are listed to show the differences between them.

Complex stress functions and plastic zone sizes for notch cracks subjected to various loading conditions

In this paper, the conformal mapping method is used to solve the plane problem of an infinite plate containing a central lip-shaped notch subjected to biaxial loading at a remote boundary or a surface uniform pressure on the notch. The stress intensity factors KI and KII are obtained by the derived complex stress functions. The simple analytical expressions can be applied to the situation of cracks originating from a circular or an elliptical notch. The plastic zone sizes for such notch cracks are subsequently evaluated in light of the Dugdale strip yield concept. The results are consistent with available numerical data.

Inhomogeneous isoparametric elements for stress analysis of composites

The concept of inhomogeneous element is proposed and the formulations of the inhomogeneous isoparametric elements for stress analysis of four kinds of problems are derived. As an example of applications of the inhomogeneous elements, the stress distribution in a cone-like composite syntheticrope termination is calculated.

3D Finite Volume Scheme for Czochralski Single Crystal Growth

Czochralski (Cz) technique, which is used for growing single crystals, has dominated the production of single crystals for electronic applications. The Cz growth process involves multiple phases, moving interface and three-dimensional behavior. Much has been done to study these phenomena by means of numerical methods as well as experimental observations. A three-dimensional curvilinear finite volume based algorithm has been developed to model the Cz process. A body-fitted transformation based approach is adopted in conjunction with a multizone adaptive grid generation (MAGG) technique to accurately handle the three-dimensional problems of phase-change in irregular geometries with free and moving surfaces. The multizone adaptive model is used to perform a three-dimensional simulation of the Cz growth of silicon single crystals.Since the phase change interface are irregular in shape and they move in response to the solution, accurate treatment of these interfaces is important from numerical accuracy point of view. The multizone adaptive grid generation (MAGG) is the appropriate scheme for this purpose. Another challenge encountered is the moving and periodic boundary conditions, which is essential to the numerical solution of the governing equations. Special treatments are implemented to impose the periodic boundary condition in a particular direction and to determine the internal boundary position and shape varying with the combination of ambient physicochemical transport process and interfacial dynamics. As indicated above that the applications and processes characterized by multi-phase, moving interfaces and irregular shape render the associated physical phenomena three-dimensional and unsteady. Therefore a generalized 3D model rather than a 2D simulation, in which the governing equations are solved in a general non-orthogonal coordinate system, is constructed to describe and capture the features of the growth process. All this has been implemented and validated by using it to model the low pressure Cz growth of silicon. Accuracy of this scheme is demonstrated by agreement of simulation data with available experimental data. Using the quasi-steady state approximation, it is shown that the flow and temperature fields in the melt under certain operating conditions become asymmetric and unsteady even in the absence of extrinsic sources of asymmetry. Asymmetry in the flow and temperature fields, caused by high shear initiated phenomena, affects the interface shape in the azimuthal direction thus results in the thermal stress distribution in the vicinity, which has serious implications from crystal quality point of view.