Shape control of thin rigid inclusions and cracks in elastic bodies

The paper is concerned with a control of thin rigid inclusion and crack shapes in elastic bodies. It is assumed that rigid inclusions are delaminated; thus, cracks are located on the boundary of inclusions as well as outside of inclusions. We provide the problem formulations and analyze the shape sensitivity with respect to geometrical perturbations in the frame of free boundary models. Inequality type boundary conditions are considered at the crack faces to guarantee a mutual non-penetration between crack faces. Inclusion and crack shapes are considered as control functions. 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 of the structure. We prove an existence of optimal shapes in the problems considered.     

Timoshenko inclusions in elastic bodies crossing an external boundary at zero angle

The paper concerns an analysis of equilibrium problems for 2D elastic bodies with a thin Timoshenko inclusion crossing an external boundary at zero angle. The inclusion is assumed to be delaminated, thus forming a crack between the inclusion and the body. We consider elastic inclusions as well as rigid inclusions. To prevent a mutual penetration between the crack faces, inequality type boundary conditions are imposed at the crack faces.Theorems of existence and uniqueness are established. Passages to limits are investigated as a rigidity parameter of the elastic inclusion going to infinity.     

Photoelastic investigation of crack-inclusion interaction

An experimental procedure is presented for determining the mode I stress-intensity factor of an edge crack with a nearby rigid elliptical inclusion in a finite plate loaded in uniform tension. The rigid inclusion was modeled by bonding two identical steel inclusions on to the faces of a thin plate of polycarbonate. Models were constructed with edge cracks and various inclusion geometries so that the effect of parameters such as inclusion shape, orientation, and cracktip position on the stress-intensity factors of the crack could be determined. Photoelasticity experiments were used for this investigation and the results were compared to the results of a similar theoretical analysis of the interaction between a crack and an inclusion in an infinite plate. A good correlation was found between the experimental and theoretical models indicating that the results may help provide a better understanding of the toughening mechanisms in materials such as short-fiber-reinforced composites and ceramics. This experimental method is relatively easy to use making it an attractive candidate to be applied to similar problems involving cracks and inhomogeneities.     

Debonded arc-shaped interface conducting rigid line inclusions in piezoelectric composites

We consider an arc-shaped conducting rigid line inclusion located at the interface between a circular piezoelectric inhomogeneity and an unbounded piezoelectric matrix subjected to remote uniform anti-plane shear stresses and in-plane electric fields. Moreover, one side of the rigid line inclusion has become fully debonded from the matrix or the inhomogeneity leading to the formation of an insulating crack. After the introduction of two sectionally holomorphic vector functions, the problem is reduced to a vector Riemann–Hilbert problem, which can be decoupled sequentially by repeated application of the orthogonality relations between the eigenvectors for two corresponding generalized eigenvalue problems.     

Asymptotic behavior of the energy functional for a three-dimensional body with a rigid inclusion and a crack

A model of a three-dimensional elastic body containing a rigid inclusion and a crack located on the interface between the inclusion and the body is considered. Natural boundary conditions are imposed on the crack. A derivative of the energy functional with respect to the perturbation parameter is derived for an arbitrary, rather smooth perturbation of the domain, in particular, the Griffith formula is obtained.     

Interaction between a rigid inclusion and a line crack under uniform heat flux

The thermoelastic displacement boundary value problem for a rigid inclusion interacting with a line crack in an infinite plane subjected to a uniform heat flux is studied, in which the rigid body rotation of the inclusion is considered. To solve the prescribed problem, we use the principle of superposition to decompose it into two groups of problems, which are further reduced to several basic subproblems including Green’s functions of edge dislocation and heat source couple, as well as the problem of a plane containing the inclusion under uniform heat flux and the problem of the inclusion subjected to a small rotation. The problems are solved using the complex variable method along with the rational mapping function technique. The variations of the stress intensity factors at the crack tips and the rigid body rotation angles with various crack lengths and heat flux angles are shown. The effects of the inclusion shape and size are also investigated.     

Three-dimensional asymptotic stress field in the vicinity of the circumference of a penny shaped discontinuity

