沿两条内部平行裂纹受压时材料的断裂问题

顺着平面内的平行裂纹受压材料的断裂问题并不能在线性断裂力学的框架内进行描述，Grif-fith－Irvin型或COC型的断裂判据，虽然可以用来处理经典的线性断裂力学，但对本题则完全不适用，这是因为这些压力对应力强度系数没有影响，与裂缝孔隙值也没有关系[1，2]，这一类问题只能采用新的方法，本文的第一作者曾建议过一种新方法，在这一方法中仍然使用了线性关系，但这种线性关系是从变形固体力学中的非线性方程导出的[3，4，5]．这里必须指出，这种方法曾在变形体稳定性问题中广泛地采用过。作为断裂开始的判据，我们采用了裂缝缺陷附近的局部失稳的判据，在这类情况下，我们认为是从失稳过程引发断裂过程的。     

The solution of a problem in contact fracture mechanics on the nucleation and development of a bridged crack in the hub of a friction pair

The contact deformation of the hub of a plunger pair is considered. It is assumed that, during the repeated reciprocating motion of the plunger, a crack is initiated and fracture of the materials of the elements of the contact pair occurs. The problem of the equilibrium of the hub of a friction pair with a crack nucleus reduces to solving a system of non-linear singular integrodifferential equations with a Cauchy-type kernel. The normal and shear forces in the zone where the crack originates are found from the solution of this system of equations. The condition for the appearance of a crack is formulated, taking account of the criterion of the limit traction of the bonds in the material. A problem for the plunger of a friction pair as applied to a borehole sucker rod pump is considered as an example. In conclusion, the case when there are several arbitrarily distributed rectilinear bridged cracks, with bonds between the crack faces in the end zone, close to the contact surface of the hub is investigated.     

Equilibrium cracks in composites reinforced with unidirectional fibres

A discrete-continuum approach, proposed by Novozhilov for analysing the equilibrium states of a brittle of crack in an isotropic body, is applied to a penny-shaped crack situated in a fibre-reinforced composite perpendicular to the fibres. The structural non-uniformity of the material is taken into account by the presence of unbroken fibres in the narrow part of the crack, adjoining the edge, and the different effect of the strength properties of the fibre and matrix on the limit state of the crack. Using this model, the range of dimensions of equilibrium cracks is established and an estimate is given of the critical size of the bridged part of the crack, corresponding to the onset of catastrophic fracture. It is shown that this dimension has the same value for a penny-shaped crack and for a crack under plane strain, does not depend on the form of the load and, under the condition of its smallness, is a brittle fracture characteristic of a fibre-reinforced material. The possibility of using this fracture model for two types of ceramics is analysed on the basis of experimental data.     

Phase field modeling of complex crack patterns in multi-physics problems

