Shape and topology sensitivity analysis for cracks in elastic bodies on boundaries of rigid inclusions

We consider an elastic body with a rigid inclusion and a crack located at the boundary of the inclusion. It is assumed that nonpenetration conditions are imposed at the crack faces which do not allow the opposite crack faces to penetrate each other. We analyze the variational formulation of the problem and provide shape and topology sensitivity analysis of the solution in two and three spatial dimensions. The differentiability of the energy with respect to the crack length, for the crack located at the boundary of rigid inclusion, is established.     

Problem of a crack on the boundary of a rigid inclusion in an elastic plate

The problem of equilibrium of a thin elastic plate containing a rigid inclusion is considered. On part of the interface between the elastic plate and the rigid inclusion, there is a vertical crack. It is assumed that, on both crack edges, the boundary conditions are given as inequalities describing the mutual impenetrability of the edges. The solvability of the problem is proven and the character of satisfaction of the boundary conditions is described. It is also shown that the problem is the limit problem for a family of other problems posed for a wider region and describing equilibrium of elastic plates with a vertical crack as the rigidity parameter tends to infinity.     

The effect of an elastic triangular inclusion on a crack

IntroductionWiththedevelopmentofparticleandfiberreinforcedcomposites,theinclusion_crackinteractionproblemisbecominganimportantfieldbeingstudied .Andasamodel,itisalsousedtostudytheeffectsofmaterialdefectsonthestrengthandfractureofengineeringstructure.TheinterationbetweencircularinclusionandcrackwasstudiedinRefs.[1 -6 ] ;InRefs.[7-1 2 ] ,theinterationbetweenlineinclusionandcrackswasdiscussed ;TheinterationbetweenellipticalinclusionandcrackwasstudiedinRefs.[1 3,1 4] .However,withthedevelopmento…     

Interaction between a compliant disk-shaped inclusion and a crack upon incidence of an elastic wave

The propagation of harmonic elastic wave in an infinite three-dimensional matrix containing an interacting low-rigidity disk-shaped inclusion and a crack. The problem is reduced to a system of boundary integral equations for functions that characterize jumps of displacements on the inclusion and crack. The unknown functions are determined by numerical solution of the system of boundary integral equations. For the symmetric problem, graphs are given of the dynamic stress intensity factors in the vicinity of the circular inclusion and the crack on the wavenumber for different distances between them and different compliance parameters of the inclusion.     

A crack along the interface of a circular inclusion embedded in an infinite solid

The two-dimensional problem of an arc shaped crack lying along the interface of a circular elastic inclusion embedded in an infinite matrix with different elastic constants is considered. Based on the complex variable method of Muskhelishvili, closed-form solutions for the stresses and the displacements around the crack are obtained when general biaxial loads are applied at infinity. These solutions are then combined with A.A. Griffith's virtual work argument to give a criterion of crack extension, namely the de-bonding of the interface. The critical applied loads are expressed explicitly in terms of a function of the inclusion radius and the central angle subtended by the crack arc. In the case of simple tension the critical load is inversely proportional to the square-root of the inclusion radius. By analyzing the variation of the cleavage stress near the crack tip, the deviation of the crack into the matrix is discussed. The case of uniaxial tension is worked out in detail.     

The Method of R-Functions in the Solution of Elastic Problems on the Basis of Reissner's Mixed Variational Principle

A method is presented for solving boundary-value elastic problems on the basis of the variational–structural method of R-functions and Reissner's mixed variational principle. A mathematical formulation is given to problems on the deformation of elastic bodies under mixed boundary conditions and bodies interacting with smooth rigid dies. Solutions satisfying all the boundary conditions are proposed. For undetermined components of these solutions, the resolving equations are derived and their properties are studied. A posteriori estimation of numerical solutions is made. As examples, solutions are found to a problem on the stress–strain state of a short cylinder and to a contact problem on a cylinder interacting with a smooth die. A numerical method of solving such problems is analyzed for convergence, and the accuracy of the solutions is estimated.     

Strain gradient solution for a finite-domain Eshelby-type anti-plane strain inclusion problem

