Next Article in Journal
Surface Roughening Behavior of the 6063-T4 Aluminum Alloy during Quasi-in Situ Uniaxial Stretching
Previous Article in Journal
Electrical and Optoelectronic Properties Enhancement of n-ZnO/p-GaAs Heterojunction Solar Cells via an Optimized Design for Higher Efficiency
Previous Article in Special Issue
Tensile-Tearing Fracture Analysis of U-Notched Spruce Samples
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Application of the Generalized Method of Moving Coordinates to Calculating Stress Fields near an Elliptical Hole

by
Sergei Alexandrov
1,2,
Marina Rynkovskaya
2 and
Shang-Nan Tsai
3,*
1
Ishlinsky Institute for Problems in Mechanics RAS, 101-1 Prospect Vernadskogo, 119526 Moscow, Russia
2
Department of Civil Engineering, Peoples’ Friendship University of Russia (RUDN University), 6 Miklukho-Maklaya St, 117198 Moscow, Russia
3
Department of Mechanical and Electro-Mechanical Engineering, National Sun Yat-sen University, No. 70 Lien-Hai Rd., Kaohsiung 804, Taiwan
*
Author to whom correspondence should be addressed.
Materials 2022, 15(18), 6266; https://doi.org/10.3390/ma15186266
Submission received: 30 July 2022 / Revised: 2 September 2022 / Accepted: 6 September 2022 / Published: 9 September 2022
(This article belongs to the Special Issue Fatigue and Fracture of Materials)

Abstract

:
The distribution of stresses near holes is of great importance in fracture mechanics and material modeling. The present paper provides a general stress solution near a traction-free surface for an arbitrary piecewise linear yield criterion, assuming plane-strain conditions. The generalized method of moving coordinates is proven efficient in this case. In particular, the solution reduces to evaluating one ordinary integral. The boundary value problem solved is a Cauchy problem for a hyperbolic system of equations. Therefore, the stress solution in the plastic region is independent of other boundary conditions, though the occurrence of plastic yielding at a specific point is path-dependent. The general solution applies to calculating the stress field near an elliptic hole. It is shown that the parameter that controls the pressure-dependency of the yield criterion affects the stress field significantly. The aspect ratio is less significant as compared to that parameter. However, for a given material, the aspect ratio should also be considered to predict the stress field accurately, especially in the near vicinity of the hole. The solution reduces to an available solution for the pressure-independent yield criterion, which is a particular yield criterion of the considered class of yield criteria.

1. Introduction

The distribution of stresses near holes and cracks is of great importance in fracture mechanics and material modeling. A vast amount of literature is devoted to developing methods for finding such distributions using the theories of elasticity and plasticity. The linearly elastic stress field at the base of a sharp stationary crack has been found in [1]. This analysis has been extended to the plastic range in [2,3,4].
The importance of stress solutions near holes is twofold. Firstly, such solutions are required to understand the fracture process. Paper [5] has presented a stress solution in the plastic region around an elliptic hole. This solution can be used in conjunction with the ductile fracture model proposed in [6]. Exact solutions for an elliptic hole embedded in a thermoelectric material have been presented in [7,8]. Paper [9] has found the stress distribution around an elliptic hole for non-local elastic material models. An elastic solution for a plate subject to biaxial loading and containing an inclined elliptic defect has been provided in [10]. The stress concentration generated by an elliptic hole in a functionally graded panel subjected to uniform tension has been analyzed in [11]. The effect of elastic anisotropy on the stress distribution around an elliptical crack has been revealed in [12], assuming that the medium is subjected to general loading at infinity. The influence of multiple pre-existing holes on the fracture process in granite subject to compressive loading has been studied in [13]. The stress concentration due to a circular hole in a plate subject to the uniaxial extension has been investigated in [14], using new constitutive equations with density-dependent elasticity moduli. This result has been extended to bi-axial deformation in [15]. A perturbation technique has been adopted in [16] for deriving analytical solutions for an elastic/plastic plate with a hole, including an elliptic hole. In particular, the elastic/plastic boundary has been found. Another analytical method for determining the elastic/plastic boundary in the vicinity of holes has been proposed in [17].
Another area of applications of the stress solutions near elliptic holes is that structural members contain such holes. Examples of such members have been considered in [18,19].
The present paper extends the solution [5] to an arbitrary piecewise linear yield criterion. Such yield criteria are widely used for describing metallic and nonmetallic materials [20,21,22]. The generalized method of moving coordinates is developed and used. The solution is semi-analytical. One needs to evaluate one ordinary integral to calculate the principal stresses at any point of the plastic region.

2. Basic Equations

