 Open Access
 Authors : Ebikiton Ndiwari , Zuonaki Ongodiebi
 Paper ID : IJERTV9IS070581
 Volume & Issue : Volume 09, Issue 07 (July 2020)
 Published (First Online): 04082020
 ISSN (Online) : 22780181
 Publisher Name : IJERT
 License: This work is licensed under a Creative Commons Attribution 4.0 International License
Biharmonic Solution for the Forcing Term in a NonHomogeous Equation of Statics in the Theory of Elastic Mixture
Ebikiton Ndiwari1*, Zuonaki Ongodiebi2
*1, 2 Department of Mathematics and Computer Science, Niger Delta University, Nigeria
Abstract: In this paper, we analysed the problem of plane elasticity for a doubly connected body with outer and inner boundaries in the form of a regular polygon with a common centre and parallel sides. The sides of the polygon are exposed to external gravitational force and the biharmonic solution of the forcing term is determined. This is achieved by defining the forcing term in a non homogeneous equation of statics, using complex variable theory. The forces are analysed under twodimensional stress function, from which the equilibrium and compatibility equations were derived. Using the compatibility equation and the stressstrain relations, we derived our biharmonic equation. Our results show that the theoretical frame work of the forcing term is consistent with the previously existing results.
Keywords: Elasticity, biharmonic, forcing term, compactibility
1 INTRODUCTION
The application of the methods of conformal mappings and boundary value problems of analytic functions has proved to be the most effective way of solving boundary value problems of elasticity and plate bending. However, for a simplyconnected domain, these methods yield effective results (especially for domains mapped onto the circle by rational functions). However, these methods still remain poorly [1]. Nevertheless, for some practically important classes of doubly connected domains bounded by polygons including the polygonal domain with a curvilinear 2gonal holes, we may succeed in constructing effectively (in the analytical form) functions conformally mapping this domain onto the circular ring [2]. In addition to this, the KolosovMuskhelishvili methods make it possible to decompose these problems, into two RiemannHilbert problems, for the circular ring and by solving the latter problem, we can construct the sought complex potentials in analytic form.
Theories of mixtures in the framework of rational continuum thermodynamics have been developed throughout the sixties and seventies, and subsequent development in various constitutive theories and thermodynamic analysis are too numerous to document [3]. Boundaryvalue problems for a finite domain with part of its boundary being unknown and the other a polygonal line were solved in [4]. A similar boundaryvalue problems of plane elasticity for infinite plates weakened by unknown full strength holes with normal stresses on their boundaries and forces applied at infinity were analysed in [5, 6, 7].
A mixed problem of elasticity was solved in [8, 9] for a convex polygon and for a doubly connected domain with a polygonal boundary. Also, linear and nonlinear static boundaryvalue problems for doubly or multiply connected isotropic and anisotropic elastic bodies ( plates and shell) were solved by Maksimyuk and Chernyshenko [10], Liu,Ishih [7] discussed the entropy flux of transversely isotropic elastic bodies of homogeneous type [11], while [12, 13, 14] gave a solution of a non classical problem of oscillation of two component mixtures. A fundamental solution of the system of differential equations of stationary oscillations of twotemperature elastic mixtures theory was provided by [15].
The problem of plane elasticity for a doubly connected body with outer and inner boundaries in the form of a regular polygon with a common centre and parallel sides had been reported in literature; thereby motivating this research.
2 MATHEMATICAL FORMULATION
We consider a homogenous isotropic elastic body occupying a doubly connected domain on the complex plane
z x iy . Its outer and inner boundaries are L0 and L1 respectively. The body is a rectangle whose centre is at
z = 0. The neighborhood of the vertices of the inner rectangle symmetric angles of equidistance from the centre is shown in Fig 1.
Y
F F
1
F F
LOO
F F
Figure 1: isotropic Elastic Body
We assumed that the edges of the isotropic elastic body is exposed to an external force in the form of load. Our aim is to determine the biharmonic solution for the force (F), and the stress state.
3. METHOD OF SOLUTION
We use nonhomogeneous equation in the theory of elastic mixtures as our governing equation to define the forcing term (F) [1618]. It is shown that the displacement vector components are represented in this theory by means of four arbitrary analytic functions. In the twodimensional case, the basic nonhomogenous governing equation [16] of the theory of elastic mixture has the form:
au+ bgraddivu + cu + dgraddivu = F (1)
Where is the twodimensional laplacian, grad and div are the principal operators of the field theory, is the partial density (positive constant of the mixture), F is the mass force, u = w and u = w are the displacement vectors, a, b, c, and d are combination of constitutive constants, characterizing the physical properties of the mixture[19].

