二元复合材料板的自由振动: 半解析法

二元复合材料板是超材料板结构中常见的单元之一. 针对由材料参数相差两个量级的基体和嵌入体组成的二元复合材料板, 提出结构自由振动的半解析模型, 并对其振动特性进行了研究. 基于区域分解法和二元材料的分布, 将二维平板分解成两个子区域. 通过在振型函数中附加区域试函数, 来描述复合材料板面内刚度突变引起局部位移和转角的非光滑性. 基于二元复合材料板的基本边界条件和两子区连接处的变形协调条件, 构造了新的振型函数. 基于经典薄板理论, 利用带特殊试函数的里兹法, 求得不同几何构型下二元复合材料板的固有频率和振型, 并研究了嵌入体的尺寸和位置对结构振动特性的影响规律. 通过收敛分析并与有限元仿真结果对比, 验证了本文方法的准确性. 研究结果表明: 传统的全局试函数在分析具有振动局部化的模态时会得到不准确的结果, 而附加区域试函数可以显著提高里兹法的收敛速度以及结果的准确性; 嵌入体位置对低阶固有频率的作用不明显, 却能显著改变低阶振型节线的分布和振动局部化发生的区域.      

An efficient method for vibration and stability analysis of rectangular plates axially moving in fluid

An efficient method is developed to investigate the vibration and stability of moving plates immersed in fluid by applying the Kirchhoff plate theory and finite element method.The fluid is considered as an ideal fluid and is described with Bernoulli’s equation and the linear potential flow theory.Hamilton’s principle is used to acquire the dynamic equations of the immersed moving plate.The mass matrix,stiffness matrix,and gyroscopic inertia matrix are determined by the exact analytical integration.The numerical results show that the fundamental natural frequency of the submersed moving plates gradually decreases to zero with an increase in the axial speed,and consequently,the coupling phenomenon occurs between the first-and second-order modes.It is also found that the natural frequency of the submersed moving plates reduces with an increase in the fluid density or the immersion level.Moreover,the natural frequency will drop obviously if the plate is located near the rigid wall.In addition,the developed method has been verified in comparison with available results for special cases.     

Vibrations and stability of a periodically supported rectangular plate immersed in axial flow

Vibrations and stability of a thin rectangular plate, infinitely long and wide, periodically supported in both directions (so that it is composed by an infinite number of supported rectangular plates with slope continuity at the edges) and immersed in axial liquid flow on its upper side is studied theoretically. The flow is bounded by a rigid wall and the model is based on potential flow theory. The Galerkin method is applied to determine the expression of the flow perturbation potential. Then the Rayleigh–Ritz method is used to discretize the system. The stability of the coupled system is analyzed by solving the eigenvalue problem as a function of the flow velocity; divergence instability is detected. The convergence analysis is presented to determine the accuracy of the computed eigenfrequencies and stability limits. Finally, the effects of the plate aspect ratio and of the channel height ratio on the critical velocity giving divergence instability and vibration frequencies are investigated.     

非局部杆纵向振动的特性研究

基于非局部理论，研究弹性杆在任意边界约束条件下的纵向振动特性．根据Chebyshev 谱级数建立非局部弹性杆的纵向位移形式．在杆的两端引入纵向约束弹簧，通过设置弹簧刚度系数，模拟经典边界及弹性边界．建立非局部杆的能量表达式，由瑞利-里兹法得到齐次线性方程组，求解对应的矩阵特征值与特征向量问题获得非局部杆的固有频率和振型．通过数值仿真计算，研究非局部特征系数与边界约束条件对非局部杆振动频率的影响．结果表明本文方法合理简便，具有良好的精度，且适用于任意弹性边界条件．     

On the immersed boundary‐lattice Boltzmann simulations of incompressible flows with freely moving objects

