Jtam-A4.dvi


JOURNAL OF THEORETICAL

AND APPLIED MECHANICS

50, 4, pp. 1097-1108, Warsaw 2012

50th Anniversary of JTAM

MATHEMATICAL MODELLING OF A DISC WEAKENED BY

AN ECCENTRIC CIRCULAR HOLE

Milan Bižić, Dragan Petrović

University of Kragujevac, Faculty of Mechanical Engineering Kraljevo, Serbia

e-mail: bizic.m@mfkv.kg.ac.rs; petrovic.d@mfkv.kg.ac.rs

Dragan Pančić

Wagon Factory Kraljevo, Serbia; e-mail: pancic.dr@gmail.com

Zoran –Dinović

Vienna University of Technology, Institute of Sensor and Actuator Systems, Austria;

Integrated Microsystems Austria GmbH, Austria; e-mail: djinovic@ima-mst.at

The task of this paper is identification of stresses in a homogeneous isotropic disc weakened
with an eccentric circular holewhich is loadedbypressure in the internal contour of the hole.
By application of a complex variablemethod, themathematicalmodel that allows complete
analytical solution of stresses of the disc is formed. The methodology can be applied for
the solution of any disc weakenedwith an eccentric circular hole. The comparative analysis
has shown a high accuracy of analytically obtained results with FEM results obtained by
calculations inANSYS 12 software package.The application of the results of this paper is of
great importance for quality design and optimization of thin-walled structures of disc type
weakened by a circular hole.

Key words: modelling, disc, eccentric hole

1. Introduction

The papers by Timoshenko and Goodier (1951), Timoshenko and Woinowsky-Kreiger (1959)
provide methods for solving some typical problems of the theory of elasticity. These methods
are applied broadly to some practical engineering problems, particularly those with lower ap-
proximation, and give reasonably accurate results. However, when there is a case where the
problem solving significantly differs from the given theoretical model, it is necessary to make
quite a rough approximation,which is the reasonwhy the obtained results are not sufficiently ac-
curate. A particular problem is the calculation of machine elements and structures with sudden
change of geometry, the existence of openings, sharp corners, etc. At these places, there is the
phenomenon of stress concentration which occupies a significant place in the study of problems
of the theory of elasticity. The data from technical practice show that frequent breakdowns and
accidents happen on such machines or structural elements and that the consequences are often
tragic human victims and enormousmaterial damage. Inmany cases, the cause is an incomplete
and inaccurate identification of stresses and strains.

In studiesof thisphenomenon, inmost cases, the thinplateweakenedbyacircularhole loaded
with a certain type of loading is analysed (Bakhshandeh et al., 2008; Bizic and Petrovic, 2011;
Bojic et al., 2010; Chandrashekhara and Muthanna, 1978; Chen and Archer, 1989; Mizushima
andHamada, 1983; Troyani et al., 2002;Wang, 2004; Yang et al., 2010; Zhang and Shen, 2011).
Also, the analysis may include determination of stresses in a thin plate in the case when there
are two adjacent circular holes or more holes (Arshadnejad et al., 2009; Chen et al., 2000;



1098 M. Bižić et al.

Kratochvil andBecker, 2011;Wu andMarkenscoff, 1996). Such cases are very frequently present
in constructionof railwayand roadvehicles, vessels, aircrafts, civil engineeringmachinery,mining
and transportation machinery, cranes, tooling machines, steel structures andmany others. The
holes exist on them for various reasons such as construction requirements, optimization of the
structure, reduction of self-weight, esthetic reasons, and so on. In all cases, the problem is the
accurate determination of stresses in some points of the element in loaded plates weakened by
circular holes, in which the stress concentration is present. One of directions for solving this
problem, whose theoretical basis are defined in Bower (2010), Lu (1995), Muskhelishvili (1963),
Savin (1961), is based on the use of complex functions and complex analysis (Mitrinovic, 1981).
The basic equations of the theory of elasticity and stress functions are expressed in a complex
form,and for their solution themostwidelyusedmethod is themethodof conformalmapping.By
applying the complex variable method, it is possible to determine theoretically the exact stress
state of the observed element in a plate weakened by a hole (Huan-chun et al., 1987; Simha
and Mohapatra, 1998). In the design and calculation phase of the mentioned machinery and
structures, this allows very accurate theoretical determination of stresses, which are impossible
to be found by conventional procedures. The papers by Batra andNie (2010). Radi and Strozzi
(2009) deal with the analysis of mechanical elements similar to the disk weakened with an
eccentric circular hole. In line with these researches, a very interesting problem ismathematical
modelling and identification of stresses in a discweakenedwith an eccentric circular holewhen it
is loaded by the pressure in the internal contour of the hole. Thiswasmotivation for the research
published in this paper,where themathematicalmodel that allows analytical solution of stresses
is obtained by using the complex variable method (CVM). In order to verify analytically the
obtained results, the stresses are also determined by the finite element method (FEM).

