Published: 30 December 2023

A third-order shear deformation plate bending formulation for thick plates: first principles derivation and applications

Charles Chinwuba Ike1
1Department of Civil Engineering, Enugu State University of Science and Technology, Agbani, 402004 Enugu State, Nigeria
Views 209
Reads 122
Downloads 231

Abstract

A third-order shear deformation plate bending formulation is presented in this study from the first principles. The derivation assumed a displacement field constructed using third-order polynomial function of the transverse (z) coordinate; and made to apriori satisfy the linear three-dimensional (3D) kinematics relations as well as the transverse shear stress free boundary conditions at the top and bottom plate surfaces. The formulation thus has no need for shear stress correction factors of the first-order shear deformation plate theories. The domain equations of equilibrium are obtained as a set of three coupled differential equations in terms of three unknown displacements. The system of coupled equations is solved for simply supported rectangular and square plates subjected to four cases of loading distributions: sinusoidal loading, uniformly distributed loading, linearly distributed loading and point load at the plate center. Navier’s double trigonometric series method is used to construct trial solutions for the three displacement functions such that the boundary conditions are satisfied identically. The integration problem is thus reduced to an algebraic problem and is solved for each considered loading. It is found that the present formulation gives exact results for the normal stresses σxx for sinusoidal and uniformly distributed loads. The study further showed that the results for deflection and stresses agreed with Krishna Murty’s higher order shear deformation plate theory results. The present formulation gave accurate results because of the inclusion of transverse normal strain effects in the formulation. The formulation gives a quadratic variation of the transverse shear stresses across the thickness in consonance with the theory of elasticity method.

1. Introduction

Plates are three-dimensional (3D) structural members characterized by in-plane dimensions and a transverse dimension that is usually smaller than the in-plane dimensions. They are widely used in civil, structural, mechanical, naval, marine engineering to carry static, dynamic or compression loadings. They may be made of homogeneous or non-homogeneous materials, isotropic, orthotropic, laminated or composite materials, and hence can be classified accordingly based on their material makeup. They are also classified based on their thickness relative to the in-plane dimensions as thin, moderately thick or thick plates.

Depending on the nature of applied loadings, plates may be subject to static flexure, dynamic flexure or buckling.

Generally, plate problems are three-dimensional (3D) problems of elasticity theory. However, a rigorous 3D elasticity formulations of plate behaviour entails complicated mathematical problems even for simple considerations of isotropic, homogeneous materials.

Problems of plates have therefore been expressed by simplifying the 3D elasticity formulations to one-dimensional (1D) and two-dimensional (2D) approximations; and this has been the main objective of plate research.

The classical thin plate theory (CPT) commonly called the Kirchhoff plate theory (KPT) assumed the Navier-Kirchhoff hypothesis [1-10].

– orthogonality of the cross-sectional planes to the middle surface prior to flexural deformation and after bending deformation.

– the invariance of thickness during bending deformation.

– the following are the notable advantages of the KPT [1-16]:

– the equation of equilibrium is linear and contains only one unknown function – deflection of the middle surface.

– the internal force resultants are expressible in terms of the transverse deflection function.

– it gives parabolic variation of shear stresses over the thickness in consonance with mechanics of structures.

The main disadvantage of the KPT is the inability to consider transverse shear deformations, thus limiting the scope of applicability to thin plates where transverse shear deformations are ignorable without significant errors [17-22]. Despite the disadvantages, KPT has been found to be satisfactory for thin plates and various methods for solving KPT are found in references [1-22] and [23-24].

Research efforts to develop improved formulations and postulations to consider shear deformation effects and thus extend the scope of plate theory to moderately thick and thick plates led to the derivations by Reissner [25], [26]; Mindlin [27]; Krishna Murty [28], Srinivas and Rao [29], Shimpi and Patel [30]; Sayyad and Shinde [31]; Ghugal and Gajbiye [32], Sayyad and Shinde [33]; Bathini and Reddy [34], [35]; Bathini et al. [36]; Eipakchi and Moshir [37]; Zagaripoor et al. [38], Raissi et al. [39] and Rodrigues et al. [40].

2. Review of previous works

Ike [20], [22] studied Mindlin’s first order shear deformable plates. Nwoji et al [21] obtained satisfactory solutions for the flexural analysis of simply supported rectangular Mindlin plates subjected to sinusoidal transverse load distribution using the Navier’s double trigonometric series method.

Ike [41] used Fourier series method to find the stresses and deflections in thick beams. Ike et al [42] used least squares method solve the stresses in rectangular plates under parabolic edge loads. Onah et al [43] derived stress function for solving elastostatic problems of thick circular plates.

Onyeka et al [44-46] and Onyeka and Okeke [47] used polynomial displacement function in an energy formulation to solve the flexural problem of thick plate with simply supported, free and clamped boundaries. Their work considered shear deformation and did not need shear correction factors. They used 3D kinematics and constitutive relations in formulating the energy functional and minimization procedure for the equilibrium equations. This solution for center deflections were in error by 2.9 %-3.7 % compared with the exact solution.

Onyeka and Mama [48] and Onyeka et al [49] have used a trigonometric displacement function in an energy functional minimization method to obtain satisfactory solutions respectively for the bending and stability of thick plates.

Rouzegar and Abdoli Sharifpoor [50] developed a finite element formulation based on two-variable refined plate theory (RPT) for the flexural analysis of isotropic and orthotropic plates. The RPT used is applicable to thin and thick plates and gives parabolic transverse shear stress variation across the plate thickness, thus obviating the need for shear correction as boundary conditions are satisfied. They used variational principles to obtain the equilibrium and weak form equations; and considered a 4-node rectangular plate element with six nodal degrees of freedom. They found satisfactory solutions using MATLAB software on the resulting algebraic problem.

Gajbhiye et al [51] used a 5th order shear and normal deformation theory to satisfactorily solve the eigenvalue problem of free vibration of simply supported thick isotropic square plates.

Gajbhiye et al [52] used a quasi-three-dimensional theory that considered shear and normal deformation to solve the bending problem of simply supported sandwich plates. They considered shear stress free boundary conditions, thus obviating the need for stress correction. They used virtual work principle for the domain equations, and Navier’s method for satisfactory solutions for sinusoidal and uniform transverse loadings.

Ghugal and Sayyad [53], [54] used a trigonometric shear deformation theory (TSDT) to obtain satisfactory solutions to the elasticity problem of thick laminated plates under transverse loading.

Ghugal and Gajbhiye [32] developed a 5th order shear deformation theory that considered transverse shear deformation and transverse normal strain effects, and applied it to the bending analysis of thick plates. They used virtual work principles to obtain field equations and boundary conditions and Navier’s series method to solve the resulting boundary value problem (BVP).

Zargaripoor et al [55] used exact wave propagation approach for the first time to obtain free vibration and buckling solutions for thick rectangular plates modelled using third-order shear deformation plate theory. They considered plates with opposite simply supported edges while the other edges may be clamped or simply supported. They derived the matrices of wave propagation and reflection for the plate problem and by superposition, obtained the characteristic equation, which was solved for the dimensionless frequencies and buckling loads for the different boundary conditions studied.

Kumar et al [56] have used radial basis function based meshfree methods for the analysis of thick plates using higher order shear deformation theory. Makvandi et al [57] studied the behaviour of moderately thick plates under compressive load. Civalek and Ulker [58] used the harmonic differential quadratic method for the flexural analysis of thin isotropic circular plates. Civalek [59] used the discrete singular convolution (DSC) method to solve bending problems of thick rectangular plates. Other seminal works with significant insight to the thick plate problem are found in references [60-83].

In this study a third-order shear deformation plate bending formulation is presented from first principles for the modelling and solution of thick plate bending problems under transverse loadings.

2.1. Novelty of the study

The novelty of the study is the first principle approach adopted in the paper for the formulation of the thick plate bending problem using a third order shear deformation plate theory that satisfies the transverse shear stress boundary conditions. The study presents a systematic study that uses equilibrium method to derive the governing differential equations of equilibrium. The Navier’s double trigonometric series method is used in a first principles, systematic way to derive solutions for:

(i) sinusoidal distribution of transverse load,

(ii) uniformly distributed transverse load,

(iii) linearly varying transverse load,

(iv) transverse point load acting at any point on the plate domain.

2.2. Theoretical framework

