半无穷大裂纹端部粘聚力分析

准脆性材料裂纹端部断裂过程区粘聚力是导致非线性断裂特性的重要原因,根据准脆性材料的断裂特性,对存在粘聚力分布的半无穷大裂纹力学分析模型,由变形叠加原理得到以该粘聚应力分布为未知函数的积分方程,通过对积分方程的分析推证,得到了该分布函数解的数学结构和级数型表达式;提出了由实际裂纹张开位移,确定裂纹端部粘聚力分布函数的两种方法:其一由连续的裂纹张开位移通过积分变换求解未知函数级数展开项的系数,其二是由离散的裂纹张开位移数据通过最小二乘法确定该函数;推导出了相应方法求解未知量的代数方程,并且给出了适当的算例和讨论。     

幂律非线性粘弹性材料中的裂纹扩展*

对蠕变不可压幂律非线性粘弹性材料中裂纹的蠕变扩展进行了分析,为描述银纹带中的力学行为,假设在裂纹尖端邻域中断裂过程区中分布着阻抗裂纹张开的粘聚应力б_f,.通过对均匀应力参考状态平凡解的摄动,将非线性粘弹性问题化成线性问题处理.对于幂指数.n≌1的弱非线性情况得到了应力与位移表达式.提出断裂过程区局域能量判据,导出了裂纹孕育时间t*与蠕变扩展率a的预测公式.     

Subcritical and critical states of a crack with failure zones

The problem on the stress–strain state near a mode I crack in an infinite plate is solved in the frame of a cohesive zone model. The complex variable method of Muskhelishvili is used to obtain the crack opening displacements caused by the cohesive traction, which models the failure zone at the crack tip, as well as by the external load. The finite stress condition and logarithmic singularity of the derivative of the separation with respect to the coordinate at the tip of a physical crack are taken into account.The cohesive traction distribution is sought in a piecewise linear form, nodal values of which are being numerically chosen to satisfy the traction-separation law. According to this law, the cohesive traction is coupled with the corresponding separation and fracture toughness. The tips of the physical crack and cohesive zone (geometric variables) along with the discrete cohesive traction are used as the problem parameters determining the stress-strain state. If the crack length is included in the set, then the critical crack size can be found for the given loading intensity.The obtained determining system of equations is solved numerically. To find the initial point for a standard numerical algorithm, the asymptotic determining system is derived. In this system, the geometric variables can be easily eliminated, which make it possible to linearize the system.In the numerical examples, the one-parameter traction-separation laws are used. Influence of the shape parameters of the law on the critical crack size and the corresponding cohesive length is studied. The possibility of using asymptotic solutions for determining the critical parameters is analysed. It is established that the critical crack length slightly depends on the shape parameter, while the cohesive length shows a strong dependence on the shape of cohesive laws.     

Simulation of crack growth using cohesive crack method

We present an effective cohesive discrete crack method in the context of the Reproducing Kernel Particle Method (RKPM) in order to study fracture of concrete structures. The discrete crack approach is based on the visibility method and a simple node splitting scheme. We also present an effective implementation of the visibility method and an iteration free algorithm by including the cohesive force term directly into the stiffness equations. The crack is represented by straight-line segments and the cohesive zone model is employed to model the post-localization behavior of concrete. The method is applied to several examples involving mode I and mixed-mode fracture. These results are compared to experimental data and show good agreement.     

Interlaminar Crack Propagation Analysis of ENF Specimens Based on the Cohesive Zone Model北大核心CSCD

Based on the classical laminated plate theory and the cohesive zone model, a theoretical model for general delamination cracked laminates was established for crack propagation of pure mode Ⅱ ENF specimens. Compared with the conventional beam theory, the proposed model fully considered the softening process of the cohesive zone and introduced the nonlinear behavior of ENF specimens before failure. The predicted failure load is smaller than that under the beam theory and closer to the experimental data in literatures. Compared with the beam theory with only fracture toughness considered, the proposed model can simultaneously analyze the influences of the interface strength, the fracture toughness and the initial interface stiffness on the load-displacement curves in ENF tests. The results show that, the interface strength mainly affects the mechanical behavior of specimens before failure, but has no influence on crack propagation. The fracture toughness is the main parameter affecting crack propagation, and the initial interface stiffness only affects the linear elastic loading stage. The cohesive zone length increases with the fracture toughness and decreases with the interface strength. The effect of the interface strength on the cohesive zone length is more obvious than that of the fracture toughness. When the adhesive zone tip reaches the half length of the specimen, the adhesive zone length will decrease to a certain extent. Copyright ©2022 Applied Mathematics and Mechanics. All rights reserved.     

