Vibration of triangular plates of variable thickness

The problem of the natural frequencies and mode shapes of cantilevered triangular plates with variable thickness and arbitrary planform is solved using the finite element technique. This is done for various combinations of four non-dimensional geometric parameters, namely, the aspect ratio, the two thickness ratios along the two coordinate directions and the sweepback angle. The frequencies for the various cases are tabulated and a few typical mode shapes have been presented graphically.     

Finite element buckling analysis of stiffened plates

An isoparametric stiffened plate bending element for the buckling analysis of stiffened plates has been presented. In the present approach, the stiffener can be positioned anywhere within the plate element and need not necessarily be placed on the nodal lines. The element, being isoparametric quadratic, can readily accommodate curved boundaries, laminated materials and transverse shear deformation. The formulation is applicable to thin as well as thick plates. The buckling loads for various rectangular and skew stiffened plates with varying skew angles and stiffness parameters have been indicated. The results show good agreement with those published.     

Linear free flexural vibration of cracked functionally graded plates in thermal environment

In this paper, the linear free flexural vibrations of functionally graded material plates with a through center crack is studied using an 8-noded shear flexible element. The material properties are assumed to be temperature dependent and graded in the thickness direction. The effective material properties are estimated using the Mori–Tanaka homogenization scheme. The formulation is developed based on first-order shear deformation theory. The shear correction factors are evaluated employing the energy equivalence principle. The variation of the plates natural frequency is studied considering various parameters such as the crack length, plate aspect ratio, skew angle, temperature, thickness and boundary conditions. The results obtained here reveal that the natural frequency of the plate decreases with increase in temperature gradient, crack length and gradient index.     

Design of laminated composite plates for optimal dynamic characteristics using a constrained global optimization technique

The lamination arrangements of moderately thick laminated composite plates for optimal dynamic characteristics are studied via a constrained multi-start global optimization technique. In the optimization process, the dynamical analysis of laminated composite plates is accomplished by utilizing a shear deformable laminated composite finite element, in which the exact expressions for determining shear correction factors were adopted and the modal damping model constructed based on an energy concept. The optimal layups of laminated composite plates with maximum fundamental frequency or modal damping are then designed by maximizing the frequency or modal damping capacity of the plate via the multi-start global optimization technique. The effects of length-to-thickness ratio, aspect ratio and number of layer groups upon the optimum fiber orientations or layer group thicknesses are investigated by means of a number of examples of the design of symmetrically laminated composite plates.     

An equilibrium method for prediction of transverse shear stresses in a thick laminated plate

First two equations of equilibrium are utilized to compute the transverse shear stress variation through thickness of a thick laminated plate after in-plane stresses have been computed using an assumed quadratic displacement triangular element based on transverse inextensibility and layerwise constant shear angle theory (LCST). Centroid of the triangle is the point of exceptional accuracy for transverse shear stresses. Numerical results indicate close agreement with elasticity theory. An interesting comparison between the present theory and that based on assumed stress hybrid finite element approach suggests that the latter does not satisfy the condition of free normal traction at the edge. Comparison with numerical results obtained by using constant shear angle theory suggests that LCST is close to the elasticity solution while the CST is closer to classical (CLT) solution. It is also demonstrated that the reduced integration gives faster convergence when the present theory is applied to a thin plate.     

A simple displacement-based 3-node triangular element for linear and geometrically nonlinear analysis of laminated composite plates

A simple displacement-based 3-node, 18-degree-of-freedom flat triangular plate/shell element LDT18 is proposed in this paper for linear and geometrically nonlinear finite element analysis of thin and thick laminated composite plates. The presented element is based on the first-order shear deformation theory (FSDT), and the total Lagrangian approach is employed to formulate the element for geometrically nonlinear analysis. The deflection and rotation functions of the element boundary are obtained from the Timoshenko’s laminated composite beam functions, hence convergence to the thin plate solution can be achieved theoretically and shear-locking problem is avoided naturally. The plane displacement interpolation functions of the Airman’s triangular membrane element with drilling degrees of freedom are taken as the in-plane displacements of the element. Numerical examples demonstrate that the present element is accurate and efficient for linear and geometrically nonlinear analysis of thin to moderately thick laminated composite plates.     

Nonlinear analysis of skew plates using the finite element method

In this paper finite element analysis of the large deflection behaviour of skew plates has been done. A high precision conforming triangular plate bending element has been used. The central deflection, bending and membrane stresses have been reported for simply supported and clamped rhombic plates. The variations of these quantities have been studied for different skew angles.     

