Implementation of a constrained Ritz series modeling technique for acoustic cavity-structural systems

A prior study [Ginsberg, J. H. (2010b). J. Acoust. Soc. Am. 127, 2749-2758] used Ritz series in conjunction with Hamilton's principle to derive general equations describing the time domain response of an acoustic cavity bounded by an elastic structure. The equations of motion are supplemented by constraint equations that explicitly enforce velocity continuity at the cavity's surface. These constraints are imposed by the surface traction, which is represented by unknown coefficients of Ritz-type series. The resulting set of equations are differential-algebraic type. Three methods are presented to convert the governing equations to forms that are familiar to structural acoustics, including one that transforms them from differential-algebraic type to the standard ordinary differential equations associated with linear multi-degree-of-freedom vibratory systems. In cases where only the structure is excited, the formulation offers options as to how displacement/velocity boundary conditions on the nonstructural boundary are enforced, as well as whether zero pressure boundary conditions are enforced at all. An example of a one-dimensional waveguide that is closed at one end by an oscillator is used to explore the quality of solutions obtained from each of these options. Results for natural frequencies and mode functions are examined for accuracy and convergence.     

球状物体作微形变时弹性势能的研究

研究一个弹性球在保持体积不变的前提下,形变为长旋转椭球时的弹性势能.首先通过最小势能原理和物理条件确定位移函数,然后根据弹性势能密度-应变-位移三者之间的关系,通过积分计算得到弹性势能,并且进行了有益的讨论.     

An alternative approach for the analysis of nonlinear vibrations of pipes conveying fluid

Based on a novel extended version of the Lagrange equations for systems containing non-material volumes, the nonlinear equations of motion for cantilever pipe systems conveying fluid are deduced. An alternative to existing methods utilizing Newtonian balance equations or Hamilton's principle is thus provided. The application of the extended Lagrange equations in combination with a Ritz method directly results in a set of nonlinear ordinary differential equations of motion, as opposed to the methods of derivation previously published, which result in partial differential equations. The pipe is modeled as a Euler elastica, where large deflections are considered without order-of-magnitude assumptions. For the equations of motion, a dimensional reduction with arbitrary order of approximation is introduced afterwards and compared with existing lower-order formulations from the literature. The effects of nonlinearities in the equations of motion are studied numerically. The numerical solutions of the extended Lagrange equations of the cantilever pipe system are compared with a second approach based on discrete masses and modeled in the framework of the multibody software HOTINT/MBS. Instability phenomena for an increasing number of discrete masses are presented and convergence towards the solution for pipes conveying fluid is shown.     

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.     

A modified Fourier series solution for vibration analysis of truncated conical shells with general boundary conditions

Free vibration analysis of truncated conical shells with general elastic boundary conditions is presented in this paper. An accurate modified Fourier series solution is developed, in which, regardless of the boundary conditions, each displacement of the conical shell is invariantly expressed as a new form of improved series expansions composed of a standard Fourier series and closed-form auxiliary functions introduced to ensure and accelerate the convergence of the series expansion. All the expansion coefficients are treated as the generalized coordinates and determined using the Rayleigh–Ritz method. By using the present method, conical shells with arbitrary boundary conditions including all classical and elastic end restraints can be solved in a unified form. The accuracy and convergence of the current approach are validated by numerical examples and comparison with FEM results and those from the literature, and excellent accuracy is demonstrated. Comprehensive studies on the effects of elastic restraint parameters, semi-vertex angle and the ratio of length to radius are also reported. Some new results are presented for cases with elastic boundary restraints which may serve as benchmark solution for future researches.     

Stability of an axially accelerating string subjected to frictional guiding forces

The dynamic response of an axially translating continuum subjected to the combined effects of a pair of spring supported frictional guides and axial acceleration is investigated; such systems are both non-conservative and gyroscopic. The continuum is modeled as a tensioned string translating between two rigid supports with a time-dependent velocity profile. The equations of motion are derived with the extended Hamilton's principle and discretized in the space domain with the finite element method. The stability of the system is analyzed with the Floquet theory for cases where the transport velocity is a periodic function of time. Direct time integration using an adaptive step Runge-Kutta algorithm is used to verify the results of the Floquet theory. This approach can also be employed in the general case of arbitrary time-varying velocity. Results are given in the form of time history diagrams and instability point grids for different sets of parameters such as the location of the stationary load, the stiffness of the elastic support, and the values of initial tension. This work showed that presence of friction adversely affects stability, but using non-zero spring stiffness on the guiding force has a stabilizing effect. This work also showed that the use of the finite element method and Floquet theory is an effective combination to analyze stability in gyroscopic systems with stationary friction loads.     

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.     

Vibration of electro-elastic versus magneto-elastic circular/annular plates using the Chebyshev–Ritz method