Considered is plane-strain deformation of a material obeying an arbitrary piecewise linear yield criterion. Such yield criteria are represented by a linear function of the two principal stresses in a typical plane of flow. These stresses are denoted as σ 1 and σ 2 . The yield criterion can be represented as:
λ p + q = 1 .
Here, λ is a constitutive parameter and
p = σ 1 + σ 2 2   and   q = σ 1 σ 2 2 .
Note that Equation (1) represents the totality of piecewise linear yield criteria if the stresses are non-dimensionalized using an appropriate reference stress. For example, a standard form of the Mohr–Coulomb yield criterion is [23]:
p sin χ + q = k cos χ
where the cohesion k and angle of internal friction χ are constants. Equation (3) reduces to Equation (1) if λ = sin χ and the reference stress is k cos χ . It is also possible to choose without loss of generality that σ 1 > σ 2 . Then, q > 0 .
Let ψ be the inclination of the direction of the stress σ 1 relative to the x 1 axis of a Cartesian coordinate system ( x 1 ,   y 1 ) , measured anticlockwise. Then, the stress components referred to the Cartesian coordinates are:
σ x x = ( σ 1 + σ 2 ) 2 + ( σ 1 σ 2 ) 2 cos 2 ψ , σ y y = ( σ 1 + σ 2 ) 2 ( σ 1 σ 2 ) 2 cos 2 ψ ,   σ x y = ( σ 1 σ 2 ) 2 sin 2 ψ .
Equations (1), (2) and (4) combine to give:
σ x x = p ( 1 + λ cos 2 ψ ) + cos 2 ψ , σ y y = p ( 1 + λ cos 2 ψ ) cos 2 ψ ,   σ x y = ( 1 λ p ) sin 2 ψ .
Equation (1) and the equilibrium equations comprise a statically determinate system. This system, without the assumption that the yield criterion is piecewise linear, has been analyzed in [24]. In the case of Equation (1), this analysis provides the following characteristics:
d y 1 d x 1 = tan ( ψ ϕ )   and   d y 1 d x 1 = tan ( ψ + ϕ ) ,
which are termed the α and β lines, respectively, where
tan ϕ = 1 λ 1 + λ .
The relations along the characteristics are:
1 λ 2 d p + 2 ( 1 λ p ) d ψ = 0   and   1 λ 2 d p 2 ( 1 λ p ) d ψ = 0
on the α and β lines, respectively. It is seen from Equations (6) and (7) that the angle between the direction of the stress σ 1 and each characteristic line is constant. Therefore, Equation (8) can be rewritten as:
1 λ 2 d p + 2 ( 1 λ p ) d φ = 0   and   1 λ 2 d p 2 ( 1 λ p ) d φ = 0
where
φ = ψ ϕ
is the angle between the x axis and the α characteristic line measured from the axis anticlockwise. If both families of characteristics are curved, the equations in (9) can be immediately integrated to give:
P φ = 2 α cos ϕ φ 0   and   P + φ = 2 β cos ϕ + φ 0
where
P = 1 λ 2 2 λ ln ( 1 λ p 1 λ p 0 ) .
The constants φ 0 and p 0 have been introduced for further convenience. Note that Equation (12) reduces to P = ( p p 0 ) / 2 as λ 0 . Solving the equations in (11) for P and φ yields:
P = ( β α ) sin 2 ϕ   and   φ = ( β + α ) sin 2 ϕ + φ 0 .

3. Generalized Moving Coordinates

