Analytical Behavior Prediction for Skewed Thick Plates on Elastic Foundation

This paper presents analytical solutions for the problem of skewed thick plates under transverse load on a Winkler foundation, which has not been reported in the literature. The thick plate solution is obtained by using a framework of an oblique coordinate system. First, the governing differential equation in that system is derived, and the solution is obtained using deflection and rotation as derivatives of the potential function developed here. This method is applicable for arbitrary loading conditions, boundary conditions, and materials. The solution technique is applied to two illustrative application examples, and the results are compared with numerical solutions. The two approaches yielded results in good agreement.


Introduction
Plates resting on an elastic foundation are important structural elements that are used in a wide range of applications including foundation structures and concrete pavements of roads.In the past few decades, there have been efforts to analytically investigate the behavior of plates resting on an elastic foundation.Most of the work has been performed on rectangular plates and circular plates.For example, Timoshenko and Woinowsky-Krieger 1 provided an analytical solution for rectangular thin plates, and Celep 2 presented the behavior of circular thin plates on an elastic foundation.
Due to its mathematical complexity, there are very few studies which have considered skewed plates on an elastic foundation.Kennedy 3 was the first to develop an analytical solution as polynomial series for skewed isotropic thin plates with four edges clamped on an elastic foundation.Ng and Das 4 and Chell et al. 5 developed an analytical solution as trigonometric series for skewed sandwich thin plates on an elastic foundation.Though these studies are informative, they have limitations in terms of loading and Mathematical Problems in Engineering In addition, there is another issue; these methods are applicable only to thin plates of which the thickness-to-side ratio is less than 1/20 because Kirchhoff plate theory has been employed to develop the governing equation in these studies.Though the Kirchhoff theory is widely used in plate analysis, it suffers from under predicting deflections when thick plates are analyzed because it neglects the effect of the transverse shear deformation 6 .To resolve this issue, the Mindlin theory 7, 8 , which relaxes the perpendicular restriction for transverse normals, was developed so that the effect of transverse shear deformation could be considered.Several studies which have employed the Mindlin theory for analysis are worth mentioning.For example, Kobayashi and Sonoda 9 developed an analytical solution for rectangular thick plates with two opposite edges, simply supported.In their research, Levytype trigonometric series solutions were derived.Liew et al. 10 and Liu 11 solved the rectangular thick plate problem by the differential quadrature method and offered a number of numerical results under several boundary conditions.Though there are some studies dealing with rectangular thick plates, there has been no research analyzing skewed thick plates.
In the present manuscript, analytical solutions for skewed thick plates on an elastic foundation with arbitrary boundary and loading conditions are newly developed based on the Mindlin theory.The Winkler foundation, which is the simplest model to represent the elastic foundation, is employed here.The present method is new in the following two aspects.First of all, the present method has developed analytical solutions based on the Mindlin theory for skewed thick plates, and such solutions have never previously been developed in the literature.Second, the present method allows arbitrary loading conditions, boundary conditions, and materials.We first show the derivation of a governing equation in a framework of an oblique coordinate system, and then it is solved using trigonometric series.The governing equation in this system is derived and the solution is obtained using the deflection and rotation as derivatives of a potential function.The solution technique is applied to an illustrative application example with different skew angles and modulus of foundations, and the results are compared with the commercial finite element package ANSYS 12 .The present approach is in reasonable agreement with the FEM solution.

Governing Equation under the Oblique Coordinate System
When a plate's boundary profile is a parallelogram, the oblique Cartesian coordinate system can be advantageous.We first present the concept of the oblique coordinate system and then show the governing differential equation of skewed thick plates on a Winkler foundation based on the Mindlin theory.The authors have developed an analytical solution for skewed thick plates without an elastic foundation in 13 , and the formulas developed in this Mathematical Problems in Engineering chapter have utilized the relationships developed in the research.Figure 1 shows the oblique coordinate system spanned by the X and Y axes, along with the rectangular coordinate system x and y, with angle XOY denoted as skew angle α.Parallelogram ABCD in Figure 1 represents the skewed plate of interest, and the edge lengths CD and AD are 2a and 2b, respectively.
The relationship between the rectangular and oblique coordinate systems can be written as follows 14, 15 : where φ X and φ Y are, respectively, the rotations normal to the X and Y axes.The corresponding relationships for stress, strain, moment, and shear force under the two coordinate systems are also available in 14 .They are readily derivable according to the principles of statics and continuum mechanics.Their corresponding quantities in the oblique system are derived in this research effort.The resulting flexural stiffness matrix under the oblique coordinate system D o is related to its counterpart in the rectangular system D r as follows 16 :

