Plane-stress crack-tip fields for power-law hardening orthotropic materials

Plane stress mode I near-tip fields in orthotropic materials are examined. Plastic orthotropy is described by Hill's quadratic yield function and the strain hardening behavior is given by an appropriate generalization of a uniaxial tensile power-law stress-strain relation. Pronounced changes in the pattern of the angular variations of crack-tip fields have been observed with the degree of plastic orthotropy and the amount of strain hardening. Possible shapes and sizes of plastic zones (as inferred from effective stress contours) are presented for high- and low-hardening materials and a wide range of plastic orthotropy. The shape of the plastic zone for a particular case of plastic orthotropy agreed remarkably well with the zone of intense straining induced by an appropriately orientated crack within a graphite/epoxy laminate.

---

Résumé On examine les champs de contraintes planes selon un mode I au voisinage de l'extrémité d'une fissure dans des matériaux orthotropes. L'orthotrope plastique est décrite par la fonction quadratique de plastification de Hill, et le comportement à l'écrouissage est donné par une généralisation adéquate d'une relation tensioncilatation de forme parabolique, sous traction mono-axiale. On a observé des modifications profondes dans l'aspect des variations angulaires des champs d'extrémité de fissure, selon le degré d'orthotropie plastique et infensité de l'écrouissage. Pour des matériaux très sujets ou peu sujets à l'écrouissage, et pour une large gamme d'orthotropies plastiques, on présente les formes et dimensions possibles des zones plastiques, telles qu'elles se deduisent des contours effectifs de contraintes. La forme de la zone plastique correspondant au cas particulier d'une orthotropie plastique s'accorde remarquablement bien à la zone de dilatation importante créée par une fissure d'orientation appropriée, dans une plaque de graphite-epoxy.

     

Application of an orthotropy rescaling method to edge cracks and kinked cracks in orthotropic two-dimensional solids

Oblique edge cracks and kinked cracks in orthotropic materials with inclined principal material directions under inplane loadings are investigated. The Stroh formalism is modified by introducing new complex functions, which recovers a classical solution for a degenerate orthotropic material with multiple characteristic roots. An orthotropy rescaling technique is presented based on the modified Stroh formalism. Stress intensity factors for edge cracks as well as kinked cracks are obtained in terms of solutions for a material with cubic symmetry by applying the orthotropy rescaling method. Explicit expressions of the stress intensity factors for a degenerate orthotropic material are obtained in terms of solutions for an isotropic material. The effects of orthotropic parameter, material orientation, and crack angle on the stress intensity factors for the degenerate orthotropic material are discussed. The stress intensity factors for cubic symmetry materials are calculated from finite element analyses, which can be used to evaluate the stress intensity factors for orthotropic materials. The energy release rate for the kinked crack in an orthotropic material is also obtained.     

Crack growth in an orthotropic medium

Abstract

The elastostatic problem of an orthotropic body having a central inclined crack and subjected at infinity to a uniform biaxial load is considered. It is assumed that the crack line does not coincide with an axis of elastic symmetry of the body. The problem must be considered as one of general orthotropy, due in particular to the fact that the elastic coefficients of the material change with rotation of the reference system. The stress fields at the crack tip are reported and the presence of the non‐singular terms underlined. The Strain Energy Density Theory is extended to orthotropic materials. It is assumed that the Critical Strain Energy Density Factor has a polar variation. The crack initiation is determined via minimization of the ratio of the strain energy density over the material critical strain energy density, pointing out the effects of orthotropy and load biaxiality. The effects of the non‐singular terms on crack growth for different orthotropic materials is also studied, underling the relation between orthotropy and non‐singular terms.     

Plastic limit analysis of cylindrically orthotropic circular plates

