Large rotations in first-order shear deformation FE analysis of laminated shells

The paper deals with the geometrically non-linear analysis of laminated composite beams, plates and shells in the framework of the first-order transverse shear deformation (FOSD) theory. A central point of the present paper is the discussion of the relevance of five- and six-parameter variants, respectively, of the FOSD hypothesis for large rotation plate and shell problems. In particular, it is shown that the assumption of constant through-thickness distribution of the transverse normal displacements is acceptable only for small and moderate rotation problems. Implications inherent in this assumption that are incompatible with large rotations are discussed from the point of view of the transverse normal strain-displacement relations as well as in the light of an enhanced, accurate large rotation formulation based on the use of Euler angles. The latter one is implemented as an updating process within a Total Lagrangian formulation of the six-parameter FOSD large rotation plate and shell theory. Numerical solutions are obtained by using isoparametric eight-node Serendipity-type shell finite elements with reduced integration. The Riks-Wempner-Ramm arc-length control method is used to trace primary and secondary equilibrium paths in the pre- and post-buckling range of deformation. A number of sample problems of non-linear, large rotation response of composite laminated plate and shell structures are presented including symmetric and asymmetric snap-through and snap-back problems.     

板壳问题的三维无网格伽辽金直接分析法

对于板壳问题,共有三种数值模拟方案:线性或非线性的板壳理论、退化连续体方案和直接三维连续体方案。无网格法近似函数可具有C1甚至更高的连续性,便于在K irchhoff-Love理论中应用。但当各种无网格法用于M ind lin-R e issner板理论时,会遇到数值锁死的困扰。对比之下,三维连续体方案是最简单,最精确但并不常用的一种方案。无网格法近似函数具有高度光滑性,在板壳的厚度方向仅布置2~5层点就可以很好地捕捉此方向场的梯度,同时还可以在一定参数范围内避免剪切和体积锁死,在处理复杂本构关系、非线性板壳等问题中更是具有很大优势。本文采用无网格伽辽金法(EFG)和三维连续体方案分析了线性板壳问题,与有限单元法做了对比,并讨论了数值锁死等问题。     

Natural vibrations of a cantilever thin elastic orthotropic cylindrical shell

The natural vibrations of a cantilever thin elastic orthotropic circular cylindrical shell are studied. Dispersion equations for the determination of possible natural frequencies of cantilever closed shells and open shells with Navier hinged boundary conditions at the longitudinal edges are derived from the classical dynamic theory of orthotropic cylindrical shells. It is proved that there are asymptotic relationships between these problems and the problems for a cantilever orthotropic strip plate and for a cantilever rectangular plate and the eigenvalue problem for a semi-infinite closed orthotropic cylindrical shell with free end and for the same but open shell with Navier hinged boundary conditions at the longitudinal edges. A procedure to identify types of vibrations is presented. Orthotropic cylindrical shells with different radii and lengths are used as an example to find approximate values of the dimensionless natural frequency and damping factor for vibration modes __________ Translated from Prikladnaya Mekhanika, Vol. 44, No. 5, pp. 68–91, May 2008.     

Vibrations of multilayer plates under the effect ofimpulse loads. Three-dimensional theory

A new analytical–numerical approach to investigation of the response of multilayerplates to impulse loading is described in this paper. The plates behaviour is described by theequations of the three-dimensional elasticity theory. According to the approach being proposed,the sought for functions included in the system of equations and the boundary and initialconditions are presented as Fourier series expansions in the tangential directions. The derivativesof these functions in the transverse direction are replaced by their finite-difference presentations.As a result of such transforms, the problem of vibration of a multilayer plate is reduced tointegration of a system of ordinary differential equations with constant coefficients. Integration isperformed by expansion into the Taylors series. The possibilities of the approach proposed andthe validity of results obtained is illustrated by several examples of calculating vibration processesand the processes of propagation of elastic waves. A comparison of the results obtained on thebasis of other approaches has been performed.     

Direct approach-based analysis of plates composed of functionally graded materials