For simulating freely moving problems, conventional immersed boundary‐lattice Boltzmann methods encounter two major difficulties of an extremely large flow domain and the incompressible limit. To remove these two difficulties, this work proposes an immersed boundary‐lattice Boltzmann flux solver (IB‐LBFS) in the arbitrary Lagragian–Eulerian (ALE) coordinates and establishes a dynamic similarity theory. In the ALE‐based IB‐LBFS, the flow filed is obtained by using the LBFS on a moving Cartesian mesh, and the no‐slip boundary condition is implemented by using the boundary condition‐enforced immersed boundary method. The velocity of the Cartesian mesh is set the same as the translational velocity of the freely moving object so that there is no relative motion between the plate center and the mesh. This enables the ALE‐based IB‐LBFS to study flows with a freely moving object in a large open flow domain. By normalizing the governing equations for the flow domain and the motion of rigid body, six non‐dimensional parameters are derived and maintained to be the same in both physical systems and the lattice Boltzmann framework. This similarity algorithm enables the lattice Boltzmann equation‐based solver to study a general freely moving problem within the incompressible limit. The proposed solver and dynamic similarity theory have been successfully validated by simulating the flow around an in‐line oscillating cylinder, single particle sedimentation, and flows with a freely falling plate. The obtained results agree well with both numerical and experimental data. Copyright © 2016 John Wiley & Sons, Ltd.     

轴向运动导电板磁弹性非线性动力学及分岔特性

研究磁场环境中轴向运动导电薄板磁弹性动力学及分岔特性。考虑几何非线性因素,在给出薄板运动的动能、应变能及外力虚功的基础上,应用哈密顿变分原理,得到磁场中轴向运动薄板的非线性磁弹性振动方程,并给出洛伦兹电磁力的确定形式。针对横向磁场环境中条形板共振特性进行分析,应用多尺度法和奇异性理论,得到稳态运动下的分岔响应方程以及普适开折对应的转迁集。通过算例,分别得到以磁感应强度、轴向运动速度和激励力为分岔控制参数的分岔图、最大李雅普诺夫指数图和庞加莱映射图等计算结果,讨论不同分岔参数对系统呈现的倍周期和混沌运动的影响。结果表明,通过相应参数的改变可实现对系统复杂动力学行为的控制。     

三种典型轴向运动结构的振动特性对比

轴向运动结构的横向振动一直是动力学领域的研究热点之一.目前大多数的文献只涉及对一种模型的研究,而针对几种模型的对比分析较少.本文对3种典型轴向运动结构(Euler梁、窄板和对边简支对边自由的板)的振动特性进行了对比分析.针对工程中不同的结构参数,本文为其理论研究中选择更加合理的模型提供了参考.通过复模态方法求解了3种模型的控制方程,给出了其相应的固有频率及模态函数.对于板模型,同时考虑了其自由边界的两种刚体位移以及弯扭耦合振动3种情况.通过数值算例给出了3种模型的前四阶固有频率随轴速和长宽比的变化情况,并应用微分求积法对复模态方法得到的解析解进行验证.特别采用三维图的形式分析了不同的轴速、阻尼、刚度和长宽比等参数混合时对3种模型第一阶固有频率的影响,着重研究了窄板和梁的不同的长宽比和轴速混合时对两者的第一阶固有频率的相对误差的影响.结果表明:随着轴速的增大,3种模型的固有频率逐渐减小. 窄板是板的一种简化模型.在各参数值发生变化时,阻尼对第一阶固有频率的影响最小.长宽比很大,轴速很小或为零时,复杂模型可以简化为简单模型.      

Principal parametric resonance of axially accelerating rectangular thin plate in magnetic field

Nonlinear parametric vibration and stability is investigated for an axially accelerating rectangular thin plate subjected to parametric excitations resulting from the axial time-varying tension and axial time-varying speed in the magnetic field. Consid- ering geometric nonlinearity, based on the expressions of total kinetic energy, potential energy, and electromagnetic force, the nonlinear magneto-elastic vibration equations of axially moving rectangular thin plate are derived by using the Hamilton principle. Based on displacement mode hypothesis, by using the Galerkin method, the nonlinear para- metric oscillation equation of the axially moving rectangular thin plate with four simply supported edges in the transverse magnetic field is obtained. The nonlinear principal parametric resonance amplitude-frequency equation is further derived by means of the multiple-scale method. The stability of the steady-state solution is also discussed, and the critical condition of stability is determined. As numerical examples for an axially moving rectangular thin plate, the influences of the detuning parameter, axial speed, axial tension, and magnetic induction intensity on the principal parametric resonance behavior are investigated.     

AN EFFICIENT DYNAMIC MODELING METHOD OF AN AXIALLY MOVING CANTILEVER BEAM AND FREQUENCY RESPONSE ANALYSIS

