数学网格和物理网格分离的有限单元法(II):粘聚裂纹扩展问题中的应用

强化有限单元法将物理网格与数学网格分离开来,可以方便地描述非连续变形;粘聚区域模型是模拟断裂过程区作用最简单有效的方法,且可以避免裂纹尖端的应力奇异性.本文以平面问题为例,将强化有限单元法与粘聚区域模型相结合,利用富集数学节点描述任意粘聚裂纹扩展过程中的非连续变形问题,提出了裂纹扩展过程中数学节点富集和数学单元定义的方法.本文还导出了与平面4～8节点平面等参单元对应的8～16节点粘聚裂纹单元列武.最后,通过三点弯梁的裂纹扩展过程模拟验证了本文提出的粘聚裂纹扩展模拟方法的有效性.     

双K断裂模型粘聚韧度KcIc实用插值计算方法

为了配合《水工混凝土断裂试验规程》在工程界的推广使用,提出了双K断裂模型中在临界状态时分布在虚拟裂缝上的粘聚力σ(x)产生的粘聚韧度KcIc的二元拉格朗日插值计算方法,进而由准脆性材料韧度三参数定律计算了混凝土起裂韧度KiInci。与粘聚韧度KcIc积分计算和简化计算方法的比较,表明此插值计算方法具有很高的精度,计算过程简化,使用计算器即可实现,便于工程人员实际使用。     

混凝土断裂过程及尺寸效应分析

为研究混凝土裂纹断裂过程和最大承载力计算方法,通过实验机对四种不同尺寸混凝土紧凑拉伸断裂试件进行了加载过程实验,并对其中一个试件进行应变片跟踪测试.由实验结果分析得到了一系列关系曲线,如试件的载荷-加载点位移关系曲线,断裂损伤区变形随载荷变化曲线;并且计算了不同尺寸断裂试件的应力强度因子.结合计算粘聚裂纹应力强度因子的公式与断裂准则,完成了对承载力理论值的计算,并将其与实验峰值平均值进行对比,其结果是两者相比误差较小,表明此种计算裂纹结构最大承载力方法是可行的.     

整体加筋壁板裂纹扩展轨迹模拟及控制分析

应用有限元方法并结合裂纹扩展方向判据分析模拟了飞机整体加筋壁板裂纹的扩展轨迹,与实验结果对比证明此方法的有效性;并通过壳单元模型和粘聚区渐进损伤模型详细讨论了影响整体加筋壁板裂纹扩展的筋条参数,研究了壁板筋条对裂纹止裂的控制,并对不同参数下裂纹的应力强度因子和扩展情况作对比分析.文章所述方法和结论对于飞机整体加筋壁板的损伤容限设计具有指导意义.     

双K断裂模型粘聚韧度K_(Ic)~c实用插值计算方法

为了配合<水工混凝土断裂试验规程>在工程界的推广使用,提出了双K断裂模型中在临界状态时分布在虚拟裂缝上的粘聚力σ(x)产生的粘聚韧度K_(Ic)~c的二元拉格朗日插值计算方法,进而由准脆性材料韧度三参数定律计算了混凝土起裂韧度K_K~(ini).与粘聚韧度K_(Ic)~c积分计算和简化计算方法的比较,表明此插值计算方法具有很高的精度,计算过程简化,使用计算器即可实现,便于工程人员实际使用.     

压电复合材料粘接界面断裂有限元模拟

根据数字化FRMM(Fix-Ratio Mix-Mode)断裂试验,得到了压电复合材料试件的断裂韧性和位移及应变场。本文在试验的基础上,通过非线性有限元软件ABAQUS及用户子程序UMAT进行了模拟分析,采用基于损伤力学的粘聚区模型(CZM)对压电复合材料界面的起裂和脱胶扩展进行了分析,并与VCCT方法进行了比较。计算得到的荷载位移曲线更接近于试验结果,但在裂纹扩展路径上的吻合需要对粘聚区法则进一步修正。通过进一步对CZM参数进行分析,表明界面粘结强度和界面刚度对计算结果的影响很大。研究结果表明,粘聚区模型可以很好地表征压电复合材料弱粘接界面脱胶断裂问题。     

基于扩展有限元的地质聚合物混凝土断裂过程研究

扩展有限元法是基于常规有限元框架分析裂纹等不连续力学问题的一种有效数值方法,在常规的有限元位移表达式中,增加了能够反映位移不连续性的跳跃函数和渐进缝尖位移场函数来对不连续结构附近的节点自由度进行局部加强。本文介绍了扩展有限元法及粘聚力模型的基本原理,给出了基于扩展有限元法的地质聚合物混凝土断裂过程分析方法。分别采用四种不同的软化曲线对I型缺口地质聚合物混凝土梁从裂纹萌生、扩展直至断裂破坏的全过程进行了模拟,并基于双K断裂准则分析了其断裂韧性。结果表明,Petersson模型与试验结果吻合较好,最后基于模拟结果进一步揭示了断裂过程区的演化过程。     

