CSM玻璃钢复合型断裂的有限元分析

对切短玻璃纤维毡增强聚脂层板的复合型断裂进行了有限元分析。采用八节点四边形等参元的正规型式计算应力分布与损伤区扩展;而用坍塌(collapsed)三角形四分之一点(quarter-point)奇异元计算应力强度因子K_Ⅰ与K_Ⅱ。用节点位移约束与释放技术计算了裂纹扩展过程。对决定应力强度因子K_Ⅰ与K_Ⅱ的三种方法进行了对比。对施加于裂纹顶点节点约束条件的影响进行了评价。最后基于实验测得的断裂载荷与临界裂纹长度,估算了材料在这种受力条件下的复合型临界应力强度因子K_ⅠC和K_ⅡC。     

任意形状孔口单边裂纹问题的边界配位解法

本文提出了一组复应力函数,采用边界配位方法对不同形状孔口(包括圆、椭圆、矩形及菱形孔口)的单边裂纹平板的应力强度因子进行了计算.计算结果表明,对长度和宽度远大于孔口和裂纹几何尺寸的试件,配位法与用其他方法所得的无限大板含圆或椭圆孔边裂纹问题的解符合得很好.同时,对其他孔口问题,特别是有限大板情形,本文给出了一系列计算结果.本文所提出的函数及计算过程可以应用于任意形状孔口单边裂纹平板的计算.     

二维多裂纹弹性体的有效弹性模量

采用严格考虑微裂纹相互作用的分析方法 ,针对微裂纹随机分布和平行分布两种情况 ,计算了无限大体中代表性体元 (RVE)的有效弹性模量 ,并与其他细观损伤理论所得结果和试验结果进行了比较 .数值计算结果表明 ,所用方法简单精确 ,对求解多裂纹体问题非常有效     

分析复合材料层合板弯曲和振动的一种有效无网格方法

基于高阶剪切法向变形板理论（HOSNDPT）利用无网格方法对层合板弯曲和振动问题进行数值分析.在通常的径向点插值法（RPIM）中对每个Gauss(高斯)点或计算点需要求矩矩阵的逆,且受到影响域半径大小的限制.而在加权节点径向点插值法（WN-RPIM） 近似中,求解系统矩阵的逆的数量等于问题域中的节点数量,它远远小于Gauss点的数目,可以大大减少矩矩阵求逆的计算量,且克服了RPIM中影响域半径大小的限制.首先,将三维板位移分解成厚度和面内位移的乘积,在厚度方向使用正交Legendre多项式作为基函数,在板的面内使用WN RPIM来构造形函数.然后,通过对层合板的弯曲问题进行数值计算表明WN-RPIM的计算精度和稳定性.最后,将该方法推广到对不同边界条件、不同厚跨比、不同铺设方式的层合板振动问题的数值计算,数值结果表明了本文提供方法的适用性和有效性.     

外部压力作用下复合材料层合圆柱壳的脱层扩展分析

含脱层的层合圆柱壳在承受外部压力作用时,容易发生脱层扩展,进而引起结构失效．基于可动边界变分原理并考虑脱层间的接触效应对圆柱壳在外部压力作用下的脱层扩展进行了分析,同时应用Griffith准则,导出了脱层前缘各点处的能量释放率表达式,且以轴对称脱层层合圆柱壳为例进行了数值计算,讨论了脱层大小、脱层深度、几何尺寸、材料性质及纤维铺层方式等因素对脱层扩展的影响．     

任意形状孔口双边裂纹平板的应力强度因子计算

本文采用Muskhclishvili弹性力学的复变函数和边界配位方法对不同形状孔口双边裂纹问题进行了研究,计算了圆孔、椭圆孔、矩形孔、菱形孔等不同形状孔口双边裂纹,以及Ⅰ型和复合型等不同类型断裂试件的应力强度因子,本文方法简单方便,精度较高,与某些已有计算结果的问题比较,本文方法所得的结果是令人满意的.同时,本方法可以应用于不同几何形状和加载条件下的孔口双边裂纹有限大板的计算,是解这一类问题的一致有效方法.     

非线性塑性材料内部空洞闭合模型的研究

