`ISRN Civil EngineeringVolume 2013 (2013), Article ID 562482, 8 pageshttp://dx.doi.org/10.1155/2013/562482`
Research Article

## Poisson’s Theory for Analysis of Bending of Isotropic and Anisotropic Plates

Department of Aerospace Engineering, Indian Institute of Science, Bangalore 560 012, India

Received 30 May 2013; Accepted 19 June 2013

Academic Editors: M. Garg and D. Huang

Copyright © 2013 K. Vijayakumar. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

#### Abstract

Sixteen-decade-old problem of Poisson-Kirchhoff’s boundary conditions paradox is resolved in the case of isotropic plates through a theory designated as “Poisson’s theory of plates in bending.” It is based on “assuming” zero transverse shear stresses instead of strains. Reactive (statically equivalent) transverse shear stresses are gradients of a function (in place of in-plane displacements as gradients of vertical deflection) so that reactive transverse stresses are independent of material constants in the preliminary solution. Equations governing in-plane displacements are independent of the vertical (transverse) deflection . Coupling of these equations with is the root cause for the boundary conditions paradox. Edge support condition on does not play any role in obtaining in-plane displacements. Normally, solutions to the displacements are obtained from governing equations based on the stationary property of relevant total potential and reactive transverse shear stresses are expressed in terms of these displacements. In the present study, a reverse process in obtaining preliminary solution is adapted in which reactive transverse stresses are determined first and displacements are obtained in terms of these stresses. Equations governing second-order corrections to preliminary solutions of bending of anisotropic plates are derived through application of an iterative method used earlier for the analysis of bending of isotropic plates.

#### 1. Introduction

Kirchhoff’s theory [1] and first-order shear deformation theory based on Hencky’s work [2] abbreviated as FSDT of plates in bending are simple theories and continuously used to obtain design information. Kirchhoff’s theory consists of a single variable model in which in-plane displacements are expressed in terms of gradients of vertical deflection so that zero face shear conditions are satisfied. is governed by a fourth-order equation associated with two edge conditions instead of three edge conditions required in a 3D problem. Consequence of this lacuna is the well-known Poisson-Kirchhoff boundary conditions paradox (see Reissner’s article [3]).

Assumption of zero transverse shear strains is discarded in FSDT forming a three-variable model. Vertical deflection and in-plane displacements are coupled in the governing differential equations and boundary conditions. Reactive (statically equivalent) transverse shears are combined with in-plane shear resulting in approximation of associated torsion problem instead of flexure problem. In Kirchhoff's theory, in-plane shear is combined with transverse shear implied in Kelvin and Tait’s physical interpretation of contracted boundary condition [4]. In fact, torsion problem is associated with flexure problem whereas flexure problem (unlike directly or indirectly implied in energy methods) is independent of torsion problem. Reissner [3] in his article felt that inclusion of transverse shear deformation effects was presumed to be the key in resolving Poisson-Kirchhoff’s boundary conditions paradox. Recently, it is shown that the second-order correction to by either Reissner’s theory or FSDT corresponds to approximate solution of a torsion problem [5]. Coupling between bending and associated torsion problems is eliminated earlier [6] by using zero rotation about the vertical axis. The proposal of zero does not, however, resolve the paradox in a satisfactory manner. In the reference [6], face deflection in (8) is from zero face shear conditions whereas neutral plane deflectio n from edge support condition is associated with self-equilibrating transverse shear stresses.

The condition decoupling the bending and torsion problems is satisfied in Kirchhoff’s theory. If this condition is imposed in FSDT, sum of the strains () in the isotropic plate is governed by a second-order equation due to applied transverse loads. Reactive transverse shear stresses are in terms of gradients of () uncoupled from . Thickness-wise linear distribution of is zero from the equilibrium equation governing transverse stresses. Normal strain from constitutive relation is linear in in terms of (). Reactive transverse shear stresses and thickness-wise linear strain form the basis for resolving the paradox and for obtaining higher order corrections to the displacements. The theory thus developed is designated as “Poisson’s theory of plates in bending”.

