The stress concentration near a rigid line inclusion in a prestressed, elastic material. Part II.: Implications on shear band nucleation, growth and energy release rate

The full-field and asymptotic solutions derived in Part I of this article (for a lamellar rigid inclusion, embedded in a uniformly prestressed, incompressible and orthotropic elastic sheet, subject to a far-field deformation increment) are employed to analyse shear band formation, as promoted by the near-tip stress singularity. Since these solutions involve the prestress as a parameter, stress and deformation fields can be investigated near the boundary of ellipticity loss (but still within the elliptic range). In the vicinity of this boundary, the incremental stress and displacement fields evidence localized deformations with patterns organized into shear bands, evidencing inclinations corresponding to those predicted at ellipticity loss. These localized deformation patterns are shown to explain experimental results on highly deformed soft materials containing thin, stiff inclusions. Finally, the incremental energy release rate and incremental J-integral are derived, related to a reduction (or growth) of the stiffener. It is shown that this is always positive (or negative), but tends to zero approaching the Ellipticity boundary, which implies that reduction of the lamellar inclusion dies out and, simultaneously, shear bands develop.     

Effects of pre-stress on crack-tip fields in elastic,incompressible solids

A closed-form asymptotic solution is provided for velocity fields and the nominal stress rates near the tip of a stationary crack in a homogeneously pre-stressed configuration of a nonlinear elastic, incompressible material. In particular, a biaxial pre-stress is assumed with stress axes parallel and orthogonal to the crack faces. Two boundary conditions are considered on the crack faces, namely a constant pressure or a constant dead loading, both preserving an homogeneous ground state. Starting from this configuration, small superimposed Mode I or Mode II deformations are solved, in the framework of Biot's incremental theory of elasticity. In this way a definition of an incremental stress intensity factor is introduced, slightly different for pressure or dead loading conditions on crack faces. Specific examples are finally developed for various hyperelastic materials, including the J₂-deformation theory of plasticity. The presence of pre-stress is shown to strongly influence the angular variation of the asymptotic crack-tip fields, even if the nominal stress rate displays a square root singularity as in the infinitesimal theory. Relationships between the solution with shear band formation at the crack tip and instability of the crack surfaces are given in evidence.     

Finite strains at the tip of a crack in a sheet of hyperelastic material: III. General bimaterial case

In this last in a series of three papers, we summarize an asymptotic analysis of the near-tip stress and deformation fields for an interface crack between two sheets of Generalized Neo-Hookean materials. This investigation, which is consistent with the nonlinear elastostatic theory of plane stress, allows for an arbitrary choice, on both sides of the three parameters characterizing this class of hyperelastic materials. The first three terms of the approximation series are obtained, showing the existence of a non-oscillatory and contact-free solution to the interface crack problem. The analytical results are compared with a full-field solution obtained numerically using the finite element method.     

Dynamics of a prestressed stiff layer on an elastic half space: filtering and band gap characteristics of periodic structural models derived from long-wave asymptotics

Anti-plane and plane-strain, time-harmonic, small-amplitude vibrations of an elastic layer on an elastic half space are considered, superimposed upon a state of finite, uniform stress and strain. A (compressible) elastic material is considered, orthotropic with orthotropy axes aligned parallel and orthogonal both to the layer and the prestress principal directions. A non-uniform mass density is assumed in the layer. A formal long-wave asymptotic solution is derived under the assumptions of high contrast between the stiffnesses of the layer and the half space and between certain prestress components and the current elastic shear modulus.It is shown that (i) the layer asymptotically behaves as a beam subject to transversal and axial vibrations; (ii) the response of the half space can be found in a closed-form, under the assumption of plane wave motion (which becomes consistent when the density of the layer is uniform), otherwise it is represented by a hypersingular integral equation; (iii) if the nonlocality introduced by the hypersingular integral equation is restricted to an influence area of finite extent, the integral can be analytically approximated, so that a Winkler-type spring model representing the half space is rigorously derived. For uniform density of the layer, the constants defining the spring model are given as functions of the prestress and anisotropy parameters of the half space; and, finally, (iv) the asymptotic solution provides new analytical expressions for incremental displacement of the layer, which, compared to the exact numerical solution (also included), are shown to perform quite well, even for values of parameters much beyond the limits imposed by the asymptotic analysis.The asymptotic analysis allows us to explore, for the first time, dynamic properties of a periodic layer bonded to an elastic half space and subject to a uniform prestress state. We find that the system exhibits band gaps (ranges of forbidden frequencies) and that the prestress can be used as a parameter tuning the filtering properties of the structure, an effect which may have important consequences in the design of resonant devices.     

Asymptotic fields around an interfacial crack with a cohesive zone ahead of the crack tip