基于典型体元(RVE)模型和Rayleigh-Ritz法,对材料内部空洞从球形闭合成裂纹的过程进行了定量研究.基体材料的本构关系采用幂次粘性方程.通过研究材料变形过程中内部球形空洞和圆形裂纹的演化规律,得到了各自的体积应变率的表达式.采用插值近似,建立了空洞闭合的解析模型,发现空洞变形的主要机理来自于空洞周围基体材料的塑性流动.空洞闭合模型反映了材料属性、远场应力三轴度、远场等效应变对空洞闭合的定量规律.空洞闭合模型的预测结果与文献中的数值结果和有限元计算结果相吻合.空洞闭合模型与CAE(computer-aided engineering)技术相结合可对材料加工工艺进行优化设计,为消除材料内部空洞提供了一条有应用前景的新途径.     

覆冰输电线路脱冰跳跃高度理论算法

建立连续档线路脱冰理论分析简化模型,基于导线脱冰跳跃过程中的能量关系、应力弧垂关系、变形关系和平衡关系,给出求解导线最大脱冰跳跃高度的理论计算方法.利用算例与数值模拟方法和Morgan理论算法进行比较,验证了方法的正确性.进而针对典型输电线路的脱冰跳跃问题,利用获得的理论算法、数值模拟方法与现行输电线路设计规程中的脱冰跳跃高度计算公式的计算结果进行比较分析,讨论了该理论算法的实用性.采用所获得的冰跳高度理论算法,可在冰区输电线路的设计中快速确定导线冰跳时的绝缘间隙.     

考虑表面效应时孔边均布径向多裂纹Ⅲ型断裂力学分析

理论研究了纳米尺度孔边均布径向多裂纹的Ⅲ型断裂性能.基于Gurtin-Murdoch表面弹性理论和保角映射技术,获得了孔和裂纹应力场的解析解,给出了裂纹尖端应力强度因子的闭合解.基于解答分析了应力强度因子的尺寸效应,讨论了裂纹数量、裂纹/孔径比和缺陷表面性能对应力强度因子的影响.结果表明:当孔和裂纹尺寸在纳米量级时,无量纲应力强度因子具有显著的尺寸效应;应力强度因子随裂纹数量的变化规律受裂纹/孔径比的影响;裂纹/孔径比对应力强度因子的影响受到缺陷表面性能的制约,同时表面性能对应力强度因子的影响也受限于裂纹/孔径比;表面效应对应力强度因子的影响与裂纹数量无关.     

带复合型摩擦裂纹的圆盘动态断裂实验有限元分析

为验证考虑裂纹面接触和动态荷载时,中心裂纹巴西圆盘（CCBD）试件用于分离式Hopkinson压杆（SHPB）系统中测量脆性材料复合型动态断裂韧度的可行性,以及研究裂纹面接触对动态断裂韧度实验结果的影响.通过有限元法建立SHPB CCBD三维有限元模型,计算了不同加载条件下CCBD试件的动态应力强度因子（DSIF）.结果表明:在实验中,将考虑裂纹面接触的应力强度因子(SIF)准静态公式推广为动态公式,需要判定断裂时间是否达到应力平衡的时间条件;压剪复合型加载时,裂纹面接触导致裂纹面应力变化,会对Ⅱ型裂纹的DSIF产生显著影响,不考虑裂纹面接触的影响将会导致Ⅱ型DSIF的测试值偏大.     

压电复合材料圆形夹杂中螺型位错和界面裂纹的干涉分析

研究了无穷远纵向剪切和面内电场共同作用下,压电复合材料圆形夹杂中螺型位错与界面裂纹的电弹耦合干涉作用．运用Riemann-Schwarz 对称原理,并结合复变函数奇性主部分析方法,获得了该问题的一般解答．作为典型算例,求出了界面含一条裂纹时基体和夹杂区域复势函数和电弹性场的封闭形式解．应用广义Peach-Koehler公式,导出了位错力的解析表达式．分析了裂纹几何参数和材料的电弹性常数对位错力的影响规律．结果表明,界面裂纹对位错力和位错平衡位置有很强的扰动效应,当界面裂纹长度达到临界值时,可以改变位错力的方向．该结果可以作为格林函数研究圆形夹杂内裂纹和界面裂纹的干涉效应．其公式的退化结果与已有文献完全一致．     

