基于XFEM与自适应网格的非均质材料裂纹扩展模拟

针对含有间断的非均匀材料的断裂问题，本文将虚节点多边形单元的形函数引入到扩展有限元（XFEM）中，提出了一种基于四叉树结构的动态网格细化方法，该方法可对间断面附近单元实现可调控的多层级细化，特别是对于裂纹扩展问题，可实现裂尖附近单元的动态网格细化与粗化。基于以上网格细化方法，本文提出了针对非均匀材质裂纹扩展问题的计算方法VP-XFEM。为验证算法的准确性与计算效率，针对含有孔洞及材料界面的断裂问题，本文给出了相应的算例。结果显示，与传统的一致性网格的XFEM相比，VP-XFEM能够明显改善计算精度与计算效率。     

Mechanics of crack propagation in materials with initial (residual) stresses (review)

Major results on the mechanics of crack propagation in materials with initial (residual) stresses are analyzed. The case of straight cracks of constant width that propagate at a constant speed in a material with initial (residual) stresses acting along the cracks is examined. The results were obtained, based on linearized solid mechanics, in a universal form for isotropic and orthotropic, compressible and incompressible elastic materials with an arbitrary elastic potential in the cases of finite (large) and small initial strains. The stresses and displacements in the linearized theory are expressed in terms of analytical functions of complex variables when solving dynamic plane and antiplane problems. These complex variables depend on the crack propagation rate and the material properties. The exact solutions analyzed were obtained for growing (mode I, II, III) cracks and the case of wedging by using methods of complex variable theory, such as Riemann–Hilbert problem methods and the Keldysh–Sedov formula. As the initial (residual) stresses tend to zero, these exact solutions of linearized solid mechanics transform into the respective exact solutions of classical linear solid mechanics based on the Muskhelishvili, Lekhnitskii, and Galin complex representations. New mechanical effects in the dynamic problems under consideration are analyzed. The influence of initial (residual) stresses and crack propagation rate is established. In addition, the following two related problems are briefly analyzed within the framework of linearized solid mechanics: growing cracks at the interface of two materials with initial (residual) stresses and brittle fracture under compression along cracks     

Numerical analysis of quasi-static crack branching in brittle solids by a modified displacement discontinuity method

Mechanism of quasi-static crack branching in brittle solids has been analyzed by a modified displacement discontinuity method. It has been assumed that the pre-existing cracks in brittle solids may propagate at the crack tips due to the initiation and propagation of the kink (or wing) cracks. The originated wing cracks will act as new cracks and can be further propagated from their tips according to the linear elastic fracture mechanics (LEFM) theory. The kink displacement discontinuity formulations (considering the linear and quadratic interpolation functions) are specially developed to calculate the displacement discontinuities for the left and right sides of a kink point so that the first and second mode kink stress intensity factors can be estimated. The crack tips are also treated by boundary displacement collocation technique considering the singularity variation of the displacements and stresses near the crack tip. The propagating direction of the secondary cracks can be predicted by using the maximum tangential stress criterion. An iterative algorithm is used to predict the crack propagating path assuming an incremental increase of the crack length in the predicted direction (straight and curved cracks have been treated). The same approach has been used for estimating the crack propagating direction and path of the original and wing cracks considering the special crack tip elements. Some example problems are numerically solved assuming quasi-static conditions. These results are compared with the corresponding experimental and numerical results given in the literature. This comparison validates the accuracy and applicability of the proposed method.     

裂纹体分析的权函数理论与应用: 回顾和展望

