A NEW MATRIX PERTURBATION METHOD FOR ANALYTICAL SOLUTION OF THE COMPLEX MODAL EIGENVALUE PROBLEM OF VISCOUSLY DAMPED LINEAR VIBRATION SYSTEMS

A new matrix perturbation analysis method is presented for efficient approximate solution of the complex modal quadratic generalized eigenvalue problem of viscously damped linear vibration systems. First, the damping matrix is decomposed into the sum of a proportional-and a nonproportional-damping parts, and the solutions of the real modal eigenproblem with the proportional dampings are determined, which are a set of initial approximate solutions of the complex modal eigenproblem. Second, by taking the nonproportional-damping part as a small modification to the proportional one and using the matrix perturbation analysis method, a set of approximate solutions of the complex modal eigenvalue problem can be obtained analytically. The result is quite simple. The new method is applicable to the systems with viscous dampings-which do not deviate far away from the proportional-damping case. It is particularly important that the solution technique be also effective to the systems with heavy, but not over, dampings. The solution formulas of complex modal eigenvlaues and eigenvectors are derived up to second-order perturbation terms. The effectiveness of the perturbation algorithm is illustrated by an exemplar numerical problem with heavy dampings. In addition, the practicability of approximately estimating the complex modal eigenvalues, under the proportional-damping hypothesis, of damped vibration systems is discussed by several numerical examples.

A symmetric substructuring method for analyzing the natural frequencies of conical origami structures

Conical origami structures are characterized by their substantial out-of-plane stiffness and energy-absorptioncapacity.Previous investigations have commonly focused on the static characteristics of these lightweight struc-tures.However,the efficient analysis of the natural vibrations of these structures is pivotal for designing conicalorigami structures with programmable stiffness and mass.In this paper,we propose a novel method to analyzethe natural vibrations of such structures by combining a symmetric substructuring method(SSM)and a gener-alized eigenvalue analysis.SSM exploits the inherent symmetry of the structure to decompose it into a finiteset of repetitive substructures.In doing so,we reduce the dimensions of matrices and improve computationalefficiency by adopting the stiffness and mass matrices of the substructures in the generalized eigenvalue analysis.Finite element simulations of pin-jointed models are used to validate the computational results of the proposedapproach.Moreover,the parametric analysis of the structures demonstrates the influences of the number of seg-ments along the circumference and the radius of the cone on the structural mass and natural frequencies of thestructures.Furthermore,we present a comparison between six-fold and four-fold conical origami structures anddiscuss the influence of various geometric parameters on their natural frequencies.This study provides a strategyfor efficiently analyzing the natural vibration of symmetric origami structures and has the potential to contributeto the efficient design and customization of origami metastructures with programmable stiffness.

Inverse Eigenvalue Problem in Structural Dynamics Design

A kind of inverse eigenvalue problem in structural dynamics design is considered. The problem is formulated as an optimization problem. The properties of this problem are analyzed, and the existence of the optimum solution is proved. The directional derivative of the objective function is obtained and a necessary condition for a point to be a local minimum point is given. Then a numerical algorithm for solving the problem is presented and a plane-truss problem is discussed to show the applications of the theories and the algorithm.

Time Collocation Method for Structural Dynamic Problems

In order to achieve highly accurate and efficient numerical calculations of structural dynamics, time collocation method is presented. For a given time interval, the numerical solution of the method is approximated by a polynomial. The polynomial coefficients are evaluated by solving algebraic equation. Once the polynomial coefficients are evaluated, the numerical solutions at any time in the interval can be easily calculated. New formulae are derived for the polynomial coefficients,which are more practical and succinct than those previously given. Two structural dynamic equations are calculated by the proposed method. The numerical solutions are compared with the traditional fourth-order Runge-Kutta method. The results show that the method proposed is highly accurate and computationally efficient. In addition, an important advantage of the method is the simplicity in software programming.

Mixed Finite Volume Element Method for Vibration Equations of Beam with Structural Damping