The plastic limit analysis of cylindrically orthotropic circular plates is developed using a piecewise linear orthotropic yield criterion. The yield criterion is a modification of an isotropic formulation that consists of a series of weighted piecewise linear components. The piecewise linear yield criterion enables an analytical solution for the plastic limit load of cylindrically orthotropic circular plates. Plastic limit analysis for both simply supported and clamped circular plates under uniformly distributed load are carried out. Parametric studies are conducted to investigate the sensitivity of the plastic limit loads to material orthotropy and influences of orthotropic ratio and chosen yield criteria on the plastic limit loads of the circular plates are discussed. It is found that the plastic limit loads of the orthotropic circular plates are affected significantly by the orthotropic ratio. Enhancement of the circumferential yield moment will increase dramatically the plastic limit load of the plates. Moment and velocity fields of the plates in plastic limit state are also derived and discussed. The results obtained from the present study are helpful in understanding the failure mechanism of orthotropic circular plates and is useful for design.     

The asymptotic structure of small-scale yielding interfacial free-edge joint and crack-tip fields

Summary The problem of the small-scale yielding (SSY) plane-strain asymptotic fields for the interfacial free-edge joint singularity is examined in detail, and comparisons are made with the interfacial crack tip. The geometries are idealized as isotropic elasto-plastic materials with Ramberg-Osgood power-law hardening properties bonded to a rigid elastic substrate. The resulting fields are shown to be singular and are presented in terms of radial and angular distributions of stress and displacement, and as idealized plastic slip-line sectors. A fourth-order Runge-Kutta numerical method provides solutions to fundamental equations of equilibrium and compatibility that are verified with those of a highly focused finite element (FE) analysis. It is shown that, as in the case of the crack, the asymptotic singular fields are only dependent on the hardening parameter and only a small range of interfacial mode-mix ratios are permitted. The order for the stress singularity may be formulated in terms of the hardening parameter and the elastic solution for incompressible material. The rigid-slip-line field for the interfacial free-edge joint is presented, and it is shown that there is some significant similarity between the asymptotic fields of the deviatoric polar stresses for the joint and the crack-tip having an elastic wedge sector.     

Influence of anisotropy on the limit load of bi-material welded cracked joints subject to tension

A plane-strain upper bound limit load solution for bi-material welded joints subject tension is found. It is assumed that each material obeys Hill's orthotropic yield criterion and one of the principal axes of orthotropy coincides with the tensile direction. A crack of arbitrary length is located at the bi-material interface. The solution is based on a simple discontinuous kinematically admissible velocity field and is an extension of the corresponding solution for the specimen made of isotropic materials. The main purpose of the paper is to demonstrate that the influence of anisotropy on the magnitude of limit loads may be much more significant than other effects.     

Plane-stress crack-tip fields for perfectly plastic orthotropic materials

Plane stress mode I crack-tip fields for perfectly plastic orthotropic materials are studied. Plastic orthotropy is described by Hill's quadratic yield function. The construction of crack-tip fields is based on the general crack-tip field analysis for elastic perfectly plastic materials given by Rice [1] and guided by the corresponding low-hardening power-law solutions. Two very different types of plane-stress crack-tip fields emerge as plastic orthotropy is varied. The first one consists of a centered fan sector in front of the crack tip and two neighboring constant stress sectors. The second one consists of a constant stress sector in front of the crack tip, a constant stress sector bordering the crack face, and a centered fan sector between the two constant stress sectors. All the perfectly plastic crack-tip solutons are verified by the corresponding low-hardening power-law solutions. General trends of crack-tip stress solutions as functions of plastic orthotropy and implications of these solutions to the design of ductile composite materials are discussed.

---

