Research Papers

Analytical Solutions and Stress Concentration Factors for Annuli With Inhomogeneous Boundary Conditions

[+] Author and Article Information
S. Shahzad

Department of Civil Engineering,
Aalto University,
Rakentajanaukio 4,
Espoo 02150, Finland

J. Niiranen

Department of Civil Engineering,
Aalto University,
Rakentajanaukio 4,
Espoo 02150, Finland

Contributed by the Applied Mechanics Division of ASME for publication in the JOURNAL OF APPLIED MECHANICS. Manuscript received December 20, 2017; final manuscript received April 18, 2018; published online May 10, 2018. Assoc. Editor: Thomas Siegmund.

J. Appl. Mech 85(7), 071008 (May 10, 2018) (14 pages) Paper No: JAM-17-1686; doi: 10.1115/1.4040079 History: Received December 20, 2017; Revised April 18, 2018

Analytical displacement and stress fields with stress concentration factors (SCFs) are derived for linearly elastic annular regions subject to inhomogeneous boundary conditions: an infinite class of the mth order polynomial antiplane tractions or displacements. The solution of the Laplace equation governing the out-of-plane problem covers both rigid and void circular inclusions forming the core of the annulus. The results show first that the SCF and the loading order are inversely proportional. In particular, the SCF approaches value 2 when either the outer boundary of the annulus tends to infinity or the order of the polynomial loading increases. Second, the number of peculiar points on the inner contour having null stress increases with the increasing loading order. The analytical solution is confirmed and extended to noncircular enclosures via finite element analysis by exploiting the heat-stress analogy. The results show that the closed-form solution for a circular annulus can be used as an accurate approximation for noncircular enclosures. Altogether, the results shown can be exploited for analyzing complex loading conditions and/or multiple rigid or void inclusions for enhancing the design of hollow and reinforced composites materials.

Copyright © 2018 by ASME
Your Session has timed out. Please sign back in to continue.


