Two Round Holes in Infinite Plate under General In-plane Load

Stresses around two interacting round holes of diameters D1 and D2, centers s apart on the x axis, in an infinite plate under any uniform in-plane far-field state given as stress resultants Nxx, Nyy, Nxy and the plate thickness t (far-field stresses N/t). The field is the superposition of the two single-hole Laurent expansions carried to convergence (the problem of Ling 1948 and Haddon 1967): the stress components, principal stresses and maximum shear at a user input point about hole 1, the hoop-stress extremes on each hole with Kt referred to the maximum far-field principal stress S1, the ligament stresses, and the stress field around the holes.

  • Two round holes of any diameters, centers s apart, in an infinite plate under Nxx, Nyy, Nxy (biaxial tension or compression and shear)

  • Stress components, principal stresses and principal angle at a user input point about hole 1

  • Hoop-stress extremes on each hole, Ktg = σedge,max / S1 for each hole and the ligament stresses between them

  • Fringe plots and line plots of the stress field around both holes

Loads and Hole Geometry


Unit System: /

Convention

Two round holes on the x axis, centers s apart, in an infinite plate under N_xx, N_yy, N_xy; point P at (r, theta) about hole 1, theta from the y axis

Far-field Loads and Thickness





Two Holes




Stress State Location



×

Two Round Holes in Infinite Plate under General In-plane Load

Two Round Holes in an Infinite Plate under General In-plane Load

Problem Statement

An infinite thin plate of thickness t contains two round holes of diameters D1 and D2 (radii a1, a2) whose centers lie on the x axis a distance s apart, hole 1 at the origin and hole 2 at x = s. Far from the holes the plate carries a uniform in-plane stress state given by the resultants Nxx, Nyy, Nxy, so that σxx∞ = Nxx/t acts along the line of centers, σyy∞ = Nyy/t across it and τxy∞ = Nxy/t in shear. The holes interact: under a load across the line of centers the ligament between them crowds the stress and the peaks on the x axis of each hole rise above Kirsch's 3σ (at the outer points for moderate spacing, where the other hole's σxx lowers the facing point, and at the facing points once the holes are close); under a load along the line each hole shields the other and the peaks fall below 3σ and move away from the ligament; a small hole beside a large one sees the large hole's own concentration as its far field and can reach twice Kirsch's factor. The questions are the stress at a point P(r, θ) about hole 1, θ measured from the y axis towards x, the peak hoop stress on each hole, the gross factors Ktg = σmax/S1 referred to the maximum far-field principal stress, and the stresses in the ligament.

Two round holes on the x axis under N_xx, N_yy, N_xy; P at (r, theta) about hole 1

Figure A: Convention. Hole 1 is at the origin and hole 2 at x = s; θ is measured from the y axis towards x about each hole's own center, the same convention as the other hole tools, so θ = 90° on hole 1 and θ = −90° on hole 2 are the two points facing each other.

Key Assumptions

Infinite plate, two round holes, linear elastic, plane stress, uniform far field. Both holes are traction-free and carry no resultant force; the solution does not depend on the elastic constants (it is the same in plane strain). The material stays elastic at the peak. The series is carried to convergence as long as the gap between the edges is at least a fifth of the larger diameter (equal holes: D/s ≤ 0.833). A finite plate, a third hole or a row of holes are outside the tool (the row has its own tool).

Solution by Superposed Single-hole Expansions

With Muskhelishvili's potentials Φ(z) and Ψ(z) the stresses are σxx + σyy = 4 Re Φ and σyy − σxx + 2iτxy = 2[z̅ Φ′ + Ψ]. Each traction-free hole without a resultant force contributes the Laurent terms αn (z − c)−n to Φ and βn (z − c)−n to Ψ, n ≥ 2 (Kirsch's solution is n = 2, 4), about its own center c. For two holes the two expansions are added to the far field,

Φ(z) = Γ + ∑j=1,2 ∑n αn(j) (z − cj)−n, z̅ Φ′ + Ψ = Γ′ + ∑j ∑n [ βn(j) − n αn(j) e−2iφj ] (z − cj)−n,