In this paper, for the initial and boundary value problem of beams with structural damping, by introducing intermediate variables, the original fourth-order problem is transformed into second-order partial differential equations, and the mixed finite volume element scheme is constructed, and the existence, uniqueness and convergence of the scheme are analyzed. Numerical examples are provided to confirm the theoretical results. In the end, we test the value of δ to observe its influence on the model.

Fully coupled flow-induced vibration of structures under small deformation with GMRES method

Lagrangian-Eulerian formulations based on a generalized variational principle of fluid-solid coupling dynamics are established to describe flow-induced vibration of a structure under small deformation in an incompressible viscous fluid flow. The spatial discretization of the formulations is based on the multi-linear interpolating functions by using the finite element method for both the fluid and solid structures. The generalized trapezoidal rule is used to obtain apparently non-symmetric linear equations in an incremental form for the variables of the flow and vibration. The nonlinear convective term and time factors are contained in the non-symmetric coefficient matrix of the equations. The generalized minimum residual （GMRES） method is used to solve the incremental equations. A new stable algorithm of GMRES-Hughes-Newmark is developed to deal with the flow-induced vibration with dynamical fluid-structure interaction in complex geometries. Good agreement between the simulations and laboratory measurements of the pressure and blade vibration accelerations in a hydro turbine passage was obtained, indicating that the GiViRES-Hughes-Newmark algorithm presented in this paper is suitable for dealing with the flow-induced vibration of structures under small deformation.

Eigenvalue Analysis of Thin Plate with Complicated Shapes By a Novel Infinite Element Method

A novel infinite element method(IEM)is presented for solving plate vibration problems in this paper.In the proposed IEM,the substructure domain is partitioned into multiple layers of geometrically similar finite elements which use only the data of the boundary nodes.A convergence criterion based on the trace of the mass matrix is used to determine the number of layers in the IE model partitioning process.Furthermore,in implementing the Craig-Bampton(CB)reduction method,the inversion of the global stiffness matrix is calculated using only the stiffness matrix of the first element layer.The validity and performance of the proposed method are investigated by means of four illustrative problems.The first example considers the case of a simple clamped rectangular plate.It is observed that the IEM results are consistent with the theoretical results for first six natural frequencies.The second example considers the frequency response of a clamped rectangular plate with a crack.The main feature of IEM is that a very fine and good quality virtual mesh can be created around the crack tip.The third and fourth examples consider the natural frequency of a multiple point supported plate and a perforated plate,respectively.The results are obtained just need to adjust the reference point or boundary nodes.The parametric analyses for various geometric profiles are easy to be conducted using these numerical techniques.In general,the results presented in this study have shown that the proposed method provides a direct,convenient and accurate tool for eigenvalue analysis of thin plate structure with complicated shapes.

A NUMERICAL CALCULATION METHOD FOR EIGENVALUE PROBLEMS OF NONLINEAR INTERNAL WAVES

Generally speaking, the background shear current U（z） must be taken into account in eigenvalue problems of nonlinear internal waves in ocean, as is different from those of linear internal waves. A numerical calculation method for eigenvalue problems of nonlinear internal waves is presented in this paper on the basis of the Thompson-Haskell＇s calculation method. As an application of this method, at a station （21°N, 117°15′E） in the South China Sea, a modal structure and parameters of nonlinear internal waves are calculated, and the results closely agree with the calculated results based on observation by Yang et al..

ACG-TypeMethod for InverseQuadratic Eigenvalue Problems in Model Updating of Structural Dynamics

In this paper we first present a CG-type method for inverse eigenvalue problem of constructing real and symmetric matrices M,D and K for the quadratic pencil Q(λ)=λ^(2)M+λD+K,so that Q(λ)has a prescribed subset of eigenvalues and eigenvectors.This method can determine the solvability of the inverse eigenvalue problem automatically.We then consider the least squares model for updating a quadratic pencil Q(λ).More precisely,we update the model coefficient matrices M,C and K so that(i)the updated model reproduces the measured data,(ii)the symmetry of the original model is preserved,and(iii)the difference between the analytical triplet(M,D,K)and the updated triplet(M_(new),D_(new),K_(new))is minimized.In this paper a computationally efficient method is provided for such model updating and numerical examples are given to illustrate the effectiveness of the proposed method.

