CT试样三维疲劳裂纹扩展数值模拟

在循环载荷下疲劳裂纹的裂纹形貌在稳定扩展区近似为半椭圆形状,因此通过Paris方法根据疲劳裂纹表面尖端点应力强度因子的变化幅值(△K)得到扩展速率与真实的裂纹速率会有误差.为了更好的研究疲劳裂纹的性质,本文通过分析紧凑拉伸(CT)试样的疲劳裂纹扩展后的三维形貌,采用Jiang-Sehitoglu循环塑性模型和疲劳准则以...     

过载作用下聚乙烯疲劳裂纹扩展行为研究

论文针对中密度聚乙烯材料(MDPE)，采用平板试样进行了I型疲劳裂纹扩展和单次过载下裂纹扩展试验．发现与金属材料类似，单次拉伸过载对聚乙烯(PE)的疲劳裂纹扩展有明显的迟滞作用，降低了裂纹扩展速率．试验还通过变载荷刻线法获取疲劳裂纹扩展前缘的实际形貌和变化规律，对常规变载荷刻线方法进行了调整和验证，其修正方法对高分子材料的疲劳裂纹扩展前缘刻线具有较好的效果．通过观察发现含楔形塑性区的裂尖钝化是裂纹迟滞的主要原因．过载引入的塑性区内残余应力对裂纹迟滞也起了重要作用．论文利用Dugdale模型计算了塑性区尺寸，使用基于残余应力的Wheeler模型对过载迟滞进行了很好的拟合．     

用可动节点有限元模型模拟裂纹扩展过程

本文针对裂纹扩展有限元模拟模型存在的问题,提出一种可动节点模型.使裂尖在单元尺度内任意移动,实现了扩展模拟连续性.加进松弛扩展,突出其所占比重,使模拟更加逼真.用弹塑性有限变形杂交元自编模拟裂纹扩展的通用FORTRAN 程序,对含中心裂纹铝板进行模拟,假定材料各向同性并遵守等向强化规律,采用米赛斯屈服准则和Drucker 塑性势流动理论做增量弹塑性静力分析,得到了与实验符合较好的结果.     

考虑材料循环塑性的疲劳裂纹扩展模拟

提出了一种考虑材料循环塑性性能的研究疲劳裂纹扩展与闭合行为的有限元模拟方法．对所选用的循环塑性本构关系进行了基本实验检验．探讨了在疲劳裂纹扩展有限元分析中网格尺寸的影响，给出了网格优化准则．研究了在循环硬化条件下考虑裂纹合效应时裂纹面张开廓形、裂纹尖端应力、应变场和正反向塑性区的演变规律．对于循环硬化和不同循环应力比Ｒ等因素对裂纹张开应力水平的影响也作了考察     

CTS试件中复合型疲劳裂纹扩展

针对复合型循环载荷作用下的金属构件中的裂纹扩展问题进行了实验分析和理论建模. 首先 采用紧凑拉剪试件(CTS)和 Richard研制的复合型载荷加载装置，对承受复合型循环载荷的裂纹进行了实验研究. 实验选择了两种金属材料试件，分别承受3种形式的复合型循环载荷的作用，在裂纹尖端具 有相同的初始应力场强度的条件下考察复合型循环载荷对裂纹扩展规律的影响. 实验结果表明，疲劳裂纹的扩展速率与加载角度有关. 对于同样金属材料的试件，当裂尖处 初始应力场强度相等时，载荷越接近于II型，裂纹增长速率越快. 采用等效应力强度 因子(I型和II型应力强度因子的组合)、裂纹扩展速率及复合强度等参数，以实验数据为 基础，建立了一个疲劳裂纹扩展模型，用来预测裂纹在不同模式疲劳载荷作用下的扩展速率. 为验证其有效性，该模型被应用于钢制试件的数值模拟计算中. 实验结果与模拟计算曲线保 持一致，表明该模型可以用来估算带裂纹金属构件的寿命.     

CRACK FACE CONTACT PROBLEM ANALYSIS USING THE SCALED BOUNDARY FINITE ELEMENT METHOD