with φj the polar angle about hole j, Γ = (σxx∞ + σyy∞)/4 and Γ′ = (σyy∞ − σxx∞)/2 + iτxy∞. The geometry is symmetric about the x axis, so the normal loads give real coefficients and the shear load imaginary ones; the two are solved separately and superposed. Each hole's terms leave tractions on the other hole, so the coefficients of both holes are found together by imposing σrr = τrθ = 0 at Nj + 12 points of the upper half of each hole in a least-squares sense (Horii and Nemat-Nasser 1985), the order Nj following aj/(s − ak) (8 at D/s = 0.1, 19 at 0.5, 52 at 0.8 for equal holes; a large hole beside a small one needs the higher order on the large hole). Ling (1948) solved the equal-hole problem exactly in bipolar coordinates and Haddon (1967) the unequal-hole problem by complex potentials; the series reproduces both.

Stress Concentration Factors

The hoop stress on each edge is scanned for its extremes. The gross factors Ktg = σmax/S1 are referred to the maximum far-field principal stress, as in the other general-load tools (S2 when S1 is not tensile). In an infinite plate the ligament has no net-section stress of its own; Peterson's charts for the load across the line of centers use the convention σnom = σ √(1 − λ²)/(1 − λ), λ = D/s, which the tool prints for equal holes together with the fitted polynomials of Pilkey's tables:

tension along the line of centers (Chart 4.21): Ktg = 3.000 − 0.712λ + 0.271λ²,
tension across the line of centers (Chart 4.22): Ktn = 3.0000 − 3.0018λ + 1.0099λ²,
equal biaxial tension (Chart 4.23): Ktn = 2.000 − 2.119λ + 2.493λ² − 1.372λ³.

Table 1 compares the converged series with these fits and Table 2 with the values of Ling, Haddon and Chiang collected by Badiger and Ramakrishna (2024), for equal holes and for a hole beside a larger one. The series reproduces the exact values to every printed digit; the fits are within 1 % of it for λ ≤ 0.6 and the across-the-line fit falls below it by 2 % at λ = 0.7 and 8 % at 0.8. The tool's values are the converged series.

Table 1: Two equal holes, series against the fits (Ktn on σ√(1 − λ²)/(1 − λ), Ktg on the gross stress)
D/salong: Ktgfitacross: Ktgacross: Ktnfitbiaxial: Ktnfit
0.12.9802.9323.0012.7142.7101.8191.812
0.22.9272.8683.0042.4532.4401.6731.665
0.32.8532.8113.0152.2122.1901.5571.552
0.42.7772.7593.0351.9871.9611.4651.464
0.52.7152.7123.0661.7701.7521.3921.392
0.62.6702.6703.1161.5581.5631.3301.330
0.72.6422.6343.3811.4201.3941.2721.268
0.82.6232.6044.0341.3451.2451.2191.198
Table 2: Gross factor of the pair against the published values (gap between the edges g in units of the radius a of hole 1; hole 2 has the radius b)
b/ag/aloadseriesLing 1948Haddon 1967Chiangfinite elements
11along2.6502.6232.6502.6252.682
12along2.7152.7032.7152.7032.752
14along2.8272.8252.8272.8252.863
11across3.2643.2643.2643.2643.316
12across3.0663.0663.0663.0663.108
14across3.0203.0203.0203.0203.062
21along2.9252.9242.961
31along2.9822.9823.021
41along2.9962.9933.035
21across4.0514.0514.112
31across4.8544.8544.899
41across5.5455.5375.583

For the unequal pairs the peak of the load across the line sits on the small hole (4.051 at b/a = 2 against 3.263 on the large one), and the peak of the load along the line sits on the large hole, which is nearly undisturbed (2.925 to 2.996) while the small hole in its shadow falls to 2.0, 1.5 and 1.1.

Checks

A center distance of 2000 radii gives Kirsch's 3 and −1 on each hole; the field of two equal holes is symmetric about x = s/2 to 10−14; the root-mean-square traction left on the two edges is below 3 × 10−8 of the load (10−13 for close holes, where the order is higher); pure shear at s = 4a gives 4.55 against 4 for a single hole; the results do not change when the order is doubled.

Limitations

