三维界面裂纹的理论分析

使用界面裂纹的基本解及有限部积分的方法 ,将三维界面裂纹问题 ,归为解一组以裂纹面位移间断为未知函数的超奇异积分微分方程 ,然后为此组方程的求解作了系统的理论分析     

横观各向同性电磁弹性介质中裂纹对SH波的散射

研究横观各向同性电磁弹性介质中裂纹和反平面剪切波之间的相互作用．根据电磁弹性介质的平衡运动微分方程、电位移和磁感应强度微分方程,得到SH波传播的控制场方程．引入线性变换,将控制场方程简化为Helmholtz方程和两个Laplace方程A·D2通过Fourier变换,并采用非电磁渗透型裂面边界条件,得到了柯西奇异积分方程组．利用Chebyshev多项式求解积分方程,得到应力场、电场和磁场以及动应力强度因子的表达,并给出了数值算例．     

各向异性介质中周期性界面裂纹的弹塑性问题

应用富里叶积分变换方法将裂纹边值问题化为对偶积分方程组,再用定积分变换法将问题进一步化为奇异积分方程组,求得了双材料各向异性弹塑性介质中周期性界面裂纹反平面问题的封闭形式解,并作为特例讨论了各向同性双材料问题、各向异性单一材料问题及各向同性—各向异性双材料问题.结果表明:裂纹尖端前沿的塑性区尺寸、裂纹的张开位移(COD)均决定于两种材料流动极限中的较小者及裂纹的长度和相邻两裂纹的间距,此外,COD还与材料模量有关.     

求解平片裂纹问题的有限部积分与边界元法

本文利用位移的Somigliana公式和有限部积分的概念,导出了求解三维弹性力学中的任意形状平片裂纹问题的超奇异积分方程组,进而联合使用有限部积分法与边界元法对所得方程建立了数值法.为验证本文的方法,计算了若干数值例子的裂纹面的位移间断及裂纹前沿的应力强度因子,它们与理论值相比符合很好.     

横观各向同性电磁弹性固体耦合方程的一般解

横观各向同性电磁弹性固体的耦合特征由5个关于弹性位移、电位和磁位的二阶偏微分方程控制.基于势函数理论,耦合的方程组被简化为5个非耦合的关于势函数的广义Laplace方程.弹性场和电磁场由势函数表示,这构成了横观各向同性电磁弹性固体的一般解.     

层状横观各向同性饱和地基上圆板的非轴对称振动

研究了层状横观各向同性饱和地基上弹性圆板的非轴对称振动问题．首先,通过方位角的Fourier变换,将圆柱坐标系下横观各向同性饱和土的三维动力方程转化为一阶常微分方程组,基于径向Hankel变换,建立问题的状态方程,求解状态方程后得到传递矩阵;其次,利用传递矩阵,结合层状饱和地基的边界条件、排水条件及层间接触和连续条件,给出了任意简谐激振力作用下层状横观各向同性饱和地基动力响应的通解;然后,按混合边值问题建立层状饱和地基上弹性圆板非轴对称振动的对偶积分方程,并将对偶积分方程化为易于数值计算的第二类Fredholm积分方程,并给出了算例．     

带一类非Carleman正位移的拟线性奇异积分方程

空气动力学中出现的某些混合型偏微分方程组的边值问题往往归结为带非Carleman位移的奇异积分方程。关于带Carleman位移的线性奇异积分方程已有系统理论。可见专著和文献等。而关于带非Carleman位移的奇异积分方程的理论目前还很不完整。专著介绍了一些结果。 本文研究带位移的拟线性奇异积分方程     

三维无限接合体双材料多界面裂纹的超奇异微积分方程法

利用有限部积分的概念,导出了三维无限接合体中多个界面裂纹,在任意载荷作用下的超奇异微积分方程组.数值分析中,未知的位移间断采用基本分布函数和多项式乘积的形式来近似,其中基本分布函数是根据界而裂纹应力的振荡奇异性来选取的.作为典型算例,研究了存在两个矩形界面裂纹时,裂纹之间距离、裂纹形状及双材料弹性常数对应力强度因子的影响.计算表明,应力强度因子随裂纹间的距离的增大而减小.     

