风力发电塔的自振特性分析

将风力发电塔视为带有附加质量的变截面悬臂梁,进行横向振动的自振特性分析.采用直接模态摄动法建立风力发电塔自振特性的近似求解方法,与采用梁单元模型和壳单元模型的有限元法的计算结果相比较.数值计算结果表明直接模态摄动法具有较好的精度,形成了半解析解形式.     

模态分析与动态子结构方法新进展

综述在模态分析与动态子结构方法研究的一些最新进展.首先回顾经典的位移展开定理和模态叠加原理.为了加速经典方法的收敛速度、提高计算效率, 进一步介绍两个新的结构位移展开定理(采用固定界面模态的位移展开新定理, 给出采用低阶固定界面模态的高精度位移展开式;采用混合模态的位移展开新定理, 给出采用低阶混合模态表示的高精度位移展开式)和相应的动力学新解法.相应上述3个位移展开定理, 介绍采用解析推导的方法构造出3类动态精确子结构方法, 各种子结构模态综合法实质上都是它们的某种近似与变化形式, 从而形成系统的动态子结构分析技术.上述介绍的模态分析与动态子结构方法新进展与经典模态分析技术一起形成结构动力学分析技术的系统理论.      

求解非比例阻尼体系复模态的实模态摄动法

根据工程结构的实际情况,建立了非比例阻尼结构体系复模态特性的近似求解方 法------实模态摄动法. 这一方法以复Ritz向量展开原理为基础, 把非比例阻尼结构体系复模态特性的分析过程分解为两个基本步骤,首先以结构体系的实模 态向量构建复Ritz向量的求解子空间,然后通过非线性复代数方程组的求解代替扩阶后的复 特征值方程的求解,从而简化了计算过程.通过两个算例表明: 这一方法不仅计算简便，而且具有较高的计算精度和执行效率, 对于复杂的非比例阻尼系统是很适用的，具有一定的工程应用价值.     

多自由度系统复模态理论的摄动方法 (三)非对称系统复特征值问题的摄动迭代法

基于实模态理论和相应特征值计算方法,本文在位形空间中用摄动方法求解非对称系统的复特征值和复模态。为提高精度,增加了一个迭代格式。最后给出算例并讨论了计算精度和收敛性。     

工程结构的随机特征问题研究及其在梁结构中的应用

采用子结构模态综合和摄动随机有限元相结合求解工程结构的随机特征问题。为求出随机特征对的方差，借助于模态截断概念推出诸特征值与特征向量对随机变量的偏导数。以染结构为典型算法，定量研究了子结构动模态的选取个数与随机特征对的计算精度间关系，以梁的长细比首次确定使用Timoshenko梁和Euler-Bernoulli梁两模型求解梁类工程结构随机特征问题的适用范围。     

核反应堆主系统管道的动力特性研究

本文系统研究了核反应堆主系统管道在局部破裂工况下的动力特性及非线性动力响应。采用统一动力子结构方法建立复杂管系的动力计算模型,采用Lanczos法和Arnoldi法分别求解系统的实、复特征值问题,分析了流速大小对系统固有频率的影响,在给定破口处的流体排放力后,采用Lanczos基矢量与模态综合相结合的自由度减缩技术,兼用伪力法、摄动法等技巧,求解管道系统的非线性动力响应。计算表明,采用本文的理论方法所研制的软件对核安全设计审批具有实用价值。文中给出了主系统管道的若干算例。     

混合界面直接分支模态综合法

在直接分支模态综合法(本质上属自由界面模态综合法)的基础上,提出并深入研究了混合界面直接分支模态综合法。混合界面法将复杂庞大的结构取做主体部件,而将与此主体部件相连的其余部分视为分支部件。对主体部件采用自由界面,分支部件则采用固定界面。各子结构之间完全独立计算,并且通过计入主、从子结构的剩余模态矩阵的办法来补偿所舍去的高阶模态分量,保证了计算精度。本方法从工程应用的实际出发,通过模拟航空发动机机匣管路系统的数值算例,验证了混合界面法的正确性。     