A solution for the finite-domain Eshelby-type inclusion problem of a finite elastic body containing an anti-plane strain inclusion of arbitrary cross-sectional shape prescribed with a uniform eigenstrain and a uniform eigenstrain gradient is derived in a general form using a simplified strain gradient elasticity theory (SSGET). The formulation is facilitated by an extended Betti’s reciprocal theorem and an extended Somigliana’s identity based on the SSGET and suitable for anti-plane strain problems. The disturbed displacement field is obtained in terms of the SSGET-based Green’s function for an infinite anti-plane strain elastic body. The solution reduces to that of the infinite-domain anti-plane strain inclusion problem when the boundary effect is not considered. The problem of a circular cylindrical inclusion embedded concentrically in a finite cylindrical elastic matrix undergoing anti-plane strain deformations is analytically solved by applying the general solution, with the Eshelby tensor and its average over the circular cross section of the inclusion obtained in closed forms. This Eshelby tensor, being dependent on the position, inclusion size, matrix size, and a material length scale parameter, captures the inclusion size and boundary effects, unlike existing ones. It reduces to the classical linear elasticity-based Eshelby tensor for the circular cylindrical inclusion in an infinite matrix if both the strain gradient and boundary effects are suppressed. Numerical results quantitatively show that the inclusion size effect can be quite large when the inclusion is small and that the boundary effect can dominate when the inclusion volume fraction is high. However, the inclusion size effect is diminishing with the increase of the inclusion size, and the boundary effect is vanishing as the inclusion volume fraction becomes sufficiently low.     

基于Hamilton体系的弹性力学问题的比例边界有限元方法

比例边界有限元方法是求解偏微分方程的一种半解析半数值解法。对于弹性力学问题,可采用基于力学相似性、基于比例坐标相似变换的加权余量法和虚功原理得到以位移为未知量的系统控制方程,属于Lagrange体系。但在求解时,又引入了表面力为未知量,控制方程属于Hamilton体系。因而,本文提出在比例边界有限元离散方法的基础上,利...     

考虑晶界滑错的金属多晶体响应的自洽方法

本文首先建立含有三种介质（各向异性基体、各向异性夹杂，界面层）的平面应变夹杂模型，将基体和夹杂位移场展开为多项式级数，假设界面层很薄，运用变分原理得出这一问题的近似解。将上述夹杂问题的解和ＨＩＬＬ自洽方法相结合，给出了考虑晶界滑错效应的金属多晶体弹塑性响应。     

Interaction between rigid-disc inclusion and penny-shaped crack under elastic time-harmonic wave incidence

Influence of a rigid-disc massive inclusion on a neighboring penny-shaped crack induced by the time-harmonic wave propagation in an infinite elastic matrix is investigated by the numerical solution of associated 3D elastodynamic problem. No restrictions on the mutual orientation of interacting objects and direction of wave incidence are assumed. The inclusion is perfectly bonded with a matrix and supposes the translations and rotations, the crack faces are load-free. Frequency-domain problem is reduced to a system of boundary integral equations (BIEs) relative to the interfacial stress jumps (ISJs) on the inclusion and the crack opening displacements (CODs). The subtraction technique in conjunction with mapping technique, under taking into account the structure of solution at the fronts of inclusion and crack, is applied for regularization of BIEs obtained. A discrete analogue of equations is constructed by using the collocation scheme. Numerical calculations are carried out for the grazing incidence of a plane P-wave on the crack, where the interacting inclusion is coplanar and perpendicular to the crack, and has the same radius. The shielding and amplification effects of inclusion are assessed by the analysis of mode-I stress intensity factor (SIF) in the crack vicinity depending on the wave number, incident wave direction, position of the crack front point, inclusion mass, crack-inclusion orientation and distance.     

Effect of bulk deformation and surface tension on the growth of a coherent inclusion

The object of this paper is the study of the moving boundary problem modeling the growth of a spherical solid inclusion in an infinite solid matrix. The displacement in bulk is assumed infinitesimal, while the phases are modeled as isotropic elastic bodies, and the interface structure is described by a constant surface tension. Existence of solutions is proved, and their asymptotic behavior in time is studied, with particular attention to the competition between surface tension and bulk deformation.     

