A cohesive plastic and damage zone model for dynamic crack growth in rate-dependent materials

Mode I steady-state dynamic crack growth in rate-dependent viscoplastic solids containing damage, under small scale yielding conditions, is analyzed based on a modified cohesive zone model. A multi-scale approach is used to describe the entire non-linear zone consisting of a plastic region and a damage region, each of which has its own constitutive law. Traction in the damage region is characterized by a softening power-law, in terms of the ultimate strength, a softening index and a rate sensitivity factor. In the plastic region, the cohesive law is assumed to be both strain hardening and rate dependent. The critical crack opening displacement at the physical crack-tip controls crack growth. The governing integral equations are derived and solved by a collocation method combined with associated boundary conditions. Numerical results are presented for the traction and opening profiles along the cohesive zone, the fracture energy and lengths of the damage and non-linear zones at different crack speeds and for different material parameters. The importance of factors, such as material softening, plastic deformation, crack speed and viscosity, is identified by parametric studies. In addition, the competition of plastic flow and material damage, and its effect on crack growth, are discussed.     

Predicting ductile fracture of low carbon steel sheets: Stress-based versus mixed stress/strain-based Mohr–Coulomb model

Two distinct implementations of the Mohr–Coulomb failure model are used in conjunction with a non-associated quadratic plasticity model to describe the onset of fracture in low carbon steel sheets. The stress-based version corresponds to the original Mohr–Coulomb model in stress space. For the mixed stress/strain-based version, the Mohr–Coulomb failure criterion is first transformed into the space of stress triaxiality, Lode angle parameter and equivalent plastic strain and then used as stress-state dependent weighting function in a damage indicator model. Basic fracture experiments including tensile specimens of different notch radii and a punch test are performed to calibrate the material parameters of the respective models. Subsequently, the models are used to predict the crack initiation in a Hasek test and during the stamping of an anticlastic structure. Unlike for the calibration experiments, the loading history during stamping is highly non-linear. Both models can be calibrated with similar accuracy, but the strain-based model predicts the instant of onset of fracture with greater accuracy in the stamping experiment which is an advantage of the empirical damage accumulation rule.     

无铆连接模具参数对异种金属板料接头成形的影响

为了研究无铆连接模具参数对 304 不锈钢和 AL6061 板料接头成形的影响，基于 DEFORM-2D 建立了有限元模型，分析了凹模深度、凹槽深度和凸模直径三个主要模具参数对接头成形的影响．研究结果表明，凹模深度越大，越有利于下板料变形，反之，则有利于上板料变形．当凹模直径一定时，凹槽深度和凸模直径过小不利于形成互锁，而凹槽深度不宜过大，避免出现充不满现象，凸模直径过大导致接头颈厚值过小．通过正交试验发现，对接头成形的影响程度从大到小依次为凸模直径、凹槽深度和凹模深度，确定了最优模具参数，并通过无铆连接实验验证了正交试验的可靠性．此外，通过剪切实验获得了接头最大剪切失效载荷为 1.8 kN．     

准脆性材料断裂过程区中的圆饼微裂纹损伤计算

准脆性工程材料及结构在外力作用下，不仅引起内部缺陷变化和微裂纹的出现及发展，且使得其结构承载能力降低或性能劣化．在其材料失效过程中常存在裂缝与断裂损伤过程区．为研究材料细观缺陷或微裂纹与宏观破坏的规律，通过细观力学方法，对于代表性体积单元RVE中的圆饼型微裂纹的尺寸与密度变化，探讨其宏观断裂过程区力学参量与损伤之间的量化关系．借助宏观断裂过程区的黏聚裂纹模型，将损伤单元RVE嵌入到宏观裂缝端部的断裂过程区中，对其进行联接细观损伤到宏观破坏的力学多尺度研究．文中也通过实验数据，对其理论计算结果进行了算例的讨论与分析．     

Fatigue and fracture of bovine dentin

In this paper, the fatigue and fracture properties of bovine dentin are evaluated usingin vitro experimental analyses. Double cantilever beam (DCB) specimens were prepared from bovine maxillary molars and subjected to zeroto-tension cyclic loads. The fatigue crack growth rate was evaluated as a function of the dentin tubule orientation using the Paris law. Wedge-loaded DCB specimens were also prepared and subjected to monotonic opening loads. Moiré interferometry was used to acquire the in-plane displacement field during stable crack growth, and the instantaneous wedge load and crack length were acquired to evaluate the crack growth resistance and crack tip opening displacement (CTOD) with crack extension. The rate of fatigue crack growth was generally larger for crack propagation occurring perpendicular to the dentin tubules. The Moiré fringe fields documented during monotonic crack growth exhibited non-linear deformation occurring within a confined region adjacent to the crack tip. Both the wedge load and CTOD response provided evidence that a fracture process zone contributes to energy dissipation during crack extension and that dentin exhibits a risingR-curve behavior. Results from this preliminary investigation are being used as a guide for an evaluation of the fatigue and fracture properties of human dentin.     

