随机超载对疲劳裂纹扩展迟滞效应的模拟

考虑超载的迟滞效应,对随机超载作用下的疲劳裂纹扩展进行了模拟计算．载荷谱为在基本恒幅循环载荷基础上加入一以泊松流发生的随机超载序列,超载的大小为均匀分布．采用裂纹闭合模型考虑超载的迟滞效应,认为裂纹张开应力在超载引起的塑性区内按线性规律衰减．循环续循环模拟计算出裂纹从初始长度一直到疲劳破坏的扩展曲线．据此,计算了各种超载发生强度和大小下的疲劳裂纹扩展寿命的平均值与标准差。     

单元类型和尺寸对裂尖疲劳塑性数值模拟的影响研究

本文采用Jiang-Sehitoglu循环塑性模型和多轴疲劳准则对紧凑拉伸式样裂尖的循环塑性变形、裂纹扩展速率和残余应力进行了有限元数值模拟,着重考察了单元的类型和最小单元尺寸对裂尖循环塑性和裂纹扩展速率的影响.紧凑拉伸试样的材料为1070钢,数值模拟采用了线性单元(四节点)和二次单元(八节点)两种单元,裂尖附近有限元单元的最小尺寸从0.007mm到0.24mm不等.文中将裂纹扩展速率的预测值与实验值进行了比较,通过对裂纹扩展速率的比较,确定在疲劳塑性分析时对单元类型和尺寸进行合理选取.     

裂纹闭合和循环塑性预应变及循环载荷压缩部分对裂纹扩展速率的影响

本文针对三种国产材料 Ｌｙ１１ｃｚ、 Ｌｙ１２ｃｚ 铝合金和 １８ Ｍｎ Ｈ Ｐ钢,通过实验初步考察了循环塑性预应变和循环载荷压缩部分对疲劳裂纹扩展的影响;采用电测法,测定了两种铝合金材料疲劳裂纹扩展的张开应力和有效应力强度因子幅值比 Ｕ。结果表明：（１）材料循环塑性预应变和循环载荷压缩部分,都使疲劳裂纹扩展速率提高;（２）常幅载荷下,在疲劳裂纹稳定扩展阶段,有效应力强度因子幅值比 Ｕ 与应力比 Ｒ 有关,与裂纹长度ａ 无关,并依赖于材料的力学性能。     

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

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

弹塑性疲劳微弯裂纹尖端塑性区问题

主要研究疲劳载荷作用下弹塑性裂纹弯曲延伸扩展问题.通过分析论证,比较精确地研究了疲劳载荷作用下弯曲延伸裂纹尖端塑性区域边界上交变应力的分布状况.综合考虑了疲劳作用应力,塑性区域交变应力,利用二阶摄动方法,研究计算了弯曲延伸裂纹尖端塑性区域的范围,并预测了疲劳载荷作用下弹塑性裂纹扩展路径.     

拉压循环载荷下正交异性钢桥面板肋-面板焊缝裂纹扩展分析

采用四步法计算了考虑循环载荷中压应力影响的正交异性钢桥面板的肋-面板焊缝表面裂纹扩展。第一步是基于正交异性钢桥面板的疲劳分析模型,计算肋-面板焊缝处的应力,第二步是通过肋-面板焊缝的三维局部模型,用Schwartz-Neumann交替法计算焊缝表面裂纹的应力强度因子分布,第三步是用二维断裂力学模型和增量塑性损伤模型,计算循环载荷中的压应力对裂纹扩展的影响,第四步是用第二步中的三维裂纹分析结果和第三步中的二维断裂力学模型得到的裂纹扩展公式,计算钢桥面板的肋-面板焊缝表面裂纹扩展。计算结果表明,对应于正交异性钢桥面板肋-面板焊缝处的循环应力,本文所用模型的裂纹尖端反向塑性区导致裂纹扩展率增加50%以上。研究结果为正交异性钢桥面板肋-面板焊缝裂纹的疲劳寿命分析提供了研究基础。     

一种经验型疲劳小裂纹扩展模型

金属材料疲劳寿命由裂纹萌生和裂纹扩展寿命两部分组成，其中对于萌生寿命中的小裂纹分析是精确描述裂纹萌生寿命的关键．而小裂纹在扩展过程中由于尺寸相对较小，导致传统线弹性断裂力学预测方法失效，需要对其进行改进，考虑裂纹尖端塑性区引起的残余压应力对小裂纹扩展速度的影响．本文针对此问题进行了初步分析，通过对塑性区引起的残余应力的量化，结合小裂纹门槛值特性，提出了一种经验型修正的小裂纹扩展模型，用于定量预测裂纹的萌生寿命．使用铝合金6082-T6缺口试样进行了疲劳实验，并与理论结果进行了对比，验证了所提模型的有效性．     

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

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

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

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

滚动接触下的轨道表面三维裂纹研究方法

针对已有的求解三维裂纹的应力强度因子方法工作量大、难以实现复杂的三维模型建立的问题,提出了一种求解广义移动荷载下三维裂纹尖端应力强度因子的普遍方法。采用包络式三维实体建模方法,有效避免了三维裂纹有限元模型网格划分时因受限于拓扑误差容忍而出现网格划分失败的现象;通过MATLAB三维插值施加节点荷载,解决了不规则三维裂纹有限元模型施加复杂空间荷载的难题。计算结果表明:荷载中心距离裂纹中心约4mm时裂纹尖端应力强度因子达到最大值7.41×107Pa·m1/2;轮轨接触疲劳以张开型为主。此外,给出了最大应力强度因子条件下裂纹尖端塑性区的分布及裂纹扩展情况,可为轮轨滚动接触疲劳的研究及设计提供一定的参考。     

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.     

高分子材料平面应变拉伸变形局部化