The thick plate considered which is shown in Fig. 1 has length a, width b, and thickness h, and is defined using the Cartesian coordinates by: 0xa, 0yb, -h/2zh/2.

Fig. 1Thick plate

Thick plate

Basic assumptions of the formulation

(i) The plate material is homogeneous, isotropic and linearly elastic.

(ii The plate is only subjected to transverse load.

(iii) The body forces are ignored, but can be incorporated into the formulation by adding them to the distributed transverse load.

(iv) The in-plane components of displacement in the x and y directions and the transverse component of displacement in the z direction are small in comparison with the plate thickness.

2.3. Displacement field

Ignoring the in-plane deformations, the displacement field components for flexural deformations are:

1a
ux,y,z=zαx+z3βx=zαx+z2βx,
1b
vx,y,z=zαy+z3βy=zαy+z2βy,
1c
wx,y,z=wx,y,z=0=w0x,y,

where, ux,y,z is the displacement field component in the x coordinate direction, v(x,y,z) is the displacement field component in the y coordinate direction, w0(x,y) is the transverse (z) component of the displacement of the middle surface (z= 0), αx, αy, βx, βy are the unknown displacement variables αx(x,y) and αy(x,y) are the rotations of the normals to the middle surface of the plate about the y and x axes respectively; while βx(x,y) and βy(x,y) are called the displacement warping functions.

Five unknown functions, namely αx(x,y), αy(x,y), βx(x,y), βy(x,y) and w0x,y are used to describe the displacement field components in the Reddy plate theory (RdPT). Two of these unknown functions (αx(x,y), αy(x,y)) are encountered in the Mindlin plate theory (MdPT) which is a first order shear deformation plate theory (FSDPT). The displacement warping functions βx(x,y), βy(x,y) are derivable/expressible in terms of αx(x,y) and αy(x,y) by the imposition of the transverse shear stress free boundary conditions at the top and bottom faces of the plate.

2.4. Strain field (Kinematics)

Assuming infinitesimally small displacements, the strain fields are obtained from the strain-displacement relations of the small displacement linear elasticity as follows:

2
εxx=ux, εyy=vy, εzz=wz,γxy=uy+vx, γyz=vz+wy, γxz=uz+wx,

where εxx, εyy and εzz are the normal strains in the x, y and z Cartesian coordinate directions respectively; γxy, γyz and γxz are the shear strains.

By substitution of the displacement field components – Eq. (1) – into Eq. (2) we have:

3a
εxx=xzαx+z3βx=zαxx+z3βxx,
3b
εyy=xzαy+z3βy=zαyy+z3βyy,
3c
εzz=w0z(x,y)=0,
3d
γxy=xzαy+z3βy+yzαx+z3βx,
3e
γxy=zαxy+αyx+z3βxy+βyx,
3f
γyz=w0y+zzαy+z3βy=w0y+αy+3z2βy,
3g
γxz=w0x+zzαx+z3βx=w0x+αx+3z2βx.

2.5. Stress fields

The stress fields are determined by using the stress-strain relations of isotropic homogeneous elasticity. Thus:

4a
σxx=E1-μ2εxx+μεyy,
4b
σyy=E1-μ2εyy+μεxx,
4c
τxy=Gγxy,
4d
τyz=Gγyz,
4e
τxz=Gγxz,

where:

5
G=E2(1+μ),

and E is the Young’s modulus of elasticity, G is the shear modulus, μ is the Poisson’s ratio.

Substituting the expressions for normal and shear strains into Eqs. (4a-4e), the stress fields are found as:

6a
σxx=E1-μ2zαxx+z3βxx+μE1-μ2zαyy+z3βyy,
6b
σyy=E1-μ2zαyy+z3βyy+μE1-μ2zαxx+z3βxx,
6c
τyz=Gαy+3z2βy+w0y,
6d
τxz=Gαx+3z2βx+w0x,
6e
τxy=Gzαxy+αyx+Gz3βxy+βyx.

2.6. Enforcement of boundary conditions

The thick plate bending problem considered in this study is subjected to a distributed transverse load of intensity qx,y on the top surface (z=h/2) while the bottom surface (z=-h/2) of the plate is free of load. The top and bottom surfaces of the plate (z=±h/2) are free from shear stresses. The shear stress free boundary conditions at the top and bottom surfaces of the plate can be expressed as follows:

7a
τyzx,y,z=±h2=0,
7b
τxzx,y,z=±h2=0.

Applying the shear stress free boundary conditions on Eq. (6c) we have:

8a
Gαy+3z2βy+w0yz=±h/2=0,
8b
Gαy+3±h22βy+w0y=0.

G0 hence:

8c
αy+3h24βy+w0y=0.

Solving for βy gives:

8d
βy=-43h2αy+w0y.

Similarly, applying the boundary conditions – Eq. (7b) on Eq. (6d) gives:

9a
τxzx,y,z=±h2=Gαx+3z2βx+w0xz=±h/2=0.

G0. Hence:

9b
αx+3h24βx+w0x=0.

Solving for βx gives:

9c
βx=-43h2αx+w0x.

The relationships between βx and αx and w0 and between βy and αy and w0, thus reduce the number of unknown displacement parameters in the formulation to three, namely αx(x,y), αy(x,y) and w0(x,y).

2.7. Displacement field that satisfy shear stress free boundary conditions

The displacement field components are thus simplified and expressed in terms of the three unknown displacement parameters as follows:

10a
ux,y,z=zαx-43zh2w0x+αx,
10b
vx,y,z=zαy-43zh2w0y+αy,
10c
wx,y,z=w0x,y.

The displacement field components are now simplified and expressed in terms of only three unknown displacement parameters (αx(x,y), αy(x,y) and w0(x,y)) as in the MdPT and other FSDPTs even though a third order displacement variation has been apriori assumed in the formulation.

2.8. Strain fields that satisfy shear stress free boundary conditions

The resulting strain fields that satisfy the shear stress free boundary conditions at the top and bottom plate surfaces are now found as follows:

11a
εxx=zαxx-43zh22w0x2+αxx,
11b
εyy=zαyy-43zh22w0y2+αyy,
11c
εzz=0,
11d
γxy=zαxy-43zh22w0xy+αxy+αyx-43zh22w0xy+αyx,
11e
γyz=1-4z2h2αy+w0y,
11f
γxz=1-4z2h2αx+w0x.

2.9. Stress fields that satisfy shear stress free boundary conditions

The stress fields that satisfy the shear stress free boundary conditions on z=±h/2 surfaces are:

12a
σxx=Ez1-μ2αxx-43zh22w0x2+αxx+μEz1-μ2αyy-43zh22w0y2+αyy,
12b
σyy=Ez1-μ2αyy-43zh22w0y2+αyy+μEz1-μ2αxx-43zh22w0x2+αxx,
12c
τxy=Gzαxy+αyx-G43z3h222w0xy+αxy+αyx,
12d
τyz=G1-4z2h2αy+w0y,
12e
τzx=G1-4z2h2αx+w0x.

Simplification gives the normal stress fields as:

13a
σxx=Ez1-μ2αxx+μαyy-43z3h2E1-μ22w0x2+αxx+μ2w0y2+μαyy,
13b
σyy=Ez1-μ2αyy+μαxx-43z3h2E1-μ22w0y2+αyy+μ2w0x2+μαxx.

2.10. Internal stress resultants

The internal stress resultants are given by the following integration problems over the plate thickness:

14a
Mxx=-h/2h/2σxxz dz,
14b
Myy=-h/2h/2σyyz dz,
14c
Mxy=-h/2h/2τxyz dz,
14d
Qx=-h/2h/2τxy dz,
14e
Qy=-h/2h/2τyz dz,
14f
Px=-h/2h/2σxxz3 dz,
14g
Py=-h/2h/2z3σyy dz,
14h
Pxy=-h/2h/2z3τxy dz,
14i
Rx=-h/2h/2z2τxz dz,
14j
Ry=-h/2h/2z2τyz dz,

where Mxx, Myy are bending moments; Mxy is twisting moment; Px, Py, Pxy, Rx and Ry are the higher-order internal stress resultants.

By substitution of the expressions for the stress fields, and integration over the thickness of the plate, explicit expressions are obtained for the internal stress resultants. These expressions relate the internal stress resultants to the unknown displacements of the formulation.

Thus,

15
Mxx=-h/2h/2Ez21-μ2αxx+μαyy-43E1-μ2z4h22w0x2+μ2w0y2+αxx+μαyydz.

Simplifying:

Mxx=αxx+μαyy-h/2h/2Ez21-μ2dz
-2w0x2+μ2w0y2+αxx+μαyy-h2h/243E1-μ2z4h2dz.

Let:

16
D=-h/2h/2Ez21-μ2 dz=Eh312(1-μ2).

D is the modulus of flexural rigidity:

17
Mxx=Dαxx+μαyy-D52w0x2+μ2w0y2+αxx+μαyy,
17a
Mxx=Dαxx+μαyy-D52w0x2+μ2w0y2-D5αxx+μαyy,
17b
Mxx=4D5αxx+μαyy-D52w0x2+μ2w0y2.

Similarly:

18a
Myy=-h/2h/2Ez21-μ2αyy+μαxx-43E1-μ2z4h22w0y2+μ2w0x2+αyy+μαxxdz,
18b
Myy=4D5αyy+μαxx-D52w0y2+μ2w0x2,
19a
Mxy=-h/2h/2Gz2αxy+αyx--h/2h/243h2Gz4dz22w0xy+αxy+αyx.

Simplifying:

19b
Mxy=Gh315αxy+αyx-Gh3302w0xy,
20
Qx=-h/2h/2G1-4z2h2αx+w0xdz=2Gh3αx+w0x,
21
Qy=-h/2h/2G1-4z2h2αy+w0ydz=2Gh3αy+w0y,
22a
Pxx=-h/2h/2Ez41-μ2αxx+μαyydz--h/2h/243h2Ez61-μ2dz2w0x2+μ2w0y2+αxx+μαyy,
22b
Pxx=4Dh235αxx+μαyy-Dh2282w0x2+μ2w0y2,
23a
Pyy=-h/2h/2Ez41-μ2αyy+μαxxdz--h2h243h2Ez61-μ2dz2w0y2+μ2w0x2+αyy+μαxx,
23b
Pyy=4Dh235αyy+μαxx-Dh2282w0y2+μ2w0x2,
24a
Pxy=-h/2h/2Gz4dzαxy+αyx--h2h243h2Gz6dz22w0xy+αxy+αyx,
24b
Pxy=Gh530αxy+αyx-Gh51682w0xy,
25
Rx=-h/2h/2Gαx+w0x z21-4z2h2dz=Gh330αx+w0x,
26
Ry=-h/2h/2Gαy+w0y z21-4z2h2dz=Gh330αy+w0y.

3. Governing equations of equilibrium

The principle of minimum potential energy gives the total potential energy functional βx and Π as follows:

27
=z=(-h/2)z=h/2R212(σxxεxx+σyyεyy+τxyγxy+τyzγyz+τxzγxz)dxdydz-R2qw(x,y)dxdy,

where R2 is the two-dimensional (2D) domain of the plate (0xa,0yb), a and b are the in-plane dimensions of the plate, q(x,y) is the distribution of transverse load over the plate domain.

For equilibrium:

28
δ=0.

The equilibrium of the internal stress resultants or the principle of minimum total potential energy are used to obtain the differential equations of equilibrium of the plate. Using the equilibrium approach, the equations of equilibrium are found using:

29a
Nxxx+Nxyy=0,
29b
Nxyx+Nyyy=0,
29c
Qxx+Qyy+43h22Pxxx2+z2Pxyxy+2Pyyy2+q(x,y)=0,
29d
Mxxx+Myyy-Qx=0,
29e
Mxyx+Myyy-Qy=0,

where Nxx, Nyy, Nxy are resultant in-plane forces.

The governing differential equations of equilibrium of the shear deformable thick plate are found in terms of the three displacement parameters as follows:

30
3w0x3+3w0xy2-1742αxx2+1-μ22αxy2+1+μ22αyxy
+211-μh2αx+w0x=0,
31
3w0y3+3w0yx2-1742αyy2+1-μ22αyx2+1+μ22αxxy
+211-μh2αy+w0y=0,
32
4w0-1653αxx3+3αxxy2+3αyx2y+3αyy3-3365h2(1-μ)2w0
-3365h2(1-μ)αxx+αyy-21qx,yD=0,

where:

33
2=2x2+2y2,
34
4=22=(2)2=2x2+2y22,
35
4=4x4+24x2y2+4y4.

Alternatively, the domain equations are:

36
x(2w0)-1742αxx2+1-μ22αxy2+1+μ22αyxy+211-μh2αx+w0x=0,
37
y(2w0)-1742αyy2+1-μ22αyx2+1+μ22αxxy+211-μh2αy+w0y=0,
38
4w0(x,y)-165x(2αx)+y(2αy)-3365h2(1-μ)2w0 -3365h2(1-μ)αxx+αyy-21qx,yD=0.

3.1. Alternative exact method for the study by Srivinas and Rao [29]

Srivinas and Rao [29] presented an alternative method for formulating and solving the candidate problem of this study. They presented a unified exact analysis for the static flexural solutions of simply supported thick plates and laminated plates based on a three-dimensional linear, small deformation theory of elasticity. Their work considered orthotropic plates, but isotropic solutions were also presented as special cases of the orthotropic material. They obtained formally exact solutions which were presented as simple infinite series for stresses and displacements for static flexural solutions.

3.2. Advantages and limitations of the third order shear deformable plate formulation

These formulated third order shear deformable plate equations avoid some issues associated with the use of lower order and simpler Kirchoff plate theory (KPT) and first order shear deformable plate theories (FSDPTs) of Timoshenko and Mindlin.

In particular, the advantages are:

1) The formulation allows the determination of a more accurate stress analysis without the use of shear correction factors of the FSDPTs of Timoshenko and Mindlin.