断裂力学是工程材料和结构的疲劳与断裂分析、损伤容限设计和结构完整性评定的理论基础. 应力强度因子作为线弹性裂纹尖端奇异场的单一表征参量和裂纹扩展驱动力, 在裂纹体的断裂力学分析中发挥着关键作用. 权函数法为复杂受载裂纹体的应力强度因子求解计算提供了强有力的解析工具, 不但具有远高于各类数值解法的计算效率, 而且精度可靠, 使用方便. 本文结合笔者团队在权函数法方面的长期研究工作, 对该方法自20世纪70年代初提出至今半个世纪以来, 国际断裂界在二维和三维权函数理论与应用方面的主要研究进展作了回顾和评述, 并对其未来发展提出了展望. 主要内容涵盖: 当前国际断裂界广泛应用的3种二维裂纹解析权函数法简介和以格林函数为基准的验证评价; 三维裂纹问题的片条合成权函数法和点载荷权函数法; 权函数法在复杂受载裂纹体的应力强度因子和裂纹张开位移等关键力学参量计算、内聚力/桥连等裂纹模型分析、共线多裂纹权函数理论及其在剩余强度预测等方面的应用, 以及复杂裂纹几何的工程化权函数分析和权函数法的反向应用问题.      

ON PROBLEMS OF FRACTURE OF MATERIALS IN COMPRESSION ALONG TWO INTERNAL PARALLEL CRACKS

I.Relati0ns0fN0n-LinearThe0ryandDerivati0nofLinearizedRelati0nsinCoordinates0fN0n.DeformedStateNotationsareintroduced:x,-x'-Lagrangiancoordinateswhichinthenatural(non-def0rmed)statecoincidewithCartesiancoordinateswith0rthse,.Coordinatesx,willbeassumedtobe…     

Analysis of rotating solids with cracks by the boundary element method

The use of the boundary element method for the solution of linear elastic fracture mechanics problems, without body forces, is quite extensive since the method is intrinsically well suited to the analysis of high stress gradients associated with crack problems. The crack-tip stresses for rotating bodies are similar to the stresses for stationary bodies and therefore all the advantages of the boundary element procedure can be encompassed in the extension of the method to the solution of rotating bodies with cracks. In the present analysis, the additional volume integral that arises from the treatment of inertial body forces is eliminated by using the well-known particular integral procedure. The matrix ill-conditioning that results from the need to model co-planar crack surfaces of non-symmetrical cracks is avoided by using the multi-region approach. The accuracy of the numerical solutions is improved by utilizing quarter-point elements with traction singular enhancement at the crack-tip. The procedure is then applied to the solution of arbitrary cracks in two- and three-dimensional rotating bodies.     

Combined analysis of fracture under stresses acting along cracks

A joint approach to the study of two non-classical fracture mechanisms, namely fracture of cracked materials with initial (residual) stresses acting along the crack planes and fracture under compression along parallel cracks, is considered in the framework of three-dimensional linearized solid mechanics. Mathematical statements of problems for pre-stressed solids that contain interacting circular cracks are given. Problems for an infinite solid containing two parallel co-axial cracks and for a space with the periodical set of co-axial parallel cracks as well as for a half-space with near-the-surface crack are solved. Several patterns of loading on the crack faces (normal loading, radial shear and torsion) are considered. The effects of initial stresses on stress intensity factors are analyzed for highly elastic materials with some types of elastic potentials. Formulation of fracture criteria accounting effect of initial (residual) stresses is given. Critical parameters of fracture of solids containing interacting cracks under compression along the cracks are calculated. The influence of geometrical parameters of the problems as well as physical and mechanical properties of materials on these critical parameters is analyzed.     

Stress intensity factors for three-dimensional cracks in functionally graded materials using enriched finite elements

A three-dimensional enriched finite-element methodology is presented to compute stress intensity factors for three-dimensional cracks contained in functionally graded materials (FGMs). A general-purpose 3D finite-element based fracture analysis program, FRAC3D, is enhanced to include this capability. First, using available solutions from the literature, comparisons have been made in terms of stresses under different loading conditions, such as uniform tensile, bending and thermal loads. Mesh refinement studies are also performed. The fracture solutions are obtained for edge cracks in an FGM strip and surface cracks in a finite-thickness FGM plate and compared with existing solutions in the literature. Further analyses are performed to study the behavior of stress intensity factor near the free surface where crack front terminates. It is shown that three-dimensional enriched finite elements provide accurate and efficient fracture solutions for three-dimensional cracks contained in functionally graded materials.     