An eigenfunction expansion method is presented to obtain three-dimensional asymptotic stress fields in the vicinity of the front of a penny shaped discontinuity, e.g., crack, anticrack (infinitely rigid lamella), etc., subjected to the far-field torsion (mode III), extension/bending (mode I) and sliding shear/twisting (mode II) loadings. Five different discontinuity-surface boundary conditions are considered: (i) penny shaped crack, (ii) penny shaped anticrack or perfectly bonded thin rigid inclusion, (iii) penny shaped thin transversely rigid inclusion (frictionless planar slip permitted), (iv) penny shaped thin rigid inclusion in part perfectly bonded, the remainder with frictionless slip, and (v) penny shaped thin rigid inclusion alongside penny shaped crack. The computed stress singularity for a penny shaped anticrack is the same as that of the corresponding crack. The main difference is, however, that all the stress components at the circular tip of an anticrack depend on Poisson’s ratio under modes I and II.     

Wave momentum and scattering of elastic waves by two-dimensional thin objects

The elastic wave momentum equation is applied to scattering of dilatational and shear waves by two-dimensional thin objects. It is shown that the sources of wave momentum are located at the edges of these objects. For a stress-zero crack or for a rigid inclusion there are two sources at each edge, for a fluid-filled crack there is just one. The scattered wave is expressed in terms of these sources. This reduces the number of independent variables by a factor two. An application to inverse scattering problems is also given.     

P波对界面部分脱胶的刚性圆柱夹杂物的散射

本文求解了弹性Ｐ波对界面部分脱胶的可动刚性圆柱夹杂物的散射问题。将脱胶区看作表面不相接触的弧形界面裂纹，借助波函数展开法并利用边界条件将问题转化为一组对偶级数方程。然后通过引入裂纹面的位错密度函数，将其化为一组具有Ｈｉｌｂｅｒｔ核的第二类奇异积分方程，并进一步化为Ｃａｕｃｈｙ型奇异积分方程组，数值求解方程组可获得动应力强度因子，夹杂物刚体振动位移和散射截面等重要参量。结果显示该类结构在较低的频率上发生共振，此低频共振特性与脱胶区大小，入射波方向、材料组合等多种参数有关。与已有方法相比，本文的方法更具一般性，适用于任意材料组合。     

Stress state of a piezoceramic body with a plane crack opened by a rigid inclusion

The paper establishes a relationship between the solutions for cracks located in the isotropy plane of a transversely isotropic piezoceramic medium and opened (without friction) by rigid inclusions and the solutions for cracks in a purely elastic medium. This makes it possible to calculate the stress intensity factor (SIF) for cracks in an electroelastic medium from the SIF for an elastic isotropic material, without the need to solve the electroelastic problem. The use of the approach is exemplified by a penny-shaped crack opened by either a disk-shaped rigid inclusion of constant thickness or a rigid oblate spheroidal inclusion in an electroelastic medium __________ Translated from Prikladnaya Mekhanika, Vol. 44, No. 7, pp. 47–60, July 2008.     

The in-plane loading of rigid disc inclusion embedded in a crack

The paper examines the in-plane loading of a disc shaped rigid disc inclusion which is embedded in bonded contact with the plane surfaces of a penny-shaped crack. The mixed boundary value problem governing the elastostatic problem is reduced to the solution of a system of coupled integral equations, which are solved numerically to determine results of engineering interest. These results include the in-plane stiffness of the disc inclusion and the crack opening mode stress intensity factor at the boundary of the penny-shaped crack.     

Internally indented penny-shaped cracks: A comparison of analytical and boundary integral equation estimates

The problem of the axisymmetric internal indentation of a penny-shaped crack by a rigid circular inclusion is discussed. The paper presents a comparison of analytical and boundary integral equation results for the stress intensity factor at the boundary of the penny-shaped crack indented by a smooth inclusion. Numerical results presented in the paper examines the influence of features such as adhesion at the inclusion-elastic medium interface and finite geometry of the elastic solid containing the penny shaped crack.     

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.     

Problem of interaction of thin rigid needle-shaped inclusions located between different elastic materials

We study a piecewise-homogeneous elastic plane composed of two half-planes with different elastic parameters and two thin rigid needle-shaped inclusions located between them. One inclusion is rigidly connected with the environment, and the other inclusion is not, while contacting with it like a smooth rigid punch. We consider the plane deformed state generated by stresses given at infinity. The problem is reduced to a combination of a matrix Riemann boundary-value problem from the theory of analytic functions and a matrix Hilbert problem, which can be solved in terms of integrals through the reduction to two separate scalar Riemann boundary-value problems on a twosheeted Riemann surface.We explicitly obtain the complex potentials of the composite elastic plane, the stress intensity factors near the tips of the inclusion, and the rotation angles of the inclusions. We also present numerical examples illustrating how the stresses near the inclusions depend on the elastic and geometric parameters of the problem.     

