A simple and less-costly meshless local Petrov-Galerkin (MLPG) method for the dynamic fracture problem

A simple and less-costly MLPG method using the Heaviside step function as the test function in each sub-domain avoids the need for both a domain integral, except inertial force and body force integral in the attendant symmetric weak form, and a singular integral for analysis of elasto-dynamic deformations near a crack tip. The Newmark family of the methods is applied into the time integration scheme. A numerical example, namely, a rectangular plate with a central crack with plate edges parallel to the crack axis loaded in tension is solved by this method. The results show that the stresses near the crack tip agree well with those obtained from another MLPG method using the weight function of the moving least square approximation as a test function of the weighted residual method. Time histories of dynamic stress intensity factors (DSIF) for mode-I are determined form the computed stress fields.     

Application of the weight function method to two-dimensional elastodynamic fracture mechanics

In this paper, the weight function method is used for two-dimensional mixed-mode crack analyses of clastostatic and elastodynamic problems. By the use of the Laplace transformation method and an indirect boundary element method, the dynamic stress intensity factors for a finite sheet containing a central or an edge crack are evaluated. A Green's function method is introduced which depends on the weight function for an impulsive applied load. The Green's function can be used to determine stress intensity factors for arbitrary time dependence of the boundary conditions. The stress intensity factors obtained by the weight function method are compared where possible, with existing solutions.     

An indirect boundary element method for three-dimensional dynamic fracture mechanics

Indirect boundary element methods (fictitious load and displacement discontinuity) have been developed for the analysis of three-dimensional elastostatic and elastodynamic fracture mechanics problems. A set of boundary integral equations for fictitious loads and displacement discontinuities have been derived. The stress intensity factors were obtained by the stress equivalent method for static loading. For dynamic loading the problem was studied in Laplace transform space where the numerical calculation procedure, for the stress intensity factor K_I(p), is the same: as that for the static problem. The Durbin inversion method for Laplace transforms was used to obtain the stress intensity factors in the time domain K_I(t). Results of this analysis are presented for a square bar, with either a rectangular or a circular crack, under static and dynamic loads.     

An effective method of stress intensity factor calculation for cracks emanating from a triangular or square hole under biaxial loads

This paper is concerned with stress intensity factors for cracks emanating from a triangular or square hole under biaxial loads by means of a new boundary element method. The boundary element method consists of the constant displacement discontinuity element presented by Crouch and Starfied and the crack‐tip displacement discontinuity elements proposed by the author. In the boundary element implementation, the left or the right crack‐tip displacement discontinuity element is placed locally at the corresponding left or right crack tip on top of the constant displacement discontinuity elements that cover the entire crack surface and the other boundaries. The method is called a Hybrid Displacement Discontinuity Method (HDDM). Numerical examples are included to show that the method is very efficient and accurate for calculating stress intensity factors for plane elastic crack problems. In addition, the present numerical results can reveal the effect of the biaxial loads on stress intensity factors.     

A variational technique for boundary element analysis of 3D fracture mechanics weight functions: static

The dual boundary element method coupled with the weight function technique is developed for the analysis of three-dimensional elastostatic fracture mechanics mixed-mode problems. The weight functions used to calculate the stress intensity factors are defined by the derivatives of traction and displacement for a reference problem. A knowledge of the weight functions allows the stress intensity factors for any loading on the boundary to be calculated by means of a simple boundary integration without singularities. Values of mixed-mode stress intensity factors are presented for an edge crack in a rectangular bar and a slant circular crack embedded in a cylindrical bar, for both uniform tensile and pure bending loads applied to the ends of the bars. © 1998 John Wiley & Sons, Ltd.     

Elastodynamic stress-intensity factors for a semi-infinite 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.     

Green's functions for the stress intensity factor evolution in finite cracks in orthotropic materials

The elastodynamic response of an infinite orthotropic material with finite crack under concentrated loads is examined. Solution for the stress intensity factor history around the crack tips is found. Laplace and Fourier transforms are employed to solve the equations of motion leading to a Fredholm integral equation on the Laplace transform domain. The dynamic stress intensity factor history can be computed by numerical Laplace transform inversion of the solution of the Fredholm equation. Numerical values of the dynamic stress intensity factor history for some example materials are obtained. This solution can be used as a Green's function to solve dynamic problems involving fini te cracks.     

Dynamic propagation of a finite crack in a micropolar elastic solid

Summary The dynamic propagation of a finite crack under mode-I loading in a micropolar elastic solid is investigated. By using an integral transform method, a pair of two-dimensional singular integral equations governing stress and couple stress is formulated in terms of displacement transverse to the crack, macro and micro rotations, and microinertia. These equations are solved numerically, and solutions for dynamic stress intensity and couple stress intensity factors are obtained by utilizing the values of the strengths of the square root singularities in macrorotation and the gradient of microrotation at the crack tips.     