This paper considers an interfacial crack with a cohesive zone ahead of the crack tip in a linearly elastic isotropic bi-material and derives the mixed-mode asymptotic stress and displacement fields around the crack and cohesive zone under plane deformation conditions (plane stress or plane strain). The field solution is obtained using elliptic coordinates and complex functions and can be represented in terms of a complete set of complex eigenfunction terms. The imaginary portion of the eigenvalues is characterized by a bi-material mismatch parameter ε = arctanh(β)/π, where β is a Dundurs parameter, and the resulting fields do not contain stress singularity. The behaviors of “Mode I” type and “Mode II” type fields based on dominant eigenfunction terms are discussed in detail. For completeness, the counterpart for the Mode III solution is included in an appendix.     

Three-dimensional stress analysis of textile composites: Part I. Numerical analysis

Three-dimensional stress analysis in a unit-cell of a plain-woven composite was performed by using B-spline displacement approximation. The spline approximation provides continuity of displacement and stress components within each yarn and matrix subregion. Two types of unit-cell problems with and without inter-yarn delamination were considered. A penalty function approach along with a contact surface characteristic function was used to obtain a full-field numerical solution for the frictionless contact problem between delaminated yarn surfaces.Yarn interfaces at yarn-crossover locations represent three-material wedge-type regions resulting in singular stress behavior. In the case of unit-cells with perfect bonding between the yarn interfaces, the numerical values of the inter-yarn normal stress did not exhibit trends typical for unbounded stress behavior, whereas the inter-yarn shear stress components displayed discontinuous behavior typical for numerical results in the vicinity of the stress singularity. In the presence of the delamination, both the inter-yarn normal and shear stress components exhibited unbounded behavior near the singularity. Notably, the inter-yarn normal stress showed signs of singular behavior in both cases of open and closed delaminations. Due to the stress singularity that exists at yarn-crossover locations containing three materials (yarn–yarn–matrix) interface intersections, the full-field numerical solution, even with high-order approximation functions, was not able to capture the directional nonuniqueness of the stress values in the vicinity of the singularity, and therefore calls for incorporation of the asymptotic singular stress analysis, which will be given in a follow-on paper [Sihn and Roy, International Journal of Solids and Structures (accepted for publication)].     

An analytical and experimental stress analysis of a practical mode II fracture test specimen

A boundary-collocation method has been employed to determine the Mode II stress-intensity factors for a pair of through-the-thickness edge cracks in a finite isotropic plate. An elastostatic analysis has been carried out in terms of the complete Williams stress function employing both even and old components. The results of the numerical analysis were verified by a two-step procedure whereby the symmetric (Mode I) and antisymmetric (Mode II) portions of the solution were independently compared with existing solutions. Since no previous analytical solutions existed for the asymmetric loading of an edge-cracked plate, the complete solution was verified by comparison with a photoelastic analysis. A compact shear (CS) specimen of Hysol epoxy resin was loaded in a photoelastic experiment designed to study the isochromatic-fringe patterns resulting from the Mode II crack-tip stress distribution. The experiment verified that a pure mode II stress distribution existed in the neighborhood of the crack tips, and confirmed the accuracy of the boundary-collocation solution for the Mode II stress-intensity factors. Specimen center-line stress-distribution data were obtained photoelastically and employed to refine the boundary-collocation analysis. Agreement between the analytically and experimentally determined Mode II stress-intensity factors was excellent.     

Mode II dynamic fields near a crack tip growing in an elastic-perfectly-plastic solid

An asymptotic solution is given for Mode II dynamic fields in the neighborhood of the tip of a steadily advancing crack in an incompressible elastic—perfectly-plastic solid (plane strain). It is shown that, like for Modes I and III (Gao and Nemat-Nasser, 1983), the complete dynamic solution for Mode II predicts a logarithmic singularity for the strain field, but unlike for those modes which involve no elastic unloading, the pure Mode II solution includes two elastic sectors next to the stress-free crack surfaces. This is in contradiction to the quasi-static solution which predicts a small central plastic zone, followed by two large elastic zones, and then two very small plastic zones adjacent to the stress-free crack faces. The stress field for the complete dynamic solution varies throughout the entire crack tip neighborhood, admitting finite jumps at two shock fronts within the central plastic sector. This dynamic stress field is consistent with that of the stationary crack solution, and indeed reduces to it as the crack growth speed becomes zero.     

High-order asymptotics and perturbation problems for 3D interfacial cracks