横观各向同性饱和地基的三维动力响应

首先引入位移函数,将直角坐标系下横观各向同性饱和土Biot波动方程转化为2个解耦的六阶和二阶控制方程;然后基于双重Fourier变换,求解了Biot波动方程,得到以土骨架位移和孔隙水压力为基本未知量的积分形式的一般解,并用一般解给出了饱和土总应力分量的表达式．在此基础上系统研究了横观各向同性饱和半空间体的稳态动力响应问题,考虑表面排水和不排水两种情况,得到了半空间体在任意分布的表面谐振荷载作用下,表面位移的稳态动力响应,文末给出了算例．     

横观各向同性饱和地基与中厚圆板的非轴对称动力相互作用

研究横观各向同性饱和土地基上中厚弹性圆板的非轴对称振动问题,即首先利用Fourier展开和Hankel变换技术,求解了简谐激励下横观各向同性饱和土地基的非轴对称Biot波动方程,然后按混合边值问题建立地基与弹性中厚圆板非轴对称动力相互作用的对偶积分方程,并将对偶积分方程化为易于数值计算的第二类Fredholm积分方程．文末给出了算例．数值结果表明,在一定频率范围内,地基表面的位移幅值随激振频率增加而增大,随距离的增大以振荡形式衰减变化．     

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.     

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.     

3D modeling of crack growth in electro-magneto-thermo-elastic coupled viscoplastic multiphase composites

This work presents a time-domain hypersingular integral equation (TD-HIE) method for modeling 3D crack growth in electro-magneto-thermo-elastic coupled viscoplastic multiphase composites (EMTE-CVP-MCs) under extended incremental loads rate through intricate theoretical analysis and numerical simulations. Using Green’s functions, the extended general incremental displacement rate solutions are obtained by time-domain boundary element method. Three-dimensional arbitrary crack growth problem in EMTE-CVP-MCs is reduced to solving a set of TD-HIEs coupled with boundary integral equations, in which the unknown functions are the extended incremental displacement discontinuities gradient. Then, the behavior of the extended incremental displacement discontinuities gradient around the crack front terminating at the interface is analyzed by the time-domain main-part analysis method of TD-HIE. Also, analytical solutions of the extended singular incremental stresses gradient and extended incremental integral near the crack fronts in EMTE-CVP-MCs are provided. In addition, a numerical method of the TD-HIE for a 3D crack subjected to extended incremental loads rate is put forward with the extended incremental displacement discontinuities gradient approximated by the product of time-domain basic density functions and polynomials. Finally, examples are presented to demonstrate the application of the proposed method.     

Analysis of anti-plane interface cracks in one-dimensional hexagonal quasicrystal coating

The displacement discontinuity method is extended to study the fracture behavior of interface cracks in one-dimensional hexagonal quasicrystal coating subjected to anti-plane loading. The Fredholm integral equation of the first kind is established in terms of displacement discontinuities. The fundamental solution for anti-plane displacement discontinuity is derived by the Fourier transform method. The singularity of stress near the crack front is analyzed, and Chebyshev polynomials of the second kind are numerically adopted to solve the integral equations. The displacement discontinuities across crack faces, the stress intensity factors, and the energy release rate are calculated from the coefficients of Chebyshev polynomials. In combination with numerical simulations, a comprehensive study of influencing factors on the fracture behavior is conducted.     

A new boundary element method for the solution of plane steady-state thermoelastic fracture mechanics problems

An efficient numerical technique is developed for plane, homogeneous, isotropic, steady-state thermoelasticity problems involving arbitrary internal smooth and/or kinkedcracks. The thermal stress intensity factors and relative crack surface displacements due to steady-state temperature distributions are determined and compared to available solutions obtained by other methods. In these analyses the thermal boundary conditions across the crack surface are assumed to be insulated. The present approach involves coupling the direct boundary integral equations to newly developed crack integral equations.     

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.     

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

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