Résumé On étudie les champs de contraintes planes de mode I à l'extrémité d'une fissure, dans les matériaux orthotropiques parfaitement plastique. L'orthotropie plastique est décrite par la fonction quadratique de plasticité de Hill. On base les constructions des champs de constraintes sur l'analyse générale des constraintes à l'extrémité d'une fissure fournie par Rice pour les matériaux élastiques parfaitement plastiques, que l'on règle par les lois paraboliques caractérisant un faible écrouissage. Lorsque l'on modifie l'orthotropie plastique, il apparaît deux types de champs de contraintes à l'extrémité de la fissure très différents. Le premier comporte un secteur en éventail centré sur le front de fissure, et deux secteurs voisins à contraintes constantes. Le second consiste en une secteur à contrainte au bord de la surface de la fissure, et un secteur en éventail centré sur les deux secteurs à contraintes constantes. Toutes les solutions relatives à une extrémité de fissures parfaitement plastique sont vérifiées par les fonctions paraboliques d'écrouissage faible correspondantes. On discute des tendances générales que suivent les solutions pour les contraintes en extrémité de fissure selon l'orthrotropie plastique, et des implications que comportent ces solutions dans la conception de matériaux composites ductiles.

     

On the effect of orthotropy in a cracked cylindrical shell

A pressurized cylindrical shell containing a longitudinal crack is considered. It is assumed that the material has a special orthotropy, namely that the shear modulus of the sheet may be evaluated from the measured Young's moduli and the Poisson's ratios rather than being an independent material constant. Two examples, one for a mildly orthotropic (titanium) and the other for a strongly orthotropic (graphite) material approximately satisfying the condition of special orthotropy are given. The results show that the stress intensity factors are rather strongly dependent on the degree of orthotropy.     

Verification of a non-local stress criterion for mixed mode fracture in wood

An extension of a non-local stress fracture criterion to orthotropic materials based on the damage model of an elastic solid containing growing microcracks was presented in this paper. By taking this approach, a new fracture condition expressed in terms of the mixed mode stress intensity factors for orthotropic materials was proposed and its applicability to predict of a crack initiation and propagation in wood was validated. Predicted values of the stress intensity factors at failure were compared to experimental observations carried out on wood specimens for cracks arbitrarily oriented with respect to the orthotropy axes. Special considerations were applied to the comparison of the non-local stress fracture criterion with some classical fracture criteria for orthotropic materials.     

Stresses in an orthotropic strip containing A Griffith Crack

The problem of determining the stress and displacement fields in an orthotropic elastic strip containing a Griffith crack situated symmetrically and oriented in a direction normal to the edges of the strip is considered. A general solution in terms of two potential functions is presented. The mixed boundary conditions lead to dual integral equations, which are reduced to Fredholm integral equation of second kind and are solved by the use of Gaussian quadrature formula. Numerical solutions for a fiber-reinforced composite material and some isotropic materials are carried out and the effect of orthotropy on various quantities of physical interest, in fracture mechanics, is discussed.     

Effect of orthotropic material on finite element modeling of completely dentate mandible

The objective of this investigation is to construct a high quality complete dentate mandible model with detailed biological structures, and assign mandibular bone with inherent orthotropic material characteristics. Three different types of scan data are used to elaborate detailed mandibular structures, including the cortical and cancellous bone, tooth enamel, dentin, periodontal ligament, temporal fossa, TMJ articular disk, temporal cartilage, and condylar cartilage. In addition, an extended orthotropic material assignment methodology based on harmonic fields is used to handle the alveolar ridge region of dentate mandible, to generate compatible orthotropic axes fields. The influence of orthotropic material on the biomechanical behavior of complete dentate mandible is analyzed compared with commonly used isotropic model. The result revealed that the orthotropic model would induce higher stress values and more well-distributed stress pattern than the isotropic model, especially for the cancellous bone. And the orthotropic model would induce lower volumetric strain values than the isotropic model on the cortical bone. It was concluded that elastic orthotropy had a significant effect on the simulated stress value and distribution pattern, as well as the volumetric strain, and demonstrated the mechanical optimality of the mandible.     

Deep penetration analysis with dynamic cylindrical cavitation fields

