Crack nucleation sensitivity analysis

A simple analytical expression for crack nucleation sensitivity analysis is proposed relying on the concept of topological derivative and applied within a two‐dimensional linear elastic fracture mechanics theory (LEFM). In particular, the topological asymptotic expansion of the total potential energy together with a Griffith‐type energy of an elastic cracked body is calculated. As a main result, we derive a crack nucleation criterion based on the topological derivative and a criterion for determining the direction of crack growth based on the topological gradient. The proposed methodology leads to an axiomatic approach of crack nucleation sensitivity analysis. Copyright © 2010 John Wiley & Sons, Ltd.     

Derivatives of the energy functional for 2D‐problems with a crack under Signorini and friction conditions

We consider the two‐dimensional elasticity problem for an elastic body with a crack under unilateral constraints imposed at the crack. We assume that both the Signorini condition for non‐penetration of the crack faces and the condition of given friction between them are fulfilled. The problem is non‐linear and can be described by a variational inequality. Varying the shape of the crack by a local coordinate transformation of the domain, the first derivative of the energy functional to the problem with respect to the crack length is obtained, which gives the criterion for the crack growing. The regularity of the solution is discussed and the singular solution is performed. Copyright © 2000 John Wiley & Sons, Ltd.     

Initiation of Cracks in Griffith’s Theory: An Argument of Continuity in Favor of Global Minimization

The initiation of a crack in a sound body is a real issue in the setting of Griffith’s theory of brittle fracture. If one uses the concept of critical energy release rate (Griffith’s criterion), it is in general impossible to initiate a crack. On the other hand, if we replace it by a least energy principle (Francfort–Marigo’s criterion), it becomes possible to predict the onset of cracking in any circumstance. However this latter criterion can appear too strong. We propose here to reinforce its interest by an argument of continuity. Specifically, we consider the issue of the initiation of a crack at a notch whose angle ω is considered as a parameter. The result predicted by the Griffith criterion is not continuous with respect to ω, since no initiation occurs when ω>0 while a crack initiates when ω=0. In contrast, the Francfort–Marigo’s criterion delivers a response which is continuous with respect to ω, even though the onset of cracking is necessarily brutal when ω>0. The theoretical analysis is illustrated by numerical computations.     

含有界面裂纹的粘弹性层合板应变能释放率的计算

将裂纹扩展所对应的能量释放率定义为同一时刻,同样载荷条件下两种状态的能量之差.一是裂纹长度为a时,系统内能,第二状态是指裂纹长度为a+Δa时系统内能.这样,所定义的能量释放率相当于在无限短时间内,裂纹从a扩展到a+Δa所释放的能量.通过计算发现,对于给定的加载历史,应变能释放率是时间的函数,它的最大值相对应于层间开裂临界状态.在William工作的基础上,根据经典梁的理论求得双悬臂梁结构的应变能释放率的显函表达式.     

Shape sensitivity of a plane crack front

The 3D‐elasticity model of a solid with a plane crack under the stress‐free boundary conditions at the crack is considered. We investigate variations of a solution and of energy functionals with respect to perturbations of the crack front in the plane. The corresponding expansions at least up to the second‐order terms are obtained. The strong derivatives of the solution are constructed as an iterative solution of the same elasticity problem with specified right‐hand sides. Using the expansion of the potential and surface energy, we consider an approximate quadratic form for local shape optimization of the crack front defined by the Griffith criterion. To specify its properties, a procedure of discrete optimization is proposed, which reduces to a matrix variational inequality. At least for a small load we prove its solvability and find a quasi‐static model of the crack growth depending on the loading parameter. Copyright © 2003 John Wiley & Sons, Ltd.     

A computational framework of three dimensional configurational-force-driven crack propagation