A weight function method for mixed modes hole‐edge cracks

A weight function approach is proposed to calculate the stress intensity factor and crack opening displacement for cracks emanating from a circular hole in an infinite sheet subjected to mixed modes load. The weight function for a pure mode II hole‐edge crack is given in this paper. The stress intensity factors for a mixed modes hole‐edge crack are obtained by using the present mode II weight function and existing mode I Green (weight) function for a hole‐edge crack. Without complex derivation, the weight functions for a single hole‐edge crack and a centre crack in infinite sheets are used to study 2 unequal‐length hole‐edge cracks. The stress intensity factor and crack opening displacement obtained from the present weight function method are compared well with available results from literature and finite element analysis. Compared with the alternative methods, the present weight function approach is simple, accurate, efficient, and versatile in calculating the stress intensity factor and crack opening displacement.     

Dynamic stress intensity factor due to concentrated loads on a propagating semi-infinite crack in orthotropic materials

The elastodynamic response of an infinite orthotropic material with a semi-infinite crack propagating at constant speed under the action of concentrated loads on the crack faces is examined. Solution for the stress intensity factor history around the crack tip is found for the loading modes I and II. Laplace and Fourier transforms along with the Wiener-Hopf technique are employed to solve the equations of motion. The asymptotic expression for the stress near the crack tip is analyzed which lead to a closed-form solution of the dynamic stress intensity factor. It is found that the stress intensity factor for the propagating crack is proportional to the stress intensity factor for a stationary crack by a factor similar to the universal function k(v) from the isotropic case. Results are presented for orthotropic materials as well as for the isotropic case.     

Three-dimensional crack problem analysis using boundary element method with finite-part integrals

A set of hypersingular integral equations of a three-dimensional finite elastic solid with an embedded planar crack subjected to arbitrary loads is derived. Then a new numerical method for these equations is proposed by using the boundary element method combined with the finite-part integral method. According to the analytical theory of the hypersingular integral equations of planar crack problems, the square root models of the displacement discontinuities in elements near the crack front are applied, and thus the stress intensity factors can be directly calculated from these. Finally, the stress intensity factor solutions to several typical planar crack problems in a finite body are evaluated. This revised version was published online in August 2006 with corrections to the Cover Date.     

Weight functions for bimaterial interface cracks

An efficient approach using the analytically decoupled near-tip displacement solution for bimaterial interface cracks presented in this paper involves: (1) the calculation of the decoupled strain energy release rates G _I and G _II associated respectively with the decoupled stress intensity factors K _I and K _II and (2) the extension of Rice's displacement derivative representation of Bueckner's weight function vectors beyond the homogeneous media. It is shown that the stress intensity factors for a bimaterial interface crack predicted by the present approach agree very well with those solutions available in the literature. The computational efficiency is enhanced through the use of singular elements in the crack-tip neighborhood.As reported in the homogeneous case, the calculated weight function for a bimaterial interface crack is load-independent but depends strongly on geometry and constraint conditions. Due to the coupling nature of the stress intensity factors of a bimaterial interface crack, the invariant characteristics of the dimensionless weight function vectors are different from those of a crack in homogeneous material. In addition, the elastic constants of two constituents can significantly alter the weight function behavior for a cracked bimaterial medium.Due to the load-independent characteristic of the weight functions, the stress intensity factors for a bimaterial interface crack can be obtained accurately and inexpensively by performing the sum of worklike products between the applied loads and the weight functions for the cracked bimaterial body under any loading conditions once the weight functions are explicitly predetermined. The same calculation can also be applied for the identical cracked bimaterial medium with different constraint conditions by including the self-equilibrium forces that contain both the external loads and the reaction forces induced at the constraint locations. Moreover, the physical interpretation of the weight functions can provide a guidance for damage tolerant design application.     

Weight function for an external circumferential semielliptical crack in a cylinder

In this paper, the weight function was extracted at the deepest point of a semielliptical circumferential crack. The crack is assumed to exist on the outer surface of the cylinder. For this purpose, the three‐dimensional finite element method was accomplished to specify two reference loads, which are indispensable for determining the weight function. The verification study confirms the accuracy of the derived weight function under prescribed mechanical loading on the crack surfaces. There is consistency among the solution results compared with those in the literature. The second part describes the application of the weight function for the thermal boundary conditions. Steady‐state thermal stress intensity factors are demonstrated using the weight function and presented as a closed‐form solution. The results were compared with the finite element data on the special case of thermal loading, and good agreement is obtained.     

Transient analysis of a piezoelectric strip with a permeable crack under anti-plane impact loads