The singular elastostatic field due to a crack in rubberlike materials

Within the framework of finite-strain elastostatics an asymptotic analysis is carried out in order to calculate the singular field near the crack tip in a slab under conditions of plane deformation. A class of Ogden-Ball hyperelastic rubberlike materials and general loading conditions ensuring vanishing tractions on the crack faces near the crack tip are considered. It is shown that the singular deformation field near the crack tip can be specified by applying a rigid body rotation with a subsequent parallel translation to a so-called canonical field. The adjective canonical is adopted here to denote the field with symmetrically opening crack faces, just resembling the displacement field of the symmetric mode in linear elastic fracture mechanics. No analogy with the antisymmetric mode is possible, and the crack equilibrium criterion requires only one stress intensity factor to be determined.     

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.     

A Westergaard-type stress function for line inclusion problems

A direct correspondence is shown to exist between the Westergaard stress function for crack problems and a newly-introduced stress function for rigid line inclusion problems. A correspondence is also shown to exist between stress intensities in the crack and the inclusion problems, as well as between the opening displacement of the crack and the axial force in the inclusion. A scheme is presented to modify these rigid-inclusion solutions to account approximately for non-zero compliance of real fibers in a composite material.     

Interaction between annular crack and rigid disc inclusion in a transversely isotropic solid

An analytical approach has been presented to analyze the asymmetric and axisymmetric interactions between an annular crack and a rigid disc inclusion embedded in a transversely isotropic full-space. With the aid of a method of potential functions, Hankel and Abel transforms, the solution of the problems is reduced to a system of Fredholm integral equations, which are solved by using a numerical method. In each case, the stiffness of the disc inclusion and the stress intensity factor at the tips of the annular crack for different degrees of material anisotropy and different ratios of the inner and outer radius of the crack are illustrated graphically. Several limiting cases such as penny-shaped crack and external crack along with some exact solutions are presented to demonstrate the efficiency of the method.     

压电材料反平面应变状态的任意形状夹杂问题

应用复函数的Ｆａｂｅｒ级数展开方法，分析了含任意形状夹杂的压电材料反平面应变问题，给出了问题的复势函数解。利用这个解，具体讨论了椭圆形夹杂及其极限（几何方面与物理方面）问题。并给出了三角形、正方形夹杂的近似结果。其特例结果与早期工作一致     

Role of inclusion stiffness and interfacial strength on dynamic matrix crack growth: An experimental study

Experimental simulations of dynamic crack growth past inclusions of two different elastic moduli, stiff (glass) and compliant (polyurethane) relative to the matrix (epoxy), are carried out in a 2D setting. Full-field surface deformations are mapped in the crack–inclusion vicinity optically. The crack growth behavior as a function of inclusion–matrix interfacial strength and the inclusion location relative to the crack is studied under stress-wave loading conditions. An ultra high-speed rotating mirror-type digital camera is used to record random speckle patterns in the crack–inclusion vicinity to quantify in-plane displacement fields. The crack-tip deformation histories from the time of impact until complete fracture are mapped and fracture parameters are extracted. The crack front is arrested by the symmetrically located compliant inclusion for about half the duration needed for complete fracture event. The dynamically propagating crack is attracted and trapped by the weakly bonded inclusion interface for both stiff and compliant symmetrically located inclusion cases, whereas it is deflected away by the strongly bonded stiff inclusion and attracted by strongly bonded compliant inclusion when located eccentrically. The crack is arrested by a strongly bonded compliant inclusion for a significant fraction of the total dynamic event and is longer than the one for the weakly bonded counterpart. The compliant inclusion cases show higher fracture toughness than the stiff inclusion cases. Measured crack-tip mode-mixities correlate well with the observed crack attraction and repulsion mechanisms. Macroscopic examination of fracture surfaces reveals much higher surface roughness and ruggedness after crack–inclusion interaction for compliant inclusion than the stiff one. Implications of these observations on the dynamic fracture behavior of micron size A-glass and polyamide (PA6) particle filled epoxy is demonstrated. Filled-epoxy with 3% V_f of PA6 filler is shown to produce the same dynamic fracture toughness enhancement as the one due to 10% V_f glass.