Stress-strain state of an anisotropic plate with curved cracks and thin rigid inclusions

We solve the bending problem for an anisotropic plate with flaws like smooth curved nonoverlapping through cracks and rigid inclusions. The problem is solved by the method of Lekhnitskii complex potentials specified as Cauchy type integrals over the flaw contours with an unknown integrand density function. We use the Sokhotskii—Plemelj formulas to reduce the boundary-value problem to a system of singular integral equations with the additional conditions that the displacements in the plate are single-valued when going around the cut contours and the equilibrium conditions for stress-free rigid inclusions. After the singular integrals are approximated by the Gauss-Chebyshev quadrature formulas, the problem is reduced to solving a system of linear algebraic equations. We study the local stress distribution near flaw tips. We analyze the mutual influence of flaws on the stress distribution character near their vertices and compare the well-known solutions for isotropic plates with the solutions obtained by passing to the limit in the anisotropy parameters (“weakly anisotropic material”) and by using the method proposed here.     

Bueckner's work conjugate integrals and weight functions for a crack in anisotropic solids

The Bueckner work conjugate integrals are studied for cracks in anisotropic elastic solids. The difficulties in separating Lekhnitskii's two complex arguments involved in the integrals are overcome and explicit functional representations of the integrals are given for several typical cases. It is found that the pseudoorthogonal property of the eigenfunction expansion forms presented previously for isotropic cases, isotropic bimaterials, and orthotropic cases, are proved to be also valid in the present case of anisotropic material. Finally, Some useful path-independent integrals and weight functions are proposed. The project supported by the National Natural Science Foundation of China (19891180) and Doctorate Foundation of Xi'an Jiaotong University     

Analysis of three-dimensional edge cracks under tensile loading

Three-dimensional edge cracks are analyzed using the Self-Similar Crack Expansion (SSCE) method with a boundary integral equation technique. The boundary integral equations for surface cracks in a half space are presented based on a half space Green's function (Mindlin, 1936). By using the SSCE method, the stress intensity factors are determined by the crack-opening displacement over the crack surface. In discrete boundary integral equations, the regular and singular integrals on the crack surface elements are evaluated by an analytical method, and the closed form expressions of the integrals are given for subsurface cracks and edge crakcs. This globally numerical and locally analytical method improves the solution accuracy and computational effort. Numerical results for edge cracks under tensile loading with various geometries, such as rectangular cracks, elliptical cracks, and semi-circular cracks, are presented using the SSCE method. Results for stress intensity factors of those surface breaking cracks are in good agreement with other numerical and analytical solutions.     

A continuum damage model for piezoelectric materials

In this paper, a constitutive model is proposed for piezoelectric material solids containing distributed cracks. The model is formulated in a framework of continuum damage mechanics using second rank tensors as internal variables. The Helrnhotlz free energy of piezoelectric mate- rials with damage is then expressed as a polynomial including the transformed strains, the electric field vector and the tensorial damage variables by using the integrity bases restricted by the initial orthotropic symmetry of the material. By using the Talreja＇s tensor valued internal state damage variables as well as the Helrnhotlz free energy of the piezoelectric material, the constitutive relations of piezoelectric materials with damage are derived. The model is applied to a special case of piezoelectric plate with transverse matrix cracks. With the Kirchhoff hypothesis of plate, the free vibration equations of the piezoelectric rectangular plate considering damage is established. By using Galerkin method, the equations are solved. Numerical results show the effect of the damage on the free vibration of the piezoelectric plate under the close-circuit condition, and the present results are compared with those of the three-dimensional theory.     

A Self-Similar Crack Expansion method for three-dimensional cracks in an infinite or semi-infinite medium

The Self-Similar Crack Expansion (SSCE) method is used to calculate stress intensity factors for three-dimensional cracks in an infinite medium or semi-infinite medium by the boundary integral element technique, whereby, the stress intensity factors at crack tips are determined by calculating the crack-opening displacements over the crack surface. For elements on the crack surface, regular integrals and singular integrals are precisely evaluated based on closed form expressions, which improves the accuracy. Examples show that this method yields very accurate results for stress intensity factors of penny-shaped cracks and elliptical cracks in the full space, with errors of less than 1% as compared with analytical solutions. The stress intensity factors of subsurface cracks are in good agreement with other analytical solutions.     

Further investigation of subinterface cracks

Summary Further investigation of subinterface cracks in bimaterial solids is presented. The traction method is proposed permitting easy calculation of the stress intensity factors of the cracks. The elastic T-terms of the cracks are determined. The J ₁ and J ₂ integrals are analyzed for a contour enclosing all the cracks in the global coordinate system x,y, where the x-axis is parallel to the interface. Numerical examples are given, and the results are presented for two kinds of material combinations, Cu/Al₂O₃ and Ni/MgO. Accepted for publication 24 November 1996     

SOLUTION OF MULTIPLE CRACKS IN A FINITE PLATE OF AN ELASTIC ISOTROPIC MATERIAL WITH THE DISTRIBUTED DISLOCATION METHOD