幂硬化粘塑性材料动态扩展裂纹尖端场的渐近解

采用弹性－粘塑性本构模型，对幂硬化粘塑性介质中反平面剪切动态扩展裂纹尖端的应力，应变场进行了渐近分析，给出了反平面剪切动态扩展纹尖端场的渐进方程。分析结果表明，在裂纹法端应力具有（ｌｎＲ／ｒ）１／ｎ－１的奇异性，应变具有（ｌｎＲ／Ｒｎ／ｎ－１的奇异性。从而提示了幂硬化粘塑材料反平面剪动态扩展裂纹尖端场的渐近行为。     

混凝土双K断裂参数计算的半解析有限元法

混凝土裂缝扩展的双$K$断裂准则，用于描述混凝土结构裂缝的起裂、稳定扩 展和失稳断裂. 其相应的双$K$断裂参数(起裂断裂韧度$K_{\rm IC}^{\rm ini} $和失 稳断裂韧度$K_{\rm IC}^{\rm un}$)一般通过简便的试验和基于虚拟裂缝扩展粘 聚力的解析方法确定. 利用平面扇形域哈 密顿体系的方程，通过分离变量法及共轭辛本征函数向量展开法，以解析的方法推导出基于混 凝土虚拟裂缝扩展线性粘聚力模型的平面裂缝解析元列式. 将该解析元与有限元相结合，构成 半解析的有限元法，可求解任意结构几何形状的混凝土平面裂缝双$K$断裂参数的计算问题. 数值计算结果表明半解析有限元法对该类问题的求解是十分有效的.     

各向同性材料接头和界面相交裂纹应力奇异性特征分析

提出了用插值矩阵法分析各向同性材料接头以及与界面相交的平面裂纹应力奇异性。基于接头和裂纹端部附近区域位移场渐近展开,将位移场的渐近展开式的典型项代入线弹性力学基本方程,得到关于平面内各向同性材料接头以及与两相材料界面相交裂纹应力奇异性指数的一组非线性常微分方程的特征值问题,运用插值矩阵法求解,获得了两相材料平面接头端部应力奇异性指数以及与界面以任意角相交的裂纹尖端的应力奇异性指数随裂纹角的变化规律,数值计算结果与已有结果比较表明,本文方法具有很高的精度和效率。     

内聚力模型裂纹问题分析的解析奇异单元

内聚力模型已经被广泛应用于需要考虑断裂过程区的裂纹问题当中,然而常用的数值方法应用于分析内聚力模型裂纹问题时还存在着一些不足,比如不能准确的给出断裂过程区的长度、需要网格加密等。为了克服这些缺点,论文构造了一个新型的解析奇异单元,并将之应用于基于内聚力模型的裂纹分析当中。首先将虚拟裂纹表面处的内聚力用拉格拉日插值的方法近似表示为多项式的形式,而多项式表示的内聚力所对应的特解可以被解析地给出。然后利用一个简单的迭代分析,基于内聚力模型的裂纹问题就可以被模拟出来了。最后,给出二个数值算例来证明本文方法的有效性。     

混凝土黏聚开裂模型若干进展

黏聚模型是用来描述混凝土断裂行为的基本模型, 首先介绍了混凝土的黏聚开裂模型的基本概念,总结了确定黏聚区的本构方程的各种方法,即直接单轴拉伸测试、J积分方法、R曲线法、柔度法和逆推法.然后介绍了黏聚模型在I型和复合型裂纹问题、疲劳断裂问题中的应用以及黏聚模型与混凝土尺寸效应的关系.最后对黏聚开裂模型与桥联模型、带状裂缝模型进行了比较和总结, 指出了该模型存在的问题, 并对其以后的发展方向提出了建议.      

Analysis on the cohesive stress at half infinite crack tip

The nonlinear fracture behavior of quasi-brittle materials is closely related with the cohesive force distribution of fracture process zone at crack tip. Based on fracture character of quasi-brittle materials, a mechanical analysis model of half infinite crack with cohesive stress is presented. A pair of integral equations is established according to the superposition principle of crack opening displacement in solids, and the fictitious adhesive stress is unknown function . The properties of integral equations are analyzed, and the series function expression of cohesive stress is certified. By means of the data of actual crack opening displacement, two approaches to gain the cohesive stress distribution are proposed through resolving algebra equation. They are the integral transformation method for continuous displacement of actual crack opening, and the least square method for the discrete data of crack opening displacement. The calculation examples of two approaches and associated discussions are give     

Fracture process zone analysis of concrete using moiré interferometry

