Kinked crack analysis by a hybridized boundary element/boundary collocation method

In this research a two dimensional displacement discontinuity method (which is a kind of indirect boundary element method) using higher order elements (i.e. a source element with a cubic variation of displacement discontinuities having four sub-elements) is used to obtain the displacement discontinuities along each boundary element. In this paper, three kinds of the higher order boundary elements are used: the ordinary elements, the kink elements and the special crack tip elements.The boundary collocation technique is used for the calculation of the displacement discontinuities at the center of each sub-elements. Again a special boundary collocation technique is used to treat the kinked source elements occur in the crack analysis. Considering the two source elements (each having four sub-elements) joined at a corner (kink point). The collocation points in the cubic element model which are outside of the kink point are moved to the crack kink then the displacement discontinuities on the left and right sides of the kink are calculated. The displacement discontinuities of the kink point are obtained by averaging the corresponding values of its left and right sides. The special crack tip elements are also treated by the boundary displacement collocation technique considering the singularity variation of the displacements and stresses near the crack tip. Some simple example problems are solved numerically by the proposed method. The numerical results are compared with the corresponding results obtained by the previous methods cited in the literature. This comparison shows a very good agreement between the results and verify the accuracy and validity of the proposed method.     

Stress intensity factors of a plate with two cracks emanating from an arbitrary hole

In this paper, Muskhelishvili complex function theory and boundary collocation method are used to calculate the stress intensity factors (SIF) of a plate with two cracks emanating from an arbitrary hole. The calculated examples include a circular, elliptical, rectangular, or rhombic hole in a plate. The principle and procedure by the method is not only rather simple, but also has good accuracy. The SIF values calculated compare very favorably with the existing solutions. At the same time,the method can be used for different finite plate with two cracks emanating from a hole with more complex geometrical and loading conditions. It is an effective unified approach for this kind of fracture problems.     

Torsion of a thick-walled cylinder with an external crack: boundary collocation method

The boundary collocation method is used to obtain the torsional rigidity and Mode III stress intensity factor of a thick-walled cylinder with an external radial crack. When the internal radius of the cylinder is very small, the results agree well with those obtained previously from other methods for an edge crack in a solid cylindrical bar. The present method is shown to be expedient when applied to obtain results for different ratios of the internal and external radii of the cracked cylinder.     

The anti-plane solution for the edge cracks originating from an arbitrary hole in a piezoelectric material

In this paper, a general and simple way was found to solve the problem of an arbitrary hole with edge cracks in transversely isotropic piezoelectric materials based on the complex variable method and the method of numerical conformal mapping. Firstly, the approximate mapping function which maps the outside of the arbitrary hole and the cracks into the outside of a circular hole is derived after a series of conformal mapping process. Secondly, based on the assumption that the surface of the cracks and hole is electrically impermeable and traction-free, the approximate expressions for the complex potential, fields intensity factors and energy release rates are presented, respectively. Thirdly, under the in-plane electric loading together with the out-plane mechanical loading, the influences of the hole size, crack length and mechanical/electric loading on the fields intensity factors and energy release rates are analyzed. Finally, some particular holes with edge cracks are studied in numerical analysis. The result shows that, the mechanical loading always promotes crack growth, while the electric loading may retard crack growth.     

The mode III cracks emanating from an elliptical hole in the piezo-electro-magneto-elastic materials

The mode III crack problem in a medium possessing coupled electro-magneto-elastic is considered. Two asymmetrical edge cracks emanate from an elliptical hole. Combined stress, electric and magnetic loads are considered. The elliptical hole and cracks are assumed to be either magneto-electrically impermeable or permeable. The closed-form solution for stress, electric and magnetic intensity factors are derived explicitly. Also the solution for energy release rate is given in closed form. The solution is based on the complex variable method combined with the method of conformal mapping. Numerical computations are given to illustrate the effect of variable geometrical and material parameters on stress, electric and magnetic intensity factors and energy release rate.     

The behavior of two non-symmetrical permeable cracks emanating from an elliptical hole in a piezoelectric solid

Based on the Stroh-type formalism, we present a concise analytic method to solve the problem of complicated defects in piezoelectric materials. Using this method and the technique of conformal mapping, the problem of two non-symmetrical collinear cracks emanating from an elliptical hole in a piezoelectric solid is investigated under remotely uniform in-plane electric loading and anti-plane mechanical loading. The exact solutions of the field intensity factors and the energy release rate are presented in closed-form under the permeable electric boundary condition. With the variation of the geometrical parameters, the present results can be reduced to the well-known results of a mode-III crack in piezoelectric materials. Moreover, new special models used for simulating more practical defects in a piezoelectric solid are obtained, such as two symmetrical edge cracks and single edge crack emanating from an elliptical hole or circular hole, T-shaped crack, cross-shaped crack, and semi-infinite plane with an edge crack. Numerical results are then presented to reveal the effects of geometrical parameters and the applied mechanical loading on the field intensity factors and the energy release rate.     

Fracture calculation of bending plates by boundary collocation method

