Surface initiated crack growth simulation in moving lubricated contact

The paper describes a two-dimensional computational model for simulating surface initiated crack growth in the lubricated contact area that leads to surface pitting of mechanical components. The model assumes size and orientation of the initial crack which is subjected to contact loading conditions, accounting for the elasto-hydrodynamic-lubrication effects and tangential loading due to sliding. The influence of a lubricating fluid, driven into the crack by hydraulic mechanism, is also considered. The minimum strain energy density criterion is used to analyze crack propagation with the aim of the finite element analysis. The model is applied to a real pitting problem of a gear. The results for pit sizes correlate well with those observed in experimental testing.     

A cyclic steady-state method for fatigue crack propagation: Evaluation of plasticity-induced crack closure in 3D

The numerical study of plasticity-induced crack closure using the node-release technique presents many difficulties widely studied in literature. For instance various rules, proposed for overcoming mesh sensitivity, are challenged by more recent studies. This paper intends to propose and evaluate a numerical method for the investigation of crack propagation under fatigue loading, and particularly for the assessment of plasticity-induced crack closure in three-dimension. The method is an extension of the “steady-state method” to cyclic loadings. The steady-state method allows a direct computation (on a fixed mesh, without releasing nodes) of stress and strain fields around the crack tip and in the wake for a steady crack growth. The method is extended to simulate crack propagation under fatigue loading. Therefore it constitutes a valuable numerical tool for gaining insight into the physics of crack propagation, as it provides accurate mechanical fields around the crack tip and their relation with crack growth rate, various loading modes and parameters. The proposed method is also compared with the classical node-release technique. A very good agreement between the two methods is found. However the steady-state method needs much less mesh refinement and computational time. Following an analysis of some features of the fatigue crack, a discussion on a crack closure criterion is opened, and a reliable criterion for the determination of local crack closure is proposed.     

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

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

Elastic-viscoplastic field at mixed-mode interface crack-tip under compression and shear

For a compression-shear mixed mode interface crack, it is difficult to solve the stress and strain fields considering the material viscosity, the crack-tip singularity, the frictional effect, and the mixed loading level. In this paper, a mechanical model of the dynamic propagation interface crack for the compression-shear mixed mode is proposed using an elastic-viscoplastic constitutive model. The governing equations of propagation crack interface at the crack-tip are given. The numerical analysis is performed for the interface crack of the compression-shear mixed mode by introducing a displacement function and some boundary conditions. The distributed regularities of stress field of the interface crack-tip are discussed with several special parameters. The final results show that the viscosity effect and the frictional contact effect on the crack surface and the mixed-load parameter are important factors in studying the mixed mode interface crack- tip fields. These fields are controlled by the viscosity coefficient, the Mach number, and the singularity exponent.     

Analysis of crack moving and curving in anisotropic solids

The general equations for a dynamically curved crack in an anisotropic solid are derived, and the asymptotic fields of a moving crack under arbitrary distributed loading on the crack surface are calculated from them. For a moving crack under mixed-mode loading conditions a general Muskhelishvili type approach is proposed to calculate intensity factors due to crack surface loading in anisotropic materials. The kinking and curving caused by dynamic loading in anisotropic materials are calculated using the maximum normal stress ratio criterion. The results show that cracks in anisotropic solids may deviate from the straight path and approach a direction parallel to the stiff axis even under symmetric loading and that a crack will tend to deviate more from the crack path to the direction of the stiff axis as the crack speed becomes higher.     

Sufficient criterion of fracture in the case with a complex stress state and non-proportional deformation of the material in the pre-fracture zone

A general case of proportional loading with a complex stress state of the material in the pre-fracture zone, which is typical for polycrystalline solids with plastic deformation, is considered. A sufficient criterion of fracture is proposed for the case of a complex stress state with non-proportional deformation of the material in the pre-fracture zone. Critical parameters of fracture (pre-fracture zone length and load) for cracks propagating in quasi-brittle materials are obtained with the use of a modified Leonov-Panasyuk-Dugdale model. The pre-fracture zone width is determined by solving the problem of the plasticity theory in the vicinity of the crack tip. The proposed modification of the Leonov-Panasyuk-Dugdale model makes it possible to estimate the critical opening of the crack and the critical displacement of the crack flanks. Inequalities that describe different mechanisms of material fracture under proportional loading (predominantly shear fracture mechanism and fracture mechanism through cleavage) are derived.     

