含裂纹金属薄板胶接补强的应力分析

本文应用弹塑性有限元分析胶接补强的含裂纹薄板结构,剪切单元被改进使之可连结平面等参单元并用于分析胶层应力。计算了补强后裂纹板的应力强度因子,并将计算结果与实验数据作了比较,分析表明,胶接补强可显著降低裂尖的应力集中,使裂纹板内应力分布趋于均匀。     

含裂纹粘接结构的边界元应力分析

本文采用子域法和直接应力奇异单元法求解二维粘接结构中的裂纹问题。子域法把粘接结构划分为三个子域,根据每个子域的边界积分方程和子域间的界面条件,可以建立粘接结构的边界积分方程组。直接应力奇异单元能够在整个单元长度上反映裂纹端部的1/r~(1/2)奇异性,在计算时可以通过坐标变换消除奇异单元积分中的奇异性,直接计算出应力强度因子。含裂纹多层结构的数值示例和粘接补强单边裂纹板的应力测试和疲劳试验结果证实了本文方法的有效性。     

胶补加筋板应力强度因子及胶层剪应力计算

 在铆钉力法的基础上，提出计算胶接修补加筋板结构应力强度因子及胶层剪应力的解析法, 对分析胶接修补的加筋板强度具有现实意义. 利用影响系数，可以得到筋条、补片及蒙皮位 移的表达式，根据位移协调条件，建立联立方程，从而得到胶层剪应力及铆钉力. 通过解析 式，计算裂纹尖端的应力强度因子，分析补片对结构修补的效果，综合考虑了筋条及补片对 应力强度因子的影响.     

胶补加筋板应力强度因子及胶层剪应力计算

在铆钉力法的基础上，提出计算胶接修补加筋板结构应力强度因子及胶层剪应力的解析法，对分析胶接修补的加筋板强度具有现实意义．利用影响系数，可以得到筋条、补片及蒙皮位移的表达式，根据位移协调条件，建立联立方程，从而得到胶层剪应力及铆钉力．通过解析式，计算裂纹尖端的应力强度因子，分析补片对结构修补的效果，综合考虑了筋条及补片对应力强度因子的影响．     

含多源损伤的胶接修补结构有限元分析方法

采用有限元法研究含多源损伤结构的胶接修补问题,利用二维三层有限元模型对损伤区进行了数值模拟,并选取典型多源损伤情况中含共线双裂纹铝板结构为算例,详细分析了含多裂纹胶接修补结构中两裂纹相对位置、补片尺寸、铺层和厚度对应力强度因子的影响.结果表明,复合材料胶接修补可明显降低含共线双裂纹母板的应力强度因子;对于确定的裂纹和应力场,应对复合材料补片长度和厚度等参数进行优化设计,以获得最佳的修补效果.     

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

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

基于有限元方法对人工径向裂纹薄圆板振动辐射声场的研究

在有限元模态分析的基础上，提取了含人工径向裂纹的周界固定薄圆板各单元的相关参数，并将瑞利积分离散化后进而计算了裂纹圆板辐射声场的轴向声压分布和r=0．5m球面上声场指向性。结果表明，径向裂纹不仅使模态裂解为关于裂纹的对称模态和反对称模态，而且使声场分布有显著的变化。通过完整圆板辐射声场计算验证了本文提出方法的正确性。可用于求解任意形状穿透裂纹薄圆板辐射声场的计算。     

含界面裂纹Reissner板弯曲问题分析的奇异单元

首先,采用特征函数渐近展开法,推导了Reissner板弯曲界面裂纹尖端附近位移场渐近展开的前两阶显式表达式,并利用所获得的位移场渐近表达式构造了一种可用于Reissner板弯曲界面裂纹分析的奇异单元。然后,将该奇异单元与外部的常规有限单元相结合,开展了含界面裂纹Reissner板弯曲断裂问题的数值分析。奇异单元可以较好地描述裂纹尖端附近的内力场与位移场,其优势是它与常规单元进行连接时不需要使用过渡单元,并且可以直接给出应力强度因子等断裂参数的高精度数值结果。最后,通过两个数值算例验证了本文方法的有效性。     

FRP加固金属裂纹板的断裂力学分析

传统的金属结构加固方法会形成新的疲劳源,而粘贴FRP加固则具有明显的优势．提出了“三维实体-弹簧-壳元”有限元模型,金属板采用三维实体单元, FRP采用壳单元,用弹簧单元来模拟FRP与金属板之间的胶层,对金属裂纹板粘贴FRP加固后的性能进行了线弹性断裂力学分析,并对影响金属板裂纹前缘应力强度因子的参数进行了讨论．分析结果表明,采用高弹性模量的FRP和增加FRP的厚度对改善加固效果非常明显．     

自然单元法在线弹性断裂力学中的应用研究

采用Laplace插值构造近似场函数,代入弹性力学边值问题的变分弱形式,由变分原理推导出自然单元法的离散控制方程。运用线弹性断裂力学理论,分析了有限板单边Ⅰ型裂纹的应力强度因子和有限板拉剪复合裂纹的扩展,给出了分析线弹性断裂力学问题的自然单元法。数值计算结果表明了方法的正确性和有效性。     

单搭接胶接接头拉伸剪切性能的数值模拟与实验研究

采用数值模拟和光测技术对单向拉伸载荷作用下单搭接胶接接头中的剪切性能进行分析,研究了不同厚度胶层中切应力的变化规律。用有限元方法(FEM)对不同胶层厚度的试件进行建模,得到了拉伸载荷下胶粘剂中的切应力分布及其统计参数。利用数字图像相关(digitalimage correlation,DIC)方法对试件的变形场进行测量。结果表明,当胶粘剂的厚度较小时,胶粘剂中的切应力的分布统计参数随着其厚度的增加会有显著的变化,但是当厚度超过一定的数值时,统计参数对厚度的变化不再敏感。     