Crack propagation calculation using trial cracks

We present a method to perform finite element calculations for crack propagation problems with arbitrary crack directions in two dimensions. The crack direction (angle of propagation) is determined by inserting small “trial cracks” at the crack tip. For each trial crack, the domain is remeshed to allow crack propagation between elements. The trial cracks are then opened and the energy release rate is measured. The optimum crack direction (i.e., the crack direction with maximum energy release) is determined by an optimisation procedure. Although the method is computationally expensive due to the need to perform several calculations for each crack increment, it has the advantage that the energy release rate can be calculated even in cases where other methods fail. After explaining the method, it is applied to some test examples. (© 2008 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

点群6一维六方准晶椭圆孔口与半无限裂纹反平面弹性问题的解析解

导出了点群6-维六方准晶反平面弹性问题的控制方程.利用复变方法,给出了点群6-维六方准晶在周期平面内的反平面弹性问题的应力分量以及边界条件的复变表示,通过引入适当的保角变换,研究了点群6-维六方准晶中带有椭圆孔口与半无限裂纹的反平面弹性问题,得到了椭圆孔口问题应力场的解析解,给出了半无限裂纹问题在裂纹尖端处的应力强度因子的解析解.在极限情形下,椭圆孔口转化为Griffith裂纹,并得到该裂纹在裂尖处的应力强度因子的解析解.当点群6-维六方准晶体的对称性增加时,其椭圆孔口与半无限裂纹的反平面弹性问题的解退化为点群6mm-维六方准晶带有椭圆孔口与半无限裂纹的反平面弹性问题的解。     

Structural Crack Identification Based on the Variational Mode Decomposition北大核心CSCD

In order to enrich the bridge damage detection method and further improve the accuracy of bridge damage identification, a detection method for simply supported beams with cracks under dynamic loads was proposed not based on the complete finite element model. Under the premise of not blocking traffic, the method only needs to analyze and deal with the acceleration responses of the simply supported beam span, which reduces the mounting, dismounting and maintenance of sensors in practical engineering. At the same time, based on the model, an analytical formula of the acceleration at the midspan of the simply supported cracked beam was derived. Based on the theoretical derivation, the instantaneous energy and the mean energy difference were constructed through the variational mode decomposition and the Hilbert transform, and these 2 crack identification indexes were used to effectively identify small cracks with a crack depth ratio of only 5%. Then the influences of different wheel loads, environmental noises and damage degrees on detection results were studied. The results show that: ① the instantaneous frequency has a better recognition effect for crack positions; ② the mean energy difference is sensitive to crack depth ratio δ and the wheel load magnitude; ③ this method has strong noise robustness. © 2022 Editorial Office of Applied Mathematics and Mechanics. All rights reserved.     

PROCRACK-A Software Tool for Finite Element Simulation of Three-Dimensional Fatigue Crack Growth

The lifetime of structural components is limited by fatigue cracks. After initiation from an existing defect the crack grows subcritically until it reaches a critical size. At this point it becomes unstable and the structural component fails. PROCRACK is a powerful tool that enables the commercial finite element software Abaqus to calculate the crack propagation in pre-cracked components. The complete capability of Abaqus can be used to simulate nearly arbitrary geometries. Abaqus/CAE is used for the three-dimensional modeling of the initial crack at a geometrical level by means of points, lines and triangles. The numerical analysis is performed by Abaqus/Standard. PROCRACK automatically generates a tube-shaped submodel around the crack front to calculate the stress intensity factors with high accuracy. The Paris law or the NASGRO equation can be used to calculate the incremental crack growth. The shape of the crack and the finite element mesh are updated in every crack growth step. (© 2014 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

带裂纹三点弯曲试样的动态应力强度因子分析

提出了计算带单边裂纹三点弯曲试样动态应力强度因子的新方法．首先由权函数的普遍形式和两种参考载荷下的应力强度因子,得到了带单边裂纹三点弯曲试样的权函数,然后考虑试样的转动惯性和剪切变形,根据振动理论推导出无裂纹梁内的动应力响应和分布,最后由权函数的思想推导出了带裂纹三点弯曲试样动态应力强度因子公式．通过有限元数值计算,验证了该方法的正确性,结果比较表明公式具有较高的精度．另外,还研究了冲击载荷下三点弯曲试样的动态应力强度因子随裂纹长度和加载速率的变化规律．     

带裂纹十次对称二维准晶平面弹性的无摩擦接触问题

借助经典平面弹性复变函数方法，研究了单个刚性凸基底压头作用下，带任意形状裂纹十次对称二维准晶半平面弹性的无摩擦接触问题.利用十次对称二维准晶位移、应力的复变函数表达式, 带任意形状裂纹的准晶半平面弹性无摩擦接触问题被转换为可解的解析函数复合边值问题，进而简化成一类可解的Riemann边值问题.通过求解Riemann边值问题，得到了应力函数的封闭解, 并给出了裂纹端点处应力强度因子和压头下方准晶体表面任意点处接触应力的显式表达式.从压头下方接触应力的表达式可以看出, 接触应力在压头边缘和裂纹端点处具有奇异性.当忽略相位子场影响时, 该文所得结论与弹性材料对应结果一致.数值算例分别给出了单个平底刚性压头无摩擦压入带单个垂直裂纹和水平裂纹的十次对称二维准晶下半平面的结果.该文所得结论为准晶材料的应用提供了理论参考.     

A Variational–Based Formulation of Regularized Brittle Fracture

The paper performs a comparative study of variational-based brittle fracture with its gradient-type regularization, and outlines aspects of the numerical implementations of both approaches. The latter smoothes out sharp displacement discontinuities of cracks. On the side of discrete crack modeling, we propose a variational framework of configurational-force-driven crack propagation. The latter provides the basis for the computation of material nodal forces and drives the crack propagation in our proposed finite element framework with adaptive nodal doubling. Such a formulation is of limited applicability for the modeling of crack inititation in homogeneous bodies without defects and in situations with complex crack branching. This can be overcome by a regularized crack modeling. Here, an elliptic approximation of the crack surface term yields a regularized two field functional, where an additional scalar field approximates the set of discontinuities. We provide a simple interpretation of such a transition from the sharp crack to the regularized setting. It results in a smooth continuum-damage-type theory with a specific gradient-damage and hardening terms, depending on a length scale that represents the width of a zone that surrounds the crack. Such a variational framework is implemented by a coupled two-field finite element framework in a staightforward manner. We compare representative numerical results obtained by both methods for selected crack patterns and highlight the pro and contra of both meshes. (© 2009 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Detection of matrix cracks in composite laminates by using the modal strain energy method

As the foremost mechanism of damage development, matrix cracking is the critical damage found in the early stage of structural failure of composites. This study aims to nondestructively detect matrix cracks in composite laminates by using an experimental modal analysis (EMA). An AS4/PEEK composite was used to fabricate cross-ply [0₂/90₁₂/0₂] and quasi-isotropic [(±45/0/90)₄] _s laminates. The damage in the form of a matrix crack in the laminates was created by using a tensile load. The EMA was conducted on the laminates to obtain the modal displacements before and after damage. The displacements were then employed to compute the modal strain energy and to define the damage index, which is used for detecting matrix cracks. Limited by the mesh points of measurements, we used the differential quadrature method to calculate the partial differentials in the strain energy formula. The results obtained were validated by using the X-ray radiography method and three-point bending tests. The experimental results showed that the damage index well identified the location of breadthwise matrix cracks inside the laminates. However, the resolution of the damage index became poor if the spans of matrix cracks were short or the matrix cracks were scattered over the laminates.     

On the calculation of energy release rate for curved cracks of rubbery materials

The energy release rate G for 2-D rubbery material problems with curved cracks is calculated with finite element solutions in this paper. Two approaches, a generalized domain integral method and the virtual crack extension method, are investigated. The generalized domain integral method is demonstrated to be “patch independent” and therefore a complicated finite element mesh adjacent to the crack tip area is not required.