Theory of Complex Variables
We solve equation (1), using complex variable as follows:
= + with its conjugate given as = (2) which in addition results to
2 = + (3)
Expressing equation (3) in partial differential equation
2
= +
(4)
and by subtraction, equation (2) becomes,
2 = (5)
Expressing equation (5) in partial differential equation
2
=
(6)
Adding equation (4) and (6); we have,
2(
+ ) = 2
(7)
Subtracting equation (6) from (4); we have,
2(
) = 2
(8)
Multiplying equation (7) and (8); we obtain,
2
2
4
= 4(
+ ) + 4i(
) (9)
2
2
2 2
2 2
2
2
2
2
Equating the two right terms of equation (9) to the real and imaginary parts of equation (4) and (6) respectively; to have,
2
2
2
2
4(
+ ) = 4 (
+ ) (10i)
2
2
2
2
2
2
4i(
) = 4 (
) (10ii)
Replacing the two right hand terms of equation (9) by the two right hand terms of equation (10i) and (10ii); to have,
2
2
4
= 4 (
+ ) 4 (
) (11)
Let the displacement vectors u and u in their complex form be represented by
= u1 + iu2 and = u1 + iu2 (12)
Introducing the displacement vectors into equation (11); we have,
2
2
4
= 4 (
+ ) 4 ( ) (13)
Note: The displacement vector () depends on the elastic (z) and plastic () regions respectively. In order make equation (1) solvable, we adopt the method used in [1] and let
and
= 4
2
2
= 4
2
2
(14)
+
= 2 ( + ) = 2 = 2 (15)
Substituting equation (14) and (16) for in equation (1) ; we have,
2
2
<>4
+ 4 + 2 + 2 = (16)
2
2
Note: Our laplacian here is define as
= . =
where
.
= =
Thus, substituting for grad with in equation (16); we have,
2 2
2 2
4 + 4 + 2 + 2
= (17)
(4 + 4 + 2 + 2) = (18)
(19)
(19)
(4 + 4 + 2 + 2) =
4 + 4 + 2 + 2 = (20)
From Pompeiu formula, the integral [11] of (19); gives,
4 + 4 + 2 + 2 =
(21)
where is the analytic nonhomogeneous terms, define as,
= u + iv is the displacement vector component at the transformed state, as a result of the contact with external force.