Zou, W. N. , Lee, Y. G. , and He, Q. C. , 2017, “Inclusions Inside a Bounded Elastic Body Undergoing Anti-Plane Shear,” Math. Mech. Solids, 23(4), pp. 588–605. [CrossRef]
Zappalorto, M. , Lazzarin, P. , and Yates, J. R. , 2008, “Elastic Stress Distributions for Hyperbolic and Parabolic Notches in round Shafts Under Torsion and Uniform Antiplane Shear Loadings,” Int. J. Solids Struct., 45(18–19), pp. 4879–4901. [CrossRef]
Shahzad, S. , and Dal Corso, F. , 2018, “Torsion Elastic Solids with Sparsely Distributed Collinear Voids,” (in preparation).
Radaj, D. , 1990, Design and Analysis of Fatigue Resistant Welded Structures, Elsevier, Cambridge, UK. [CrossRef]
Savin, G. N. , 1961, Stress Concentration around Holes, Pergamon Press, London.
Muskhelishvili, N. I. , 1963, Some Basic Problems of the Mathematical Theory of Elasticity, P. Noordhoff Ltd., Groningen, The Netherlands.
Slaughter, W. S. , 2002, The Linearized Theory of Elasticity, Springer Science & Business Media, New York. [CrossRef]
Toubal, L. , Karama, M. , and Lorrain, B. , 2005, “Stress Concentration in a Circular Hole in Composite Plate,” Compos. Struct., 68(1), pp. 31–36. [CrossRef]
Givoli, D. , and Elishakoff, I. , 1992, “Stress Concentration at a Nearly Circular Hole With Uncertain Irregularities,” ASME. J. Appl. Mech., 59(2S), pp. 65–71. [CrossRef]
Hoff, N. J. , 1981, “Stress Concentrations in Cylindrically Orthotropic Composite Plates With a Circular Hole,” ASME. J. Appl. Mech., 48(3), pp. 563–569. [CrossRef]
Lubarda, V. A. , 2003, “Circular Inclusions in Anti-Plane Strain Couple Stress Elasticity,” Int. J. Solids Struct., 40(15), pp. 3827–3851. [CrossRef]
Meguid, S. A. , and GongStress, S. X. , 1993, “Stress Concentration Around Interacting Circular Holes: A Comparison Between Theory and Experiments,” Eng. Fract. Mech., 44(2), pp. 247–256. [CrossRef]
Misseroni, D. , Dal Corso, F. , Shahzad, S. , and Bigoni, D. , 2014, “Stress Concentration Near Stiff Inclusions: Validation of Rigid Inclusion Model and Boundary Layers by Means of Photoelasticity,” Eng. Fract. Mech., 121–122, pp. 87–97. [CrossRef]
Dal Corso, F. , Bigoni, D. , Noselli, G. , Misseroni, D. , and Shahzad, S. , 2014, “Rigid Inclusions: Stress Singularity, Inclusion Neutrality and Shear Bands,” Third International Conference on Fracture, Fatigue and Wear, Kitakyushu, Japan, Sept. 1–3, pp. 40–42.
Setiawan, Y. , Gan, B. S. , and Han, A. L. , 2017, “Modeling of the ITZ Zone in Concrete: Experiment and Numerical Simulation,” Comput. Concr., 19(6), pp. 647–655.
Ru, C. Q. , and Schiavone, P. , 1997, “A Circular Inhomogeneity With Circumferentially Inhomogeneous Interface in Antiplane Shear,” Proc. R. Soc. A, 453(1967), pp. 2551–2572. [CrossRef]
Zappalorto, M. , Lazzarin, P. , and Berto, F. , 2009, “Stress Concentration Near Holes in the Elastic Plane Subjected to Antiplane Deformation,” Mater. Sci., 48(4), pp. 415–426.
Sih, G. C. , 1965, “Stress Distribution Near Internal Crack Tips for Longitudinal Shear Problems,” ASME J. Appl. Mech., 32(1), pp. 51–58. [CrossRef]
Kohno, Y. , and Ishikawa, H. , 1995, “Singularities and Stress Intensities at the Corner Point of a Polygonal Hole and Rigid Polygonal Inclusion Under Antiplane Shear,” Int. J. Eng. Sci., 33(11), pp. 1547–1560. [CrossRef]
Seweryn, A. , and Molski, K. , 1996, “Elastic Stress Singularities and Corresponding Generalized Stress Intensity Factors for Angular Corners Under Various Boundary Conditions,” Eng. Fract. Mech., 55(4), pp. 529–556. [CrossRef]
Lubarda, V. A. , 2015, “On the Circumferential Shear Stress Around Circular and Elliptical Holes,” Arch. Appl. Mech., 85(2), pp. 223–235. [CrossRef]
Movchan, A. B. , Movchan, N. V. , and Poulton, C. G. , 2002, Asymptotic Models of Fields in Dilute and Densely Packed Composites, Imperial College Press, London. [CrossRef]
Chen, J. , and Wu, A. , 2006, “Null-Field Approach for the Multi-Inclusion Problem Under Antiplane Shears,” ASME. J. Appl. Mech., 74(3), pp. 469–487. [CrossRef]
Sendeckyj, G. P. , 1971, “Multiple Circular Inclusion Problems in Longitudinal Shear Deformation,” J. Elasticity, 1(1), pp. 83–86. [CrossRef]
Gong, S. X. , 1995, “Antiplane Interaction Among Multiple Circular Inclusions,” Mech. Res. Commun., 22(3), pp. 257–262. [CrossRef]
Bacca, M. , Dal Corso, F. , Veber, D. , and Bigoni, D. , 2013, “Anisotropic Effective Higher-Order Response of Heterogeneous Cauchy Elastic Materials,” Mech. Res. Commun., 54, pp. 63–71. [CrossRef]
Bigoni, D. , and Drugan, W. J. , 2007, “Analytical Derivation of Cosserat Moduli Via Homogenization of Heterogeneous Elastic Materials,” ASME J. Appl. Mech., 74(4), pp. 741–753. [CrossRef]
Schiavone, P. , 2003, “Neutrality of the Elliptic Inhomogeneity in the Case of Non-Uniform Loading,” Int. J. Eng. Sci., 41(18), pp. 2081–2090. [CrossRef]
Van Vliet, D. , Schiavone, P. , and Mioduchowski, A. , 2003, “On the Design of Neutral Elastic Inhomogeneities in the Case of Non-Uniform Loading,” Math. Mech. Solids, 41(2), pp. 2081–2090.
Vasudevan, M. , and Schiavone, P. , 2006, “New Results Concerning the Identification of Neutral Inhomogeneities in Plane Elasticity,” Arch. Mech., 58(1), pp. 45–58. http://am.ippt.pan.pl/am/article/view/v58p45
Dal Corso, F. , Shahzad, S. , and Bigoni, D. , 2016, “Isotoxal Star-Shaped Polygonal Voids and Rigid Inclusions in Nonuniform Antiplane Shear Fields—Part I: Formulation and Full-Field Solution,” Int. J. Solids Struct., 85–86, pp. 67–75. [CrossRef]
Dal Corso, F. , Shahzad, S. , and Bigoni, D. , 2016, “Isotoxal Star-Shaped Polygonal Voids and Rigid Inclusions in Nonuniform Antiplane Shear Fields—Part II: Singularities, Annihilation and Invisibility,” Int. J. Solids Struct., 85–86, pp. 76–88. [CrossRef]
Shahzad, S. , 2016, “Stress Singularities, Annihilations and Invisibilities Induced by Polygonal Inclusions in Linear Elasticity,” Ph.D. thesis, University of Trento, Trento, Italy. http://eprints-phd.biblio.unitn.it/1769/
Shahzad, S. , Dal Corso, F. , and Bigoni, D. , 2016, “Hypocycloidal Inclusions in Nonuniform Out-of-Plane Elasticity: Stress Singularity vs Stress Reduction,” J. Elasticity, 126(2), pp. 215–229. [CrossRef]
Shahzad, S. , and Niiranen, J. , 2018, “Analytical Solution With Validity Analysis for an Elliptical Void and Rigid Inclusion Under Uniform and Nonuniform Antiplane Loading,” (in preparation).
Polyanin, A. D. , 2002, Handbook of Linear Partial Differential Equations for Engineers and Scientists, Chapman and Hall/CRC, Boca Raton, FL.
Morrison, N. , 1994, Introduction to Fourier Analysis, Wiley, New York.
Carslaw, H. S. , and Jaeger, J. C. , 1959, Conduction of Heat in Solids, Oxford University Press, London.
Chen, Z. , Zhu, W. , Di, Q. , and Wang, W. , 2015, “Burst Pressure Analysis of Pipes With Geometric Eccentricity and Small Thickness-to-Diameter Ratio,” J. Pet. Sci. Eng., 127, pp. 452–458. [CrossRef]