The moving coordinates ( x ¯ ,   y ¯ ) are used for solving boundary value problems in pressure-independent plasticity. This method has been described in many monographs, for example, [24]. In short, x ¯ and y ¯ are the coordinates of the point C under consideration, referred to axes passing through the origin O and parallel to the characteristic directions at C. The geometric interpretation of the moving coordinates is shown in Figure 1a. In the case of pressure-independent plasticity, the characteristic directions are orthogonal. Therefore, the moving coordinates are orthogonal. It is seen from Equations (6) and (7) that the characteristic directions are not orthogonal if λ 0 .
The geometric interpretation of the generalized moving coordinates, which are also denoted as ( x , y ) , is shown in Figure 1b. In this figure, ( x ,   y ) is a skew rectilinear coordinate system with the origin at O. The x axis is tangent to the base α characteristic line passing through O, and the y axis is tangent to the base β characteristic line passing through O. Let e x and e y be the unit base vectors of the ( x ,   y ) coordinate system. The unit vectors e α and e β are directed along the α and β characteristic lines at the point C under consideration. Therefore, the angle between the vectors e x and e α is φ . Then (Figure 1b):
e α e x = cos φ , e α e y = cos ( φ 2 ϕ ) , e β e x = cos ( φ + 2 ϕ ) , e β e y = cos φ , e y e x = cos 2 ϕ .
The position vector of any point can be represented as
R = x e x + y e y = x ¯ e α + y ¯ e β .
Using Equations (14) and (15), one can find the scalar products, R e x and R e y , as:
R e x = x + y cos 2 ϕ = x ¯ cos φ + y ¯ cos ( φ + 2 ϕ ) , R e y = x cos 2 ϕ + y = x ¯ cos ( φ 2 ϕ ) + y ¯ cos φ .
Solving these equations for x ¯ and y ¯ yields:
x ¯ = x [ cos φ cos 2 ϕ cos ( φ + 2 ϕ ) ] sin 2 2 ϕ + y [ cos φ cos 2 ϕ cos ( φ + 2 ϕ ) ] sin 2 2 ϕ , y ¯ = x [ cos φ cos 2 ϕ cos ( φ 2 ϕ ) ] sin 2 2 ϕ + y [ cos φ cos 2 ϕ cos ( φ 2 ϕ ) ] sin 2 2 ϕ .
Alternatively, the equations in (16) can be solved for x and y . As a result,
x = x ¯ [ cos φ cos 2 ϕ cos ( φ 2 ϕ ) ] sin 2 2 ϕ + y ¯ [ cos ( φ + 2 ϕ ) cos φ cos 2 ϕ ] sin 2 2 ϕ , y = x ¯ [ cos ( φ 2 ϕ ) cos φ cos 2 ϕ ] sin 2 2 ϕ + y ¯ [ cos φ cos 2 ϕ cos ( φ + 2 ϕ ) ] sin 2 2 ϕ .
Differentiating the first equation in (17) with respect to β and the second with respect to α leads to:
sin 2 2 ϕ x ¯ β = [ cos φ cos 2 ϕ cos ( φ + 2 ϕ ) ] x β + [ cos φ cos 2 ϕ cos ( φ + 2 ϕ ) ] y β { x [ sin φ cos 2 ϕ sin ( φ + 2 ϕ ) ] + y [ sin φ cos 2 ϕ sin ( φ + 2 ϕ ) ] } φ β , sin 2 2 ϕ y ¯ α = [ cos φ cos 2 ϕ cos ( φ 2 ϕ ) ] x α + [ cos φ cos 2 ϕ cos ( φ 2 ϕ ) ] y α { x [ sin φ cos 2 ϕ sin ( φ 2 ϕ ) ] + y [ sin φ cos 2 ϕ sin ( φ 2 ϕ ) ] } φ α .
It is seen from Equation (13) that φ / α = φ / β = sin 2 ϕ . Therefore, Equation (19) becomes:
x ¯ β = [ cos φ cos 2 ϕ cos ( φ + 2 ϕ ) ] sin 2 2 ϕ x β + [ cos φ cos 2 ϕ cos ( φ + 2 ϕ ) ] sin 2 2 ϕ y β { x [ sin φ cos 2 ϕ sin ( φ + 2 ϕ ) ] + y [ sin φ cos 2 ϕ sin ( φ + 2 ϕ ) ] } csc 2 ϕ , y ¯ α = [ cos φ cos 2 ϕ cos ( φ 2 ϕ ) ] sin 2 2 ϕ x α + [ cos φ cos 2 ϕ cos ( φ 2 ϕ ) ] sin 2 2 ϕ y α { x [ sin φ cos 2 ϕ sin ( φ 2 ϕ ) ] + y [ sin φ cos 2 ϕ sin ( φ 2 ϕ ) ] } csc 2 ϕ .
Using Equation (18), one can get after some trigonometry:
{ x [ sin φ cos 2 ϕ sin ( φ + 2 ϕ ) ] + y [ sin φ cos 2 ϕ sin ( φ + 2 ϕ ) ] } csc 2 ϕ = ( y ¯ + x ¯ cos 2 ϕ ) , { x [ sin φ cos 2 ϕ sin ( φ 2 ϕ ) ] + y [ sin φ cos 2 ϕ sin ( φ 2 ϕ ) ] } csc 2 ϕ = y ¯ cos 2 ϕ + x ¯ .
Equation (6) can be rewritten using Equation (10) as:
y 1 α = tan φ x 1 α   and   y 1 β = tan ( φ + 2 ϕ ) x 1 β .
One can situate the origin of the ( x 1 , y 1 ) coordinate system at the origin of the ( x , y ) coordinate system and direct the x 1 axis along the x axis (Figure 2). Then,
x = x 1 y 1 cot 2 ϕ   and   y = y 1 sin 2 ϕ .
Differentiating, one arrives at:
x α = x 1 α + y 1 α cot 2 ϕ , x β = x 1 β + y 1 β cot 2 ϕ , y α = 1 sin 2 ϕ y 1 α ,   y β = 1 sin 2 ϕ y 1 β .
Consider the terms of the equations in (20) containing the derivatives x / α , y / α , x / β , and y / β . Using Equations (22) and (24), one can represent these terms as:
[ cos φ cos 2 ϕ cos ( φ 2 ϕ ) ] sin 2 2 ϕ x α + [ cos φ cos 2 ϕ cos ( φ 2 ϕ ) ] sin 2 2 ϕ y α = { [ cos φ cos 2 ϕ cos ( φ 2 ϕ ) ] ( cot φ cos 2 ϕ ) + [ cos φ cos 2 ϕ cos ( φ 2 ϕ ) ] csc 2 ϕ } csc 2 2 ϕ y 1 α , [ cos φ cos 2 ϕ cos ( φ + 2 ϕ ) ] sin 2 2 ϕ x β + [ cos φ cos 2 ϕ cos ( φ + 2 ϕ ) ] sin 2 2 ϕ y β = { [ cos φ cos 2 ϕ cos ( φ + 2 ϕ ) ] [ cot ( φ + 2 ϕ ) cot 2 ϕ ] + [ cos φ cos 2 ϕ cos ( φ + 2 ϕ ) ] csc 2 ϕ } csc 2 2 ϕ y 1 α
Using trigonometric identities, one can show that the coefficient of y 1 / α vanishes in each of these equations. Therefore, Equations (20) and (22) combine to give:
x ¯ β ( y ¯ + x ¯ cos 2 ϕ ) = 0   and   y ¯ α + y ¯ cos 2 ϕ + x ¯ = 0 .
Introduce new variables X ¯ and Y ¯ as
x ¯ = X ¯ exp [ ( β α ) cos 2 ϕ ]   and   y ¯ = Y ¯ exp [ ( β α ) cos 2 ϕ ] .
Substituting Equation (27) into Equation (26) yields:
X ¯ β Y ¯ = 0   and   Y ¯ α + X ¯ = 0 .
It is seen from this equation that the quantities X ¯ and Y ¯ separately satisfy the equation of telegraphy. The latter is integrated by the method of Riemann.

