改进频率相关黏性阻尼模型的时程计算方法

孙攀旭 杨红 赵雯桐 刘庆林

摘   要：復阻尼模型的时域计算结果不能稳定收敛. 迟滞阻尼模型存在能量耗散与实际不符以及非线性的缺陷. 针对复阻尼模型和迟滞阻尼模型的缺陷，本文依据频域转化原则得到了频率相关黏性阻尼模型. 为实现结构体系的时程计算，基于加速度与位移的关系假定，进一步得到了改进频率相关黏性阻尼模型.改进频率相关黏性模型保留了结构每周期耗散能量与外激励频率无关的优点，同时克服了迟滞阻尼模型中能量耗散与实际不符的缺陷，还保证了单一振动频率下单自由度结构的线性特征.假定时间步长内结构处于单一频率的简谐振动，引入常平均加速度法，提出了单自由度体系的时程计算方法. 在此基础上，结合模态叠加法，推导了多自由度体系的时程计算公式. 算例结果表明，改进频率相关黏性阻尼模型可克服复阻尼模型频域法的缺陷，同时有效避免复阻尼模型时域法计算结果的发散现象.

关键词：复阻尼;频率相关黏性阻尼;时程计算;稳定收敛;线性特征

中图分类号：TU311.3                           文献标志码：A

文章编号：1674—2974（2020）09—0113—07

Abstract：There is divergent phenomenon existing in time-domain calculation results based on complex damping model. Hysteretic damping model has the shortcomings that energy dissipation is not consistent with the practical case and a defect of nonlinearity in linear elastic stage. To overcome the above shortcomings of complex damping model and hysteretic damping model， a frequency dependent viscous damping model is obtained based on the principles of frequency domain transformation in this paper. For the realization of structural time-history calculation method， an improved frequency dependent viscous damping model is obtained based on the assumption of the relationship  between acceleration and displacement. The improved frequency dependent viscous damping model has the advantage that energy consumption is  independent of external excitation frequencies. At the same time， the energy consumption in  the proposed model is consistent with the practical case， and it maintains linear characteristic of single-degree-of-freedom structure with a single vibration frequency. It is assumed that structural response is harmonic vibration response in every time step. It contains the single frequency. By introducing the constant average acceleration method， a time-history calculation method of single degree of freedom system can be put forward. On this basis， combined with the modal superposition method， the time-history calculation formulas of multi-degree of freedom system are obtained. The analysis results of the cases show that improved frequency viscous damping model can overcome the shortcoming of the frequency-domain method based on complex damping model. It can also avoid the divergent phenomenon in calculation results of time-domain method based on the complex damping model.

Key words：complex damping;frequency dependent viscous damping;time-history calculation;stable convergence;linear characteristic

目前应用最为广泛的阻尼模型是黏性阻尼模型，黏性阻尼模型由于其数学处理上的简易性，在结构动力计算中得到了广泛应用，但其存在每周期耗散能量与外激励频率相关的特点[1-4]，这与大部分材料在实验中每周期耗散能量与外频率无关的现象不一致[5].

能够更真实地描述实验现象的复阻尼模型具有体系每周期耗散能量与外激励频率无关的优点[6]，但其自由振动运动方程的通解中存在发散项，造成复阻尼模型时程迭代计算结果不收敛[7]. 潘玉华等[8]提出了复阻尼运动方程的高斯精细积分法，吴泽玉等[9]提出了复阻尼运动方程的增维精细积分法，Yuan等[10]提出了复阻尼模型的迭代更新法，但上述方法均无法避免复阻尼模型的时域发散现象.针对复阻尼模型的缺陷，周正华等[11]依据线弹性体的时域本构关系，采用最小二乘法使其近似等于频域内的复阻尼本构关系，得到一种时域复阻尼本构方程，但该方程仅适用于给定的频率范围，且用于多自由体系计算时，其误差有待进一步分析. Reggio等[12]采用Maxwell-Wiechert本构模型，在频域范围内近似等效于复阻尼本构模型，得到在时域内稳定收敛的运动方程，但计算过程过于复杂. Wang[13]在频域内采用Rayleigh阻尼矩阵等效复阻尼矩阵，进而避免发散现象，但存在计算结果不唯一、合理性不易判定的缺点.

