Nonlinear Oscillations and Stability of Parametrically Excited Cylindrical Shells

Based on Donnell shallow shell equations, the nonlinear vibrations and dynamic instability of axially loaded circular cylindrical shells under both static and harmonic forces is theoretically analyzed. First the problem is reduced to a finite degree-of-freedom one by using the Galerkin method; then the resulting set of coupled nonlinear ordinary differential equations of motion are solved by the Runge–Kutta method. To study the nonlinear behavior of the shell, several numerical strategies were used to obtain Poincaré maps, Lyapunov exponents, stable and unstable fixed points, bifurcation diagrams, and basins of attraction. Particular attention is paid to two dynamic instability phenomena that may arise under these loading conditions: parametric excitation of flexural modes and escape from the pre-buckling potential well. Calculations are carried out for the principal and secondary instability regions associated with the lowest natural frequency of the shell. Special attention is given to the determination of the instability boundaries in control space and the identification of the bifurcational events connected with these boundaries. The results clarify the importance of modal coupling in the post-buckling solution and the strong role of nonlinearities on the dynamics of cylindrical shells.     

Free vibration analysis of a hung clamped-free cylindrical shell partially submerged in fluid

The free flexural vibration of a hung clamped-free cylindrical shell partially submerged in a fluid is investigated. The fluid is assumed to be inviscid and irrotational. The cylindrical shell is modelled by using the Rayleigh-Ritz method based on Sanders’ shell theory. The kinetic energy of the fluid is derived by solving the boundary-value problem related to the fluid motion. The natural vibration characteristics of the partially submerged cylindrical shell are discussed with respect to the added virtual mass approach. In this study, the nondimensionalized added virtual mass incremental factor for the partially submerged finite shell is derived. This factor can be readily used to estimate the change in the natural frequency of the shell due to the presence of a fluid.     

Nonlinear buckling and postbuckling of heated functionally graded cylindrical shells under combined axial compression and radial pressure

The nonlinear large deflection theory of cylindrical shells is extended to discuss nonlinear buckling and postbuckling behaviors of functionally graded (FG) cylindrical shells which are synchronously subjected to axial compression and lateral loads. In this analysis, the non-linear strain-displacement relations of large deformation and the Ritz energy method are used. The material properties of the shells vary smoothly through the shell thickness according to a power law distribution of the volume fraction of the constituent materials. Meanwhile, by taking the temperature-dependent material properties into account, various effects of external thermal environment are also investigated. The non-linear critical condition is found by defining the possible lowest point of external force. Numerical results show various effects of the inhomogeneous parameter, dimensional parameters and external thermal environments on non-linear buckling behaviors of combine-loaded FG cylindrical shells. In addition, the postbuckling equilibrium paths are also plotted for axially loaded pre-pressured FG cylindrical shells and there is an interesting mode jump exhibited.     

径向余弦冲击的长圆柱壳弹性脉冲屈曲

根据已有的实验现象假定变形模态，采用Ｌａｇｒａｎｇｅ方法分析了单位长度薄圆环在余弦冲击载荷作用下的弹性脉冲动屈曲。导出了动力屈曲方程组，借助数值方法求解方程，并与已有的计算进行了比较，得出结论     

Coupled thermoelasticity of functionally graded cylindrical shells

The coupled thermoelastic response of a functionally graded circular cylindrical shell is studied. The coupled thermoelastic and the energy equations are simultaneously solved for a functionally graded axisymmetric cylindrical shell subjected to thermal shock load. A second-order shear deformation shell theory that accounts for the transverse shear strains and rotations is considered. Including the thermo-mechanical coupling and rotary inertia, a Galerkin finite element formulation in space domain and the Laplace transform in time domain are used to formulate the problem. The inverse Laplace transform is obtained using a numerical algorithm. The shell is graded through the thickness assuming a volume fraction of metal and ceramic, using a power law distribution. The results are validated with the known data in the literature.     

圆柱壳撞水时的弹塑性动力屈曲研究

本文研究圆柱壳在流固中冲击载荷下的弹塑性动力曲问题。建立了液体－气体－固体三相的数学模型。其中结构部分控制方程由弹塑性力学中关于加速度的最小原理获得，本构关系采用增量理论，液体假设为不可压缩，空气层认为最等熵压缩。分别讨论了不同冲击高度时压力变化规律，屈曲对载荷的影响。屈曲沿壳长的分布及发展规律等。     

Nonlinear vibrations of a rotating shaft with broadband random variations of internal damping

A simple Jeffcott rotor is considered with broadband temporal random variations of internal damping which are described using the theory of Markov processes. Transverse response of the rotor with stiffening nonlinearity either in external damping or in restoring force is studied by stochastic averaging method. This method reduces the problems to stochastic differential equations (SDEs) for which analytical solutions are obtained for the Fokker–Planck–Kolmogorov (FPK) equations for stationary probability density functions (PDFs) of the squared whirl radius of the shaft. These PDFs do exist beyond the dynamic instability threshold and they correspond to forward whirl of the rotor. At rotation speeds just slightly above the instability threshold, the response PDF has integrable singularity at zero which corresponds to intermittency in the response.     