比例边界有限元侧面上有任意荷载时,将侧面载荷分解成关于径向方向局部坐标的多项式函数的和,推导给出了考虑侧面载荷存在的新型形函数,并基于该形函数推导了刚度矩阵和等效节点载荷列阵.首次对比例边界有限元法求解裂纹面接触问题进行了研究,运用Lagrange乘子引入接触界面约束条件,推导给出了比例边界有限元求解裂纹面接触问题的控制方程.将裂纹面单元分为非裂尖单元和含有侧面的裂尖单元.在非裂尖单元中的裂纹面,裂纹面作为多边形单元的边界,边界上的接触力可等效到节点上,通过在节点上构造Lagrange乘子,采用点对点接触约束进行处理.对于含有侧面的裂尖单元,在整个侧面上构造Lagrange乘子的插值场,采用边对边接触约束进行处理.对三个不同的接触约束状态下的算例进行了数值计算,通过与解析解及有限元软件ABAQUS计算结果的对比,验证了本文提出的比例边界有限元点对点和边对边接触求解裂纹面接触问题的精确性与有效性.     

比例边界有限元法求解裂纹面接触问题

比例边界有限元侧面上有任意荷载时,将侧面载荷分解成关于径向方向局部坐标的多项式函数的和,推导给出了考虑侧面载荷存在的新型形函数,并基于该形函数推导了刚度矩阵和等效节点载荷列阵.首次对比例边界有限元法求解裂纹面接触问题进行了研究,运用Lagrange乘子引入接触界面约束条件,推导给出了比例边界有限元求解裂纹面接触问题的控制方程.将裂纹面单元分为非裂尖单元和含有侧面的裂尖单元.在非裂尖单元中的裂纹面,裂纹面作为多边形单元的边界,边界上的接触力可等效到节点上,通过在节点上构造Lagrange乘子,采用点对点接触约束进行处理.对于含有侧面的裂尖单元,在整个侧面上构造Lagrange乘子的插值场,采用边对边接触约束进行处理.对三个不同的接触约束状态下的算例进行了数值计算,通过与解析解及有限元软件ABAQUS计算结果的对比,验证了本文提出的比例边界有限元点对点和边对边接触求解裂纹面接触问题的精确性与有效性.     

正交异性钢桥面板肋-面板焊缝疲劳分析

基于焊缝的局部三维断裂力学模型和超重多轴货车的载荷谱,进行正交异性钢桥面板的肋-面板焊缝表面裂纹的疲劳寿命分析。采用Schwartz-Neuman交替法计算肋-面板焊缝处半椭圆表面裂纹的应力强度因子,基于裂尖反向塑性区模型考虑循环载荷中压应力对疲劳裂纹扩展的作用。正交异性钢桥面板的肋-面板焊缝的应力计算结果表明:超载货车作用下肋-面板焊缝处的横向应力峰值和应力幅都有明显增加;相比于标准疲劳荷载车,超载货车作用下肋-面板焊缝处半椭圆表面裂纹的裂纹扩展率增大了6.1倍;对应于正交异性钢桥面板的肋-面板焊缝处的拉-压循环应力,平面应变状态下的裂尖反向塑性区使裂纹扩展率增加了3.7倍;基于所得裂纹扩展速率,本研究给出仅在严重超载的五轴和六轴货车作用下正交异性钢桥面板的肋-面板焊缝疲劳寿命不足20年,远远低于桥梁的设计寿命。因此,考虑超载多轴货车的载荷谱和循环载荷中的压应力对肋-面板焊缝疲劳裂纹扩展的影响十分重要。     

含缺口试件疲劳全寿命预测一种建议的方法

传统的研究含缺口构件的疲劳的方法是将疲劳启裂和疲劳裂纹扩展两个过程完全独立起来,用不同的方法来模拟,相互间并没有定量的关系。本文是基于最新发展的多轴疲劳损伤理论,建立了一种适用于各种载荷条件下的疲劳启裂和裂纹扩展的普适方法。根据从弹塑性分析中得到的应力应变,确定疲劳损伤模型,建立能够预测疲劳启裂、裂纹扩展速率和扩展方向的新方法。整个模拟可以分为两步:弹-塑性应力分析得到材料的应力应变分布;再运用一个通用的疲劳准则预测疲劳裂纹启裂和裂纹扩展。通过对1070号钢含缺口试件的疲劳全寿命预测,得到了与实验非常吻合的模拟结果。     

疲劳裂纹塑性闭合特性的数值分析方法