The problem of a through permeable crack situated in the mid-plane of a piezoelectric strip is considered under anti-plane impact loads for two cases. The first is that the strip boundaries are free of stresses and of electric displacements, and the second is that the strip boundaries are clamped rigid electrodes. The method adopted is to reduce the mixed initial-boundary value problem, by using integral transform techniques, to dual integral equations, which are further transformed into a Fredholm integral equation of the second kind by introducing an auxiliary function. The dynamic stress intensity factor and energy release rate in the Laplace transform domain are obtained in explicit form in terms of the auxiliary function. Some numerical results for the dynamic stress intensity factor are presented graphically in the physical space by using numerical techniques for solving the resulting Fredholm integral equation and inverting Laplace transform.     

A numerical analysis of cracks emanating from a square hole in a rectangular plate under biaxial loads

This note concerns with stress intensity factors of cracks emanating from a square hole in rectangular plate under biaxial loads by means of the boundary element method which consists of the non-singular displacement discontinuity element presented by Crouch and Starfied and the crack tip displacement discontinuity elements proposed by the author. In the boundary element implementation the left or the right crack tip displacement discontinuity element is placed locally at corresponding left or right crack tip on top of the constant displacement discontinuity elements that cover the entire crack surface and the other boundary. The present numerical results illustrate that the present approach is very effective and accurate for calculating stress intensity factors of complicated cracks in a finite plate and can reveal the effect of the biaxial load and the cracked body geometry on stress intensity factors.     

Numerical simulations of cracked plate using XIGA under different loads and boundary conditions

This article deals with the numerical simulation of cracked plate using extended isogeometric analysis (XIGA) under different loads and boundary conditions. The plate formulation is done using first-order shear deformation theory. The crack faces are modeled by the Heaviside function, whereas the singularity in stress field at the crack tip is modeled by crack tip enrichment functions. The stress intensity factors for the cracked plate are numerically computed using a domain-based interaction integral. The results obtained by XIGA for the center and edge crack plate are compared with extended finite element method and/or literature results for different types of loads and boundary conditions.     

动荷载作用下含裂缝公路结构体的应力强度因子

以沈阳-大连高速公路为工程背景，基于弹性动力学理论，采用平面应变有限单元法，分析了车辆荷载对含裂缝路面体的动态作用，分析过程中，车辆荷载简化为正弦分布柔性荷载；路面结构体计算模型抽象为平面应变模型；路面结构体为弹性的连续介质，为了反映裂尖应力，位移场的奇异性和减少模型网格数，在裂尖环向设置了奇异单元。通过计算得到裂尖的位移场，由位移外插得到I-型应力强度因子随加载时间的变化规律。同时探讨了初始裂缝长度和公路结构材料阻尼比的变化对I-型应力强度因子分布规律的影响，为路面体的动态破坏研究提供了一定的理论参考。     

Closed‐form thermal stress intensity factors for an internal circumferential crack in a thick‐walled cylinder

In this paper the method of weight functions is employed to calculate the stress intensity factors for an internal circumferential crack in a thick‐walled cylinder. The pressurized cylinder is also subjected to convection cooling on the inner surface. Finite element method is used to determine an accurate weight function for the crack and a closed‐form thermal stress intensity factor with the aid of the weight function method is extracted. The influence of crack parameter and the heat transfer coefficient on the stress intensity factors are determined. Comparison of the results in the special cases with those cited in the literature and the finite element data shows that the results are in very good agreement.     

Weight function,stress intensity factor and crack opening displacement solutions to periodic collinear edge hole cracks

Periodic collinear edge hole cracks and arbitrary small cracks emanating from collinear holes, which are two typical multiple site damages occurred in the aircraft structures, are studied by using the weigh function method. An explicit closed form weight function for periodic edge hole cracks in an infinite sheet is obtained and further used to calculate the stress intensity factor and crack opening displacement for various loading cases. Compared to finite element method, the present weight function is accurate and highly efficient. The interactions of the holes and cracks on the stress intensity factor and crack opening displacement are quantitatively determined by using the present weight function. An approximate weight function method is also proposed for arbitrary small cracks emanating from multiple collinear holes. This method is very useful for calculating the stress intensity factor for arbitrary small cracks.     

Calculation of stress intensity factors in three dimensions by finite element methods

Finite element methods are used to calculate the stress intensity factors for three-dimensional geometries containing a number of depths of crack subjected to various loads. Special elements are used at the tip to represent the variation of the displacement with respect to the square root of the distance from the tip. The stress intensity factors are determined by comparison of the displacements in the special elements, by a method of virtual crack extensions, and, in one case, by an integral around the tip. With meshes containing between 50 and 100 quadratic isoparametric elements, results accurate to within 1 per cent or 4 per cent (depending on case) of known solutions are demonstrated.