A frequency-domain BEM for 3D non-synchronous crack interaction analysis in elastic solids

A frequency-domain boundary element method (BEM) is presented for non-synchronous crack interaction analysis in three-dimensional (3D), infinite, isotropic and linear elastic solids with multiple coplanar cracks. The cracks are subjected to non-synchronous time-harmonic crack-surface loading with contrast frequencies. Hypersingular frequency-domain traction boundary integral equations (BIEs) are applied to solve the boundary value problem. A collocation method is adopted for solving the BIEs numerically. The local square-root behavior of the crack-opening-displacements at the crack-front is taken into account in the present method. For two coplanar penny-shaped cracks of equal radius subjected to non-synchronous time-harmonic crack-surface loading, numerical results for the dynamic stress intensity factors are presented and discussed.     

Strength distributions and size effects for 2D and 3D composites with Weibull fibers in an elastic matrix

Monte Carlo simulation and theoretical modeling are used to study the statistical failure modes in unidirectional composites consisting of elastic fibers in an elastic matrix. Both linear and hexagonal fiber arrays are considered, forming 2D and 3D composites, respectively. Failure is idealized using the chain-of-bundles model in terms of -bundles of length , which is the length-scale of fiber load transfer. Within each -bundle, fiber load redistribution is determined by local load-sharing models that approximate the in-plane fiber load redistribution from planar break clusters, as predicted from 2D and 3D shear-lag models. As a result the -bundle failure models are 1D and 2D, respectively. Fiber elements have random strengths following either a Weibull or a power-law distribution with shape and scale parameters and , respectively. Under Weibull fiber strength, failure simulations for 2D -bundles, reveal two regimes: When fiber strength variability is low (roughly >2) the dominant failure mode is by growing clusters of fiber breaks, one of which becomes catastrophic. When this variability is high (roughly 0<<2) cluster formation is suppressed by a dispersed failure mode due to the blocking effects of a few strong fibers. For 1D -bundles or for 2D -bundles under power-law fiber strength, the transitional value of drops to 1 or lower, and overall, it may slowly decrease with increasing bundle size. For the two regimes, closed-form approximations to the distribution of -bundle strength are developed under the local load-sharing model and an equal load-sharing model of Daniels, respectively. The results compare favorably with simulations on -bundles with up to 1500 fibers.     

Stability analysis of an explicit finite element scheme for plane wave motions in elastic solids

 Expressions for critical timesteps are provided for an explicit finite element method for plane elastodynamic problems in isotropic, linear elastic solids. Both 4-node and 8-node quadrilateral elements are considered. The method involves solving for the eigenvalues directly from the eigenvalue problem at the element level. The characteristic polynomial is of order 8 for 4-node elements and 16 for 8-node elements. Due to the complexity of these equations, direct solution of these polynomials had not been attempted previously. The commonly used critical time-step estimates in the literature were obtained by reducing the characteristic equation for 4-node elements to a second-order equation involving only the normal strain modes of deformation. Furthermore, the results appear to be valid only for lumped-mass 4-node elements. In this paper, the characteristic equations are solved directly for the eigenvalues using <ty>Mathematica<ty> and critical time-step estimates are provided for both lumped and consistent mass matrix formulations. For lumped-mass method, both full and reduced integration are considered. In each case, the natural modes of deformation are obtained and it is shown that when Poisson's ratio is below a certain transition value, either shear-mode or hourglass mode of deformation dominates depending on the formulation. And when Poisson's ratio is above the transition value, in all the cases, the uniform normal strain mode dominates. Consequently, depending on Poisson's ratio the critical time-step also assumes two different expressions. The approach used in this work also has a definite pedagogical merit as the same approach is used in obtaining time-step estimates for simpler problems such as rod and beam elements. Received: 8 January 2002 / Accepted: 12 July 2002 The support of NSF under grant number DMI-9820880 is gratefully acknowledged.     

Cohesive zone representations of failure between elastic or rigid solids and ductile solids