塑性诱导裂纹闭合是导致裂纹闭合发生的主要机理之一.利用弹塑性有限元法模拟含中心裂纹矩形板试件的疲劳裂纹扩展,并确定疲劳裂纹张开、闭合应力水平.通过计算,考察应力比R、裂纹长度、最大应力强度因子Kmax等对疲劳裂纹张开闭合应力的影响规律.论文阐述了所采用的裂纹扩展模拟方法及确定裂纹张开和闭合应力的原理.采用了等K加载方式,即在裂纹扩展中裂尖应力强度因子的最大值Kmax保持不变(给定R比,最大应力σmax随裂纹长度变化).分析了两种Kmax水平下R比分别为0.3,0,-0.5和-1.0共8种载荷工况.结果表明,对各个载荷工况,用瞬时最大应力σmax正则化的裂纹张开、闭合应力水平σop/σmax和σcl/σmax与裂纹长度无关.等K循环加载比等幅循环加载更有利于分析影响裂纹闭合水平的因素和闭合效应对疲劳裂纹扩展的影响规律,为建立基于裂纹闭合效应的疲劳裂纹扩展规律模型提供了一种新的思路.     

Modeling of fatigue crack growth from a notch

When a fatigue crack is nucleated and propagates into the vicinity of the notch, the crack growth rate is generally higher than that can be expected by using the stress intensity factor concept. The current study attempted to describe the crack growth at notches quantitatively with a detailed consideration of the cyclic plasticity of the material. An elastic–plastic finite element analysis was conducted to obtain the stress and strain histories of the notched component. A single multiaxial fatigue criterion was used to determine the crack initiation from the notch and the subsequent crack growth. Round compact specimens made of 1070 steel were subjected to Mode I cyclic loading with different R-ratios at room temperature. The approach developed was able to quantitatively capture the crack growth behavior near the notch. When the R-ratio was positive, the crack growth near a notch was mainly influenced by the plasticity created by the notch and the resulted fatigue damage during crack initiation. When the R-ratio was negative, the contact of the cracked surfaces during a part of a loading cycle reduced the cyclic plasticity of the material near the crack tip. The combined effect of notch plasticity and possible contact of cracked surface were responsible for the observed crack growth phenomenon near a notch.     

A reexamination of plasticity-induced crack closure in fatigue crack propagation

The crack closure concept is often used to consider the R-ratio and overload effects on fatigue crack growth. The presumption is that when the crack is closed, the external load produces negligible fatigue damage in the cracked component. The current investigation provides a reassessment of the frequently used concept with an emphasis on the plasticity-induced crack closure. A center cracked specimen made of 1070 steel was investigated. The specimen was subjected to plane-stress mode I loading. An elastic–plastic stress analysis was conducted for the cracked specimens using the finite element method. By applying the commonly used one-node-per-cycle debonding scheme for the crack closure simulations, it was shown that the predicted crack opening load did not stabilize when the extended crack was less than four times of the plastic zone size. The predicted opening load was strongly influenced by the plasticity model used. When the elastic–perfectly plastic (EPP) stress–strain relationship was used together with the kinematic hardening plasticity theory, the predicted crack opening load was found to be critically dependent on the element size of the finite element mesh model. For R = 0, the predicted crack opening load was greatly reduced when the finite element size became very fine. The kinematic hardening rule with the bilinear (BL) stress–strain relationship predicted crack closure with less dependence on the element size. When a recently developed cyclic plasticity model was used, the element size effect on the predicted crack opening level was insignificant. While crack closure may occur, it was demonstrated that cyclic plasticity persisted in the material near the crack tip. The cyclic plasticity was reduced but not negligible when the crack was closed. The traditional approaches may have overestimated the effect of crack closure in fatigue crack growth predictions.     

A stress criterion for crack growth

A stress criterion for crack growth was developed from test results with 7075-T6 aluminum-sheet specimens containing transverse machined cracks. Stress distributions near the crack tip were obtained using strain gages and by reducing the strain data to stresses with the aid of Reuss plasticity theory. These distributions indicated the biaxial nature of stress at the crack tip, the high stress gradients a short distance from the tip, and the variation in stress-concentration factor with crack length. Crack growth was found to occur when the effective stress at the crack root reached the engineering ultimate strength.     

Analysis of dynamic growth of a tensile crack in an elastic-plastic material

The influence of inertia on the stress and deformation fields near the tip of a crack growing in an elastic-plastic material is studied. The material is characterized by the von Mises yield criterion and J₂ flow theory of plasticity. The crack grows steadily under plane strain conditions in the tensile opening mode. Features of the stress and deformation state at points near the moving crack tip are described for elastic-perfectly plastic response and for several crack propagation speeds. It is found that inertia has a significant effect on the elastic-plastic response of material particles near the crack tip, and that elastic unloading may occur behind the crack tip for higher speeds. The relationship between the applied crack driving force, represented by a remote stress intensity factor, and the crack tip speed is examined on the basis of a critical crack tip opening angle growth criterion. The calculated result is compared with dynamic fracture toughness versus crack speed data for a 4340 steel.     

Effect of yield surface vertex on crack-tip fields in mode III