This paper presents a numerical solution to model multiple cracks in a finite plate of an elastic isotropic material. Both the boundaries and the cracks are modeled by distributed dislocations. This method results in a system of singular integral equations with Cauchy kernels which can be solved by Gauss-Chebyshev quadrature method. Four examples are provided to assess the capability of this method.     

Surface cracks in a plate of finite width under extension or bending

In this paper the problem of a finite plate containing collinear surface cracks is considered. The problem is solved by using the line spring model with plane elasticity and Reissner's plate theory. The main purpose of the study is to investigate the effect of interaction between two cracks or between cracks and stress-free plate boundaries on the stress intensity factors and to provide extensive numerical results which may be useful in applications. First, some sample results are obtained and are compared with the existing finite element results. Then the problem is solved for a single (internal) crack, two collinear cracks and two corner cracks for wide range of relative dimensions. Particularly in corner cracks the agreement with the finite element solution is surprisingly very good. The results are obtained for semielliptic and rectangular crack profiles which may, in practice, correspond to two limiting cases of the actual profile of a subcritically growing surface crack.     

THERMAL FRACTURE OF FUNCTIONALLY GRADED PLATE WITH PARALLEL SURFACE CRACKS

This work examines the fracture behavior of a functionally graded material （FGM） plate containing parallel surface cracks with alternating lengths subjected to a thermal shock. The thermal stress intensity factors （TSIFs） at the tips of long and short cracks are calculated using a singular integral equation technique. The critical thermal shock △Tc that causes crack initiation is calculated using a stress intensity factor criterion. Numerical examples of TSIFs and △Tc for an Al2O3/Si3N4 FGM plate are presented to illustrate the effects of thermal property gradation, crack spacing and crack length ratio on the TSIFs and △Tc. It is found that for a given crack length ratio, the TSIFs at the tips of both long and short cracks can be reduced significantly and △Tc can be enhanced by introducing appropriate material gradation. The TSIFs also decrease dramatically with a decrease in crack spacing. The TSIF at the tips of short cracks may be higher than that for the long cracks under certain crack geometry conditions. Hence, the short cracks instead of long cracks may first start to grow under the thermal shock loading.     

Natural Vibrations of a Composite Plate Strip with Macrocracks

The fundamental frequency of a clamped composite plate strip with either one free macrocrack or two parallel free cracks is determined using the exact equations of the theory of elasticity. The plate is assumed made of an anisotropic orthotropic material with effective mechanical properties. The problem is solved by the finite-element method. Numerical examples are given to demonstrate how the problem parameters influence the fundamental frequency.     

非均质中厚板的无网格LRPIM动力学分析

用局部加权残值法建立了非均质中厚板的局部径向点插值离散系统方程,采用无网格局部径向点插值法分析了非均质中厚板的自由振动和强迫振动问题。用径向基函数耦合多项式基函数来近似试函数,用四次样条函数做为加权残值法中的权函数。所构造的形函数具有Kronecker delta性质,可以很方便地施加本质边界条件。该方法不需要任何形式的网格划分,所有的积分都在规则形状的子域及其边界上进行。在计算过程中,取积分中的高斯点的材料参数来模拟问题域材料特性的变化。计算结果表明,利用该方法计算非均质中厚板的自由振动和强迫振动问题可以得到具有较高精度的解。     

Dynamic stress intensity factors around two rectangular cracks in an infinite elastic plate under impact load

A dynamic problem for two equal rectangular cracks in an infinite elastic plate is considered. The two cracks are placed perpendicular to the plane surfaces of the plate. An incoming shock tensile stress is returned by the cracks. In the Laplace transform domain, the boundary conditions at the two sides of the plate are satisfied using the Fourier transform technique. The mixed boundary conditions are reduced to dual integral equations. Crack displacement is expanded in a series of functions which are zero outside of the cracks. The unknown coefficients in the series are determined by the Schmidt method. The stress intensity factors are defined in the Laplace transform domain and these are inverted using a numerical method.     

On Cauchy singular integrals and stress intensity factors for 2D mode I cracks in an infinite solid

This paper comments on some of the different numerical techniques commonly employed in evaluating Cauchy singular integrals of the first kind; e.g. as pertaining to 2D through cracks in a brittle material undergoing Mode I loading. In addition, a different more direct method is proposed here. Also, two different ways to calculate the stress intensity factor (K_I) are contrasted. The accuracy attained by the different methods in calculating K_I, and the factors affecting the calculation, are compared. Finally, comments on calculating the stress field of a 2D crack and important considerations are presented.     

Electromagnetoelastic problem for an infinite plate with known electrical potentials at hole boundaries

The plane electromagnetoelastic problem for an infinite plate with cracks and holes for known electrical potentials at the boundaries of the cracks and holes is solved using generalized complex potentials and conformal mapping, series, and least squares methods. The electromagnetoelastic state and variation of the stress intensity coefficient versus the geometric and physicomechanical properties of the material was studied for plates with two circular holes, a crack and a circular hole, or two cracks.     