4. Solution near a Traction-Free Surface

One of the principal stresses vanishes on any traction-free surface. The other principal stress and p are determined from Equations (1) and (2). It is necessary to consider two cases. One of these cases demands:
σ 1 = 0 ,   σ 2 = 2 1 + λ ,   and   p = 1 1 + λ
on a traction-free surface. The other case demands:
σ 2 = 0 ,   σ 1 = 2 1 λ ,   and   p = 1 1 λ
on a traction-free surface.
Consider Equation (29). The general structure of the solution is shown in Figure 3. It is required to find X ¯ and Y ¯ at a point C. BC is an α characteristic line, and AC is a β characteristic line. Riemann’s method for the telegraph equation results in
B C A ( G f α f G α ) d α + ( f G β G f β ) d β = 0 .
Here f is X ¯ or Y ¯ and.
G ( α C , β C , α , β ) J 0 [ 2 ( α C α ) ( β C β ) ] .
Also, J 0 ( z ) is the Bessel function of the first kind of zero order, α C is the value of α at the point C, and β C is the value of β at the point C. Note that
J 0 ( 0 ) = 1 .   and   J 1 ( 0 ) = 0
where J 1 ( z ) = d J 0 ( z ) / d z is the Bessel function of the first kind of first order.
Since d β = 0 on BC and d α = 0 on AC, Equation (31) can be rewritten as:
B C A ( G f α f G α ) d α + ( f G β G f β ) d β = α B α C ( G f α f G α ) d α + β C β A ( f G β G f β ) d β + A B ( G f α f G α ) d α + ( f G β G f β ) d β = 0 .
Moreover, β = β C on BC and α = α C on AC. Therefore, z = 0 on these lines. Then, it follows from Equations (33) and (34) that:
f C = f B + f A 2 1 2 A B ( G f α f G α ) d α + ( f G β G f β ) d β .
Here f A , f B , and f C are the value of f at A, B, and C, respectively.
Put p 0 = 1 / ( 1 + λ ) in (12). Then, it follows from Equations (12), (13), and (29) that:
P = 0 ,   d P = 0 ,   d β = d φ 2 sin 2 ϕ ,   and   d α = d φ 2 sin 2 ϕ
on AB.
Using Equation (36), one can rewrite the integral in Equation (35) as
A B ( G f α f G α ) d α + ( f G β G f β ) d β = = 1 2 sin 2 ϕ A B ( G f α f G α + f G β G f β ) d φ .
It follows from Equation (13) that:
α = sin 2 ϕ ( φ P )   and   β = sin 2 ϕ ( φ + P ) .
Equations (37) and (38) combine to give:
A B ( G f α f G α ) d α + ( f G β G f β ) d β = A B ( f G P G f P ) d φ .
The derivative f / P can be eliminated using Equation (28). In particular, it follows from Equations (28) and (38) that:
Y ¯ P = Y ¯ φ csc 2 ϕ Y ¯ α = Y ¯ φ + X ¯ csc 2 ϕ , X ¯ P = csc 2 ϕ X ¯ β X ¯ φ = Y ¯ csc 2 ϕ X ¯ φ .
Moreover, by replacing α and β in Equation (32) with φ and P using Equation (13), one gets:
G ( φ C , P C , φ , P ) J 0 ( Z )
where φ C is the value of φ at C, P C is the value of P at C, and
Z = csc 2 ϕ ( φ φ C ) 2 ( P P C ) 2 .
It follows from Equations (41) and (42) that:
G P = J 1 ( Z ) ( P P C ) Z sin 2 2 ϕ   and   G φ = J 1 ( Z ) ( φ φ C ) Z sin 2 2 ϕ .
Substituting Equation (40) into Equation (39) and integrating by parts, one arrives at:
A B ( X ¯ G P G X ¯ P ) d φ = A B [ X ¯ ( G P G φ ) G Y ¯ csc 2 ϕ ] d φ + X ¯ B X ¯ A , A B ( Y ¯ G P G Y ¯ P ) d φ = A B [ Y ¯ ( G P + G φ ) G X ¯ csc 2 ϕ ] d φ Y ¯ B + Y ¯ A .
Eliminating here the derivatives G / P and G / φ using Equation (43) yields:
A B ( X ¯ G P G X ¯ P ) d φ = csc 2 ϕ A B [ X ¯ J 1 ( Z ) ( φ φ C P C ) Z sin 2 ϕ J 0 ( Z ) Y ¯ ] d φ + X ¯ B X ¯ A , A B ( Y ¯ G P G Y ¯ P ) d φ = csc 2 ϕ A B [ Y ¯ J 1 ( Z ) ( φ φ C + P C ) Z sin 2 ϕ + J 0 ( Z ) X ¯ ] d φ Y ¯ B + Y ¯ A .
Substituting Equation (45) into Equation (35) leads to:
X ¯ C = X ¯ A csc 2 ϕ 2 φ A φ B [ X ¯ J 1 ( Z ) ( φ φ C P C ) Z sin 2 ϕ J 0 ( Z ) Y ¯ ] d φ , Y ¯ C = Y ¯ B + csc 2 ϕ 2 φ A φ B [ Y ¯ J 1 ( Z ) ( φ φ C + P C ) Z sin 2 ϕ + J 0 ( Z ) X ¯ ] d φ .
The coefficients of X ¯ and Y ¯ in the integrands are known functions of φ due to Equations (36) and (42).
Consider Equation (30). In this case, AC is an α characteristic line, and BC is a β characteristic line (Figure 3). Putting p 0 = 1 / ( 1 + λ ) in Equation (12) and repeating the line of reasoning that has led to Equation (46) in the previous case, one gets:
X ¯ C = X ¯ B + csc 2 ϕ 2 φ A φ B [ X ¯ J 1 ( Z ) ( φ φ C P C ) Z sin 2 ϕ J 0 ( Z ) Y ¯ ] d φ , Y ¯ C = Y ¯ A csc 2 ϕ 2 φ A φ B [ Y ¯ J 1 ( Z ) ( φ φ C + P C ) Z sin 2 ϕ + J 0 ( Z ) X ¯ ] d φ .
It follows from Equation (11) that:
P C = φ A φ B 2   and   φ C = φ B + φ A 2
if Equation (29) is valid and
P C = φ B φ A 2   and   φ C = φ B + φ A 2
Equation (30) is valid. Using Equations (48) and (49), one can eliminate φ C and P C in Equations (46) and (47).
The free surface geometry determines the quantities X ¯ and Y ¯ involved in Equations (46) and (47). Since P = 0 on the free surface, it is seen from Equations (13) and (27) that:
X ¯ = x ¯   and   Y ¯ = y ¯
on AB. One of the principal stress directions is orthogonal to the free surface. Therefore, replacing φ with ψ in Equations (46) and (47) is convenient. It is seen from Equation (10) that d ψ = d φ . Then, using Equations (48)–(50), one rewrites Equations (46) and (47) as:
X ¯ C = x ¯ A csc 2 ϕ 2 ψ A ψ B [ x ¯ J 1 ( Z ) ( ψ ψ B ) Z sin 2 ϕ J 0 ( Z ) y ¯ ] d ψ , Y ¯ C = y ¯ B + csc 2 ϕ 2 ψ A ψ B [ y ¯ J 1 ( Z ) ( ψ ψ A ) Z sin 2 ϕ + J 0 ( Z ) x ¯ ] d ψ ,
and
X ¯ C = x ¯ B + csc 2 ϕ 2 ψ A ψ B [ x ¯ J 1 ( Z ) ( ψ ψ B ) Z sin 2 ϕ J 0 ( Z ) y ¯ ] d ψ , Y ¯ C = y ¯ A csc 2 ϕ 2 ψ A ψ B [ y ¯ J 1 ( Z ) ( ψ ψ A ) Z sin 2 ϕ + J 0 ( Z ) x ¯ ] d ψ ,
respectively.
Equations (51) and (52) coincide with the solution [5] for pressure-independent material (i.e., ϕ = π / 4 ).