Grahic Jump Location
Fig. 1

Annular cross sections composed of an inclusion having a circular contour of radius a embedded in a matrix with elliptical contour (defined by the semimajor axis b and the semiminor axis R̂) and of regular polygons circumscribed by circle of radius R̂. In particular, the extreme cases of circular inclusions defined by the Greek letter χ are taken into consideration, see Eq. (5).

Grahic Jump Location
Fig. 2

Contour plots under 3D shaded surface plots for a circular annulus undergoing nonuniform polynomial antiplane displacement fields w(r, θ) corresponding to uniform, linear, and quadratic (m = 0, 1, 2) antiplane shearing loads (b0(m)=0)

Grahic Jump Location
Fig. 3

Cases of dimensionless radial shear stress τrz(m)(r,0)/(amF(m)(θ)) (flat curves between the stress values of 2 and 1) and circumferential shear stress component τθz(m)(r,0)/(am∂F(m)(θ)/∂θ) (exploding curves between the stress values of 2 and ) as a function of the volume fractions f undergoing different loading orders (m = 0, 1, 2, 5, 10)

Grahic Jump Location
Fig. 4

Dimensionless shear stress τ(m)(r,0)/τ∞(m)(a,0) along the horizontal line bisecting the circular annulus with different values of the volume fractions f as a function of dimensionless radial distance r/a is reported. In particular, cases of uniform (m = 0) and linear (m = 1) antiplane loading conditions activating stress concentrations for void (b0(m)=0) and for rigid inclusion (c0(m)=0) are considered. The results show that for thin annulus, the void becomes very dangerous from the failure point of view due to the high stress concentration with respect to the rigid inclusion. If inclusion becomes very small, then the SCF value 2 is reached for both the boundary value problems but the void approaches the value from above and the rigid inclusion from below.