Crack initiation at contact surface

By using the S-theory, the crack initiation angle from the contact surface of rectangular rigid punch and flat-surface substrate has been investigated. The coefficient of friction at the edge of contact, which characterizes the asymptotic stress field, is considered as the controlling parameter in the analysis. The predicted results are in agreement with the experimental observations. The information gained may lend insight into the different stages of damage associated with the complex process of fretting fatigue.     

Experimental determination of damage parameters in uniaxially-loaded metal-matrix composites using the overall approach

A recently developed damage-plasticity model for metal matrix composites is examined using the overall approach to damage. In this approach, the damage mechanisms in the composite system are reflected through a single overall damage tensor assuming an elasto-plastic matrix with elastic fibers. It is shown how the model is applied to uniaxially loaded metal-matrix specimens. Each specimen is a symmetric laminate composed of four plies with different orientations. A numerical procedure is formulated and implemented to calculate the stresses, strains, and damage variables in these specimens.

In order to test the validity of the model, a series of uniaxial test experiments are conducted for the damage characterization of the metal-matrix composite. The experiments are conducted on a titanium aluminide SiC-reinforced laminate. Laminate layups of (0/90)s and (± 45)s are tested under uniaxial tension. In order to investigate and model damage evolution, each specimen is loaded to different load levels ranging from the fracture load to 70% of fracture load. Scanning electron microscopy is performed on representative cross-sections of all specimens. This information is then used to quantify damage in the composite system. Damage variables are defined based on the crack density and correlated with the theoretical model. The experimental and numerical results are then compared. Certain recommendations are finally made concerning the determination of damage parameters in metal-matrix composites.     

A rate-dependent damage model for brittle materials based on the dominant crack

A rate-dependent, continuum damage model is developed for brittle materials under dynamic loading. This model improves on the approach (ISOSCM) of [Addessio, F.L., Johnson, J.N., 1990. A constitutive model for the dynamic response of brittle materials. Journal of Applied Physics 67, 3275–3286] in several respects. (1) A new damage surface is found by applying the generalized Griffith instability criterion to the dominant crack (having the most unstable orientation), rather than by averaging the instability condition over all crack orientations as done previously. The new surface removes a discontinuity in the damage surface in ISOSCM when the pressure changes sign. (2) The strain due to crack opening is more consistent with crack mechanics, with only the tensile principal stresses contributing to the crack opening strain. This is achieved by incorporating a projection operator in the equation for the crack opening strain. One consequence of incorporating the projection operator is a prediction of shear dilatancy, which is not accounted for in ISOSCM. (3) The evolution of damage, which is based on the energy-release rate for the dominant crack, has a physical basis, whereas in the previous approach the damage growth rate was assumed to be an exponential function of the distance from the stress state to the damage surface without specific physical justification.An implicit algorithm has been developed so that a larger time step can be used than with the explicit algorithm used in ISOSCM. The numerical results of a silicon carbide (SiC) ceramic under several loading paths (hydrostatic tension/compression, uniaxial strain, uniaxial stress, and shear) and strain rates are presented to illustrate the main features of the model.     

微小缺口/孔洞的激光衍射无损检测

基于裂纹和孔洞的小尺寸特征,提出了微小缺口/孔洞的激光衍射无损探测技术,给出了解析表达式.通过这一技术,对单向拉伸试件中所含单边缺口和中心孔洞在外载作用下的演化过程进行了实时原位检测,获得了缺口/孔洞孔径随载荷的变化曲线及模拟裂纹时裂纹的张开位移、裂纹开裂长度及应力强度因子等一系列断裂参数.     

Forming limit diagram prediction of tailor welded blank by modified M–K model

The Marciniak and Kuczynski (M–K) model for necking prediction in sheet metal forming was based on the in-plane forming. Bending which was resulted from out-of-plane forming was not considered in the M–K model. Whereas most of the sheet metal forming processes and also standard test of hemispherical punch for forming limit diagram are out-of-plane forming, it is important to consider bending effect in the M–K model. Therefore, in this study bending strain is added to stretching strain of M–K model and a new model is presented for forming limit diagram (FLD) prediction. This modified M–K (MM–K) model is written in the python programming language and it is used as a post-processing criterion for FLD prediction in the commercial software Abaqus. The MM–K model was used to predict FLD and weld line movement in the tailor welded blank forming. It was found that the predicted results by MM–K model are in a good agreement with experimental data.     