5. Numerical Example

The distributions of x ¯ and y ¯ along AB are required to evaluate the integrals in Equations (51) and (52). The free surface shape determines these distributions. The numerical solution below is for an elliptic hole. Its equation can be written as:
y 1 2 + x 1 2 a 2 = 1   or   ( sin 2 θ + cos 2 θ a 2 ) ρ 2 = 1
where x 1 = ρ cos θ and y 1 = ρ sin θ . Equations (17) and (18) have been derived upon the assumption that the x axis coincides with the x 1 axis. Therefore, the origin of the ( x , y ) coordinate system should be situated at a point of the ellipse where the tangent to the α is parallel to the x 1 axis. This condition is equivalent to:
ψ = ϕ
at the origin of the ( x , y ) coordinate system. Let θ 0 be the value of θ at which this condition is satisfied. The components of the unit vector orthogonal to the ellipse are readily determined from Equation (53). As a result,
a 2 tan θ = tan ψ .
Then, the condition Equation (54) yields:
a 2 tan θ 0 = tan ϕ .
The ellipse, ( x ,   y ) coordinate system, and ( x 1 ,   y 1 ) coordinate system are shown in Figure 4. The unit base vectors of the ( x ,   y ) coordinate system are denoted as e x and e y , and the unit base vectors of the ( x 1 ,   y 1 ) coordinate system as i and j. The position vector of a generic point is:
R = R 0 + r .
Using Equation (56), one can represent the vector R 0 as:
R 0 = x 0 i + y 0 j , x 0 = a 2 cos ϕ 1 ( 1 a 2 ) cos 2 ϕ , y 0 = sin ϕ 1 ( 1 a 2 ) cos 2 ϕ .
It follows from the geometry of Figure 4 that:
e x i = 1 , e x j = 0 , e y i = cos 2 ϕ , e y j = sin 2 ϕ .
Using Equations (58) and (59), one can find the scalar products R i and R j from Equation (57) as:
x 1 = x 0 + x + y cos 2 ϕ   and   y 1 = y 0 + y sin 2 ϕ .
Solving these equations for x and y yields:
x = x 1 y 1 cot 2 ϕ x 0 + y 0 cot 2 ϕ   and   y = y 1 y 0 sin 2 ϕ .
It is convenient to rewrite these equations in terms of ρ and θ as
x = ρ ( cos θ sin θ cot 2 ϕ ) x 0 + y 0 cot 2 ϕ   and   y = ρ sin θ y 0 sin 2 ϕ .
One can eliminate r in these equations using Equation (54). As a result,
x = a ( cos θ sin θ cot 2 ϕ ) a 2 sin 2 θ + cos 2 θ x 0 + y 0 cot 2 ϕ , y = a sin θ sin 2 ϕ a 2 sin 2 θ + cos 2 θ y 0 sin 2 ϕ .
Equation (55) allows for x and y to be calculated from Equation (63) as functions of ψ . Then, Equations (10) and (17) supply x ¯ and y ¯ as functions of ψ . These functions have been used in Equation (51) to calculate X ¯ C and Y ¯ C . Having X ¯ C and Y ¯ C , one can calculate the x 1 and y 1 coordinates of point C using Equations (18), (27), and (60). The principal stresses at this point are found from Equations (1), (2), (12), and (48) (or (49)). The numerical results below are for Equation (29). Figure 5 depicts the variation of the stress σ 1 along the x 1 axis for several values of λ and a. Figure 6 shows the variation of the stress σ 2 . In these figures, s is the distance from the hole surface, s = x 1 a . Figure 5 and Figure 6 reveal a significant effect of λ on the stresses. In particular, the effect of a is invisible in these figures. However, the hole shape does affect the stresses near the surface. Figure 7 depicts the variation of the principal stresses for several a-values at λ = 1 / 4 . It is seen from this figure that the value of | σ 1 | near a circular hole ( a = 1 ) is larger than its value near an elliptic hole. However, elliptical holes produce higher values of | σ 1 | than the corresponding circular hole at some distance from the hole’s surface. The same tendency occurs with the stress | σ 2 | .