Grahic Jump Location
Fig. 5

Loading order m = 0 remains the same, while aspect ratio f varies. Contour plots of the dimensionless shear stress modulus τ(0)(r,θ)/τ∞(0)(r,θ) (top row) and the values of SCF τ(0)(a,θ)/τ∞(0)(a,θ) (bottom row) along the inner contour subject to uniform antiplane shear (b0(0)=0) are reported. Results show that, when f → 1, then SCF tends to infinity, while for f → 0 the SCF value tends to 2 (as shown in Refs. [19], [21], and [32]). Notice that the analytical solution is confirmed through the comparison of SCF analytical curves (solid line) and numerical curves (dashed line) obtained via finite element analysis discussed comprehensively in Sec. 4.

Grahic Jump Location
Fig. 6

The aspect ratio of the annulus remains constant, while the loading order m varies. Contour plots of the dimensionless shear stress modulus τ(m)(r,θ)/τ∞(m)(a,θ) (top row) and the values of SCFs τ(m)(a,θ)/τ∞(m)(a,θ) (bottom row) along the inner contour of the circular hole subject to uniform (m = 0), linear (m = 1) and quadratic (m = 2) antiplane loading (b0(m)=0) are reported. Results show that the value of SCF along the boundary of circular hole is inversely proportional to the order of the polynomial loading. Whereas the number of annihilation points (with null stress) occurring at the inner contour is directly proportional to the loading order. The number of annihilation points is given by 2(m + 1).

Grahic Jump Location
Fig. 7

A rectangular (I) and an elliptical enclosure (II) modeling hollow core slab and symmetrical hollow airfoil, respectively, are meshed at two levels with free triangular elements

Grahic Jump Location
Fig. 8

Contour lines of the dimensionless shear stress modulus τ(m)(r,θ)/τ∞(m)(r,θ) for the value of volume fraction f = 0.04 (i.e., a/R = 1/5, top row) and f = 0.25 (i.e., a/R = 1/2, bottom row) and the values of SCF (middle column) along the inner contour of a circular annulus subject to uniform (b0(0)=0) show that, even a small amount of relative eccentricity d¯=d/a leads to the higher SCF

Grahic Jump Location
Fig. 9

Values of the SCF are shown for noncircular annulus (circular void and square enclosure) subject to uniform (top row) and linear (bottom row) antiplane loading condition (b0(m)=0). Different aspect ratios are considered for the noncircular annulus to evaluate the influence on the SCF value. Results show that, for uniform loading, the error is greater than 1%, while already for the linear loading, the analytical solution is very close to the numerical one.

Grahic Jump Location
Fig. 10

Cracking of mortar specimens containing cylindrical inclusion. Source: Courtesy by Han and co-workers [15].

Grahic Jump Location
Fig. 11

Simulation of a hollow core slab subject to uniform (m = 0) antiplane shear stress (b0(0)=0) is addressed showing that the interaction between the inclusions does not have significant influence on the stress concentration values. While the horizontal vicinity of the external boundary increases significantly the SCF value.

Grahic Jump Location
Fig. 12

An elliptical enclosure approximating a symmetrical airfoil subject to uniform (m = 0) antiplane shear stress (b0(0)=0) is reported and the interaction between the inclusions does not have significant influence on the stress concentration values

Grahic Jump Location
Fig. 13

A quarter of a square prism with collinear circular void subject to uniform antiplane shear (m = 0): dimensionless shear stress field (top left), cut plane (bottom left), and dimensionless stress curves on the perimeter p of the circular arc of the cut plane




Some tools below are only available to our subscribers or users with an online account.

Related Content

Customize your page view by dragging and repositioning the boxes below.

Related Journal Articles

Sorry! You do not have access to this content. For assistance or to subscribe, please contact us:

  • TELEPHONE: 1-800-843-2763 (Toll-free in the USA)
  • EMAIL: asmedigitalcollection@asme.org
Sign In