2.2
The flexural stiffness matrices relate the moments to the curvatures in the respective coordinate systems.For example, D r in the rectangular coordinate system for isotropic material is 1 : where E is Young's modulus, ν is Poisson's ratio, and t is the thickness of the plate.Note that 2.2 is also applicable for other, more complex situations, such as orthotropic or anisotropic materials.By using the flexural stiffness matrices, the moment-strain relationship of the rectangular and oblique coordinate systems can be described as in the following 2.4 and 2.5 , respectively 6 : In this derivation, the Mindlin theory is applied, which assumes that strains are linearly distributed along the thickness direction in the plate's cross-section.Similarly, the extensional stiffness matrix in the oblique coordinate system A O is related to that of the rectangular system A r as follows 13 : The extensional stiffness matrix relates the shear forces to the shear strains.For example, A r for isotropic material is Mathematical Problems in Engineering By using the extensional stiffness matrices, the relationships between the shear force and deflection, rotation angle of the rectangular and oblique coordinate system can be described as in the following 2.8 and 2.9 , respectively, 6 : where w is the transverse deformation perpendicular to the plane of the plate, and K s is the shear correction factor to account for nonuniform transverse shear distribution.Hereafter, the components in D o and A o are referred to using their respective elements D 11 to D 33 and A 44 to A 55 as follows: where the diagonal components of D O relate the moments to the curvatures in the same directions.The off-diagonal terms relate the same moments to the curvatures in other directions due to the Poisson's effect and coordinate system obliquity.Similarly, the diagonal components of A o relate the shear forces to the shear strains in the same directions, and off-diagonal terms to the shear strains in other directions due to obliquity.The following 2.11 to 2.13 are equilibrium conditions of the skewed plates shown in Figure 1.
Equilibrium of force in the z direction perpendicular to the plane of the plate Equilibrium of moments along the X-axis Equilibrium of moments along the Y-axis  where Q in 2.11 is the load normal to the upper surface of the plate, and k is the modulus of the foundation.By substituting 2.5 , 2.9 , and 2.10 into 2.11 to 2.13 , the following 2.14 to 2.16 are obtained in the oblique system:

2.16
To make the solution process simpler, a new function ψ is introduced below to represent the condition of the skewed thick plate.We assume that w consists of terms up to the 4th

2.19
w has only even derivatives and φ X and φ Y have only odd derivatives though it is assumed that w, φ X , and φ Y may have both even and odd derivatives.This is because the coefficients of odd derivatives of w and the coefficients of even derivatives of φ X and φ Y should be 0 to satisfy 2.14 to 2.16 .By substituting these relations into 2.14 , the governing equation of the Mindlin skewed thick plate is then formulated as a 6th-order partial differential equation as follows: where L is a linear differential operator in the oblique coordinate system

Analytical Solution in Series Form
In this section, a general solution to the governing differential equation 2.11 is developed as the sum of a fundamental homogeneous and a particular nonhomogeneous solution, detailed separately, as follows.

Homogeneous Solution
The homogeneous solution ψ h is the solution to 2.20 for Q 0, obtained as a trigonometric series in 3.1 below 3.1 where i √ −1 is the imaginary unit, and C eXf , C eY f , S eXf , and S eY f trigonometric functions are as follows: where the bar on λ denotes the conjugate of λ. λ 1X , λ 2X , λ 3X , λ 1Y , and λ 2Y , and λ 3Y are the eigenvalues to be obtained by satisfying L ψ h 0. For example, λ eX is derived by solving the following equation: The trigonometric function ψ h in 3.1 has 24l unknowns A h , B h , C h , . .., and X h h 1, 2, 3, . . ., l , with l being the number of the trigonometric terms needed for convergence.In addition, as for the coefficients of real parts A h , C h , . .., and W h , l 1 terms are considered so that 12 more unknowns A l 1 , C l 1 , . .., and W l 1 can be determined.Therefore, the homogeneous solution ψ h has 24l 12 unknowns, and these will be determined according to the boundary conditions discussed below.

Particular Solution
For a particular solution in the series form, the transverse load Q in 2.20 is expanded to a trigonometric series as follows:

3.4
Mathematical Problems in Engineering  where K jk and L jk are to be determined to satisfy 2.20 and 3.4 , and m is the number of the trigonometric terms needed for convergence.The general solution for ψ is derived as the sum of the homogeneous solution and the particular solution as ψ ψ h ψ p .

3.6
Since no unknowns exist in the particular solution, the total number of unknowns in the general solution is still 24l 12, as in the homogeneous solution.