Dynamic cylindrical cavitation fields are studied for a family of plastic orthotropic solids with arbitrary strain hardening response. Analysis is within the framework of plane-strain, steady state flow theory of associated plasticity. New formulae for cavitation pressure are validated against accurate numerical analysis and contact is made with existing studies. A uniform procedure is presented for estimating penetration depth of rigid axisymmetric projectiles at normal impact. Comparison with available experimental data reveals a very good agreement for both spherical and cylindrical dynamic cavitation models. Quasi-static cavitation pressure formulae can predict penetration depth with an appropriate scaling of the yield stress. The scaling factors appear to be independent of material properties but reflect the shape of head profile.     

Quasi-statically growing crack-tip fields in elastic perfectly plastic pressure-sensitive materials under plane strain conditions

Quasi-statically growing crack-tip fields in elastic perfectly plastic pressure-sensitive materials under plane strain conditions are investigated in this paper. The materials are assumed to follow the Drucker-Prager yield criterion and the normality flow rule. The asymptotic mode I crack-tip fields are assumed to follow the five-sector assembly of Drugan et al. (1982) for Mises materials. The crack-tip sectors, in turns, from the front of the crack tip are a constant stress sector, a centered fan sector, a non-singular plastic sector, an elastic sector and finally a trailing non-singular plastic sector bordering the crack face. The results of the asymptotic analysis show that as the pressure sensitivity increases, the plastic deformation shifts to the front of the tip, the angular span of the elastic unloading sector increases, and the angular span of the trailing non-singular plastic sector bordering the crack surface decreases. As the pressure sensitivity increases to about 0.6, the angular span of the trailing non-singular plastic sector almost vanishes. The effects of the border conditions between the centered fan sector and the first non-singular plastic sector on the solutions of the crack-tip fields for both Mises and pressure-sensitive materials are investigated in details. This revised version was published online in August 2006 with corrections to the Cover Date.     

Cracked shells under skew-symmetric loading

In this paper, the general problem of a shell containing a through crack in one of the principal planes of curvature and under general skew-symmetric loading is considered. By employing a Reissner type shell theory which takes into account the effect of transverse shear strains, all boundary conditions on the crack surfaces are satisfied separately. Consequently, unlike those obtained from the classical shell theory, the angular distributions of the stress components around the crack tips are shown to be identical to the distributions obtained from the plane and anti-plane elasticity solutions. Extensive results are given for axially and circumferentially cracked cylindrical shells, spherical shells, and toroidal shells under uniform in-plane shearing, out of plane shearing, and torsion. Taking advantage of the fact that the problem is formulated for “specially” orthotropic materials, the effect of orthotropy on the results is also studied in some detail.     

Studies on Elastic Coupling Parameters of Fracture Modes in Orthotropic Materials

The characterisation of cracks is usually done using the well known three basic fracture modes, namely opening, shearing and tearing modes. In isotropic materials these modes are uncoupled and provide a convenient way to define the fracture parameters. It is well known that these fracture modes are coupled in anisotropic materials. In the case of orthotropic materials also, coupling exists between the fracture modes, unless the crack plane coincides with one of the axes of orthotropy. The strength of coupling depends upon the orientation of the axes of orthotropy with respect to the crack plane and so the energy release rate components associated with each of the modes vary with crack orientation. The variation, of these energy release rate components with the crack orientation with respect to orthotropic axes, is analyzed in this paper. Results indicate that in addition to the orthotropic planes there exists other planes with reference to which fracture modes are uncoupled. This revised version was published online in July 2006 with corrections to the Cover Date.     

Near-tip dynamic fields for a crack advancing in a power-law elastic-plastic material: Modes I, II and III

Near-tip dynamic asymptotic stress fields of a crack advancing in an incompressible power-law elastic-plastic material are presented. It is shown that the stress- and strain-singularity are, respectively, of the order (In(R₀/r))^1/(n−1) and (In(R₀/r))^n/(n−1), where R₀ is a length parameter, r measures distance from the crack tip, and n is the power-law exponent. The angular variations of these fields are identical with those corresponding to dynamic crack growth in an elastic-perfectly-plastic material (Gao and Nemat-Nasser, 1983a,b).     

The crack problem in a specially orthotropic shell with double curvature

