双材料中平片裂纹问题的超奇异积分方程解法

利用三维断裂力学的超奇异积分方程方法,对双材料空间中重直于界面的平片裂纹Ⅰ型问题进行了研究。首先根据双材料空间的弹性力学基本解,使用边界积分方程方法,在有限部积分的意义下导出了以裂纹面位罗间断为未知函数的超奇异积分方程,并为其建立了数值法。在此基础上,讨论了用裂纹面位移问题计算应力强度因子的方法。最后用此计算了几个典型的Ⅰ型下片裂纹问题的应力强度因子,其数值结果令人满意。     

压电材料中三维I型断裂力学分析

讨论了不可导通情况下三维横观各向同性压电材料中受拉伸和电载荷作用的平片裂纹I型断裂力学问题. 使用有限部积分概念，从三维线性压电理论出发，严格得到了一组以裂纹面位移间断和电势间断为未知变量的超奇异积分方程组；应用二维超奇异积分的主部分析法，从理论上分析得到了裂纹前沿应力和电势奇性指数以及应力和电位移奇性场，从而找到了以裂纹面位移间断和电势间断表示的应力和电位移强度因子、能量释放率表达式；为所得到的超奇异积分方程组建立了数值法，并用此计算了若干典型的平片裂纹问题，数值结果令人满意.     

压电材料中三维I型断裂力学分析

讨论了不可导通情况下三维横观各向刚性压电材料中受拉伸和电载荷作用的平片裂纹Ⅰ型断裂力学问题．使用自限部分概念，从二维线性压电理论出发，严格得到了一组以裂纹面位移间断和电势间断为未知变量的超奇异积分方程组；应用二维超奇异积分的主部分析法，从理论上分析得到了裂纹前沿应力和电势奇性指数以及应力和电位移奇性场，从而找到了以裂纹面位移间断和电势间断表示的应力和电位移强度因子、能量释放率表达式；为所得到的超奇异积分方程组建立了数值法，并用此计算了若干典型的平片裂纹问题，数值结果令人满意．     

三维界面裂纹的奇性应力场和应力强度因子分析

使用有限部积分概念和极限方法得到了三维平片界面裂纹的超奇异积分－微分方程组后,进一步利用二维超奇异积分主部分析方法,对裂纹前沿的应力场作了理论分析,并获得了其奇性应力场和裂纹面位移间断表示复位应力强度因子的精确表达式,为三维平片界面裂纹的超奇异积分－微分方程组的求解建立了数值方法,并分析了界面椭圆平片裂纹问题,和现有解比较,所得数值结果令人满意。     

三维有限体中两平行裂纹干扰问题的有限部积分与边界元法

利用Ｓｏｍｉｇｌｉａｎａ公式及有限部积分的概念，导出含两平行平片裂纹三维有限体裂纹干扰问题的超奇异积分方程组，联合使用有限部积分与边界元法，建立了数值求解方法，为提高数值计算结果的精度，在裂纹前疝附近单元，采用平方根位移模型，并在此基础雌出直接计算应力强度因子的公式，最后计算若干典型例子裂纹前沿的应力强度因子。     

三维断裂力学的超奇异积分方程方法

本文利用有限部积分的概念和方法,严格地证明了三维弹性体中受任意载荷作用的平片裂纹问题的超奇异积分方程组,并对未知解的性态作了理论分析,得到了性态指数,在此基础上通过主部分析,精确地求得了裂纹前沿光滑点附近的奇性应力场,从而找到了以裂纹面位移间断(位错)表示的应力强度因子表达式,最后对所得的超奇异积分方程组建立了数值法,并用此计算了若干典型的平片裂纹问题,数值结果令人满意。     

椭圆平片裂纹前沿应力强度因子解析解的超奇异积分方程方法