轴向移动梁动力学问题具有广泛的工程应用背景,如：机械手、机床主轴、武器身管等.计算轴向移动梁动力学响应是评估结构动力学性能以及最终指导结构设计的一个重要手段.采用Rayleigh-Ritz法、拉格朗日方程推导了轴向移动悬臂梁时变动力学方程.选取幂级数函数构造试函数对轴向移动系统动力问题进行求解.幂级数函数良好的积分与微分性能,使得推导容易以矩阵的形式快速进行,便于符号运算软件直接生成MATLAB程序.由于MATLAB基本数据单位为矩阵,符号软件生成的程序只需经过简单修改便可进行动力学计算.大大缩短了轴向移动梁从建模到动力学分析的时间,过程十分高效.通过四组算例,将本文方法计算得到的动力学响应与文献数据进行对比,对该方法准确性进行了验证,并给出了幂级数函数拟合阶数的选取原则.以此为基础,研究了轴向移动梁的频率响应特性.分为考虑重力与忽略重力两种情况,讨论了轴向振动幅度对其频率响应特性的影响.     

STUDY ON THE DYNAMIC BEHAVIOR OF AXIALLY MOVING RECTANGULAR PLATES PARTIALLY SUBMERSED IN FLUID

The natural frequencies, complex modes and critical speeds of an axially moving rectangular plate, which is partially immersed in a fluid and subjected to a pretension, are investigated. The effects of free surface waves, compressibility and viscidity of the fluid are neglected in the analysis. The subsection functions are used to describe the discontinuous characteristics of the system due to partial immersion. The classical thin plate theory is adopted to formulate the equations of motion of a vibrating plate. The velocity potential and Bernoulli's equation are used to describe the fluid pressure acting on the moving plate. The effect of fluid on the vibrations of the plate may be equivalent to the added mass on the plate. The effects of distance ratio, moving speed, immersed-depth ratio, boundary conditions, stiffness ratio and aspect ratio of the plate as well as the fluid-plate density ratios on the free vibrations of the moving plate-fluid system are investigated.     

HYDROELASTIC VIBRATION OF A CIRCULAR CONTAINER BOTTOM PLATE USING THE GALERKIN METHOD

The free vibration of a flexible thin plate placed into a circular hole and elastically connected to the rigid bottom slab of a circular cylindrical container filled with fluid having a free surface is studied. The liquid is assumed to be incompressible, inviscid and irrotational. The effect of the free surface wave is also taken into account in the analysis. First of all, the exact expression of velocity potential of the liquid movement is derived by a combination of the superposition method and the method of separation of variables. With the help of the Fourier–Bessel series expansion, part of the unknown coefficients in the solution is determined by the consistency condition between the liquid movement and the plate vibration, in the form of integrals associated with the dynamic deflection of the plate. Then, the Galerkin method is applied to derive the eigenfrequency equation of the fluid–plate interaction. Finally, the effects of various parameters and the free surface wave on eigenfrequencies of the fluid–plate system are discussed. As a consequence, the accuracy of the nondimensional added virtual mass incremental (NAVMI) factor solution has also been evaluated by comparing with the more accurate Galerkin solution. It is shown that the proposed method is also applicable to the vibration analysis of circular plates in contact with an infinite liquid by only taking a finite but larger size of liquid to replace the infinite liquid in the computation.     

Dynamical analysis of axially moving plate by finite difference method

The complex natural frequencies for linear free vibrations and bifurcation and chaos for forced nonlinear vibration of axially moving viscoelastic plate are investigated in this paper. The governing partial differential equation of out-of-plane motion of the plate is derived by Newton’s second law. The finite difference method in spatial field is applied to the differential equation to study the instability due to flutter and divergence. The finite difference method in both spatial and temporal field is used in the analysis of a nonlinear partial differential equation to detect bifurcations and chaos of a nonlinear forced vibration of the system. Numerical results show that, with the increasing axially moving speed, the increasing excitation amplitude, and the decreasing viscosity coefficient, the equilibrium loses its stability and bifurcates into periodic motion, and then the periodic motion becomes chaotic motion by period-doubling bifurcation.     

Moving least squares reconstruction for sharp interface immersed boundary methods