Analytical and experimental modal analysis of nonuniformly ring-stiffened cylindrical shells

In this research, the free vibration analysis of cylindrical shells with circumferential stiffeners, i.e., rings with nonuniform stiffener eccentricity and unequal stiffener spacing, is investigated using analytical and experimental methods. The Ritz method is applied in analytical solution, while stiffeners are treated as discrete elements. The polynomial functions are used for Ritz functions. The effects of nonuniformity of stiffener distribution on natural frequencies are considered for free–free boundary conditions. Results show that, at constant stiffener mass, significant increments in natural frequencies can be achieved using nonuniform stiffener distribution. In experimental method, modal testing is performed to obtain modal parameters, including natural frequencies, mode shapes, and damping in each mode. Analytical results are compared with experimental ones, showing good agreement. Because of insufficient experimental modal data for nonuniform stiffener distribution, the results of modal testing obtained in this study could be a useful reference for validating the accuracy of other analytical and numerical methods for free vibration analysis.     

Stability of shells of revolution made of a damageable nonlinear elastic material

The problem of bifurcation instability of shells of revolution made of a nonlinear elastic material progressively damaged under loading is formulated and solved __________ Translated from Prikladnaya Mekhanika, Vol. 43, No. 5, pp. 35–45, May 2007.     

Plate/shell structure topology optimization of orthotropic material for buckling problem based on independent continuous topological variables

The purpose of the present work is to study the buckling problem with plate/shell topology optimiza-tion of orthotropic material.A model of buckling topology optimization is established based on the independent,con-tinuous, and mapping method, which considers structural mass as objective and buckling critical loads as constraints. Firstly, composite exponential function (CEF) and power function(PF)as filter functions are introduced to recognize the element mass,the element stiffness matrix,and the ele-ment geometric stiffness matrix.The filter functions of the orthotropic material stiffness are deduced. Then these fil-ter functions are put into buckling topology optimization of a differential equation to analyze the design sensitiv-ity.Furthermore,the buckling constraints are approximately expressed as explicit functions with respect to the design vari-ables based on the first-order Taylor expansion.The objective function is standardized based on the second-order Taylor expansion. Therefore,the optimization model is translated into a quadratic program.Finally,the dual sequence quadratic programming(DSQP)algorithm and the global convergence method of moving asymptotes algorithm with two different filter functions(CEF and PF)are applied to solve the opti-mal model.Three numerical results show that DSQP&CEF has the best performance in the view of structural mass and discretion.     

Stability of cylindrical shells made of a damageable physically nonlinear material

The bifurcation-instability problem for cylindrical shells made of physically nonlinear materials progressively damaged with increasing load is formulated and solved __________ Translated from Prikladnaya Mekhanika, Vol. 42, No. 10, pp. 56–66, October 2006.     

Buckling analysis of multi-walled carbon nanotubes: a continuum model accounting for van der Waals interaction

Explicit formulas are derived for the van der Waals (vdW) interaction between any two layers of a multi-walled carbon nanotube (CNT). Based on the derived formulas, an efficient algorithm is established for the buckling analysis of multi-walled CNTs, in which individual tubes are modeled as a continuum cylindrical shell. The explicit expressions are also derived for the buckling of double-walled CNTs. In previous studies by Ru (J. Appl. Phys. 87 (2000b) 7227) and Wang et al. (Int. J. Solids Struct. 40 (2003) 3893), only the vdW interaction between adjacent two layers was considered and the vdW interaction between the other two layers was neglected. The vdW interaction coefficient was treated as a constant that was not dependent on the radii of the tubes. However, the formulas derived herein reveal that the vdW interaction coefficients are dependent on the change of interlayer spacing and the radii of the tubes. With the increase of radii, the coefficients approach constants, and the constants between two adjacent layers are about 10% higher than those reported by Wang et al. (Int. J. Solids. Struct. 40 (2003) 3893). In addition, the numerical results show that the vdW interaction will lead to a higher critical buckling load in multi-walled CNTs. The effect of the tube radius on the critical buckling load of a multi-walled CNT is also examined.     

正交各向异性功能梯度圆柱壳的自由振动

利用从弹性力学的三维基本方程导得的状态方程，结合层合近似模型，求得了柱型正交各向异性功能梯度圆柱壳的自由振动频率，给出了数值算例，并与前人结果作了比较，文末讨论了相关参数的影响。     

Free vibration of circular cylindrical shell with constrained layer damping

Free vibration characteristics of circular cylindrical shell with passive constrained layer damping (PCLD) are presented. Wave propagation approach rather than finite element method, transfer matrix method, and Rayleigh-Ritz method is used to solve the problem of vibration of PCLD circular cylindrical shell under a simply supported boundary condition at two ends. The governing equations of motion for the orthotropic cylindrical shell with PCLD are derived on the base of Sanders’ thin shell theory. Numerical results show that the present method is more effective in comparison with other methods. The effects of the thickness of viscoelastic core and constrained layer, the elastic modulus ratio of orthotropic constrained layer, the complex shear modulus of viscoelastic core on frequency parameter, and the loss factor are discussed.     

