Vibration analysis of circular cylindrical double-shell structures under general coupling and end boundary conditions

In this paper, the free and forced vibration analysis of circular cylindrical double-shell structures under arbitrary boundary conditions is presented. This is achieved by employing the improved Fourier series method based on Hamilton’s principle. In the formulation, each displacement component of the cylindrical shells and annular plates is invariantly expanded as the superposition of a standard Fourier series with several supplementary functions introduced to remove the potential discontinuities of the original displacement and its derives at the boundaries. With the introduction of four sets of boundary springs at the coupling interfaces and end boundaries of the shell–plate combination, both elastic and rigid coupling and end boundary conditions can be easily obtained by assigning the stiffnesses of the artificial springs to certain values. The natural frequencies and mode shapes of the structures as well as frequency responses under forced vibration are obtained with the Rayleigh–Ritz procedure. The convergence of the method is validated by comparing the present results with those obtained by the finite element method. Several numerical results including natural frequencies and mode shapes are presented to demonstrate the excellent accuracy and reliability of the current method. Finally, a number of parameter studies concerning various end and coupling boundary conditions, different dimensions of shells and annular plates are also performed.     

Free vibration analysis of functionally graded conical shells and annular plates using the Haar wavelet method

The main aim of this paper is to provide a simple yet efficient solution for the free vibration analysis of functionally graded (FG) conical shells and annular plates. A solution approach based on Haar wavelet is introduced and the first-order shear deformation shell theory is adopted to formulate the theoretical model. The material properties of the shells are assumed to vary continuously in the thickness direction according to general four-parameter power-law distributions in terms of volume fractions of the constituents. The separation of variables is first performed; then Haar wavelet discretization is applied with respect to the axial direction and Fourier series is assumed with respect to the circumferential direction. The constants appearing from the integrating process are determined by boundary conditions, and thus the partial differential equations are transformed into algebraic equations. Then natural frequencies of the FG shells are obtained by solving algebraic equations. Accuracy and reliability of the current method are validated by comparing the present results with the existing solutions. Effects of some geometrical and material parameters on the natural frequencies of shells are discussed and some selected mode shapes are given for illustrative purposes. It’s found that accurate frequencies can be obtained by using a small number of collocation points and boundary conditions can be easily achieved. The advantages of this current solution method consist in its simplicity, fast convergence and excellent accuracy.     

Acoustic analysis of a rectangular cavity with general impedance boundary conditions

A Fourier series method is proposed for the acoustic analysis of a rectangular cavity with impedance boundary conditions arbitrarily specified on any of the walls. The sound pressure is expressed as the combination of a three-dimensional Fourier cosine series and six supplementary two-dimensional expansions introduced to ensure (accelerate) the uniform and absolute convergence (rate) of the series representation in the cavity including the boundary surfaces. The expansion coefficients are determined using the Rayleigh-Ritz method. Since the pressure field is constructed adequately smooth throughout the entire solution domain, the Rayleigh-Ritz solution is mathematically equivalent to what is obtained from a strong formulation based on directly solving the governing equations and the boundary conditions. To unify the treatments of arbitrary nonuniform impedance boundary conditions, the impedance distribution function on each specified surface is invariantly expressed as a double Fourier series expansion so that all the relevant integrals can be calculated analytically. The modal parameters for the acoustic cavity can be simultaneously obtained from solving a standard matrix eigenvalue problem instead of iteratively solving a nonlinear transcendental equation as in the existing methods. Several numerical examples are presented to demonstrate the effectiveness and reliability of the current method for various impedance boundary conditions, including nonuniform impedance distributions.     

A unified approach for predicting sound radiation from baffled rectangular plates with arbitrary boundary conditions

This paper discusses sound radiation from a baffled rectangular plate with each of its edges arbitrarily supported in the form of elastic restraints. The plate displacement function is universally expressed as a 2-D Fourier cosine series supplemented by several 1-D series. The unknown Fourier expansion coefficients are then determined by using the Rayleigh-Ritz procedure. Once the vibration field is solved, the displacement function is further simplified to a single standard 2-D Fourier cosine series in the subsequent acoustic analysis. Thus, the sound radiation from a rectangular plate can always be obtained from the radiation resistance matrix for an invariant set of cosine functions, regardless of its actual dimensions and boundary conditions. Further, this radiation resistance matrix, unlike the traditional ones for modal functions, only needs to be calculated once for all plates with the same aspect ratio. In order to determine the radiation resistance matrix effectively, an analytical formula is derived in the form of a power series of the non-dimensional acoustic wavenumber; the formula is mathematically valid and accurate for any wavenumber. Several numerical examples are presented to validate the formulations and show the effect of the boundary conditions on the radiation behavior of planar sources.     