一种计算非比例阻尼结构地震响应的新方法

在非比例阻尼结构地震响应分析中, 直接积分法和强迫解耦法均具有鲜明的优缺点, 考虑到计算精度和计算效率的均衡, 提出了一种可用于非比例阻尼结构地震响应计算的新方法---多自由度模态方程方法. 同时, 通过推导指出了直接积分法和强迫解耦法是所提出方法的两种特殊形式, 从而构建了非比例阻尼结构地震响应计算方法的完整理论体系. 最后采用数值算例验证了多自由度模态方程方法在非比例阻尼结构地震响应计算中的有效性, 并说明其可以通过合理划分结构分区来调节计算精度和计算效率.      

基于键合图法的多个子结构自由界面模态综合法研究和应用

将键合图方法用于动态子结构研究,提出了一种求解多个子结构自由界面模态综合法的新方法。通过一系列物理和数学上的分析,详细推导出多个子结构自由界面模态综合法的计算过程。在本文给出的算例中,基于键合图方法的自由界面模态综合法,通过建立各个子结构的状态空间方程,计算子结构相应的特征值矩阵和振型矩阵;这些子结构在进行模态综合后,获得整体结构的特征值矩阵和振型矩阵,该结果和原来整体结构完全相同,这进一步说明了本文提出的新方法的正确性。运用本文提出的方法建立状态空间方程,在子结构的模态综合过程中,不会产生系统特征值的增失根问题,确保了子结构在综合后其整体结构模态信息的完整性。     

子结构Ritz向量综合法

讨论了在动态子结构法中应用静态Ｒｉｔｚ向量替代模态向量时所应考虑的基本问题，即如何合理地定义荷载空间分布向量｛ｆ（ｓ）｝，以避免在子结构中出现低阶模态不连续分布的情况。文中给出了在自由子结构和约束子结构内进行Ｒｉｔｚ变换的一般形式及修正公式。算例表明：应用子结构Ｒｉｔｚ向量综合技术能有效地提高线性结构模态分析的计算效率     

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 approximatesolution of the complex modal quadratic generalized eigenvalue problem of viscouslydamped linear vibration systems.First,the damping matrix is decomposed into the sum of aproportional-and a nonproportional-damping parts,and the solutions of the real modaleigenproblem with the proportional dampings are determined,which are a set of initialapproximate solutions of the complex modal eigenproblem.Second,by taking thenonproportional-damping part as a small modification to the proportional one and using thematrix perturbation analysis method,a set of approximate solutions of the complex modaleigenvalue problem can be obtained analytically.The result is quite simple.The new methodis applicable to the systems with viscous dampings-which do not deviate far away from theproportional-damping case.It is particularly important that the solution technique be alsoeffective to the systems with heavy,but not over,dampings.The solution formul     

A simplified frequency equation and its approximate solution of a beam with an incipient crack from a wave perspective

This paper presents a simplified frequency equation and its approximate solution to predict the modal frequencies of a beam with an incipient crack. The physical implication of the simplified frequency equation is fully described from a wave perspective for the cracked beam with arbitrary support conditions. The approximate solution of the proposed frequency equation is derived from a wave perspective as well. The asymptotic equivalence is demonstrated between the approximate solution and that obtained by the first order perturbation method as the mode number increases. The validity of the proposed approach is demonstrated through comparison to both numerical results from finite element analysis and experimental data.     

计算特征向量摄动量的混合基展开法

在结构修改和模型校正中，模态展开法是计算特征向量摄动量的常用方法之一，但当高阶模态被截断时，它会带来很大的截断误差。本文利用已知的有限阶模态，构造了Ｎ维欧氏空间的一个新基－混合基，并将特征向量的摄动量在新基上展开来计算特征向量的一、二阶 摄动量。该方法使得不管截模态个数的多少，其精度总与全模态展开法相同，且计算量都远少于全模记展开法；与改进的部分民开法相比，本方法不要求所截留的模态边连续的低阶模态     