The classical plate theory can be applied to thin plates made of classical materials like steel. The first theory allowing the analysis of such plates was elaborated by Kirchhoff. But this approach was connected with various limitations (e.g., constant material properties in the thickness direction). In addition, some mathematical inconsistencies like the order of the governing equation and the number of boundary conditions exist. During the last century many suggestions for improvements of the classical plate theory were made. The engineering direction of improvements was ruled by applications (e.g., the use of laminates or sandwiches as the plate material), and so new hypotheses for the derivation of the governing equations were introduced. In addition, some mathematical approaches like power series expansions or asymptotic integration techniques were applied. A conceptional different direction is connected with the direct approach in the plate theory. This paper presents the extension of Zhilin’s direct approach to plates made of functionally graded materials. The second author was supported by DFG grant 436RUS17/21/07.     

A variational differential quadrature solution to finite deformation problems of hyperelastic shell-type structures: a two-point formulation in Cartesian coordinates

A new numerical approach is presented to compute the large deformations of shell-type structures made of the Saint Venant-Kirchhoff and Neo-Hookean materials based on the seven-parameter shell theory. A work conjugate pair of the first Piola Kirchhoff stress tensor and deformation gradient tensor is considered for the stress and strain measures in the paper. Through introducing the displacement vector, the deformation gradient, and the stress tensor in the Cartesian coordinate system and by means of the chain rule for taking derivative of tensors, the difficulties in using the curvilinear coordinate system are bypassed. The variational differential quadrature (VDQ) method as a pointwise numerical method is also used to discretize the weak form of the governing equations. Being locking-free, the simple implementation, computational efficiency, and fast convergence rate are the main features of the proposed numerical approach. Some well-known benchmark problems are solved to assess the approach. The results indicate that it is capable of addressing the large deformation problems of elastic and hyperelastic shell-type structures efficiently.     

A polar theory for vibrations of thin elastic shells

In relation to a polar continuum, this paper presents a 2-D shear deformable theory for the high frequency vibrations of a thin elastic shell. To begin with, the 3-D fundamental equations of the micropolar elastic continuum are expressed as the Euler–Lagrange equations of a unified variational principle. Next, the kinematic variables of the shell are represented by the power series expansions in its thickness coordinate, and then, they are used to establish the 2-D theory by means of the variational principle. The 2-D theory is derived in invariant variational and differential forms and governs all the types of vibrations of the functionally graded micropolar shell. Lastly, the uniqueness is investigated in solutions of the initial mixed boundary value problems defined by the 2-D theory, and some of special cases are indicated in the theory.     

梁板壳的几何大变形--从近似的非线性理论到有限变形理论

对梁板壳的线性理论、近似几何非线性理论与有限变形理论作了比较,介绍了有限转动理论,指出了应用有限变形理论求解梁板壳的大变形问题的高效率、高精度的巨大优越性。     

Axial compression of circular cylindrical bar with coaxial subsurface cylindrical crack of finite length

A fracture stability of a circular cylindrical bar with a coaxial surface cylindrical crack subjected to an axial compression is considered. A state of subcritical initial strain is assumed. A non-classical fracture criterion is based on a local stability loss near the defect. The theory of integral Fourier transforms and series expansions are used to reduce these problems to a system of paired integral equations and then to a system of linear algebraic equations with respect to the contraction parameter.     

Large deflection analysis of unsymmetrically laminated composite plates: analytical-numerical type approach

Large deflection analysis of laminated composite plates is considered. The Galerkin method along with Newton-Raphson method is applied to large deflection analysis of laminated composite plates with various edge conditions. The von Kármán plate theory is utilized and the governing differential equations are solved by choosing suitable polynomials as trial functions to approximate the plate displacement functions. The solutions are compared to that of Dynamic Relaxation and finite elements. A very close agreement has been observed with these approximating methods. In the solution process, analytical computation has been done wherever it is possible, and analytical-numerical type approach has been made for all problems.     

Integral formulation for a circular cylindrical cavity in infinite solid and a finite length coaxial cylindrical crack compressed axially

A method for solving problems of fracture of an infinite solid with a circular cylindrical cavity and a coaxial cylindrical crack near the surface under an uniform axial compression is proposed using a non-classical criterial approach associated with a mechanism of a local stability loss near the defect. The theory of integral Fourier transforms and series expansions are used to reduce these problems to a system of paired integral equations and then to a system of linear algebraic equations with respect to the contraction parameter.     

Solution of contact problems on the basis of a refined theory of plates and shells