基于非齐次复合Poisson过程的正交异性钢桥面板细节疲劳裂纹随机扩展研究

传统的正交异性钢桥面板疲劳损伤评估常采用确定性和可靠性分析方法,忽略了疲劳裂纹扩展的随机性影响,针对这一问题,提出钢桥面板细节疲劳随机扩展分析方法。本文以南溪长江大桥为工程背景,基于长期车辆荷载监测数据,建立了车辆荷载非齐次复合Poisson过程模型。建立钢桥面板有限元模型,采用瞬态分析方法将随机车辆荷载转化成细节疲劳应力,基于线弹性断裂力学理论推导U肋-顶板焊接细节疲劳裂纹扩展时变微分方程,实现宏观关系式疲劳应力幅次数-疲劳损伤至微观表达式应力时间序列-疲劳损伤转换,讨论了车载次序及超载对疲劳裂纹扩展的影响。研究结果表明,非齐次复合泊松过程模型能够较好描述随机车流运营状态,车辆荷载的次序对疲劳裂纹扩展速率的影响不可忽略,重车排序靠前时能够促使疲劳裂纹扩展增速,南溪长江大桥细节点的车辆超载迟滞效应修正系数取值0.804。     

STUDY ON THE SURFACE CRACK GROWTH BEHAVIOR IN 14MnNbq BRIDGE STEEL

Three-dimensional crack closure correction methods are investigated in this paper.The fatigue crack growth tests of surface cracks in 14MnNbq steel for bridge plate subjected to tensile and bending loadings are systematically conducted.The experimentally measured fatigue crack growth rates of surface cracks are compared with those of through-thickness cracks in detail.It is found that the crack growth rates of surface cracks are lower than those of through-thickness cracks.In order to correct their differences in fatigue crack growth rates, a dimensionless crack closure correction model is proposed.Although this correction model is determined only by the experimental data of surface cracks under tensile loading with a constant ratio R=0.05, it can correlate the surface crack growth rates with reasonable accuracy under tensile and bending loadings with various stress ratios ranging from 0 to 0.5.Furthermore, predictions of fatigue life and crack aspect ratio for surface cracks are discussed, and the predicted results are also compared with those obtained from other prediction approaches.Comparison results show that the proposed crack closure correction model gives better prediction of fatigue life than other models.     

Bayesian Statistic Based Multivariate Gaussian Process Approach for Offline/Online Fatigue Crack Growth Prediction

Offline and online fatigue crack growth prediction of Aluminum 2024 compact-tension (CT) specimens under variable loading has been modeled, using multivariate Gaussian Process (GP) technique. The GP model is a Bayesian statistic stochastic model that projects the input space to an output space by probabilistically inferring the underlying nonlinear function. For the offline prediction, the input space of the model is trained with parameters that affect fatigue crack growth, such as the number of fatigue cycles, minimum load, maximum load, and load ratio. For the online prediction, the model input space is trained using piezoelectric sensor signal features rather than training the input space with loading parameters, which are difficult to measure in a real time scenario. Principal Component Analysis (PCA) is used to extract the principal features from sensor signals. In both the offline and online case, the output space is trained with known associated crack lengths or crack growth rates. Once the GP model is trained, a new output space for which the corresponding crack length or crack growth rate is not known, is predicted using the trained GP model. The models are validated through several numerical examples.     

Dynamic fracture of concrete – compact tension specimen

The behavior of concrete structures is strongly influenced by the loading rate. Compared to quasi-static loading concrete loaded by impact loading acts in a different way. First, there is a strain-rate influence on strength, stiffness, and ductility, and, second, there are inertia forces activated. Both influences are clearly demonstrated in experiments. Moreover, for concrete structures, which exhibit damage and fracture phenomena, the failure mode and cracking pattern depend on loading rate. In general, there is a tendency that with the increase of loading rate the failure mode changes from mode-I to mixed mode. Furthermore, theoretical and experimental investigations indicate that after the crack reaches critical speed of propagation there is crack branching. The present paper focuses on 3D finite-element study of the crack propagation of the concrete compact tension specimen. The rate sensitive microplane model is used as a constitutive law for concrete. The strain-rate influence is captured by the activation energy theory. Inertia forces are implicitly accounted for through dynamic finite element analysis. The results of the study show that the fracture of the specimen strongly depends on the loading rate. For relatively low loading rates there is a single crack due to the mode-I fracture. However, with the increase of loading rate crack branching is observed. Up to certain threshold (critical) loading rate the maximal crack velocity increases with increase of loading rate, however, for higher loading rates maximal velocity of the crack propagation becomes independent of the loading rate. The critical crack velocity at the onset of crack branching is found to be approximately 500 m/s.     

