Linear and non-linear deflection analysis of thick rectangular plates—II. Numerical applications

Abstract

Variational methods are widely used for the solution of complex differential equations in mechanics for which exact solutions are not possible. The finite difference method, although well known as an efficient numerical method, was applied in the past only for the analysis of linear and non-linear thin plates. In this paper the suitability of the method for the analysis of non-linear deflection of thick plates is studied for the first time. While there are major differences between small deflection and large deflection plate theories, the former can be treated as a particular case of the latter, when the centre deflection of the plate is less than or equal to 0.2–0.25 of the thickness of the plate. The finite difference method as applied here is a modified finite difference approach to the ordinary finite difference method generally used for the solution of thin plate problems. In this analysis thin plates are treated as a particular case of the corresponding thick plate when the boundary conditions of the plates are taken into account. The method is first applied to investigate the deflection behaviour of clamped and simply supported square isotropic thick plates. After the validity of the method is established, it is then extended to the solution of rectangular thick plates of various aspect ratios and thicknesses. Generally, beginning with the use of a limited number of mesh sizes for a given plate aspect ratio and boundary conditions, a general solution of the problem including the investigation of accuracy and convergence was extended to rectangular thick plates by providing more detailed functions satisfying the rectangular mesh sizes generated automatically by the program. Whenever possible results obtained by the present method are compared with existing solutions in the technical literature obtained by much more laborious methods and close agreements are found. The significant number of results presented here are not currently available in the technical literature. The submatrices involved in the formation of the finite difference equations from the governing differential equations are generated directly by the computer program. The subroutine SOLINV using the change of variable technique illustrated elsewhere takes care of the solution of the general system. Simplicity in formulation and quick convergence are the obvious advantages of the finite difference formulation developed to compute small and large deflection analysis of thick plates in comparison with other numerical methods requiring extensive computer facilities.