Triangular grid method for stress-wave propagation in 2-D orthotropic materials

A triangular grid method is presented to calculate propagation problems of elastic stress waves in 2-D orthotropic materials. This method is based on the dynamic equilibrium equations of the computational cells formed among the auxiliary triangular grids. The solution is obtained by calculating alternately the nodal displacements and the central point stresses of the spatial grids. The numerical results are compared with the corresponding solutions of the finite element method. Comparisons show that the triangular grid method yields a higher calculational speed than the finite element method. The stress concentrations are investigated from wave-field analyses when the stress wave propagates within an orthotropic plate with a hole. Finally, the presented numerical method is used to study the features of wave propagation and diffraction in a square orthotropic plate with a hole when an impact load is applied to the top of the plate.This work was supported by National Natural Science Foundation of China (Nos. 10025212 and 10232040) and Natural Science Foundation of Liaoning province (No. 20021070).     

FREE VIBRATION OF PIEZOELECTRIC CYLINDRICAL SHELLS

Three displacement functions are introduced to represent each mechanical displacementaccording to the 3-D theory in this paper,By expanding the displacement functions and the electric po-tential in orthogonal series,the free vibration equation of piezoelectric cylindrical shells satisfying SS3boundary conditions can be obtained.The equation was solved by utilizing Bessel functions with com-plex arguments.Results are presented graphically as well as in table,and compared with those of otherreferences.Some frequencies that were missing in Ref.[9]are discovered.     

Large amplitude free flexural vibration of rings using finite element approach

Here, the large amplitude free flexural vibrations of isotropic/laminated orthotropic rings are investigated, using a shear flexible curved beam element based on field consistency principle. A laminated refined beam theory is introduced for developing the element, which satisfies the interface transverse shear stress and displacement continuity, and has a vanishing shear stress on the inner and outer surfaces of the beam. The formulation includes in-plane and rotary inertia effects, and the non-linearity due to the finite deformation of the ring. The governing equations obtained using Lagrange's equations of motion are solved through the direct integration technique. Amplitude-frequency relationships evaluated from the dynamic response history are examined. Detailed numerical results are presented considering various parameters such as radius-to-thickness ratio, circumferential wave number and ovality for isotropic and laminated orthotropic rings. The nature and degree of the participation of various modes in non-linear asymmetric vibration of oval ring brought out through the present study are useful for accurate modelling of the closed non-circular structures.     

Dynamics of the geometrically non-linear analysis of anisotropic laminated cylindrical shells

A general approach, based on shearable shell theory, to predict the influence of geometric non-linearities on the natural frequencies of an elastic anisotropic laminated cylindrical shell incorporating large displacements and rotations is presented in this paper. The effects of shear deformations and rotary inertia are taken into account in the equations of motion. The hybrid finite element approach and shearable shell theory are used to determine the shape function matrix. The analytical solution is divided into two parts. In part one, the displacement functions are obtained by the exact solution of the equilibrium equations of a cylindrical shell based on shearable shell theory instead of the usually used and more arbitrary interpolating polynomials. The mass and linear stiffness matrices are derived by exact analytical integration. In part two, the modal coefficients are obtained, using Green's exact strain-displacement relations, for these displacement functions. The second- and third-order non-linear stiffness matrices are then calculated by precise analytical integration and superimposed on the linear part of equations to establish the non-linear modal equations. Comparison with available results is satisfactorily good.     

The energy transfer mechanism of a perturbed solid-body rotation flow in a rotating pipe

Three-dimensional direct numerical simulations of a solid-body rotation superposed on a uniform axial flow entering a rotating constant-area pipe of finite length are presented. Steady in time profiles of the radial, axial, and circumferential velocities are imposed at the pipe inlet. Convective boundary conditions are imposed at the pipe outlet.The Wang and Rusak(Phys. Fluids 8:1007–1016, 1996.doi:10.1063/1.86882) axisymmetric instability mechanism is retrieved at certain operational conditions in terms of incoming flow swirl levels and the Reynolds number. However, at other operational conditions there exists a dominant,three-dimensional spiral type of instability mode that is consistent with the linear stability theory of Wang et al.(J. Fluid Mech. 797: 284–321, 2016). The growth of this mode leads to a spiral type of flow roll-up that subsequently nonlinearly saturates on a large amplitude rotating spiral wave. The energy transfer mechanism between the bulk of the flow and the perturbations is studied by the Reynolds-Orr equation. The production or loss of the perturbation kinetic energy is combined of three components: the viscous loss, the convective loss at the pipe outlet, and the gain of energy at the outlet through the work done by the pressure perturbation. Theenergy transfer in the nonlinear stage is shown to be a natural extension of the linear stage with a nonlinear saturated process.