An identification method of generalized modal parameters of linear vibrating defective systems

Based on the generalized mode theory of linear vibrating defective systems, an identification method of generalized modal parameters is presented in this paper. By the use of this method, which combines the direct method with the iteration method in frequency domain, all the generalized modal parameters can be identified without any initial value. It is shown that the present method is effective and useful.     

线性振动亏损系统广义模态参数的识别方法

基于线性振动亏损系统的广义模态理论,本文提出了一种识别亏损系统广义模态参数的频域方法。该方法将直接法与迭代法相结合,无需人工初值,可分步识别出亏损系统的全部广义模态参数。模拟识别结果表明,本文方法是有效且可行的。     

局部损伤梁动力问题的近似计算方法

应用模态摄动法求解局部损伤梁动力特性。这一方法是在无损伤梁的模态子空间内实施,将含局部损伤梁的微分方程的求解转化为代数方程求解。通过算例,表明这一方法可有效的简化计算,同时计算结果具有较高的精度。     

Iterative methods for scattering problems in isotropic or anisotropic elastic waveguides

We consider the time-harmonic problem of the diffraction of an incident propagative mode by a localized defect, in an infinite elastic waveguide. We propose several iterative algorithms to compute an approximate solution of the problem, using a classical finite element discretization in a small area around the perturbation, and a modal expansion in the unbounded straight parts of the guide. Each algorithm can be related to a so-called domain decomposition method, with an overlap between the domains. Specific transmission conditions are used, so that at each step of the algorithm only the sparse finite element matrix has to be inverted, the modal expansion being obtained by a simple projection, using a bi-orthogonality relation. The benefit of using an overlap between the finite element domain and the modal domain is emphasized. An original choice of transmission conditions is proposed which enhances the effect of the overlap and allows us to handle arbitrary anisotropic materials. As a by-product, we derive transparent boundary conditions for an arbitrary anisotropic waveguide. The transparency of these new boundary conditions is checked for two- and three-dimensional anisotropic waveguides. Finally, in the isotropic case, numerical validation for two- and three-dimensional waveguides illustrates the efficiency of the new approach, compared to other existing methods, in terms of number of iterations and CPU time.     

基于精细积分技术的非线性动力学方程的同伦摄动法

将精细积分技术（PIM）和同伦摄动方法（HPM）相结合,给出了一种求解非线性动力学方程的新的渐近数值方法。采用精细积分法求解非线性问题时,需要将非线性项对时间参数按Taylor级数展开,在展开项少时,计算精度对时间步长敏感;随着展开项的增加,计算格式会变得越来越复杂。采用同伦摄动法,则具有相对筒单的计算格式,但计算精度较差,应用范围也限于低维非线性微分方程。将这两种方法相结合得到的新的渐近数值方法则同时具备了两者的优点,既使同伦摄动方法的应用范围推广到高维非线性动力学方程的求解,又使精细积分方法在求解非线性问题时具有较简单的计算格式。数值算例表明,该方法具有较高的数值精度和计算效率。     

Static and Dynamic Analysis of Non-Linear Uncertain Structures

The procedures usually adopted in the evaluation of the stochastic response of structures with uncertain parameters, the so-called stochastic structures, are affected by some limits. Namely, the major drawbacks are: the conspicuous computational time required for the analysis of many degree of freedom systems and the loss of accuracy in the case of large uncertainty in the parameters. A method able to reduce the previous inconveniences was recently introduced in the study of statically loaded linear structures. The method, named improved perturbation method, was also extended to the field of linear dynamics. In both cases, the improved perturbation method provides a good approximation, even in the case of moderately large deviation of the uncertain parameters, and the computational time required is comparable to conventional first order perturbation. The present paper intends to apply the improved perturbation method in the second order analysis of geometrically non-linear uncertain systems subjected to static and dynamic deterministic forces.