Strain gradient solution for a finite-domain Eshelby-type plane strain inclusion problem and Eshelby’s tensor for a cylindrical inclusion in a finite elastic matrix

A solution for the finite-domain Eshelby-type inclusion problem of a finite elastic body containing a plane strain inclusion prescribed with a uniform eigenstrain and a uniform eigenstrain gradient is derived in a general form using a simplified strain gradient elasticity theory (SSGET). The formulation is facilitated by an extended Betti’s reciprocal theorem and an extended Somigliana’s identity based on the SSGET and suitable for plane strain problems. The disturbed displacement field is obtained in terms of the SSGET-based Green’s function for an infinite plane strain elastic body, which differs from that in earlier studies using the three-dimensional Green’s function. The solution reduces to that of the infinite-domain inclusion problem when the boundary effect is suppressed. The problem of a cylindrical inclusion embedded concentrically in a finite plane strain cylindrical elastic matrix of an enhanced continuum is analytically solved for the first time by applying the general solution, with the Eshelby tensor and its average over the circular cross section of the inclusion obtained in closed forms. This Eshelby tensor, being dependent on the position, inclusion size, matrix size, and a material length scale parameter, captures the inclusion size and boundary effects, unlike existing ones. It reduces to the classical elasticity-based Eshelby tensor for the cylindrical inclusion in an infinite matrix if both the strain gradient and boundary effects are not considered. Numerical results quantitatively show that the inclusion size effect can be quite large when the inclusion is very small and that the boundary effect can dominate when the inclusion volume fraction is very high. However, the inclusion size effect is diminishing with the increase of the inclusion size, and the boundary effect is vanishing as the inclusion volume fraction becomes sufficiently low.     

A new integral equation formulation of two-dimensional inclusion–crack problems

A new integral equation formulation of two-dimensional infinite isotropic medium (matrix) with various inclusions and cracks is presented in this paper. The proposed integral formulation only contains the unknown displacements on the inclusion–matrix interfaces and the discontinuous displacements over the cracks. In order to solve the inclusion–crack problems, the displacement integral equation is used when the source points are acting on the inclusion–matrix interfaces, whilst the stress integral equation is adopted when the source points are being on the crack surfaces. Thus, the resulting system of equations can be formulated so that the displacements on the inclusion–matrix interfaces and the discontinuous displacements over the cracks can be obtained. Based on one point formulation, the stress intensity factors at the crack tips can be achieved. Numerical results from the present method are in excellent agreement with those from the conventional boundary element method.     

GURTIN VARIATIONAL PRINCIPLE AND FINITE ELEMENT SIMULATION FOR DYNAMICAL PROBLEMS OF FLUID-SATURATED POROUS MEDIA

Based on the theory of porous media,a general Gurtin variational principle for theinitial boundary value problem of dynamical response of fluid-saturated elastic porous media isdeveloped by assuming infinitesimal deformation and incompressible constituents of the solid andfluid phase.The finite element formulation based on this variational principle is also derived.Asthe functional of the variational principle is a spatial integral of the convolution formulation,thegeneral finite element discretization in space results in symmetrical differential-integral equationsin the time domain.In some situations,the differential-integral equations can be reduced to sym-metrical differential equations and,as a numerical example,it is employed to analyze the reflectionof one-dimensional longitudinal wave in a fluid-saturated porous solid.The numerical results canprovide further understanding of the wave propagation in porous media.     

Boundary element analysis of interaction between an elastic rectangular inclusion and a crack

The interaction between an elastic rectangular inclusion and a kinked crack inan infinite elastic body was considered by using boundary element method. The new complexboundary integral equations were derived. By introducing a complex unknown function H(t)related to the interface displacement density and traction and applying integration by parts,the traction continuous condition was satisfied automatically. Only one complex boundaryintegral equation was obtained on interface and involves only singularity of order l/ r. Toverify the validity and effectiveness of the present boundary element method, some typicalexamples were calculated. The obtained results show that the crack stress intensity factorsdecrease as the shear modulus of inclusion increases. Thus, the crack propagation is easiernear a softer inclusion and the harder inclusion is helpful for crack arrest.     

A rigid triangular inclusion partially bonded in an elastic matrix