2) The resulting displacement field for the present theory determines a quadratic variation of shear strains and shear stresses across the thickness, with the shear strains and shear stresses vanishing at the top and bottom surfaces of the plate and attaining the maximum values at the middle plane. This leads to the satisfaction of the shear stress-free boundary conditions at the top and bottom plate surfaces z=±h/2, and the obviation of the need for shear correction factors.

3) The governing equations formulated take due consideration of transverse shear deformation in deriving the equations of equilibrium, and hence make the resulting formulation applicable to moderately thick and thick plates where transverse shear deformation is vital to their behaviour.

The disadvantages are:

1) The formulation assumes that the plate material is homogenous and hence the equations would not apply to non-homogenous plates.

2) The assumption of isotropic plate material also restricts the applicability to isotropic plates.

3) The formulation applies to linear, small displacement static flexure problems and is inapplicable to large displacement, nonlinear flexural problems.

4) The formulation contains three unknown displacements, which increase the difficulty in obtaining solutions to the governing equations as the solutions involve prior determination of the three unknown displacement parameters.

4. Application to the flexural analysis of simply supported thick plate

A rectangular thick plate with simply supported edges x= 0, x=a, y= 0, y=b is considered. The plate has in-plane dimensions a×b and thickness h, and is subjected to a distributed transverse load of intensity q(x,y) over the domain 0xa,0yb.

The boundary conditions along the simply supported edges are given by along x=0:

39a
w(x=0,y)=0,
39b
Pxx(x=0,y)=0,
39c
Mxx(x=0,y)=0,
39d
αy(x=0,y)=0.

Similarly, along x=a:

40a
w(x=a,y)=0,
40b
Pxx(x=a,y)=0,
40c
Mxx(x=a,y)=0,
40d
αy(x=a,y)=0.

Along y=0:

41a
w(x,y=0)=0,
41b
Pyy(x,y=0)=0,
41c
Myy(x,y=0)=0,
41d
αx(x,y=0)=0.

Along y=b:

42a
w(x,y=b)=0,
42b
Pyy(x,y=b)=0,
42c
Myy(x,y=b)=0,
42d
αx(x,y=b)=0.

4.1. Navier’s double trigonometric series method of solution

Using Navier’s double trigonometric series method and considering the displacement shape functions that satisfy apriori the boundary conditions along the four edges, w(x,y), αx(x,y) and αy(x,y) are expressed in terms of double series of infinite terms as:

43
wx,y=m=1n=1Amnsinλmxsinγny,
44
αxx,y=m=1n=1Bmncosλmxsinγny,
45
αyx,y=m=1n=1Cmnsinλmxcosγny,

where Amn, Bmn and Cmn are the unknown parameters of w(x,y), αx(x,y), and αy(x,y) respectively, and λm and γn are defined in terms of m, n, a and b as:

