Computation of certain double integrals that are characteristic of dynamic problems of the theory of cracks in a semi-infinite body

We propose a method of computing certain double integrals that are encountered in dynamic problems of the theory of elasticity for semi-infinite bodies with cracks. Application of this method makes it possible to reduce a system of boundary integral equations for the interaction of a crack with the boundary of the half-space to a form that contains only integrals over the regions occupied by the cracks. Translated fromMatematichni Metody i Fiziko-Mekhanichni Polya, Vol. 38, 1995.     

Elastic stress concentration near boundary defects in the case of a longitudinal shear deformation

We study the concentration of stresses due to two boundary defects located in an elastic half-space with stress-free boundary loaded at infinity with a constant shear load. The problem is reduced to solving singular integral equations for the cases in which the half-space contains defects of different types: two cracks, two inclusions, and a crack and an inclusion, whose solutions are sought by the method of mechanical quadratures. The interaction of the defects as they approach each other and the influence of their relative sizes are studied numerically. Translated fromDinamicheskie Sistemy, Vol. 11, 1992.     

Stress analysis in a sheet with multiple cracks

The stress field due to the presence of a Volterra dislocation in an isotropic elastic sheet is obtained. The stress components exhibit the familiar Cauchy type singularity at dislocation location. The solution is utilized to construct integral equations for elastic sheets weakened by multiple embedded or edge cracks. The cracks are perpendicular to the sheet boundary and applied traction is such that crack closing may not occur. The integral equations are solved numerically and stress intensity factors (SIFs) are determined on a crack edges.     

多裂纹问题计算分析的本征COD边界积分方程方法

针对多裂纹问题，若采用常规的数值求解技术，计算效率较低.为实现多裂纹问题的大规模数值模拟，建立了本征裂纹张开位移（crack opening displacement, COD）边界积分方程及其迭代算法，并引入Eshelby矩阵的定义，将多裂纹分为近场裂纹和远场裂纹来处理裂纹间的相互影响.以采用常单元作为离散单元的快速多极边界元法为参照，对提出的计算模型和迭代算法进行了数值验证.结果表明，本征COD边界积分方程方法在处理多裂纹问题时取得较大的改进，其计算效率显著高于传统的边界元法和快速多极边界元法.     

含曲线裂纹圆柱扭转问题的新边界元法

研究含曲线裂纹圆柱的Saint-Venant扭转,将问题化归为裂纹上边界积分方程的求解．利用裂纹尖端的奇异元和线性元插值模型,给出了扭转刚度和应力强度因子的边界元计算公式．对圆弧裂纹、曲折裂纹以及直线裂纹的典型问题进行了数值计算,并与用Gauss-Chebyshev求积法计算的直裂纹情形结果进行了比较,证明了方法的有效性和正确性．     

A study of the influence of the orientation of a planar surface crack in the shape of a limaçon of pascal on the stress concentration in a half-space

We study the stress state in the vicinity of a planar surface crack whose boundary is described by the limaçon of pascal. The problem is solved by a conformal mapping of the region occupied by the crack onto part of a disk in the plane. This makes it possible to apply a numerical-analytic method for solving systems of double singular integral equations of the mathematical theory of cracks. We present the graphs of the dependence of the stress intensity factor on the angular coordinate for cracks that are part of a limaçon of Pascal.     

Influence of initial stresses on fracture of composite materials containing interacting cracks

Considered in this study are the axially-symmetric problems of fracture of composite materials with interacting cracks, which are subjected to initial (residual) stresses acting along the cracks planes. An analytical approach within the framework of three-dimensional linearized mechanics of solids is used. Two geometric schemes of cracks location are studied: a circular crack is located parallel to the surface of a semi-infinite composite with initial stresses, and two parallel co-axial penny-shaped cracks are contained in an infinite composite material with initial stresses. The cracks are assumed to be under a normal or a radial shear load. Analysis involves reducing the problems to systems of second-kind Fredholm integral equations, where the solutions are identified with harmonic potential functions. Representations of the stress intensity factors near the cracks edges are obtained. These stress intensity factors are influenced by the initial stresses. The presence of the free boundary and the interaction between cracks has a significant effect on the stress intensity factors as well. The parameters of fracture for two types of composites (a laminar composite made of aluminum/boron/silicate glass with epoxy-maleic resin and a carbon/plastic composite with stochastic reinforcement by short ellipsoidal carbon fibers) are analyzed numerically. The dependence of the stress intensity factors on the initial stresses, physical-mechanical parameters of the composites, and the geometric parameters of the problem are investigated.     

压电材料中两个非对称平行裂纹的基本解

采用Schmidt方法分析压电材料中非对称平行的双可导通裂纹的断裂性能．利用Fourier变换使问题的求解转换为求解两对以裂纹面位移之差为未知变量的对偶积分方程．为了求解对偶积分方程,直接把裂纹面位移差函数展开成Jacobi多项式形式．最终得到了裂纹的应力强度因子与电位移强度因子之间的关系．数值结果表明,应力强度因子和电位移强度因子与裂纹间的距离、裂纹的几何尺寸有关;与不可导通裂纹有关结果相比,可导通裂纹的电位移强度因子远小于相应问题不可导通裂纹的电位移强度因子．同时可以发现裂纹间的“屏蔽”效应也在压电材料中出现．     

