基于微分求积法求解非粘滞阻尼系统的时程响应

研究一种含有指数型非粘滞阻尼线性多自由度振动系统的时程分析问题。该非粘滞阻尼模型假设阻尼力与质点速度的时间历程相关,数学表述为质点速度与核函数的卷积。由于阻尼模型的改变,常用的数值积分方法（如Newmark-β法、Wilson-θ法）不能直接应用于这种非粘滞阻尼系统。基于一种无条件稳定的微分求积方法,给出了这种非粘滞阻...     

非粘滞阻尼系统时程响应分析的精细积分方法

考虑一个具有非粘滞阻尼特性的多自由度系统响应的时程分析问题.该非粘滞阻尼模型假设阻尼力与质点速度的时间历程相关,数学表达式体现为阻尼力等于质点速度与某一核函数的卷积.在利用状态空间方法将系统运动方程转换成一阶的状态方程的基础上,采用精细积分方法对状态方程进行数值求解,得到一种求解该阻尼系统时程响应的精确、高效的计算方法.通过两个数值算例表明,采用该方法得到几乎精确的数值计算结果,而且计算效率有成数量级的提高.     

基于Volterra积分方程Taylor展开的卷积型非粘滞阻尼系统响应的时程分析方法

考虑一个具有卷积型非粘滞阻尼特性的多自由度系统响应的时程分析问题。非粘滞阻尼模型假设阻尼力与质点速度的时间历程相关,数学表达式为阻尼力等于质点速度与某一核函数的卷积,该模型为常用的粘滞阻尼模型的一般化形式。以一种在特定区间内求解第二类Volterra方程的Taylor展开法为基础,对所分析时段中各时间点的响应函数逐步作Taylor展开,代入卷积核来消去运动方程的积分项,通过求解推导出的时变线性方程组完成对卷积型阻尼模型系统的时程响应分析。数值算例验证了该方法的有效性。该方法增大时间步长可以大幅减少计算量,计算精度有所下降。     

一种结构动力响应分析的快速高阶算法

为了提高基于高阶格式的结构动力响应微分求积分析方法的计算效率,发展了一种求解动力方程的快速算法．利用微分求积原理将结构动力方程转化为标准Sylvester方程的形式,通过对系数矩阵进行矩阵分解,进而将动力响应Sylvester方程化为一系列标准线性方程组,采用相关成熟算法求解这些线性方程组后即可获得结构动力时程响应的全部解答．结构动力响应微分求积分析方法为高阶数值方法,一步计算可以获得多个时点处的动力响应．基于本文快速算法,不必直接对矩阵方程进行求解．数值算例表明,本文快速算法能够准确地计算出结构动力响应,具有数值精度高、收敛性好的优点．     

基于Kanai-Tajimi谱的非黏滞阻尼结构地震动响应的简明闭式解

针对非黏滞阻尼结构基于Kanai-Tajimi谱的卷积-微分混合动力方程解法较繁琐的问题,提出了 一种新的简明封闭解法.非黏滞阻尼模型能较好地模拟实际工程材料的阻尼特征,常以指数型核函数的卷积形式表示,给出其对应的微分型本构关系.Kanai-Tajimi谱随机地震动模型能较好地描述场地的随机地震动特性,其工程应用时所获...     

附加自由阻尼结构系统的有限元分析和拓扑优化设计方法研究

给出了一种适用于车身等复杂附加阻尼结构系统的有限元建模和动态特性分析方法,包括附加自由阻尼薄壁结构和附加阻尼材料粘弹性特性的有限元建模方法以及附加阻尼结构动态特性的有限元分析方法。在此基础上,又以强迫振动响应最小为优化设计目标,给出了一个附加阻尼结构的拓扑优化设计方法,包括优化设计问题的列式和适用的求解算法。通过对一个简单薄壁构件和一个车身地板上的附加阻尼结构的拓扑优化设计,验证了提出的有限元分析和拓扑优化设计方法的有效性。     

结构阻尼时域本构模型及其应用

阻尼合金作为一种新型结构功能材料在不少领域已获应用，由于其阻尼值较大且随频率呈复杂交化关系，传统的线性粘性阻尼理论或经典的非频变结构阻尼理论难以精确地描述其耗能行为。本文应用粘弹性阻尼理论，根据阻尼合金储能模量和损耗因子在频域的实测数据．应用最优化方法拟合出标准线性体模型中的本构参数；根据积分形式的三参量本构关系和变形体虚功原理，推导出了有限元形式的动力学方程；讨论了三参数初值的选取；对包含卷积积分的有限元动力学方程通过数学推导将其化为三阶线性微分方程组，再转化为标准状态变量方程，应用数值求解。数值计算实例证明了所提方法的正确和有效性。     

结构阻尼模型及系统时域动响应

振动界对经典的非频变结构阻尼模型一直有种种争议,本文检查了它的前提并指出将该模型推广到系统时域动响应分析时导致的悖论.为了提供结构阻尼系统的时域动响应实用分析方法,文中分别提出了基于响应优势谱成份的粘性阻尼近似模型和基于粘弹性理论的三参数频变结构阻尼模型;它们既可逼近结构阻尼缓频变的特性,又可保证系统的时域动响应具有因果性.     

