• Magazine of the Korean Society of Agricultural Engineers
• /
• v.22 no.2
• /
• pp.83-100
• /
• 1980

This study was intended to computerize the structural analysis of two-dimensional continuum by finite element method, and to provide a preparatory basis for more sophisticated and more generalized computer programs of this kind. A computer program, applicable to any shape of two-dimensional continuum, was formulated on the basis of 16-degree-of- freedom rectangular element. Various computational aspects pertaining to the implementation of finite element method were reviewed and settled in the course of programming. The validity of the program was checked through several case studies. To assess the accuracy and the convergence characteristics of the method, the results computed by the program were compared with solutions by other methods, namely the analytical Navier's method and the framework method. Through actual programming and analysis of the computed results, the following facts were recognized; 1) The stiffness matrix should necessarily be assembled in a condensed form in order to make it possible to discretize the continuum into practically adequate number of elements without using back-up storage. 2) For minimization of solution time, in-core solution of the equilibrium equation is essential. LDLT decomposition is recommended for stiffness matrices condensed by the compacted column storage scheme. 3) As for rectangular plates, the finite element method shows better performances both in the accuracy and in the rate of convergence than the framework method. As the number of elements increases, the error of the finite element method approaches around 1%. 4) Regardless of the structural shape, there is a uniform tendency in convergence characteristics dependent on the shape of element. Square elements show the best performance. 5) The accuracy of computation is independent of the interpolation function selected.

• Structural Engineering and Mechanics
• /
• v.74 no.3
• /
• pp.365-380
• /
• 2020

The focus of this paper is to develop an analytical approach based on an efficient shear deformation theory with stretching effect for bending stress analysis of cross-ply laminated composite plates subjected to transverse parabolic load and line load by using a new kinematic model, in which the axial displacements involve an undetermined integral component in order to reduce the number of unknowns and a sinusoidal function in terms of the thickness coordinate to include the effect of transverse shear deformation. The present theory contains only five unknowns and satisfies the zero shear stress conditions on the top and bottom surfaces of the plate without using any shear correction factors. The governing differential equations and its boundary conditions are derived by employing the static version of principle of virtual work. Closed-form solutions for simply supported cross-ply laminated plates are obtained applying Navier's solution technique, and the numerical case studies are compared with the theoretical results to verify the utility of the proposed model. Lastly, it can be seen that the present outlined theory is more accurate and useful than some higher-order shear deformation theories developed previously to study the static flexure of laminated composite plates.