Analysis solution method for 3D planar crack problems of two-dimensional hexagonal quasicrystals with thermal effects

An analysis solution method (ASM) is proposed for analyzing arbitrarily shaped planar cracks in two-dimensional (2D) hexagonal quasicrystal (QC) media. The extended displacement discontinuity (EDD) boundary integral equations governing three-dimensional (3D) crack problems are transferred to simplified integral-differential forms by introducing some complex quantities. The proposed ASM is based on the analogy between these EDD boundary equations for 3D planar cracks problems of 2D hexagonal QCs and those in isotropic thermoelastic materials. Mixed model crack problems under combined normal, tangential and thermal loadings are considered in 2D hexagonal QC media. By virtue of ASM, the solutions to 3D planar crack problems under various types of loadings for 2D hexagonal QCs are formulated through comparison to the corresponding solutions of isotropic thermoelastic materials which have been studied intensively and extensively. As an application, analytical solutions of a penny-shaped crack subjected uniform distributed combined loadings are obtained. Especially, the analytical solutions to a penny-shaped crack subjected to the anti-symmetric uniform thermal loading are first derived for 2D hexagonal QCs. Numerical solutions obtained by EDD boundary element method provide a way to verify the validity of the presented formulation. The influences of phonon-phason coupling effect on fracture parameters of 2D hexagonal QCs are assessed.     

Fundamental two-dimensional solutions of electroelasticity for compound piezoceramic bodies

We construct fundamental solutions of the two-dimensional equations of electroelasticity for antiplanar strain of a piezoceramic space with an interphase defect. We study the orders of the powers of the singularities at the vertices of the defect for two cases: an interphase crack and a stiff fiber continuously joined to the upper half-space which has flaked off the lower half-space. The solution is constructed in closed form. Bibliography: 4 titles. Translated fromTeoreticheskaya i Prikladnaya Mekhanika, No. 28, 1998, pp. 62–67.     

Isoperimetric estimates of the solutions of a class of pseudodifferential equations and their application to crack problems

Isoperimetric estimates are obtained of solutions of boundary value problems for a class of pseudodifferential equations. This class of equations includes the equation of problems on plane normal discontinuity cracks located in a homogeneous linearly elastic space and an inhomogeneous space whose Young's modulus has a power-law dependence on the distance to the plane of the crack. As it applies to crack problems, the established inequalities yield, in particular, isoperimetric estimates of the maximum opening of the crack and its volume under arbitrary loads.     

Transient response of dissimilar piezoelectric layers with multiple interacting interface cracks

In this study, we examine the dynamic behavior of two bonded dissimilar piezoelectric layers containing multiple interfacial cracks subjected to electro-mechanical impact loading. The problem was formulated through Fourier transformation into singular integral equations in which the unknown variables are the jumps of displacement and electric potential across the crack surface in the Laplace transform domain. The resulting integral equations together with the corresponding single-valued conditions are solved numerically for the densities of electro-elastic dislocations on a crack surface. The dynamic field intensity factors and dynamic energy release rate (DERR) history are obtained for both permeable and impermeable crack. The stress field is also obtained for the interface crack under impact loads. The results show that the field intensity factors at the crack tips and dynamic energy release rate depend on the interfacial crack geometry, electromechanical coupling and the electric boundary conditions on the crack surface.     

夹层压电材料中垂直于界面的共线双裂纹动力学问题分析

采用Schmidt方法分析了在简谐反平面剪切波作用下,两个半空间夹层压电材料中的共线裂纹的动力学行为．压电材料层内裂纹垂直于界面,电边界条件假设为可导通．通过Fourier变换,使问题的求解转换为两对三重积分对偶方程．通过数值计算,给出了裂纹的几何尺寸、压电材料常数、入射波频率等对于应力强度因子的影响．结果表明,在不同的入射波频率范围,动力场将阻碍或促使压电材料内裂纹的扩展．与不可导通电边界条件相比,导通裂纹表面的电位移强度因子比不可导通裂纹的电位移强度因子要小许多．     

Stress concentration in an anisotropic half-plane with holes and cracks

On the basis of general representations of the generalized complex potentials for a multiconnected half-plane, which the authors have obtained, we solve problems for a multiconnected half-plane with holes and cracks when external forces or dies act on the boundary of the half-plane. Using conformal mapping for an ellipse and the method of least squares, we reduce these problems to solving a system of linear algebraic equations. For different anisotropic materials we give the results of studies of the stress distributions and the variation of the stress intensity factors for a half-plane with a crack in the case of tension at infinity, internal pressure on the edges of the crack, and the action of normal forces on the rectilinear boundary. Two figures, 2 tables. Bibliography: 2 titles. Translated fromTeoreticheskaya i Prikladnaya Mekhanika, No. 27, 1997, pp. 63–72.     