A cohesive zone model that describes tangential separation as well as normal separation along an interface is reviewed. The model is based on nonlinear traction-separation relations between the normal and tangential components of the interface tractions and relative displacements. To illustrate the application of the cohesive zone model in studies of material failure or crack growth, analyses of matrix-fibre debonding in metal matrix composites are presented, taking into account effects of residual stresses or of nonlocal plasticity for the matrix. Also studies of interface crack growth under mixed mode conditions are discussed.     

Spin diffusion in 2D and 3D quantum solids

A moment method is used to compute the anisotropic spin diffusion constant in two-dimensional (2D) adsorbed and bulk (3D) quantum solids in which the spin motion is induced by an exchange Hamiltonian. Computations are carried out in 2D for the square and triangular lattices and in 3D for the hcp lattice. It is assumed that there are pair and three-particle exchange processes only. Since, in hcp³He, exchange processes out of the basal planeJ may occur at a different rate from processes in the planeJ, comparison with experimental results on single crystals should allow the determination ofJ andJ. Our results are given as functions of the ratioy=J/J and of the angle betweenc axis and field gradient. The 2D triangular lattice is shown to correspond to the special casey=0 of the hcp lattice. Our square-lattice result compares well with that of Morita (who used a different technique), supporting the validity of our method.     

Computation of stresses in non-homogeneous elastic solids by local integral equation method: a comparative study

The paper deals with the numerical implementation of local integral equations for solution of boundary value problems and interior computations of displacements and their gradients in functionally graded elastic solids. Two kinds of meshless approximations and one element based approximation are employed in various formulations. The numerical stability, accuracy, convergence of accuracy and cost efficiency are investigated in numerous test examples with exact benchmark solutions.     

On 3D Cracks Intersecting a Free Surface in Laminated Composites

In two- and three-dimensional linear elasticity, the singularities together with matched asymptotic expansions allow to extend the brittle fracture mechanics. Although there exist some difference between 2D and 3D approaches, the usual crack tip singularity exponent remains in both cases the hinge value between strong and weak singularities. In 3D, 0 was expected to be also a hinge, but it seems difficult to exhibit solutions with a negative exponent. One aim of this paper is to investigate numerically such specific cases and to derive some asymptotics of the classical stress intensity factors. The second part is dedicated to prove that, in case of a small linear ligament, negative exponents cannot exist. This revised version was published online in July 2006 with corrections to the Cover Date.     

Microcrack nucleation and growth in elastic lamellar solids

The nucleation and growth of microcracks in elastic lamellar microstructures is studied numerically. The analyses are carried out within a framework where the continuum is characterized by two constitutive relations: one relating the stress and strain in the bulk material and the other relating the traction and separation across a specified set of cohesive surfaces. In such a framework, fracture initiation and crack growth, including micro-crack nucleation ahead of the main crack, arise naturally as a consequence of the imposed loading, without any additional assumptions concerning criteria for crack growth, crack path selection or micro-crack nucleation. Full transient analyses are carried out and plane strain conditions are assumed. The specific problem analyzed is a compact tension specimen with two regions of differing lamellar orientation separated by a fracture resistant layer of finite width d, which is small compared to the physical dimensions of the specimen. An initial crack, normal to the applied loading, is assumed to exist in the first region whose lamellar orientation is fixed. The lamellar orientation of the second region, , is varied, as is the thickness of the fracture resistant layer. It is found that microcrack nucleation in the second region is highly sensitive to the lamellar orientation in that region for small values of d. However, microcrack nucleation becomes rather insensitive to with increasing d. It is also shown that a linear elastic fracture mechanics model with one adjustable parameter gives good agreement with the numerical results for fracture initiation.     

A numerical analysis for cracks emanating from a surface semi-spherical cavity in an infinite elastic body by FRANC3D

This note is specifically concerned with cracks emanating from a surface semi-spherical cavity in an infinite body (see Fig. 1) by using the boundary element software FRANC3D developed by a fracture mechanics investigation group of Cornell University. The numerical results can reveal the effect of the geometry of the surface cavity on the stress intensity factors.     