We propose a new approach for reconstructing velocity boundary conditions in sharp-inerface immersed boundary (IB) methods based on the moving least squares (MLS) interpolation method. The MLS is employed to not only reconstruct velocity boundary conditions but also to calculate the pressure and velocity gradients in the vicinity of the immersed body, which are required in fluid structure interaction problems to obtain the force exerted by the fluid on the structure. To extend the method to arbitrarily complex geometries with nonconvex shaped boundaries, the visibility method is combined with the MLS method. The performance of the proposed curvilinear IB MLS (CURVIB-MLS) is demonstrated by systematic grid-refinement studies for two- and three-dimensional tests and compared with the standard CURVIB method employing standard wall-normal interpolation for reconstructing boundary conditions. The test problems are flow in a lid-driven cavity with a sphere, uniform flow over a sphere, flow on a NACA0018 airfoil at incidence, and vortex-induced vibration of an elastically-mounted cylinder. We show that the CURVIB-MLS formulation yields a method that is easier to implement in complex geometries and exhibits higher accuracy and rate of convergence relative to the standard CURVIB method. The MLS approach is also shown to dramatically improve the accuracy of calculating the pressure and viscous forces imparted by the flow on the body and improve the overall accuracy of FSI simulations. Finally, the CURVIB-MLS approach is able to qualitatively capture on relatively coarse grids important features of complex separated flows that the standard CURVIB method is able to capture only on finer grids.     

采用压电材料提高超声速飞行器壁板结构的颤振特性

采用压电材料对结构进行振动主动控制已经进行了广泛研究,论文进一步采用压电材料改进超声速壁板结构的气动弹性颤振特性,研究中考虑压电材料力电耦合效应的影响.采用Hamilton原理和Rayleigh-Ritz方法建立壁板及压电材料整体结构的运动方程,采用超声速活塞理论模拟气动力,利用加速度反馈控制策略对压电材料施加外电压,获得结构的主动质量.求解运动方程的特征值问题获得固有频率,进而确定气动弹性颤振边界,分析了反馈控制增益对超声速飞行器壁板结构主动颤振特性的影响,研究表明,采用压电材料可以提高超声速壁板结构的气动弹性颤振特性.     

An immersed boundary method for incompressible flows using volume of body function

A simple and effective immersed boundary method using volume of body (VOB) function is implemented on unstructured Cartesian meshes. The flow solver is a second‐order accurate implicit pressure‐correction method for the incompressible Navier–Stokes equations. The domain inside the immersed body is viewed as being occupied by the same fluid as outside with a prescribed divergence‐free velocity field. Under this view a fluid–body interface is similar to a fluid–fluid interface encountered in the volume of fluid (VOF) method for the two‐fluid flow problems. The body can thus be identified by the VOB function similar to the VOF function. In fluid–body interface cells the velocity is obtained by a volume‐averaged mixture of body and fluid velocities. The pressure inside the immersed body satisfies the same pressure Poisson equation as outside. To enhance stability and convergence, multigrid methods are developed to solve the difference equations for both pressure and velocity. Various steady and unsteady flows with stationary and moving bodies are computed to validate and to demonstrate the capability of the current method. Copyright © 2005 John Wiley & Sons, Ltd.     

Generation of Transient Waves by Impulsive Disturbances in an Inviscid Fluid with an Ice-Cover

The dynamic responses of an ice-covered fluid to impulsive disturbances are analytically investigated for two- and three-dimensional cases. The initially quiescent fluid of infinite depth is assumed to be inviscid, incompressible and homogenous. The thin ice-cover is modelled as a homogenous elastic plate with negligible inertia. Four types of impulsive concentrated disturbances are considered, namely an instantaneous mass source immersed in the fluid, an instantaneously dynamic load on the plate, an initial impulse on the surface of the fluid, and an initial displacement of the ice plate. The linearized initial-boundary-value problem is formulated within the framework of potential flow. The solutions in integral form for the vertical deflexions at the ice-water interface are obtained by means of a joint Laplace-Fourier transform. The asymptotic representations of the wave motions for large time with a fixed distance-to-time ratio are derived by making use of the method of stationary phase. It is found that there exists a minimal group velocity and the wave system observed depends on the moving speed of the observer. For an observer moving with the speed larger than the minimal group velocity, there exist two trains of waves, namely the long gravity waves and the short flexural waves, the latter riding on the former. Moreover, the deflexions of the ice-plate for an observer moving with a speed near the minimal group velocity are expressed in terms of the Airy functions. The effects of the presence of an ice-cover on the resultant wave amplitudes, the wavelengths and periods are discussed in detail. The explicit expressions for the free-surface gravity waves can readily be recovered by the present results as the thickness of ice-plate tends to zero.     