An asymptotic crack-tip analysis of stress and strain fields is carried out for an antiplane shear crack (Mode III) based on a corner theory of plasticity. Because of the nonproportional loading history experienced by a material element near the crack tip in stable crack growth, classical flow theory may predict an overly stiff response of the elastic plastic solid, as is the case in plastic buckling problems. The corner theory used here accounts for this anomalous behavior. The results are compared with those of a similar analysis based on the J₂ flow theory of plasticity.     

Stress and energy density fields near a blunted crack front in power hardening material

Using Jaumann and Dienes rates of Euler stress in elastic-plastic constitutive equations of finite deformation, plane strain finite element analysis for a compact tension specimen with a blunted crack front is made. The Euler stress, Kirchhoff stress and volume strain energy density near a blunted crack tip are computed. Constitutive relations with different deformation rates affect the the near crack tip solution in a region within an order of magnitude of the crack opening displacement. The results differed from the corresponding solution of deformation plasticity (or nonlinear elasticity) with increasing deformation. They are smaller in a local region of about 2 to 10 times of the crack opening distance.The volume energy density near the crack tip is computed, the stationary values of which determine the locations of extensive yielding and possible sites of crack initiation. It remained nearly constant with increasing deformation. Such a character tends to support the volume energy density criterion as a means for quantifying the ductile fracture behavior of metals.     

Analyses of steady quasistatic crack growth in plane strain tension in elastic-plastic materials with non-isotropic hardening

Steady state crack propagation problems of elastic-plastic materials in Mode I, plane strain under small scale yielding conditions were investigated with the aid of the finite element method. The elastic-perfectly plastic solution shows that elastic unloading wedges subtended by the crack tip in the plastic wake region do exist and that the stress state around the crack tip is similar to the modified Prandtl fan solution. To demonstrate the effects of a vertex on the yield surface, the small strain version of a phenomenological J₂, corner theory of plasticity (Christoffersen, J. and Hutchinson, J. W. J. Mech. Phys. Solids,27, 465 C 1979) with a power law stress strain relation was used to govern the strain hardening of the material. The results are compared with the conventional J₂ incremental plasticity solution. To take account of Bauschinger like effects caused by the stress history near the crack tip, a simple kinematic hardening rule with a bilinear stress strain relation was also studied. The results are again compared with the smooth yield surface isotropic hardening solution for the same stress strain curve. There appears to be more potential for steady state crack growth in the conventional J₂ incremental plasticity material than in the other two plasticity laws considered here if a crack opening displacement fracture criterion is used. However, a fracture criterion dependent on both stress and strain could lead to a contrary prediction.     

NiTi形状记忆合金裂纹尖端热力耦合行为研究

本文对NiTi形状记忆合金I型裂纹尖端热力耦合行为进行了数值仿真分析和实验验证。建立了包含相变和热力耦合的本构模型,通过有限元计算得到了裂纹尖端附近的纵向应变、马氏体体积分数和温度场分布,依据马氏体相变情况对裂纹尖端有效应力强度因子进行了修正,揭示了加载速率对形状记忆合金裂纹尖端有效应力强度影子的影响规律。参数研究表明,随着加载频率的增加,裂纹尖端附近温度逐渐升高,马氏体相变区域逐渐缩小,有效应力强度因子呈下降趋势,形状记忆合金表现出增韧效应,有助于减缓裂纹扩展。本研究结果对于揭示热力耦合作用下超弹性形状记忆合金疲劳裂纹扩展规律具有重要参考意义。     

孔洞对平面裂纹应力强度因子影响分析的广义参数Williams单元

本文采用圆形奇异区广义参数Williams单元（W单元）建立了中心裂纹与圆孔共存的平面应力模型,奇异区外围利用ABAQUS有限元软件自动网格离散技术与FORTRAN95编程前处理相结合,克服了自主编程中网格离散的局限性．算例分析了圆孔位置和几何参数对I-II混合型裂纹尖端应力强度因子(SIFs)的影响,并与扩展有限元法(XFEM)计算结果进行比较．结果表明：靠近圆孔一侧的裂尖SIFs大于远离圆孔一侧的裂尖SIFs;控制圆孔左边缘到裂纹中心的距离,则两侧裂尖SIFs随圆孔半径的增大而增大;圆孔中心与裂纹中心水平距离越远,圆孔对裂纹扩展的影响越小．同时,基于圆形奇异区的W单元直接计算得到的裂尖SIFs与扩展有限元法得到的解吻合较好,证明了W单元对奇异区离散形状不敏感,且具有高效率和高精度．