基于Dirichlet-to-Neumann映射计算太赫兹频段旋电光子晶体带隙结构

用扩展Dirichlet-to-Neumann(DtN)映射方法计算太赫兹频段二维旋电光子晶体的带隙结构,由于DtN映射方法仅需在单元晶格的边界上进行离散,避免在其内部离散,因此减少了未知数的个数,使计算速度大幅度加快.首先,构造旋电光子晶体单元晶格的DtN映射;其次,将光子晶体的带隙结构问题转化为矩阵的特征值问题进行求解;最后,数值模拟验证用DtN映射方法计算旋电光子晶体带隙结构的有效性.

含内支承的弹性边界约束梁振动功率流研究

为了有效解决弯矩、剪力等动力学参量在支承处跳跃问题,本文建立弹性边界条件下含内部支承梁结构振动功率流分析模型,将梁结构在支承位置分割为子梁I-耦合弹簧-子梁II的形式。梁结构振动位移采用边界光滑傅里叶级数进行分段展开,根据能量原理建立考虑内支承刚度的任意边界约束梁结构系统拉格朗日函数,结合瑞利-里兹方法对模型求解获得系统特征方程。在此基础上,根据振动功率流理论计算结构振动系统功率流分布与振动响应特性。通过与传统方法、数值求解对比,对本文模型的正确性与可靠性给予验证。基于所建立模型,研究横向支承刚度、旋转约束刚度、内部支承位置等参数对任意弹性边界梁结构振动功率流、支承传递损失、系统总输出功率流的影响。本文模型解决了传统方法求解突变位置功率流的数值不稳定问题,且系统边界变化时无需对模型进行重建,具有高效、高精度等优点,为弹性边界条件内部支承梁结构振动功率流特性分析提供了有效手段。

基于Dirichlet-to-Neumann映射分析计算旋电光子晶体波导结构

拓展了Dirichlet-to-Neumann(DtN)映射方法,将其应用于分析二维旋电光子晶体波导结构.首先基于旋电光子晶体单元晶格的DtN映射,构造出波导结构超级晶格的DtN映射,然后在超级晶格的边界上建立起线性特征值问题进行求解.由于该方法只需要在晶格的边界上进行离散,避免了在计算区域内部的离散,故所涉及的矩阵规模相对较小,极大地提高了计算速度.文章最后用数值算例验证了该方法的有效性.

任意凸四边形区域上二阶椭圆特征值问题基于高斯点的一种有效的谱配置法

作者提出了任意凸四边形区域上基于高斯点的二阶椭圆特征值问题的一种有效谱配置法.该方法首先利用等参变换将任意凸四边形区域上的函数转化为[-1,1]×[-1,1]上的函数,然后根据边界条件,可根据Legendre多项式的正交性构造一组有效的基函数,将数值解表示为这组基函数的展开组合.再次,通过编程计算出每个基函数在这些高斯点上的节点值,将离散格式推导为一个线性的矩阵特征系统.最后,给出了一些数值算例来表明算法的正确性和有效性,数值结果表明了该方法是有效的和收敛的.

非轴对称条件下用三维模式匹配法计算电阻率测井响应

采用模式匹配法研究了非轴对称条件下普通电阻率测井响应.根据非轴对称的地层模型,建立了新型的坐标系,确定在井轴与地层法线所形成的XOY平面上用数值方法,与此面垂直的方向上用解析方法.在XOY平面上的数值分析中,构建了广义特征值问题;在Z轴方向的解析部分中,人为地划分了一些介质层,并根据电磁场在层界面的连续性条件,推导出了各层之间场量的传递关系,从而摒弃了传统上二维模式匹配法中的复杂的反射阵和透射阵理论,完成了三维模式匹配理论.多种地层模型的数值结果对比表明,该方法在轴对称条件下的计算结果和二维模式匹配法的计算结果具有很好的一致性.此外,应用此算法还研究了倾斜地层等各种非轴对称的地层模型电阻率测井响应,系统地考察了井斜角对普通电阻率测井响应的影响.

快速子空间迭代法、迭代Ritz向量法与迭代Lanczos法的比较