We present an asymptotic algorithm for analysis of a singularly perturbed problem in a domain containing an interfacial crack. The crack is assumed to be flat and its front, initially straight, is perturbed in the plane containing the crack. The aim of the work is to determine the asymptotic representation of the stress-intensity factors near the edge of the crack. Mathematically, the limit problem is reduced to the analysis of a matrix, 3×3, Wiener-Hopf problem, and its solution generates the “weight matrix-function” characterised by a special singular solution near the crack edge. The two-term asymptotic representation for the weight function components is required by the asymptotic algorithm, together with two-term asymptotics for stress components associated with the physical fields near the edge of the crack. The particular feature of the solution is the coupling between the normal opening mode (Mode-I), and the shear modes (Mode-II and Mode-III), and the oscillatory behaviour of certain stress components near the crack edge. Explicit asymptotic formulae for the stress-intensity factors are obtained at the edge of a “wavy crack” at an interface.     

On the lamb solution and Rayleigh-wave-induced cracking

This paper examines the extension of surface microcracks induced by a surface or Rayleigh wave (R-wave). This problem is examined both theoretically and experimentally. The theoretical approach involves a full-field reappraisal of the Lamb solution for a surface wave propagating in a homogeneous, isotropic, elastic, two-dimensional material for the cases of plane strain and plane stress. Using the Griffith-Irwin energy-release-rate fracture criterion for cracks under combined Mode I and Mode II loading, a prediction is made of the path and final length of the surface microcrack extension produced by the R-wave. Predictions of the crack-extension direction are also obtained using the maximum normal-stress fracture criterion. The experimental approach uses dynamic photoelasticity to observe the isochromatic patterns associated with an R-wave propagating along the narrow edge of a transparent birefringent plate, examining in detail the process of crack extension. When the theoretically and experimentally obtained results are compared, reasonable agreement is obtained.     

Stress-strain state in the vicinity of a crack tip under mixed loading

A method is proposed to calculate the eigenvalues of the class of nonlinear eigenvalue problems resulting from the problem of determining the stress-strain state in the vicinity of a crack tip in power-law materials over the entire range of mixed modes of deformation, from the opening mode to pure shear. The proposed approach was used to found eigenvalues of the problem that differ from the well-known eigenvalue corresponding to the Hutchinson-Rice-Rosengren solution. The resulting asymptotic form of the stress field is a self-similar intermediate asymptotic solution of the problem of a crack in a damaged medium under mixed loading. Using the new asymptotic form of the stress field and introducing a self-similar variable, we obtained an asymptotic solution of the problem of a crack in a damaged medium and constructed the regions of dispersed material near the crack.     

On the effects of characteristic lengths in bending and torsion on Mode III crack in couple stress elasticity

The problem of a stationary semi-infinite crack in an elastic solid with microstructures subject to remote classical K_III field is investigated in the present work. The material behavior is described by the indeterminate theory of couple stress elasticity developed by Koiter. This constitutive model includes the characteristic lengths in bending and torsion and thus it is able to account for the underlying microstructure of the material as well as for the strong size effects arising at small scales. The stress and displacement fields turn out to be strongly influenced by the ratio between the characteristic lengths. Moreover, the symmetric stress field turns out to be finite at the crack tip, whereas the skew-symmetric stress field displays a strong singularity. Ahead of the crack tip within a zone smaller than the characteristic length in torsion, the total shear stress and reduced tractions occur with the opposite sign with respect to the classical LEFM solution, due to the relative rotation of the microstructural particles currently at the crack tip. The asymptotic fields dominate within this zone, which however has limited physical relevance and becomes vanishing small for a characteristic length in torsion of zero. In this limiting case the full-field solution recovers the classical K_III field with square-root stress singularity. Outside the zone where the total shear stress is negative, the full-field solution exhibits a bounded maximum for the total shear stress ahead of the crack tip, whose magnitude can be adopted as a measure of the critical stress level for crack advancing. The corresponding fracture criterion defines a critical stress intensity factor, which increases with the characteristic length in torsion. Moreover, the occurrence of a sharp crack profile denotes that the crack becomes stiffer with respect to the classical elastic response, thus revealing that the presence of microstructures may shield the crack tip from fracture.     

Regular Perturbation Methods for a Region with a Crack

The paper considers a model problem for Poisson's equation for a region containing a crack or a set of cracks under arbitrary linear perturbation. Variational formulation of the problem using smooth mapping of regions yields a complete asymptotic expansion of the solution in the perturbation parameter, which is a generalized shape derivative. This global asymptotic expansion of the solution was used to derive representations of arbitraryorder derivatives for the potential energy function, stress intensity factors, and invariant energy integrals in general form and for basis perturbations of the region (shear, tension, and rotation). The problem of the local growth of a branching crack for the Griffith fracture criterion and the linearized problem of optimal location of a rectilinear crack in a body with the energy function as a cost functional were formulated.     

On the diffusion of load from a stiffener into an infinite wedge-shaped plate

