A new insight in stress recovery tecniques and hessian reconstruction methods

Stress recovery techniques have been an active research topic in the last few years since, in 1987, Zienkiewicz and Zhu proposed a procedure called Superconvergent Patch Recovery (SPR). This procedure is a last-squares fit of stresses at super-convergent points over patches of elements and it leads to enhanced stress fields that can be used for evaluating finite element discretization errors. In subsequent years, numerous improved forms of this procedure have been proposed attempting to add equilibrium constraints to improve its performances. Later, another superconvergent technique, called Recovery by Equilibrium in Patches (REP), has been proposed. In this case the idea is to impose equilibrium in a weak form over patches and solve the resultant equations by a last-square scheme. In recent years another procedure, based on minimization of complementary energy, called Recovery by Compatibility in Patches (RCP) has been proposed in. This procedure, in many ways, can be seen as the dual form of REP as it substantially imposes compatibility in a weak form among a set of self-equilibrated stress fields. In this thesis a new insight in RCP is presented and the procedure is improved aiming at obtaining convergent second order derivatives of the stress resultants. In order to achieve this result, two different strategies and their combination have been tested. The first one is to consider larger patches in the spirit of what proposed in [4] and the second one is to perform a second recovery on the recovered stresses. Some numerical tests in plane stress conditions are presented, showing the effectiveness of these procedures. Afterwards, a new recovery technique called Last Square Displacements (LSD) is introduced. This new procedure is based on last square interpolation of nodal displacements resulting from the finite element solution. In fact, it has been observed that the major part of the error affecting stress resultants is introduced when shape functions are derived in order to obtain strains components from displacements. This procedure shows to be ultraconvergent and is extremely cost effective, as it needs in input only nodal displacements directly coming from finite element solution, avoiding any other post-processing in order to obtain stress resultants using the traditional method. Numerical tests in plane stress conditions are than presented showing that the procedure is ultraconvergent and leads to convergent first and second order derivatives of stress resultants. In the end, transverse stress profiles reconstruction using First-order Shear Deformation Theory for laminated plates and three dimensional equilibrium equations is presented. It can be seen that accuracy of this reconstruction depends on accuracy of first and second derivatives of stress resultants, which is not guaranteed by most of available low order plate finite elements. RCP and LSD procedures are than used to compute convergent first and second order derivatives of stress resultants ensuring convergence of reconstructed transverse shear and normal stress profiles respectively. Numerical tests are presented and discussed showing the effectiveness of both procedures.

Transverse vibrations of hybrid laminated plates

In this paper, an attempt is made to obtain the free vibration response of hybrid, laminated rectangular and skew plates. The Galerkin technique is employed to obtain an approximate solution of the governing differential equations. It is found that this technique is well suited for the study of such problems. Results are presented in a graphical form for plates with one pair of opposite edges simply supported and the other two edges clamped. The method is quite general and can be applied to any other boundary conditions.

On Uniform Approximate Solutions in Bending of Symmetric Laminated Plates

A layer-wise theory with the analysis of face ply independent of lamination is used in the bending of symmetric laminates with anisotropic plies. More realistic and practical edge conditions as in Kirchhoff's theory are considered. An iterative procedure based on point-wise equilibrium equations is adapted. The necessity of a solution of an auxiliary problem in the interior plies is explained and used in the generation of proper sequence of two dimensional problems. Displacements are expanded in terms of polynomials in thickness coordinate such that continuity of transverse stresses across interfaces is assured. Solution of a fourth order system of a supplementary problem in the face ply is necessary to ensure the continuity of in-plane displacements across interfaces and to rectify inadequacies of these polynomial expansions in the interior distribution of approximate solutions. Vertical deflection does not play any role in obtaining all six stress components and two in-plane displacements. In overcoming lacuna in Kirchhoff's theory, widely used first order shear deformation theory and other sixth and higher order theories based on energy principles at laminate level in smeared laminate theories and at ply level in layer-wise theories are not useful in the generation of a proper sequence of 2-D problems converging to 3-D problems. Relevance of present analysis is demonstrated through solutions in a simple text book problem of simply supported square plate under doubly sinusoidal load.

Lamination-dependent shear deformation models for bending of angle-ply laminated plates

Lamination-dependent shear corrective terms in the analysis of bending of laminated plates are derived from a priori assumed linear thicknesswise distributions for gradients of transverse shear stresses by using CLPT inplane stresses in the two in-plane equilibrium equations of elasticity in each ply. In the development of a general model for angle-ply laminated plates, special cases like cylindrical bending of laminates in either direction, symmetric laminates, cross-ply laminates, antisymmetric angle-ply laminates, homogeneous plates are taken into consideration. Adding these corrective terms to the assumed displacements in (i) Classical Laminate Plate Theory (CLPT) and (ii) Classical Laminate Shear Deformation Theory (CLSDT), two new refined lamination-dependent shear deformation models are developed. Closed form solutions from these models are obtained for antisymmetric angle-ply laminates under sinusoidal load for a type of simply supported boundary conditions. Results obtained from the present models and also from Ren's model (1987) are compared with each other.

Vibration analysis of laminated plates using a refined shear deformation-theory

A previously published discrete-layer shear deformation theory is used to analyze free vibration of laminated plates. The theory includes the assumption that the transverse shear strains across any two layers are linearly dependent on each other. The theory has the same dependent variables as first order shear deformation theory, but the set of governing differential equations is of twelfth order. No shear correction factors are required. Free vibration of simply supported symmetric and antisymmetric cross-ply plates is calculated. The numerical results are in good agreement with those from three-dimensional elasticity theory.

Static analysis of functionally graded cylindrical and conical shells or panels using the generalized unconstrained third order theory coupled with the stress recovery