6. Conclusions

From this work, the following conclusions can be drawn:
(1)
The generalized method of moving coordinates is efficient for determining stress fields near traction-free surfaces. In particular, determining the stress components at any point of the plastic region requires evaluating one ordinary integral in Equations (51) or (52).
(2)
The solutions in Equations (51) and (52) are independent of other boundary conditions, though the shape of the plastic region is.
(3)
The solution found is identical to the solution in [5] if λ = 0 in Equation (1).
(4)
The parameter λ involved in Equation (1) has a much greater effect on the stress field than the aspect ratio (Figure 5 and Figure 6). However, it is a consequence of the yield criterion. For a given material, the aspect ratio affects the stress distribution in a small region near the hole (Figure 7).
(5)
The solution found is for perfectly plastic material. However, slip-line solutions for the class of problems considered agree well with experimental results for strain-hardening materials [25].

Author Contributions

Conceptualization, S.A. and S.-N.T.; formal analysis, M.R.; writing—review and editing, S.A. and S.-N.T. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by Department of Mechanical and Electro-Mechanical Engineering, National Sun Yat-sen University, Taiwan, under project No. 11TA303. The APC was also funded by the same department and university.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

This paper has been supported by the RUDN University Strategic Academic Leadership Program.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Williams, M.L. On the stress distribution at the base of a stationary crack. J. Appl. Mech. 1956, 24, 109–114. [Google Scholar] [CrossRef]
  2. Rice, J.R. A path independent integral and the approximate analysis of strain concentration by notches and cracks. J. Appl. Mech. 1968, 35, 379–386. [Google Scholar] [CrossRef]
  3. Li, J.; Hancock, J.W. Mode I and mixed mode fields with incomplete crack tip plasticity. Int. J. Solids Struct. 1999, 36, 711–725. [Google Scholar] [CrossRef]
  4. Papanastasiou, P.; Durban, D. Singular plastic fields in non-associative pressure sensitive solids. Int. J. Solids Struct. 2001, 38, 1539–1550. [Google Scholar] [CrossRef]
  5. Thomason, P.F. Riemann-Integral solutions for the plastic slip-line fields around elliptical holes. J. Appl. Mech. 1978, 45, 678–679. [Google Scholar] [CrossRef]
  6. Thomason, P.F. Ductile fracture and the stability of incompressible plasticity in the presence of microvoids. Acta Metall. 1981, 29, 763–777. [Google Scholar] [CrossRef]
  7. Zhang, A.B.; Wang, B.L. Explicit solutions of an elliptic hole or a crack problem in thermoelectric materials. Eng. Fract. Mech. 2016, 151, 11–21. [Google Scholar] [CrossRef]
  8. Zhang, A.B.; Wang, B.L.; Wang, J.; Du, J.K. Two-dimensional problem of thermoelectric materials with an elliptic hole or a rigid inclusion. Int. J. Therm. Sci. 2017, 117, 184–195. [Google Scholar] [CrossRef]
  9. Tuna, M.; Trovalusci, P. Stress distribution around an elliptic hole in a plate with ‘implicit’ and ‘explicit’ non-local models. Compos. Struct. 2021, 256, 113003. [Google Scholar] [CrossRef]
  10. Ostsemin, A.A.; Utkin, P.B. Stress-strain state and stress intensity factors of an inclined elliptic defect in a plate under biaxial loading. Mech. Solids 2010, 45, 214–225. [Google Scholar] [CrossRef]
  11. Wang, W.; Yuan, H.; Li, X.; Shi, P. Stress concentration and damage factor due to central elliptical hole in functionally graded panels subjected to uniform tensile traction. Materials 2019, 12, 422. [Google Scholar] [CrossRef]
  12. Willis, J.R. The stress field around an elliptical crack in an anisotropic elastic medium. Int. J. Eng. Sci. 1968, 6, 253–263. [Google Scholar] [CrossRef]
  13. Huang, Y.-H.; Yang, S.-Q.; Tian, W.-L. Cracking process of a granite specimen that contains multiple pre-existing holes under uniaxial compression. Fatigue Fract. Eng. Mater. Struct. 2019, 42, 1341–1356. [Google Scholar] [CrossRef]
  14. Murru, P.T.; Rajagopal, K.R. Stress concentration due to the presence of a hole within the context of elastic bodies. Mat. Des. Process. Commun. 2021, 3, e219. [Google Scholar] [CrossRef]
  15. Murru, P.T.; Rajagopal, K.R. Stress concentration due to the bi-axial deformation of a plate of a porous elastic body with a hole. Z. Angew. Math. Mech. 2021, 101, e202100103. [Google Scholar] [CrossRef]
  16. Meyer, M.; Sayir, M.B. The elasto-plastic plate with a hole: Analytical solutions derived by singular perturbations. Z. Angew. Math. Phys. 1995, 46, 427–445. [Google Scholar]
  17. Yaocai Ma, Y.; Lu, A.; Cai, H. Analytical method for determining the elastoplastic interface of a circular hole subjected to biaxial tension-compression loads. Mech. Based Des. Struct. Mach. 2022, 50, 3206–3223. [Google Scholar] [CrossRef]
  18. Li, Z.; Mantell, S.C.; Davidson, J.H. Mechanical analysis of streamlined tubes with non-uniform wall thickness for heat exchangers. J. Strain Anal. Eng. Des. 2005, 40, 275–285. [Google Scholar] [CrossRef]
  19. Jafari Fesharaki, J.; Roghani, M. Mechanical behavior and optimization of functionally graded hollow cylinder with an elliptic hole. Mech. Adv. Compos. Struct. 2020, 7, 189–201. [Google Scholar] [CrossRef]
  20. Druyanov, B. Technological Mechanics of Porous Bodies; Clarendon Press: New York, NY, USA, 1993. [Google Scholar]
  21. Cox, G.M.; Thamwattana, N.; McCue, S.W.; Hill, J.M. Coulomb–Mohr Granular Materials: Quasi-static flows and the highly frictional limit. ASME. Appl. Mech. Rev. 2008, 61, 060802. [Google Scholar] [CrossRef]
  22. Burenin, A.A.; Tkacheva, A.V. Piecewise linear plastic potentials as a tool for calculating plane transient temperature stresses. Mech. Solids 2020, 55, 791–799. [Google Scholar] [CrossRef]
  23. Spencer, A.J.M. A theory of the kinematics of ideal soils under plane strain conditions. J. Mech. Phys. Solids 1964, 12, 337–351. [Google Scholar] [CrossRef]
  24. Hill, R. The Mathematical Theory of Plasticity; Clarendon Press: Oxford, UK, 1950. [Google Scholar]
  25. Bastawros, A.; Kyung-Suk, K. Experimental analysis of near-crack-tip plastic flow and deformation characteristics (I): Polycrystalline aluminum. J. Mech. Phys. Solids 2000, 48, 67–98. [Google Scholar] [CrossRef]