Global approaches for accurate loading analyses and crack path predictions at single and two-cracks systems

Based on the work by Eshelby, the path-independent J_k-, M-, L- and interaction- or I_k-integrals were introduced and applied to cracks for the accurate calculation of crack tip loading quantities. Applying the FE-method to solve boundary value problems with cracks, numerically inaccurate values are observed within the crack tip region affecting the accuracy of local approaches. Simulating crack paths, local approaches face further problems as cracks are running towards interfaces, internal boundaries or other crack faces. Within global approaches, path-independent integrals are calculated along remote contours far from the crack tip, essentially exploiting numerically reliable data requiring special treatment only for the near-tip crack faces. To provide path-independence, additional integrals along interfaces, internal boundaries and crack faces are necessary. In this paper, new global approaches of path-independent integrals are presented and applied to the calculation of crack paths at two-cracks systems. A second focus is directed to the accurate loading analysis and crack path prediction considering anisotropic properties and material interfaces. The numerical model provides crack paths which are in good agreement with those obtained from crack growth experiments. (© 2015 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Analysis of arbitrarily shaped planar cracks in two-dimensional hexagonal quasicrystals with thermal effects. Part I: Theoretical solutions

An extended displacement discontinuity (EDD) boundary integral equation method is proposed for analysis of arbitrarily shaped planar cracks in two-dimensional (2D) hexagonal quasicrystals (QCs) with thermal effects. The EDDs include the phonon and phason displacement discontinuities and the temperature discontinuity on the crack surface. Green's functions for unit point EDDs in an infinite three-dimensional medium of 2D hexagonal QC are derived using the Hankel transform method. Based on the Green's functions and the superposition theorem, the EDD boundary integral equations for an arbitrarily shaped planar crack in an infinite 2D hexagonal QC body are established. Using the EDD boundary integral equation method, the asymptotic behavior along the crack front is studied and the classical singular index of 1/2 is obtained at the crack edge. The extended stress intensity factors are expressed in terms of the EDDs across crack surfaces. Finally, the energy release rate is obtained using the definitions of the stress intensity factors.     

三维横观各向同性介质界面裂纹的边界积分方程方法

基于两相三维横观各向同性介质的基本解和Somigliana恒等式,对三维横观各向同性介质中的任意形状的平片界面裂纹,以裂纹面上的不连续位移为待求参量建立了超奇异积分_微分方程,界面平行于横观各向同性面.根据发散积分的有限部积分理论,应用积分方程方法研究得到裂纹前沿的位移和应力场的表达式、奇性指数以及应力强度因子的不连续位移表达式.在非震荡情形下,超奇异积分_微分方程退化为超奇异积分方程,与均匀介质的超奇异积分方程形式完全相同.     

A modified scaled boundary finite element method for dynamic response of a discontinuous layered half-space

In this paper, a modified scaled boundary finite element method is proposed to deal with the dynamic analysis of a discontinuous layered half-space. In order to describe the geometry of discontinuous layered half-space exactly, splicing lines, rather than a point, are chosen as the scaling center. Based on the modified scaled boundary transformation of the geometry, the Galerkin's weighted residual technique is applied to obtain the corresponding scaled boundary finite element equations in displacement. Then a modified version of dimensionless frequency is defined, and the governing first-order partial differential equations in dynamic stiffness with respect to the excitation frequency are obtained. The global stiffness is obtained by adding the dynamic stiffness of the interior domain calculated by a standard finite element method, and the dynamic stiffness of far field is calculated by the proposed method. The comparison of two existing solutions for a horizontal layered half-space confirms the accuracy and efficiency of the proposed approach. Finally, the dynamic response of a discontinuous layered half-space due to vertical uniform strip loadings is investigated.     

Parabolic Equations with Partially BMO Coefficients and Boundary Value Problems in Sobolev Spaces with Mixed Norms

Second order parabolic equations in Sobolev spaces with mixed norms are studied. The leading coefficients (except a ¹¹) are measurable in time and one spatial variable, and have small BMO (bounded mean oscillation) semi-norms as functions of the other spatial variables. The coefficient a ¹¹ is measurable in time and have a small BMO semi-norm in the spatial variables. The unique solvability of equations in the whole space is established. Then this result is applied to solving Dirichlet and Neumann boundary value problems for parabolic equations defined on a half-space or on a bounded domain.     

Radiation of short waves by a pair of plates backed by a semi-infinite set of ribbed stiffeners

The following acoustic diffraction problem is considered. The upper half-space is filled by an acoustic medium. Two semi-infinite thin plates are situated on the boundary of the upper half-space. One of the plates is backed by a semi-infinite periodical set of ribbed stiffeners. The source of an acoustic field is positioned on the other plate. The problem is reduced to an infinite system of linear algebraic equations. Such a system can be solved in a shortwave approximation. An expression for the acoustic field at a large distance from the junction of the plates is obtained. Translated fromZapiski Nauchnykh Seminarov POMI, Vol. 210, 1994, pp. 38–46. Translated by A. V. Badanin.