Implementation of a cohesive zone model for the investigation of the dynamic behavior of a rotating shaft with a transverse crack

The influence of a transverse crack on the vibration of a rotating shaft has been at the focus of attention of many researchers. The knowledge of the dynamic behavior of cracked shaft has helped in predicting the presence of a crack in a rotor. Here, the changing stiffness of the cracked shaft is investigated based on a cohesive zone model. This model is developed for mode-I plane strain and accounts for triaxiality of the stress state explicitly by using basic elastic-plastic constitutive relations. Then, the proposed numerical solution is compared to the switching crack model, which is based on linear elastic fracture mechanics. The cohesive zone model is implemented in finite element techniques to predict and to analyse the dynamic behavior of cracked rotor system. Timoshenko beam theory is used to model the discrete shaft under the effect of gravity, unbalance force and gyroscopic effect. The analysis includes the cohesive function for describing the breathing crack and the reduction of the second moment of area of the element at the location of the crack. (© 2010 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

求解层间界面反平面剪切破坏的剪切梁模型(Ⅰ)——基本特性

在反平面剪切载荷及侧压力共同作用下引起的裂纹及裂纹扩展导致的层间界面失效,是岩土工程层间界面及砌体结构中界面层上典型的失效方式.运用弹性力学和断裂力学的理论原理,提出了能够反映上述层间界面断裂失效问题力学特性的剪切梁模型.文中采用具有应力软化特性的“粘性裂纹”(内聚力裂纹)模型来表述层间裂纹前方损伤过程区的本构行为.对通过粘性层结合在一起的两个弹性板,在反平面剪切载荷及侧压力共同作用下的力学行为作了解析分析计算,研究了层间界面裂纹扩展规律.     

Implementation of a cohesive zone model for crack closure analysis of rotors

The presence of a crack in a rotor introduces a local flexibility which affects its dynamic response. Moreover, the crack may open and close during the vibration period. The crack status is a function of time and also depends on the rotational speed and the vibration amplitude of the rotor. This nonlinear case is still a challenging research topic especially in the field of closing crack in the rotating shaft. A cohesive zone model is developed in order to analyze the stiffness of a crack in a rotating shaft. The proposed expression will be compared to three different crack models, namely, a breathing crack model, a switching crack model and an open crack model. Moreover, a cohesive law to predict and to analyse the stress at the crack tip is presented. The numerical model is implemented using a finite element formulation. (© 2009 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Assessment of asymmetric crack initiation in open-hole plates using a coupled stress and energy criterion

An application of the finite fracture mechanics concept to open-hole plates subject to combined tensile and bending loading is presented. In finite fracture mechanics, the simultaneous satisfaction of both, a stress and an energy criterion, is enforced as a condition for crack initiation. Efficient modeling and closed-form expressions for the dependence of the stress and energy quantities on governing structural and material parameters allow for a comprehensive numerical analysis of the onset of asymmetric crack patterns. The obtained failure load predictions are found to agree well with a cohesive zone model and experimental data from literature. (© 2016 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Crack opening displacement for two unequal straight cracks with coalesced plastic zones – A modified Dugdale model

The problem investigated is of an infinite plate weakened by two collinear unequal hairline straight quasi-static cracks. Uniform constant tension is applied at infinity in a direction perpendicular to the rims of the cracks. Consequently the rims of the cracks open in Mode I type deformation. The tension at infinity is increased to the limit such that the plastic zones developed at the two adjacent interior tips of cracks get coalesced. To arrest the crack from further opening normal cohesive variable stress distribution is applied on the rims of the plastic zones. Closed form analytic expressions are obtained for load bearing capacity and crack opening displacement (COD). An illustrative case is discussed to study the behavior of load bearing capacity and crack opening displacement with respect to affecting parameters viz. crack length, plastic zone length and inter crack distance between the two cracks. Results obtained are reported graphically and analyzed.     

A coupled stress and energy failure model for crack initiation in arbitrarily shaped adhesive lap joints

In this work, crack initiation in adhesive lap joints of arbitrary joint configuration is studied by means of a finite fracture mechanics approach. The analysis is based on a general stress solution for adhesive joints combined with a coupled stress and energy criterion. The instantaneous formation of a crack of finite size is predicted if a stress and energy criterion are satisfied simultaneously. The closed-form analytical solution of the stress field allows for an efficient evaluation of the crack initiation load and corresponding finite crack length. A comparison to experimental results from literature and to numerical results obtained with a cohesive zone model approach shows a good agreement. (© 2016 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Material state change relationships to fracture path development for large-strain fatigue of composite materials

The long-term performance of engineering structures is typically discussed in terms of such concepts as structural integrity, durability, damage tolerance, fracture toughness, etc. These familiar concepts are usually addressed by considering balance equations, crack growth relationships, constitutive equations with constant material properties, and constant or cyclically applied load conditions. The loading histories are represented by changing stress (or strain) states only. For many situations, especially for those associated with high-performance engineering structures, the local state of the material may also change during service, so that the properties used in the equations are functions of time and history of applied conditions. For example, the local values of stiffness, strength, and conductivity are altered by material degradation to create "property fields" that replace the global constants, and introduce time and history into the governing equations. The present paper will examine a small set of such problems, which involve the accumulation of distributed damage and the development of an eventual fracture path leading to failure. Specifically, the paper discusses this problem in the context of material state changes measured by impedance variations as a method of following the details of fracture path development. An analysis and interpretations of observations will be presented, and limitations and opportunities associated with this general concept will be discussed.     

On the computation of material forces in continuum‐atomistic modelling

Material forces govern the behavior and evolution of defects in solids. In elastic materials these forces which are associated with the Eshelby stress tensor are used to describe fracture sensitivities and can be employed to compute the J‐integral [2]. Since crack propagation begins with a variety of fundamental processes which occur within highly localized ultra–fine volume of material that constitute the fracture process zone surrounding a crack tip [3], the question of appropriate growth criteria, i.e. how far and in which direction a crack will glide under a certain loading condition is implied by the material force. (© 2004 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Cohesive zone model for a transverse breathing crack in a rotor

The breathing mechanism of a transversely cracked rotor and its influence on a rotor system that appears due to the shaft weight is studied. This breathing mechanism is based on experimental and simulation result for the crack shape reported in the literature. If the crack depth is small, the crack closure line is a straight line while for larger crack depths the crack closure becomes more curved. For both cases, a method is proposed for the evaluation of the stiffness losses in the cross section that contains the crack. This method is based on a cohesive zone model (CZM) instead of linear elastic fracture mechanics (LEFM) approach, because LEFM is valid only for the fully open crack and cannot be extended to other intermediate situations. As the crack is closed, the stress intensity factor (SIF) will not appear at the boundary between the closed cracked areas and the open cracked areas. The CZM is developed for mode-I plane strain conditions and accounts explicitly for triaxiality of the stress state by using constitutive relations. The proposed model gives more realistic results than models based on LEFM for the stiffness losses of the crack rotor system for a wide range of the crack depth. (© 2011 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

基于裂纹尖端二阶弹性解的断裂过程区尺度

基于Westergaard应力函数裂纹尖端二阶弹性解,推导了裂纹尖端微裂区的轮廓线和特征尺寸的解析表达式;采用幂函数模型描述的拉应变软化模型,确定了在最大拉应力强度理论和最大拉应变强度理论下断裂过程区（FPZ）临界值的解析表达式;将基于Westergaard应力函数一阶弹性解及二阶弹性解、Muskhelishvili应力函数和Duan-Nakagawa模型确定的FPZ临界值进行了比较．结果表明裂纹尖端微裂区和FPZ临界值随着Poisson比的减小而增加并逐渐趋近于应用最大拉应力强度理论确定的结果;二阶弹性解确定的裂纹尖端微裂区和FPZ临界值大于一阶弹性解的值;FPZ临界值随着拉应变软化指数的增加而增加;二阶弹性解确定的FPZ临界值的精度远高于一阶弹性解确定的值.     

Simulating quasi-brittle failures including damage-induced softening based on the mechanism of stress redistribution

There have been extensive analytical and numerical studies on crack propagation and softening behavior in quasi-brittle materials, such as concrete, rock and bone, based on the smeared crack model or the discrete crack model, each of which has advantages and disadvantages. It is well known that stress redistribution has an important effect on cracking behavior and may be taken as the physics mechanism. Using the principle of superposition, the theoretical explanation is provided for stress redistribution during the progressive fracture process. Although this explanation deals with problems of linear elastic fracture mechanics, it is a clue of modeling softening behavior. Then, it is shown that the above stress redistribution based method satisfies the necessary condition that the specimen returns to the unstressed state whenever all the external forces are removed. Finally, a new model for simulating quasi-brittle softening is established via a kind of "sub-element" methodology. Numerical examples show the effectiveness of the above subelement method in modeling the behavior of block media with quasi-brittle block-block interfaces.     

功能梯度条硬币型裂纹扭转冲击响应

研究非均匀条中硬币型裂纹的扭转冲击问题．材料的剪切模量假定按特定的梯度变化．采用Laplace 和Hankel 变换将问题化为求解Fredholm积分方程,通过将Bessel函数渐进展开获得裂纹尖端动态应力场．考查非均匀参数和功能梯度条高度对裂尖动态断裂行为的影响．动应力强度因子和能量密度因子的清晰表达式表明,作为裂纹扩展力,对于这里所研究的问题,二者是等价的．动应力强度因子的数值结果显示,增加剪切模量的非均匀参数可以抑制动应力强度因子的幅度,而条形域的高度对动态断裂特性的影响较小．     

Inplane Stress Singularities at the Corner of an Anisotropic Material Discontinuity

In the mechanics of composite laminates the local mechanical inplane fields at corners of anisotropic material discontinuities are of particular interest since they can have singular behavior. In the present study, the stress and strain fields in the local near field of such corners are investigated by an asymptotic analysis. The order of the singularity of these mechanical inplane fields are determined in closed‐form manner by use of the complex potential method based on Lekhnitskii's approach. Various different geometrical setups and material combinations of corners with material discontinuities are investigated with regard to their effect on the singular behavior of the mechanical fields present. These examples show that the order of singularity considered is clearly weaker than the typical crack tip singularity in fracture mechanics. Nevertheless, it may render the corner a critical location for the onset of failure. (© 2004 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

利用内聚力模型(CZM)模拟弹粘塑性多晶体的裂纹扩展

采用内聚力模型（CZM）,模拟多晶体中起裂于晶界的二维平面应变裂纹扩展．结果表明,弹粘塑性体中,初始裂纹尖端不会最先开裂．晶体本构的率敏感指数表征了塑性变形和内聚力区耗散两种机制的相互竞争．率敏感指数越大,塑性耗散能越大,内聚力区粘着能越小,使材料的塑性变形越容易,内聚力区诱发的破坏越不易;率敏感指数越小,材料响应越接近弹塑性性质,塑性耗散能减小,粘着能增大,外力功易转化为内聚力区的粘着能,使内聚力单元更易分离．增大内聚力区结合强度或临界张开位移使晶内和晶界的三轴应力度减小,即提高内聚力区韧性也使基体材料抗孔洞损伤能力提高．     

The effect of microdamage on fatigue-crack growth

The theory of fatigue-crack growth, based on a synthesis of fracture mechanics and continuum mechanics of microdamage accumulation, is applied to the problem of crack growth under cyclic loading, taking into account the plastic deformations in the tip zones. The model of a thin plastic zone, which is a region of considerable plastic deformations at the crack tip, is supplemented by taking into account the effect of microdamage on the value of the specific fracture work and the limit stresses in the tip zone. Governing equations which describe fatigue-crack growth taking these factors into account are derived. The effect of the material characteristics and the load parameters on the growth rate and the distribution of microdamage in the tip zone and on its extensions is investigated by a computational experiment. Particular attention is given to the initial stage when crack growth may occur abruptly and the growth rate depends substantially on the initial conditions