The two-dimensional problem of a rigid rounded-off angle triangular inclusion partially bonded in an infinite elastic plate is studied. The unbonded part of the inclusion boundary forms an interfacial crack. Based on the complex variable method for curvilinear boundaries, the problem is reduced to a non-homogeneous Hilbert problem and the stress and displacement fields in the plate are obtained in closed form. Special attention is paid in the investigation of the stress field in the vicinity of the crack tip. It is found that the stresses present an oscillatory singularity and the general equations for the local stresses are derived. The singular stress field is coupled with the maximum circumferential stress and the minimum strain energy density criteria to study the fracture characteristics of the composite plate. Results are given for the complex stress intensity factors, the local stresses, the crack extension angles and the critical applied loads for unstable crack growth from its more vulnerable tip or two types of interfacial cracks along the inclusion boundary.     

摩擦约束塑性力学变分不等原理的半反推法

带摩擦约束的弹塑性接触问题,由于摩擦约束条件是一种判别性的条件,它的变分问题的逆问题的研究比较困难。本文对弹塑性接触力学中的变分不等问题的逆问题进行了研究,改进了半反推法并将其应用到弹塑性变分不等原理的研究中,导出了摩擦约束弹塑性增量广义变分不等原理中的能量泛函,消除了用拉氏乘子法可能产生的临界变分现象,在证明中,巧妙地处理了增量表示的接触摩擦边界条件,避免了使用非线性泛函分析和凸分析,简化了证明。     

Nonlinear buckling of a spherical shell embedded in an elastic medium with imperfect interface

This work analyzes nonlinear buckling of a single spherical shell imperfectly bonded to an infinite elastic matrix under a compressive remote load. The inclusion is modeled using a nonlinear shell formulation and the matrix is treated as a linear elastic body. Imperfect bonding conditions are realized through a linear spring interface model. A variational method is used to derive the governing differential equations, which are cast into a tractable set of nonlinear algebraic equations using the Galerkin method. An incremental iterative technique based on the modified Newton–Raphson method is employed to find the critical load of the system. The accuracy and convergence properties of the proposed method are validated through finite element analysis. The study is relevant to the analysis of compressive failure of syntactic foams used in marine and aerospace applications. Results are specialized to glass particle-vinyl ester matrix syntactic foams to test the hypothesis as to whether microballoons’ buckling is a dominant failure mechanism in such composites under compression. Parametric studies are conducted to understand the effect of interfacial properties and inclusion wall thickness on the overall mechanical behavior of the composite. Comparisons between analytical findings and experimental results on compressive response of syntactic foams and isolated microballoons indicate that inclusion buckling is unlikely a determinant of compressive failure in vinyl ester-glass systems. In particular, the matrix is found to exert a beneficial stabilizing effect on the inclusions, which fail under brittle fracture before the onset of buckling.     

Compliance and Hill polarization tensor of a crack in an anisotropic matrix

This work aims at developing an efficient method to compute the compliance due to a crack modeled as a flat ellipsoid of any shape in an infinite elastic matrix of arbitrary anisotropy (Eshelby problem) when no closed-form solution seems currently available. Whereas the solution of this problem usually requires the calculation of the so-called fourth-order Hill polarization tensor if the ellipsoid is not singular, it is shown that the crack compliance can be derived from the first-order term in the Taylor expansion of the Hill tensor with respect to the smallest aspect ratio of the ellipsoidal inclusion. For a 3D ellipsoidal crack model, this first-order term is expressed as a simple integral thanks to the Cauchy residue theorem. A similar method allows to express the same term in the case of a cylindrical crack model without any integral. A numerical example is finally treated.     

An analytical approach for the calculation of stress-intensity factors in transformation-toughened ceramics

Stress-induced transformation toughening in Zirconia-containing ceramics is described analytically by means of a quantitative model: A Griffith crack which interacts with a transformed, circular Zirconia inclusion. Due to its volume expansion, a ZrO₂-particle compresses its flanks, whereas a particle in front of the crack opens the flanks such that the crack will be attracted and finally absorbed. Erdogan's integral equation technique is applied to calculate the dislocation functions and the stress-intensity-factors which correspond to these situations. In order to derive analytical expressions, the elastic constants of the inclusion and the matrix are assumed to be equal.