A variational formulation of quasi-static brittle fracture is considered and a new finite-element-based computational framework is developed for propagation of cracks in three-dimensional bodies. We outline a consistent thermodynamical framework for crack propagation in elastic solids and show that the crack propagation direction associated with the classical Griffith criterion is identified by the material configurational force which maximizes the local dissipation at the crack front. The evolving crack discontinuity is realized by the doubling of critical nodes and triangular interface facets of the tetrahedral mesh. The crucial step for the success of the procedure is its embedding into an r-adaptive crack-facet reorientation procedure based on configurational-force-based indicators in conjunction with crack front constraints. We further propose a staggered algorithm which minimizes the stored energy at frozen crack state followed by the successive crack releases at frozen deformation. This constitutes a sequence of positive definite subproblems with successively decreasing overall stiffness, providing a very robust algorithmic setting in the postcritical range. (© 2008 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

On crack propagation shapes in elastic bodies

We consider a two dimensional elastic isotropic body with a curvilinear crack. The formula for the derivative of the energy functional with respect to the crack length is discussed. It is proved that this derivative is independent of the crack path provided that we consider quite smooth crack propagation shapes. An estimate for the derivative of the energy functional being uniform with respect to the crack propagation shape is derived.     

On crack propagation shapes in elastic bodies

We consider a two dimensional elastic isotropic body with a curvilinear crack. The formula for the derivative of the energy functional with respect to the crack length is discussed. It is proved that this derivative is independent of the crack path provided that we consider quite smooth crack propagation shapes. An estimate for the derivative of the energy functional being uniform with respect to the crack propagation shape is derived.     

On subcritical development of a shear crack in a composite with viscoelastic components

We present results of an investigation of the development of a transverse shear crack in a composite material with linearly viscoelastic components under external shear load. The solution is divided into the following two main stages: determination of the time dependence of the crack tip opening displacement and determination of the crack-growth kinetics as a result of the solution of integral equations. In the first stage, we use the solution of the corresponding elastic problem of determination of the crack opening displacement and the problem of determination of the effective moduli of the composite reinforced with unidirectional discrete fibers. Using the theoretically proved principle of elasto-viscoelastic analogy and the method of Laplace inverse transformation, we obtain a solution in a time domain. In the second stage, using the criterion of critical crack opening displacement for a transverse shear crack and an equation for the viscoelastic crack opening displacement of this crack, we construct an equation of crack growth. We present results of the numerical solution, which illustrate the influence of relations between the relaxation parameters of the materials of the components on the durability of the body with a crack.     

On elastic bodies with thin rigid inclusions and cracks

This paper is concerned with the analysis of equilibrium problems for two‐dimensional elastic bodies with thin rigid inclusions and cracks. Inequality‐type boundary conditions are imposed at the crack faces providing a mutual non‐penetration between the crack faces. A rigid inclusion may have a delamination, thus forming a crack with non‐penetration between the opposite faces. We analyze variational and differential problem formulations. Different geometrical situations are considered, in particular, a crack may be parallel to the inclusion as well as the crack may cross the inclusion, and also a deviation of the crack from the rigid inclusion is considered. We obtain a formula for the derivative of the energy functional with respect to the crack length for considering this derivative as a cost functional. An optimal control problem is analyzed to control the crack growth. Copyright © 2010 John Wiley & Sons, Ltd.     

Interface cracks in anisotropic composites

The linear model equations of elasticity often give rise to oscillatory solutions in some vicinity of interface crack fronts. In this paper we apply the Wiener–Hopf method which yields the asymptotic behaviour of the elastic fields and, in addition, criteria to prevent oscillatory solutions. The exponents of the asymptotic expansions are found as eigenvalues of the symbol of corresponding boundary pseudodifferential equations. The method works for three‐dimensional anisotropic bodies and we demonstrate it for the example of two anisotropic bodies, one of which is bounded and the other one is its exterior complement. The common boundary is a smooth surface. On one part of this surface, called the interface, the bodies are bonded, while on the complementary part there is a crack. By applying the potential method, the problem is reduced to an equivalent system of Boundary Pseudodifferential Equations (BPE) on the interface with the stress vector as the unknown. The BPEs are defined via Poincaré–Steklov operators. We prove the unique solvability of these BPEs and obtain the full asymptotic expansion of the solution near the crack front. As a special case we consider the interface crack between two different isotropic materials and derive an explicit criterion which prevents oscillatory solutions. Copyright © 1999 John Wiley & Sons, Ltd.     

Blister Patterns and Energy Minimization in Compressed Thin Films on Compliant Substrates

This paper is motivated by the complex blister patterns sometimes seen in thin elastic films on thick, compliant substrates. These patterns are often induced by an elastic misfit that compresses the film. Blistering permits the film to expand locally, reducing the elastic energy of the system. It is therefore natural to ask: what is the minimum elastic energy achievable by blistering on a fixed area fraction of the substrate? This is a variational problem involving both the elastic deformation of the film and substrate and the geometry of the blistered region. It involves three small parameters: the nondimensionalized thickness of the film, the compliance ratio of the film/substrate pair, and the mismatch strain. In formulating the problem, we use a small‐slope (Föppl–von Kármán) approximation for the elastic energy of the film, and a local approximation for the elastic energy of the substrate. For a one‐dimensional version of the problem, we obtain “matching” upper and lower bounds on the minimum energy, in the sense that both bounds have the same scaling behavior with respect to the small parameters. The upper bound is straightforward and familiar: it is achieved by periodic blistering on a specific length scale. The lower bound is more subtle, since it must be proved without any assumption on the geometry of the blistered region. For a two‐dimensional version of the problem, our results are less complete. Our upper and lower bounds only “match” in their scaling with respect to the nondimensionalized thickness, not in the dependence on the compliance ratio and the mismatch strain. The lower bound is an easy consequence of our one‐dimensional analysis. The upper bound considers a two‐dimensional lattice of blisters and uses ideas from the literature on the folding or “crumpling” of a confined elastic sheet. Our main two‐dimensional result is that in a certain parameter regime, the elastic energy of this lattice is significantly lower than that of a few large blisters. © 2015 Wiley Periodicals, Inc.     

Two arbitrarily located normal forces and a penny‐shaped crack: a complete solution

The following problem is considered: a penny‐shaped crack is located in the plane z=0 of a transversely isotropic elastic space and interacts with two equal and opposite normal forces, which are located arbitrarily, but symmetrically with respect to the plane of the crack. An exact closed‐form solution is obtained and expressed in terms of elementary functions for the fields of stresses and displacements in the whole space. This kind of problem deemed to be intractable by the methods of contemporary mathematical analysis, and has never been attempted before, even in the case of an isotropic body. Copyright © 1999 John Wiley & Sons, Ltd.     

Mathematical modelling and the study of the influence of initial stresses on the SIF and ERR at the crack tips in a plate-strip of orthotropic material

A plate-strip fabricated from the orthotropic material and containing a crack whose edges are parallel to the face planes of the plate is considered. It is assumed that the strip is stretched (or compressed) initially along the crack edges by uniformly distributed external normal forces acting on the simply supported ends of the plate-strip. After this initial stretching (or compression) the crack edges are loaded by additional uniformly distributed normal (opening) forces. As a result of the action of these additional forces the stress concentration characterized by the stress intensity factor (SIF) of mode I or by the energy release rate (ERR) of mixed mode arises at the crack tips. In this paper, the influence of the initial stresses on the SIF or ERR is modelled mathematically by the use of the three-dimensional linearized theory of elasticity. The aim of the present investigations is to study the effect of the mechanical–orthotropic properties of the plate-strip material on this influence by the use of the finite element method (FEM) modelling of the corresponding boundary-value problem.     

Statistical methods for automatic crack detection based on vibrothermography sequence‐of‐images data

Vibrothermography is a relatively new nondestructive evaluation technique for finding cracks through frictional heat generated from crack surface vibrations under external excitations. The vibrothermography inspection method provides a sequence of infrared images as the output. We use a matched filter technique to increase the signal‐to‐noise ratio of the sequence‐of‐images data. An automatic crack detection criterion based on the features extracted from the matched filter output greatly increases the sensitivity of the vibrothermography inspection method. In this paper, we develop a three‐dimensional matched filter for the sequence‐of‐images data, which presents the statistical analysis for the matched filter output, and evaluate the probability of detection. Our results show the crack detection criterion based on the matched filter output provides an improved detection capability. Copyright © 2010 John Wiley & Sons, Ltd.     

一维六方准晶双材料中圆孔边共线界面裂纹的反平面问题

研究了一维六方准晶双材料中圆孔边不对称共线界面裂纹的反平面问题。利用Stroh公式和复变函数方法得到了声子场和相位子场耦合作用下的复势函数,给出了裂纹尖端应力强度因子和能量释放率的解析表达式。通过数值算例,讨论了圆孔半径和裂纹长度对应力强度因子的影响,以及耦合系数、声子场应力和相位子场应力对能量释放率的影响。结果表明:当圆孔半径不变时,应力强度因子随右裂纹长度的增大趋向稳定值。当相位子场应力取一定值时,能量释放率达到最小值,说明特定的相位子场应力可以抑制裂纹的扩展。     

Uniform stabilization of a one‐dimensional hybrid thermo‐elastic structure

This paper is concerned with the stabilization of a one‐dimensional hybrid thermo‐elastic structure consisting of an extensible thermo‐elastic beam which is hinged at one end with a rigid body attached to its free end. The model takes account of the effect of stretching on bending and rotational inertia. The property of uniform stability of the energy associated with the model is asserted by constructing an appropriate Lyapunov functional for an abstract second order evolution problem. Critical use is made of a multiplier of an operator theoretic nature, which involves the fractional power A^?1/2 of the bi‐harmonic operator pair A acting in the abstract evolution problem. An explicit decay rate of the energy is obtained. Copyright © 2003 John Wiley & Sons, Ltd.     

Shape optimization of rigid inclusions for elastic plates with cracks

In the paper, we consider an optimal control problem of finding the most safe rigid inclusion shapes in elastic plates with cracks from the viewpoint of the Griffith rupture criterion. We make use of a general Kirchhoff–Love plate model with both vertical and horizontal displacements, and nonpenetration conditions are fulfilled on the crack faces. The dependence of the first derivative of the energy functional with respect to the crack length on regular shape perturbations of the rigid inclusion is analyzed. It is shown that there exists a solution of the optimal control problem.     

Junction problem for elastic and rigid inclusions in elastic bodies

An equilibrium problem for an elastic body is considered. It is assumed that the body has a thin elastic inclusion and a thin rigid inclusion. We analyze a junction problem assuming that the inclusions have a joint point. Different equivalent problem formulations are discussed, and existence of solutions is proved. A set of junction conditions is found. We investigate a convergence to infinity and to zero of a rigidity parameter of the elastic inclusion. A delamination of the elastic inclusion is also investigated. In this case, inequality‐type boundary conditions are imposed at the crack faces to prevent a mutual penetration between crack faces. Copyright © 2015 John Wiley & Sons, Ltd.     

Optimal Control of Inclusion and Crack Shapes in Elastic Bodies

The paper is concerned with the control of the shape of rigid and elastic inclusions and crack paths in elastic bodies. We provide the corresponding problem formulations and analyze the shape sensitivity of such inclusions and cracks with respect to different perturbations. Inequality type boundary conditions are imposed at the crack faces to provide a mutual nonpenetration between crack faces. Inclusion and crack shapes are considered as control functions and control objectives, respectively. The cost functional, which is based on the Griffith rupture criterion, characterizes the energy release rate and provides the shape sensitivity with respect to a change of the geometry. We prove an existence of optimal solutions.