2. Theoretical formulation

For obtaining an analytical solution of the problem of stresses in a disc weakened with an
eccentric circular hole, it is necessary to start with the theoretical formulation of the plane
stress condition. If the thin plate is loaded by forces that are evenly spaced along its thickness
and act in parallel to its base, there is a plane stress condition. The plane stress condition is
defined by the following four groups of equations (Timoshenko andGoodier, 1951).

The first group of equations link the stresses and volume forces

∂σx
∂x
+
∂τxy
∂y
+X =0

∂τxy
∂x
+
∂σy
∂y
+Y =0 (2.1)

where: σx,σy – components of normal stresses, τxy – shear stress, X, Y – components of volume
forces.

The second group of equations is the relation between stress and strain

σx = λε+2µ
∂u

∂x
σy = λε+2µ

∂v

∂y
τxy = µ

(∂u

∂y
+
∂v

∂x

)

(2.2)

where: ε – surface deformation, λ, µ – Lame’s constants, that are

λ =
νE

(1−2ν)(1+ν)
µ =

E

2(1+ν)

and E – modulus of elasticity, ν – Poisson’s ratio, u, v – displacements in the direction of the
coordinate axes x and y, respectively.



Mathematical modelling of a disc weakened... 1099

The third group of equations are compatibility conditions, where solutions have a unique ID
only when they satisfy the conditions of compatibility

∂2εx
∂y2
+
∂2εy
∂x2
=
∂2γxy
∂x∂y

(2.3)

where: εx, εy – deformational components, γxy =(∂u/∂y)+(∂v)/(∂x) – slip component.
Finally, the fourth group of equation are boundary conditions

Xn = lσx+mτxy Yn = lτxy+mσy (2.4)

where: Xn, Yn – components of the vector of external forces, l = cos(n,x), m = cos(n,y) –
direction cosine, n – vector of the external normal to the contour.
Therefore, the plane stress condition is completely defined by equations (2.1)-(2.4). In one

point with the stresses σx, σy and τxy, it is always possible to find a coordinate system in
which the normal stresses have extreme values and inwhich the shear stresses equal zero. These
stresses are named the principal stresses and are defined by the following equations

σ1,2 =
σx+σy
2
±

√

(σx−σy
2

)2

+ τ2xy (2.5)

2.1. Plane stress condition expressed through complex potentials

The problem solving in the plane stress condition is reduced to finding a stress function
that will uniquely determine stresses and deformations, while equations (2.1)-(2.4) will match
with it. Determination of stresses in the disc weakened by an eccentric circular hole is this type
of problem in which the volume forces can be neglected. As a consequence, there is always a
function U(x,y) throughwhich it is possible to express the stresses by the following expressions

σx =
∂2U(x,y)

∂2y
σy =

∂2U(x,y)

∂2x
τxy =−

∂2U(x,y)

∂x∂y
(2.6)

The function U(x,y) is called the stress function and it has to match with the following bihar-
monic equation

∂4U

∂x4
+2

∂4U

∂x2∂y2
+
∂4U

∂y4
=0 (2.7)

