Stress state of a polygonal plate having a central circular hole with two linear cracks in the small physical non-linear condition

Eyvazov, Etimad Bayramoglu; Öge, Mecit

doi:https://doi.org/10.5194/ms-10-553-2019

Articles | Volume 10, issue 2

https://doi.org/10.5194/ms-10-553-2019

Articles | Volume 10, issue 2

Research article

18 Nov 2019

Research article |

| 18 Nov 2019

Stress state of a polygonal plate having a central circular hole with two linear cracks in the small physical non-linear condition

Etimad Bayramoglu Eyvazov and Mecit Öge

Abstract

In this study, the stress state of a polygonal plate having a central circular hole with two linear cracks was considered. External forces were applied to the contours of polygonal plate and internal forces were applied to the contours of the circular hole. Conformal mapping function, initially developed by Kuliyev, was used to the make stress analysis of the polygonal plate in physical non-linear statement. As a result, stress concentrations were determined by using conformal mapping function at the end points of cracks. Critical loads where fracture began were defined for different loads. Analytical results of solutions were compared and found to be in agreement with the numerical results in the literature.

This paper has been retracted.

Download & links

Article (PDF, 482 KB)

Retraction notice
This paper has been retracted. Please read the editorial note.
Article (482 KB)

Editorial note

Download & links

This paper has been retracted. Please read the editorial note.

How to cite

How to cite.

Dates

Received: 16 May 2019 – Accepted: 07 Oct 2019 – Published: 18 Nov 2019

1 Introduction

In recent years, solution of physical non-linear problems in the fields of physics and engineering that cannot be solved through linearization of assumptions and principal equations (such as stress-strain relationship), has necessitated making improvements in classical elasticity theory, which in turn has led to the development of non-linear elasticity theory.

