﻿ An Infinite Elastic Plate Weakened by a Generalized Curvilinear Hole and Goursat Functions

Applied Mathematics
Vol.5 No.4(2014), Article ID:43794,16 pages DOI:10.4236/am.2014.54070

An Infinite Elastic Plate Weakened by a Generalized Curvilinear Hole and Goursat Functions

Mohamed Abdellah Abdou1, Azhar Rashad Jan2

1Department of Mathematics, Alexandria University, Alexandria, Egypt

2Department of Mathematics, Umm Al-Qura University, Mecca, KSA

Email: abdella_777@yahoo.com, azhaarjaan@ymail.com

Received 23 December 2013; revised 23 January 2014; accepted 2 February 2014

ABSTRACT

Complex variables method has been used to solve the first and second fundamental problems for an infinite plate weakened by a generalized curvilinear hole C. The curvilinear hole is conformally mapped on the domain outside or inside a unit circle γ using a general rational mapping function with complex constants. Many special and new cases are derived from this work. Some of the work of the previous authors in this domain will be considered as special cases of this paper. Also the interesting cases when the shape of the hole takes different famous shapes are included. The components of stresses for some examples are obtained.

Keywords:Goursat Functions; Conformal Mapping; Curvilinear Hole; Stress Components

1. Introduction

The boundary value problems for isotropic homogeneous performed infinite plates have been discussed by several authors: see Colton and Kress [1] , Popov [2] , Noda et al. [3] and Schinzinger and Laura [4] . Some authors used Laurent’s theorem to express the solution in the series form, see England [5] , Parkus [6] and Kalandiya [7] . Others used complex variables method of Cauchy integrals to express the solution of the boundary value problems in the form of two complex potential functions, Goursat functions, by using many rational mappings, see Muskhelishvili [8] , El-Sirafy and Abdou [9] , Abdou and Khar-Eldin [10] , Abdou and Khamis [11] , Abdou [12] and Abdou et al. [13] . In all previous works, the coefficients of the rational mappings were real.

It is worth mentioning that Exadaktylos and Stavropoulou [14] and Exadaktylos et al. [15] considered rational mapping functions with complex constants that conformally maps the holes inside a unit circle, using Laurent’s method. Also Abdou and Asseri [16] [17] considered more general rational mapping functions with complex constants that conformally maps the holes outside and inside a unit circle, using Cauchy singular method. All the previous four works will be considered as special cases of this work.

It is known that, see Muskhelishvili [8] , the first and second fundamental problems in the plane theory of elasticity are equivalent to finding two analytic functions and of one complex argument. These analytic potential functions, Goursat functions, must satisfy the boundary conditions

, (1)

where, and is a given function of stress, for the first fundamental problem. While

, and is a given function of the displacement for the second fundamental problem; l, m are called the Lame’s constants and t denotes the affix of a point on the boundary.

In terms of the rational mapping function does not vanish or become infinite for, the infinite region outside a closed contour conformally mapped outside the unit circle. The two complex potential functions, , in this case, take the forms

(2)

, (3)

where X, Y are the components of the resultant vector of all external forces acting on the boundary and, are complex constants. The two complex functions and are single valued analytic functions within the region outside the unit circle and. For the first fundamental problem, we have, and.

In the absence of body forces, Muskhelishvili [8] has considered the stress components in the plane theory of elasticity in the form

(4)

(5)

In this work, the complex variables method will be applied to solve the first and second fundamental problems for an infinite plate with a generalized curvilinear hole C conformally mapped on the domain outside a unit circle by the generalized rational mapping function

(6)

where; and is a parameter restricted such that does not vanish or become infinite outside the unit circle. The holes take different famous shapes which make these studies applicable for many phenomena throughout the nature like tunnels, caves, excavations in soil or rock, etc. Moreover, the results of Goursat functions when the transformation mapping (6) is conformally mapped inside the unit circle are discussed and obtained.