Initiation and stable growth of penny-shaped crack in aging viscoelastic transversally isotropic material

Considered is the long-term cracking of an aging transversally isotropic material containing a Mode I penny-shaped crack under remotely applied tensile stress. The aging material properties are described by the Boltzmann–Volterra’s linear theory for integral operators with non-difference kernels. It applied to wood, concrete, some polymers and rocks. Only the symmetric case is considered where the crack lies in the plane of isotropy. The modified Leonov–Panasyuk–Dugdale’s crack model is used with a constant process zone assuming that the critical opening displacement is the fracture criterion. Volterra’s principle is applied to derive the equations of subcritical crack growth. Numerical calculations are made for subcritical crack growth for the specific example of transversally isotropic material simulating the behavior of reinforced concrete.     

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

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

Dynamic fracture behaviour in fibre-reinforced cementitious composites

The object of this work is to simulate the dynamic fracture propagation in fibre-reinforced cementitious composites, in particular, in steel fibre reinforced concrete (SFRC). Beams loaded in a three-point bend configuration through a drop-weight impact device are considered. A single cohesive crack is assumed to propagate at the middle section; the opening of this crack is governed by a rate-dependent cohesive law; the fibres around the fracture plane are explicitly represented through truss elements. The fibre pull-out behaviour is depicted by an equivalent constitutive law, which is obtained from an analytical load–slip curve. The obtained load–displacement curves and crack propagation velocities are compared with their experimental counterparts. The good agreement with experimental data testifies to the feasibility of the proposed methodology and paves the way to its application in a multi-scale framework.     

A cohesive zone model for fatigue crack growth allowing for crack retardation

A damage-based cohesive model is developed for simulating crack growth due to fatigue loading. The cohesive model follows a linear damage-dependent traction–separation relation coupled with a damage evolution equation. The rate of damage evolution is characterized by three material parameters corresponding to common features of fatigue behavior captured by the model, namely, damage accumulation, crack retardation and stress threshold. Good agreement is obtained between finite element solutions using the model and fatigue test results for an aluminum alloy under different load ratios and for the overload effect on ductile 316 L steel.     

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.     

Application of microline/grid methods to temperature dependent microscopic deformation and damage in CFRP laminates

Tensile tests of CFRP symmetric cross-ply laminates are carried out in a scanning electron microscope (SEM), and microscopic interlaminar deformation and damage near the transverse crack tip are visualized by microlines or microgrids printed on the specimen edge surface. The local deformation around the transverse crack tip is observed at 20°C, 80°C, 120°C and 160°C to evaluate the effect of thermal residual stress on the microscopic deformation and damage in the interlaminar region near the transverse crack tip. Temperature dependence of the axial crack opening displacement (COD) is also measured. The displacement field of the specimen edge surface obtained from these experimental results is compared with the theoretical model proposed by Lee, Allen and Harris. The analysis is modified to consider the temperature effect, and a good agreement is obtained between the modified theoretical predictions and experimental results.     

Modeling hydrogen attack effect on creep fracture toughness

The effect of high temperature hydrogen attack on creep crack growth rates in steels is studied by modeling the interaction between creep deformation and gaseous pressures generated by hydrogen and methane. The equilibrium methane pressure as a function of hydrogen pressure, temperature and carbide types for carbon steels and Cr–Mo steels is calculated. This gaseous driving force is incorporated into a micromechanics model for void growth along grain boundaries of a creeping solid. Growth and coalescence of voids along grain boundaries is modeled by a microporous strip of cell elements, referred to as the fracture process zone. The cell elements are governed by a nonlinear viscous constitutive relation for a voided material. Two rate sensitivities as well as two types of grain boundaries are considered in this computational study. Simulations of creep crack growth accelerated by gaseous pressures are performed under conditions of small-scale and extensive creep. The computed crack growth rates at elevated temperatures are able to reproduce the trends of experimental results.     

刻槽式MEFP的成形及侵彻钢靶模式

多爆炸成形弹丸(MEFP)能有效提高炸药利用率和命中概率。基于LS-DYNA软件对刻槽式MEFP战斗部成形及侵彻双层无间隔钢靶模式进行了研究,得到了刻槽式MEFP战斗部成形及对靶板侵彻的数值计算结果。结果表明:刻槽式MEFP成形过程要经过翻转、头部挤压断裂、尾部拉伸断裂、交叉飞散、稳定飞行5个阶段, 聚能爆轰波对药型罩的径向力是其交叉飞散的主要原因;侵彻过程要经过开坑、联合侵彻、贯穿3个阶段。模拟结果与实验结果基本吻合。