Vibration characteristics of thin rotating cylindrical shells with various boundary conditions

An analysis is presented for the vibration characteristics of thin rotating cylindrical shells with various boundary conditions by use of Fourier series expansion method. Based on Sanders’ shell equations, the governing equations of motion which take into account the effects of centrifugal and Coriolis forces as well as the initial hoop tension due to rotating are derived. The displacement field is expressed as a product of Fourier series expressions which represents the axial modal displacements and trigonometric functions which represents the circumferential modal displacements. Stokes’ transformation is employed to derive the derivatives of the Fourier series expressions. Then, through the process of formula derivation, an explicit expression of the exact frequency equation can be obtained for a thin rotating cylinder with classical boundary conditions of any type. Once the frequency equation has been determined, the frequencies are calculated numerically. To validate the present analysis, comparisons between the results of the present method and previous studies are performed and very good agreement is achieved. Finally, the method is applied to investigate the vibration characteristics of thin rotating cylindrical shells under various boundaries, and the results are presented.     

Free vibrational characteristics of isotropic coupled cylindrical-conical shells

This paper presents the free vibrational characteristics of isotropic coupled conical-cylindrical shells. The equations of motion for the cylindrical and conical shells are solved using two different methods. A wave solution is used to describe the displacements of the cylindrical shell, while the displacements of the conical sections are solved using a power series solution. Both Donnell-Mushtari and Flügge equations of motion are used and the limitations associated with each thin shell theory are discussed. Natural frequencies are presented for different boundary conditions. The effect of the boundary conditions and the influence of the semi-vertex cone angle are described. The results from the theoretical model presented here are compared with those obtained by previous researchers and from a finite element model.     

Vibrations of rectangular plates reinforced by any number of beams of arbitrary lengths and placement angles

This paper presents an analytical method for the vibration analysis of plates reinforced by any number of beams of arbitrary lengths and placement angles. Both the plate and stiffening beams are generally modeled as three-dimensional (3-D) structures having six displacement components at a point, and the coupling at an interface is generically described by a set of distributed elastic springs. Each of the displacement functions is here invariably expressed as a modified Fourier series, which consists of a standard Fourier cosine series plus several supplementary series/functions used to ensure and improve uniform convergence of the series representation. Unlike most existing techniques, the current method offers a unified solution to the vibration problems for a wide spectrum of stiffened plates, regardless of their boundary conditions, coupling conditions, and reinforcement configurations. Several numerical examples are presented to validate the methodology and demonstrate the effect on modal parameters for a stiffened plate with various boundary conditions, coupling conditions, and reinforcement configurations.     

Free vibrations of cylindrical shells with elastic-support boundary conditions

This paper concerns the free vibrations of cylindrical shells with elastic boundary conditions. Based on the Flügge classical thin shell theory, the equations of motion for the cylindrical shells are solved by the method of wave propagations. The wave numbers are obtained by directly solving an eighth order equation. The elastic-support boundary conditions can be arbitrarily specified in terms of 8 independent sets of distributed springs. All the classical homogeneous boundary conditions can be considered as the special cases when the stiffness for each set of springs is equal to either infinity or zero. The present solutions are validated by the results previously given by other researchers and/or obtained using finite element models. The effects on the frequency parameters of elastic restraints are investigated for shells of different geometrical characteristics.     

Influence of boundary restraint on sound attenuation performance of a duct-membrane silencer

Low frequency noise in duct is a challenge for the traditional passive noise control techniques. Recently, a so-called duct-membrane silencer has attracted much research attention due to its simple configuration and potential application, however, the current studies are merely limited to the cases in which just the classical boundary conditions are considered. Actually, as an important factor affecting the modal characteristics of the membrane, and the existing studies are not enough to fully understand the vibro-acoustic characteristics of such silencer with complicated boundary conditions. Motivated by this, in this paper, the structural–acoustic coupling model of duct-membrane system is established by a modified Fourier series method in combination with Rayleigh–Ritz procedure, in which the transverse elastic boundary restraints are taken into account. Energy principle is formulated for the vibro-acoustic coupling of such duct-membrane silencer to obtain the system matrix equation. Numerical results are then presented to validate the proposed model, and the influence of boundary restraining stiffness on sound attenuation performance is also studied. To the best of authors’ knowledge, this work represents the first time that the elastic boundary restraints have been considered for such duct-membrane silencing system.     

Free vibration of rectangular plates with edges having different degrees of rotational restraint

A numerical method developed by the author has been used as a basis for determining natural frequencies of rectangular plates possessing different degrees of elastic restraints along the edges. The basic functions satisfying the boundary conditions along two opposite edges for such cases have been derived. Comparison of results with others that are available indicates excellent accuracy. Many new results have been presented.     