Also, many applications for the first and second fundamental problems are considered and the components of stress and strain have been obtained and plotted to investigate their physical meaning. Moreover, computer work using maple 9.5 has been used in applications to give the shapes of holes and curves of stresses with some calculations of stresses at their important points.

2. The Rational Mapping

The physical interest of the mapping (6) comes from its special cases and its different shapes of holes that can be obtained, see Figures 1-6.

From the rational mapping we can discuss the following:

1) The number of the holes corners is subjected to’s values. There are given by.

2) The shape of the hole depending on the values of n’s and m’s.

Figure 1. ℓ = 2, n = (0, 0), m = (1, 2), d = (2, 6).

Figure 2. ℓ = 2, n = (0.1, 0.1), m = (0.2, −0.2), d = (1, 1).

Figure 3. ℓ = 3, n = (0.001, 0.001), m = (0.1, −0.05), d = (1, 1).

Figure 4. ℓ = 5, n = (0.001, 0.001), m = (0.1, −0.05), d = (1, 1).

Figure 5. ℓ = 7, n = (0.001, 0.001), m = (0.1, −0.05), d = (1, 1).

Figure 6. ℓ = 10, n = (0.009, 0.007), m = (0.1, −0.05), d = (2, 2).

3) Entering none zero values of the complex constants m and d never gives symmetric graphs. While, entering zero values for all imaginary parts of both m and d, we get symmetric shapes around the x-axis. On the other hand, entering zero values for all real parts of both m and d, we get symmetric shape around the y-axis.

4) The complex constant m works on circling the shape from the symmetry situation and the circling angle is given by,. Positive values of q means that the circling will be in the positive direction i.e. in the anti clockwise direction and for negative values the circling will be in the negative direction i.e. in clockwise direction.

5) Using the rational mapping function, enables us to enter none zero vales of d, m and n complex constants and getting symmetric shapes around the x-axis. But by substituting zero values for real parts, we get the same shapes that have been gotten by using the rational mapping, and the invested shapes of it will be given by substituting zero values for the imaginary parts of the complex constants d, m and n .

6) The complex constant d works on expanding the corners of the hole shape.

3. Goursat Functions

In this section, we use the transformation mapping (6) in the boundary conditions (1), and complex variables method, Cauchy method, to obtain a closed form expression for the Goursat functions and respectively. Therefore, we write

, (7)

where

, (8)

(9)

and is a regular function for.

Using (7) in the boundary conditions (1) and on, we get

(10)

where

(11)

(12)

(13)

and

. (14)

The function with its derivatives must satisfy the Hölder condition. Multiplying both sides of (10) by and integrating with respect to s on, we have

(15)

where

(16)

and the complex constant b, will be determined, is given by

(17)

Differentiating (15) with respect to z, then using the result in (17), the complex constant b takes the form

, (18)

where

,

Also, the function can be determined from (1) in the form

(19)

where

The two formulas (15) and (19) are representing the Goursat functions for the first and second fundamental problems for an infinite elastic plate weakened by The two formulas (15) and (19) are representing the Goursat functions for the first and second fundamental problems for an infinite elastic plate weakened by a generalized curvilinear hole C, that can be transformed outside a unit circle g by the rational mapping (6).

An important new case for discussion is using the transformation mapping

(20)

This mapping function, when;, transforms the points in the z-plane inside the unit circle g in z-plane. And, in this case, the Goursat functions, become

(21)

(22)

4. Special Cases

Here, we discuss the following:

1) By considering the reality of the constants of the mapping (1.6), the Goursat functions, in this case, are agree with work of Abdou and Khar-Eldin [10] of Equations (15) and (19), on notation the difference in notation.

2) When and and for finite expansion, the transformation mapping (6), in this case, becomes

, d, nJ are complex constants                      (23)

The Goursat functions, in this case, become

(24)

(25)

The results of the two formulas (24) and (25) are in agreement with the work of Abdou and Asseri [16] , on noting the difference in notation.