对无限大三维均质弹性体中任意平片裂纹的超奇异积分方程,巧妙地引入椭球坐标系和利用裂纹表面位移间断人有平方根的特性,获得了受任意方向均布压力作用下椭圆平片裂纹问题的超奇异积分方程的解析解。运用这些解析解和应力强度因子的定义,得到了裂纹前沿Ⅰ型、Ⅱ型和Ⅲ型和应力强度因子的精确表达式,所得结果与现有精确解完全一致。     

三维体内含垂直于双材料界面混合型裂纹的超奇异积分解法

利用双材料位移基本解和Somigliana公式,将三维体内含垂直于双材料界面混合型裂纹问题归结为求解一组超奇异积分方程。使用主部分析法,通过对裂纹前沿应力奇性的分析,得到用裂纹面位移间断表示的应力强度因子的计算公式,进而利用超奇异积分方程未知解的理论分析结果和有限部积分理论,给出了超奇异积分方程的数值求解方法。最后,对典型算例的应力强度因子做了计算,并讨论了应力强度因子数值结果的收敛性及其随各参数变化的规律。     

磁电热弹耦合材料三维多裂纹超奇异积分法

用超奇异积分方程法将多场耦合载荷作用下磁电热弹耦合材料内含任意形状和位置三维多裂纹问题转化为求解一以广义位移间断为未知函数的超奇异积分方程组问题,退化得到内含任意形状平行三维多裂纹问题的超奇异积分方程组;推导出平行三维多裂纹问题的裂纹前沿广义奇异应力场解析表达式、定义了广义(应力、应变能)强度因子和广义能量释放率;应用有限部积分概念及体积力法,为超奇异积分方程组建立了数值求解方法,编制了FORTRAN程序,以平行双裂纹为例,通过典型算例,研究了广义(应力、应变能)强度因子随裂纹位置、裂纹形状及材料参数变化规律,得到裂纹断裂评定准则. 最后,分析了裂纹间干扰、屏蔽作用及其在工程实际中的应用.      

横观各向同性压电材料中裂纹问题的边界元法

采用Somigiliana公式给出了三维横观各向同性压电材料中的非渗漏裂纹问题的一般解和超奇异积分方程，其中未知函数为裂纹面上的位移间断和电势间断．在此基础上，使用有限部积分和边界元结合的方法，建立了超奇异积分方程的数值求解方法，并给出了一些典型数值算例的应力强度因子和电位移强度因子的数值结果，结果令人满意．     

APPLICATION OF THE HYPERSINGULAR INTEGRAL EQUATIONS METHOD TO THE INTERACTION BETWEEN TWO PARALLEL PLANAR CRACKS IN 3-D ELASTICITY

By using the analytic theory of hypersingular integral equations in three-dimensional fracture mechanics, the interactions between two parallel planar cracksunder arbitrary loads are investigated. According to the concepts and method of finite-part integrals, a set of hypersingular integral equations is derived, in which theunknown functions are the displacement discontinuities of the crack surfaces. Then itsnumerical method is proposed by combining the finite-part integral method with theboundary element method. Based on the above results, the method for calculating thestress intensity factors with the displacement discontinuities of the crack surfaces ispresented. Finally, several typical examples are calculated and the numerical resultsare satisfactory.     

Finite-part integral and boundary element method to solve flat crack problems

Using the Somigliana formula and the concepts of finite-part integral, a set of hypersingular integral equations to solve the arbitrary fiat crack in three-dimensional elasticity is derived and its numerical method is then proposed by combining the finitepart integral method with boundary element method. In order to verify the method, several numerical examples are carried out. The results of the displacement discontinuities of the crack surface and the stress intensity factors at the crack front are in good agrernent with the' theoretical solutions.     

Finite-part integral and boundary element method for interaction between two parallel planar cracks in three-dimensional finite bodies