Solving equation (2.7) and determining the stress function in real form is often very complex,
and for a number of problems is practically impossible. One method of solving this problem
is the transition into the complex area and solution of the problem in a complex form. The
stress function, which is a function of two independent variables x and y, is expressed through
two functions of one complex variable. In this way, the problem of solving a single function of
two independent variables reduces the problem of two complex functions of one independent
variable. These complex functions are φ(z) and ψ(z), and they are called complex potentials.
Equations of the plane stress condition expressed trough complex potentials are (Muskheli-

shvili, 1963)

σx+σy =2[φ(z)+φ(z)] = 4Reφ(z)

σy −σx+2iτxy =2[zφ
′(z)+ψ(z)]

(2.8)

where: Re – real part, z – conjugated complex number.
The complex potentials are determined through the conditions made at the contour. In

solving the problem of determining the stresses of the disc weakened by an eccentric circular
hole, several methods for determining the complex potentials can be applied. In this paper, the
method of power series is applied, but before that the conformal mapping is defined.



1100 M. Bižić et al.

2.2. Conformal mapping

If z = x + iy and ζ = ζ + iη are two complex variables that are linked with the relation
z =ω(ζ), where ω(ζ) is an unambiguous analytical function in the area of Σ in the plane area
of change ζ, then eachmapping (Fig. 1) that is applied by using these functions in which values
of the angles are preserved, is called a conformal mapping (Muskhelishvili, 1963).

Fig. 1. Conformalmapping

In other words, the task is to find a function for themapping inwhich the angle between the
two curves in the plane z will be copiedwithout changes in the angle between the corresponding
curves in the plane ζ (α = β). There are advanced methods for formation of the conformal
mapping function ω(ζ). In this paper, the ready-made function of conformal mapping is used.

2.3. Stresses in a mirrored area

Some area S in the plane z is mapped into a circular area Σ in the plane ζ (Fig. 1). In a
new field Σ in the plane ζ, the polar coordinates θ and ρ are introduced by using the following
relation

ζ = ρeiθ (2.9)

The circles ρ = const and θ = const of the field Σ in the plane ζ match curves which are also
marked with ρ = const and θ = const. Themapping is done by using the following analytical
function

x+iy = ω(ρeiθ) (2.10)

So, the lines ρ = const and θ = const are coordinate lines in the mirrored area Σ, and in
the plane z they intersect at right angles. If through the point on the plane z pass the curves
ρ = const and θ = const (Fig. 2), which are mutually orthogonal, then the angle between
the tangent to the curve θ = const in the direction of increasing ρ and x axis can be marked
with α.
The link between the stress ratio in Cartesian and polar coordinates is as follows:

σρ+σθ = σx+σy

σθ −σρ+2iτρθ =e
2iα(σy−σx+2iτxy)

(2.11)

where: σρ – normal component of stress on the curve ρ = const, σθ – normal component of
stress on the curve θ = const, τρθ – tangential component of stress on both curves.
The stresses that are expressed in the polar coordinate system through the complex potential

in the newmirrored area are defined by the expressions

σρ+σθ =4Reφ(z)

σθ −σρ+2iτρθ =
2ζ2

ρ2ω′(ζ)
[ω(ζ)φ′(ζ)+ω′(ζ)ψ(ζ)]

(2.12)



Mathematical modelling of a disc weakened... 1101

Fig. 2. The link between stress in Cartesian and polar coordinates in the mirrored area

3. Determination of complex potentials for a disc weakened with an eccentric

circular hole

The disc is limited by two circles whose centers are shifted by the eccentricity e, as shown in
Fig. 3.

Fig. 3. The disc weakened by an eccentric circular hole

The condition is that themainmoment and themain vector on the contours equal zero. The
function of conformal mapping by which the disc mirrors on a circular ring is

z = ω(ζ)= c
ζ +1

ζ −1
(3.1)

where: c – real constant determined from the expression

c =

√

e4+R41+R
4
2−2(R1R2)

2−2e2(R21+R
2
2)

2e
(3.2)

The outer contour of the disc of radius R1 is mapped to the inner contour of the circular ring
of radius ρ0, according to the equation

ρ0 =
c+
√

c2+R21

R1
(3.3)

The inner contour of the disc of radius R2 is mapped to the outer contour of the circular ring
of radius R1, according to the equation

