The development of triangular and quadrilateral membrane, plate and shell elements in Java using the object oriented programming technique is presented. The membrane elements developed are the constant strain triangle (CST) element and the four node isoparametric quadrilateral membrane element. The plate bending elements developed are the discrete Kirchoff triangular (DKT) element and discrete Kirchoff quadrilateral (DKQ) element. The flat shell elements are developed by super imposing the stiffness of the membrane element and plate bending element. A finite element analysis program is also developed to check the accuracy of the developed elements. The program is developed using the object oriented programming approach as an alternative to traditional procedural programming. Several test structures are analyzed using the developed program for each developed element and the results are compared with those obtained from the commercial finite element analysis program SAP 2000. The results indicate that all elements give accurate displacements. However, there were significant differences in stresses for the shell elements, which can be attributable to the approximate approach in these elements to model the drilling degree of freedom. / Master of Science

Finite element analysis of laminated shells using a three-noded flat triangular shell element is presented. The flat shell element is obtained by combining the Discrete Kirchhoff Theory (DKT) plate bending element and a membrane element similar to the Allman element, but derived from the Linear Strain Triangular (LST) element. Though this combination has been employed in the literature for linear static analysis of laminated plates, the results presented are not adequate to ascertain that the element would perform well in the case of static and dynamic analysis of general shells. The element is first thoroughly tested for linear static analysis of laminated plates and shells and is extended for free vibration, thermal, and geometrically nonlinear analysis. The major drawback of the DKT plate bending element is that the transverse displacement is not explicitly defined within the interior of the element. Hence obtaining the consistent mass matrix or the derivatives of the transverse displacement that are required for forming the geometric stiffness matrix is not straight forward. This problem is alleviated by borrowing shape functions from other similar elements or using simple displacement fields. In the present research, free vibration analysis is performed both by using a lumped mass matrix and a so called consistent mass matrix, obtained by borrowing shape functions from an existing element, in order to compare the performance of the two methods. The geometrically nonlinear analysis is performed using an updated Lagrangian formulation employing Green strain and Second Piola-Kirchhoff (PK2) stress measures. A linear displacement field is used for the transverse displacement in order to compute the derivatives of the transverse displacement that are required to compute the geometric stiffness or the initial stress matrix. Several numerical examples are solved to demonstrate the accuracy of the formulation for both small and large rotation analysis of laminated plates and shells. The results are compared with those available in the existing literature and those obtained using the commercial finite element package ABAQUS and are found to be in good agreement. The element is employed for two main applications involving large flexible structures. The first application is the control of thermal deformations of a spherical mirror segment, which is a segment of a multi-segmented primary mirror used in a space telescope. The feasibility of controlling the surface distortions of the mirror segment due to arbitrary thermal fields, using discrete and distributed actuators, is studied. This kind of study was required for the design of a multi-segmented primary mirror of a next generation space telescope. The second application is the analysis of an inflatable structure, being considered by the US Army for housing vehicles and personnel. The tent structure is made up of membranes supported by arches stiffened by internal pressure. The updated Lagrangian formulation of the flat shell element has been developed primarily for the nonlinear analysis of the tent structure, since such a structure is expected to undergo large deformations and rotations under the action of environmental loads like the wind and snow loads. The wind load is modeled as a nonuniform pressure load and the snow load as lumped concentrated loads. Since the direction of the pressure load is assumed to be normal to the current configuration of the structure, it changes as the structure undergoes deformation. This is called the follower action. As a result, the pressure load is a function of the displacements and hence contributes to the tangent stiffness matrix in the case of geometrically nonlinear analysis. The thermal load also contributes to the system tangent stiffness matrix. In the case of the thermal load this contribution is similar to the initial stress matrix and hence no additional effort is required to compute this contribution. In the case of the pressure load, this contribution (called the pressure stiffness) is in general unsymmetric but can be systematically derived from the principle of virtual work. The follower effects of the pressure load have been included in the updated Lagrangian formulation of the flat shell element and have been validated using standard examples in the literature involving deformation-dependent pressure loads. The element can be used to obtain the nonlinear response of the tent structure under wind and snow loads. / Ph. D.

Experimental Spatial Dynamics Modeling (ESDM) is the new process of constructing a three dimensional, complex-valued dynamic model of a harmonically vibrating structure using numerical models and laser-based experimental data obtained from a Scanning Laser Doppler Vibrometer (SLDV). In ESDM process, a finite element formulation is used to construct a numerical model of a structure. A conventional finite element such as rod, beam, or plate element, can be used to construct the numerical model of a structure from its mid-plane. In this research, the damped flat shell element is developed to construct the numerical models of a cantilever beam and a simply supported flat plate. The velocity transformation function developed in this research will make possible to use the FE model, constructed by the damped flat shell element, and the laser-based experimental data within a framework of ESDM in the consistent manner. / Master of Science