Figure 1. Moving coordinates: (a) pressure-independent yield criterion, (b) pressure-dependent yield criterion.
Figure 1. Moving coordinates: (a) pressure-independent yield criterion, (b) pressure-dependent yield criterion.
Materials 15 06266 g001
Figure 2. ( x 1 ,   y 1 ) coordinate system chosen for deriving Equation (28).
Figure 2. ( x 1 ,   y 1 ) coordinate system chosen for deriving Equation (28).
Materials 15 06266 g002
Figure 3. Illustration of the general boundary value problem.
Figure 3. Illustration of the general boundary value problem.
Materials 15 06266 g003
Figure 4. Illustration of the boundary value problem for an elliptic hole.
Figure 4. Illustration of the boundary value problem for an elliptic hole.
Materials 15 06266 g004
Figure 5. Dependence of the principal stress σ1 on the distance from the free surface along the x1-axis for several λ-values: (a) a = 1 (circular hole), (b) a = 1/2, (c) a = 1/4, (d) a = 1/6.
Figure 5. Dependence of the principal stress σ1 on the distance from the free surface along the x1-axis for several λ-values: (a) a = 1 (circular hole), (b) a = 1/2, (c) a = 1/4, (d) a = 1/6.
Materials 15 06266 g005
Figure 6. Dependence of the principal stress σ2 on the distance from the free surface along the x1-axis for several λ-values: (a) a = 1 (circular hole), (b) a = 1/2, (c) a = 1/4, (d) a = 1/6.
Figure 6. Dependence of the principal stress σ2 on the distance from the free surface along the x1-axis for several λ-values: (a) a = 1 (circular hole), (b) a = 1/2, (c) a = 1/4, (d) a = 1/6.
Materials 15 06266 g006
Figure 7. Dependence of the principal stress on the distance from the free surface along the x1-axis for several a-values at λ = −1/4: (a) σ1-stress, (b) σ2-stress.
Figure 7. Dependence of the principal stress on the distance from the free surface along the x1-axis for several a-values at λ = −1/4: (a) σ1-stress, (b) σ2-stress.
Materials 15 06266 g007
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Alexandrov, S.; Rynkovskaya, M.; Tsai, S.-N. Application of the Generalized Method of Moving Coordinates to Calculating Stress Fields near an Elliptical Hole. Materials 2022, 15, 6266. https://doi.org/10.3390/ma15186266

AMA Style

Alexandrov S, Rynkovskaya M, Tsai S-N. Application of the Generalized Method of Moving Coordinates to Calculating Stress Fields near an Elliptical Hole. Materials. 2022; 15(18):6266. https://doi.org/10.3390/ma15186266

Chicago/Turabian Style

Alexandrov, Sergei, Marina Rynkovskaya, and Shang-Nan Tsai. 2022. "Application of the Generalized Method of Moving Coordinates to Calculating Stress Fields near an Elliptical Hole" Materials 15, no. 18: 6266. https://doi.org/10.3390/ma15186266

APA Style

Alexandrov, S., Rynkovskaya, M., & Tsai, S. -N. (2022). Application of the Generalized Method of Moving Coordinates to Calculating Stress Fields near an Elliptical Hole. Materials, 15(18), 6266. https://doi.org/10.3390/ma15186266

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop