Plate Bending

Plate Bending

1. Introduction

A plate is a flat body whose thickness is much smaller than its other dimensions. Plate carries a lateral load by bending. A plate develops bending moments in two directions and a twisting moment. Two plate theories are well known depending on whether transverse shear deformation is considered or not considered.

1.1 Thin plate theory/ Kirchhoff plate theory

The basic assumption for the classical Kirchhoff plate bending theory is that a straight line normal to the midplane of the plate before deformation remains normal even after deformation. In this theory transverse shear deformation is neglected.

1.2 Mindlin theory/ Mindlin—Reissner theory

The plane normal to the midplane before deformation will not remain normal any longer after the deformation. In this theory transverse shear deformation is accounted.

Figure 1: Plate, geometry and degrees of freedom

Figure 1 shows the plate, its geometry and degrees of freedom. Plate is placed in xy plane. θx and θy represent rotations along Y-axes and X-axes respectively. A plate has a thickness t and has a midsurface at a distance t/2 from each lateral surface. We locate plate midsurface along the xy plane, z = 0 is the plate mid surface. Length of the plate is ‘a’ and its breadth ‘b’.

2. Kirchhoff plate theory

In this theory the transverse deformation is neglected. The in-plane displacement u, v are expresses as

[2.1]

Where x and y are the in-plane axes located at the mid-plane of the plate, and z is along the thickness of the plate as shown in Figure 1. Thus u, v are the displacements along the x-axes and y-axes respectively, and w is the transverse displacement or deflection along the z-axes.

The inplane strains in terms of the displacements are written as

[2.2]

As the transverse deformation is not considered, we have

[2.3]

For a thin plate assuming the plane stress condition, we have

[2.4]

The above stresses give rise to the following moments

[2.5]

Thus the moment-curvature relations for a homogeneous and isotropic Kirchhoff plate are

[2.6]

Where D is flexural rigidity, given by

[2.7]

3. Mindlin plate theory

It assumes that transverse deformation occurs in the plate. The deformations and strains are given by

[3.1]

The moment curvature relations will be for bending and shear deformations given as follows

[3.2]

For bending and

[3.3]

For the transverse shear deformation. Where Qx and Qy are transverse shear force.

Stress on cross section of the plate is shown in Figure 2.

Figure 2: Stresses and distributed lateral forces on plate

4. Finite Element Method for Plate

Mindlin plate theory is considered in finite elements of plates. The displacement interpolation for plates is given by

[4.1]

Where Ni is the shape function associated with the plate. Most common shape functions used for the plate are Q4 and Q8.

5. Example Problem

A code PlateBending is written in MATLAB to calculate the bending displacement of a thin plate under uniform transverse pressure. Two different boundary conditions are considered. Plate is dicretized using isoparametric Q4 elements. The values obtained with present code are compared with standard FEM software and are in good agreement.

A simple supported and a clamped square plate under uniform transverse pressure is considered. The following is the data used in the code.

Young’s Modulus, E = 10920 Nm/Kg^2

Plate thickness, t = 0.1 m

Poisson’s ratio, υ = 0.3

Length of the plate, a = 1 m

Breadth of the plate, b = 1m

Pressure = 1 unit

Four noded Isoparametric elements are used. Plate is meshed with 200 elements. Meshing of the plate and the transverse load applied on the plate are shown in the Figure 3. Static analysis is done and the values obtained with present code are compared with standard FEM software. To show the accuracy of the code, values of the present code are compared with standard FEM software. Table 1 shows the transverse displacement ‘w’ obtained using present code and FEM software. The transverse displacement is written in the non-dimensional from as given in Eq (5.1)

Figure 3: Plate under Uniform transverse pressure

[5.1]

Table 1 : Transverse displacement of the plate

Figure 4 shows the deformation of the simply supported plate under transverse pressure.

Figure 4: Plate deformation under transverse pressure

Figure 5 shows the contour plot of transverse displacement of the simply supported plate from FEM software and Figure 6 shows the contour plot from the present code. Both the plots are in good agreement. Figure 7 shows the contour plot transverse displacement on the deformed plate.

Figure 5: Contour plot of w from FEM software

Figure 6: Contour plot of w from present code

Figure 7: Profile of w on deformed plate

Code PlateBending can obtained from the following link:

Plate Bending:

http://www.mathworks.de/matlabcentral/fileexchange/32029-plate-bending

Plane Stress << Prev Next >> Stress Recovery