Nonhomogeneous Part
Equating the nonhomogeneous parts of equations (21) and (1) That is,
= = (F
+ iF ) (21*)
So that
= (+)
1
= F
2
+ iF
(+) 1 2
which gives
+ +
= F1 + iF2 (22)
Equating coefficient in equation (22); to have,
= F1 (23)
+
= F2 (24)
Introducing new variables [15]; such that,
=
+
=
+
Substituting for u and v in equation (22) and (23) gives
( + ) ( + ) = F1
which implies
2 + 2
= F
2
2 1
2
2
(
2
2
+ ) = F
2
2 1
2 = F1 (25)
a similar procedure gives:
2 = F2 (26)
hence, from equation (25) and (26); we have,
2 + 2 = F1 + F2 (27)
2( + ) = (F1 + F2)
2 = F (28)
Note: The classical field theory describing gravity is the Newtonian gravitation, which describes the gravitational force F, as a mutual interaction between two masses M, and m. That is,
F =
2
where G is the earth gravitational constant, and r is the radius of the earth.
(29)
The massive body M has a gravitational field g. Since the gravitational force is conservative, the field g, can be written as a gradient of gravitational scalar potential ; that is,
G = (30)
Also, in the case of gravitational field due to an attracting massive object of density , Gauss law for gravity in differential form can be used to obtain the corresponding Poisson equation for gravity[14] as,
. = 4 (31)
Substituting for g in equation (31) using equation ((30)
. () = 4 (32)
2 = 4 (33)
Equation (33) is called Poisson equation for gravity[2]. Hence, equation (28) is equivalent to equation (33); because, they both involve the mutual interaction between two masses (M and m). That is
2 = F = 2 = 4 F = 4
F = 4 (34)
Thus, our forcing term (F) is the gravitational force, which is Poisson in nature; as such, it is restricted to a plane.

Biharmonic Solution for the Forcing Term (F)
We determine the stress state of the forcing term, by considering Euler Cauchy Stress Principle [20], which states that upon any surface (real or imaginary) that divides a body, the action of one part of the body on another is equivalent to the system of distributed forces, and is represented by a field Tn, called stress vector, defined on the surface A, assumed to depend continuously on the unit vector n [21].
This is mathematically as:
= ndA (35)
Thus,
dF = TndA
= n (36)
where Tn is the resultant stress vector, components n and n are the normal and shear stresses respectively. Hence, the stress distribution diagram becomes
ABBBBBBB
ABBBBBBB
A
A
Figure 2: Stress distribution on a plane
where, x and y are the normal stresses in the x and y directions.
xy = yx, is the symmetric shear stress in the x and y directions Hence, we deduce the equilibrium equation from figure 2, as:
+
= 0 (37i)
+ = 0 (37ii)
equations (37) is called the equilibrium equation.
Equation (37), involves two equations with three unknown,( , ). For compatibility, we introduce an extra equation from the straindisplacement relation of the deformation process [22].
=
(38ii)
(38i)
=
+
(38iii)
=
2
2
where u = u(, ) and v = v(, ) are the displacement vectors in the transformed state. Differentiating equation (38) twice with respect to , , respectively, result to
2 .
=
(39i)
2
2
2 .
2
=
(39ii)
2
2
2
2 2 2 2
2 2 2 2
= . + . = +
(39iii)
2
2
2
2
Hence, adding equations (39i), (39ii) and (39iii); we have,
2 + 2
2
=
(40)
2
2
equation (40) is called compatibility equation.
To solve equation (40), we use stressstrain relationship [22] for stress:
= 1 ( ) (41i)
= 1( ) (41ii)
= 2 (1 + )
= –
(41iii)
where:
= Poisson ratio
E = Young modulus
G = Modulus of rigidity
Substituting for , into equation (40); we have,
1 2
2
2 2
[ ( ) + ( )] =2 2
(1 + )
2 (
) + (
2
) = 2(1 + )
(42)
2
2
2
2
To eliminate the shearing term () in equation (42), we differentiate the equilibrium equation (37i) and (37ii) w.r.t.
respectively. That is,
2 + 2
= 0 (43i)
2
2 + 2
= 0 (43ii)
2
adding equation (43i) and (43ii); we have,
2 + 2 + 2 2
= 0 (44)
2
2
2
2 2
2 2
= 1 ( + ) (45)
2 2
2
2
2
Substituting for in equation (42); give,
2
2
1 2
2
2 ( ) +
2 ( ) = 2(1 + ) 2 ( 2
2 )
2
2
2
2
2 ( ) + 2 ( ) = (1 + ) ( 2
2 )
Expanding and simplifying the equation above; we have,
2
2
2
2
2 +
2 = 2
2
2 + 2 + 2 + 2 = 0 (46)
2
2
2
2
The solution to equation (46), can be obtained by introducing a new function , called Airys stress function [22]. For the case under consideration, we can define ; so that,
=
2
2 , =
2
2 , =
2
Substituting for , in equation (46); that is,
2 2 2 2 2 2 2 2
2 2 + 2 2 + 2 2 + 2 2 = 0
4
4
4
4 + 4 + 2 22 = 0
4
4
4
4 + 2 22 +
4 = 0
2 2 2 2
[( + ) ( + )] = 02 2 2 2
2. 2 = 0
4 = 0 (47)
equaion (47) is called biharmonic equation. Hence, our forcing term in equation (34) is,
= 4 = 4 = 0.

Stress State of the Forcing Term on the Plane
Stress is force per unit area; that is,
= Hence,
Force = Area Ã— Stress
With the above definition, figures 36 illustrate the stress distribution on the plane.
A
Figure 3
A
Figure 3
Â°
Â°
Figure 4
Â°
Â°
Deducing from figure 4, the force impact on the plane in xdirection is,
A. sin
A.
cos A
A. y
cos
A. sin
Figure 5: Stress distribution in the xdirection. Figure 6: Stress distribution in ydirection.
= ( cos + sin ) (48)
= cos +
sin (49)
Also, deducing from figure 6, the force impact on the plane in the ydirection is,
= sin + cos (50)
= ( sin +
cos ) (51)
Thus, total stress impact on the plane is,
= =
cos sin
= ( )
sin cos
Recall from our main result,
= 4 = 4 = 0
0 =
cos sin )
(
(
sin cos
0 = 2() 2 ()
2() = 2 ()
2 = = 2() (52)
2
Recall that is our resultant stress vector earlier mentioned in equation (36)
2 = 1 = 2() (53)
So that for
2 = 1
= 11 = 45Â° 2 = 90Â°
For
2 = 2()
= (54)
With the result of equation (54), we can generate the following table to ascertain the stress state of the following angles
Table 1: Relationship between and ()
Â°
0
30
60
90
120
150
()
0.00
0.58
1.73
0.00
0.58
1.73
Tn
2 *
*
*
1 *
0 30 60
1
2
*90
o
120 150
*
*
Figure 7: Graph of Stress() against Angle().

Discussion of Result
The graph in Fig.7, show the relationship between the angle () and the stress state of the isotropic elastic body. It is shown that the body attain stability at angle 90Â°; where, the resultant stresses resolved to zero. At angle 30Â° 0 60,Â°the body undergoes an elastic motion; while, at angle 120Â° 150Â°, the body suffers a plastic deformation.

Conclusion
In this paper, we considered the biharmonic solution for the forcing term (F) in a nonhomogeneous equation of statics in the theory of elastic mixture. It was found that our theoretical solution for the forcing term of the isotropic elastic body (Fig.1) is consistent with the existing result of [23]. In this paper, we analysed the problem of plane elasticity for a doubly connected body with outer and inner boundaries in the form of a regular polygon with a common centre and parallel sides. The sides of the polygon are exposed to external gravitational force and the biharmonic solution of the forcing term is determined. This is achieved by defining the forcing term in a nonhomogeneous equation of statics, using complex variable theory. The forces are analysed under twodimensional stress function, from which the equilibrium and compatibility equations were derived. Using the compatibility equation and the stressstrain relations, we derived our biharmonic equation. Our results show that the theoretical frame work of the forcing term is consistent with the previously existing results.
REFERENCES

Kapanadze G, Gulua B. About one problem of plane elasticity for a polygonal domain with a curvilinear hole. AMIM. 2016; 21:2129

Kapanadze G. A. On conformal mapping of doublyconnected domain bounded by convex polygon with linear section on circular ring. Bull. Georgian Acad. Sci. 2000; 61(2): 192194.

Rajagopal K. R, Tao L. Mechanics of mixtures. World Scientific Singapore; 1995.

Bantsuri R. D, Isakhanov AS. Some inverse problems of elasticity. Tbilisi Institute.1987;(87):320.

Bantsuri R. D. Some inverse problems of plane elasticity and bending of thin plates. Proc International symposium Dedicated to the centenary of Acad. N. Muskhelishvili, Tbilisi. 1993;(6):100 106.

Banichuk N. V. Shape of optimization for elastic bodies. Moscow: Nauka; 1980.

Lui, IShih. Entropy flux relation for viscoelastic bodies. Journal of Elasticity. 2008; 90:259 270.

Liu, Ishih. On entropy flux of transversely isotropic elastic bodies. Journal of Elasticity. 2009 96:97104.

Odishelidze N. T, Kriado F. F. A mixed problem of plane elasticity for a domain with unknown boundary. Internationally Applied Mechanics. 2006; 42(3): 342349.

Maksimyuk V. A, Chernyshenko IS. Mixed functional in the theory of non linearly elastic shells International Applied Mechanics. 2004; 40(11):12261262.

Kaloerov S. A, Boronenko OI. Twodimensional magnetoelastic problems for a connected piezomagnetic body. International Applied Mechanics. 2005; 41(10): 11381148.

Giorgashvili L, Skhvitaridze K. Problems of statics of twocomponent elastic mixtures. Math. J. 2005; 12(4): 619635. MathSciNet, Google Scholar.

Giorgashvilis L, Karseladze G, Sadunishvili G. Solution of a boundary value problem of statics twocomponent elastic mixtures for a space with two nonintersecting spherical cavities. Men. Differ. Equ. Math.Phys. 2008; 45:85115. MATHmathSciNet, Google Scholar.

Kapanadze G. A. On a problem of the bending of a plate for doubly connected domain bounded polygons. (Russian). Prikl Mat. Mekh. 66 no.4, 616620; Translation In J.Appl. Math. Mech 2002; 66(4): 601604.

Basheleishvili M. Analogies of the kolosovMuskheleishvili general representation formulas CauchyRiemann conditions in the theory of elastic mixtures. Georgian Mathematical Journal 1997;4(3):223242

Skhvitaridze K, Kharashvili M, Burchuladze D. Boundaryvalue problems of statics in the two temperature elastic mixture theory for a halfspace. Journal of Mathematical Sciences.2015; 206(4): 445456.

Giorgashvili L, Karseladze G, Kharashvili M. Basic boundaryvalue problems for a solid with double porosity and two nonintersecting spherical cavities. Journal of Mathematical Sciences; 2015.

Meladze R, Kharashvilin M, Skhvitaridz K. On a twocomponent elastic mixture with temperature values in a scalar field. Journal of Mathematical Sciences; 2015.

Natela Zirakashvili. On the numerical solution of some twodimensional boundary delocalization problems. Meccanica. 2013; 48(7): 17911804.

Truesdell AI, Toupin P. G. Classical and computational solid mechanics. Berlin World Scientific Springer Verlag. 2001; 209501.

Peter C. Continuum mechanics: Concise theory and problem. (New York). Dover series. 1999; pp.4750.

Chou T, Pagano G. Two and three dimensional stress function. 2nd ed. Aeche, Germany Institute of general mechanics. 2001; pp. 117.

Udoh P. J, Ndiwari E. Boundaryvalued equations for force term in non homogeneous of statics in the theory of elastic mixture. Asian Research Journal of Mathematics. 2018; 8(1 111.