3) When , the transformation mapping (6) becomes

(d, m, n are complex constants).                     (26)

The Goursat functions, in this case, of the two formulas (15) and (19) agree with the all results of Abdou and Asseri [17] .

4) In the mapping function (20) if we let m = 0, then for finite expansion, we will have the following mapping function

. (27)

with the corresponding Goursat functions

(28)

(29)

The three Formulas (27)-(29) are equivalent to those derived by Exadaktylos and Stavropoulou [14] , where they used Laurent's theorem, after considering in (27)-(29) the following special cases: and. The constant l is called the situ stress relief factor i.e. for l = 0 no excavation has been occurred and for l = 1, the tunnels is fully excavated. More information and physical meaning for l can be found in the work of Exadaktylos and Stavropoulou [14] .

5) Also, in (27), if we allow the index inside the summation sign to take the form, in addition to the consideration of the reality all constants. Then using in (28) and (29) the following:, , and; P is the intensity of the tensile stress and z is the above rational mapping, the results will agree with the work of Exadaktylos et al. [15] .

5. Applications

1) For and , we have the Goursat functions in the form

(30)

(31)

The complex constant b has been determined by Equation (18) and its value was calculated by using Maple 9.5. Here, we have the Goursat functions for an infinite plate weakened by a curvilinear hole C which is free from stresses. The plate stretched at infinity by the application of a uniform tensile stress of intensity P, making an angle q with the x-axis.

For, and the relation between the stress components sxx, syy, sxy and the angle q are considered in Figures 7-9.

2) For and f = Pt, P is a real constant, we have

, (32)

(33)

Thus, (32) and (33) give the solution of the first fundamental problem for an isotropic infinite plate with a curvilinear hole, when there are no external forces and the edge of the hole is subject to a uniform pressure P.

If in application (2) we write, we have the case of the first fundamental problem, when the edge of the hole is subject to uniform tangential stress T. The Goursat functions, in this case, are obtained directly from the two formulas (32) and (33) by putting –iT instead of P.

For and the relation between the stress components sxx ,syy ,sxy and the angle q, using Maple 9.5 are considered in Figures 10-12.

Figure 7. Max. and Min. values of σxx are [2.09232, [θ = 3.7365]], [−0.49508, [θ = 1.82742]].

Figure 8. Max. and Min. values of σyy are [1.24503, [θ = 1.84222]], [−2.63598, [θ = 0.73521]].

Figure 9. Maximum value of σxy is [0.25045, [θ = 2.75306]], Minimum value of σxy is [−7.61349, [θ = 3.65448]].

Figure 10. Max. and Min. values of σxx are [7.26497, [θ = 3.66684]], [−0.74520, [θ = 2.95389]].

Figure 11. Max. and Min. values of σyy are [0.45410, [θ = 1.87226]], [−6.94161, [θ = 3.66649]].

Figure 12. Maximum value of σxy is [1.99460, [θ = 3.77003]], Minimum value of σxy is [−0.63053, [θ = 0.63788]].

3) For and, we have

, (34)

, (35)

where

Here, we have the case of uni-directional tension of an infinite plate with a rigid curvilinear centre. The constant e, which represents the angle of rotation, can be determined from the condition that the resultant moment of the forces, acting on the curvilinear centre from the surrounding material, must vanish i.e.

(36)

Hence, we have

(37)

(38)

where

(39)

For θand m = 0.6 the relation between the stress components sxx, syy, sxy and the angle q are considered in Figures 13-15.

From the previous results, we can establish the following Case (1): In the case of Bi-axial tension, we have

Hence, we get

, (40)

, (41)

where

.

Figure 13. Max. and Min. values of σxx are [19.94724, [θ = 3.72527]], [−31.09023, [θ = 1.48899]].

Figure 14. Max. and Min. values of σyy are [43.21619, [θ = 1.52536]], [−28.86760, [θ = 3.68310]].

