An exact linear elastic solution is derived for a pair of opposing point loads, which act in the middle of a straight slot having rounded ends. This problem is analogous to the mode III crack problem where two concentrated forces act in opposite directions to open the surfaces of a finite-length crack. The corresponding path independent J integral for this slot problem is also determined.
Prior to this analysis, no exact solution has been found for a linear slot with rounded ends subject to two opposing antiplane point loads. Perhaps the closest analogy that one may find in the literature is for two opposing point loads acting on the surface of an elliptical hole . However, experience has shown that for remotely applied shear tractions, large differences exist in stress concentration factors for solutions of a slot having this particular geometry and those of an ellipse with identical aspect ratio of slot length to slot width [2, 3]. This observation suggests that it is questionable in general to use the solution of an elliptical hole problem as a close approximation of a solution for a slot.
The shape of the slot analyzed here comes from the free-streamline problem of an ideal fluid flowing past two side-by-side flat plates  where the flow at infinity is perpendicular to the plates. The analogy between this solution and that of an antiplane slot problem with a uniform shear loading at infinity was explored in detail in [5, 6]. Further, linear elastic solutions in the plane for remotely applied tensile loadings of a slot of a similar shape was investigated in [7, 8] using numerical methods of analysis.
where a is a scaling factor with units of length and z is the complex variable
On the unit circle γ, the parameters of (4) are related to the angle θ by the following relationships
where ±θ0 and ±π ∓ θ0 are the four locations on the slot boundary (ρ = 1)where the flat surfaces meet the rounded surfaces. In Figure 1, point A marks one these locations in the first quadrant of the xy plane.
2 Stress analysis in cartesian coordinates
In the complex plane, a Westergaard stress function  is often used to solve mode III crack problems, which are a class of antiplane linear elastic problems for isotropic materials. Accordingly, the Westergaard stress function and its related functions have the following relationships
where τx (x, y) and τy (x, y) are antiplane shear stresses, φ (x, y) is a real-valued stress function, and δ (x, y) is a real-valued out-of-plane displacement. In (6), the parameter G is the shear modulus, which is assumed constant. The stresses (6) are the only stresses present in this class of problems. Similarly, δ (x, y) is the only component of displacement in an antiplane problem, which acts in a direction perpendicular to the xy plane.
Now equilibrium requires 
which is automatically satisfied by the stress function φ (x, y). The compatibility equation for strain further requires 
where γx and γy are the engineering shear strains in the x and y directions respectively. As the shear strains are related to the displacement as follows
the compatibility equations (8) are automatically satisfied by the displacement δ (x, y).
However, Hooke’s law imposes the following additional relationships between stress and strain
for equilibrium to be satisfied.
Similarly, by substituting the relationships for the shear stresses involving partial derivatives of φ from (6) into (10) to obtain their relationships with strain, one concludes by substituting these strains into (8) that φ must also satisfy Laplace’s equation in order that the strain compatibility equation be satisfied
Therefore, as relationships (6) indicate, φ and Gδ are conjugate harmonic functions. Consequently, their relationships are consistent with the general properties of the real and imaginary parts of a complex function of z , which in this case is the Westergaard stress function.
3 Stress analysis in curvilinear coordinates
Often, the solution of a problem with complicated boundaries can be simplified by the adoption of an appropriate orthogonal curvilinear coordinate system that accommodates the boundaries naturally. Let us suppose that it is useful to define such a system as follows
where f is an arbitrary function and u and v are real functions that constitute an orthogonal curvilinear coordinate system by virtue of the properties of a function of a complex variable .
It is the goal here to develop a connection between the curvilinear coordinate system defined by (13) and the previously defined functions of (6). Using the chain rule of partial differentiation, one finds that for the stress function φ
Certain geometric relationships exist between the antiplane shear stresses in the (u, v) directions respectively τu and τv and their Cartesian counterparts. By resolution of the in-plane stress vectors, one finds
where α is the angle between the x axis and the direction of a unit vector in the u direction.
At this point it is helpful to introduce the metric coefficient
Using relationships similar to those derived in 
Because of (12), the function φ must also satisfy Laplace’s equation, which in the (u, v) system assumes the form
4 Elliptical hole analysis
In order to gain insight into the slot solution, let us first explore the solution of the analogous elliptical hole point load problem, which has a known solution . The mapping function of the boundary of an elliptical hole in an infinite plate onto the surface and interior of a unit circle has the form 
In (25), a is the semimajor axis and b is the semiminor axis of the ellipse. The mapping function of the elliptical hole on the unit circle will have the same form as (3) for the associated polar coordinates.
Let us now define the following function, which has logarithmic singularities at ζ = ±i,
where A is a parameter to be determine from equilibrium in terms of the point loads P. By differentiation of (26) with respect to ζ , one finds that the function has 1/r singularities at ζ = ±i. These properties characterize stresses that are typically associated with point loads in linear elasticity theory. The derivative of (26) with respect to ζ is explicitly
Now assume a Westergaard stress function of the form
which has the following alternative representation in terms of ζ using the chain rule of differential calculus
where ω' (ζ) is the first derivative of (24) with respect to ζ,
Now to determine the constant A of (31) one needs to evaluate the integral appearing in the following expression
This relationship reflects a symmetry requirement that one half of the magnitude of the load P must be carried along the positive x axis from the tip of the semimajor axis to infinity. In terms of the coordinates defined by (3), the relationships to the right of the integral sign of (32) are
where the sign in (34) is determined by the directions of the loads.
Let us now introduce cylindrical elliptical coordinates (u, v) in terms of complex variable notation as follows
The equation of equilibrium in cylindrical elliptical co-ordinates for the antiplane problem is derived from a more general relationship given in  as
where u0 is the value of the coordinate u of the elliptical hole.
A plot of the stress function φ is shown in Figure 2 for c = 1, u0 = 0.2, and P = 1. Note that the stress function is constant along the elliptical hole u0 except for a jump that occurs at v = ±π/2, where the point loads are applied.
On the elliptical hole u = u0, one determines from (41) that the boundary condition on traction is satisfied as τu is everywhere zero except at the locations of the concentrated loads, which are singularities. This behavior is depicted in Figure 3.
A local maximum value of the shear stress, away from the point loads, is consistent with the value determined previously in 
5 Slot analysis
Now that experience has been gained by working through the simpler elliptical hole problem, the analogous slot solution of Figure 1 will be derived.
By taking the first derivative of (1) with respect to ζ , one finds
The functions analogous to those of (26)–(29) assume similar forms to those given previously for the elliptical hole. They provide upon substitution of (43) for ω' (ζ) in (29) the following form of the Westergaard function for the slot
The relationship between A and P is determined for the slot problem from equilibrium considerations in much the same way as for the elliptical hole. However, direct integration of the expressions involved are much more complicated than they were for the elliptical hole. See Appendix A for details of the evaluation of the integral indirectly. The result of this integration is identical to the value of A found previously for the elliptical hole in terms of the concentrated load P for the directions shown in Figure 1,
Along the slot surface u = 0 in this coordinate system. Now the real-valued stress function assumes the following form for the slot problem in these coordinates
The shear stresses follow from (22) as
The value of a local maximum of the shear stress, away from the point loads, is found at the four locations where the flat surfaces of the slot meet the rounded ends of the slot along its boundary. This value at point A of Figure 1 is
The location where the maximum stress occurs in a slot of this geometry can differ in how the load is applied. For example, in the case of uniform tensile loadings at infinity, the largest principal stress occurs at a point analogous to point M of Figure 1, provided the direction of the applied load is perpendicular to the long axis of the slot . Similarly, in the case of uniform biaxial tensile loadings of the slot at infinity, the maximum value of the principal stress falls somewhere between points M and A of Figure 1 along the curved portion of the slot depending on the value of m1.
6 Path independent integral for slot problem
The path-independent J integral has shown to be a useful tool in the analysis of crack growth problems. The geometry of the slot of Figure 1 can serve as a crack having a finite width with rounded tips. The J integral as defined in  reduces to the following form for this problem
where W is the strain energy density and Γ is the path of integration. The J integral for this particular slot problem can be reduced to the evaluation of this integral along a path limited to the rounded portion of the left hand side of the slot. This evaluation will be performed using the coordinate θ of (3). In , the angle θ was related to a conformal mapping of the region exterior to the slot onto the surface and interior of a unit circle. The positive direction for the angle θ on the unit circle is in the counterclockwise direction. In the physical plane, the angle θ = 0 defines a ray that begins at the origin of the xy coordinate system and continues along the negative x axis. A point on the surface of the slot traverses a path in the clockwise direction for a change of the angle θ in the counterclockwise direction.
The form of W for the antiplane problem addressed here is
Along the rounded portion of the of the left hand slot tip
Taking the derivative of y of (60) with respect to θ gives
Now, using properties of trigonometric functions, one can show
For the purposes of evaluation, this integral is further simplified by making the following substitutions into it
where the factor of two before the integral results from the symmetry of the integrand and an adjustment of the limits of integration. The integral in (65) is readily evaluated using the symbolic computer program Mathematica® as
The limit of J as θ0→ 0 corresponds to the value of a mode III crack subject to a pair of antiplane splitting forces applied at the middle of the crack
This is consistent with the value determined by substituting the stress intensity factor KIII for a pair of opposed point loads P applied at the middle of a line crack of length 2a , i.e.,
into a general relationship for the J integral  for mode III crack problems
The range of J/J0 varies monotonically from a value of one for θ0 = 0 to infinity for θ0 = π/2.
A prototype problem for the loading of a slot by forces applied near the center of its internal boundary was proposed and solved in this article The singularities that appear at the points of application of the concentrated loads were anticipated and consequently they are not particularly noteworthy. However, of considerable interest is how this symmetrical loading near the center of a slot generates states of stress at other locations along the slot boundary. For example, it was not known beforehand if the applied loads would generate larger states of stress at points A or M of Figure 1.
One notes that the linear elastic point load problem presented here is only the second analytical solution known for a slot with rounded ends. The first [2, 3, 5, 6] was found by analogy with the free streamline problem of Riabouchinsky . That particular slot solution assumes a load of uniform antiplane shear tractions at infinity rather than point loads applied along its interior boundary. In both cases, however, one finds large stresses that are proportional to 1/m11/4 at the slot tips. Consequently, as the geometry of the slot approaches that of a line crack (m1→ 0), the stresses at the tips become singular in a similar mathematical fashion with m1 for the two different types of loading.
Further, the J integral for the slot was found here in a straightforward manner that is impossible to duplicate with the analogous loading of an elliptical hole. If one attempts to calculate the J integral along the periphery of an elliptical hole, one encounters a divergent integral because of the presence of the singularities generated by the point loads. This issue is not encountered in the slot problem because the path of integration is limited to the rounded portions of the slot. Along the flat surfaces of the slot, where the point loads are applied, dy is zero, which negates the effect of the singularities on the convergence of the integral. In a similar light, the calculation of the J integral for the line crack cannot be determined by following a path along the surfaces of the crack as dy equals zero along them and there is no smooth curve connecting the surfaces. Nevertheless, one may determine the value of the J integral for the line crack (67) subject to opposing point loads from the solution for the slot (66). One simply allows the parameter m1 to approach zero, which causes the flat surfaces of the slot to converge to a line crack of length 2a with zero width.
 Neuber H. Kerbspannungslehre. 4th ed. Berlin: Springer; 2001:36. Search in Google Scholar
 Unger DJ. Linear elastic solutions for slotted plates. J Elast. 2012;108(1):67-82. Search in Google Scholar
 Unger DJ. Erratum to: linear elastic solutions for slotted plates. J Elast. 2012;108(1):83. Search in Google Scholar
 Riabouchinsky D., On steady fluid motions with free surfaces, Proc. Lond. Math. Soc. (2) 1921;19(1):206-215 Search in Google Scholar
 Unger DJ. Free streamline hydrodynamic analogy for a linear elastic antiplane slot problem with perfectly plastic ligaments at its ends. J Elast. 2018;132(2):261-70. Search in Google Scholar
 Unger DJ. Visualizing the crack driving force through fluid analogy. J Mech Behav Mater. 2019;28(1):89-94. Search in Google Scholar
 Unger DJ. Linear elastic solutions for slotted plates revisited, Proceedings 24th International Congress of Theoretical and Applied Mechanics, ed. JM Floryan, ICTAM, Montreal; 2017:2042-2043. Search in Google Scholar
 Unger DJ. Comparison of series and finite difference solutions to remote tensile loadings of a plate having a linear slot with rounded ends. J Mech Mater Struct. 2020;15(3):361-78. Search in Google Scholar
 Abramowitz M, Stegun IA. Handbook of Mathematical Functions with Formulas, Graphs and Mathematical Tables. National Bureau of Standards, Series 55, U.S. Printing Office, Washington, DC, 1964; 81:589-628. Search in Google Scholar
 Unger DJ. Analytical Fracture Mechanics. Dover: Mineola; 2011:29-34. Search in Google Scholar
 Churchill RV. Complex Variables and Applications. New York: McGraw-Hill. 1960. Search in Google Scholar
 Kachanov LM. Fundamentals of the Theory of Plasticity. Dover: Mineola; 2004:469-70. Search in Google Scholar
 Timoshenko SP, Goodier JN. Theory of Elasticity. 3rd ed. New York: McGraw-Hill; 1934:184. Search in Google Scholar
 Spiegel MR, Liu J. Mathematical Handbook of Formulas and Tables, 2nd ed., Shaum’s Outline, McGraw-Hill, New York; 1999:127. Search in Google Scholar
 Sokolnikoff IS. Mathematical Theory of Elasticity, 2nd ed., McGraw-Hill, New York; 1956;181:292. Search in Google Scholar
 Rice JR. A path independent integral and approximate analysis of strain concentrations by notches and cracks. J Appl Mech. 1968;35(2):379-86. Search in Google Scholar
 Tada H, Paris PC, Irwin GR. The Stress Analysis of Cracks Handbook. 3rd ed. New York: ASME; 2000, 138. Search in Google Scholar
 Broberg KB. Cracks and Fracture. San Diego: Academic Press; 1999, 77 Search in Google Scholar
Evaluation of parameter A
In this appendix, the parameter A of (27) will be evaluated from equilibrium considerations for the slot problem.
Upon integrating both sides of (A1), one finds between limits of integration
By substituting the following expressions into (A2) and expanding
Now restricting attention to the positive x axis, one may set both dy = 0 and ImZIII equal to zero in (A4). The latter is valid because it is equal to τx, which produces no traction along the x axis in the y direction. What remains of (A4) reflects the load that must be carried along the right hand side of the plane. For the direction of loading shown in Figure 1, one has
The complex variables w1 and w2 found in (A7) are defined in terms of the real-valued coordinates (u, v) by
© 2020 D. J. Unger, published by De Gruyter
This work is licensed under the Creative Commons Attribution 4.0 International License.