The fracture process zone (FPZ) ahead of a crack tip in concrete and mortar beams subjected to threepoint bending was studied using moiré interferometry. A large FPZ can occur in concrete before the external load reaches its maximum value. Comparing the experimental results between concrete and mortar suggests that the aggregate contributes to the formation of the large FPZ in concrete. The formation of this large FPZ makes concrete less brittle than mortar. The effect of the FPZ on the fracture property, such as stress intensity factor, is investigated by combining moiré interferometry measured displacements with the smoothing FEM method. The study shows that a large FPZ significantly affects the value of the stress intensity factor.     

Dynamic toughness in elastic nonlinear viscous solids

This work addresses the interrelationship among dissipative mechanisms—material separation in the fracture process zone (FPZ), nonelastic deformation in the surrounding background material and kinetic energy—and how they affect the macroscopic dynamic fracture toughness as well as the limiting crack speed in strain rate sensitive materials. To this end, a micromechanics-based model for void growth in a nonlinear viscous solid is incorporated into a microporous strip of cell elements that forms the FPZ. The latter is surrounded by background material described by conventional constitutive relations. In the first part of the paper, the background material is assumed to be purely elastic. Here, the computed dynamic fracture toughness is a convex function of crack velocity. In the second part, the background material as well as the FPZ are described by similar rate-sensitivity parameters. Voids grow in the strain rate strengthened FPZ as the crack velocity increases. Consequently, the work of separation increases. By contrast, the inelastic dissipation in the background material appears to be a concave function of crack velocity. At the lower crack velocity regime, where dissipation in the background material is dominant, toughness decreases as crack velocity increases. At high crack velocities, inelastic deformation enhanced by the inertial force can cause a sharp increase in toughness. Here, the computed toughness increases rapidly with crack velocity. There exist regimes where the toughness is a non-monotonic function of the crack velocity. Two length scales—the width of the FPZ and the ratio of the shear wave speed to the reference strain rate—can be shown to strongly affect the dynamic fracture toughness. Our computational simulations can predict experimental data for fracture toughness vs. crack velocity reported in several studies for amorphous polymeric materials.     

Crack propagation in the power-law nonlinear viscoelastic material

I.IntroductionCrazingdamageisacommonphenon1enonoffractureofpolymericmaterials.Theformationofcrazezoneisamid-stateinthefractureprocessofthematerialsfromperfectstatetofaiIurc.Microscopically,inthisregionthereexistssomefibrilslinkingthetwocracksurfacesandres…     

Analytical representation of the non-square-root singular stress field at a finite angle sharp notch

The stress field near the tip of a finite angle sharp notch is singular. However, unlike a crack, the order of the singularity at the notch tip is less than one-half. Under tensile loading, such a singularity is characterized by a generalized stress intensity factor which is analogous to the mode I stress intensity factor used in fracture mechanics, but which has order less than one-half. By using a cohesive zone model for a notional crack emanating from the notch tip, we relate the critical value of the generalized stress intensity factor to the fracture toughness. The results show that this relation depends not only on the notch angle, but also on the maximum stress of the cohesive zone model. As expected the dependence on that maximum stress vanishes as the notch angle approaches zero. The results of this analysis compare very well with a numerical (finite element) analysis in the literature. For mixed-mode loading the limits of applicability of using a mode I failure criterion are explored.     

Interface crack problems for mode Ⅱ of double dissimilar orthotropic composite materials

The fracture problems near the interface crack tip for mode Ⅱ of double dissimilar orthotropic composite materials are studied. The mechanical models of interface crack for mode Ⅱ are given. By translating the governing equations into the generalized bi-harmonic equations,the stress functions containing two stress singularity exponents are derived with the help of a complex function method. Based on the boundary conditions,a system of non-homogeneous linear equations is found. Two real stress singularity exponents are determined be solving this system under appropriate conditions about himaterial engineering parameters. According to the uniqueness theorem of limit,both the formulae of stress intensity factors and theoretical solutions of stress field near the interface crack tip are derived. When the two orthotropic materials are the same,the stress singularity exponents,stress intensity factors and stresses for mode Ⅱ crack of the orthotropic single material are obtained.     

I-II复合型裂缝断裂角的分层剪滞模型与分析

采用修正的剪滞理论建立了岩石、混凝土等准脆性材料的I-II复合型裂缝在单向拉伸荷载作用下的计算模型,得到了与实验相吻合且优于传统S判据的断裂角。通过对远场应力、斜裂缝区应力以及子层位移的合理简化,得到了求解剪滞分析模型的边界条件,进而得到了含斜裂缝的各子层位移分布函数。引入最大应力集中因子,对I-II复合型裂缝前缘应力场进行简化;基于斜裂缝沿最大应力集中因子方向扩展,得到裂缝的断裂角。根据斜裂缝的应力分布,设置不同的子层分区,得到了更为细化的位移分布模式。通过对计算数据的分析,针对单向拉伸荷载作用下的I-II复合型裂缝,建立了按应力场分区设置子层的分层剪滞模型,得到更为精确的斜裂缝断裂角。