Figure 15. Maximum value of σxy is [27.14849, [θ = 4.28087]], Minimum value of σxy is [−18.54277, [θ = 3.15247]].

The complex constant b has been determined by Equation (18) and its value was calculated by using Maple 12. For n = 0.1 + 0.1i, m = 0.2 − 0.2i, d = 1 + i, P = 1/4 and c = 2. The relation θ between the stress components sxx, syy, sxy and the angle q are considered in Figures 16-18.

Case (2): When the curvilinear centre not allowed to rotate, i.e. when and. This means that the rigid curvilinear kernel in restrained in its original position by a couple which is not sufficient to rotate. The Goursat functions, in this case, become

, (42)

(43)

For , , P = 1/4 and, p »3.14159 the relation between the stress components sxx ,syy , sxy and the angle q , using Maple 12 are considered in Figures 19-21.

4) When the force acts on the centre of the curvilinear kernel and the stresses vanish at infinity. In this case the kernel can not be rotate and it remains in its original position. Hence, we get

(44)

, (45)

where

.

Therefore, we have the solution of the second fundamental problem in the case when (X, Y) acts on the centre of curvilinear hole.

For, n = 0.1 + 0.1i, m = 0.2 − 0.2i, d = 1 + i, , k = χ and X = Y = 10 the relation between the stress components sxx,syy, sxy and the angle q, using Maple 9.5 are considered in Figures 22-24.

6. Conclusion and Discussion

From the previous work the following discussion and results can be concluded 1) In the theory of two-dimensional linear elasticity one of the most useful techniques for the solution of the boundary value problem for a region weakened by a curvilinear hole is to transform the region into a simpler shape to get the solution directly without difficulties.

2) The transformation mapping, , c is a complex or real constant, transforms the domain of the infinite plate with a curvilinear hole onto the domain outside (when) or into the domain inside a unit circle.

3) The physical interest of the using mapping transform comes from its different shapes of holes it treats and different directions it takes. This mapping function deals with famous shapes of tunnels, thereon it is useful in studying the stresses around tunnels. In underground engineering the tunnel is assumed to be driven in a homogeneous, isotropic, linear elastic and pre-stressed geometrical situation. Also, the tunnel is considered to be deep enough such that the stress distribution before excavation is homogeneous. Excavating underground openings in soils and rocks is done for several purposes and in multi-sizes. At least, excavation of the opening will cause the soil or rock to deform elastically. The excavation in soil or rock is a complicated, dangerous and expensive process. The mechanics of this can be very complex. However, the use of conformal mapping that allows us to study stresses and strains around a unit circle makes it useful for engineers and easier for mathematicians.

Figure 16. Max. and Min. values of σxx are [0.19217, [θ = 2.88015]], [−0.21424, [θ = 4.21641]].

Figure 17. Max. and Min. values of σyy are [0.20395, [θ = 4.29250]], [−0.21998, [θ = 2.93249]].

Figure 18. Maximum value of σxy is [0.21042, [θ = 3.54667]], Minimum value of σxy is [−0.13798, [θ = 2.47646]].

Figure 19. Max. and Min. values of σxx are [0.42251, [θ = 3.40048]], [−0.32725, [θ = 2.60422]].

Figure 20. Max. and Min. values of σyy are [0.29017, [θ = 4.09478]], [−0.47644, [θ = 2.04499]].

Figure 21. Maximum value of σxy is [0.24900, [θ = 3.77559]], Minimum value of σxy is [−0.22325, [θ = 1.68785]].

Figure 22. Max. and Min. values of σxx are [0.06313, [θ = 6.28319]], [−0.02863, [θ = 2.37434]].

Figure 23. Max. and Min. values of σyy are [0.04161, [θ = 2.45943]], [−0.02730, [θ = 5.01912]].

Figure 24. Maximum value of σxy is [−0.00237, [θ = 3.58676]], Minimum value of σxy is [−0.04160, [θ = 1.06868]].