46
λm=mπa,
47
γn=nπb.

The distributed transverse load q(x,y) is represented as the double Fourier sine series expansion:

48
qx,y=m=1n=1qmnsinλmxsinγny,

where qmn are the Fourier series coefficients of the load.

Then by Fourier series theory:

49
qmn=4ab0a0bq(x,y)sinλmxsinγny dxdy.

Substitution of Eqs. (43-45) and Eq. (49) into the governing PDEs and simplification gives the following system of equations:

50
-λm(λm2+γn2)+21(1-μ)h2λmAmn+174λm2+1-μ2γn2+21(1-μ)h2Bmn +178(1+μ)λmγnCmn=0,
51
γn21(1-μ)h2-(λm2+γn2)Amn+178(1+μ)λmγn Bmn
+174γn2+1-μ2λm2+21(1-μ)h2Cmn=0,
52
(λm2+γn2)2-336(1-μ)5h2(λm2+γn2)Amn+336(1-μ)5h2λm-165λm(λm2+γn2)Bmn
+336(1-μ)5h2γn-165γn(λm2+γn2)Cmn=21Dqmn.

4.2. Matrix representation of the algebraic problem

The resulting algebraic problem is represented in matrix format as follows:

53
a11a12a13a21a22a23a31a32a33AmnBmnCmn=00Qmn,

wherein:

54
a11=21(1-μ)h2λm-λmλm2+γn2,
a12=174λm2+1-μ2γn2+21(1-μ)h2,
a13=1781+μλmγn,
a21=γn21(1-μ)h2-(λm2+γn2),
a22=1781+μλmγn=a13,
a23=174γn2+1-μ2λm2+21(1-μ)h2,
a31=(λm2+γn2)2+3361-μ5h2λm2+γn2,
α32=336(1-μ)5h2λm2-165λmλm2+γn2,
α33=336(1-μ)5h2γn-165γnλm2+γn2,
Qmn=21qmnD.

Case 1: Sinusoidal distribution of transverse load:

55
qx,y=q0sinπxasinπyb.

Then:

56
qmn=4ab0a0bq0sinπxasinπybsinmπxasinnπyb dxdy,
57
qmn=q0,m=n=1,0,m>1, n>1.

Case 2: Uniformly distributed transverse load q0 on the surface z=-h/2 acting in the z direction:

58
qx,y=q0, 0xa, 0yb,
59a
qmax=4ab0a0bq0sinmπxasinnπyb dxdy=16q0mnπ2, m=1,3,5,7,9, n=1,3,5,7,9,
59b
qmn=0, m=2,4,6,8,10, n=2,4,6,8,10.

Case 3: Linearly varying transverse load:

60
qx,y=q0xa.

Then:

61a
qmn=4ab0a0bq0xasinmπxasinnπyb dxdy,
61b
qmn=8q0mnπ2cosmπ.

Case 4: Transverse point load P0 acting at x=x-, y=y-, then:

62
qx,y=P0δx=x-, y=y-,
63
qmn=4ab0a0bP0δ(x=x-,y=y-)sinmπxasinnπyb dxdy,
64
qmn=4P0absinmπx-asinnπy-b.

For center point loads, x-=a2, y-=b2:

65
qmn=4P0absinmπ2sinnπ2.

5. Results and discussion

For each loading case considered, non-dimensional transverse displacement w- and the normal stresses (σ-xx,σ-yy) and shear stresses (τ-xy,τ-zx,τ-yz) are calculated at defined points on the plate as follows:

66
w-=100Ewqhs, σ-xx=σxxqs2, σ-yy=σyyqs2, τxy=τxyqs2, τ-zx=τzxqs, τyz=τyzqs,

where s=a/h.

The results for the dimensionless displacements and stresses are shown in Tables 1-4 for the four cases of load distributions considered in the study. In the tables, τ-zx and τ-yz were evaluated using two procedures namely constitutive relations, which is denoted by CR and equilibrium equations denoted by EE as superscripts for each of the stress notations. The percentage error in the results obtained and other results from literature sources are calculated by comparison with the exact solution as follows:

67
% error=value of result-exact resultexact result×100 %.

The values enclosed in brackets in Tables 1-4 are the % errors for the displacements and stresses relative to the exact solutions computed by Srivinas and Rao [29] as indicated in the corresponding table.

Table 1 presents the results for dimensionless displacement w- and stresses for the case of sinusoidal load distribution on a square thick plate for h/a=0.10 and for h/a=0.25. Table 1 shows that the present results for w- at the center (for h/a=0.10) is 0.62 % greater than the exact result obtained by Srivinas and Rao [29]. Table 1 further shows that Krishna Murty’s HSDPT solutions gave a better result for w- at the center as the result coincides with the exact result. However, the present result gave better result than the Kirchhoff CPT result which showed a percentage error of –4.76 %.

Table 1Dimensionless deflection w- at (x=a/2, y=b/2, z=0), Dimensionless in-plane normal stress σ-xx at (x=a/2, y=b/2, z=h/2), σ-yy at (x=a/2, y=b/2, z=h/2) in-plane shear stress τ-xy at (x=0, y=0, z=h/2) and transverse shear stress τ-zx at (x=0, y=bx, z=0) in square thick plate subjected to single sine load

ab
ha
Theory
w-
σ-xx
σ-yy
τ-xy
τ-zxCR
τ-zxEE
1
0.1
Present work
2.960
(0.62 %)
0.289
(0.00 %)
0.2890
(0.00 %)
0.1070
0.2380
(0.00 %)
0.228
Krishna Murty [28] HSDT
2.942
(0.00%)
0.290
(0.34 %)
0.290
(0.34 %)
0.107
0.238
(0.00 %)
0.228
Kirchhoff [6] (CPT)
2.802
(–4.76 %)
0.287
(–1.50%)
0.287
(–1.50 %)
0.106
0.238
(Exact) Srivinas et al [29]
2.942
0.289
0.289
0.238
0.25
Present Study
3.7810
(3.39%)
0.2090
(2.45 %)
0.2090
(2.45%)
0.1120
0.2370
(0.38 %)
Sayyad and Ghugal [31]
3.7480
(2.32 %)
0.2130
(4.41%)
0.2130
(4.41 %)
0.1140
0.2380
(0.80 %)
Kirchhoff [6]
2.8030
(–23.48 %)
0.1970
(–3.43 %)
0.1970
(–3.43 %)
0.1060
(Exact) Srivinas and Rao [29]
3.6630
0.2040
0.2040
0.2361

Table 2Non-dimensionless transverse deflection w- at (x=a/2, y=b/2, z=0), Dimensionless in-plane stresses σ-xx at (x=a/2, y=b/2, z=h/2),σ-yy at (x=a/2, y=b/2, z=h/2) in-plane shear stress τ-xy at (x=0, y=0, z=h/2) and transverse shear stress τ-zx at (x=0,y=b/2,z=0) in rectangular thick plate subjected to uniformly distributed loading

ab
ha
Theory / Method / Procedure
w-
σ-xx
σ-yy
τ-xy
τ-zxCR
τ-zxEE
0.5
0.10
Present Study
11.420
(0.39 %)
0.612
(0.00 %)
0.278
(–1.06 %)
0.280
0.679
(0.00 %)
0.6776
(–0.206 %)
Krishna Murty [28] HSDT
11.310
(–0.57 %)
0.613
(0.163 %)
0.310
(10.32 %)
0.278
0.682
(0.441 %)
0.667
(–1.76 %)
Kirchhoff [6] (CPT)
11.060
(–2.78 %)
0.610
(–0.32 %)
0.278
(–1.06 %)
0.277
0.686
(1.03 %)
0.6865
(1.104 %)
(Exact) Srivinas and Rao [29]
11.375
0.612
0.281
0.679

Table 1 also reveals that the present results for σ-xx and σ-yy (at x=a/2, y=b/2, z=h/2) coincide with the exact result obtained by Srivinas et al [29] while Krishna Murty’s results for σ-xx and σ-yy are 0.34 % in error and CPT of Kirchhoff is –1.50 % in error.

For the case of sinusoidal load distribution for h/a=0.25, the present result for w- at the center is 3.39 % different from the exact result given by Srivinas et al [29], while the exponential shear deformation plate theory solution by Sayyad and Ghugal [31] gave better result for w- at the center as the error is 2.32 %. The error of the CPT is however –23.48 % for w- at the center for h/a=0.25.