IntroductionPlatesarecommonstructuresinaerospace ,civil,mechanicalandchemicalengineering .Duringtheproductionanduseofthestructures,flawsandcracksmayappearintheplates.Inordertoassurethesafetyofplatestructures,toavoidcleavagefractureandfatiguefailure,andtodesignsafetyandmakerationaluseofmaterials,itisnecessarytodofractureanalysisfortheplatestructures.Fractureproblemofabendingplateismorecomplexthanthatoftheplaneelasticity[1,2 ].Firstofall,therearedifferenttheoriesforplatebending ,suchasKirchhoff…     

THE BOUNDARY INTEGRAL-VARIATIONAL THEOREMS OF AN ARBITRARY ELEMENT IN THE SOLID—COMPUTE THE ENERGY RELEASE RATE OF AN ARBITRARY CRACK EXTENSION

Based on the solid mechanics of the discrete form and its variational principles proposed by Niu^[1–2], this paper puts forward four kinds of boundary integral-variational theorems of an arbitrary element. In the course of the fracture analysis, they can be used to compute the energy release rate along the normal direction of the crack boundary. When there is a hole in the solid, and whether there are given surface forces on the hole boundary or not they can be used to compute the variation of the energy along the normal direction of the hole boundary. In the course of the discrete analysis, they can be used to establish the discrete equations, so that the values of the unknown functions are solved. At the same time, from this paper we know that the J-integral proposed by Rice^[3] represents an integral to be independent of a path imperfectly.     

A boundary collocation method for the solution of a flow problem in a complex three-dimensional porous medium

A flow problem in a complex three-dimensional domain with a free surface and mixed-type boundary conditions is solved by the boundary collocation method. The solution is expressed as a combination of source functions distributed all around the domain close to the boundary, plus a special basis function to take care of a corner singularity. The resulting procedure is compared with the boundary integral elements method and is found to be simpler and more flexible to implement and faster to compute.     

A numerical analysis of cracks emanating from an elliptical hole in a 2-D elasticity plate

This paper is concerned with the stress intensity factors (SIFs) of cracks emanating from an elliptical hole in an infinite or a finite plate under biaxial loads by using a boundary element method, which consists of the non-singular displacement discontinuity element presented by Crouch and Starfield and the crack-tip displacement discontinuity elements due to the author. In the boundary element implementation the left or the right crack-tip 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. A few numerical examples are included to show that the present approach is very efficient and accurate for the calculating the SIFs of crack problems in an infinite or a finite plate. The present numerical results of cracks emanating from an elliptical hole under biaxial loads can reveal the effect of the elliptical aspect ratio and the transverse load on the SIFs.     

A boundary integral equation method for the solution of a class of crack problems

A boundary integral equation method for the solution of an important class of crack problems in elasticity is outlined. The method is applicable for deformations in which the crack faces remain in contact. Specific numerical examples are considered to illustrate the application of the method.     

Fatigue growth modeling of cracks emanating from a circular hole in infinite plate

This paper presents a numerical approach of fatigue growth analysis of cracks emanating from a hole in infinite elastic plate subjected to remote loads. It involves a generation of Bueckner’s principle and a hybrid displacement discontinuity method (a boundary element method) proposed recently by the senior author of the paper. Because of an intrinsic feature of the boundary element method, a general crack growth problem can be solved in a single region formulation. In the numerical simulation, for each increment of crack extension, remeshing of existing boundaries is not necessary. Crack extension is modeled conveniently by adding new boundary elements on the incremental crack extension to the previous crack boundaries. As an example, fatigue growth process of an inclined crack in an infinite plate under uniaxial cycle load is modeled to illustrate the effectiveness of the numerical approach. In addition, fatigue growth of cracks emanating from a circular hole in infinite elastic plate subjected to remote loads is investigated by using the numerical approach. Many numerical results are given     

Complex variable function method for the scattering of plane waves by an arbitrary hole in a porous medium

Based on the complex variable function method, a new approach for solving the scattering of plane elastic waves by a hole with an arbitrary configuration embedded in an infinite poroelastic medium is developed in the paper. The poroelastic medium is described by Biot's theory. By introducing three potentials, the governing equations for Biot's theory are reduced to three Helmholtz equations for the three potentials. The series solutions of the Helmholtz equations are obtained by the wave function expansion method. Through the conformal mapping method, the arbitrary hole in the physical plane is mapped into a unit circle in the image plane. Integration of the boundary conditions along the unit circle in the image plane yields the algebraic equations for the coefficients of the series solutions. Numerical solution of the resulting algebraic equations yields the displacements, the stresses and the pore pressure for the porous medium. In order to demonstrate the proposed approach, some numerical results are given in the paper.     

Fracture instability of reinforced panel with cracks emanating from hole