假定阻尼力与结构体系的位移成正比，且与速度的方向相反，可进一步得到迟滞阻尼模型[14，15]. 迟滞阻尼模型保留了耗散能量与外激励频率无关的优点，同时时域计算结果稳定收敛，但其对应的阻尼体系在线弹性阶段具有非线性的缺陷[16，17].

如何解决上述阻尼模型的缺点，是提高结构动力计算结果可靠性的难点之一. 针对复阻尼模型，本文采用频域转化原则得到频率相关黏性阻尼模型，同时基于加速度与位移的关系假定，进一步将其改进为适用于迭代计算的改进频率相关黏性阻尼模型.在此基础上，结合常平均加速度法和模态叠加法，提出了基于改进频率相关黏性阻尼模型的多自由度体系时程计算方法.

1   基于加速度与位移关系的改进频率相关黏

性阻尼模型

1.1   运动方程的构建

单自由度体系的复阻尼运动方程为

方程（10）为改进频率相关黏性阻尼模型的时域运动方程.

1.2   不同阻尼模型的对比分析

迟滞阻尼模型[14，15]（以下称PVS）可克服复阻尼模型的时域发散现象，其阻尼力为

由式（12）可知，PVS具有阻尼力每周期耗散的能量与外激励频率无关的优点，但稳态反应时一个周期内阻尼力做功与实际耗散能量不相等，PVS存在阻尼力消耗能量与实际耗散能量不符的缺陷.

改进频率相关黏性阻尼模型（以下称PVJ）下，由方程（10）可得阻尼力为

由式（14）可知，PVJ在稳态反应时，阻尼力每周期耗散的能量与外激励频率无关，且一个周期内阻尼力消耗能量与实际耗散能量相等.因此PVJ可有效克服PVS阻尼力消耗能量与实际耗散能量不符的缺陷.

当结构处于简谐振动时，其位移响应如式（6）所示，可计算出两种阻尼模型下阻尼力与位移响应的变化关系如图1所示. PVS中阻尼力与位移的关系呈现非连续变化（图1（a）），结构体系具有非线性特点;PVJ中阻尼力随位移连续变化（图1（b）），结构体系具有线性特点. 线弹性状态下，单一振动频率的单自由度结构具有线性特征[19]，但仅PVJ呈现出线性特点，PVS则存在理论误差.

综上，PVJ不仅克服了PVS中阻尼力消耗能量与实际不符的缺陷，还保证了单一振动频率下单自由度结构的线性特征.

2   基于改进频率相关黏性阻尼模型的时程计算方法

2.1   单自由度体系的时程计算方法

按照时间步长Δt对时间进行离散，任意时刻可表示为tk = kΔt（k = 0，1，2…）. 时间步长较小时，假定tk时刻到tk + 1时刻体系振动为简谐振动，位移响应为

式中：I（t）为瞬时振幅，θk为tk时刻到tk + 1时刻的瞬时频率，φk为tk时刻的瞬时相位.

2.1.1   单自由度体系运动方程初值

假定初始时刻t0时，体系的振动频率为有阻尼自由振动频率，初值的确定需要首先计算结构的有阻尼自由振动频率.

频率相关黏性阻尼模型下单自由度体系运动方程对应的特征方程为

2.1.2   基于常平均加速度法的时程计算方法

2.2   多自由度体系的时程计算方法

2.2.1   基于模态叠加法的时程计算方法

多自由度体系的PVJ运动方程为

对于单一材料体系，ηK为比例矩阵，满足经典阻尼条件，方程（26）可直接采用模态叠加法[21]进行计算.

2.2.2   多自由度体系运动方程初值的确定

3   算例分析

3.1   单自由度体系

由式（18）可知，频率相关黏性阻尼运动方程两个特征根的实部均为负实数，因此其通解中仅含衰减项，而没有发散项，自由振动过程将是穩定收敛的.

以质量为1 000 kg，刚度为16 000 N/m，损耗因子为0.1的单自由体系为例，其初始位移为5 cm，初始速度为8 cm/s，采用PVJ计算体系的自由振动响应，所得结果如图2所示，该算例计算结果表明PVJ计算的位移响应稳定收敛，可有效解决复阻尼模型不能计算结构自由振动响应的问题.