以高效的细胞稀疏直接快速解法为核心步骤,实现了快速的固有振动广义特征值问题解法,并在相同的允许模态误差的意义下检验了三种常用的大型矩阵特征模态算法——子空间迭代法、迭代Ritz向量法和迭代Lanc-zos法的计算效率。迭代Ritz向量法平均最快,子空间迭代法最慢,三种解法效率相差不是太大。与ANSYS的子空间迭代和Lanczos法相比,本文的子空间迭代比ANSYS的效率高很多,Lanczos法和ANSYS的效率差不多。大量较大规模的例题显示,本文对特征值算法的改进是十分有效的,算法的健壮性,通用性都达到了高水平。

考虑索局部振动的斜拉桥动力特性研究——杭州湾跨海大桥北航道桥动力特性分析

以主跨为448m的杭州湾跨海大桥北航道桥为例,研究了考虑拉索局部振动的斜拉桥动力特性.利用子结构方法,把索的振型作为广义自由度来考虑拉索局部振动.该方法与把拉索划分为多个单元的链杆法相比,具有自由度少并且计算精度不受影响的优点.研究结果表明,当考虑拉索局部振动时,会出现拉索与主梁或桥塔的耦合振型,这种索桥耦合振动现象在以往把拉索作为单一单元的斜拉桥动力特性分析中不会出现.索的局部振动虽然对结构主要振型的频率影响不大,但对振型能量分布有比较明显的影响.对发生索桥共振的振型,如不考虑索的局部振动,通常会使索的阻尼估计过高,而索的振动能量估计偏低,不能很好地反映大跨度斜拉桥的动力特性.

地震作用下PTMD对高层钢结构的减振效果

针对传统调谐质量阻尼器(TMD)的不足,本文研究了粘弹性碰撞调谐质量阻尼器(viscoelastic pounding tuned mass damper,PTMD)对高层钢结构的减振问题,采用数值仿真对PTMD进行了参数研究和优化。碰撞材料和预留间隙对PTMD的减振效果至关重要,若碰撞材料等效弹性模量较低,减振率将随着间隙的增大而减小,最后趋于恒定;碰撞材料等效弹性模量较高,减振率随预留间隙的增大,先上升后下降,最后稳定在一个定值。PTMD在不同地震强度和不同地震作用下,对高层钢结构均有较好的减振效果。从减振效果、空间要求和附加质量块重量三方面对比了TMD与PTMD,结果表明:PTMD可以大幅降低阻尼的功率要求、安装空间限制和质量块重量,比TMD更加经济有效。

变截面梁式动力吸振器的宽带吸振机理

通过对变截面梁-板耦合结构的功率流研究,提出用变截面梁吸收弹性体的振动,论证了变截面梁式动力吸振器的宽带吸振性能,分析了梁的结构参数、梁在板上的连接位置等对吸振效果的影响,为设计适合于弹性主振系的新型梁式宽带动力吸振器奠定了基础。仿真结果表明梁式动力吸振器吸振频带宽、吸振效果好,加之其结构简单,安装方便,有一定的工程应用前景。

输液管道流固耦合振动的数值分析

采用特征线法对输液管道流固耦合振动问题进行了数值分析研究。数值分析采用14-方程模型,为了提高数值分析精度,对传统特征线法进行了改进;给出了时间步长Δt和管道分段数N的确定条件,提出了自适应网格加密的方法来修正空间步长Δz以便更好地捕捉边界处液体压力的波动。对一典型液压系统的耦合振动响应进行了数值仿真,数值结果与实验结果吻合较好,证明了该方法的正确性和有效性。

梁的离散模型的模态反问题

采用集中质量法或有限差分法对梁进行离散,得到横向振动梁的弹簧-质点-刚杆模型,其质量矩阵为对角矩阵而刚度矩阵为对称五对角矩阵。已知(λ,x)、(μ,y)为梁的两个特征对,W为梁的总质量,考虑由这些数据构造离散梁的物理参数,得出构造具有正的质量和刚度的真实系统的充要条件。如果这些条件得到满足,则系统的构造是唯一的。因为构造梁的离散系统需要的数据可由测试得到,所以其结果适用于模态分析应用。