Higher-order plate bending analysis using the modified Razzaque-Irons triangle

The Razzaque-Irons (RI) triangular, Kirchhoff plate bending element is modified using a higher-order shear deformation plate bending theory which was introduced by Levinson, Murthy and Reddy. This theory allows a quadratic distribution of transverse shear strains through the plate thickness. Examples are provided to demonstrate the accuracy of this new element.     

A three-dimensional nonlinear analysis of cross-ply rectangular composite plates

The results of a three-dimensional, geometrically nonlinear, finite-element analysis of the bending of cross-ply laminated anisotropie composite plates are presented. Individual laminae are assumed to be homogeneous, orthotropic and linearly elastic. A fully three-dimensional isoparametric finite element with eight nodes (i.e. linear element) and 24 degrees of freedom (three displacement components per node) is used to model the laminated plate. The finite element results of the linear analysis are found to agree very well with the exact solutions of cross-ply laminated rectangular plates under sinusiodal loading. The finite element results of the three-dimensional, geometrically nonlinear analysis are compared with those obtained by using a shear deformable, geometrically nonlinear, plate theory. It is found that the deflections predicted by the shear deformable plate theory are in fair agreement with those predicted by three-dimensional elasticity theory; however stresses were found to be not in good agreement     

A study of moderately thick quadrilateral plate elements based on the absolute nodal coordinate formulation

Finite element analysis using plate elements based on the absolute nodal coordinate formulation (ANCF) can predict the behaviors of moderately thick plates subject to large deformation. However, the formulation is subject to numerical locking, which compromises results. This study was designed to investigate and develop techniques to prevent or mitigate numerical locking phenomena. Three different ANCF plate element types were examined. The first is the original fully parameterized quadrilateral ANCF plate element. The second is an update to this element that linearly interpolates transverse shear strains to overcome slow convergence due to transverse shear locking. Finally, the third is based on a new higher order ANCF plate element that is being introduced here. The higher order plate element makes it possible to describe a higher than first-order transverse displacement field to prevent Poisson thickness locking. The term “higher order” is used, because some nodal coordinates of the new plate element are defined by higher order derivatives. The performance of each plate element type was tested by (1) solving a comprehensive set of small deformation static problems, (2) carrying out eigenfrequency analyses, and (3) analyzing a typical dynamic scenario. The numerical calculations were made using MATLAB. The results of the static and eigenfrequency analyses were benchmarked to reference solutions provided by the commercially available finite element software ANSYS. The results show that shear locking is strongly dependent on material thickness. Poisson thickness locking is independent of thickness, but strongly depends on the Poisson effect. Poisson thickness locking becomes a problem for both of the fully parameterized element types implemented with full 3-D elasticity. Their converged results differ by about 18 % from the ANSYS results. Corresponding results for the new higher order ANCF plate element agree with the benchmark. ANCF plate elements can describe the trapezoidal mode; therefore, they do not suffer from Poisson locking, a reported problem for fully parameterized ANCF beam elements. For cases with shear deformation loading, shear locking slows solution convergence for models based on either the original fully parameterized plate element or the newly introduced higher order plate element.     

Large deflection analysis of skew plates by lumped triangular element formulation

A finite difference scheme with triangular mesh is presented for the analysis of skew plate problems with large deflections. The suggested formulation is independent of the boundary condition and uses energy principles to derive a set of nonlinear algebraic equations which are solved by using Newton-Raphson iterative method with incremental loading. The investigation is concerned with the behaviour of constant thickness clamped and simply supported isotropic skew plates with immovable edges and subjected to uniformly distributed transverse load. The effects of skew on plates with large deflections are investigated and comparisons are made with existing results; good agreement is shown.     

Vibration of stiffened skew plates by using B-spline functions

This paper presents a general procedure for calculating the free vibration of stiffened skew plates by the Rayleigh-Ritz method with B-spline functions as coordinate functions. The stiffened skew plates are modelled as the skew plate with a number of stiffening beams.The results are compared with existing values based on other numerical methods. Vibration characteristics of stiffened skew plates are also studied with changing the arrangements of stiffening beams, the stiffness parameters of beams, skew angle and aspect ratio.     

A refined higher-order C° plate bending element