The aim of this paper is to analyze three-dimensional free vibration of magneto-elastic/electro-elastic circular/annular plates with different boundary conditions using the Chebyshev–Ritz method, in which a set of duplicate Chebyshev polynomial series multiplied by the boundary function satisfying the boundary conditions are chosen as the trial functions of the displacement components, the electric potential and the magnetic potential. Convergence of the method is checked using various Chebyshev polynomial terms. The effect of geometrical parameters and material properties of magneto-elastic/electro-elastic circular/annular plates on the eigenfrequencies of free vibration is considered.     

Vibration of a two-member open frame

The dynamics of a two-member open frame structure undergoing both in- and out-of-plane motion is examined. The frames are modelled using the Euler-Bernoulli beam theory and are further generalized by permitting an arbitrary angle between the beams and the attachment of a payload at the end of the second beam. The equations of motion are derived using Hamilton's principle and the orthogonality conditions are presented. It is shown that the in- and out-of-plane motions can be decoupled by including the axial deformation components in the assumed displacement fields. The natural frequencies of the system and the contribution of each member into the system potential energy are examined via numerical examples.     

封闭声腔的结构-声耦合的瞬态响应及固有特性分析的时域有限体积法

基于波动理论,采用时域有限体积法(TDFVM)研究封闭声腔结构-声耦合问题的瞬态响应及固有特性。该方法在结构与流体区域分别求解结构动力学方程与非均匀介质中的声波动方程,根据流固分界面上的力平衡与质点振速连续条件实现结构-声耦合。通过数值算例验证方法的正确性和精确性,在此基础上研究封闭声腔结构与空气耦合的瞬态响应及固有特性,分析水深变化对耦合系统声振特性和固有特性的影响。结果表明,随着水深的增加,结构与水的耦合会更加强烈,导致耦合系统的特征频率降低,同时空气腔深度的减小,导致声腔的部分固有频率增加。该法对计算机内存要求低,且可以考虑含有非均匀流体的结构-声耦合。      

Natural vibrations and stability of shells of revolution interacting with an internal fluid flow

A finite-element algorithm is proposed to investigate the dynamic behavior of elastic shells of revolution containing a quiescent or a flowing inviscid fluid in the framework of linear theory. The fluid behavior is described using the perturbed velocity potential. The shell behavior is treated in the framework of the classical shell theory and variational principle of virtual displacements incorporating a linearized Bernoulli equation for calculation of hydrodynamic pressure acting on the shell. The problem reduces to evaluation and analysis of the eigenvalues in the connected system of equations obtained by coupling the equations for velocity perturbations with the equations for shell displacements. For cylindrical shells, the results of numerical simulations are compared with recently published experimental, analytical and numerical data. The paper also reports the results of studying the dynamic behavior of shells under various boundary conditions for the perturbed velocity potential. The investigation made for conical shells has shown that under certain conditions an increase in the cone angle can change a divergent type of instability to a flutter type.     

NON-LINEAR FORCED VIBRATIONS OF PLATES BY AN ASYMPTOTIC-NUMERICAL METHOD

Non-linear forced vibrations of thin elastic plates have been investigated by an asymptotic-numerical method (ANM). Various types of harmonic excitation forces such as distributed and concentrated are considered. Using the harmonic balance method and Hamilton's principle, the equation of motion is converted into an operational formulation. Based on the finite element method a starting point corresponding to a non-linear solution associated to a given frequency and amplitude of excitation is computed. Applying perturbation techniques in the vicinity of this solution, the non-linear governing equation obtained is transformed into a sequence of linear problems having the same stiffness matrix. Employing one matrix inversion, a large number of terms of the perturbation series of the displacement and frequency can be easily computed with a small computation time. Iterations of this method lead to a powerful path-following technique. Comprehensive numerical tests for forced vibrations of plates subjected to time-harmonic lateral excitations are reported.     

Numerical solution for curved crack problem in elastic half-plane using hypersingular integral equation

A hypersingular integral equation for the curved crack problems of an elastic half-plane is introduced. Formulation of the equation is based on the usage of a modified complex potential. The potential is generally expressed in the form of a Cauchy-type integral. The modified complex potential is composed of the principal part and the complementary part. The principal part of the complex potential is actually equivalent to the original complex potential for the curved crack in an infinite plate. The role of the complementary part is to eliminate the boundary traction along the boundary of the half-plane caused by the principal part. From the assumed boundary traction condition, a hypersingular integral equation is obtained for the curved crack problems of an elastic half-plane. The curve length coordinate method is used to obtain a final solution. Several numerical examples are presented that prove the efficiency of the suggested method.     

Analysis of transient responses in a laminated piezoelectric cylindrical shell