Summary The stress-distribution in a wedge-shaped plate with a stiffener upon one of the edges is considered. The stiffener is loaded by an axial force. The problem leads to the solution of a biharmonic equation with one mixed boundary condition. The problem is reduced to the standard problem of the stress-distribution in a wedge. The reduction has been executed by the solution of a difference equation for the transform of the shear-stress along the stiffened edge. For this solution we give two representations: one by means of an infinite product and one by means of an integral. Full discussion is given on asymptotic behaviour and on the numerical aspects.     

水下爆炸载荷作用下加筋板的毁伤模式

为研究加筋板结构在水下爆炸载荷作用下的破坏模式,进行了一系列模型试验。通过对试验结果的分析,将加筋板的毁伤模式总体分为塑性大变形（模式I）、拉伸破坏（模式II）和剪切破坏（模式III）等3大类。又根据载荷强弱和加强筋的相对强度,将每种破坏模式细分为3种子模式。提出了一个量纲一的损伤因子,作为预报毁伤模式的判据。     

A MODE Ⅱ CRACK IN A TWO-DIMENSIONAL OCTAGONAL QUASICRYSTALS

The present study develops the fracture theory for a two-dimensional octagonal quasicrystals. The exact analytic solution of a Mode Ⅱ Griffith crack in the material was obtained by using the Fourier transform and dual integral equations theory, then the displacement and stress fields, stress intensity factor and strain energy release rate were determined, the physical sense of the results relative to phason and the difference between mechanical behaviors of the crack problem in crystal and quasicrystal were figured out. These provide important information for studying the deformation and fracture of the new solid phase.     

Fracture mechanics of plates and shells applied to fail-safe analysis of fuselage Part I: Theory

In this paper, a general theory on the asymptotic field near the crack tip for plates and shells with and without shear deformation effect is established. It is found that four stress intensity factors, two for symmetrical and antisymmetrical stretching and two for symmetrical and antisymmetrical bending, are required to describe arbitrary asymptotic fields near the crack tip for plates without shear deformation. An additional stress intensity factor is required for the transverse shearing force induced by antisymmetrical bending when the shear deformation is included in the analysis. It is also proven by means of the complex variable technique that for problems of plates with shear deformation, there exist similarities in the asymptotic expressions of moments and membrane forces and also in the asymptotic expressions of in-plane displacements and rotations of the mid-surface. The energy release rate associated with crack growth in the direction of the crack line can be expressed in terms of stress intensity factors by means of Irwin's method of work and energy associated with a virtual crack extension. A combined stress intensity factor can be defined through the total energy release rate. The theory of the fracture of plates is generalized and applied to the study of problems in the fracture of shells. An example of an infinitely long cylindrical shell with a circumferential crack subjected to remote axial tension is given to demonstrate the application of the theory and to test the accuracy of the numerical analysis used for the problem.     

Finite strains at the tip of a crack in a sheet of hyperelastic material: II. Special bimaterial cases

An asymptotic analysis of the strain and stress near-tip fields for a crack in a sheet of Generalized Neo-Hookean materials is presented in this second in a series of three papers. The analysis is based on the nonlinear plane stress theory of elasticity and concerns two special cases of the interface crack problem: in the first situation both components have the same hardening behavior; next, we investigate the particular case of a sheet of Generalized Neo-Hookean material bonded to a rigid substrate. The transition between the two special cases is studied in detail. The analytical results are also compared with a full-field finite element solution.     

Analytical model of thermal striping for a micro-cracked solid

Using the formal asymptotic approximation of the Mode I stress intensity factor for an edge crack in a thermoelastic half plane containing several small voids obtained in [Nieves, M.J., Movchan, A.B., and Jones, I.S., 2011. Asymptotic study of a thermoelastic problem in a semi-infinite body containing a surface-breaking crack and small perforations. QJMAM 64 (3), 349–369] we investigate the effect of micro-cracks on this stress intensity factor. In numerical examples, we show how the behaviour of the stress intensity factor as a function of crack depth is affected by micro-cracks of different orientations occurring in the half space.     

Experimental investigation of viscoelastic drop deformation in Newtonian matrix at high capillary number under simple shear flow

The steady state deformation of a viscoelastic drop (Boger fluid) in a Newtonian liquid at high capillary number under simple shear flow is investigated by direct visualization using a specially designed Couette apparatus which enables visualization from two perpendicular directions. Two drop deformation modes are found: (1) Mode I – drop deformation in the flow direction and (2) Mode II – drop deformation in the vorticity direction. The drop deformation mode depends on the relative strength of the elastic contribution to viscous contribution. If the elastic contribution is weak compared to the viscous contribution, the drop elongates in the flow direction via Mode I. If the elastic contribution is strong, the drop elongates in the vorticity direction via Mode II. The drop size also affects the drop deformation. At the same capillary number, bigger drops have larger deformations than smaller drops.