By using the Somigliana representation and the concepts of finite-part integrals, a set of hypersingular integral equations of the interaction between two parallel planar cracks in a three-dimensional finite body subjected to arbitrary loads is derived, and then its numerical method is proposed by the finite-part integral method combined with the boundary element method. According to the analytic theory of hypersingular integral equations, the square root models of displacement discontinuities in the elements near the crack front are applied, and thus the computational precision is raised. Based on this, the stress intensity factors can be directly calculated. Finally, the stress intensity factors of several typical interaction problems are calculated.     

Finite-part integral and boundary element method to solve three-dimensional crack problems in piezoelectric materials

Using the hypersingular integral equation method based on body force method, a planar crack in a three-dimensional transversely isotropic piezoelectric solid under mechanical and electrical loads is analyzed. This crack problem is reduced to solve a set of hypersingular integral equations. Compare with the crack problems in elastic isotropic materials, it is shown that for the impermeable crack, the intensity factors for piezoelectric materials can be obtained from those for elastic isotropic materials. Based on the exact analytical solution of the singular stresses and electrical displacements near the crack front, the numerical method of the hypersingular integral equation is proposed by the finite-part integral method and boundary element method, which the square root models of the displacement and electric potential discontinuities in elements near the crack front are applied. Finally, the numerical solutions of the stress and electric field intensity factors of some examples are given.     

Analysis of multiple interfacial cracks in three-dimensional bimaterials using hypersingular integro-differential equation method

By using the concept of finite-part integral, a set of hypersingular integro-differential equations for multiple interracial cracks in a three-dimensional infinite bimaterial subjected to arbitrary loads is derived. In the numerical analysis, unknown displacement discontinuities are approximated with the products of the fundamental density functions and power series. The fundamental functions are chosen to express a two-dimensional interface crack rigorously. As illustrative examples, the stress intensity factors for two rectangular interface cracks are calculated for various spacing, crack shape and elastic constants. It is shown that the stress intensity factors decrease with the crack spacing.     

横观各向同性材料的三维断裂力学问题

从三维横观各向同性材料弹性力学理论出发, 使用Hadamard有限部积分概念, 导出了三维状态下单位位移间断(位错)集度的基 本解. 在此基础上, 进一步运用极限理论, 将任意载荷作用下, 三维无限大横观各向 同性材料弹性体中, 含有一个位于弹性对称面内的任意形状的片状裂纹问题, 归结为求 解一组超奇异积分方程的问题. 通过二维超奇异积分的主部分析方法, 精确地求得了裂纹前沿光滑点附近的应力奇异指数和奇异应力场, 从而找到了以裂纹表面位移间断表示的应力强度因子表达式及裂纹局部扩展所提供 的能量释放率. 作为以上理论的实际应用，最后给出了一个圆形片状裂纹问题 的精确解例和一个正方形片状裂纹问题的数值解例. 对受轴对称法向均布载荷作用下圆形片状裂纹问题, 讨论了超奇异积分方程的精确求解方法, 并获得了位移间断和应力强度因子的封闭解, 此结果与现有理论解完全一致.     

Application of hypersingular integral equation method to three-dimensional crack in electromagnetothermoelastic multiphase composites

A three-dimensional crack problem in electromagnetothermoelastic multiphase composites (EMTE-MCs) under extended loads is investigated in this paper. Using Green’s functions, the extended general displacement solutions are obtained by the boundary element method. This crack problem is reduced to solving a set of hypersingular integral equations coupled with boundary integral equations, in which the unknown functions are the extended displacement discontinuities. Then, the behavior of the extended displacement discontinuities around the crack front terminating at the interface is analyzed by the main-part analysis method of hypersingular integral equations. Analytical solutions for the extended singular stresses, the extended stress intensity factors (SIFs) and the extended energy release rate near the crack front in EMTE-MCs are provided. Also, a numerical method of the hypersingular integral equations for a rectangular crack subjected to extended loads is put forward with the extended displacement discontinuities approximated by the product of basic density functions and polynomials. In addition, distributions of extended SIFs varying with the shape of the crack are presented. The results show that the present method accurately yields smooth variations of extended SIFs along the crack front.