将基于应变软化玻璃状高分子材料微观特征建立的ＢＰＡ８－链分子网络模型引入ＵｐｄａｔｉｎｇＬａｇｒａｎｇｅ有限元方法，建立了适于变形局部化分析的大变形弹塑性有限元驱动应力法．在此基础上，数值模拟了初始各向同性高分子材料平面应变拉伸变形局部化的传播过程．探讨了ＢＰＡ模型对具有加工硬化特性的结晶性高分子材料变形分析的适应性；分析了局部化传播过程中颈缩截面的非均匀应力三轴效应；最后，讨论了网格尺寸以及初始几何不均匀性对颈缩扩散以及应力三轴效应的影响     

Finite element analysis of grain-by-grain deformation by crystal plasticity with couple stress

Rigid–plastic crystal plasticity with the rate-sensitive constitutive behavior of a slip system has been formulated within the framework of a two-dimensional finite element method to predict the grain-by-grain deformation of single- and polycrystalline FCC metals. For that purpose, individual grains are represented by several numbers of finite elements to describe the sub-grain deformation behavior, and couple stress has been introduced into the equilibrium equation to be able to describe the size effect as well as to prevent mesh-dependent predictions. A modified virtual work-rate principle with an approximate interface constraint has been suggested to use a C ⁰-continuous element in the finite element implementation, and the couple stress work-rate has been formulated on the basis of an assumed constitutive behavior. Simulated plane-strain compressions of a single crystal cube show that the shearing and the deformation load are closely related to the imbedded lattice orientation of the crystal grain, and that the sub-grain deformation and the load magnitude can be controlled by the couple stress hardening. It is also confirmed that almost the same predictions are obtained for different mesh systems by considering the couple stress hardening. Simulated plane-strain compressions of a bi-crystal show considerably curved grain-by-grain surface profiles after large reduction for several combinations of the imbedded lattice orientation. The high couple stress hardening predicted around grain boundaries is supposed to be related to the grain size effect. It is also supposed that consideration of couple stress is necessary to predict the sub-grain or the grain-by-grain deformation, and the couple stress hardening may be used to describe the state of microstructures in grain.     

自增强厚壁圆筒疲劳裂纹扩展寿命的可靠性分析

本文提供了自增强厚壁圆筒疲劳裂纹扩展寿命的可靠性分析的方法和公式,自增强残余应力用符合厚壁筒用钢具有强化和包辛格效应实际性能推得的公式进行计算,残余应力对应的应力强度因子的计算采用了有限元法,得到了工作内压与自增强残余应力共同作用下的厚壁筒应力强度因子公式,通过实验测定了厚壁筒用钢的断裂韧度和疲劳裂纹扩展速率等性能。     

Effects of specimen thickness, hardening and crack closure for the plastic strip model

The triaxial stress constraint T_z and the effective yield stress distribution in the plastic zone for strain hardening materials are considered for the strip where the thickness effect is investigated by introducing a plastic constraint factor α. This factor depends on the specimen thickness, crack length, load level and hardening exponent. A simple expression of the plastic zone length and an expression involving α are obtained. Application of the strip model to crack closure shows that the specimen thickness has an influence and the results are compared with those found finite element.     

A WEIGHTED RESIDUAL METHOD FOR ELASTIC-PLASTIC ANALYSIS NEAR A CRACK TIP AND THE CALCULATION OF THE PLASTIC STRESS INTENSITY FACTORS

In this paper,a weighted residual method for the elastic-plastic analysis near a crack tip is systematically given by taking the model of power-law hardening under plane strain condition as a sample.The elastic-plastic solutions of the crack tip field and an approach based on the superposition of the nonlinear finite element method on the complete solution in the whole crack body field,to calculate the plastic stress intensity factors,are also developed.Therefore,a complete analysis based on the calculation both for the crack tip field and for the whole crack body field is provided.     

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.     

Finite element analysis of mode III steady crack growth in anisotropic hardening materials

In this paper, the steady crack growth of mode III under small scale yielding conditions is investigated for anisotropic hardening materials by the finite element method. The elastic-plastic stiffness matrix for anisotropic materials is given. The results show the significant influences of anisotropic hardening behaviour on the shape and size of plastic zone and deformation field near the crack tip. With a COD fracture criterion, the ratio of stress intensity factorsk _ss/k_c varies appreciably with the anisotropic hardening parameterM and the hardening exponentN.     

循环硬化材料本构模型的隐式应力积分和有限元实现

针对新发展的、能够描述循环硬化行为应变幅值依赖性的粘塑性本构模型，讨论了它的数值实现方法。首先，为了能够对材料的循环棘轮行为（Ratcheting）和循环应力松弛现象进行描述，对已有的本构模型进行了改进；然后，在改进模型的基础上，建立了一个新的、全隐式应力积分算法，进而推导了相应的一致切线刚度（Consistent Tangent Modulus）矩阵的表达式；最后，通过ABAQUS用户材料子程序UMAT将上述本构模型进行了有限元实现，并通过一些算例对一些构件的循环变形行为进行了有限元数值模拟，讨论了该类本构模型有限元实现的必要性和合理性。     

Effect of plastic zone on the fatigue crack propagation behavior between two fatigue cracks

In this study, the fatigue crack propagation behavior in the stress interaction field between two different fatigue cracks is studied by experiment and finite element analysis. In the experiment, the offset distance between two cracks and the applied stress are varied to create different stress interaction fields. The size of the plastic zone area is used to examine the crack propagation path and rate. Three types of crack propagation in the interaction field were found by experiment, and the crack propagation behavior of two cracks was significantly changed as different stresses were applied. The size of the plastic zone obtained by finite element analysis can be used to explain crack propagation behavior qualitatively.