串联电气设备支架隔震体系地震响应半解析法

通过并联橡胶隔震支座,建立串联高架电气设备支架隔震体系及力学模型,应用分布参数梁振动理论, 通过边界条件引入集中参数,推导其频率方程,并用数值方法求得频率及振型. 应用Betti定律,推导具有集中分布参数柔性节点的多节电瓷型高压电气设备的正交条件,得到广义质量及广义载荷. 将该串联隔震体系的非比例阻尼分解为比例阻尼部分和非比例阻尼部分,应用Hamilton原理推导出非比例阻尼部分等效振型阻尼比,实现串联电气设备支架隔震体系振动方程的解耦,然后通过振型叠加法求得结构的地震响应. 最后应用该半解析法与有限元数值积分法对一330kV电压互感器地震响应进行分析. 结果表明:该隔震体系能显著减小电气设备的地震响应,同时该半解析法求解的地震响应与有限元数值积分求解的结果相一致,说明该隔震体系的有效性与该半解析方法求解非比例阻尼串联电气设备支架隔震体系地震响应的正确性.      

液体边界层对平板结构振动阻尼效应分析

本文分析了液体边界层对平板振动的阻尼效应，基于液体的Stokes边界层理论，建立了考虑液体边界层阻尼效应的平板振动单自由度系统运动方程，利用Laplace变换方法求解了运动方程，得到了平板结构的位移响应解答，制作了一个单自由度质量-弹簧系统来测量水体粘滞阻尼效应，并将理论与试验结果进行了对比分析，表明Stokes边界层理论能较好模拟水体的粘滞阻尼特性，边界层阻尼明显减小了结构的动力响应，本文最后讨论了所提的分析方法的适用范围．     

航天器噪声试验中结构振动响应预示方法研究

航天器在随运载火箭发射过程中要承受严酷的噪声环境,需通过噪声试验来检验航天器承受噪声环境并能正常工作的能力.航天器噪声试验中结构振动的响应特性是结构强度设计应该考虑的因素之一,更是制定器上组件随机振动试验条件的重要依据,因此有必要在航天器研制初期对噪声载荷作用下的结构振动进行响应预示.文章应用商用有限元分析软件MSC.Patran和MSC.Nastran建立了某型号航天器结构舱板的有限元模型,将噪声载荷声压谱转换为脉动压力功率谱密度,进而采用模态法分析结构在噪声载荷作用下的随机振动响应,并将仿真预示结果与试验结果进行对比研究,在仿真分析中考虑阻尼参数模型和流场附加质量效应等因素的影响;通过研究表明:采用阻尼比随频率提高而减小的经验阻尼参数模型可以较好地反映中高频响应特性、得到较为准确的总均方根响应分析结果,进一步采用虚拟质量法考虑流场附加质量效应可以得到较为准确的功率谱密度响应分析结果.文章提出的仿真分析方法建模简便、计算成本低,适用于在航天器研制初期对航天器噪声试验中的结构振动进行响应预示.      

有阻尼体系的时域边界元法

针对时域边界元方法的特点，提出了一种新的阻尼模型（称为比例时耗阻尼），并成功地应用于时域边界元方法中．这是在时域边界元方法中首次考虑阻尼这一影响结构动力响应的重要因素，为时域边界元在实际工程中的应用解决了一个难题．通过对简单问题的分析和计算，验证了本文模型的正确性．     

Unified two-phase nonlocal formulation for vibration of functionally graded beams resting on nonlocal viscoelastic Winkler-Pasternak foundation

A nonlocal study of the vibration responses of functionally graded (FG) beams supported by a viscoelastic Winkler-Pasternak foundation is presented. The damping responses of both the Winkler and Pasternak layers of the foundation are considered in the formulation, which were not considered in most literature on this subject, and the bending deformation of the beams and the elastic and damping responses of the foundation as nonlocal by uniting the equivalently differential formulation of well-posed strain-driven (ε-D) and stress-driven (σ-D) two-phase local/nonlocal integral models with constitutive constraints are comprehensively considered, which can address both the stiffness softening and toughing effects due to scale reduction. The generalized differential quadrature method (GDQM) is used to solve the complex eigenvalue problem. After verifying the solution procedure, a series of benchmark results for the vibration frequency of different bounded FG beams supported by the foundation are obtained. Subsequently, the effects of the nonlocality of the foundation on the undamped/damping vibration frequency of the beams are examined.     

Structural damping identification method using normal FRFs