Recently developed continuum phase field models for brittle fracture show excellent modeling capability in situations with complex crack topologies including branching in the small and large strain applications. This work presents a generalization towards fully coupled multi-physics problems at large strains. A modular concept is outlined for the linking of the diffusive crack modeling with complex multi field material response, where the focus is put on the model problem of finite thermo-elasticity. This concerns a generalization of crack driving forces from the energetic definitions towards stress-based criteria, the constitutive modeling of degradation of non-mechanical fluxes on generated crack faces. Particular assumptions are made on the generation of convective heat exchanges approximating surface load integrals of the sharp crack approach by distinct volume integrals. The coupling effect is also shown in generation of cracks due to thermally induced stress states. We finally demonstrate the performance of the phase field formulation of fracture at large strains by means of representative numerical examples. (© 2014 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Contact Interaction of Faces of the Plane Crack under Harmonic Loading

The present work is devoted to the solution of the three-dimensional fracture mechanics problem for a linear elastic, homogeneous and isotropic solid with a stationary plane crack under normal time-harmonic loading. The problem has been solved by the method of boundary integral equations with the allowance for the contact interaction of the opposite faces of the crack. (© 2005 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Thermo-mechanical modeling of turbine blades taking into account structural defects

In turbine blades of aero-engines typical defects are cracks due to high mechanical and thermal loads. The extended finite element method (XFEM) is used for simulations of fracture mechanics problems with cracks. Discontinuities in the displacement and temperature field are allowed and the crack opening displacement and crack tip stress field are reproduced accurately. Since crack closure and non-physical penetration of the crack surfaces may occur under certain load conditions, it becomes necessary to enforce the non-penetration condition for crack surfaces. This contact formulation is assumed to be frictionless. The node-to-segment approach proposed in [3] is extended to ten-node tetrahedral elements with quadratic shape functions. (© 2017 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

A phase field model for fracture

The variational formulation of brittle fracture as formulated for example by Francfort and Marigo in [1], where the total energy is minimized with respect to any admissible crack set and displacement field, allows the identification of crack paths, branching of preexisting cracks and even crack initiation without additional criteria. For its numerical treatment a continuous approximation of the model in the sense of Γ-convergence has been presented by Bourdin in [2]. In the regularized Francfort–Marigo model cracks are represented by an additional field variable (secondary variable) s∈[0,1] which is 0 if the material is cracked and 1 if it is undamaged. In this work, we reinterpret the crack variable as a phase field order parameter and address cracking as a phase transition problem. The crack growth is governed by the evolution equation of the order parameter which resembles the Ginzburg–Landau equation. The numerical treatment is done by finite elements combined with an implicit Euler scheme for the time integration. (© 2008 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Transverse cracks in cross-ply laminates. 1. Stress analysis

During service loading of cross-ply laminates, transverse cracks occur in plies. The cracks parallel to the fiber direction are extended over the full thickness of transverse plies and often cross the entire test specimen width. It is widely recognized that the changes of laminate thermomechanical constants, caused by the transverse cracking of composite laminates, can be significant. Theoretical stress analysis in the cross-ply laminates in the vicinity of cracks is performed using numerical (FE) and analytical methods. The effect of transverse cracks on the degradation of elastic properties will be discussed in Part 2 [1]. Approximate analytical micromechanical models based on shear lag predictions, variational analysis, and numerical 2D finite element calculations were verified in their predictive abilities. The three variational models used are based on the principle of minimum complementary energy and have different degrees of accuracy with respect to the stress assumptions used (Hashin's, 2D 0° and 2D 0°/90° models). Using FEM, the plane stress and strain state were analyzed. The effect of material properties and layer thickness on the stress distribution in a 90° layer was evaluated by varying the crack spacing. The crack opening displacement (COD), normalized with respect to the far field strain, is proposed as a measure of reduction of the mechanical properties. Since the CODs are rather insensitive to the crack spacing (crack density) in a wide region, they will be used in modeling the stiffness reduction in these laminates [1].Translated from Mekhanika Kompozitnykh Materialov, Vol. 33, No. 6, pp. 796–820, November–December, 1997.     

The contact problem of elasticity theory for bodies with cracks

The plane contact problem of a stamp impressed into an elastic half-plane containing arbitrarily arranged rectilinear subsurface cracks is formulated and investigated by asymptotic methods. Partial or total overlapping of the crack edges is assumed. The problem reduces to a system of linear singular integrodifferential equations with side conditions in the form of equalities and inequalities. An analytic solution of the problem is obtained in the form of asymptotic power series /1/ in the relative dimension of the greatest crack. Dependences of the first terms of the asymptotic expansions of the desired functions on the mutual location of the cracks and the contact domains, the pressure and friction stress distributions, and the crack size and orientation are determined. Numerical results are presented.

Analysis of the influence of the stress-free boundary of the half-plane on the state of stress and strain of the elastic material near the cracks is presented in /2, 3/. The problem of a crack in an elastic plane whose edges overlap partially is also examined in /3/ by numerical methods.     

压电复合材料圆形夹杂中螺型位错和界面裂纹的干涉分析

研究了无穷远纵向剪切和面内电场共同作用下,压电复合材料圆形夹杂中螺型位错与界面裂纹的电弹耦合干涉作用．运用Riemann-Schwarz 对称原理,并结合复变函数奇性主部分析方法,获得了该问题的一般解答．作为典型算例,求出了界面含一条裂纹时基体和夹杂区域复势函数和电弹性场的封闭形式解．应用广义Peach-Koehler公式,导出了位错力的解析表达式．分析了裂纹几何参数和材料的电弹性常数对位错力的影响规律．结果表明,界面裂纹对位错力和位错平衡位置有很强的扰动效应,当界面裂纹长度达到临界值时,可以改变位错力的方向．该结果可以作为格林函数研究圆形夹杂内裂纹和界面裂纹的干涉效应．其公式的退化结果与已有文献完全一致．     

Influence of initial stresses on fracture of composite materials containing interacting cracks

Considered in this study are the axially-symmetric problems of fracture of composite materials with interacting cracks, which are subjected to initial (residual) stresses acting along the cracks planes. An analytical approach within the framework of three-dimensional linearized mechanics of solids is used. Two geometric schemes of cracks location are studied: a circular crack is located parallel to the surface of a semi-infinite composite with initial stresses, and two parallel co-axial penny-shaped cracks are contained in an infinite composite material with initial stresses. The cracks are assumed to be under a normal or a radial shear load. Analysis involves reducing the problems to systems of second-kind Fredholm integral equations, where the solutions are identified with harmonic potential functions. Representations of the stress intensity factors near the cracks edges are obtained. These stress intensity factors are influenced by the initial stresses. The presence of the free boundary and the interaction between cracks has a significant effect on the stress intensity factors as well. The parameters of fracture for two types of composites (a laminar composite made of aluminum/boron/silicate glass with epoxy-maleic resin and a carbon/plastic composite with stochastic reinforcement by short ellipsoidal carbon fibers) are analyzed numerically. The dependence of the stress intensity factors on the initial stresses, physical-mechanical parameters of the composites, and the geometric parameters of the problem are investigated.     

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)     

Elastodynamic Contact Problems for Interface Cracks under Harmonic Loading

2-D fracture dynamics' problems for elastic bimaterials with cracks located at the bonding interface under the oblique time harmonic wave are considered in the study. The system of boundary integral equations for displacements and tractions is derived from Somigliana identity taking the contact interaction of the opposite crack faces into account. For the numerical solution the collocation method with piecewise constant approximation is used. The numerical results are obtained for various values of the angle of the wave incidence and the wave frequency taking the friction effects into account. (© 2011 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Towards Phase Field Modeling of Ductile Fracture in Gradient-Extended Elastic-Plastic Solids

The modeling of failure in ductile metals must account for complex phenomena at a micro-scale as well as the final rupture at the macro-scale. Within a top-down viewpoint, this can be achieved by the combination of a micro-structure-informed elastic-plastic model with a concept for the modeling of macroscopic crack discontinuities. In this context, it is important to account for material length scales and thermo-mechanical coupling effects due to dissipative heating. This can be achieved by the construction of non-standard, gradient-enhanced models of plasticity with a full embedding into continuum thermodynamics [1,2]. The modeling of macroscopic cracks can be achieved in a convenient way by recently developed continuum phase field approaches to fracture based on regularized crack discontinuities. This avoids the use of complex discretization methods for crack discontinuities, and can account for complex crack patterns within a pure continuum formulation. Moreover, the phase field modeling of fracture is related to gradient theories of continuum damage mechanics, and fits nicely the structure of constitutive models for gradient plasticity. The main focus of this work is the extensions to gradient thermoplasticity and phase field formulation of ductile fracture, conceptually in line with the work [3]. (© 2014 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Ductile failure with gradient plasticity coupled to the phase-field fracture at finite strains

Most metals fail in a ductile fashion, i.e, fracture is preceded by significant plastic deformation. The modeling of failure in ductile metals must account for complex phenomena at micro-scale, such as nucleation, growth and coalescence of micro-voids. In this work, we start with von-Mises plasticity model without considering void generation. The modeling of macroscopic cracks can be achieved in a convenient way by the continuum phase field approaches to fracture, which are based on the regularization of sharp crack discontinuities [1]. This avoids the use of complex discretization methods for crack discontinuities and can account for complex crack patterns. The key aspect of this work is the extension of the energetic and the stress-based phase field driving force function in brittle fracture to account for a coupled elasto-plastic response in line with our recent work [3]. We develop a new theoretical and computational framework for the phase field modeling of ductile fracture in elastic-plastic solids. To account for large strains, the constitutive model is constructed in the logarithmic strain space, which simplify the model equations and results in a formulation similar to small strains. We demonstrate the modeling capabilities and algorithmic performance of the proposed formulation by representative simulations of ductile failure mechanisms in metals. (© 2015 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Delamination buckling of a rectangular orthotropic composite plate containing a band crack

The delamination buckling problem for a rectangular plate made of an orthotropic composite material is studied. The plate contains a band crack whose faces have an initial infinitesimal imperfection. The subsequent development of this imperfection due to an external compressive load acting along the crack is studied through the use of the three-dimensional geometrically nonlinear field equations of elasticity theory for anisotropic bodies. A criterion of initial imperfection is used in determining the critical forces. The corresponding boundary-value problems are solved by employing the boundary-form perturbation technique and the FEM. Numerical results for the critical force are presented.     

Evolutional contact with Coulomb's friction on a periodic microstructure

We consider the elasticity problem in a heterogeneous domain with an ε-periodic micro-structure, ε ≪ 1, including a multiple micro-contact in a simply connected matrix domain with inclusions completely surrounded by cracks, which do not connect the boundary, or a textile-like material. The contact is described by the Signorini and Coulomb-friction contact conditions. In the case of the Coulomb friction, the dissipative functional is state dependent, like in [2]. A time discretization scheme from [2] reduces the contact problem to the Tresca one (with prescribed frictional traction or state independent dissipation) on each time-increment. We further look for the spatial homogenization. The limiting energy and the dissipation term in the stability condition were obtained for the contact with Tresca's friction law in [4] for closed cracks and can be extended to textile-like materials. Using these results and the concept of energetic solutions for evolutional quasi-variational problems from [2], for a uniform time-step partition, the existence can be proved for the solution of the continuous problem and a subsequence of incremental solutions weakly converging to the continuous one uniformly in time. Furthermore, the irreversible frictional displacements at micro-level lead to a kind of an evolutional plastic behavior of the homogenized medium. (© 2013 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Crack – interface crack interactions in functionally graded/homogeneous composite bimaterials subjected to a heat flux

A bimaterial containing an interface crack and consisting of a homogeneous material and a functionally graded material (FGM) with a system of small internal cracks is considered. The thermal fracture of the biomaterial under the action of a heat flux applied to it at infinity is investigated. The problem is studied in the case where the interface crack is much larger than the internal ones. It is assumed that the thermal properties of the FGM are continuous functions of the thickness coordinate. Asymptotic analytical formulas for the thermal stress intensity factors (TSIFs) at the tips of the interface crack are obtained as series in a small parameter (the ratio between sizes of the internal and interface cracks). Then, the critical heat fluxes and the angles of propagation direction of the interface crack are calculated using the criterion of maximum circumferential stress. A parametric analysis shows that the propagation direction of the interface crack depends on the location and orientation of the system of internal cracks. The parameters of inhomogeneity of the FGM affect the value of TSIFs and, hence, the deflection angle of the interface crack.     

Unsteady propagation of longitudinal shear cracks

In all problems of unsteady crack propagation which have been solved to date [1 to 3], it has been assumed that the crack propagates at a constant speed. This assumption was not prompted by physical considerations of the problem, but by the methods of solution, therefore, the applicability of the results is limited. It would be more realistic to consider the speed of crack propagation as a function of time based on explicit physical hypotheses. Unfortunately, the general case of the resultant problem cannot be solved by existing methods. However, the problem of longitudinal shear cracks i.e. the plane problem in which the displacement is parallel to the crack boundary, may be solved for an arbitrary given variation in crack propagation speed, utilizing the method developed in connection with the theory of supersonic flows [4 and 5].

Note that equilibrium problems of longitudinal shear cracks have been studied in [6 and 7].