Methods to compute 3D residual stress distributions in hyperelastic tubes with application to arterial walls

In this paper the problem of modeling three-dimensional residual stress distributions in hyperelastic tubes is addressed. First, the problem of a radially opened straight and bent tube, where the opening angle depends on the axial position, is explored with the semi-inverse method. As a result a rather complicated system of nonlinear partial differential equations is achieved which is difficult to solve. Second, a different approximate method considers the tube as a composition of two, three, four or more rings in the axial direction. Also here the opening angle of the tube depends on the axial position. Some numerical solutions for the stress components in the radial, circumferential and axial directions are analyzed in more detail. Third, the tube wall is divided into a number of radial layers, with different mechanical properties, and an approximate method to treat that problems is presented. It is emphasized that the proposed approach can also be used to compute 3D residual stress distributions in arterial walls. A final conclusion points to possible future research directions.     

Crack-tip fields for matrix cracks between dissimilar elastic and creeping materials

The stress fields near the tip of a matrix crack terminating at and perpendicular to a planar interface under symmetric in-plane loading in plane strain are investigated. The bimaterial interface is formed by a linearly elastic material and an elastic power-law creeping material in which the crack is located. Using generalized expansions at the crack tip in each region and matching the stresses and displacements across the interface in an asymptotic sense, a series asymptotic solution is constructed for the stresses and strain rates near the crack tip. It is found that the stress singularities, to the leading order, are the same in each material; the stress exponent is real. The oscillatory higher-order terms exist in both regions and stress higher-order term with the order of O(r°) appears in the elastic material. The stress exponents and the angular distributions for singular terms and higher order terms are obtained for different creep exponents and material properties in each region. A full agreement between asymptotic solutions and the full-field finite element results for a set of test examples with different times has been obtained.     

Boundary-only element solutions of 2D and 3D nonlinear and nonhomogeneous elastic problems

This paper presents a robust boundary element method (BEM) that can be used to solve elastic problems with nonlinearly varying material parameters, such as the functionally graded material (FGM) and damage mechanics problems. The main feature of this method is that no internal cells are required to evaluate domain integrals appearing in the conventional integral equations derived for these problems, and very few internal points are needed to improve the computational accuracy. In addition, one of the basic field quantities used in the boundary integral equations is normalized by the material parameter. As a result, no gradients of the field quantities are involved in the integral equations. Another advantage of using the normalized quantities is that no material parameters are included in the boundary integrals, so that a unified equation form can be established for multi-region problems which have different material parameters. This is very efficient for solving composite structural problems.     

3D modelling of strong discontinuities in elastoplastic solids: fixed and rotating localization formulations

This paper is concerned with the incorporation of displacement discontinuities into a continuum theory of elastoplasticity for the modelling of localization processes such as cracking in brittle materials. Based on the strong discontinuity approach (SDA) (Computational Mechanics 1993; 12: 277–296) mesh objective 2D and 3D finite element formulations are developed using linear and quadratic 2D elements as well as 8‐noded 3D elements. In the formulation of the finite‐element model proposed in the paper, the analogy with standard formulations is emphasized. The parameter defining the amplitude of the displacement jump within the finite element is condensed out at the material level without employing the standard static condensation technique. This approach results in linearized constitutive equations formally identical to continuum models. Therefore, the standard return mapping algorithm is used to solve the non‐linear equations. In analogy to concepts used in continuum smeared crack models, a rotating formulation of the SDA is proposed in addition to the standard concept of fixed discontinuities. It is shown that the rotating localization approach reduces locking effects observed in analyses based on fixed localization directions. The applicability of the proposed SDA finite‐element model as well as its numerical performance is investigated by means of a three‐dimensional ultimate load analysis of a steel anchor embedded in a concrete block subjected to a shear force. Copyright © 2003 John Wiley & Sons, Ltd.     

Critical field and impurity scattering: Application to layered and 3D superconducting systems