Thermoelastoplastic deformation of a thick-walled cylinder with a radial crack

The thermoelastoplastic fracture mechanics problem of a thick-walled cylinder subjected to internal pressure and a nonuniform temperature field is solved by the method of elastic solutions combined with the finite-element method. The correctness of the solution is provided by using the Barenblatt crack model, in which the stress and strain fields are regular. The elastoplastic problem of a cracked cylinder subjected to internal pressure and a nonuniform temperature field are solved. The calculation results are compared with available data. __________ Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 49, No. 3, pp. 173–183, May–June, 2008.     

Inverse problem of fracture mechanics for a compound cylinder

We consider the plane problem of fracture mechanics for a compound cylinder. We assume that the sleeve (internal cylinder) is reinforced with negative allowance by the external cylinder and that near the sleeve surface there are N arbitrarily located rectilinear cracks of length 2l _k (k = 1, 2, ..., N). A minimax criterion is used to determine the negative allowance in the junction theoretically minimizing the fracture parameters (the stress intensity coefficients) of the compound cylinder. A simplified method for minimizing the fracture parameters of the compound cylinder is considered separately.     

Three-dimensional elasticity problems related to cracks: Exact solution by inversion transformation

Prior to determining the conditions of brittle and quasi-brittle fracture of elastic solids with cracks, it is necessary to solve the corresponding three-dimensional elasticity problem. Since analytical solutions are known only for certain simple configurations such as an infinite space containing a plane, penny-shaped, or elliptical crack, and a half-plane or a strip crack [1–3], numerical procedure may have to be used.This paper is concerned with a tensile crack (or Mode I crack) on a plane domain in infinite elastic medium. A technique is proposed for constructing exact solutions for crack configurations that are obtainable by inversion transformation from the canonical contours such as a circle, ellipse or half-plane, for which exact solutions are known. Simple formulae are derived which yield the stress intensity factor distribution over the transformed crack boundary with loads appropriately adjusted according to the initial crack region. Hence, it is not necessary to find the complete solution. Solutions obtained by this method are presented for cracks bounded by convex as well as nonconvex contour, i.e., oval and sickle-shape contours.     

Effect of contact interaction of the crack faces for a crack under harmonic loading

This paper is concerned with dynamic problems in fracture mechanics for elastic solids having cracks with contacting faces. The contact problem for a penny-shaped crack with a nonzero initial opening under normally incident harmonic wave is solved by the method of boundary integral equations. The solutions are compared with those that neglect the contact interaction of the crack faces. Results are presented for different values of the initial crack opening Presented at the 6th International Conference on Modern Practice in Stress and Vibration Analysis (Bath, United Kingdom, September 5–7, 2006). Published in Prikladnaya Mekhanika, Vol. 43, No. 7, pp. 125–131, July 2007.     

Crack kinking in sandwich structures under three-point bending

Sandwich beams under three-point bending containing cracks in the core material very close to the upper skin interface are investigated. The cracks considered parallel or with an imperceptible inclination to the longitudinal beam axis and at different distances from the upper skin interface, are analyzed with static non-linear elastic two and three dimensional finite element analyses. From the proposed analyses stress intensity factors are calculated using linear elastic fracture mechanics. It is demonstrated that the crack propagation on the compression side of the core is mainly subjected in shear. The strain energy density criterion is used in order to determine the angle of kinking of the crack into the core.     

Pressure effects and fracture of a rubbery particulate composite

The effects of pressure on the fracture properties of a rubbery particulate composite are studied. The results indicate elevated fracture toughness values and lower crack growth rates relative to ambient pressure conditions. The regression method employed for determining toughness establishes a minimum thickness requirement for application of linear elastic fracture mechanics concepts. Also, it is demonstrated that single edge notched tension specimens can be successfully used to predict the fracture initiation of specimens with semielliptical surface cracks.     

Longitudinal and Bending Stresses in a Beam with Saw Cut Edge Cracks∗

Abstract

The Griffith-Irwin theory of brittle fracture of elastic solids predicts the propagation of cracks on the basis of the energy release rate. This depends upon the stress intensity factors for a given crack configuration. The present paper provides these informations for the problem of an infinite number of periodic, non-coplanar, parallel edge cracks in a strip. Two types of crack configurations, namely, periodic cracks of equal length starting from one edge and a set of two coplanar symmetrical edge cracks of equal length are solved for constant and linearly varying pressure distributions. These problems arise naturally in structural mechanics while investigating stresses in extension and bending of cracked strips. Final results are obtained from the numerical solution of certain Fredholm integral equations of the second kind derived from a dual series of Papkovich-Fadle eigenfunctions     

表面裂纹疲劳扩展的数值模拟

建立了一种无形状约束的模拟表面裂纹在线弹性断裂力学条件下疲劳扩展的数值方法,并研究了表面疲劳裂纹形状演化和裂纹尖端应力强度因子(SIF)的分布特征。该方法以三维有限单元技术和Paris疲劳裂纹扩展规律为基础,并在裂纹扩展增量计算中考虑了裂纹闭合影响。本文第一部分主要介绍模拟三维疲劳裂纹扩展的数值方法的理论背景和相关的技术细节。着重分析和讨论基于三维有限单元法计算裂纹SIF所涉及的几个主要问题:裂纹尖端单元网格密度对估算精度的影响;自由表面的影响及其修正方法;裂纹尖端非正交单元网格的影响及修正方法。     

多裂纹扩展的数值流形法

数值流形法的求解体系建立在两套覆盖(包括数学覆盖和物理覆盖) 和接触环路的基础之上,实现了对连续和非连续问题的统一求解. 在处理裂纹问题时,数学覆盖无需与裂纹重合,方便岩体破坏过程的模拟. 通过在裂纹尖端影响区域内的物理片上增加用于模拟应力奇异性的增强位移函数,发展了扩展的数值流形法. 在此基础上,提出一种多裂纹扩展的控制算法,并给出了裂纹扩展过程中材料体的整体响应. 针对典型的线弹性断裂力学问题, 给出的数值算例表明所建议的方法是正确有效的.      

Combination of fracture and damage mechanics for numerical failure analysis

A new approach for the analysis of crack propagation in brittle materials is proposed, which is based on a combination of fracture mechanics and continuum damage mechanics within the context of the finite element method. The approach combines the accuracy of singular crack-tip elements from fracture mechanics theories with the flexibility of crack representation by softening zones in damage mechanics formulations. A super element is constructed in which the typical elements are joined together. The crack propagation is decided on either of two fracture criteria; one criterion is based on the energy release rate or the J-integral, the other on the largest principal stress in the crack-tip region. Contrary to many damage mechanics methods, the combined fracture⧹damage approach is not sensitive to variations in the finite element division. Applications to situations of mixed-mode crack propagation in both two- and three-dimensional problems reveal that the calculated crack paths are independent of the element size and the element orientation and are accurate within one element from the theoretical (curvilinear) crack paths.     

On discrete models in fracture mechanics

In two papers published in 1969, where the foundations of the theory of brittle fracture were laid, Valentin Valentinovich Novozhilov introduced a necessary and sufficient criterion for estimating the strength of an elastic body weakened by a cut and formulated the concept of brittle fracture. The criterion is based on taking account of the body discrete structure, and the fracture process is treated as loss of stability of elastic equilibrium. These two facts, fracture discreteness and loss of stability in the deformation of the breaking bond, underlie numerous phenomena, which could not be discovered in the framework of continuum mechanics. In the present paper, the following effects of taking these facts into account are briefly discussed: an unstable leg of the stress-strain diagram and the capture of part of the fracture energy by the lattice, the role of the lattice dimension, radiation in quasi-static and dynamic crack propagation including the radiation structure, the role of the dynamic factor, irregularities in the crack growth, crack propagation under the action of high-frequency waves including steady-state modes with periodically varying speed, and cracks in continuum discrete models.