ρ1 =
c+
√

c2+R22

R2
(3.4)



1102 M. Bižić et al.

The conformal mapping is defined by the following relation

x2+
(

y− c
ζζ +1

ζζ −1

)2

=

∣

∣

∣

∣

2c
√

ζζ

ζζ −1

∣

∣

∣

∣

2

(3.5)

From relation (3.5), it is noticed that the concentric circles ζζ = ρ2 correspond to the circles in
S area with the center on the y axis shifted for the value d of radius r. The two parameters
(d and r) are concluded from relation (3.5) as follows

d = c
ρ2+1

ρ2−1
r =

2cρ

ρ2−1
(3.6)

The expressions for the stresses in polar coordinates are obtained by solving the system of
equations (2.12)

σθ =2Reφ(ζ)+
1

2
ReΩ(ζ) σρ =2Reφ(ζ)−

1

2
ReΩ(ζ)

τρθ =
1

2
ImΩ(ζ)

(3.7)

In equations (3.7), Ω(ζ) is a complex function that is defined by the following equation

Ω(ζ)=
2ζ2

ρ2ω′(ζ)
[ω(ζ)φ′(ζ)+ω′(ζ)ψ(ζ)] (3.8)

The complex potentials are

φ(ζ)= C0+(ζ −1)
2
(

C1+
C2
ζ2

)

ψ(ζ)=
ρ20(ζ −1)

2

(ρ20− ζ)
2

[

φ(ζ)+C0+
(ρ20
ζ
−1
)2(

C1+
C2ζ

2

ρ40

)]

+
(ρ20+ ζ)(ζ −1)

2

2(ρ20− ζ)
φ′(ζ)

(3.9)

The real constants are determined in the following way

C0 =
1

2
C1
[

2(ρ0ρ1)
2+2(ρ20+ρ

2
1)−2(ρ

2
0+ρ

2
1)
2
−

pρ21
2(ρ21−ρ

2
0)

]

C1 =−
pρ21ρ

2
0

(ρ21−ρ
2
0)(ρ

2
1ρ
2
0)(ρ

2
0+ρ

2
1)−4(ρ0ρ1)

2+ρ20ρ
2
1

C2 =−C1(ρ0ρ1)
2

(3.10)

In relations (3.10), p is the pressure that operates in the inner part of the contour of the
eccentric hole.

4. Determination of stress using CVM

The previously defined theoretical equations were applied to the concrete example of the disc of
thickness h =5mm, radius R1 =50mm, weakened by an eccentric circular hole of eccentricity
e =30mm, radius R2 =10mm, loaded on the inside contourwith the pressure p =1.0kN/cm

2.
By variation of the polar coordinates ρ and θ in the mirrored area, it is possible to determine
stress at any point of the disc. In this paper using the CVM, the specific numerical values of
stresses σρ and σθ and principal stresses σ1 and σ2 on the outer contour of the disc and at
the contour of the eccentric hole, as well as the stresses σx and σy at the intersection along to



Mathematical modelling of a disc weakened... 1103

the vertical axis of symmetry of the disc are determined. The numerical values of stresses σρ
and σθ are determined in function of the coordinates ρ1 and ρ0 in the outer contour of the disc
and the contour of the eccentric hole in themapped field calculated according to equations (3.3)
and (3.4), for 36 points that are defined by the angle θ that ranges from 0◦-360◦ with a step
of 10◦. The numerical values of stresses σx and σy are calculated in 22 points on the interval
from ρ0 to ρ1 (y =0.5cm to y =10.5cm) with step

∆ =
ρ1−ρ0
10

(4.1)

Therefore, the stresses σx and σy are determined for values of the coordinate ρ with the step ∆
according to equation (4.1) and for the angles θ = 0◦ and θ = 180◦. This means that the
numerical values of stresses σx and σy are determined for the intersection that matches the
y axis.
Basedon thepreviouslydefinedequations, analgorithm fornumerical solutionof theproblem