In this paper the crack problem of a shallow shell with two nonzero curvatures is considered. It is assumed that the crack lies in one of the principal planes of curvature and the shell is under Mode I loading condition. The material is assumed to be specially orthotropic. After giving the general formulation of the problem the asymptotic behavior of the stress state around the crack tip is examined. The analysis is based on Reissner's transverse shear theory. Thus, as in the bending of cracked plates, the asymptotic results are shown to be consistent with that obtained from the plane elasticity solution of crack problems. Rather extensive numerical results are obtained which show the effect of material orthotropy on the stress intensity factors in cylindrical and spherical shells and in shells with double curvature. Other results include the stress intensity factors in isotropic toroidal shells with positive or negative curvature ratio, the distribution of the membrane stress resultant outside the crack and the influence of the material orthotropy on the angular distribution of the stresses around the crack tip.     

Plane-strain mixed-mode near-tip fields in elastic perfectly plastic solids under small-scale yielding conditions

Within the context of the small-strain approach, plane-strain mixed-mode near-tip fields of a stationary crack in an elastic perfectly plastic Mises solid under small-scale yielding conditions are examined by finite element methods. Steady-state stress fields in the immediate vicinity of the crack tip are obtained as the remote loading of the elastic K-field increases. Asymptotic crack-tip solutions consisting of constant stress sectors, centered fan sectors, and an elastic sector are then constructed accordingly. The asymptotic crack-tip stress solutions agree well with the numerical results for a whole spectrum of mixed-mode loadings. Our mixed-mode near-tip solution with an elastic sector differs from that of Saka et al. by one (plastic) constant stress sector situated between the elastic sector and the neighbouring fan sector. The effect of the existence of the elastic sector on the near-tip fields is discussed in the light of the computational results. The plastic mixity factor of the near-tip field is given as a function of the elastic mixity factor of the prescribed K-field. This function is well bounded by that of the perfectly plastic limit of the corresponding solutions for power-law hardening materials given by Shih. Some issues pertaining to the numerical procedures such as the implementation of the small-scale yielding assumption are also addressed.     

On the inclined crack problem in an orthotropic medium under biaxial loading

The problem of an inclined crack in an orthotropic medium under biaxial loading is analyzed. A suitable coordinate transformation is introduced and two decoupled systems of the Cauchy–Riemann type are obtained in terms of complex potentials. The crack problem is solved by using the method of analytic continuation and closed form expressions of the near tip stress and displacement fields are derived. The influence of load biaxiality on the stress intensity factors, as well as on the local stress components is studied and graphically represented. Moreover, the action of material orthotropy on various quantities describing the crack characteristic is pointed out.     

Effects of non-singular stresses on crack-tip fields for pressure-sensitive materials,Part 1: Plane strain case

Mode I near-tip stress fields for elastic perfectly plastic pressure-sensitive materials under plane strain and small-scale yielding conditions are presented. A Coulomb-type yield criterion described by a linear combination of the effective stress and the hydrostatic stress is adopted in the analysis. The finite element computational results sampled at the distance of a few crack opening displacements from the tip show that, as the pressure sensitivity increases, the magnitudes of the normalized radial and hoop stress ahead of the tip decrease, the total angular span of the singular plastic sectors decreases, and the angular span of the elastic sectors bordering the crack surfaces increases. When non-singular T stresses are considered along the boundary layer of the small-scale yielding model, the near-tip stresses decrease as the T stress decreases. The plastic zone shifts toward the crack surfaces as the T stress increases. When the discontinuities of the radial stress and the out-of-plane normal stress along the border between the plastic sector and the elastic sector are allowed, the angular variations of the asymptotic crack-tip fields agree well with those of the finite element computations. Variation of the Q stresses for pressure-sensitive materials can be found from the asymptotic solutions when the plastic zone size ahead of the tip is relatively larger than the crack opening displacement. In addition the T stress is shown to have strong effects on the plastic zone sizes and shapes which could affect the toughening of pressure-sensitive materials.