A hybrid numerical method is proposed for analysis of transient responses in a multilayered piezoelectric cylindrical shell.In the present method,the associated equations of the displacement field and the electro-potential field are developed using an analytical-numerical method.The piezoelectric cylindrical shell is discretized into layered annular elements along the wall thickness direction.The governing equations are determined by Hamilton's Principle considering the coupling between the elastic and elec...     

Free-vibration acoustic resonance of a nonlinear elastic bar

Free-vibration acoustic resonance of a one-dimensional nonlinear elastic bar was investigated by direct analysis in the calculus of variations. The Lagrangian density of the bar includes a cubic term of the deformation gradient, which is responsible for both geometric and constitutive nonlinearities. By expanding the deformation function into a complex Fourier series, we derived the action integral in an analytic form and evaluated its stationary conditions numerically with the Ritz method for the first three resonant vibration modes. This revealed that the bar shows the following prominent nonlinear features: (i) amplitude dependence of the resonance frequency; (ii) symmetry breaking in the vibration pattern; and (iii) excitation of the high-frequency mode around nodal-like points. Stability of the resonant vibrations was also addressed in terms of a convex condition on the strain energy density.     

A variational formulation for the finite element analysis of the vibratory and acoustical response of high speed machinery

The principle of virtual velocities is employed herein to develop a variational theorem for determining the vibratory and acoustical response of high speed mechanisms immersed in a perfect fluid (air). Both the solid and fluid media are modeled as continua and the basic functional expression is generalized by using Lagrange multipliers to incorporate field equations, boundary conditions and a kinematic constraint on the intetface region between the two types of continua. The resulting mixed variational equation of motion then provides the basis for finite element analyses of these acousto-mechanical systems by mixed and displacement (velocity) formulations.     

各向异性立方晶体的弹性格林函数及其应用

本文计算了各向异性立方晶体的弹性格林函数的级数展开式,给出了直到二级近似的展开式的系数。将所得结果应用于弹性偶极子模型,给出了对称中心所产生的位移场及两对称中心间的弹性相互作用的表示式。应用于Cu,K等强各向异性立方晶体,虽然级数的收敛较慢,但所得关于对称中心的位移场,及二对称中心间的互作用能的数值结果,竟与基于点阵的不连续性作出的点阵静力学计算所得的结果基本一致。从而表明,本文给出的直到二级近似的弹性格林函数的解析表示提供了一个可以普遍应用的简便的方法。它可以较准确地描述立方晶体的某些力学行为。 关键词：     

Distribution functions for the Brownian motion of particles in a periodic potential driven by an external force

The stationary distribution functions for the Brownian motion of particles driven by an external force are calculated by expanding the velocity part into Hermite functions and the space part into a Fourier series. Insertion into the Fokker-Planck equation leads to a matrix continued fraction for the lowest two coefficients of the Hermite functions. Higher order terms are found by reverse iteration. Results are shown for a cosine potential. The good convergence allows the calculation in the full range of damping constants. For small friction the distribution function is in good agreement with previous results and the maxima are given by the solutions without noise.     

Vibration analysis of a partially treated multi-layer beam with magnetorheological fluid

Semi-active vibration control based on magnetorheological (MR) materials offers excellent potential for high bandwidth control through rapid variations in the rheological properties of the fluid under varying magnetic field. Such fluids may be conveniently applied to partial or more critical components of a large structure to realize more efficient and compact vibration control mechanism with variable damping. This study investigates the properties and vibration responses of a partially treated multi-layer MR fluid beam. The governing equations of a partially treated multi-layered MR beam are formulated using finite element method and Ritz formulation. The validity of the proposed finite element formulations is demonstrated by comparing the results with those obtained from the Ritz formulation and the experimental measurements. The properties of different configurations of a partially treated MR-fluid beam are evaluated to investigate the influences of the location and length of the MR-fluid for different boundary conditions. The properties in terms of natural frequencies and loss factors corresponding to various modes are evaluated under different magnetic field intensities and compared with those of the fully treated beams. The effect of location of the fluid treatment on deflection mode shapes is also investigated. The forced vibration responses of the various configurations of partially treated MR sandwich beam are also evaluated under harmonic force excitations. The results suggest that the natural frequencies and transverse displacement response of the partially treated MR beams are strongly influenced not only by the intensity of the applied magnetic field, but also by the location and the length of the fluid pocket. The application of partial treatment could also alter the deflection pattern of the beam, particularly the location of the peak deflection.     

On a more precise statement of Hamilton's principle

It has been recognized in the literature of the calculus of variations that the classical statement of the principle of least action (Hamilton's principle for conservative systems) is not strictly correct. Recently, mathematical proofs have been offered for what is claimed to be a more precise statement of Hamilton's principle for conservative systems. According to a widely publicized version of this more precise statement, the action integral for conservative systems is a minimum for discrete systems for small time intervals only and is never minimum for continuous systems. In this paper, two contradictions to this more precise statement are demonstrated, one for a discrete system and one for a continuous system.