基于紧凑拉伸剪切结构的复合型疲劳裂纹扩展研究

针对含I-II型复合裂纹的紧凑拉伸剪切(CTS)试样,研究了不同加载角度下的裂纹扩展路径及裂纹扩展寿命,通过实验数据给出了适宜于CTS试样的等效应力强度因子关系式,并基于此提出了一种新的I-II型复合裂纹扩展模型。研究表明,CTS试样的裂纹沿与加载方向近垂直的方向扩展,基于Tanaka公式的等效应力强度因子更适合于本文CTS试件的裂纹扩展寿命评估。当加载角度处于0°~45°之间时,提出的复合型裂纹扩展模型预测误差控制在5.49%之内,验证了分析模型的可行性和准确性。     

循环加载条件下混凝土Ⅰ型裂缝扩展模拟的接触算法

沿裂缝可能开展路径设置接触点对,把求解接触问题的有限元混合法进行扩展,以实现循环加载条件下混凝土Ⅰ型裂缝扩展的数值模拟.根据荷载过程状态和缝面接触状态的不同,将循环荷载作用下的混凝土本构关系抽象为力学模型中的六种不同接触状态,不同的接触状态对应着不同的位移-应力曲线关系.以缝面张开位移和接触应力作为接触状态转变的判断参量,并给出了各状态的转变关系和数值判断条件,以接触算法实现了循环加载条件下混凝土Ⅰ型断裂扩展的数值模拟.首先给出了有限元混合法求解接触问题的基本思路,然后引入循环荷载下的混凝土本构关系,再对Ⅰ型裂缝扩展的数值实现方法进行了阐述,最后通过数值算例说明了数值实现方法的正确性和有效性.     

基于有限断裂法和比例边界有限元法的裂纹扩展模拟

基于有限断裂法和比例边界有限元法提出了一种裂缝开裂过程模拟的数值模型。采用基于有限断裂法的混合断裂准则作为起裂及扩展的判断标准,当最大环向应力和能量释放率同时达到其临界值时,裂缝扩展。结合多边形比例边界有限元法,可以半解析地求解裂尖区域附近的应力场和位移场,在裂尖附近无需富集即可获得高精度的解。计算能量释放率时,只需将裂尖多边形内的裂尖位置局部调整,无需改变整体网格的分布,网格重剖分的工作量降至最少。裂缝扩展步长通过混合断裂准则确定,避免了人为假设的随意性,并可以实现裂缝变步长扩展的模拟,更符合实际情况。通过对四点剪切梁的复合型裂缝扩展过程的模拟,对本文模型进行了验证,并应用于重力坝模型的裂缝扩展模拟,计算结果表明,本文提出的模型简单易行且精度较高。     

Prediction of fatigue crack growth life under variable-amplitude loading using finite element analysis

This article presents an elastic-plastic study aiming at predicting the fatigue crack growth (FCG) of 2024-T3 aluminum alloys under variable-amplitude loading. The proposed analysis needs the estimation of the residual stress distribution ahead of the crack tip during propagation. An elastic-plastic FE analysis has been implemented for modeling FCG using Chaboche's model. The FE study has been carried out through consideration of the loading history effect using the memory rules. Three different loading spectra have been applied in this work. The obtained results have been compared to the experimental ones and it has been proved that the suggested model has a better prediction of the FCG lives of cracked 2024-T3 aluminum alloy structures subjected to variable-amplitude loading.     

Frictional energy dissipation in materials containing cracks