The present method gave more accurate results as the error is 2.45 % compared with 4.41 % error of the ESDPT of Sayyad and Ghugal and –3.43 % error of the CPT.

In the both cases of h/a=0.1 and h/a=0.25, the τ-xywas calculated by the present method, but was not determined by the exact results of Srivinas et al [29]. In both cases, the present study gave satisfactory results for τ-zx.

Table 2 presents non-dimensional values of w- and stresses for rectangular thick plate (with a/b= 0.5, h/a=0.10) subjected to uniformly distributed loading. Table 2 shows that the present results for w- at the center has a relative difference of 0.39 % as compared with the exact solution by Srivinas et al [29]. The present result for w- is more accurate than the result by Krishna Murty’s HSDPT which has an error of –0.57 % and the CPT which has an error of –2.78 %. Table 2 further shows that the present results for σ-xx is exact as it coincides with the exact solution by Srivinas et al [29], while Krishna Murty’s result gave an error of 0.163% and CPT gave an error of –0.32 %.

The present result for σ-yy at the point (x=a/2, y=b/2, z=h/2) gave an error of –1.06 % compared with the exact result by Srivinas et al [29]. The present result is more accurate than Krishna Murty’s HSDPT result which has an error of 10.32 %. The present result for τ-zxCR gave an exact result compared with Krishna Murty’s HSDPT result which has an error of 0.441 % and CPT with an error of 1.03 %.

Table 3Dimensionless deflection w- at (x=a/2, y=b/2, z=0), in-plane normal stress σ-xx at (x=a/2, y=b/2, z=h/2), σ-yy at (x=a/2, y=b/2, z=h/2) in-plane shear stress τ-xy at (x=0, y=0, z=h/2) and transverse shear stress τ-zx at (x=0,y=b/2,z=0) in square and rectangular thick isotropic plate subjected to linearly distributed loading

ab
ha
Theory / Method /
Reference
w-
σ-xx
σ-yy
τ-xy
τ-zxCR
τ-zxEE
0.5
0.1
Present work
5.7100
(0.395 %)
0.3060
(0.00 %)
0.139
(–1.068 %)
0.1400
0.3395
(0.00 %)
0.3388
(–0.206 %)
Krishna Murty
[28] HSDT
5.6878
(0.0017 %)
0.3067
(0.228 %)
0.155
(10.32 %)
0.1390
0.341
(0.441 %)
0.3375
(–0.589 %)
Kirchhoff [6] (CPT)
5.5300
(–2.772 %)
0.3048
(–0.39 %)
0.139
(–1.068 %)
0.1385
0.343
(1.03 %)
0.3433
(1.104 %)
(Exact) Srivinas and Rao [29]
5.6875
0.3060
0.1405
0.3395
1
0.1
Present Study
2.3325
(0.560 %)
0.1445
(0.00 %)
0.1445
(0.00 %)
0.0995
0.2463
(0.984 %)
0.243
(–0.20 %)
Krishna Murty [28]
HSDT
2.3199
(0.017 %)
0.1453
(0.553 %)
0.1453
(0.553 %)
0.0975
0.2454
(0.615 %)
0.240
(–1.44 %)
Kirchhoff [6]
2.2180
(–4.375 %)
0.1435
(–0.007 %)
0.1435
(–0.007 %)
0.0975
0.247
(1.64 %)
(Exact) Srivinas
and Rao [29]
2.3195
0.1445
0.1445
0.2439

Table 3 presents the results for w- and stresses for square and rectangular thick plates subjected to linearly distributed loading. For rectangular thick plate h/a=0.1, a/b=0.5, the present study gave w- that is 0.395 % in error compared with the exact result by Srivinas et al [29]. Krishna Murty’s HSDPT result with an error of 0.0017 % gave better accuracy while the present study is more accurate than the CPT which has –2.772 % error. Table 3 also shows that the present method gave accurate results for σ-xx and σ-yy than the Krishna Murty’s HSDPT and Kirchhoff’s CPT results. For σ-xx the present result is identical with the exact result while Krishna Murty had an error of 0.228 % and CPT an error of –0.39 %. The present method gave exact results for τ-zxCR while Krishna Murty’s HSDPT result had an error of 0.441 % and CPT had an error of 1.03 %.

For square thick plate presented in Table 3, the present study gave an error of 0.56 % for w- at the center compared with 0.017 % error of the Krishna Murty HSDPT and –4.375 % error of the CPT. Table 3 also shows that the present results for σ-xx and σ-yy are exact compared with the Srivinas et al [29] results; and are more accurate than Krishna Murty’s and CPT results.

Table 4 presents the results for w- and stresses for thick square plates for h/a=0.10 subjected to center point load. Table 4 shows that the present study gave an error of 2.701 % for w- at the center which is more accurate than the CPT results with an error of –10.173 %. Moreover, Table 4 shows that Krishna Murty’s HSDPT gave more accurate result for w- at the center as the error is –0.063 %. However, for σ-xx, the present results gave an error of –2.481 % showing better accuracy than the Krishna Murty’s result (which has –6.609 % error) and CPT results (which has –27.03 % error).

Table 4Dimensionless transverse deflection w- at (x=a/2, y=b/2, z=0), normal stress σ-xx at (x=a/2, y=b/2, z=h/2), and transverse shear stress τ-zx at (x=0,y=b/2,z=0) in square thick isotropic plate subjected to center point load

ab
ha
Method / Reference
w-
σ-xx
τ-zxCR
1.5
0.10
Present study
14.4717
(2.701 %)
2.4956
(–2.481 %)
0.9174
(24.66 %)
Krishna Murty [28] HSDT
14.0821
(–0.063 %)
2.3902
(–6.609 %)
0.7936
(7.840 %)
Kirchhoff [6] (CPT)
12.6575
(–10.173 %)
1.8672
(–27.03 %)
(Exact) Srivinas and Rao [29]
14.0910
2.5591
0.7359

6. Conclusions

The study has presented from first principles, a third-order shear deformation plate bending formulation for thick plates.

1) The formulation is displacement based and includes the effect of transverse shear deformations thus making it applicable to thick plates.

2) The formulation is based on infinitesimally small displacement assumptions, and linear 3D elasticity kinematics are used.

3) The resulting formulation satisfies the transverse shear stress free boundary conditions at the top and bottom surfaces of the plate, and hence transverse shear stress correction factors of the Mindlin first-order shear deformation plate theories are not needed in the present formulation.

4) The resulting governing differential equations of equilibrium are a set of three coupled differential equations in terms of three unknown displacements, namely w0(x,y), αx(x,y) and αy(x,y).

5) The resulting formulation is solved for simply supported rectangular and square thick plates subjected to four cases of loading distributions:

– sinusoidal loading

– uniformly distributed loading

– linearly distributed loading

– point load at the plate center.

6) For simply supported boundaries, Navier’s double trigonometric series method is used to construct trial solutions for w(x,y), αx(x,y) and αy(x,y) that apriori satisfy all the boundary conditions. The trial solutions are expressed in terms of unknown displacement parameters Amn, Bmn, Cmn.

7) Navier’s double trigonometric series method reduces the integration problem over the 2D domain of the plate to an algebraic problem where the unknown displacement parameters are Amn, Bmn, Cmn for each load case studied.

8) Unlike the Kirchhoff plate theory (KPT) the presented formulation can be used for both thin and thick plates and yields parabolic (quadratic) variation of transverse shear stresses across the plate thickness, in agreement with theory of elasticity methods.

9) The presented formulation gives exact results for σxx (at x=a/2, y=b/2, z=h/2) for sinusoidal loading (for b/a=0.1, a/b=1); uniformly distributed and linearly distributed loads (for h/a=0.1 and a/b=0.5).

10) The displacements and stresses obtained by the present formulations are in good agreement with the Krishna Murty’s HSDPT results.

11) The present formulation gave accurate results for transverse deflection because of the inclusion of transverse normal strain effects in the formulation.