The previously mentioned Poisson’s theory is applied to the analysis of anisotropic plates using reactive transverse shear stresses as gradients of a function. Normally, solutions to the displacements are obtained from governing equations based on stationary property of relevant total potential and the reactive transverse shear stresses are expressed in terms of these displacements. In the present work, reverse process is adapted in which reactive transverse stresses are determined first and the displacements are obtained in terms of these stresses. Equations governing second-order corrections to the preliminary solutions are derived through the application of an iterative method used earlier [7] for the analysis of bending of isotropic plates.

#### 2. Equations of Equilibrium and Edge Conditions

For simplicity in presentation, a rectilinear domain bounded by , , , with reference to Cartesian coordinate system () is considered. For convenience, coordinates and displacements in nondimensional form , , , , and half-thickness ratio with reference to a characteristic length in - plane are used ( is defined such that mod of and are equal to or less than 1). With the previous notation and indicating interchange, equilibrium equations in terms of stress components are in which suffix after “,” denotes partial derivative operator.

Edge conditions are prescribed such that () are even in , and () are odd in . In the primary flexure problem, the plate is subjected to asymmetric load and zero shear stresses along faces. Three conditions to be satisfied along (and ) constant edges are prescribed in the form

#### 3. Stress-Strain and Strain-Displacement Relations

In displacement based models, stress components are expressed in terms of displacements, via six stress-strain constitutive relations and six strain-displacement relations. In the present study, these relations are confined to the classical small deformation theory of elasticity.

It is convenient to denote displacement and stress components of anisotropic plates as

Strain-stress relations in terms of compliances [] with the usual summation convention are We have from semi-inverted strain-stress relations with [] denoting inverse of []: Normal strain is given by Transverse shear strains with [] denoting inverse of [] are Strains from strain-displacement relations are

#### 4. Functions and Their Use

We use thickness-wise distribution functions generated from recurrence relations [7] with , , such that . They are (up to ) Displacements, strains, and stresses are expressed in the form (sum ) (Note that ,  ) In the previous equations, variables associated with functions are functions of only.

#### 5. Sequence of Trivially Known Steps without Any Assumptions in Preliminary Analysis

(a) Due to prescribed zero along faces of the plate, in the plate.

(b) from constitutive relations (5).

(c) Static equilibrium equation (2) gives .

(d) In the absence of (), from strain-displacement relations (8).

(e) from constitutive relations (6).

(f) from relation (7).

(g) from thickness-wise integration of .

(h) From steps (b) and (g), integration of in the thickness direction gives .

In the present study as in the earlier work [6], linear in are unknown functions as in FSDT. From integration of in the face plane instead of thickness direction, one gets vertical deflection in the form Since is from satisfaction of zero face shear conditions, it can be considered as face deflection though it is same for all face parallel planes (note that prescribed zero along a segment of the neutral plane implies zero along the corresponding segment of the intersection of face plane with wall of the plate since is independent of and vice versa). It is analytic in the domain of the plate if is zero. In such a case, we note from strain-displacement relations (8) that

(i) In the absence of applied transverse loads, integrations of (1), (2) using Kirchhoff’s displacements give a homogeneous fourth-order equation governing . In the present analysis as in FSDT, static equations governing are, however, given by They are subjected to the edge conditions (3a), (3b) along - (and -) constant edges.

Equations (15) governing are uncoupled due to relations (14) in the case of isotropic plates whereas they are coupled without cross-derivatives in the case of anisotropic plates. In the absence of prescribed bending stress all along the closed boundary of the plate, in-plane distribution corresponds to pure torsion where as it is different in the presence of bending load. These distributions are discussed later in Section 6.2.