Kachanov's simplified model of microcrack interaction is applied to an investigation of the behaviour of a cracked body under predominantly compressive periodic loading, so that the cracks experience periods of closure and slip, with frictional dissipation. The model is shown to be equivalent to a discrete elastic frictional system with each crack representing one node. Theorems and algorithms from such systems are applied to determine the conditions under which the system shakes down to a state with no slip and hence no energy dissipation in friction. For conditions not too far beyond the shakedown state, the dissipation is significantly affected by the initial conditions, but with larger oscillating loads, it becomes a unique and increasing function of load amplitude. The effect of crack interaction is assessed by comparison with an uncoupled model, for which the dissipation is obtained as a summation of closed form expressions over the crack population. For small numbers of cracks, the results are significantly dependent on the randomly chosen crack locations and sizes, but with larger populations, a statistically significant decrease in dissipation is observed with increasing interaction terms.     

基于能量耗散分析的混凝土随机疲劳损伤模型

基于微-细观随机断裂模型(MMSF)发展了一类混凝土随机疲劳损伤模型。该模型将MMSF中的微弹簧视为一能量耗散单元并考察了其跨尺度的能量耗散过程。在纳观尺度,引入速率过程理论描述裂纹的扩展速度,并基于裂纹层级模型和自相似假定完成从纳观尺度到微观尺度的过渡。由此建立了微弹簧在疲劳荷载作用下的多尺度耗能描述。此外,为了考虑疲劳加载中混凝土多条裂纹的相互影响,引入了包含损伤扩展效应和损伤愈合效应的疲劳损伤因子来修正耗能表达。通过与相关试验的对比,证明了该模型能反映疲劳荷载作用下混凝土的主要力学行为,如疲劳损伤三阶段特点、疲劳寿命的离散性和疲劳寿命随加载频率的变化趋势等。     

Effect of pure mode I,II or III loading or mode mixity on crack growth in a homogeneous solid

Crack growth is analyzed numerically under combined mode I, II and III loading, or under loading in one of these modes alone. The solid is a ductile metal modelled as elastic–plastic, and the fracture process is represented in terms of a cohesive zone model. The analyses are carried out for conditions of small-scale yielding, with the elastic mixed mode solution applied as boundary conditions on the outer edge of the region analyzed. For pure mode I loading crack growth continued far beyond the maximum fracture toughness shows that the predicted subsequent steady-state toughness is well below the maximum. The reason for this is discussed in terms of the local stress and strain fields around the tip. For pure mode II or mode III loading it is shown that there is no maximum before the steady-state. Also results for different mixed mode conditions are presented and discussed in relation to the results for loading in only one mode. Most of the results are based on assuming that the peak tractions for tangential separation are equal to that for normal separation, but it is shown that a relatively smaller peak traction for tangential separation may significantly affect the predictions.     

Repeated impact on fatigue fracture of carbon spring En-42J steel

The influence of repeated impact loading on fatigue fracture of carbon spring En-42J steel is studied. A highly localized plastic deformation is observed at the time of impact. The circular disc specimen is failed by impact fatigue resulting in small chips and ordinary fatigue failure by slow crack growth. The crack growth behavior depends on specimen geometry, load spectra and material type. The resultant of applied and residual stress was highest at the edge of the circular disc where the crack initiated. The fractured surface showed typical cleavage facets and dimples indicating acicular martensite. A model is developed to explain fatigue crack initiation and propagation behavior under impact loading. Stress calculations are made to determine catastrophic failure. Crack initiation occupied the major part of impact fatigue life while crack growth was fast and led quickly to sudden fracture.     

Damage accumulation for growth of anti-plane crack in work-hardening material

Elastic–plastic solutions of an anti-plane crack in an infinite body are used in conjunction with a continuum damage model to describe the conditions necessary for the onset of crack instability, fatigue crack propagation due to cyclic loading, and rates of crack growth due to time dependent events. A power law relates the stress to the strain of the material. The damage, which invokes nucleation, growth and coalescence of microvoids due to elevated strain, is confined to the plastic zone surrounding the crack tip. For applied loading below the yield stress, the small-scale and large-scale yielding solutions are used to determine the influence of strain hardening on crack instability and failure. Crack growth due to cyclic loading and time-dependent deformations are studied using the small-scale yielding solution of the deformation theory of plasticity.     

短裂纹闭合的尺寸效应分析

对循环荷作用下缺口根部静态短裂纹和扩展短裂纹的闭合现象进行了弹塑性有限元分析。考虑不同应力水平和应力比，尤其是零－压载荷对短裂纹的作用。在此基础上提出分析长，短裂纹闭合尺寸效应的修正模型。与实验结果相比，本文模型对短裂纹的闭合预测甚好。