Determination of Unknown Constants for Series Solution
Using the Mindlin theory, the boundary conditions for various edges are given below for determining the unknown constants in the homogeneous solution.The normal and tangential directions to the edge are denoted here using the subscripts n and s, respectively.The moments on the edges are noted using these subscripts in a way that is consistent with Mathematical Problems in Engineering

X/2a
Analytical the directions of the stresses thereby induced.Namely, M n is for the moment causing normal stresses and M s is the torsional moment inducing shear stresses.3 Hard simply supported SS2 : w 0, M n 0, M s 0. 4.3 4 Free F : Note that the Kirchhoff theory treats SS1 and SS2 in 4.2 and 4.3 as having the same boundary condition.The difference between them is explained graphically in Figure 2. The boundary condition of SS1 restricts the tangential rotation by supporting two points in the cross-section, thereby generating a nonzero torsional moment.In contrast, the boundary condition of SS2 supports the plate only at one point in the cross-section, allowing a tangential rotation and generating no twisting moment.The boundary conditions in 4.1 to 4.4 can be unified as follows: where Γ 1 X, Y to Γ 12 X, Y represent the left hand side of 4.1 to 4.4 .Γ 1 X, Y to Γ 12 X, Y are expanded as a Fourier series as follows for the solution method pursued in this paper: where coefficients a 0d , a cd , and b cd are Fourier coefficients for the boundary condition Γ d X, Y .Because of the need for truncation, l terms are kept for each of the 12 boundary conditions so that a total of 12 2l 1 equations are established for the same number of Fourier coefficients that have the 24l 12 unknowns included to be solved for.

Application Examples
In this section, two application examples are presented using the developed analytical solution for skewed thick plates.They are also compared with solutions published in 10 and the FEM analysis result obtained using the commercial package ANSYS 11.In the analysis by ANSYS, 2D 8-node quadrilateral plate elements SHELL99 appropriate for thick plate analysis are used for the skewed plates with various skewed angles.The numbers of nodes and elements are 4961 and 1600, respectively.The convergence figure of transverse displacements along with the number of nodes is provided for both of the application examples.In addition, the effect on convergence of the number of terms l and m in the fundamental and particular solutions is studied.In the following examples, the shear correction factor K s is taken to be 5/6 since this factor is commonly used in plate analyses.

Isotropic Clamped Skewed Thick Plate under Uniformly Distributed Load
Isotropic skewed thick plates are analyzed here.As an external force, uniformly distributed load q 0 is applied, and, as skewed angles, α 30 • , 60 As a first step, convergence of the finite element model is checked.Figure 3 shows the convergence figure of transverse displacements w at the center of the plates when α is 30 • .It can be concluded that the finite element model converges adequately when the number of elements 1600.Similar results are confirmed for other α.
Then, the numbers of terms in the series solution l and m are determined.The deflection w at the center of the plates for l, m 5, 55 and 7, 35 differs by less than 0.1% from that of l, m 7, 55 as in Figure 4 when α 30 • and K 7. It can be concluded that the solution is already convergent while truncated at l, m 7, 55 , and similar results are observed for other α and K. Therefore, l 7 and m 55 are employed in this example.
Figures 5, 6, and 7 give a comparison between the present method and FEM analysis using ANSYS for the deflection w for α 30 • , 60 • , and 90 • along line HF defined in Figure 1 and on the top of the plate.Note that the deflection is expressed in a dimensionless form as 1000wD/16q 0 a 4 .For the plate of α 90 • , results in the literature 10 are also shown in the figure.In addition, Table 1 also shows the comparison results between the present solution, FEM, and literature at the center of the plate.The results show that our analytical and numerical solutions agree with each other very well.

Orthotropic Skewed Thick Plate under Concentrated Load
Orthotropic thick skewed plates are analyzed in this example, with the following material and geometrical properties: To show the applicability of this method for the arbitrary boundary condition, the following complex boundary condition is applied in this example: AB in Figure 1 is soft, simply supported condition SS1 , BC is clamped C , CD is free F , and DA is a hard, simply supported condition SS2 .As with the previous example, the convergence of the finite element model should be checked.Figure 8 shows the convergence study of transverse displacements w at the center of the plates when α is 30 • .It can be concluded that the finite element model converges adequately when the number of elements 1600.Similar results are confirmed for other α.
In addition, the numbers of terms including l and m in 3.1 and 3.5 need to be determined.Figure 9 shows the deflection w at the point where a concentrated load is applied for skew angle α 30 • and K 7, as one of the cases considered, for various l and m values.It is seen that the deflection at l, m 9, 115 is well converged.Similar results are observed for other cases.Therefore l, m 9, 115 is employed here and is also used as the reference for comparison.
For this example, because no previous work in the literature has been found reporting a similar experience, only FEM analysis results by ANSYS are employed for comparison with our analytical solution results.Table 2 shows comparison results of the deflection at the center of the plate X, Y 0, 0 .The deflection is expressed in a dimensionless form as 100wD 11 /4Q 0 a 2 .
Figures 10, 11, and 12 indicate the results for deflection along line HF in Figure 1.The results indicate that the response behavior for this case is much more complex than the examples above, due to nonsymmetric loading and boundary conditions.Due to the oblique coordinate system, the load at X, Y −a/2, b/2 has different relative relations to the interested responses on Y 0, along with different skew angles.For example, this causes the peak responses of the deflection to move towards the center of the plate with the skew angle decreasing from 90 to 30 degrees, and it makes the deflection of α 30 • larger than that of α 90 • .