体系初始时刻处于静止状态，分别采用PVJ和复阻尼模型时域计算方法（FZ）计算在迁安波地震作用下的位移响应时程，并与复阻尼模型的频域计算方法（FFZ）的计算结果进行对比（如图3所示），其中FFZ的计算结果可视为精确解[8]. 当地震作用持续时间大于15 s时，FZ计算的位移响应开始明显发散，因此基于复阻尼模型的时域计算方法仅能适用于地震作用持续时间较小的情况.当地震作用持续时间小于12 s时，PVJ、FZ和FFZ计算的位移响应近似相等.整个时程过程中，PVJ计算的位移响应不受地震持时的限制，位移响应一直稳定收敛.PVJ的位移在8.15 s处达到最大值，FFZ的位移在8.18 s处达到最大值，且相对误差为6.25%（见表1）.因此，PVJ与FFZ的计算结果近似相等，证明了PVJ的正确性.

FFZ为频域计算方法，仅能计算初始时刻为静止状态的结构动力响应.相比FFZ，PVJ为时域计算方法，可计算结构自由振动响应和地震作用下非静止初始状态的结构动力响应，克服了FFZ的缺陷.

3.2   多自由度体系

如图4所示，以3层剪切型钢筋混凝土框架结构为例，其损耗因子为0.1.

框架结构的初始处于静止状态，分别采用PVJ、FZ和FFZ计算结构在El Centro波和天津波作用下的地震反应，所得位移时程如图5所示. 在图5（a）、图5（b）中，当地震作用持续时间分别小于12 s、15 s时，PVJ的计算结果与FZ的计算结果一致.但是，随着地震作用持时增加，FZ的数值计算结果出现发散现象，这是复阻尼模型自由振动方程通解中包含发散解所引起的.相比之下，PVJ一直具有稳定收敛的优点. El Centro波作用下，PVJ和FFZ的位移均在4.2 s处达到最大值，且相对误差为6.60%;PVJ的位移在8.77 s处达到最大值，FFZ的位移在8.79 s处达到最大值，且相对误差为5.20%（见表2）. PVJ与FFZ的计算结果近似相等，进一步证明了本文提出的改进频率相关黏性阻尼时程计算方法的正确性.

4   结   论

经理论推导和算例分析，得到以下结论：

1）引入加速度与位移关系假定，提出了基于复阻尼模型的改进频率相关黏性阻尼模型，可有效克服黏性阻尼模型中耗散能量与外激励频率相关的缺陷.

2）与迟滞阻尼模型相比，基于加速度和位移关系的改进频率相关黏性阻尼模型克服了迟滞阻尼模型中能量耗散与实际不符的缺陷，还保证了单一振动频率下单自由度结构的线性特征.

3）结合常平均加速度法和模态叠加法，推导了改进频率相关黏性阻尼运动方程时程计算的相关公式，算例分析表明，改进频率相关黏性阻尼模型可克服复阻尼模型频域法的缺陷，同时有效避免复阻尼模型时域法计算结果的发散现象.

参考文献

[1]    FERIANI A，PEROTTI F. The formation of viscous damping matrices for the dynamic analysis of MDOF systems [J]. Earthquake Engineering and Structural Dynamics，1996，25（7）：689—709.

[2]    MASTRODDI F，EUGENI F，ERBA F. On the modal diagonalization of viscoelastic mechanical systems[J]. Mechanical Systems and Signal Processing，2017，96：159—175.

[3]    朱鏡清. 结构抗震分析原理 [M]. 北京：地震出版社，2002：62—66.

ZHU J Q. Seismic analysis of structures[M]. Beijing：Earthquake Press，2002：62—66. （In Chinese）

[4]    董云，楼梦麟. 基于结构基频确定Rayleigh阻尼系数的优化方法及其讨论 [J]. 湖南大学学报（自然科学版），2014，41（2）：8—13.

DONG Y，LOU M L. An optimization solution for Rayleigh damping coefficients based on the fundamental frequency of structure [J]. Journal of Hunan University（Natural Sciences），2014，41（2）：8—13. （In Chinese）

[5]    BERT C W. Material damping：an introductory review of mathematic measures and experimental technique [J]. Journal of Sound and Vibration，1973，29（2）：129—153.

[6]    张辉东，王元丰. 复阻尼模型结构地震时程响应研究[J]. 工程力学，2010，27（1）：109—115.