Free vibration and stability of an axially moving thin circular cylindrical shell using multiple scales method

This paper investigates the linear free vibration of axially moving simply supported thin circular cylindrical shells with constant and time-dependent velocity considering the effect of viscous structure damping. Classical shell theory is employed to express strain-displacement relation. Linear elasticity theory is used to write stress–strain relation considering Hook’s Law. Governing equations in cylindrical coordinates are derived using the Hamilton principle. Equilibrium equations are rewritten with the help of Donnell–Mushtari shell theory simplification assumptions. Motion equations for displacements in axial and circumferential directions are solved analytically concerning to displacement in the radial direction. As the displacement in the radial direction is the combination of driven and companion modes, the third motion equation is discretized using the Galerkin method. The set of ordinary differential equation obtained from the Galerkin method is solved using the steady-state method, which in practice leads to the prediction of the exact frequencies of vibration. By employing multiple scale method the critical speed values of a circular cylindrical shell and several types of instabilities are discussed. The numerical results show that by increasing the mean velocity, the system always loses stability by the divergence instability in different modes, and the critical speed values of lower modes are higher than those of higher modes. As well as the unstable regions for the resonances between velocity function fluctuation frequencies and the linear combination of natural frequencies is gained from the solvability condition of second order multiple scale method. The accuracy of the method is checked against the available data.

     

A stream function–vorticity formulation‐based immersed boundary method and its applications

A new stream function–vorticity formulation‐based immersed boundary method is presented in this paper. Different from the conventional immersed boundary method, the main feature of the present model is to accurately satisfy both governing equations and boundary conditions through velocity correction and vorticity correction procedures. The velocity correction process is performed implicitly based on the requirement that velocity at the immersed boundary interpolated from the corrected velocity field accurately satisfies the nonslip boundary condition. The vorticity correction is made through the stream function formulation rather than the vorticity transport equation. It is evaluated from the firstorder derivatives of velocity correction. Two simple and efficient ways are presented for approximation of velocity‐correction derivatives. One is based on finite difference approximation, while the other is based on derivative expressions of Dirac delta function and velocity correction. It was found that both ways can work very well. The main advantage of the proposed method lies in its simple concept, easy implementation, and robustness in stability. Numerical experiments for both stationary and moving boundary problems were conducted to validate the capability and efficiency of the present method. Good agreements with available data in the literature were achieved. Copyright © 2011 John Wiley & Sons, Ltd.     

An investigation into non-linear vibrations of thin plates

The variational and modified forms of the von Kármán-type non-linear plate equations are considered in the context of the Rayleigh-Ritz and Galerkin methods. An approximate analysis of the non-linear vibrations of thin elastic plates including inplane inertia is presented. The quantitative study confirms that the inplane inertia effects are negligible for thin plates provided the non-linearity is not too large. It is observed that the non-linear inertia terms in the transverse equation of motion should be retained in any such study. The analysis is simplified by neglecting the inplane inertia and applied to constrained and unconstrained plates. A different type of inplane boundary condition termed ‘the partially constrained’ is studied, and the inadequacy of replacing the unconstrained condition by means of an average-zero stress condition is clearly demonstrated. It is observed that in most of the cases considered the Galerkin method yields lower bounds for the non-linear coefficient of the modal equation. In all cases the Galerkin results yield less stiff models than the Rayleigh-Ritz method. The general significance of the convergence of the two methods beyond the scope of the title problem is highlighted.     

功能梯度旋转圆板的热弹性固有振动

在考虑热因素及旋转运动条件下，针对金属－陶瓷功能梯度圆板的固有振动问题进行研究．给出随温度变化且材料组分沿厚度方向按幂律分布的材料物性参数，依据热弹性理论得到圆板的能量关系式．基于哈密顿原理建立旋转金属－陶瓷功能梯度圆板热弹性动力学方程．采用伽辽金法得到边界约束下圆板的自由振动方程，确定了静挠度及固有振动频率．基于数值计算，得到系统固有频率值随体积分数指数、转速和温度等参量的变化曲线，讨论了静挠度变化规律及动力系统的奇点稳定性问题．结果表明，固有频率随体积分数指数、材料表面温度以及转速的增加而减小．