is created.The fact that isworthmentioning is that theprogramming is donewith complex sizes.
Based on such obtained numerical values of stresses, appropriate diagrams of stresses change
at the outer contour of the disc, and at the contour of the eccentric hole are created as given
in Section 5. Analysis of the data for the stress value σρ at the contour of the eccentric hole
obtained byusing theoreticalmodels has shown that the stresses are constant at any point of the
eccentric hole, which is logical because the first boundary condition was made on that contour
which implied that it was loaded with a constant pressure p = 1.0kN/cm2. Also, at the outer
contour of the disc, the values of stress σρ equal zero. The values of principal stresses σ1 and σ2
at any point of the disc with the eccentric circular hole are obtained when the values of stresses
σρ, σθ and τρθ are replaced in the following equations

σ1 =
σρ+σθ
2
+

√

(σρ−σθ
2

)2

+ τ2
ρθ

σ2 =
σρ+σθ
2
−

√

(σρ−σθ
2

)2

+ τ2
ρθ

(4.2)

At the outer contour of the disc, the values of stress σ2 equal zero at any point. The stresses
values σx, σy and τxy at any point of the disc with the eccentric hole can be obtained when σρ,
σθ and τρθ are replaced in the following equations

σx =
σρ+σθ
2
+
σρ−σθ
2
cos2θ− τρθ sin2θ

σy =
σρ+σθ
2
−
σρ−σθ
2
cos2θ+τρθ sin2θ

τxy =
σρ−σθ
2
sin2θ+ τρθ cos2θ

(4.3)

5. Determination of stress using FEM

TheFEM is based on physical discretization of the considered continuumwith elements of finite
dimensions and simple shape. In the study, a numerical model was created. It was a steel disc
weakened byan eccentric circular holewhosedimensions and load are identical to thedimensions
and load used in the previous analysis by theCVM.Also, a homogeneous isotropic disc has been
considered, while thematerial of the disc is steel withmodulus of elasticity E =21000kN/cm2,
and Poisson’s ratio ν = 0.33. The calculation was carried out by using ANSYS 12 software
package and the finite elements such as thin plates were applied. The FEM model consists
of 10492 nodes and 1816 finite elements. It is important to note that the input data for the
spatial discretization and mesh generation were not previously adjusted, but a mesh that is
generated automatically by the programANSYS12was used.The discwas loadedwith internal



1104 M. Bižić et al.

Fig. 4. Equivalent stress

Fig. 5. Normal stress σx (a) and σy (b)

Fig. 6. Shear stress τxy



Mathematical modelling of a disc weakened... 1105

pressure equal to p =1.0kN/cm2, acting within the internal contours of the eccentric hole. As
a consequence of the given discretization and load, the corresponding stresses in the disc were
obtained as shown in Figs. 4-6.

In order to compare the results, the numerical stresses obtained by the FEM were read in
specific points of the disc. Locations of these points at the outer contour of the disc and at the
contour of the eccentric circular holewere the same as 36 points thatwere definedby the angle θ
ranging from 0◦-360◦, with a step of 10◦, defined in Section 3. Locations of points along the
y axis for reading the stresses values were also the same as those in Section 3.

The principal stresses σ1 and σ2 were determined according to equation (2.5), and also the
diagrams of principal stresses at the outer contour of the disc and at the contour of the eccentric
holewere formed.Also, the diagrams of the stresses σx and σy for the intersection thatmatches
the y axis were formed.

6. Comparison of results obtained by CVM and FEM

The values of the principal stress σ1 at the outer contour of the disc (Fig. 7a) coincide only at
the maximum stress values for the angle θ = 180◦ (y = 0.5cm). At other points of the outer
contour, the stresses differ by a greater extent. The trend of stress distribution in both cases is
approximately the same. By the FEM, some very small and negligible values of the principal
stress σ2 at the outer contour of the disc (Fig. 7b) are obtained, while according to theoretically
obtained equations (4.2) derived by theCVM, such stress is equal to zero. Themaximumvalues
of the principal stress σ1 at the contour of the eccentric hole (Fig. 8a) calculated by the CVM
are obtained for the angles θ = 100◦ and θ = 260◦, while by using the FEM are obtained for
the angles θ =110◦ and θ =240◦. The deviations are very small and, in this case, the trend of
stress distribution is also approximately the same. The values of the principal stress σ2 at the
contour of the eccentric hole (Fig. 8b) calculated by the CVMare constant σ2 =−1.0kN/cm

