K.P. Rao

A doubly curved quadrilateral finite element for the analysis of laminated anisotropic thin shells of revolution

Abstract The details of development of the stiffness matrix for a doubly curved quadrilateral element suited for static and dynamic analysis of laminated anisotropic thin shells of revolution are reported. Expressing the assumed displacement state over the middle surface of the shell as products of one-dimensional first order Hermite polynomials, it is possible to ensure that the displacement state for the assembled set of such elements, is geometrically admissible. Monotonic convergence of total potential energy is therefore possible as the modelling is successively refined. Systematic evaluation of performance of the element is conducted, considering various examples for which analytical or other solutions are available.

A laminated anisotropic curved beam and shell stiffening finite element

The details of development of the stiffness matrix of a laminated anisotropic curved beam finite element are reported. It is a 16 dof element which makes use of 1-D first order Hermite interpolation polynomials for expressing its assumed displacement state. The performance of the element is evaluated considering various examples for which analytical or other solutions are available.

Interlaminar stresses in laminated composite plates, cylindrical/spherical shell panels damaged by low-velocity impact

Prediction of damage caused by low-velocity impact in laminated composite plate cylindrical/spherical shell panels is an important problem faced by designers using composites. Not only the in-plane stresses but also the interlaminar normal and shear stresses play a role in estimating the damage caused. The work reported here is an effort in getting better predictions of damage in composite plate cylindrical/spherical shell panels subjected to low velocity impact. The low-velocity impact problem is treated as a quasi-static problem. First, the in-plane stresses are calculated by 2-D nonlinear finite element analysis using a 48 degrees of freedom laminated composite shell element. The damage analysis is then carried out using a Tsai-Wu quadratic failure criterion and a maximum stress criteria. Interlaminar normal and shear stresses are predicted after taking into account the in-plane damage caused by low velocity impact. The interlaminar stresses are obtained by integrating the 3-D equations of equilibrium through the thickness. The deformed geometry is taken into account in the third equation of equilibrium (in the thickness direction). After evaluating the formulation and the computer program developed for correctness, the interlaminar stresses are predicted for composite plates/shell panels which are damaged by low-velocity impact

Stiffened composite axisymmetric shells-optimum lay-up for buckling by ranking

Linear bifurcation buckling of FRP axisymmetric shells with fully compatible FRP meridional and hoop stiffeners is studied using the finite element method. Eccentricity of the stiffeners is taken into account. The composite shell and the stiffener are assumed to be made of a repeated sublaminate construction. This type of construction is used in industry to reduce manufacturing errors and to produce more damage-tolerant laminates. In this type of construction, the sublaminate consists of smaller number of plies and the required thickness of the laminate is obtained by repeating the sublaminate many times. This paper deals with the determination of the optimum lay-up scheme in the sublaminate of a composite axisymmetric shell with composite stiffener elements so as to achieve maximum buckling load for a given geometry, loading and boundary conditions using the finite element method. A four-noded, 48-DOF doubly curved quadrilateral laminated anisotropic thin shell finite element with fully compatible two-noded, 16-DOF meridional stiffener elements (MSE) and parallel circle stiffener elements (PCSE) is used. The buckling loads computed for several cases of shells (solid/stiffened) of positive and negative Gaussian curvatures with different applied loads and boundary conditions compare well with existing results in the literature. Subsequently the computer program has been used to find the optimum lay-up scheme of the plies in the sublaminate so as to achieve maximum buckling load for typical composite solid/stiffened shells.

Stiffened composite cylindrical panels—Optimum lay-up for buckling by ranking

Buckling of discretely stiffened composite cylindrical panels made of repeated sublaminate construction is studied using a finite element method. In repeated sublaminate construction, a full laminate is obtained by repeating a basic sublaminate, which has a smaller number of plies. This paper deals with the determination of the optimum lay-up for buckling by ranking of such stiffened (longitudinal and hoop) composite cylindrical panels. For this purpose we use the particularized form of a four-noded, 48 degrees of freedom doubly curved quadrilateral thin shell finite element together with a fully compatible two-noded, 16 degrees of freedom composite stiffener element. The computer program developed has been used, after extensive checking for correctness, to obtain an optimum orientation scheme of the plies in the sublaminate so as to achieve maximum buckling load for a specified thickness of typical stiffened composite cylindrical panels.

Large deformation finite element analysis of laminated circular composite plates