Initiation of failure by yielding and/or fracture depends on the magnitude of the distortion and dilatation of material elements. According to the strain energy density theory (SED), failure is assumed to initiate at the site of the local maximum of maxima [(dW/dV)^max_max]_L by yielding and the maximum of minima [(dW/dV)^max_min]_L by fracture. The fracture is assumed to start from point L where [(dW/dV)^max_min]_L appears and tends toward G where the global maximum of dW/dV minima appears, denoted by [(dW/dV)^max_min]_G. The distance l between L and G along the anticipated crack trajectory is an indication of failure instability of the system by fracture. If l is sufficiently large and [(dW/dV)^max_min]_L exceeds the threshold, fracture initiation could lead to global failure. The local and global failure instability of a composite structural component is studied by application of the strain energy density theory. The depicted configuration is that of a panel with a circular hole reinforced by two side strips made of different material. The case of two symmetric cracks emanating from the hole and normal to the applied uniaxial tensile stress is also analyzed. Results are displayed graphically to illustrate the geometry and dissimilar material properties influence the fracture instability behavior of the two examples.     

Numerical solution of multiple hole problem by using boundary integral equation

This paper studies a numerical solution of multiple hole problem by using a boundary integral equation. The studied problem can be considered as a supposition of many single hole problems. After considering the interaction among holes, an algebraic equation is formulated, which is then solved by using an iteration technique. The hoop stress around holes can be finally determined. One numerical example is provided to check its accuracy.     

A crack emanating from the tip of bonded dissimilar materials

1.IntroductionBondedstructurescanbewidelyfoundindifferentareas,suchasweldedpressurevessels,reinforcedconcretemembers,groutingsoftfoundationandsolidrocketpropellants.etc.Thestudy'ofthesingularityofbondedstructuresisespeciallyimp6rtantnotonlyforsafedesignbutalsoforconstruction.DeinpseyandSinclairl61studiedthegeneraleaseofN-materialcompositewedges,andprovedthat'ingeneralthereexistsstresssingularitynearthetipofthewedges,andtheorderofthesingularitydependsontheelasticconstantsandthelocalgeometry.Fo…     

Exact solutions for anti-plane problem of two asymmetrical edge cracks emanating from an elliptical hole in a piezoelectric material

Based on the complex variable method and the technique of conformal mapping, the anti-plane problem of two asymmetrical edge cracks emanating from an elliptical hole in a piezoelectric material is studied. The exact solutions of field intensity factors and energy release rate are presented in closed-form with the assumption that the surfaces of the cracks and the elliptical hole are electrically impermeable. With the variation of the hole-size and the crack length, the present results can be reduced to the cases of two symmetrical edge cracks and a single edge crack emanating from a circular hole given by Wang and Gao [Wang, Y.J., Gao, C.F., 2008. The mode III cracks originating from the edge of a circular hole in a piezoelectric solid. International Journal of Solids and Structures 45, 4590–4599]. Moreover, new models used for simulating more practical defects in a piezoelectric solid are obtained, such as two symmetrical edge cracks and a single edge crack emanating from an elliptical hole, two asymmetrical edge cracks emanating from a circular hole, T-shaped crack, cross-shaped crack and semi-infinite plane with an edge crack. Numerical examples are then conducted to reveal the effects of the hole-size and the crack length on the field intensity factors and the energy release rate.     

Kinking of a crack emanating from the free surface of a beam under pure bending

A perturbation technique developed by Karihaloo et al. is employed to obtain the stress intensity factors at the tip of a kinking crack that emanates from the free surface of a beam under pure bending. Under the condition that the kink extends in the direction of vanishing K₁₁ the crack path is obtained as well as a path stability condition. From conditions on K₁ a material parameter r^* akin to that of Ramulu and Kobayashi's r_c is obtained. By analysis of the slope of the kinking crack a stability condition is obtained corroborating the stability condition from consideration of vanishing K₁₁. It is shown that for a beam in pure bending the nonsingular remote stress term T must be greater than some positive critical value for kinking to occur confirming the results of Sayir and Schindler.     

Viscoelastic crack analysis by the boundary integral equation method

Summary The linear viscoelastic three-dimensional crack problem is analyzed by combining the correspondence principle and the boundary integral equation method. In a general crack analysis the usual boundary integral equations lead to a nonunique formulation of the problem, because they do not involve information about the loading on the crack surface. Here, the boundary integro-differential equations are applied to the numerical calculation of the crack opening displacement of a penny-shaped crack in an infinite linear viscoelastic body. Moreover, the influence of several parameters of the three-parameter viscoelastic model on the crack opening displacement and the incubation time is shown.

---

Viskoelastische Rißanalyse durch Randintegralgleichungen

Übersicht Das linear viskoelastisch räumliche Rißproblem wird mit Hilfe einer Kombination von Korrespondenzprinzip und Randintegralgleichungsverfahren gelöst. In einer allgemeinen Rißanalyse führen die üblichen Randintegrale zu einer nicht eindeutigen Formulierung dieses Problems, weil die Angaben über Belastung und Rißoberfläche fehlen. Das Randintegralgleichungsverfahren wird für die numerische Berechnung der Rißerweiterung eines münzförmigen Risses in einem unendlich linear viskoelastischen Körper angewendet. Weiterhin wird der Einfluß von verschiedenen Parametern des räumlich viskoelastischen Modells auf die Rißerweiterung und die Inkubationszeit gezeigt.