The tool stops when the gap between the edges is less than a fifth of the larger diameter (equal holes: D/s > 0.833), where the series needs more than 60 terms per hole and the page takes several seconds; Ling's and Haddon's solutions reach closer spacings. It does not cover a finite plate, a hole near a straight edge, three or more holes, or holes of different shape. The hoop stress in Figure 1 is referred to the nearer hole center, so the color field has a seam where the reference center changes.




References

Ling, C. B. (1948). "On the Stresses in a Plate Containing Two Circular Holes". Journal of Applied Physics, Vol. 19, pp. 77–82. Exact solution for two equal holes in bipolar coordinates under tension along and across the line of centers and under equal biaxial tension; the tabulated factors (2.623 / 2.703 / 2.825 along the line, 3.264 / 3.066 / 3.020 across it for gaps of 1, 2 and 4 radii) that the tool reproduces.

Haddon, R. A. W. (1967). "Stresses in an Infinite Plate with Two Unequal Circular Holes". Quarterly Journal of Mechanics and Applied Mathematics, Vol. 20, No. 3, pp. 277–291. Two unequal holes by complex potentials; the source of Peterson's Charts 4.24 and 4.25 and of the unequal-hole values (4.051, 4.854, 5.537 for a hole beside one of 2, 3 and 4 times its radius at a gap of one radius) the tool reproduces.

Davies, G. A. O., and Hoddinott, R. (1963). "Stresses in a Plate Pierced by Two Unequal Circular Holes". Journal of the Royal Aeronautical Society, Vol. 67. Extension of Ling's bipolar solution to unequal holes and a check of simple superposition for a small hole beside a large one.

Pilkey, W. D., and Pilkey, D. F. (2008). Peterson's Stress Concentration Factors, 3rd ed. John Wiley & Sons, Hoboken. Section 4.3.10 and Charts 4.21–4.25: two equal holes under tension parallel and perpendicular to the line of centers and under biaxial tension, and two unequal holes; the nominal stress σnom = σ √(1 − (d/l)²)/(1 − d/l) used for the net factors.

Pilkey, W. D. (2005). Formulas for Stress, Strain, and Structural Matrices, 2nd ed. John Wiley & Sons, Hoboken. Chapter 6, Table 6-1, case 4 "Two equal circular holes in infinite plate": the fitted polynomials Kt = 3.000 − 0.712(d/L) + 0.271(d/L)² (tension along the line, gross stress), Kt = 3.0000 − 3.0018(d/L) + 1.0099(d/L)² (tension across the line, nominal stress above) and Kt = 2.000 − 2.119(d/L) + 2.493(d/L)² − 1.372(d/L)³ (equal biaxial tension), printed by the tool for comparison. Copy of the chapter filed in the site's References folder.

Badiger, S. M., and Ramakrishna, D. S. (2024). "Interaction of Two Circular Holes in an Infinite Plate by Body Force Method". Materials Physics and Mechanics, Vol. 52, No. 3, pp. 96–107. Tables 1 and 2 collect Ling's, Haddon's and Chiang's factors for equal and unequal holes with body-force and finite element values; used to check the tool (copy filed in the site's References folder).

Horii, H., and Nemat-Nasser, S. (1985). "Elastic Fields of Interacting Inhomogeneities". International Journal of Solids and Structures, Vol. 21, No. 7, pp. 731–745. Interacting holes by superposed single-hole Laurent expansions with the coefficients fixed on each boundary, the construction the tool uses.

Muskhelishvili, N. I. (1953). Some Basic Problems of the Mathematical Theory of Elasticity. Noordhoff, Groningen. The complex potentials Φ(z), Ψ(z) and the stress combinations σxx + σyy = 4 Re Φ, σyy − σxx + 2iτxy = 2[z̅Φ′ + Ψ] on which the solution is built.

Timoshenko, S. P., and Goodier, J. N. (1970). Theory of Elasticity, 3rd ed. McGraw-Hill, New York. Chapter 4, "Effect of circular holes on stress distributions in plates": the Kirsch limit Kt = 3 for the single hole.


DISCLAIMER: The methods and calculation on aNOWysis.com are provided as means to educate the reader. There is no warranty, expessed or implied, as to the relevance or accuracy of the methods or calculations. All calculations should be checked by a qualified analyst for appropriateness and accuracy. We will not be held liable for any damages, injuries or deaths caused by using methods or calculations from this site.