A 48 d.o.f., four-noded quadrilateral laminated composite shell finite element is particularised to a sector finite element and is used for the large deformation analysis of circular composite laminated plates. The strain-displacement relationships for the sector element are obtained by reducing those of the quadrilateral shell finite element by substituting proper values for the geometric parameters. Subsequently, the linear and tangent stiffness matrices are formulated using conventional methods. The Newton-Raphson method is employed as the nonlinear solution technique. The computer code developed is validated by solving an isotropic case for which results are available in the literature. The method is then applied to solve problems of cylindrically orthotropic circular plates. Some of the results of cylindrically orthotropic case are compared with those available in the literature. Subsequently, application is made to the case of laminated composite circular plates having different lay-up schemes. The computer code can handle symmetric/unsymmetric lay-up schemes. The large displacement analysis is useful in estimating the damage in composite plates caused by low-velocity impact.

Composite cylindrical panels—optimum lay-up for buckling by ranking

Fibre-reinforced composite cylindrical panels are a typical subcomponent used in various aerospace structures. Buckling is one of the important modes of failure of such a panel. Repeated sublaminate construction is often used in practice to reduce manufacturing errors and to produce more damage-tolerant laminates. In this type of construction, the full laminate is obtained by repeating the basic sublaminate, which has a smaller number of plies. This paper deals with the determination of optimum lay-up for buckling by ranking of such fibre-reinforced plastic composite cylindrical panels (which may be solid, sandwich or stiffened) using the finite element method. The particularized form of a four-noded, 48-dof doubly curved quadrilateral laminated anisotropic thin shell finite element is used. The core/stiffener can be a corrugated composite sheet or a regular grid. The computer program developed has been used, after extensive checking for correctness, to obtain the optimum orientation scheme of the plies in the sublaminate so as to achieve maximum buckling load for typical composite cylindrical panels.

Curved composite beams—optimum lay-up for buckling by ranking

Instability of laminated curved composite beams made of repeated sublaminate construction is studied using finite element method. In repeated sublaminate construction, a full laminate is obtained by repeating a basic sublaminate which has a smaller number of plies. This paper deals with the determination of optimum lay-up for buckling by ranking of such composite curved beams (which may be solid or sandwich). For this purpose, use is made of a two-noded, 16 degress of freedom curved composite beam finite element. The displacements u, v, w of the element reference axis are expressed in terms of one-dimensional first-order Hermite interpolation polynomials, and line member assumptions are invoked in formulation of the elastic stiffness matrix and geometric stiffness matrix. The nonlinear expressions for the strains, occurring in beams subjected to axial, flexural and torsional loads, are incorporated in a general instability analysis. The computer program developed has been used, after extensive checking for correctness, to obtain optimum orientation scheme of the plies in the sublaminate so as to achieve maximum buckling load for typical curved solid/sandwich composite beams.

Finite element analysis of laminated anisotropic beams of bimodulus materials

This paper presents finite element analysis of laminated anisotropic beams of bimodulus materials. The finite element has 16 d.o.f. and uses the displacement field in terms of first order Hermite interpolation polynomials. As the neutral axis position may change from point to point along the length of the beam, an iterative procedure is employed to determine the location of zero strain points along the length. Using this element some problems of laminated beams of bimodulus materials are solved for concentrated loads/moments perpendicular and parallel to the layering planes as well as combined loads.

Estimation of low-velocity impact damage in laminated composite circular plates using nonlinear finite element analysis

Nonlinear finite element analysis is used for the estimation of damage due to low-velocity impact loading of laminated composite circular plates. The impact loading is treated as an equivalent static loading by assuming the impactor to be spherical and the contact to obey Hertzian law. The stresses in the laminate are calculated using a 48 d.o.f. laminated composite sector element. Subsequently, the Tsai-Wu criterion is used to detect the zones of failure and the maximum stress criterion is used to identify the mode of failure. Then the material properties of the laminate are degraded in the failed regions. The stress analysis is performed again using the degraded properties of the plies. The iterative process is repeated until no more failure is detected in the laminate. The problem of a typical T300/N5208 composite [45 degrees/0 degrees/-45 degrees/90 degrees](s) circular plate being impacted by a spherical impactor is solved and the results are compared with experimental and analytical results available in the literature. The method proposed and the computer code developed can handle symmetric, as well as unsymmetric, laminates. It can be easily extended to cover the impact of composite rectangular plates, shell panels and shells.