Conclusions
The governing differential equation of skewed thick plates on an elastic foundation in an oblique coordinate system is formulated for the first time in this paper.This work allows derivation of an analytical solution for any boundary condition, loading condition, and material, also for the first time reported.All response quantities including shear forces, moments, stresses, strains, deflections, and rotation angles can be readily derived from the proposed potential function ψ.The two illustrative examples show that the analytical solutions are in good agreement with those reported in literature and numerical solutions by FEM.Further validation of the analytical solution developed here will be provided by an experimental test.
It is also worth noting that the approach to the governing differential equation and its analytical solution developed in this study can be used for further studies including, but not limited to, continuous plate analysis and dynamic analysis, for which only numerical solutions exist.These studies are in progress now and may contribute additional understanding of the behavior of a structure on an elastic foundation.

Figure 1 :
Figure 1: A skewed plate in an oblique coordinate system.

Figure 3 :
Figure 3: Convergence study of the center deflection of a clamped isotropic 30-degree skewed thick plate under uniform loading.

Figure 4 :
Figure 4: Truncation effect for convergence for center deflection of clamped isotropic 30-degree skewed thick plate under uniform loading.

Figure 5 :
Figure 5: Analytical and numerical solutions by ANSYS and from the literature 10 of the deflection of a clamped isotropic skewed thick plate under uniformly distributed load.Skew angle 90 • .

Figure 6 :
Figure 6: Analytical and numerical solutions by ANSYS of the deflection of a clamped isotropic skewed thick plate under uniformly distributed load.Skew angle 60 • .

Figure 7 :
Figure 7: Analytical and numerical solutions by ANSYS of the deflection of a clamped isotropic skewed thick plate under uniformly distributed load.Skew angle 30 • .1

Figure 8 :
Figure 8: Convergence study of the center deflection of an orthotropic 30-degree skewed thick plate under concentrated loading.

Figure 9 :
Figure 9: Truncation effect for convergence for maximum deflection of an orthotropic 30-degree skewed thick plate under concentrated loading.

Figure 10 :
Figure 10: Analytical and numerical solutions by ANSYS of the deflection of an orthotropic skewed thick plate under concentrated load.Skew angle 90 • .

Figure 11 :
Figure 11: Analytical and numerical solutions by ANSYS of the deflection of an orthotropic skewed thick plate under concentrated load.Skew angle 60 • .

Figure 12 :
Figure 12: Analytical and numerical solutions by ANSYS of the deflection of an orthotropic skewed thick plate under concentrated load.Skew angle 30 • .

Table 1 :
Deflection at the center of clamped skewed thick plates under uniform loading.The numbers in parentheses show the difference between these results and the analytical solution developed in this research.

Table 2 :
Deflection at the center of orthotropic skewed thick plates under concentrated loading.The numbers in parentheses show the difference between these results and the analytical solution developed in this research.
•, and 90 • are employed.The geometrical properties used in this research are a b, and t 0.2a.The clamped boundary condition is used for all edges.The modulus of foundation k is expressed in a dimensionless form as D is the bending stiffness and is expressed as D Et 3 /12 1 − ν 2 for the isotropic plate.K 3, 5, and 7 are employed here.
ν xy 0.2, a b, t 0.3a, where E x and E y are Young's modulus along the x and y directions, ν xy is the major Poisson's ratio, and G xy , G xz , and G yz are the shear modulus in the xy, xz, and yz planes.These values determine D r , D o , A r , and A o in 2.2 and 2.6 .The external transverse force is a concentrated force of Q 0 applied at X, Y −a/2, b/2 .Plates with skew angle α 30 • , 60 • , and 90 • and the modulus of foundation of K 3, 5, and 7 are analyzed here, where K is expressed in a dimensionless form as K 16a 4 k/D 11 1/4 , and D 11 is the bending stiffness which is expressed as D 11 E x t 3 /12 1 − ν xy ν yx .