2,
which is logical because there is constant pressure of the same intensity inside the contour of
the eccentric hole. In comparisonwith the FEM, these values are somewhat different, where the
deviations range from −5.9% to 2.1%. As for the values of stress components σx and σy at
the intersection which coincides with the y axis, the comparative analysis has shown that the
values largely overlap, and that overlapping ranges below 2% (Figs. 9a and 9b). Therefore, the
stresses obtained by the CVM and FEM are very similar in values and trend of distribution.

Fig. 7. Comparative diagram of the stress σ1 (a) and σ2 (b) at the outer contour of the disc



1106 M. Bižić et al.

Fig. 8. Comparative diagram of the stress σ1 (a) and σ2 (b) at the contour of the eccentric hole

Fig. 9. Comparative diagram of the stress σx (a) and σy (b) at the intersection of the y axis

7. Conclusion

The task of this paper is stress analysis in a homogeneous isotropic disc weakened with an ec-
centric circular holewhich is loaded by pressure in the internal contour of the hole. The complex
variable method (CVM)was applied which is based on the application of Muskhelishvili’s com-
plex variable function technique. The formed mathematical model allows complete analytical
solution of the stress state of the disc, especially the contour of the holewhere the stress concen-
tration is present. The technique was applied to the specific example of the disc, for which the
concrete numerical values of stresses were determined. Verification of the obtained results was
carried out by the finite element method (FEM) using the software package ANSYS 12. Com-
parative analysis has shown that stresses obtained by the CVM and FEM are very similar in
values and trend of distribution,which confirms the correctness of the establishedmathematical
model. The application of the results of this paper is of great importance for quality design and
optimization of thin-walled structures of disc type weakened by a circular hole.

References

1. Arshadnejad S., Goshtasbi K., Aghazadeh J., 2009, Stress concentration analysis between
two neighbouring circular holes under internal pressure of a non-explosive expansion material,
Yerbilimleri, 30, 3, 259-270



Mathematical modelling of a disc weakened... 1107

2. Bakhshandeh K., Rajabi I., Rahimi F., 2008, Investigation of stress concentration factor for
finite-width orthotropic rectangular plates with a circular opening using three-dimensional finite
element model, Strojniski Vestnik – Jour. of Mech. Eng., 54, 2, 140-147

3. Batra R.C., Nie G.J., 2010, Analytical solutions for functionally graded incompressible ec-
centric and non-axisymmetrically loaded circular cylinders, Composite Structures, 92, 1229-1245,
doi:10.1016/j.compstruct.2009.10.022

4. Bizic M., Petrovic D., 2011, Identification of the effects of circular hole on the stress state of
homogeneous isotropic uniaxial tensioned plate, IMK-14 – Istraživanje i Razvoj, 17, 3, 17-22

5. Bojic N., Jugovic Z., Popovic M., 2010, Impact form open when uniaxsial straining, IMK-14
– Istrazivanje i Razvoj, 16, 2, 17-22

6. Bower F.A., 2010,Applied Mechanics of Solids, CRCPress, Boca Raton

7. Chandrashekhara K., Muthanna S.K., 1978, Stress analysis of a thick plate having
a circular hole under axisymmetric radial load, Archive of Applied Mechanics, 47, 1, 1-9,
doi:10.1007/BF00536613

8. Chen K.T., Ting K., Yang W.S., 2000, Stress analysis of two-dimensional perforated pla-
tes using boundary element alternating method, Computers and Structures, 75, 5, 515-527,
doi:10.1016/S0045-7949(99)00103-0

9. Chen P.S., Archer R.R., 1989, Stress concentration factors due to the bending of a thick plate
with circular hole,Archive of Applied Mechanics, 59, 6, 401-411, doi:10.1007/BF00534307

10. Huan-Chun S., Ju-YongZ., Guo-Mei Y., Ji-Zhu Z., 1987,Complex variable functionmethod
for hole shape optimization in an elastic plane, Journal of Applied Mathematics andMechanics, 8,
2, 137-146, doi:10.1007/BF02019086