References

  • C. Ike, “Generalized integral transform method for the bending analysis of clamped rectangular thin plates,” Journal of Computational Applied Mechanics, Vol. 53, No. 4, pp. 599–625, Dec. 2022, https://doi.org/10.22059/jcamech.2022.350620.768
  • C. C. Ike, “Variational Ritz-Kantorovich-Euler Lagrange method for the elastic buckling analysis of fully clamped Kirchhoff thin plate,” ARPN Journal of Engineering and Applied Sciences, Vol. 16, No. 2, pp. 224–241, 2021.
  • C. C. Ike, “Double Fourier cosine series method for the flexural analysis of Kirchhoff plates on Winkler foundation,” Journal of Geotechnical and Transportation Engineering, Vol. 4, No. 2, pp. 30–38, 2018.
  • C. C. Ike, “Kantorovich-Euler Lagrange-Galerkin method for bending analysis of thin plates,” Nigerian Journal of Technology, Vol. 36, No. 2, pp. 351–360, 2017, https://doi.org/10.4314/nijtv36i2.5
  • C. Chinwuba Ike, “Flexural analysis of rectangular Kirchhoff plate on Winkler foundation using Galerkin-Vlasov variational method,” Mathematical Modelling of Engineering Problems, Vol. 5, No. 2, pp. 83–92, Jun. 2018, https://doi.org/10.18280/mmep.050205
  • G. Kirchhoff, “Über das Gleichgewicht und die Bewegung einer elastischen Scheibe,” Journal für die reine und angewandte Mathematik, Vol. 40, pp. 51–88, 1850.
  • S. Ullah, J. Zhang, and Y. Zhong, “New analytical solutions of buckling problem of rotationally-restrained rectangular thin plates,” International Journal of Applied Mechanics, Vol. 11, No. 10, p. 1950101, Dec. 2019, https://doi.org/10.1142/s1758825119501011
  • J. Zhang, S. Ullah, and Y. Zhong, “Accurate free vibration solutions of orthotropic rectangular thin plates by straightforward finite integral transform method,” Archive of Applied Mechanics, Vol. 90, No. 2, pp. 353–368, Feb. 2020, https://doi.org/10.1007/s00419-019-01613-1
  • V. Singh and M. H. Prashanth, “Deflection surface analysis of thin plate structures using regression technique,” in International Conference on Structural Engineering and Construction Management, 2022, https://doi.org/10.1007/978-3-031-12011-4-18
  • M. Delyavskyy and K. Rosiński, “The new approach to analysis of thin isotropic symmetrical plates,” Applied Sciences, Vol. 10, No. 17, p. 5931, Aug. 2020, https://doi.org/10.3390/app10175931
  • C. Nwoji, H. Onah, B. Mama, and C. Ike, “Ritz variational method for bending of rectangular Kirchhoff plate under transverse hydrostatic load distribution,” Mathematical Modelling of Engineering Problems, Vol. 5, No. 1, pp. 1–10, Mar. 2018, https://doi.org/10.18280/mmep.050101
  • C. C. Ike, M. E. Onyia, and E. O. Rowland-Lato, “Generalized integral transform method for bending and buckling analysis of rectangular thin plate with two opposite edges simply supported and other edges clamped,” Advances in Science, Technology and Engineering Systems Journal, Vol. 6, No. 1, pp. 283–296, Jan. 2021, https://doi.org/10.25046/aj060133
  • H. Onah, C. Nwoji, C. Ike, and B. Mama, “Elastic buckling analysis of uniaxially compressed CCSS Kirchhoff plate using single finite Fourier sine integral transform method,” Modelling, Measurement and Control B, Vol. 87, No. 2, pp. 107–111, Jun. 2018, https://doi.org/10.18280/mmc_b.870208
  • O. A. Oguaghamba and C. C. Ike, “Single finite Fourier sine integral transform method for the determination of natural frequencies of flexural vibration of Kirchhoff plates,” International Journal of Engineering Research and Technology, Vol. 13, No. 3, pp. 470–476, 2020.
  • O. A. Oguaghamba and C. C. Ike, “Galerkin-Vlasov method for the elastic buckling analysis of Kirchhoff plate with one free edge and three simply supported edges under uniform uniaxial compression,” ARPN Journal of Engineering and Applied Sciences, Vol. 15, No. 14, pp. 1574–1581, 2020.
  • M. E. Onyia, E. O. Rowland-Lato, and C. C. Ike, “Galerkin-Kantorovich method for the elastic buckling analysis of thin rectangular SCSC plates,” International Journal of Engineering Research and Technology, Vol. 13, No. 4, pp. 613–619, 2020.
  • M. E. Onyia, E. O. Rowland-Lato, and C. C. Ike, “Elastic buckling analysis of SSCF and SSSS rectangular thin plates using the single finite Fourier sine integral transform method,” International Journal of Engineering Research and Technology, Vol. 13, No. 6, pp. 1147–1158, 2020.
  • M. E. Onyia, E. O. Rowland-Lato, and C. C. Ike, “Galerkin-Vlasov variational method for the elastic buckling analysis of SSCF and SSSS rectangular plates,” International Journal of Engineering Research and Technology, Vol. 13, No. 6, pp. 1137–1146, 2020.
  • B. O. Mama, O. A. Oguaghamba, and C. C. Ike, “Single finite Fourier sine integral transform method for the flexural analysis of rectangular Kirchhoff plate with opposite edges simply supported, other edges clamped for the case of triangular load distribution,” International Journal of Engineering Research and Technology, Vol. 13, No. 7, pp. 1802–1813, 2020.
  • C. Chinwuba Ike, “Mathematical solutions for the flexural analysis of Mindlin’s first order shear deformable circular plates,” Mathematical Models in Engineering, Vol. 4, No. 2, pp. 50–72, Jun. 2018, https://doi.org/10.21595/mme.2018.19825
  • C. U. Nwoji, B. O. Mama, H. N. Onah, and C. C. Ike, “Flexural analysis of simply supported rectangular Mindlin plates under bi-sinusoidal transverse load,” ARPN Journal of Engineering and Applied Sciences, Vol. 13, No. 15, pp. 4480–4488, 2018.
  • C. C. Ike, “Equilibrium approach in the derivation of differential equations for homogeneous isotropic Mindlin plates,” Nigerian Journal of Technology, Vol. 36, No. 2, pp. 346–350, 2017, https://doi.org/10.4314/nijt.v36i2.4
  • Civalek, “Application of differential quadrature (DQ) and harmonic differential quadrature (HDQ) for buckling analysis of thin isotropic plates and elastic columns,” Engineering Structures, Vol. 26, No. 2, pp. 171–186, Jan. 2004, https://doi.org/10.1016/j.engstruct.2003.09.005
  • G. O. Alcybeev, D. P. Goloskokov, and A. V. Matrosov, “The superposition method in the problem of bending of a thin isotropic plate clamped along the contour,” Vestnik of Saint Petersburg University. Applied Mathematics. Computer Science. Control Processes, Vol. 18, No. 3, pp. 347–364, 2022, https://doi.org/10.21638/11701/spbu10.2022.305
  • E. Reisnner, “On the theory of bending of elastic plates,” Journal of Mathematics and Physics, Vol. 23, pp. 184–191, 1944.
  • E. Reissner, “The effect of transverse shear deformation on the bending of elastic plates,” Journal of Applied Mechanics, Vol. 12, No. 2, pp. A69–A77, Jun. 1945, https://doi.org/10.1115/1.4009435
  • R. D. Mindlin, “Influence of rotatory inertia and shear on flexural motions of isotropic, elastic plates,” Journal of Applied Mechanics, Vol. 18, No. 1, pp. 31–38, Mar. 1951, https://doi.org/10.1115/1.4010217
  • Murty and A. V. Krishna, “Toward a consistent plate theory,” AIAA Journal, Vol. 24, No. 6, pp. 1047–1048, Jun. 1986, https://doi.org/10.2514/3.9388
  • S. Srinivas and A. K. Rao, “Bending, vibration and buckling of simply supported thick orthotropic rectangular plates and laminates,” International Journal of Solids and Structures, Vol. 6, No. 11, pp. 1463–1481, Nov. 1970, https://doi.org/10.1016/0020-7683(70)90076-4
  • R. P. Shimpi and H. G. Patel, “A two variable refined plate theory for orthotropic plate analysis,” International Journal of Solids and Structures, Vol. 43, No. 22-23, pp. 6783–6799, Nov. 2006, https://doi.org/10.1016/j.ijsolstr.2006.02.007
  • A. S. Sayyad and Y. M. Ghugal, “Bending and free vibration analysis of thick isotropic plates by using exponential shear deformation theory,” Applied and Computational Mechanics, Vol. 6, pp. 65–82, 2012.
  • Y. M. Ghugal and P. D. Gajbhiye, “Bending analysis of thick isotropic plates by using 5th order shear deformation theory,” Journal of Applied and Computational Mechanics, Vol. 2, No. 2, pp. 80–95, May 2016, https://doi.org/10.22055/jacm.2016.12366
  • A. Sayyad and B. Shinde, “A new higher-order theory for the static and dynamic responses of sandwich FG plates,” Journal of Computational Applied Mechanics, Vol. 52, No. Online First, pp. 102–125, Nov. 2020, https://doi.org/10.22059/jcamech.2020.313152.569
  • S. Bathini and K. Vijaya Kumar Reddy, “Flexural behavior of porous functionally graded plates using a novel higher order theory,” Journal of Computational Applied Mechanics, Vol. 51, No. Online First, pp. 361–373, Mar. 2020, https://doi.org/10.22059/jcamech.2020.298540.488
  • S. Bathini and V. K. R. K., “A refined inverse hyperbolic shear deformation theory for bending analysis of functionally graded porous plates,” Journal of Computational Applied Mechanics, Vol. 51, No. Online First, pp. 417–431, Sep. 2020, https://doi.org/10.22059/jcamech.2020.308469.544
  • S. Bathini, V. K. R. K., and C. A. B., “Free vibration behavior of bi-directional functionally graded plates with porosities using a refined first order shear deformation theory,” Journal of Computational Applied Mechanics, Vol. 51, No. Online First, pp. 374–388, Jul. 2020, https://doi.org/10.22059/jcamech.2020.303046.510
  • H. Eipakchi and S. Khadem Moshir, “Dynamic response determination of viscoelastic annular plates using FSDT – perturbation approach,” Journal of Computational Applied Mechanics, Vol. 51, No. 1, pp. 98–106, Jun. 2020, https://doi.org/10.22059/jcamech.2020.283714.414
  • A. Zargaripoor, A. Daneshmehr, I. Isaac Hosseini, and A. Rajabpoor, “Free vibration analysis of nanoplates made of functionally graded materials based on nonlocal elasticity theory using finite element method,” Journal of Computational Applied Mechanics, Vol. 49, No. Online First, pp. 86–101, Feb. 2018, https://doi.org/10.22059/jcamech.2018.248906.223
  • H. Raissi, M. Shishesaz, and S. Moradi, “Applications of higher order shear deformation theories on stress distribution in a five layer sandwich plate,” Journal of Computational Applied Mechanics, No. Online First, Oct. 2017, https://doi.org/10.22059/jcamech.2017.239207.172
  • D. Rodrigues, J. Belinha, and R. Natal Jorge, “The radial point interpolation method in the bending analysis of symmetric laminates using HSDTS,” Journal of Computational Applied Mechanics, Vol. 52, No. 4, pp. 682–716, Dec. 2021, https://doi.org/10.22059/jcamech.2021.323598.616
  • C. Ike, “Fourier series method for finding displacements and stress fields in hyperbolic shear deformable thick beams subjected to distributed transverse loads,” Journal of Computational Applied Mechanics, Vol. 53, No. 1, pp. 126–141, Mar. 2022, https://doi.org/10.22059/jcamech.2022.332719.658
  • C. C. Ike, C. U. Nwoji, B. O. Mama, H. N. Onah, and M. E. Onyia, “Least squares weighted residual method for finding the elastic stress fields in rectangular plates under uniaxial parabolically distributed edge loads,” Journal of Computational Applied Mechanics, Vol. 51, No. 1, pp. 107–121, Jun. 2020, https://doi.org/10.22059/jcamech.2020.298074.484
  • C. Ike, H. Onah, M. Onyia, B. Mama, and C. Nwoji, “First principles derivation of displacement and stress function for three-dimensional elastostatic problems, and application to the flexural analysis of thick circular plates,” Journal of Computational Applied Mechanics, Vol. 51, No. Online First, pp. 184–198, Feb. 2020, https://doi.org/10.22059/jcamech.2020.295989.471
  • F. C. Onyeka, T. E. Okeke, and B. O. Mama, “Static elastic bending analysis of a three-dimensional clamped thick rectangular plate using energy method,” HighTech and Innovation Journal, Vol. 3, No. 3, pp. 267–281, Aug. 2022, https://doi.org/10.28991/hij-2022-03-03-03
  • F. Onyeka, C. Nwadavid, and B. Okwudili Mama, “Static bending solutions for an isotropic rectangular clamped/simply supported plates using 3-D plate theory,” Journal of Computational Applied Mechanics, Vol. 54, No. 1, pp. 1–18, Mar. 2023, https://doi.org/10.22059/jcamech.2022.349835.764
  • F. Onyeka, T. Edozie Okeke, C. David Nwadavid, and B. Okwudili Mama, “Analytical elasticity solution for accurate prediction of stresses in a rectangular plate bending analysis using exact 3-D theory,” Journal of Computational Applied Mechanics, Vol. 54, No. 2, pp. 167–185, Jun. 2023, https://doi.org/10.22059/jcamech.2022.351892.781
  • O. Festus and E. T. Okeke, “Analytical solution of thick rectangular plate with clamped and free support boundary condition using polynomial shear deformation theory,” Advances in Science, Technology and Engineering Systems Journal, Vol. 6, No. 1, pp. 1427–1439, Feb. 2021, https://doi.org/10.25046/aj0601162
  • F. C. Onyeka and B. O. Mama, “Analytical study of bending characteristics of an elastic rectangular plate using direct variational energy approach with trigonometric function,” Emerging Science Journal, Vol. 5, No. 6, pp. 916–928, Dec. 2021, https://doi.org/10.28991/esj-2021-01320
  • F. C. Onyeka, B. O. Mama, and T. E. Okeke, “Exact three-dimensional stability analysis of plate using a direct variational energy method,” Civil Engineering Journal, Vol. 8, No. 1, pp. 60–80, Jan. 2022, https://doi.org/10.28991/cej-2022-08-01-05
  • Jafar Rouzegar, “A finite element formulation for bending analysis of isotropic and orthotropic plates based on two-variable refined plate theory,” Scientia Iranica, Vol. 22, No. 1, 2015.
  • P. D. Gajbhiye, V. Bhaiya, and Y. M. Ghugal, “Free vibration analysis of thick isotropic plate by using 5th order shear deformation theory,” Progress in Civil and Structural Engineering, Feb. 2021, https://doi.org/10.38208/pcse.v1i1.2
  • P. D. Gajbhiye, V. Bhaiya, and Y. M. Ghugal, “Bending analysis of sandwich plates subjected to various mechanical loadings using quasi-three-dimensional theory,” Journal of Aerospace Engineering, Vol. 35, No. 4, 2022.
  • Y. M. Ghugal and A. S. Sayyad, “A static flexure of thick isotropic plates using trigonometric shear deformation theory,” Journal of Solid Mechanics, Vol. 2, No. 1, pp. 79–90, 2010.
  • Y. M. Ghugal and A. S. Sayyad, “Stress analysis of thick laminated plates using trigonometric shear deformation theory,” International Journal of Applied Mechanics, Vol. 5, No. 1, p. 1350003, Mar. 2013, https://doi.org/10.1142/s1758825113500038
  • A. Zargaripoor, A. Bahrami, and M. Nikkhah Bahrami, “Free vibration and buckling analysis of third-order shear deformation plate theory using exact wave propagation approach,” Journal of Computational Applied Mechanics, Vol. 49, No. Online First, pp. 102–124, Jan. 2018, https://doi.org/10.22059/jcamech.2018.249468.227
  • R. Kumar, M. Singh, C. Kumar, J. Singh, J. Damania, and J. Singh, “Assessment of radial basis function based meshfree method for the buckling analysis of rectangular FGM plate using HSDT and Strong form formulation,” Journal of Computational Applied Mechanics, Vol. 53, pp. 332–347, Jun. 2022, https://doi.org/10.22059/jcamech.2022.342228.716
  • H. Makvandi, S. Moradi, D. Poorveis, and K. Heydari Shirazi, “A new approach for nonlinear vibration analysis of thin and moderately thick rectangular plates under inplane compressive load,” Journal of Computational Applied Mechanics, No. Online First, Oct. 2017, https://doi.org/10.22059/jcamech.2017.240726.181
  • Ö. Civalek and M. Ülker, “Harmonic differential quadrature method (HDQ) for axisymmetric bending analysis of thin isotropic circular plates,” Structural Engineering and Mechanics, Vol. 17, No. 1, pp. 1–14, 2004.
  • Ö. Civalek, “Three-dimensional vibration, buckling and bending analyses of thick rectangular plates based on discrete singular convolution method,” International Journal of Mechanical Sciences, Vol. 49, No. 6, pp. 752–765, Jun. 2007, https://doi.org/10.1016/j.ijmecsci.2006.10.002
  • C. Ike, “Exact analytical solutions to bending problems of SFrSFr thin plates using variational Kantorovich-Vlasov method,” Journal of Computational Applied Mechanics, Vol. 54, No. 2, pp. 186–203, Jun. 2023, https://doi.org/10.22059/jcamech.2023.351563.776
  • K. Soltani et al., “A novel hyperbolic shear deformation theory, for the mechanical buckling analysis of advanced composite plates resting on elastic foundations,” Steel and Composite Structures, Vol. 30, No. 1, pp. 13–29, 2019, https://doi.org/10.12989/scs.2019.30.1.012
  • K. Nareen and R. P. Shimpi, “Refined hyperbolic shear deformation plate theory,” Proceedings of the Institution of Mechanical Engineers Part C: Journal of Mechanical Engineering Science, Vol. 229, No. 15, pp. 2675–2686, 2015, https://doi.org/10.1177/095440621456373
  • R. Kumar and J. Singh, “Assessment of higher order transverse shear deformation theories for modeling and buckling analysis of FGM plates using RBF based meshless approach,” Multidiscipline Modeling in Materials and Structures, Vol. 14, No. 5, pp. 891–907, Oct. 2018, https://doi.org/10.1108/mmms-07-2017-0069
  • Y. Bao, S. Cen, and C. Li, “Distortion-resistant and locking-free eight-node elements effectively capturing the edge effects of Mindlin-Reissner plates,” Engineering Computations, Vol. 34, No. 2, pp. 548–586, Apr. 2017, https://doi.org/10.1108/ec-04-2016-0143
  • F. Onyeka, C. Nwa-David, and T. Edozie, “Analytical solution for the static bending elastic analysis of thick rectangular plate structures using 3-D plate theory,” Engineering and Technology Journal, Vol. 40, No. 11, pp. 1548–1559, Oct. 2022, https://doi.org/10.30684/etj.2022.134687.1244
  • A. Lal, M. Vaghela, and K. Mishra, “Numerical analysis of an edge crack isotropic plate with void/inclusions under different loading by implementing XFEM,” Journal of Applied and Computational Mechanics, No. Online First, Nov. 2019, https://doi.org/10.22055/jacm.2019.31268.1848
  • T. L. Doan, P. B. Le, T. T. Tran, V. K. Trai, and Q. H. Pham, “Free vibration analysis of functionally graded porous nano-plates with different shapes resting on elastic foundation,” Journal of Applied and Computational Mechanics, Vol. 7, No. Online First, pp. 1593–1605, Feb. 2021, https://doi.org/10.22055/jacm.2021.36181.2807
  • L. Hadji and M. Avcar, “Free vibration analysis of FG porous sandwich plates under various boundary conditions,” Journal of Applied and Computational Mechanics, Vol. 7, pp. 505–519, Nov. 2020, https://doi.org/10.22055/jacm.2020.35328.2628
  • A. Assie, S. Mohamed, R. Shanab, R. Abobakr, and M. Eltaher, “Static buckling of 2D FG porous plates resting on elastic foundation based on unified shear theories,” Journal of Applied and Computational Mechanics, Vol. 9, No. Online First, pp. 239–258, Aug. 2022, https://doi.org/10.22055/jacm.2022.41265.3723
  • S. A. Salawu, G. M. Sobamowo, and O. M. Sadiq, “Nonlinear investigation of magnetic influence on dynamic behaviour of non-homogeneous varying thickness circular plates resting on elastic foundations,” Journal of Solid Mechanics, Vol. 13, No. 4, pp. 460–484, Dec. 2021, https://doi.org/10.22034/jsm.2021.1912143.1650
  • F. Boussem and L. Belounar, “A plate bending Kirchhoff element based on assumed strain functions,” Journal of Solid Mechanics, Vol. 12, No. 4, pp. 935–952, Dec. 2020, https://doi.org/10.22034/jsm.2020.1901430.1601
  • P. Roodgar Saffari, M. Fakhraie, and M. A. Roudbari, “Free vibration and transient response of heterogeneous piezoelectric sandwich annular plate using third-order shear deformation assumption,” Journal of Solid Mechanics, Vol. 12, No. 2, pp. 315–333, Jun. 2020, https://doi.org/10.22034/jsm.2019.1865985.1420
  • M. Molla-Alipour, M. Shariyat, and M. Shaban, “Free vibration analysis of bidirectional functionally graded conical/cylindrical shells and annular plates on nonlinear elastic foundations, based on a unified differential transform analytical formulation,” Journal of Solid Mechanics, Vol. 12, No. 2, pp. 385–400, Jun. 2020, https://doi.org/10.22034/jsm.2019.1869981.1450
  • N. Ahlawat and R. Lal, “Effect of Winkler foundation on radially symmetric vibrations of bi-directional FGM non-uniform Mindlin’s circular plate subjected to in-plane peripheral loading,” Journal of Solid Mechanics, Vol. 12, No. 2, pp. 455–475, Jun. 2020, https://doi.org/10.22034/jsm.2019.1873720.1466
  • M. G. Sobamowo and S. A. Salawu, “Free vibration analysis of nonlinear circular plates resting on Winkler and Pasternak foundations,” Journal of Solid Mechanics, Vol. 12, No. 1, pp. 121–135, Mar. 2020, https://doi.org/10.22034/jsm.2019.1866236.1423
  • A. Norouzzadeh, R. Ansari, and H. Rouhi, “Isogeometric analysis of Mindlin nanoplates based on the integral formulation of nonlocal elasticity,” Multidiscipline Modeling in Materials and Structures, Vol. 14, No. 5, pp. 810–827, Oct. 2018, https://doi.org/10.1108/mmms-09-2017-0109
  • T. Do, V. Pham, and H. N. Nguyen, “On the development of refined plate theory for static bending behavior of functionally graded plates,” Mathematical Problems in Engineering, Vol. 2020, pp. 1–13, Apr. 2020, https://doi.org/10.1155/2020/2836763
  • M. Malikan and B. Nguyen, “A novel one-variable first-order shear deformation theory for biaxial buckling of a size-dependent plate based on Eringen’s nonlocal differential law,” World Journal of Engineering, Vol. 15, No. 5, pp. 633–645, Oct. 2018, https://doi.org/10.1108/wje-11-2017-0357
  • V. P. Revenko, “Reduction of a three-dimensional problem of the theory of bending of thick plates to the solution of two two-dimensional problems,” Materials Science, Vol. 51, No. 6, pp. 785–792, May 2016, https://doi.org/10.1007/s11003-016-9903-7
  • Y. Zhong and Q. Xu, “Analysis bending solutions of clamped rectangular thick plate,” Mathematical Problems in Engineering, Vol. 2017, pp. 1–6, 2017, https://doi.org/10.1155/2017/7539276
  • A. J. M. Ferreira and C. M. C. Roque, “Analysis of thick plates by radial basis functions,” Acta Mechanica, Vol. 217, No. 3-4, pp. 177–190, Mar. 2011, https://doi.org/10.1007/s00707-010-0395-5
  • K. H. Lo, R. M. Christensen, and E. M. Wu, “A high-order theory of plate deformation, part 1: homogeneous plates,” ASME Journal of Applied Mechanics, Vol. 44E, No. 4, pp. 663–668, Dec. 1977, https://doi.org/10.1115/1.3424154
  • S. T. Gomaa, M. H. Baluch, H. H. Abdel‐Rahman, and A. K. Mohammed, “Finite element modelling of thick isotropic plates,” Engineering Computations, Vol. 8, No. 4, pp. 361–378, Apr. 1991, https://doi.org/10.1108/eb023845

About this article

Received
05 October 2023
Accepted
09 November 2023
Published
30 December 2023
Keywords
third-order shear deformation plate bending formulation
thick plate
Navier double trigonometric series method
transverse shear stress
transverse normal stress
Acknowledgements

The authors have not disclosed any funding.

Data Availability

The datasets generated during and/or analyzed during the current study are available from the corresponding author on reasonable request.

Conflict of interest

The authors declare that they have no conflict of interest.