Solutions of contact mixed boundary-value problems for a plate and for a cylindrical shell are given. These solutions are obtained with the use of equations for shells constructed by expanding solutions of elasticity theory equations with respect to the Legendre polynomials. Results of numerical simulations of the stress state in the vicinity of points with changing conditions on the frontal faces of the shell are presented. The results obtained are compared with analytical solutions of elasticity theory problems and with solutions obtained on the basis of the classical equations of the shell theory. __________ Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 49, No. 5, pp. 169–176, September–October, 2008.     

Non-linear vibrations of shallow circular cylindrical panels with complex geometry. Meshless discretization with the R-functions method

Geometrically non-linear forced vibrations of a shallow circular cylindrical panel with a complex shape, clamped at the edges and subjected to a radial harmonic excitation in the spectral neighborhood of the fundamental mode, are investigated. Both Donnell and the Sanders–Koiter non-linear shell theories retaining in-plane inertia are used to calculate the elastic strain energy. The discrete model of the non-linear vibrations is build using the meshfree technique based on classic approximate functions and the R-function theory, which allows for constructing the sequences of admissible functions that satisfy given boundary conditions in domains with complex geometries; Chebyshev orthogonal polynomials are used to expand shell displacements. A two-step approach is implemented in order to solve the problem: first a linear analysis is conducted to identify natural frequencies and corresponding natural modes to be used in the second step as a basis for expanding the non-linear displacements. Lagrange approach is applied to obtain a system of ordinary differential equations on both steps. Different multimodal expansions, having from 15 up to 35 generalized coordinates associated with natural modes, are used to study the convergence of the solution. The pseudo-arclength continuation method and bifurcation analysis are applied to study non-linear equations of motion. Numerical responses are obtained in the spectral neighborhood of the lowest natural frequency; results are compared to those available in the literature. Internal resonances are also detected and discussed.     

New cases of propagation of singularities along characteristic boundaries for model problems of shell theory

We continue a previous work [1] on propagation of singularities for model problems of thin shell theory in the parabolic and hyperbolic cases. The singularities along the characteristic boundaries are considered using extensions of the solutions out of the domain, adapted to either free or fixed boundaries. The corresponding transport equations are given except for the case of a characteristic fixed boundary for a hyperbolic shell, where the phenomenon is non local, but depends on the whole domain.     

Active sloshing control in a smart flexible cylindrical floating roof tank

An exact three dimensional fully-coupled hydro-elastic analysis for transient liquid sloshing in a partially-filled vertically-standing flexible circular cylindrical shell container fitted with a freely floating smart piezo-sandwich thin elastic circular plate is presented. The problem formulation is based on the linear water wave theory, the classical (Kirchhoff/Sanders) thin plate and shell models, Maxwell's equations of electrodynamics, Stokes’ transformation, and eigen-function expansions in cylindrical coordinates. The control action is achieved by combined volume displacement and volume velocity feedbacks (VDF, VVF) implemented in a second order active damping (AD) compensator via two competent evolutionary heuristic optimization techniques that systematically tune the controller gain parameters while constraining the floating panel displacement and control voltage. The uncontrolled and controlled transient responses of the coupled hydro-elastic system under various external disturbances (i.e., a harmonic base excitation, a real seismic event, a severe launch vehicle liftoff event, and a distributed impulsive transverse load on the floating panel) are calculated by means of Durbin's numerical inverse Laplace transform scheme. Moreover, the free vibration characteristics of the coupled fluid/structure interaction (FSI) system are briefly studied. The superior performance of the proposed active floating roof control configuration in effective suppression of the key hydro-elastic parameters (panel displacement, and shell displacements/stresses) is demonstrated. It is also found that, in the current FSI control problem, the Multi-objective Particle Swarm Optimization (MOPSO)-based ADC outperforms the Non-dominated Sorting Genetic Algorithm (NSGA-II)-based method, in terms of convergence rate and computational effort. Limiting cases are examined and the precision of results is verified by comparisons with the existing data as well as with the results produced by a commercial finite element package.     

Stress-strain state of non-thin plates and shells. Generalized theory (survey)