All structures exhibit some form of damping, but despite a large literature on the damping, it still remains one of the least well-understood aspects of general vibration analysis. The synthesis of damping in structural systems and machines is extremely important if a model is to be used in predicting vibration levels, transient responses, transmissibility, decay times or other characteristics in design and analysis that are dominated by energy dissipation. In this paper, new structural damping identification method using normal frequency response functions (NFRFs) which are obtained experimentally is proposed and tested with the objective that the damped finite element model is able to predict the measured FRFs accurately. The proposed structural damping identification is a direct method. In the proposed method, normal FRFs are estimated from the complex FRFs, which are obtained experimentally of the structure. The estimated normal FRFs are subsequently used for identification of general structural damping. The effectiveness of the proposed structural damping identification method is demonstrated by two numerical simulated examples and one real experimental data. Firstly, a study is performed using a lumped mass system. The lumped mass system study is followed by case involving numerical simulation of fixed–fixed beam. The effect of coordinate incompleteness and robustness of method under presence of noise is investigated. The performance of the proposed structural damping identification method is investigated for cases of light, medium, heavily and non-proportional damped structures. The numerical studies are followed by a case involving actual measured data for the case of a cantilever beam structure. The results have shown that the proposed damping identification method can be used to derive an accurate general structural damping model of the system. This is illustrated by matching the damped identified FRFs with the experimentally obtained FRFs.     

A dual-porosity model for gas reservoir flow incorporating adsorption behaviour—Part II. Numerical algorithm and example applications

In the present study, an algorithm is presented for the dual-porosity model formulated in Part I of this series. The resultant flow equation with the dual-porosity formulation is of an integro-(partial) differential equation involving differential terms for the Darcy flow in large fractures and integrals in time for diffusion within matrix blocks. The algorithm developed here to solve this equation involves a step-by-step finite difference procedure combined with a quadrature scheme. The quadrature scheme, used for the integral terms, is based on the trapezoidal method which is of second-order precision. This order of accuracy is consistent with the first- and second-order finite difference approximations used here to solve the differential terms in the derived flow equation. In an approach consistent with many petroleum reservoir and groundwater numerical flow models, the example formulation presented uses a first-order implicit algorithm. A two-dimensional example is also demonstrated, with the proposed model and numerical scheme being directly incorporated into the commercial gas reservoir simulator SIMED II that is based on a fully implicit finite difference approach. The solution procedure is applied to several problems to demonstrate its performance. Results from the derived dual-porosity formulation are also compared to the classic Warren–Root model. Whilst some of this work confirmed previous findings regarding Warren–Root inaccuracies at early times, it was also found that inaccuracy can re-enter the Warren–Root results whenever there are changes in boundary conditions leading to transient variation within the domain.     

捷联惯导系统减振设计

机抖激光陀螺的捷联惯性导航系统结构设计的关键在于惯性测量单元（IMU）的减振设计。文章描述了一种IMU减振系统设计的方法,并在实际应用当中获得良好的效果。IMU减振系统的设计方法主要分为三大内容。首先是对减振系统中减振器的布局方案进行仿真分析,建立减振系统的六自由度振动响应模型,研究减振器布局方案的频谱特性,比较得到减振效果较优的减振器方案。然后针对获得的减振器布局方案,设计IMU结构框架和各向同性的减振器,并采用有限元的方法,对IMU的减振系统进行模态分析,获得IMU减振系统的前20阶模态频率和振型,用有限元仿真分析的方法来验证减振设计。目前这种IMU减振系统的设计方法已应用到某新型捷联惯导系统当中,减少了系统陀螺精度的降低。     

圆形截面约束阻尼梁整体有限元建模与振动分析

为克服约束阻尼结构在工程应用中存在的由三维有限元建模造成的单元数目巨大的问题,本文提出了圆形实心截面梁附加约束阻尼层横向振动的整体有限元建模方法,并基于此模型,研究了材料的物理属性和几何因子对这类结构的固有频率的影响。通过与三维有限元解法相比较,证明了该方法的可行性与正确性,并给出了该方法的适用范围。     

A COMPLEX PARAMETER MODEL UPDATING METHOD FOR DAMPING CHARACTERISTICS OF FINITE ELEMENT MODELS

阻尼对于结构动力学响应具有重要的影响,但有限元模型一般很难对阻尼特性进行精确建模.基于实测频响函数,研究了一种有限元模型阻尼特性的复参数修正方法.以待修正区域各单元质量、刚度矩阵的比例修正系数为复修正参数,建立了单元矩阵比例修正的灵敏度方程直接算法,并对比分析了复修正参数与不同阻尼特性之间的数学关系.以六自由度集中参数模型和25杆平面桁架模型为例,验证了复参数修正方法在阻尼特性修正中的有效性.     

Time domain modeling of damping using anelastic displacement fields and fractional calculus

A fractional derivative model of linear viscoelasticity based on the decomposition of the displacement field into an anelastic part and elastic part is developed. The evolution equation for the anelastic part is then a differential equation of fractional order in time. By using a fractional order evolution equation for the anelastic strain the present model becomes very flexible for describing the weak frequency dependence of damping characteristics. To illustrate the modeling capability, the model parameters are fit to available frequency domain data for a high damping polymer. By studying the relaxation modulus and the relaxation spectrum the material parameters of the present viscoelastic model are given physical meaning. The use of this viscoelastic model in structural modeling is discussed and the corresponding finite element equations are outlined, including the treatment of boundary conditions. The anelastic displacement field is mathematically coupled to the total displacement field through a convolution integral with a kernel of Mittag–Leffler function type. Finally a time step algorithm for solving the finite element equations are developed and some numerical examples are presented.