An effective boundary element method for analysis of crack problems in a plane elastic plate

A simple and effective boundary element method for stress intensity factor calculation for crack problems in a plane elastic plate is presented. The boundary element method consists of the constant displacement discontinuity element presented by Crouch and Starfield and the crack-tip displacement discontinuity elements proposed by YAN Xiangqiao. In the boundary element implementation the left or the right crack-tip displacement discontinuity element was placed locally at the corresponding left or right each crack tip on top of the constant displacement discontinuity elements that cover the entire crack surface and the other boundaries. Test examples (i. e. , a center crack in an infinite plate under tension, a circular hole and a crack in an infinite plate under tension) are included to illustrate that the numerical approach is very simple and accurate for stress intensity factor calculation of plane elasticity crack problems. In addition, specifically, the stress intensity factors of branching cracks emanating from a square hole in a rectangular plate under biaxial loads were analysed. These numerical results indicate the present numerical approach is very effective for calculating stress intensity factors of complex cracks in a 2-D finite body, and are used to reveal the effect of the biaxial loads and the cracked body geometry on stress intensity factors.     

双轴载荷作用下源于椭圆孔的分支裂纹的一种边界元分析

利用一种边界元方法来研究双轴载荷作用下无限大板中源于椭圆孔的分支裂纹．该边界元方法由Crouch与Starfied建立的常位移不连续单元和笔者提出的裂尖位移不连续单元构成．在该边界元方法的实施过程中,左、右裂尖位移不连续单元分别置于裂纹的左、右裂尖处,而常位移不连续单元则分布于除了裂尖位移不连续单元占据的位置之外的整个裂纹面及其它边界,文中算例说明本数值方法对计算平面弹性裂纹的应力强度因子是非常有效的。该文对双轴载荷作用下无限大板中源于椭圆孔的分支裂纹的数值结果进一步证实本数值方法对计算复杂裂纹的应力强度因子的有效性,同时该数值结果可以揭示双轴载荷及裂纹体几何对应力强度因子的影响。     

内部压力作用下矩形板中源于椭圆孔的分支裂纹应力强度因子的一种数值分析

应用一种边界元方法来研究内部压力作用下矩形板中源于椭圆孔的分支裂纹。该边界元方法由Crouch与Starfied建立的常位移不连续单元和笔者最近提出的裂尖位移不连续单元构成。在该边界元方法的实施过程中,左、右裂尖位移不连续单元分别置于裂纹的左、右裂尖处,而常位移不连续单元则分布于除了裂尖位移不连续单元占据的位置之外的整个裂纹面及其它边界。本数值结果进一步证实这种数值方法对计算有限大板中复杂裂纹的应力强度因子的有效性,同时该数值结果可以揭示裂纹体几何对应力强度因子的影响。     

The effect of load biaxiality on crack growth in non-linear materials

The influence of load biaxiality on the stress field and fracture behavior of a cracked plate is investigated. Considered is a square plate containing a central through the thickness crack and subjected to a biaxial loading perpendicular and parallel to the crack plane. The stress field of the plate is analyzed by a finite element code based on incremental plasticity and the von Mises yield condition. A method based on the strain energy density theory is used to determine the critical stress for crack initiation. It was found that the equi-biaxial loading mode induces the smallest plastic zones, while the critical applied stress for crack initiation becomes maximum. Quite the contrary happens for the shear loading system which causes the largest plastic zones and the minimum applied stress values fro crack growth. Results showing the dependence of the above quantities on the biaxiality of the applied stress are presented in graphical form.     

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.     

Determination of stress intensity factor with direct stress approach using finite element analysis

In this article,a direct stress approach based on finite element analysis to determine the stress intensity fac-tor is improved.Firstly,by comparing the rigorous solution against the asymptotic solution for a problem of an infinite plate embedded a central crack,we found that the stresses in a restrictive interval near the crack tip given by the rigorous solution can be used to determine the stress intensity fac-tor,which is nearly equal to the stress intensity factor given by the asymptotic solution.Secondly,the crack problem is solved numerically by the finite element method.Depending on the modeling capability of the software,we designed an adaptive mesh model to simulate the stress singularity.Thus, the stress result in an appropriate interval near the crack tip is fairly approximated to the rigorous solution of the corre-sponding crack problem.Therefore,the stress intensity factor may be calculated from the stress distribution in the appro-priate interval,with a high accuracy.     

Interface crack between a compressible elastomer and a rigid substrate with finite slippage

We study the deformation of a crack between a soft elastomer and a rigid substrate with finite interfacial slippage. It is assumed that slippage occurs when the interfacial shear traction exceeds a threshold. This leads to a slip zone ahead of the crack tip where the shear traction is assumed to be equal to the constant threshold. We perform asymptotic analysis and determine closed-form solutions describing the near-tip crack opening displacement and the corresponding stress distributions. These solutions are consistent with numerical results based on finite element analysis. Our results reveal that slippage can significantly affect the deformation and stress fields near the tip of the interface crack. Specifically, depending on the direction of slippage, the crack opening profile may appear more blunted or sharpened than the parabola arising from for the case of zero interfacial shear traction or free slippage. The detailed crack opening profile is determined by the constant shear traction in the slip zone. More importantly, we find that the normal stress perpendicular to the interface can increase or decrease when slippage occurs, depending on the direction of slippage and the shear traction in the slip zone.