考虑尺寸效应的双材料轴对称界面端应力奇异性

提出了一种分析双材料轴对称界面端的应力奇异行为的特征值法.基于弹性力学空间轴对称问题的基本方程和一阶近似假设,利用分离变量形式的位移函数和无网格算法,导出了关于应力奇异性指数的离散形式的奇异性特征方程.由奇异性特征方程的特征值和特征向量,即可确定应力奇异性指数、位移角函数和应力角函数.数值求解了纤维/基体轴对称界面端模型的奇异性特征方程, 结果表明:尺寸效应参数δ(奇异点与轴对称轴的距离和应力奇异性支配区域大小的比值)影响着应力奇异性的强弱与阶次, 准一阶近似解析解只是δ>>1时的一个特例.      

双材料反平面问题界面端奇异应力场分析

利用位移函数的级数展开,对任意角度的反平面问题界面端的应力场进行了分析研究,得到了全场解。研究一阶场后发现,奇异规律与一般平面问题界面端有显著区别,在界面端关于界面对称的情况下,平角界面端(θ1 = θ2 = θ = 90°) 应力场没有奇异性,其它形状的界面端随着角度θ 从90°到180°,奇异指数也从0到0.5。当界面端是非对称时,平角界面端(θ1 θ2 = 180°)、直角界面端(θ1 = 90°,θ2 = 180°)以及其它形状界面端的奇异指数是一个与两相材料常数比Γ有关的常数。以上两种情况下的应力强度因子完全类似单相材料中裂纹尖端附近应力强度因子,故可根据定义得到     

穿过两相材料界面的曲线型刚性线

利用两相材料中集中力的基本解，建立了求解曲线型刚性线夹杂和两相材料界面相交问题的弱奇异积分方程。通过Cauchy型奇异积分方程主部分析方法，得出穿过两相材料界面的曲线型刚线性在交点处的奇性应力指数及交点处角形域内的奇性应力，并利用奇性应力定义了交点处的应力奇异因子。通过对弱奇异积分方程的数值求解，得出了刚性线端点和交点处的应力奇异因子。     

多参量断裂力学广义杂交/混合裂纹有限元法

通过研究广为人知的断裂力学单变量八节点位移裂纹QPE元和Akin族奇异单元法，本文运用经典局部裂纹解析解，与非协调假设应力杂交-混合元列式方法相结合，提出用于分层各向异性材料的多变量半解析假设应力奇异广义杂交／混合裂纹有限元法，能克服现有位移裂纹元法的域应力分布精度低和高次单元所需计算容量大的局限性，互为补充，更有利于结构裂纹扩展分析和应用研究。文中设计了一个半解析奇异裂纹平面单元，各向同性材料板算例验证了退化二次八节点协调位移裂纹元及六节点非协调奇异应力裂纹元，说明采用稀疏及加密单元网格，两类裂纹单元分别从上下逼近收敛于实验和理论参考解，可得到吻合程度较好的1／√r奇异应变和应力分量以及应力强度因子值，表明了本文奇异裂纹单元理论的优越性。     

动态断裂力学的无限相似边界元法

对弹性动力学的相似边界元法进行了进一步研究,推导了相应的计算公式,并在此基础上提出了动态断裂力学的无限相似边界元法.与传统的边界元法相比,相似边界元法由于只需在少数单元上进行数值积分,大大减少了计算量.对动态断裂力学问题,无限相似边界元法由于在裂纹尖端的边界上设置了逼近于裂纹尖端的无限个相似边界单元,可直接得到裂纹尖端具有奇异性的应力,而不需要设置奇异单元,从而突破了奇异单元对应力奇异性阶次的局限.另外,还讨论了无限相似边界元法得到的无限阶的线性代数方程组的求解方法.     

各向同性与各向异性三相材料接头应力奇异性研究

研究了各向同性与各向异性三相材料接头的应力奇性指数,通过引入奇异点附近区域位移场渐近展开的典型项,将各向同性与各向异性组合材料接头的控制方程和径向边界条件转化为变系数常微分方程的特征值问题;再利用插值矩阵法求解所建立的特征方程,得到接头端部的应力奇性指数和特征角函数。对由两个各向异性材料和一个各向同性材料以任意楔形角组成的三相接头结构的奇异性进行了研究,并比较了它们的应力奇性指数。计算结果表明:对于粘结接头,各向同性材料刚度越大应力奇异性越强;对于剥离接头,各向同性材料楔形角或材料刚度越大,第一阶应力奇异性越弱。计算结果与已有文献的结果对比吻合良好,证明了本文方法的有效性。     

不同硬化指数的幂次硬化结合材料的界面端奇异性分析

从位移匹配的观点出发,本文认为对任意结合角度的幂次硬化材料界面端弹塑性问题,如果两种材料的硬化指数不相同,则应力场的奇异次数应由硬化指数较高一方材料的材料性质和几何形状决定。进一步分析表明,奇异次数只与该材料的硬化指数n及界面端角度有关,与比例常数α等其它材料常数无关。通过边界元数值计算对上述结论进行了验证,并且发现随着硬化指数的提高,应力奇异次数降低。     