An eight noded composite plate element for the dynamic analysis of spatial mechanism systems

The development of a shear-deformable laminated plate element, based on the Mindlin plate theory, for use in large reference displacement analysis is presented. The element is sufficiently general to accept an arbitrary number of layers and an arbitrary number of orthotrophic material property sets. Coordinate mapping is utilized so that non-rectangular elements may be modeled. The Gauss quadrature method of numerical integration is utilized to evaluate volume integrals. A comparative study is done on the use of full Gauss quadrature, reduced Gauss quadrature, mixed Gauss quadrature, and closed form integration techniques for the element. Dynamic analysis is performed on the RSSR (Revolute-Spherical-Spherical-Revolute) mechanism, with the coupler modeled as a flexible plate. The results indicate the differences in the dynamic response of the transverse shear deformable eight-noded element as compared to a four-noded plate element. Dynamically induced stresses are examined, with the results indicating that the primary deformation mode of the eight-noded Mindlin plate model being bending.     

Three-dimensional mixed-mode stress intensity factors for cracks in functionally graded materials using enriched finite elements

Three-dimensional enriched finite elements are used to compute mixed-mode stress intensity factors (SIFs) for three-dimensional cracks in elastic functionally graded materials (FGMs) that are subject to general mixed-mode loading and constraint conditions. The method, which advantageously does not require special mesh configuration/modifications and post-processing of finite element results, is an enhancement of previous developments applied so far on isotropic homogeneous and isotropic interface cracks. The spatial variation of FGM material properties is taken into account at the level of element integration points. To validate the developed method, two- and three-dimensional mixed-mode fracture problems are selected from the literature for comparison. Two-dimensional cases include: inclined central crack in a large FGM medium under uniform tensile strain and stress loadings, a slanted crack in a finite-size FGM plate under exponentially varying tensile stress loading and an edge crack in a finite-size plate under shear traction load. The three-dimensional example models a deflected surface crack in a finite-size FGM plate under uniform tensile stress loading. Comparisons between current results and those from analytical and other numerical methods yield good agreement. Thus, it is concluded that the developed three-dimensional enriched finite elements are capable of accurately computing mixed-mode fracture parameters for cracks in FGMs.     

三维间断位移法及强奇异和超奇异积分的处理方法

从积分方程Somigliana等式出发,导出三维状态下单位位错集度的基本解.在此基础上,建立了边界积分方程,并给出了其离散形式.对强奇异和超奇异积分,采用了Hadamard定义的有限部分积分来处理.最后,给出了计算裂纹应力强度因子的算例,并与解析解进行了比较,证实了该方法的有效性.     

混合边界多裂隙体在简谐变温场作用下的断裂问题（张开型）

给出了无限平面作用有简谐变化的点热源时的位移场、应力场基本解、用间接法构造出混合边值多裂隙体在简谐变温场作用下的热断裂问题的边界积分方程,并离散求解.数值结果表明,该方法求解多裂隙体的简谐热断裂问题精度好,计算工作量少.文中计算了含边界裂缝的平板、含三条平行裂缝的平板在简谐变温场作用下缝端应力强度因子的变化过程,并与实验结果进行了比较,两者吻合良好.     

Size effect on cracked functional composite micro-plates by an XIGA-based effective approach

Failure of structures and their components is one of major problems in engineering. Studies on mechanical behavior of functionally graded (FG) microplates with defects or cracks by effective numerical methods are rarely reported in literature. In this paper, an effective numerical model is derived based on extended isogeometric analysis (XIGA) for assessment of vibration and buckling of FG microplates with cracks. Based on the modified couple stress theory, the non-classical theory of Reissner–Mindlin plate is extended to capture microstructure, and thus, the size effect. In such theory, possessing C¹-continuity is straightforward with the high-order continuity of non-uniform rational B-spline. Due to the use of enrichments in XIGA, crack geometry is independent of the computational mesh. Numerical examples are performed to illustrate the effects of microplate aspect ratio, crack length, internal material length scale parameter, material distribution, and boundary condition on the mechanical responses of cracked FG microplates. The obtained results are compared with reference solutions and that shows that the frequency and buckling loads increases with decreasing the size of FG microplates and crack length. The convergence of the present method is also studied.     

A simple method for measuring surface strains around cracks

A simple system has been developed to measure surface strains that occur during in situ deformation of mechanical test specimens. The system uses photolithographically deposited displacement markers and computer image recognition routines to determine in-plane displacements and strains from digital images. The strain calculating routines are integrated into a simple mouse-driven software package that facilitates the transformation from digital images to useful strain field information. Additional routines have been developed to determine crack tip stress fields and J integrals. Crack tip stress intensities have been calculated from strain maps obtained for traction-free cracks in stainless steel. The J integrals were found to be independent of contour and consistent with applied stress intensities. Crack tip stress intensities were calculated for bridged cracks in lamellar TiAl. The toughening effect of the bridging zones was determined by including the bridged region in the contours. Resistance curves generated from strain maps were consistent with those measured during mechanical testing.