(j) In FSDT, (15) correspond to neglecting shear energy from transverse shear deformations. Vertical deflection has to be obtained from (14) using the solutions for from (15). Reactive transverse stresses are expressed in terms of . They are dependent on material constants different from prescribed nonzero transverse shear stresses along the edges of the plate.

(k) In the development of Poisson’s theory, (15) are coupled with reactive transverse stresses.

#### 6. Poisson’s Theory

In the preliminary solution, steps (a)–(g) in Section 5 are unaltered. That is, transverse stresses and strains are zero and . In the absence of higher order in-plane displacement terms, reactive transverse shear stresses are parabolic and gradients of a function , that is, in-plane distributions of transverse shear stresses . Equation governing from thickness-wise integration of (2) is In the previous equation, is coefficient of . Satisfaction of the load condition along faces gives so that Edge condition on is either or its outward normal gradient is equal to the prescribed shear stress along each segment of the edge. Note that the transverse stresses thus obtained are independent of material constants. In the isotropic plate, in which and .

Integrated equilibrium equations (1) in terms of in-plane strains are (sum ): By substituting strains in terms of , (18) are two equations governing and have to be solved with conditions along - (and -) constant edge

Vertical deflection is given by (13). As mentioned earlier, cross-derivatives of do not exist in (18), (19a), and (19b) due to (14). In the isotropic case, (18) are uncoupled and are simply given by Edge conditions (19a) and (19b) are also uncoupled and they are given by The condition in solving (17) is different from the usual condition in the Kirchhoff’s theory and FSDT. The function is related to the normal strain . In the isotropic case, is proportional to . Gradients of are proportional to gradients of transverse strains in FSDT. If the plate is free of applied transverse stresses, that is, the plate is subjected to bending and twisting moments only, is proportional to . The function , thereby, is identically zero from (2) (in Kirchhoff’s theory, the tangential gradient of is proportional to the corresponding gradient of applied along the edge of the plate implied from Kelvin and Tait’s physical interpretation of the contracted transverse shear condition). This Laplace equation is not adequate to satisfy two in-plane edge conditions. One needs its conjugate harmonic functions to express in-plane displacements in the form The function was introduced earlier by Reissner [8] as a stress function in satisfying (2). Note that the two variables and are expressed in terms of gradients of two functions and .

After finding and from solving in-plane equilibrium equations (1), correction to from zero face shear conditions is given by One should note here that does not satisfy prescribed edge condition on .

In FSDT, are obtained from solving (1). There is no provision to find and it has to be obtained from zero face shear conditions in the form If the applied transverse stresses are nonhomogeneous along a segment of the edge, is coupled with in (1) and (2) and edge conditions (3a), (3b), and (3c).

With reference to the present analysis, it is relevant to note the following observation: in the preliminary solution, , , are of . As such, estimation of in-plane stresses, thereby, in-plane displacements is not dependent on . The support condition along the edge of the plate does not play any role in determining the in-plane displacements.

##### 6.1. Comparison with Analysis of Extension Problems

It is interesting to compare the previous analysis with that of extension problem. In extension problem, applied shear stresses along the faces of the plate are gradients of a potential function and they are equivalent to linearly varying body forces, independent of elastic constants, in the in-plane equilibrium equations. Rotation (but ) and equations governing displacements are coupled. Also, vertical deflection is given by thickness-wise integration of from constitutive relation and higher order corrections are dependent on this vertical deflection. In the present analysis of bending problems, in-plane distributions of transverse shear stresses are gradients of a function related to applied edge transverse shear and face normal load through equilibrium equation (2). They are also independent of elastic constants and derivative of parabolic distribution of each of them is equivalent to body force in linearly varying in-plane equilibrium equations. Vertical deflection is purely from satisfaction of zero shear stress conditions along faces of the plate. Higher order corrections are dependent on normal strains and independent of . It can be seen that steps in the analysis of bending problem are complementary to those in the classical theory of extension problem.

##### 6.2. Associated Torsion Problems