集中力拉伸楔体大变形理论分析及数值计算

对任意向集中力拉伸用下的橡胶类材料楔体尖端应力应变场进行了渐近分析,并利用有限元与无限元耦合作了数值计算,结果表明大应变情况下,楔体尖端场处于单向拉伸状态,应力主奇异项为ｒ＾－λ,ｒ为变形后至顶点距离,λ依赖于材料本构常数,数值计算与理论解吻合很好。     

三维基体裂纹与界面垂直接触的 线奇异分析1)

使用三维断裂力学的超奇异积分方程方法,对两相材料中出现的三维基体裂纹与材料分界面垂直接触的应力线奇异问题作了理论分析,求得了线奇异界面应力的表达式.     

三维基体裂纹与界面垂直接触的线奇异分析

使用三维断裂力学的超奇异积分方程方法,对两相材料中出现的三维基体裂纹与材料分界面垂直接触的应力线奇异问题作了理论分析,求得了线奇异界面应力的表达式     

On the intensity of linear elastic high order singularities ahead of cracks and re-entrant corners

The paper deals with high order elastic singular terms at cracks and re-entrant corners (sharp V-notches), which are commonly omitted in linear elastic analyses by the argument that the strain energy and displacements in the near-tip region should be bounded. The present analysis proves that these terms are fully included in the elastic part of complete elastic–plastic stress and strain solutions.The intensities of high order singular terms are found to be linked to the linear elastic stress intensity factor and the extension of the plastic zone along the crack bisector line. The smaller the plastic radius, the smaller the intensities of high order singular terms are.A physical justification of the existence of high order singular terms is provided on the basis of the strain energy density distribution detected along the crack bisector line. Finally, the influence of the V-notch opening angle is made explicit, discussing also the relationship between the singularity orders and the solution of a Williams’ type sinusoidal eigen-equation.     

A new boundary element approach of modeling singular stress fields of plane V-notch problems

In this paper, a new boundary element (BE) approach is proposed to determine the singular stress field in plane V-notch structures. The method is based on an asymptotic expansion of the stresses in a small region around a notch tip and application of the conventional BE in the remaining region of the structure. The evaluation of stress singularities at a notch tip is transformed into an eigenvalue problem of ordinary differential equations that is solved by the interpolating matrix method in order to obtain singularity orders (degrees) and associated eigen-functions of the V-notch. The combination of the eigen-analysis for the small region and the conventional BE analysis for the remaining part of the structure results in both the singular stress field near the notch tip and the notch stress intensity factors (SIFs).Examples are given for V-notch plates made of isotropic materials. Comparisons and parametric studies on stresses and notch SIFs are carried out for various V-notch plates. The studies show that the new approach is accurate and effective in simulating singular stress fields in V-notch/crack structures.     

Singular electro-mechanical fields near the apex of a piezoelectric bonded wedge under antiplane shear

This paper presents the explicit forms of singular electro-mechanical field in a piezoelectric bonded wedge subjected to antiplane shear loads. Based on the complex potential function associated with eigenfunction expansion method, the eigenvalue equations are derived analytically. Contrary to the anisotropic elastic bonded wedge, the results of this problem show that the singularity orders are single-root and may be complex. The stress intensity factors of electrical and mechanical fields are dependent. However, when the wedge angles are equal (α=β), the orders become real and double-root. The real stress intensity factors of electrical and mechanical fields are then independent. The angular functions have been validated when they are compared with the results of several degenerated cases in open literatures.     

Crack tip asymptotic character of anti-plane stress and strain rate for linear fractional constitutive relations

Linear-fractional strain rate and stress relations are used to simulate materials undergoing steady state creep. The crack tip asymptotic character of the stress and strain rate field is obtained in exact and approximate form. In the limit as the radial distance emanating from the crack tip approaches zero, the stress field corresponds to that for an ideal plastic material while the exact and approximate solutions tend to coincide. Discussed is the nonhomogeneous singular character of the strain rate field that possess different orders of singularities in a circular region around the crack tip.     

Stress singularities resulting from various boundary conditions in angular corners of plates of arbitrary thickness in extension

The stress singularities in angular corners of plates of arbitrary thickness with various boundary conditions subjected to in-plane loading are studied within the first-order plate theory. By adapting an eigenfunction expansion approach a set of characteristic equations for determining the structure and orders of singularities of the stress resultants in the vicinity of the vertex is developed. The characteristic equations derived in this paper incorporate that obtained within the classical plane theory of elasticity (M.L. Williams’ solution) and also describe the possible singular behaviour of the out-of-plane shear stress resultants induced by various boundary conditions.     

V形切口应力强度因子的一种边界元分析方法

将V形切口结构分成围绕切口尖端的小扇形和剩余结构两部分. 尖端处扇形域应力场表示成关于尖端距离$\rho$的渐近级数展开式,从线弹性理论方程推导出了一组分析平面V形切口奇异性的常微分方程特征值问题,通过求解特征方程,得到前若干个奇性指数和相应的特征向量. 再将切口尖端的位移和应力表示为有限个奇性阶和特征向量的组合. 然后用边界元法分析挖去小扇形后的剩余结构. 将位移和应力的线性组合与边界积分方程联立,求解获得切口根部区域的应力场、应力幅值系数和整体结构的位移和应力. 从而准确计算出平面V形切口的奇异应力场和应力强度因子.      