Conclusions Thus, this part of our survey has presented the main approaches that have been taken to the construction of two-dimensional (in terms of the space coordinates) equations of a generalized theory of plates and shells. The solutions of these equations represent a certain approximation of the solution of the initial three-dimensional problem. They are based on expansion of the sought functions into Fourier series in Legendre polynomials of the thickness coordinate. Studies completed on the basis of the given variants of plate and shell theory were systematized and analyzed. In terms of the method of its construction, the theory involves a regular process of replacing the solution of the three-dimensional problem by the solution (or sequence of solutions) of two-dimensional boundary-value problems or initial-boundary-value problems. Numerical results illustrating the convergence of the successive approximation were presented. It should be noted that to make comparison with the results of classical or applied theories, several of the studies cited here presented solutions of problems for thin plates and shells with allowance only for the initial terms of expansions of the stress and displacement components into base functions (Legendre polynomials).S. P. Timoshenko Institute of Mechanics, Ukrainian Academy of Sciences, Kiev. Translated from Prikladnaya Mekhanika, Vol. 29, No. 11, pp. 3–34, November, 1993.     

Multimode vibration control using several piezoelectric transducers shunted with a multiterminal network

In this paper a new approach is presented to reduce vibrations for one- and two-dimensional mechanical structures, as beam or thin plates, by means of several piezoelectric transducers shunted with a proper electric network system. The governing equations of the whole system are coupled to each other through the direct and converse piezoelectric effect. More in detail, the mechanical equations are expressed in accordance with the modal theory considering n vibration modes and the electrical equations reduce to the one-dimensional charge equation of electrostatics for each of n considered piezoelectric transducers. In this electromechanical system, a shunting electric device forms an electric subsystem working as multi degrees of freedom (dof’s) damped vibration absorber for the mechanical subsystem. Herein, it is introduced a proper transformation of the electric coordinates in order to approximate the governing equations for the whole shunted system with n uncoupled, single mode piezoelectric shunting systems that can be readily damped by the methods reported in literature. A further numerical optimisation problem on the spatial distribution of the piezoelectric elements allows to achieve a better performance. Numerical case studies of two relevant systems, a double clamped beam and a fully clamped plate, allow to take into account issues relative to the proposed approach. Laboratory experiments carried out in real time on a beam clamped at both ends consent to validate the proposed technique.     

Three dimensional degenerated shell theory and assumed strain shell elements

In this paper Reissner-Mindlin plate theory is extended to cater for curved shell structures. It can be considered as Reissner-Mindlin type shell theory. From this theory, the C(O) continuity formulation of shell elements of taking account the transverse shear deformation could be derived directly. These degenerated shell elements have been widely employed. To overcome the locking of shear and membrane and avoid zero energy modes the author proposed the formulation of the new elements with assumed strains. A wide range of numerical tests was conducted and the results illustrate that the assumed strain elements possess high accuracy and good performance.     

Asymptotic splitting in the three-dimensional problem of elasticity for non-homogeneous piezoelectric plates

A novel asymptotic approach to the theory of non-homogeneous anisotropic plates is suggested. For the problem of linear static deformations we consider solutions, which are slowly varying in the plane of the plate in comparison to the thickness direction. A small parameter is introduced in the general equations of the theory of elasticity. According to the procedure of asymptotic splitting, the principal terms of the series expansion of the solution are determined from the conditions of solvability for the minor terms. Three-dimensional conditions of compatibility make the analysis more efficient and straightforward. We obtain the system of equations of classical Kirchhoff's plate theory, including the balance equations, compatibility conditions, elastic relations and kinematic relations between the displacements and strain measures. Subsequent analysis of the edge layer near the contour of the plate is required in order to satisfy the remaining boundary conditions of the three-dimensional problem. Matching of the asymptotic expansions of the solution in the edge layer and inside the domain provides four classical plate boundary conditions. Additional effects, like electromechanical coupling for piezoelectric plates, can easily be incorporated into the model due to the modular structure of the analysis. The results of the paper constitute a sound basis to the equations of the theory of classical plates with piezoelectric effects, and provide a trustworthy algorithm for computation of the stressed state in the three-dimensional problem. Numerical and analytical studies of a sample electromechanical problem demonstrate the asymptotic nature of the present theory.     

三维退化的壳理论和假定应变的壳单元

本文将Reissner-Mindlin板理论推广到空间曲壳结构,可称为Reissner-Mindlin型壳理论。从这种理论出发,可直接导出C(0)连续的壳体单元,即考虑横向剪切变型的影向的壳体单元,这种单元在国外已被广泛地采用,为克服这种单元在应用中所出现的剪切和膜的锁制现象同时又防止出现任何零能模式,作者提出了一种采用假定应变的新的壳单元公式,并对这种单元进行了广泛的数值试验,结果表明这种单元具有较高的精度和良好的性能。