轴对称体声振耦合的边界子波谱与有限元耦合方法

探讨了子波在Helmholtz积分方程及声振耦合中的应用,在建立了求解轴对称Helmholtz积分方程的子波谱方法的基础上,构造了轴对称子波谱与轴对称有限元的耦合方法,该方法可以处理轴对称问题的任意边界条件.进行了声振耦合问题的模态分析.     

Polynomial versus trigonometric expansions for nonlinear vibrations of circular cylindrical shells with different boundary conditions

Large-amplitude (geometrically nonlinear) forced vibrations of circular cylindrical shells with different boundary conditions are investigated. The Sanders-Koiter nonlinear shell theory, which includes in-plane inertia, is used to calculate the elastic strain energy. The shell displacements (longitudinal, circumferential and radial) are expanded by means of a double mixed series: harmonic functions for the circumferential variable and three different formulations for the longitudinal variable; these three different formulations are: (a) Chebyshev orthogonal polynomials, (b) power polynomials, and (c) trigonometric functions. The same formulation is applied to study different boundary conditions; results are presented for simply supported and clamped shells. The analysis is performed in two steps: first a liner analysis is performed to identify natural modes, which are then used in the nonlinear analysis as generalized coordinates. The Lagrangian approach is applied to obtain a system of nonlinear ordinary differential equations. Different expansions involving from 14 to 34 generalized coordinates, associated with natural modes of both simply supported and clamped-clamped shells, are used to study the convergence of the solution. The nonlinear equations of motion are studied by using arclength continuation method and bifurcation analysis. Numerical responses obtained in the spectral neighborhood of the lowest natural frequency are compared with results available in literature.     

ON THE ROLES OF COMPLEMENTARY AND ADMISSIBLE BOUNDARY CONSTRAINTS IN FOURIER SOLUTIONS TO THE BOUNDARY VALUE PROBLEMS OF COMPLETELY COUPLED r TH ORDER PDES

A heretofore unavailable double Fourier series based approach, for obtaining non-separable solution to a system of completely coupled linear r th order partial differential equations with constant coefficients and subjected to general (completely coupled) boundary conditions, has been presented. The method has been successfully implemented to solve a class of hitherto unsolved boundary-value problems, pertaining to free and forced vibrations of arbitrarily laminated anisotropic doubly curved thin panels of rectangular planform, with arbitrarily prescribed (both symmetric and asymmetric with respect to the panel centerlines) admissible boundary conditions and subjected to general transverse loading.Existing solutions such as those due to Navier or Levy are based on the well-known method of separation of variables. Such solutions represent particular solutions whenever the method of separation of variables work, and when these particular solution functions fortuitously satisfy the boundary conditions. For derivation of the complementary solution, the complementary boundary constraints are introduced through boundary discontinuities of some of the particular solution functions and their partial derivatives. Such discontinuities form sets of measure zero.Various cases of lamination, geometry and dynamic response (forced and free vibrations) of a class of thin anisotropic laminated shells (curved panels) have been shown to follow from the above. Six sets of boundary conditions are used to illustrate the present method for the derivation of complementary solutions. Navier-type solutions whenever available form special cases of the present general solution.     

Free vibration analysis of circular cylindrical shells

A general analytical method is presented for evaluating the free vibration characteristics of a circular cylindrical shell with classical boundary conditions of any type. The solution is obtained through a direct solution procedure in which Sanders' shell equations are used with the axial modal displacements represented as simple Fourier series expressions. Stokes' transformation is exploited to obtain correct series expressions for the derivatives of the Fourier series. An explicit expression of the exact frequency equation can be obtained for any kind of boundary conditions. The accuracy of the method is checked against available data. The method is used to find the modal characteristics of the thermal liner model of the U.S. Fast Test Reactor (FTP). The numerical results obtained are compared with finite element method solutions.     

An Efficient Semi-Analytical Method to Compute Displacements and Stresses in an Elastic Half-Space with a Hemispherical Pit

This paper presents an efficient method to calculate the displacement and stress fields in an isotropic elastic half-space having a hemispherical pit and being subject to gravity. The method is semi-analytical and takes advantage of the axisymmetry of the problem. The Boussinesq potentials are used to obtain an analytical solution in series form, which satisfies the equilibrium equations of elastostatics, traction-free boundary conditions on the infinite plane surface and decaying conditions at infinity. The boundary conditions on the free surface of the pit are then imposed numerically, by minimising a quadratic functional of surface elastic energy. The minimisation yields a symmetric and positive definite linear system of equations for the coefficients of the series, whose particular block structure allows its solution in an efficient and robust way. The convergence of the series is verified and the obtained semi-analytical solution is then evaluated, providing numerical results. The method is validated by comparing the semi-analytical solution with the numerical results obtained using a commercial finite element software.     