11. Kratochvil J., Becker W., 2011, Asymptotic analysis of stresses in an isotropic linear ela-
stic plane or half-plane weakened by a finite number of holes, Archive of Applied Mechanics,
doi:10.1007/s00419-011-0587-z

12. Lu J.-K., 1995, Complex Variable Methods in Plane Elasticity, World Scientific Publishing, Sin-
gapore

13. Mitrinovic D.S., 1981,Complex Analysis, Gradjevinska Knjiga, Belgrade

14. Mizushima I., Hamada M., 1983, Stress analysis around circular hole in infinite plate with rigid
disk, JSME, 26, 218, 1296-1301

15. Muskhelishvili N.I., 1963, Some Basic Problems of the Mathematical Theory of Elasticity, P.
Noordhoff, Groningen

16. Radi E., Strozzi A., 2009, Jeffery solution for an elastic disk containing a sliding eccentric
circular inclusion assembled by interference fit, International Journal of Solids and Structures, 46,
4515-4526, doi:10.1016/j.ijsolstr.2009.09.016

17. Savin G.N., 1961, Stress Concentration Around Holes, PergamonPress, London

18. SimhaK.R.Y.,MohapatraS.S., 1998,Stress concentrationaround irregularholesusing complex
variable metod, Sadhana, 23, 4, 393-412

19. TimoshenkoS.,GoodierJ.N., 1951,Theory of Elasticity, 2nded.,McGraw-HillBookCompany,
NewYork

20. Timoshenko S.,Woinowsky-Kreiger S., 1959,Theory of Plates and Shells, 2nd ed.,McGraw-
Hill Book Company, NewYork

21. Troyani N., Gomes C., Sterlacci G., 2002, Theoretical stress concentration factors for short
rectangular plates with centred circular holes, Journal of Mechanical Design, 124, 1, 126-128,
doi:10.1115/1.1412849

22. Wang Q.Z., 2004, Stress concentration factors for an eccentric circular hole in a finite-width strip
or in a semi-infinite plate in tension, The Journal of Strain Analysis for Engineering Design, 39,
6, 625-630, doi:10.1243/0309324042379284



1108 M. Bižić et al.

23. Wu L., Markenscoff X., 1996, Singular stress amplification between two holes in tension,
Journal of Elasticity, 44, 2, 131-144, doi:10.1007/BF00042474

24. Yang Q., Gao C.F., Chen W., 2010, Stress analysis of a functional gradedmaterial plate with
a circular hole,Archive of Applied Mechanics, 80, 8, 895-907, doi:10.1007/s00419-009-0349-3

25. Zhang Y., Shen S., 2011, Influences of surface on the interaction between holes or edge,Archive
of Applied Mechanics, 81, 3, 385-395, doi: 10.1007/s00419-010-0416-9

Modelowanie matematyczne dysku osłabionego mimośrodowo umiejscowionym

otworem kołowym

Streszczenie

Zadaniempodjętymwartykule jest identyfikacja stanunaprężeńw jednorodnym, izotropowymdysku
osłabionymmimośrodowoumiejscowionymwycięciemkołowym,obciążonymwzdłużbrzeguwewnętrznym
naciskiem. Stosując metodę zmiennej zespolonej, w pracy sformułowano matematyczny model układu,
który pozwolił na uzyskanie w pełni analitycznego opisu rozkładu naprężeńw dysku. Zastosowanameto-
dologia, zdaniem autorów, może być użyta dla dowolnego dysku z mimośrodowym otworem o kształcie
koła. W pracy przeprowadzono ponadto analizę porównawczą z wynikami uzyskanymi numerycznie za
pomocą metody elementów skończonych z wykorzystaniem pakietu ANSYS 12. Otrzymane rezultaty
badań mogą mieć duże znaczenie praktyczne z punktu widzenia jakości projektowania i optymalizacji
cienkościennych konstrukcji zawierających elementy strukturalne z kołowymwycięciem.

Manuscript received February 24, 2012; accepted for print March 23, 2012