ZHANG H D，WANG Y F. Study on seismic time-history response of structures with complex damping [J]. Engineering Mechanics，2010，27（1）：109—115. （In Chinese）

[7]    朱敏，朱镜清. 逐步积分法求解复阻尼结构运动方程的稳定性问题[J]. 地震工程与工程振动，2001，21（4）：59—62.

ZHU M，ZHU J Q. Studies on stability of step-by-step methods under complex damping conditions [J]. Earthquake Engineering and Engineering Vibration，2001，21（4）：59—62. （In Chinese）

[8]    潘玉华，王元丰. 复阻尼结构动力方程的高斯精细时程积分法[J]. 工程力学，2012，29（2）：16—20.

PAN Y H，WANG Y F. Gauss precise time-integration of complex damping vibration systems [J]. Engineering Mechanics，2012，29（2）：16—20. （In Chinese）

[9]    吴泽玉，王东炜，李玉河. 复阻尼结构动力方程的增维精细积分法[J]. 振动与冲击，2017，36（2）：107—110.

WU Z Y，WANG D W，LI Y H. Magnified dimension precise integration method for the dynamic equations of complex damped structures [J]. Journal of vibration and Shock，2017，36（2）：107—110. （In Chinese）

[10]  YUAN Y G，ZUO K Z，ZHANG T Z. An efficient iterative updating method for hysteretic damping models [J]. Applied Mathematics and Computation，2017，320：86—98.

[11]  周正华，廖振鹏，丁海平. 一种时域复阻尼本构方程[J]. 地震工程与工程振动，1999，19（2）：37—44.

ZHOU Z H，LIAO Z P，DING H P. A time-domain complex-damping constitutive equation [J]. Earthquake Engineering and Engineering Vibration，1999，19（2）：37—44. （In Chinese）

[12]  REGGIO A，ANGELIS M D. Modelling and identification of structures with rate-independent linear damping [J]. Meccanica，2015，50（3）：617—632.

[13]  WANG J. Rayleigh coefficients for series infrastructure systems with multiple damping properties [J]. Journal of Vibration and Control，2015，21（6）：1234—1248.

[14]  CLOUGH R W，PENZIEN J. 结构动力学[M]. 王光远，译. 北京：科学出版社，1983：27—46.

CLOUGH R W，PENZIEN J. Dynamics of Structures[M]. WANG G Y （trans）. Beijing：Science Press，1983：27—46. （In Chinese）

[15]  CHEN L Y，CHEN J T，CHEN C H，et al. Free vibration of a SDOF system with hysteretic damping[J]. Mechanics Research Communications，1994，21（6）：599—604.

[16]  INAUDI J A，KELLY J M. Linear hysteretic damping and Hilbert transform[J]. Journal of Engineering Mechanics，1995，121（5）：626—632.

[17]  孙靖雅，华宏星，肖锋，等. 非线性迟滞阻尼对隔振系统力传递特性影响[J]. 振动与冲击，2014，33（10）：131—136.

SUN J Y，HUA H X，XIAO F，et al. Influence of nonlinear hysteretic damping on force transmissibility of a vibration isolation system[J]. Journal of Vibration and Shock，2014，33（10）：131—136. （In Chinese）

[18]  朱鏡清. 频率相关黏性阻尼理论及有关问题的解[J]. 振动与冲击，1992，11（4）：1—7.

ZHU J Q. Frequency dependant viscous damping theory and some related problems [J]. Journal of Vibration and Shock，1992，11（4）：1—7. （In Chinese）

[19]  刘晶波，杜修力. 结构动力学[M].北京：机械工业出版社，2005： 1—9.

LIU J B，DU X L. Dynamics of structures [M]. Beijing：China Machine Press，2005：1—9. （In Chinese）

[20]  李鸿晶，王通，廖旭. 关于 法机理的一种解释[J]. 地震工程与工程振动，2011，31（2）：55—62.

LI H J，WANG T，LIAO X. An interpretation on Newmark beta methods in mechanism of numerical analysis [J]. Journal of Earthquake Engineering and Engineering Vibration，2011，31（2）：55—62. （In Chinese）

[21]  CAUGHEY T K. Classical normal modes in damped linear dynamic systems[J]. Journal of Applied Mechanics，1960，27（3）：269—271.