All assumptions and resulting errors mentioned in elasticity theory up to date have shown that, there is insufficient information on the new variants of non-linear elasticity theory. In classic elasticity theory, linearization is mainly assumed to be in two directions. The first one is the geometry of a deformable body and the other is physical properties of material (relationship of stress-strain, i.e. the Hooke's law).

Since classical elasticity theory cannot be applied for non-linear materials, Muskhelishvili (1962) developed a method to solve the stress states of non-linear materials, in which complex variable functions are used to make stress analysis.

Stress concentration around a hole in an infinite plate subjected to a uniform load was investigated by Batista (2011) using Muskhelishvili's (1962) conformal mapping complex variable method. Pan et al. (2013) performed the stress analysis of a finite plate having a rectangular hole subjected to uniaxial tension using a modified form of Muskhelishvili's (1962) conformal mapping function. Sharma (2012) examined the stress distribution around polygonal holes in an infinite plate subjected to arbitrary biaxial in-plane loadings at infinity also by using Muskhelishvili's (1962) complex variable method. In his study, plane stress finite models were prepared in ANSYS and the results were compared with those of the complex variable method (Sharma, 2012). Sharma (2015) obtained a general solution in hypocycloidal hole with cusps in infinite anisotropic plate in order to calculate the stress intensity factors by using Muskhelishvili's (1962) complex variable method. Kuliyev (2003), on the other hand, investigated the stress-strain state of a polygonal plate weakened by a central elliptic hole with two linear cracks by using modified conformal mapping function. Kuliyev (2010) also obtained the solution for critical angular velocity in polygonal plate weakened by elliptic holes using a modified version of conformal mapping function.

In this paper, stress components of a polygonal plate having central circular hole with two linear cracks were obtained. Conformal mapping function which was developed by the author (Kuliyev, 1991, 2017) was used to determine the stress intensity factors for complex geometries. In this work, a complex geometry with additional two linear cracks is solved for the first time using Kuliyev's mapping function, and the results were compared with the literature findings which proved to be in good agreement.

2 Problem statement and solution methods

Physical law of deformability, i.e., relationship of stress-strain components is explained in two ways:

General functional relationship of stress and strain (between two different power symmetric tensor components) (Kauderer, 1961).
Inclusion of deformation energy function into elasticity potential (Globenko, 1974).

As mentioned above, determination of real physical relationships between stress-strain components and solution of some complex problems in engineering constitute the essence of non-linear elasticity theory.

As reported in the literature (Kuliyev, 1991, 2001; Amenzade, 1976; Berezhnitskii et al., 1919; Sherman, 1947, 1951) Airy Stress Function F⁽ⁿ⁾ is determined by two analytic functions, namely φ(z) and ψ(z) in linear elasticity theory (plane problems)

\begin{matrix} (1) & F^{(0)} = Re [\overline{z} \cdot φ (z) + \int ψ (z) d z] \end{matrix}

where, Re shows the real part of the expression in bracket.

The functions of elasticity theory named as φ(z) and ψ(z) in finite or infinite form of bodies (plate, cylinder, prism and etc.) and connected domains are taken as follows (Kuliyev, 1991, 2001; Amenzade, 1976; Berezhnitskii et al., 1919; Sherman, 1947, 1951).

2.1 For finite doubly – connected domains

The regular functions φ(z) and ψ(z) are as follows (Kuliyev, 1991, 2001; Amenzade, 1976; Sherman, 1947, 1951).

\begin{matrix} (2) & \begin{aligned} φ (z) & = \sum_{k = 1}^{\infty} α_{k} \cdot ξ_{1}^{- k} + \sum_{k = 0}^{\infty} A_{k} {(\frac{z}{A})}^{k} \\ ψ (z) & = \sum_{k = 1}^{\infty} β_{k} \cdot ξ_{1}^{- k} + \sum_{k = 0}^{\infty} B_{k} {(\frac{z}{A})}^{k} \end{aligned} \end{matrix}

where ξ₁ is the function of ξ=χ(z), inverse function (z=ω(ξ₁)) of which conforms L₁ inner contour around the unit circle in Fig. 1 (Kuliyev, 1991, 2001).

https://www.mech-sci.net/10/553/2019/ms-10-553-2019-f01

Figure 1A hexagonal cross-section under interior and exterior pressure.

Download

The coefficients A_k and B_k of analytic functions are as follows:

\begin{matrix} (3) & \begin{aligned} A_{k} & = \sum_{n = k}^{\infty}^{*} a_{n}^{(k)} \cdot a_{\frac{n - k}{q}}^{(n)} \\ B_{k} & = \sum_{n = k}^{\infty}^{*} b_{n} \cdot a_{\frac{n - k}{q}}^{(n)} \end{aligned} \end{matrix}

The coefficient $a_{n}^{(k)}$ in Eq. (3) is found for total contour and given in table form in (Kuliyev, 2017). The symbol ^∗ indicates the indices in series Eq. (3) with increment q. We assume that contours L₁ and L₂ are concentric in doubly connected domains.

Here, F⁽⁰⁾ is the Airy stress function in zero approximation. The function F⁽⁰⁾ is given by the known formula (the plate is under the influence of internal p₁ and external pressure p₂ equi-distributed hydrostatic pressure in Fig. 1) for polygonal doubly connected plates (Kuliyev, 1991, 2001; Muskhelishvili, 1962; Amenzade, 1976; Berezhnitskii et al., 1919).

In the considered case, the boundary conditions are as follows in zero approximation:

on the inner contour L₁ (with r radius circle and two linear cracks)
$\begin{matrix} (4) & σ_{ρ}^{(0)} |_{ρ = t_{1}} = - P_{1} \end{matrix}$
on the outer contour L₂ (equilateral polygon)
$\begin{matrix} (5) & σ_{ρ}^{(0)} |_{ρ = t_{2}} = - P_{2} \end{matrix}$

For the next approximations (the first and the second), the boundary conditions are as follows:

\begin{matrix} (6) & σ_{ρ}^{(1)} |_{ρ = t_{1}} = 0, σ_{ρ}^{(1)} |_{ρ = t_{2}} = 0 \end{matrix}

In Eqs. (4), (5), (6), t₁ and t₂ are the affixes of the points of the contours $L_{j} (j = 1; 2)$ , respectively. The function F⁽¹⁾ is found from the solution of Eq. (1).

\begin{matrix} (7) & Δ Δ F^{(2)} + L_{2} [F^{(0)}, F^{(1)}] = 0 \end{matrix}

The solution of Eq. (7) is found from the sum of general solution ΔΔF⁽¹⁾=0 of the first differential equation $(F_{gen .}^{(1)})$ and the special solution $(F_{spec .}^{(1)})$ of the mentioned differential equation.

\begin{matrix} (8) & F^{(1)} = F_{gen .}^{(1)} + F_{spec .}^{(1)} \end{matrix}

It is known that the general solution of the homogeneous differential equation ΔΔF⁽¹⁾=0 is taken in the following way (in the general case):

\begin{matrix} (9) & F_{gen .}^{(1)} = c_{1} \cdot \ln ρ + c_{2} \cdot ρ^{2} \end{matrix}

Note that, if we pass to polar (ρ,θ) coordinates, we get the following differential equation for the function F⁽¹⁾ (Kuliyev, 2017).

\begin{matrix} (10) & Δ Δ F^{(1)} = \sum_{k = 1}^{\infty} ρ^{- k} \cdot V_{1} (k) + \sum_{k = 0}^{\infty} ρ^{k} \cdot V_{2} (k) \end{matrix}

In this equation, the quantities V₁(k) and V₂(k) are algebraic sums of the same power terms of the variable p (for the fixed θ angle). The special solution $F_{spec .}^{(1)}$ is found from fourfold integration of Eq. (10) in the following way (Kuliyev, 1991, 2001):

\begin{matrix} (11) & \begin{aligned} F_{spec .}^{(1)} & = \sum_{k = 0}^{\infty} ρ^{- k + 4} \cdot α \cdot V_{1} (k) \\ + \sum_{k = 0}^{\infty} β \cdot ρ^{k + 4} \cdot V_{2} (k) \end{aligned} \end{matrix}

where

\begin{matrix} (12) & \begin{aligned} α & = {[{(4 - k)}^{2} \cdot {(2 - k)}^{2} + α_{1} + α_{2}]}^{- 1} \\ β & = {[{(k + 4)}^{2} \cdot {(2 + k)}^{2} + β_{1} + β_{2}]}^{- 1} \end{aligned} \end{matrix}

The unknown coefficients c₁ and c₂ in Eq. (9) are found from the boundary conditions of the considered problem (stress-strain state of a doubly-connected polygonal plate subjected to equi-distributed internal pressure p₁ and external pressure p₂) (Kuliyev, 1991, 2001).

If on the L₁ inner contour

ρ = t_{1} then σ_{ρ}^{(1)} = 0

If on the L₂ external contour

\begin{matrix} (13) & ρ = t_{2} = A then σ_{ρ}^{(1)} = 0 \end{matrix}

After finding the expression $F (ρ, θ, λ t) = F^{(0)} (ρ, θ) + λ F^{(1)} (ρ, θ) + λ^{2} \cdot F^{(2)} (ρ, θ) + \dots$ for the stress components σ_ρ and σ_θ we get the following expressions in the form of the following series:

\begin{matrix} (14) & \begin{aligned} σ_{ρ} & = σ_{ρ}^{(0)} + λ \cdot σ_{ρ}^{(1)} + λ^{2} \cdot σ_{ρ}^{(2)} + \dots \\ σ_{θ} & = σ_{θ}^{(0)} + λ \cdot σ_{θ}^{(1)} + λ^{2} \cdot σ_{θ}^{(2)} + \dots \\ τ_{ρ θ} & = τ_{ρ θ}^{(0)} + λ \cdot τ_{ρ θ}^{(1)} + λ^{2} \cdot τ_{ρ θ}^{(2)} + \dots \end{aligned} \end{matrix}

The obtained general solution is demonstrated by the numerical examples.

2.2 Numerical Results

2.2.1 Stress state of a cylinder made of small physical nonlinear material and possessing an annular hole and two linear cracks

In linear elasticity theory, this problem was solved completely (within Hooke's law). Therefore, here we give only the results obtained in Kuliyev (1991, 2001).

https://www.mech-sci.net/10/553/2019/ms-10-553-2019-f02

Figure 2Hollow round cylinder under interior and exterior pressure.

Download

The ratio of the main measures of the cross section is given in two variants (Fig. 2):

Variant I: $r / R = 0.5$ ; e=0.6 R; l=0.1 R
Variant II: $r / R = 0.7$ ; e=0.8 R; l=0.1 R.

The stresses $σ_{ρ}^{(0)}$ and $σ_{θ}^{(0)}$ found from characteristic points of the section (as tangential stresses τ_ρθ are too small, accepted τ_ρθ≈0) are given in Table 1. Stress intensities and critical stress values are given in Table 2.

At the most critical point of the section (at the end points of the cracks: $z = \pm e)$ , we get the following values of $10 K_{1} / p_{1} \sqrt{l}$ for the stress intensity factor K₁:

7.00 for variant I and
11.68 for variant II.

The following estimates are found for the critical value of the internal pressure p₁=p (i.e. fracture of the body-minimal stress for elongation of cracks):

For variant I: $P_{cr} = 1, 46 [σ_{b}]$
For variant II: $P_{cr} = 0, 883 [σ_{b}]$ (Table 2)

where technical ultimate stress [σ_b] (maximum stress of the material) is known from references (Kuliyev, 1991, 2001).

Table 1Stress coefficient concentrations.

Download Print Version | Download XLSX

Table 2Stress intensity (by $10 K_{1} / P \sqrt{l}$ ) and critical value of stress.

Download Print Version | Download XLSX

After finding the functions φ(z) and ψ(z), within the linear elasticity theory (in zero approximation), the stress function $F^{(0)} (x, y) = F^{(0)} (ρ, θ)$ is found by Eq. (1). Then F⁽¹⁾(ρ,θ) is found from Eq. (10) of the stress function in the first approximation. According to the found expression of the stress function F⁽¹⁾(ρ,θ), the stress σ_θ found from the characteristic points of the section by means of Eq. (14) is given in and Fig. 3.

https://www.mech-sci.net/10/553/2019/ms-10-553-2019-f03

Figure 3Hollow cylinder under interior pressure. The dotted line indicates a linear problem, the solid line a non-linear problem. (1: pure copper; 2: open heath steel).

Download

Table 3Stress concentrain of a hexogonal plate.

Download Print Version | Download XLSX

Table 4Stress intensity coefficient ( $10 K_{1} / ¶ \sqrt{l}$ ) and critical value of load (P_cr) (hexagonal cylinder).

Download Print Version | Download XLSX

2.2.2 Stress state of a hexagonal plate with annular hole and two linear cracks in physical non-linear statement

The solution of this problem in linear elasticity theory is known from references (Kuliyev, 1991, 2001). The stresses $(σ_{θ}^{(0)})$ found at any point of a hexagonal plate are calculated for different variants (for different values of the ratio of cross section measures) and given in the table. The stress intensity factor (K₁) and critical value P_cr of internal pressure are found and given for each variant at the end of Tables 3 and 4 according to the stresses found at the end points of the linear cracks. The known expression of the function is taken into account in differential Eq. (10) obtained for the solution of physical non-linear problem and accordingly, the stress for different values of internal pressure P_cr at the most critical points $F^{(0)} (x, y) = F^{(0)} (ρ, θ)$ of the section are found and given in table form.

The following expressions were considered for the stresses at the tips of the cracks.

Variant I: $r / A = 0.5$ ; $e / A = 0.6$ ; m=0.04; $b = 1, 922$ ;
$σ_{θ} / p = σ_{θ / p}^{(0)} + α σ_{θ}^{(1)} / p = 1.68 - 9.2 \cdot α \cdot p^{2}$
Variant II: $r / A = 0.7$ ; $e / A = 0.8$ ; m=0.04; b=1.37r
$σ_{θ} / p = σ_{θ / p}^{(0)} + α σ_{θ}^{(1)} / p = 3.21 - 14.62 \cdot α \cdot p^{2}$
where
$α = λ = \frac{g_{2} k}{(3 K + G) G^{2}} = 0.225 \times 10^{- 6} {cm}^{4} {kg}^{- 2}$
for pure copper.
$α = λ = 0.032 \times 10^{- 6} {cm}^{4} {kg}^{- 2}$
for open-hearth steel.

As indicated by the numerical results of both problems in physical non-linear statement, the maximal values of stresses are smaller than that of the linear statement (peak points of stress diagrams are smoothened) (Fig. 1).

2.2.3 Stress state of a thick-walled pipe (Lame problem) in small physical non-linear statement

In classic elasticity theory, the following formula is obtained for the stress function F⁽⁰⁾ (Kuliyev, 1991, 2001).

\begin{matrix} (15) & F^{(0)} (ρ, θ) = β_{1} \cdot \frac{ρ^{2}}{r^{2}} + d_{1} \cdot r \cdot \ln ρ \end{matrix}

When the pipe is under the internal pressure p₁ and external hydrostatic pressure p₂, the coefficients β₁ and d₁ are determined by the following expressions (Kuliyev, 1991, 2001).

\begin{matrix} (16) & β_{1} = \frac{r}{2} \cdot \frac{p_{1} \cdot r^{2} - p_{2} \cdot R^{2}}{R^{2} - r^{2}}; d_{1} = \frac{r^{2} R^{2} \cdot (p_{2} - p_{1})}{(R^{2} - r^{2}) \cdot r} \end{matrix}

if the Eq. (16) is written in Eq. (15);

\begin{matrix} (17) & \begin{aligned} F^{(0)} (ρ, θ) & = \frac{r^{2} R^{2} \cdot (p_{2} - p_{1})}{R^{2} - r^{2}} \cdot \ln ρ \\ + \frac{1}{2} \frac{p_{1} \cdot r^{2} - p_{2} \cdot R^{2}}{R^{2} - r^{2}} ρ^{2} \end{aligned} \end{matrix}

The obtained expression completely coincides with the expressions obtained by Tsurpal (1962). The measures of thick-walled pipes (the ratio of radius R of external circle with radius r of internal circle) were given in two variants (pure copper, aluminum, bronze and open–heath steel were chosen as pipe's materials. The following expressions are obtained for stress concentration on the inner circle for the considered problem (Lame Problem).

Variant I: $σ_{θ} / p = 1.666 - 4.444 \cdot λ \cdot p^{2}$ ; p₁=p
Variant II: $σ_{θ} / p = 4.55 - 25.3 \cdot λ \cdot p^{2}$ ; p₁=p

The stresses “σ_θ” were calculated for different values of internal pressure (p₁) ( $20 \leq p \leq 500$ ) by means of the formulas given in graphic form in both variants (Fig. 4). In the graph, the straight line represents the problem in non–linear statement and the dotted line represents the problem in linear statement.

As seen in the tables, a decline at a rate of 0.1 % occurs for the values lower than 100 MPa, resulting in a negligible difference between linear and non-linear states. As for the values higher than or equal to 100 MPa, the stress difference between linear and nonlinear states is observed to be significantly high.

https://www.mech-sci.net/10/553/2019/ms-10-553-2019-f04

Figure 4The Lame problem statement.

Download

3 Conclusion

In this work, stress concentrations were determined by using conformal mapping function at the end points of cracks. Critical loads where fracture began were defined for different loads. Analytical results of solutions were compared and found to be in agreement with the numerical results in the literature. The obtained numerical calculations show that, as the nonlinearity of the materials increase, the stress distributions at the peak points of the stresses become more regular and physical–nonlinearity of the material smoothens the peaks of stresses by decreasing maximum value of stresses and moderately increasing minimum values. Also, the physical nonlinearity of the material decreases the stress concentration, thus decreasing the maximum stress at the inner points of the cross-section.

Data availability

The underlying research data can be provided upon reasonable request by Etimad Bayramoglu Eyvazov (etimadeyvazov@karabuk.edu.tr).

Appendix A: Nomenclature

K	bulk modulus
k(s₀)	mean stress of s₀ function
$g (t_{0}^{2})$	intensity function
G	shear modulus
s₀	expression for the reduced mean stress
E	Young modulus
F	Airy stress function
I₁	strain tensor
L⁽ⁿ⁾	linear operator
λ_ij	the coefficients which determine the physical properties of the considered materials
τ₀	expression for intensity of the reduced tangential stresses
υ	Poisson ratio
χ	elongation function
Υ	shear strain function
ε_ij	strain components
δ_ij	the coefficients which determine the physical properties of the considered materials
σ_ij	stress components
Δ	Laplace operator
φ(z), ψ(z)	analitic functions
q	symmetry axis number

Author contributions

EBE conducted the numerical analyses and wrote the majority of the paper. MÖ supervised the findings, organized and structured the paper.

Competing interests

The authors declare that they have no conflict of interest.

Acknowledgements

The authors gratefully acknowledge the valuable support of Sabir A. Kuliyev, who recently passed away.

Review statement

This paper was edited by Doina Pisla and reviewed by Farshid Khosravi Maleki and one anonymous referee.

References

Amenzade, Y. A.: Theory of Elasticity, Vysshaya Shkola, Moscow, 1976 (in Russian).

Batista, M.: On the stress concentration around a hole in an infinite plate subject to a uniform load at infinity, Int. J. Mech. Sci., 53, 254–261, https://doi.org/10.1016/j.ijmecsci.2011.01.006, 2011.

Berezhnitskii, L. T., Delyavskii, M. V., and Panasyuk, V. V.: Bending of thin plates with crack type defects, Naukova Dumka, Kiew, 1919 (in Russian).

Globenko, T. M.: Influence of physical nonlinearity of material on stress concentration near square hole, Calc. Meth. Civ. Eng., 1, 39–55, 1974.

Kauderer, G.: Nonlinear mechanics, Inostrannaya Literatura, Moscow, 1961.

Kuliyev, S.: Stress state of compound polygonal plate, Mech. Res. Com., 30, 519–530, 2003.

Kuliyev, S.: Uniform rotation of a polygonal plate weakened by two linear crack holes, Mech. Res. Com., 37, 184–190, 2010.

Kuliyev, S. A.: Two-dimensional Problems of Theory of Elasticity, M. Stroyizdat, Moscow, 1991.

Kuliyev, S. A.: Some problems of elasticity theory, Azerneshr, Baku, 2001.

Kuliyev, S. A.: Conformally-mapping Functions of Complex Domains, Lap Lambert Academic Publishing, Azerbaijan, Baku, 2017.

Muskhelishvili, N. I.: Some basic problems of mathematical theory of elasticity 2nd English edn., P. Noordhooff Ltd., the Netherlands, 1962.

Pan, Z., Cheng, Y., and Liu, J.: Stress analysis of a finite platewith a rectangular hole subjected to uniaxial tension using modified stress functions, Int. J. Mech. Sci, 75, 265–277, https://doi.org/10.1016/j.ijmecsci.2013.06.014, 2013.

Sharma, D. S.: Stress distribution around polygonal holes, Int. J. Mech. Sci., 65, 115–124, https://doi.org/10.1016/j.ijmecsci.2012.09.009, 2012.

Sharma, D. S.: Stress intensity factors for hypocycloidal hole with cusps in infinite anisotropic plate, Theo. Appl. Fract. Mech., 75, 44–52, https://doi.org/10.1016/j.tafmec.2014.10.005, 2015.

Sherman, D. I.: On a solution method of elasticity theory problems for doubly-connected domains, Dokl. AN SSSR, 55, 704–710, 1947.

Sherman, D. I.: On stresses in plane weight medium with two same, symmetrically arranged annular holes, PMM, 4, 751–761, 1951.

Tsurpal, I. A.: Experimental determination of elastic constants in the nonlinear theory of elasticity, Prikl. Mekh., 8, 55–63, 1962.

Articles

Download

This paper has been retracted. Please read the editorial note.

Article (482 KB)
Full-text XML

Short summary

The stress state of a polygonal plate having a central circular hole with two linear cracks was considered. External forces were applied to the contours of polygonal plate and internal forces were applied to the contours of the circular hole. Stress concentrations were determined at the end points of cracks using conformal mapping function. Critical loads were defined for different loads. Analytical results were found to agree with the numerical literature results.