A unified definition for stress intensity factors of interface corners and cracks

Based upon linear fracture mechanics, it is well known that the singular order of stresses near the crack tip in homogeneous materials is a constant value −1/2, which is nothing to do with the material properties. For the interface cracks between two dissimilar materials, the near tip stresses are oscillatory due to the order of singularity being −1/2 ± iε and −1/2. The oscillation index ε is a constant related to the elastic properties of both materials. While for the general interface corners, their singular orders depend on the corner angle as well as the elastic properties of the materials. Owing to the difference of the singular orders of homogeneous cracks, interface cracks and interface corners, their associated stress intensity factors are usually defined separately and even not compatibly. Since homogenous cracks and interface cracks are just special cases of interface corners, in order to build a direct connection among them a unified definition for their stress intensity factors is proposed in this paper. Based upon the analytical solutions obtained previously for the multibonded anisotropic wedges, the near tip solutions for the general interface corners have been divided into five different categories depending on whether the singular order is distinct or repeated, real or complex. To provide a stable and efficient computing approach for the general mixed-mode stress intensity factors, the path-independent H-integral based on reciprocal theorem of Betti and Rayleigh is established in this paper. The complementary solutions needed for calculation of H-integral are also provided in this paper. To illustrate our results, several different kinds of examples are shown such as cracks in homogenous isotropic or anisotropic materials, central or edge notches in isotropic materials, interface cracks and interface corners between two dissimilar materials.     

穿过反平面圆形夹杂界面的曲线型刚性线

利用复变函数方法和叠加原理建立了求解刚性线夹杂问题的弱奇积分方程,利用Ｃａｕｃｈｙ型奇异积分方程主部分方法,研究了穿过反平面圆夹杂界面的曲线型刚性线夹杂在界面交点处点处的奇性应力指数以及交点处角形域内的奇性应力,并定义了交点处的应力奇性因子。利用所得的奇性应力指数,通过对弱奇异积分方程的数值求解,得出了刚性线端点和交点处的应力奇性因子。     

三维切口/裂纹结构的扩展边界元法分析

在线弹性理论中,三维 V 形切口/裂纹结构尖端区域存在多重应力奇异性,常规数值方法不易求解. 本文提出和建立了三维扩展边界元法 (XBEM),用于分析三维线弹性 V 形切口/裂纹结构完整的位移和应力场. 先将三维线弹性 V 形切口/裂纹结构分为尖端小扇形柱和挖去小扇形柱后的外围结构. 尖端小扇形柱内的位移函数采用自尖端径向距离 $r$ 的渐近级数展开式表达,其中尖端区域的应力奇异指数、位移和应力特征角函数通过插值矩阵法获得. 而级数展开式各项的幅值系数作为基本未知量. 挖去扇形域后的外围结构采用常规边界元法分析. 两者方程联立求解可获得三维 V 形切口/裂纹结构完整的位移和应力场,包括切口/裂纹尖端区域精细的应力场. 扩展边界元法具有半解析法特征,适用于一般三维 V 形切口/裂纹结构完整位移场和应力场的分析,其解可精细描述从尖端区域到整体结构区域的完整应力场. 作者研制了三维扩展边界元法程序,文中给出了两个算例,通过计算结果分析,表明了扩展边界元法求解三维 V 形切口/裂纹结构完整应力场的准确性和有效性.      

ANALYSIS OF 3-D NOTCHED/CRACKED STRUCTURES BY USING EXTENDED BOUNDARY ELEMENT METHOD 1)

在线弹性理论中,三维 V 形切口/裂纹结构尖端区域存在多重应力奇异性,常规数值方法不易求解. 本文提出和建立了三维扩展边界元法 (XBEM),用于分析三维线弹性 V 形切口/裂纹结构完整的位移和应力场. 先将三维线弹性 V 形切口/裂纹结构分为尖端小扇形柱和挖去小扇形柱后的外围结构. 尖端小扇形柱内的位移函数采用自尖端径向距离 $r$ 的渐近级数展开式表达,其中尖端区域的应力奇异指数、位移和应力特征角函数通过插值矩阵法获得. 而级数展开式各项的幅值系数作为基本未知量. 挖去扇形域后的外围结构采用常规边界元法分析. 两者方程联立求解可获得三维 V 形切口/裂纹结构完整的位移和应力场,包括切口/裂纹尖端区域精细的应力场. 扩展边界元法具有半解析法特征,适用于一般三维 V 形切口/裂纹结构完整位移场和应力场的分析,其解可精细描述从尖端区域到整体结构区域的完整应力场. 作者研制了三维扩展边界元法程序,文中给出了两个算例,通过计算结果分析,表明了扩展边界元法求解三维 V 形切口/裂纹结构完整应力场的准确性和有效性.