Computation of stress intensity factors by the sub-region mixed finite element method of lines

Based on the sub-region generalized variational principle, a sub-region mixed version of the newly-developed semi-analytical 'finite element method of lines' (FEMOL) is proposed in this paper for accurate and efficient computation of stress intensity factors (SIFs) of two-dimensional notches/cracks. The circular regions surrounding notch/crack tips are taken as the complementary energy region in which a number of leading terms of singular solutions for stresses are used, with the sought SIFs being among the unknown coefficients. The rest of the arbitrary domain is taken as the potential energy region in which FEMOL is applied to obtain approximate displacements. A mixed system of ordinary differential equations (ODEs) and algebraic equations is derived via the sub-region generalized variational principle. A singularity removal technique that eliminates the stress parameters from the mixed equation system eventually yields a standard FEMOL ODE system, the solution of which is no longer singular and is simply and efficiently obtained using a standard general-purpose ODE solver. A number of numerical examples, including bi-material notches/cracks in anti-plane and plane elasticity, are given to show the generally excellent performance of the proposed method.

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.

On the unstable stacking criterion for ideal and cracked copper-crystals

The unstable stacking criteria for an ideal copper crystal under homogeneous shearing and for a cracked copper crystal under pure mode II loading are analysed. For the ideal crystal under homogeneous shearing, the unstable stacking energy gamma(us) defined by Rice in 1992 results from shear with no relaxation in the direction normal to the slip plane. For the relaxed shear configuration, the critical condition for unstable stacking does not correspond to the relative displacement Delta = b(p)/2, where b(p) is the Burgers vector magnitude of the Shockley partial dislocation, but to the maximum shear stress. Based on this result, the unstable stacking energy Gamma(us) is defined for the relaxed lattice. For the cracked crystal under pure mode II loading, the dislocation configuration corresponding to Delta = b(p)/2 is a stable state and no instability occurs during the process of dislocation nucleation. The instability takes place at approximately Delta = 3b(p)/4. An unstable stacking energy Pi(us) is defined which corresponds to the unstable stacking state at which the dislocation emission takes place. A molecular dynamics method is applied to study this in an atomistic model and the results verify the analysis above.

The inhomogeneity of plastic-deformation in ductile single-crystals

A discrete slip model which characterizes the inhomogeneity of material properties in ductile single crystals is proposed in this paper. Based on this model rate-dependent finite element investigations are carried out which consider the finite deformation, finite rotation, latent hardening effect and elastic anisotropy. The calculation clearly exhibits the process from microscopic inhomogeneous and localized deformation to necking and the formation of LSBS and reveals several important features of shear localization. For example, the inhomogeneous deformation is influenced by the imperfections and initial non-uniformities of material properties. The inhomogeneous deformation may either induce necking which results in the lattice rotation and leads to geometrical softening, which in turn promotes the formation of CSBS, or induces heavily localized deformation. The microscopic localized deformation eventually develops into the LSBS and results in a failure. These results are in close agreement with experiment. Our calculations also find that the slip lines on the specimen's surface at necking become curved and also find that if the necking occurs before the formation of LSBS, this band must be misoriented from the operative slip systems. In this case, the formation of LSBS must involve non-crystallographic effects. These can also be indirectly confirmed by experiment. All these suggest that our present discrete slip model offers a correct description of the inhomogeneous deformation characterization in ductile crystals.

Growth of zero-mean-shear turbulent-mixed layer with suspended particles in a 2-layer fluid system

The growth behaviour of zero-mean-shear turbulent-mixed layer containing suspended solid particles has been studied experimentally and analysed theoretically in a two-layer fluid system. The potential model for estimating the turbulent entrainment rate of the mixed layer has also been suggested, including the results of the turbulent entrainment for pure two-layer fluid. The experimental results show that the entrainment behaviour of a mixed layer with the suspended particles is well described by the model. The relationship between the entrainment distance and the time, and the variation of the dimensionless entrainment rate E with the local Richardson number Ri1 for the suspended particles differ from that for the pure two-layer fluid by the factors-eta-1/5 and eta-1, respectively, where eta = 1 + sigma-0-DELTA-rho/DELTA-rho-0.

Crack extension and kinking in laminates and bicrystals

Singular fields at the tip of an interface crack in anisotropic solids are reviewed with emphasis on establishing a framework to quantify fracture resistance under mixed mode conditions. The concepts of mode mixity and surface toughness are unified by using generalized interface traction components. The similarity between the anisotropic theory and existing isotropic theory is shown. Explicit formulae are given for misoriented orthotropic bimaterials with potential applications envisioned including composite laminates and semiconductor crystals. Competition between crack extension along the interface and kinking into the substrate is investigated using a boundary layer formulation. Several case studies reveal the role of anisotropy. An explicit complex variable representation for orthotropic materials and a solution to a dislocation interacting with a crack are presented in two self-contained Appendices.

Prediction of mixed-mode cracking direction in random, short-fibre composite-materials

The prediction of cracking direction in composite materials is of significance to the design of composite structures. This paper presents several methods for predicting the cracking direction in the double grooved tension-shear specimen which gives mixed-mode cracking. Five different criteria are used in this analysis: two of them have been used by other investigators and the others are proposed by the present authors. The strain energy density criterion proposed by G.C. Sih is modified to take account of the influence of the anisotropy of the strength on the direction of crack. The two failure criteria of Tsai-Hill and Norris are extended to predict the crack orientation. The stress distributions in the near-notch zone are calculated by using the 8-node quadrilateral isoparametric finite element method. The predictions of all the criteria except one are in good agreement with the experimental measurement. In addition, on the basis of the FEM results, the size of the zone in which the singular term is dominant is estimated.

Additional response: Additional discussion of “Mixed mode ductile fracture using the strain energy density criterion”

boundary-layer flows, the skin friction and wall heat-transfer are higher and the