4) The complex variables method (Cauchy method) is considered one of the best methods for solving the integro differential equation, boundary value problem, of Equation (1) and obtaining the two complex potential functions, Goursat functions, and directly.

5) The stress is an internal force whereas positive values of it mean that stress is in the positive direction, i.e. stress acts as a tension force. On the other side, negative values of stress mean that the stress is in the negative direction, i.e. stress acts as a press force.

6) The most important issue deduced from mapping the stress components is that and vice verse.

7) When the perpendicular stress on y-axis is the maximum value and presents the body interior resistance of treatments (like rocks for example), whereas the perpendicular stress on x-axis is small according to y-axis. Thereon, it is better to treat the problem at points determined by angles that give minimum values of

.

References

1. Colton, D. and Kress, R. (1983) Integral Equation Methods in Scattering Theory. John Wiley, New York.
2. Ya. Popov, G. (1982) Contact Problems for a Linearly Deformable Functions. Kiev, Odessa.
3. Noda, N., Hentarski, R.B. and Tanigowa, Y. (2003) Thermal Stresses. Taylor and Francis, UK.
4. Schinzinger, R. and Laura, P.A. (2003) Conformal Mapping Methods and Applications. Dover Publications, New York.
5. England, A.H. (1971) Complex Variable Method in Elasticity. McGraw, London, New York.
6. Parkus, H. (1976) Thermo Elasticity. Spring-Verlag, Berlin,
7. Kalandiya, A.I. (1975) Mathematical Method of Two-Dimensional Elasticity. Mir Publishers, Moscow.
8. Muskhelishvili, N.I. (1953) Some Basic Problems of Mathematical Theory of Elasticity. Noordroof, Holland.
9. El-Sirafy, I.H. and Abdou, M.A. (1984) First and Second Fundamental Problems of Infinite Plate with a Curvilinear Hole. Journal of Mathematical and Physical Sciences, 18, 1-12.
10. Abdou, M.A. and Khar-Eldin, E.A. (1994) An Infinite Plate Weakened by a Hole Having Arbitrary Shape. Journal of Computational and Applied Mathematics, 56, 341-361. http://dx.doi.org/10.1016/0377-0427(94)90088-4
11. Abdou, M.A. and Khamis, A.K. (2000) On a Problem of an Infinite Plate with a Curvilinear Hole Having Three Poles and Arbitrary Shape. Bulletin of Calcutta Mathematical Society, 92, 309-322.
12. Abdou, M.A. (2002) Fundamental Problems for Infinite Plate with a Curvilinear Hole Having Finite Poles. Applied Mathematics and Computation, 125, 177-193. http://dx.doi.org/10.1016/S0096-3003(00)00118-1
13. Abdou, M.A., Sabbah, A.S. and Ismail, A.S. (2002) An Infinite Plate with a Curvilinear Hole and Flowing Heat. Proceeding Mathematical Physics Society Egypt, 34, 15-27.
14. Exadaktylos, G.E. and Stavropoulou, M.C. (2002) A Closed form Elastic Solution for Stress and Displacement around Tunnels. International Journal of Rock Mechanics and Mining Sciences, 39, 905-916.
15. Exadaktylos, G.E., Liolios, P.A. and Stavropoulou, M.C. (2003) A Semi-Analytical Elastic Stress-Displacement Solution for Notched Circular Openings in Rocks. International Journal of Solids and Structures, 40, 1165-1187.
16. Abdou, M.A. and Asseri, S.A. (2009) Closed Forms of Gaursat Functions in Presence of Heat for Curvilinear Holes. JJournal of Thermal Stresses, 32, 1126-1148.
17. Abdou, M.A. and Asseri, S.A. (2009) Gaursat Functions for an Infinite Plate with a Generalized Curvilinear Hole in Zeta Plane. Applied Mathematics and Computation, 212, 23-36. http://dx.doi.org/10.1016/j.amc.2009.01.079