An efficient and rigorous method of evaluating the critical fieldH _c2 is developed. The method allows one to consider in a regular way the angular dependence of the scattering cross-section. The method is applied to layered and 3D superconductors.     

Refined and generalized hybrid type quasi-3D shear deformation theory for the bending analysis of functionally graded shells

The closed-form solution of a generalized hybrid type quasi-3D higher order shear deformation theory (HSDT) for the bending analysis of functionally graded shells is presented. From the generalized quasi-3D HSDT (which involves the shear strain functions “f(ζ)” and “g(ζ)” and therefore their parameters to be selected “m” and “n”, respectively), infinite six unknowns' hybrid shear deformation theories with thickness stretching effect included, can be derived and solved in a closed-from. The generalized governing equations are also “m” and “n” parameter dependent. Navier-type closed-form solution is obtained for functionally graded shells subjected to transverse load for simply supported boundary conditions. Numerical results of new optimized hybrid type quasi-3D HSDTs are compared with the first order shear deformation theory (FSDT), and other quasi-3D HSDTs. The key conclusions that emerge from the present numerical results suggest that: (a) all non-polynomial HSDTs should be optimized in order to improve the accuracy of those theories; (b) the optimization procedure in all the cases is, in general, beneficial in terms of accuracy of the non-polynomial hybrid type quasi-3D HSDT; (c) it is possible to gain accuracy by keeping the unknowns constant; (d) there is not unique quasi-3D HSDT which performs well in any particular example problems, i.e. there exists a problem dependency matter.     

Structural, elastic and electronic properties of intermetallics in the Pt–Sn system: A density functional investigation

The structural, elastic and electronic properties of intermetallics in the Pt–Sn binary system are investigated using first-principles calculations based on density functional theory (DFT). The polycrystalline elastic properties are deduced from the calculated single-crystal elastic constants. The elastic anisotropy of these intermetallics is analyzed based on the directional dependence of the Young’s modulus and its origin explained based on the electronic nature of the crystals. All the Pt–Sn intermetallics investigated are found to be mechanically stable, ductile and metallic, and some of them show high elastic anisotropy.     

A fast and accurate analysis of the interacting cracks in linear elastic solids

This paper presents a fast and accurate solution for crack interaction problems in infinite- and half- plane solids. The new solution is based on the method of complex potentials developed by Muskhelishvili for the analysis of plane linear elasticity, and it is formulated through three steps. First, the problem is decomposed into a set of basic problems, and for each sub-problem, there is only one crack in the solid. Next, after a crack-dependent conformal mapping, the modified complex potentials associated with the sub-problems are expanded into Laurent’s series with unknown coefficients, which in turn provides a mechanism to exactly implement in the form of Fourier series the boundary condition in each sub-problem. Finally, taking into account the crack interaction via a perturbation approach, an iterative algorithm based on fast Fourier transforms (FFT) is developed to solve the unknown Fourier coefficients, and the solution of the whole problem is readily obtained with the superposition of the complex potentials in each sub-problem. The performance of the proposed method is fully investigated by comparing with benchmark results in the literatures, and superb accuracy and efficiency is observed in all situations including patterns where cracks are closely spaced. Also, the new method is able to cope with interactions among a large number of cracks, and this capability is demonstrated by a calculation of effective moduli of an elastic solid with thousands of randomly-spaced cracks.     

Stabilized finite elements for time-harmonic waves in incompressible and nearly incompressible elastic solids

The propagation of waves in elastic solids at or near the incompressible limit is of interest in many current and emerging applications. Standard low-order Galerkin finite element discretization struggles with both incompressibility and wave dispersion. Galerkin least squares stabilization is known to improve computational performance of each of these ingredients separately. A novel approach of combined pressure-curl stabilization is presented, facilitating the use of continuous, equal-order interpolation of displacements and pressure. The pressure stabilization parameter is determined by stability considerations, while the curl stabilization parameter is determined by dispersion considerations. The proposed pressure-curl–stabilized scheme provides stable and accurate results on a variety of numerical tests for incompressible and nearly incompressible elastic waves computed with linear elements.