A general finite element formulation for plate bending problem based on a higher-order displacement model and a three-dimensional state of stress and strain is attempted. The theory incorporates linear and quadratic variations of transverse normal strain and transverse shearing strains and stresses respectively through the thickness of the plate. The 9-noded quadrilateral from the family of two dimensional C° continuous isoparametric elements is then introduced and its performance is evaluated for a wide range of plates under uniformly distributed load and with different support conditions and ranging from very thick to extremely thin situations. The effect of full, reduced and selective integration schemes on the final numerical result is examined. The behaviour of this element with the present formulation is seen to be excellent under all the three integration schemes.     

Finite element flutter analysis of composite skew panels

A finite element method is used to study the dynamic instability of laminated composite skew panels subjected to supersonic flow. The FEM employs eight-noded isoparametric elements which take into account transverse shear deformation. The linearized Piston theory is applied to assess the aerodynamic loads. The effects of skew angle on critical aerodynamic parameters are investigated for different aspect ratios, boundary constraints, fibre orientations and lamination schemes. It is observed that the skew angle has a stabililzing effect on the flutter boundary, whereas couplings have a destabilizing effect. The higher aspect ratio is also found to exhibit a stabilizing effect on the flutter boundary.     

An approximate semi-analytical method for prediction of interlaminar shear stresses in an arbitrarily laminated thick plate

An approximate semi-analytical method for determination of interlaminar shear stress distribution through the thickness of an arbitrarily laminated thick plate has been presented. The method is based on the assumptions of transverse inextensibility and layerwise constant shear angle theory (LCST) and utilizes an assumed quadratic displacement potential energy based finite element method (FEM). Centroid of the triangular surface has been proved, from a rigorous methematical point of view (Aubin-Nitsche theory), to be the point of exceptional accuracy for the interlaminar shear stresses. Numerical results indicate close agreement with the available three-dimensional elasticity theory solutions. A comparison between the present theory and that due to an assumed stress hybrid FEM suggests that the (normal) traction-free-edge condition is not satisfied in the latter approach. Furthermore, the present paper is the first to present the results for interlaminar shear stresses in a two-layer thick square plate of balanced unsymmetric angle-ply construction. A comparison with the recently proposed Equilibrium Method (EM) indicates the superiority of the present method, because the latter assures faster convergence as well as simultaneous vanishing of the transverse shear stresses on both the exposed surfaces of the laminate. Superiority of the present method over the EM, in the case of a symmetric laminate, is limited to faster convergence alone. It has also been demonstrated that the combination of the present method and the reduced (quadratic order) numerical integration scheme yields convergence of the interlaminar shear stresses almost as rapidly as that of the nodal displacements, in the case of a thin plate.     

The free vibration of skew plates using the hierarchical finite element method

The hierarchical finite element method is used to determine the natural frequencies and modes of flat, isotropic skew plates. A number of such plates with different boundary conditions—including free edges and point supports—are considered in this paper. The dependence of frequency on skew angle, aspect ratio and Poisson's ratio is investigated, though succinctness prohibits a complete study exploring the full interrelation of these parameters. Extensive results are presented in diagrammatic, graphical, and tabular format; these are shown to be in very good agreement with the work of other investigators, and should prove a valuable source of data for use by engineers and scientists.     

Stability of moderately thick rectangular plates using a high precision triangular finite element

Stability of moderately thick rectangular plates is studied in this paper. A high precision triangular finite element is used in the present formulation. Explicit expressions for elastic stiffness and geometric stiffness matrices are presented; the use of numerical integration is avoided. The effect of transverse shear on the critical loads is brought out through two typical examples.     

基于修正变分法开孔板和加强板结构的自由振动分析

为研究船舶开孔板和加强板结构的振动特性,用1阶剪切变形板理论描述各向同性板的位移场,并采用修正变分原理和区域分解方法建立板的离散动力学模型.每一块子域板的位移和转角分量通过第一类切比雪夫正交多项式展开.针对加强板模型,将该方法获得的结果与已经发表的文献和有限元商用软件计算结果进行对比,验证该方法的收敛性和正确性.基于修正变分法探讨多种开孔和加强板模型的自由振动特性,充分说明该数理模型和半解析方法是一种适合处理复杂板结构问题的数值工具.     

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

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.     

Free vibration of laminated plates using a finite strip method based on a higher-order plate theory

The finite strip method based on the higher-order plate theory is developed for determining the natural frequencies of laminated plates. This method can accurately predict the through thickness effect of transverse shear deformation. Furthermore, only a few degress of freedom are required in the finite strip method. Some numerical results for various span-to-thickness ratios, material properties and stack sequences are presented for illustrative purposes. The present model provides a better way to obtain more accurate natural frequency results.