First-order shear deformation theory and higher order shear deformation theories do not provide proper corrections to initial solution (from Kirchhoff’s theory) of primary flexure problems. In these theories, corrections are due to approximate solutions of associated torsion problems. In fact, one has to have a relook at the use of shear deformation theories based on plate (instead of 3D) element equilibrium equations other than Kirchhoff’s theory in the analysis of flexure problems. In the earlier work [5], we have mentioned that this coupling of bending and torsion problems is nullified in the limit of satisfying all equations in the 3D problems. However, these higher order theories in the case of primary flexure problem defined from Kirchhoff’s theory are with reference to finding the exact solution of associated torsion problem only. In fact, one can obtain exact solution of the torsion problem (instead of using higher order polynomials in by expanding and in sine series and in cosine series. In a pure torsion problem, normal stresses and strains are zero implying that Rotation and warping displacements () are independent of vertical displacement . In FSDT, and corrections to Kirchhoff’s displacements are due to the solution of approximate torsion problem. In the isotropic rectangular plate, warping function is harmonic in --plane. It has to satisfy with prescribed along an constant edge. By expressing in product form , we have By taking where satisfying zero face shear conditions, we get If prescribed is with specified constants (), clearly there is no provision to satisfy zero in-plane shear condition along constant edges. It has to be nullified with corresponding solution from bending problem. That is why a torsion problem is associated with bending problem but not vice versa. In a corresponding bending problem, all stress components are zero except . This solution is used only in satisfying edge condition along constant edges in the presence of specified along constant edges. It shows that distribution in flexure problem is nullified in the limit in shear deformation theories due to torsion.

It is interesting to note that the previous mentioned deficiency, due to coupling with in the plate element equilibrium equations in FSDT and higher order shear deformation theories, does not exist if applied is zero all along the closed boundary of the plate. It is complementary to the fact that boundary condition paradox in Kirchhoff’s theory does not exist if tangential displacement and are zero all along the boundary of the plate.

##### 6.3. Illustrative Example: Simply Supported Isotropic Square Plate

Consider a simply supported isotropic square plate subjected to vertical load Exact values for neutral plane and face deflections obtained earlier [7] with , are In this example, Poisson-Kirchhoff’s boundary conditions paradox does not exist. Hence, Poisson’s theory gives the same solution obtained in the author’s previous mentioned work. Concerning face value in (30), Kirchhoff’s theory gives a value of 2.27 that is same for all face parallel planes. Poisson’s theory gives an additional value of 0.267 attributable to and it is same for neutral plane deflection since can be expressed as In the previous expression, is deflection of the neutral plane. Higher order correction to uncoupled from torsion is 1.262 so that total correction to the value from Kirchhoff’s theory is about 1.53 [7]. Correction due to coupling with torsion is 1.45 [5] whereas it is 1.423 from FSDT and other sixth order theories [9]. With reference to numerical values reported by Lewiński [9], the previous correction is 1.412 in Reddy’s 8th-order theory and less than 1.23 in Reissner’s 12th-order and other higher order theories. It clearly shows that these shear deformation theories do not lead to solutions of bending problems (note that the value () from FSDT corresponds to neutral plane deflection and it is in error by about 17.7% from the exact value).

In view of the previous observations, it is relevant to make the following remarks: physical validity of 3D equations is constrained by limitations of small deformation theory. For , this constraint is dependent on material constants, geometry of the domain, and applied loads. Range of validity of 2D plate theories arises for small values of . Kirchhoff’s theory gives lower bound for this . For this value of and for slightly higher values of than this lower bound, Kirchhoff’s theory is used, in spite of its known deficiencies, due to its simplicity to obtain design information. Similarly, FSDT and other sixth-order shear deformation theories are used for some additional range of values of due to the corrections over Kirchhoff’s theory though these corrections are due to coupling with approximate torsion problems. They serve only in giving guidance values for applicable design parameters for small range of values of beyond its lower bound. They do not provide initial set of equations in formulating proper sequence of sets of 2D problems converging to the 3D problem. As such, comparison of solutions from Kirchhoff theory and FSDT with the corresponding solutions in the present analysis does not serve much purpose.

Since has to be higher than , one has to obtain the correction to before finding correction to internal distribution of along with correction to . Initial correction to is obtained earlier [7] in the form It is the same for face deflection (note that in the previously mentioned author’s work is required only to obtain ). As such, is further corrected from solution of a supplementary problem. This total correction over face deflection is about 0.658 giving a value of 4.458 which is very close to the exact value 4.487 in (29). However, we note that the error in the estimated value (=3.8) of is relatively high compared to the accuracy achieved in the neutral plane deflection. It is possible to improve estimation of by including in () such that () are independent of . In the present example, correction to face deflection changes to In the previous equation, . Correction to face value changes to 1.431 giving 3.97 () for face deflection which is under 4.7% from exact value.

#### 7. Anisotropic Plate: Second-Order Corrections

We note that does not participate in the determination of and it is obtained from the constitutive relation (6) in the interior of the plate. It is zero along an edge of the plate if is specified condition. With zero transverse shear strains, specification of instead of is more appropriate along a supported edge since implies zero tangential displacement along a straight edge which is the root cause of Poisson-Kirchhoff’s boundary conditions paradox. As such, edge support condition on w does not play any role in obtaining the in-plane displacements and reactive transverse stresses.

In the first stage of iterative procedure, second-order reactive transverse stresses have to be obtained by considering higher order in-plane displacement terms and . However, one has to consider limitation in the previous Poisson’s theory like in Kirchhoff’s theory; namely, reactive is zero at locations of zeros of . It is identically zero in higher order approximations since functions are not zero along faces. Moreover, it is necessary to account for its dependence on material constants. To overcome these limitations, it is necessary to keep as a free variable by modifying in the form In the previous equation, so that . Denoting coefficient of in by , it becomes Displacements are modified such that they are corrections to face parallel plane distributions of the preliminary solution and are free to obtain reactive stresses , , and normal strain .

We have from constitutive relations Modified displacements and the corresponding derived quantities denoted with are Strain-displacement relations give In-plane and transverse shear stresses from constitutive relations are Suffix (3) in (39) is to indicate that the stresses and strains correspond to and .

We get from (2), (15), (35), (37), and (40) with reference to coefficient of : Note that is not present in the previous equation.

From integration of equilibrium equations, we have with sum on Note thatis present in the previous expression. One gets one equation governing in-plane displacements from (41), (43), noting that , in the form The second equation governing these variables is from the condition . Here, it is more convenient to express as so that and (45) becomes a fourth-order equation governing . This sixth-order system is to be solved for and with associated edge conditions along an constant edge and analogous conditions along a constant edge.

In (45) expressed in terms of and , we replace by for obtaining neutral plane deflection since contributions of and are one and the same in finding .

##### 7.1. Supplementary Problem

Because of the use of integration constant to satisfy face shear conditions, vertical deflection of neutral plane is equal to that of the face plane deflection. This is physically incorrect since face plane is bounded by elastic material on one side whereas neutral plane is bounded on both sides. This lack in interior solution is rectified by adding the solution of a supplementary problem based on leading cosine term that is enough in view of distribution of reactive shear stresses. This is based on the expansion of in sine series. Here, load condition is satisfied by the leading sine term and coefficients of all other sine terms are zero.

Corrective displacements in the supplementary problem are assumed in the form We have from integration of equilibrium equations

In-plane distributions and are added as corrections to the known and so that in the supplementary problem are We note that from integration of equilibrium equations with variables is the same as from static equilibrium equations with variables. Hence, It is convenient to use here also Correction to neutral plane deflection due to solution of supplementary problem is from

##### 7.2. Summary of Results from the First Stage of Iterative Procedure

The previous analysis gives displacements , consistent with , and corrective transverse stresses from the first stage of iteration It is to be noted that the dependence of transverse stresses on material constants is through the solution of supplementary problem. One may add terms in shear components and term in normal stress component. These components are also dependent on material constants, but they need correction from the solution of a supplementary problem. Successive application of the previous iterative procedure leads to the solution of the 3D problem in the limit.

#### 8. Concluding Remarks

Poisson’s theory developed in the present study for the analysis of bending of anisotropic plates within small deformation theory forms the basis for generation of proper sequence of 2D problems. Analysis for obtaining displacements, thereby, bending stresses along faces of the plate is different from solution of a supplementary problem in the interior of the plate. A sequence of higher order shear deformation theories lead to solution of associated torsion problem only. In the preliminary solution, reactive transverse stresses are independent of material constants. In view of layer-wise theory of symmetric laminated plates proposed by the present author [10], Poisson’s theory is useful since analysis of face plies is independent of lamination which is in confirmation of (8)-(9) in a recent NASA technical publication by Tessler et al. [11].

One significant observation is that sequence of 2D problems converging to 3D problems in the analysis of extension, bending, and torsion problems are mutually exclusive to one other.

#### Highlights

(i)Poisson-Kirchhoff boundary conditions paradox is resolved. (ii)Transverse stresses are independent of material constants in primary solution.(iii)Edge support condition on vertical deflection has no role in the analysis.(iv)Finding neutral plane deflection requires solution of a supplementary problem.(v)Twisting stress distribution in pure torsion nullifies its distribution in bending.

#### Appendix

= side length of a square plate = plate thickness = through-thickness distribution functions, = characteristic length of the plate in --plane = applied face load density = stiffness coefficients = elastic compliances = prescribed stress distributions along an constant edge = displacements in -directions, respectively = = Cartesian coordinate system = nondimensional coordinates = Laplace operator = in-plane strains = transverse strains , rotation about -axis = bending stresses = transverse stresses.

#### References

1. G. Kirchhoff, “Über das Gleichgewicht und die Bewegung einer elastischen Scheibe,” Journal für die Reine und Angewandte Mathematik, vol. 40, pp. 51–58, 1850.
2. H. Hencky, “Über die Berücksichtigung der Schubverzerrung in ebenen Platten,” Ingenieur-Archiv, vol. 16, no. 1, pp. 72–76, 1947.
3. E. Reissner, “Reflections on the theory of elastic plates,” Applied Mechanics Reviews, vol. 38, no. 11, pp. 1453–1464, 1985.
4. A. E. H. Love, A Treatise on Mathematical Theory of Elasticity, Cambridge University Press, Cambridge, UK, 4th edition, 1934.
5. K. Vijayakumar, “Modified Kirchhoff’s theory of plates including transverse shear deformations,” Mechanics Research Communications, vol. 38, pp. 211–213, 2011.
6. K. Vijayakumar, “New look at kirchhoff's theory of plates,” AIAA Journal, vol. 47, no. 4, pp. 1045–1046, 2009.
7. K. Vijayakumar, “A relook at Reissner’s theory of plates in bending,” Archive of Applied Mechanics, vol. 81, no. 11, pp. 1717–1724, 2011.
8. E. Reissner, “The effect of transverse shear deformations on the bending of elastic plates,” Journal of Applied Mechanics, vol. 12, pp. A69–A77, 1945.
9. T. Lewiński, “On the twelfth-order theory of elastic plates,” Mechanics Research Communications, vol. 17, no. 6, pp. 375–382, 1990.
10. K. Vijayakumar, “Layer-wise theory of bending of symmetric laminates with isotropic plies,” AIAA Journal, vol. 49, no. 9, pp. 2073–2076, 2011.
11. A. Tessler, M. Di Sciuva, and M. Gherlone, “Refined zigzag theory for homogeneous, laminated composite, and sandwich plates: a homogeneous limit methodology for zigzag function selection,” Tech. Rep. NASA/TP-292010216214:1, 2010.