A 2D Unconstrained Third Order Shear Deformation Theory (UTSDT) is presented for the evaluation of tangential and normal stresses in moderately thick functionally graded conical and cylindrical shells subjected to mechanical loadings. Several types of graded materials are investigated. The functionally graded material consists of ceramic and metallic constituents. A four parameter power law function is used. The UTSDT allows the presence of a finite transverse shear stress at the top and bottom surfaces of the graded shell. In addition, the initial curvature effect included in the formulation leads to the generalization of the present theory (GUTSDT). The Generalized Differential Quadrature (GDQ) method is used to discretize the derivatives in the governing equations, the external boundary conditions and the compatibility conditions. Transverse and normal stresses are also calculated by integrating the three dimensional equations of equilibrium in the thickness direction. In this way, the six components of the stress tensor at a point of the conical or cylindrical shell or panel can be given. The initial curvature effect and the role of the power law functions are shown for a wide range of functionally conical and cylindrical shells under various loading and boundary conditions. Finally, numerical examples of the available literature are worked out.

Application of the method of initial functions for the analysis of composite laminated plates

The method of initial functions has been applied for deriving higher order theories for cross-ply laminated composite thick rectangular plates. The equations of three-dimensional elasticity have been used. No a priori assumptions regarding the distribution of stresses or displacements are needed. Numerical solutions of the governing equations have been presented for simply supported edges and the results are compared with available ones.

Post-fracture behavior of laminated plates after human impact test

For safety barriers the load bearing capacity of the glass when subjected to the soft body impact should be verified. The soft body pendulum test became a testing standard to classify safety glass plates. The classification of the safety glass do not consider the structural behavior when one sheet of a laminated glass is broken; in situations when the replacement of the plate could not be very urgent, structural behavior should be evaluated. The main objective of this paper is to present the structural behavior o laminated glass plates, though modal test and human impact test, including the post fracture behavior for the laminated cases. A god reproducibility and repeatability is obtained. Two main aspects of the structural behavior can be observed: the increment of the rupture load for laminated plates after the failure of the first sheet, and some similarities with a tempered monolithic behavior of equivalent thickness.

Semi-analytical solution for nonlinear vibration of laminated FGM plates with geometric imperfections

This paper investigates the nonlinear vibration of imperfect shear deformable laminated rectangular plates comprising a homogeneous substrate and two layers of functionally graded materials (FGMs). A theoretical formulation based on Reddy's higher-order shear deformation plate theory is presented in terms of deflection, mid-plane rotations, and the stress function. A semi-analytical method, which makes use of the one-dimensional differential quadrature method, the Galerkin technique, and an iteration process, is used to obtain the vibration frequencies for plates with various boundary conditions. Material properties are assumed to be temperature-dependent. Special attention is given to the effects of sine type imperfection, localized imperfection, and global imperfection on linear and nonlinear vibration behavior. Numerical results are presented in both dimensionless tabular and graphical forms for laminated plates with graded silicon nitride/stainless steel layers. It is shown that the vibration frequencies are very much dependent on the vibration amplitude and the imperfection mode and its magnitude. While most of the imperfect laminated plates show the well-known hard-spring vibration, those with free edges can display soft-spring vibration behavior at certain imperfection levels. The influences of material composition, temperature-dependence of material properties and side-to-thickness ratio are also discussed. (C) 2004 Elsevier Ltd. All rights reserved.

Iterative modelling for stress analysis of composite laminates

The importance of interlaminar stresses has prompted a fresh look at the theory of laminated plates. An important feature in modelling such laminates is the need to provide for continuity of some strains and stresses, while at the same time allowing for the discontinuities in the others. A new modelling possibility is examined in this paper. The procedure allows for discontinuities in the in-plane stresses and transverse strains and continuity in the in-plane strains and transverse stresses. This theory is in the form of a heirarchy of formulations each representing an iterative step. Application of the theory is illustrated by considering the example of an infinite laminated strip subjected to sinusoidal loading.

An exact analysis for vibration of simply-supported homogeneous and laminated thick rectangular plates

A three-dimensional linear, small deformation theory of elasticity solution by the direct method is developed for the free vibration of simply-supported, homogeneous, isotropic, thick rectangular plates. The solution is exact and involves determining a triply infinite sequence of eigenvalues from a doubly infinite set of closed form transcendental equations. As no restrictions are placed on the thickness variation of stresses or displacements, this formulation yields a triply infinite spectrum of frequencies, instead of only one doubly infinite spectrum by thin plate theory and three doubly infinite spectra by Mindlin's thick plate theory. Further, the present analysis yields symmetric thickness modes which neither of the approximate theories can identify. Some numerical results from the two approximate theories are compared with those from the present solution and some important conclusions regarding the effect of the assumptions made in the approximate theories are drawn. The thickness variations of stresses and displacements are also discussed. The analysis is readily extended for laminated plates of isotropic materials. Numerical results are also given for three-ply laminates, and are used to assess the accuracy of thin plate theory predictions for laminates. Extension to general lateral surface conditions and forced vibrations is indicated.

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.

Dynamic-response of a laminated plate with friction damping

Recent advances of finite elements for laminated composite plates

A large amount of finite elements have been developed for finite element analysis of laminated composite plates. The laminated plate theories are reviewed and summarized in this paper. The focus of this review is on the recently developed laminated finite elements since 1990. The 2-D triangular and quadrilateral displacement-based and mixed/hybrid-based finite element models, which were developed based on the first-order shear deformation theories, the higher-order shear deformation theories, the zig-zag theories and the global-local higher-order deformation theories, and the layer-wise laminated plate theories are reviewed in this paper and also their related patents. Finally, some points on the development of the laminated finite elements are summarized.