Adaptive artificial boundary condition for the two-level Schrödinger equation with conical crossings

In this paper, we present an adaptive approach to design the artificial boundary conditions for the two-level Schrödinger equation with conical crossings on the unbounded domain. We use the windowed Fourier transform to obtain the local wave number information in the vicinity of artificial boundaries, and adopt the operator splitting method to obtain an adaptive local artificial boundary condition. Then reduce the original problem into an initial boundary value problem on the bounded computational domain, which can be solved by the finite difference method. By this numerical method, we observe the surface hopping phenomena of the two-level Schrödinger equation with conical crossings. Several numerical examples are provided to show the accuracy and convergence of the proposed method.     

Fourier spectral embedded boundary solution of the Poisson’s and Laplace equations with Dirichlet boundary conditions

A Fourier spectral embedded boundary method, for solution of the Poisson’s equation with Dirichlet boundary conditions and arbitrary forcing functions (including zero forcing function), is presented in this paper. This iterative method begins by transformation of the Dirichlet boundary conditions from the physical boundaries to some corresponding regular grid points (which are called the numerical boundaries), using a second order interpolation method. Then the transformed boundary conditions and the forcing function are extended to a square, smoothly and periodically, via multiplying them by some suitable error functions. Instead of direct solution of the resulting extended Poisson’s problem, it is suggested to define and solve an equivalent transient diffusion problem on the regular domain, until achievement of the steady solution (which is considered as the solution of the original problem). Without need of any numerical time integration method, time advancement of the solution is obtained directly, from the exact solution of the transient problem in the Fourier space. Consequently, timestep sizes can be chosen without stability limitations, which it means higher rates of convergence in comparison with the classical relaxation methods. The method is presented in details for one- and two-dimensional problems, and a new emerged phenomenon (which is called the saturation state) is illustrated both in the physical and spectral spaces. The numerical experiments have been performed on the one- and two-dimensional irregular domains to show the accuracy of the method and its superiority (from the rate of convergence viewpoint) to the other classical relaxation methods. Capability of the method, in dealing with complex geometries, and in presence of discontinuity at the boundaries, has been shown via some numerical experiments on a four-leaf shape geometry.     

Vibrations of segmented cylindrical shells by a fourier series component mode method

The vibrations of a multi-segment cylindrical shell with a common mean radius are studied. The shell is uniform for any segment but the material and geometric properties may vary from segment to segment. The solution is based on the component mode method coupled with Fourier series and Lagrange multipliers. It is shown that a single segment shell with boundary conditions of free support without tangential constraint is sufficient for an arbitrary shell with arbitrary boundary conditions. Results are presented for simply supported shells and clamped-free shells for two segments with different length and thickness.     

On the spectral accuracy of a fictitious domain method for elliptic operators in multi-dimensions

This work is a continuation of the authors efforts to develop high-order numerical methods for solving elliptic problems with complex boundaries using a fictitious domain approach. In a previous paper, a new method was proposed, based on the use of smooth forcing functions with identical shapes, mutually disjoint supports inside the fictitious domain and whose amplitudes play the role of Lagrange multipliers in relation to a discrete set of boundary constraints. For one-dimensional elliptic problems, this method shows spectral accuracy but its implementation in two dimensions seems to be limited to a fourth-order algebraic convergence rate. In this paper, a spectrally accurate formulation is presented for multi-dimensional applications. Instead of being specified locally, the forcing function is defined as a convolution of a mollifier (smooth bump function) and a Lagrange multiplier function (the amplitude of the bump). The multiplier function is then approximated by Fourier series. Using a Fourier Galerkin approximation, the spectral accuracy is demonstrated on a two-dimensional Laplacian problem and on a Stokes flow around a periodic array of cylinders. In the latter, the numerical solution achieves the same high-order accuracy as a Stokes eigenfunction expansion and is much more accurate than the solution obtained with a classical third order finite element approximation using the same number of degrees of freedom.     

Maxi-element analysis for free vibrations of orthotropic shells of revolution

The theory of free vibration for orthotropic shells of revolution with arbitrary homogeneous boundary conditions is developed. The essence of the method is to decompose the overall shell into a number of so-called cylindrical, conical, and plate “maxi-elements”. Since the eigenfunctions of each of the individual maxi-elements are analytically determined directly from the solution of the governing differential equations, the procedure has the advantage of requiring significantly fewer elements compared with the usual finite element or finite difference procedures. For the conical shell and the plate, the method of solution is novel, while the solution for the cylindrical shell has been published by the first author. The versatility and accuracy of the method is shown through the inclusion of a number of examples which present the excellent correlation with test results and other numerical schemes.