孫攀旭 楊紅 劉慶林

摘 要：混合結構的阻尼矩陣不再滿足經典阻尼條件，無法直接采用模態(tài)疊加法。遲滯阻尼模型具有每周期耗散能量與外激勵頻率無關的優(yōu)點，且時域計算結果穩(wěn)定收斂，但不滿足能量守恒原則。利用結構每周期耗散能量與阻尼做功相等，對遲滯阻尼模型的阻尼系數進行修正，得到了改進遲滯阻尼模型，并提出了對應的混合結構模態(tài)疊加法。在此基礎上，分析改進遲滯阻尼模型的反應譜特點，計算最不利地震作用效應組合，提出了基于改進遲滯阻尼理論的混合結構抗震分析SRSS法。算例分析表明，改進遲滯阻尼模型計算的地震作用效應比粘性阻尼模型計算的地震作用效應更大，其增大幅度約為15%～20%。

關鍵詞：混合結構;遲滯阻尼;模態(tài)疊加法;反應譜;SRSS

中圖分類號：TU311.3

文獻標志碼：A? 文章編號：2096-6717（2019）02-0093-06

Abstract：Damping matrix of mixed structure no longer satisfies the classical damping condition. In such a case， the modal superposition method may become invalid. The advantage of hysteretic damping model is that energy consumption is related to external excitation frequencies in each cycle. At the same time， the calculation results of time-domain are stable. However， it does not meet the principle of energy conservation. By the equivalence of energy dissipation and damping power in each cycle， damping coefficient of hysteresis damping model can be modified. Subsequently， an improved hysteretic damping model and its own modal superposition method of mixed structure are proposed in this paper. The characteristics of improved hysteretic damping model response spectrum are firstly analyzed. Then， the most unfavorable combination of seismic effects could be calculated. Based on the improved hysteretic damping model， SRSS method for seismic analysis of mixed structure is performed. Results show that seismic effects of the improved hysteretic damping model are more significant than that of viscous damping model， and the increase is about 15%～20%.

Keywords：mixed structure; hysteretic damping; modal superposition method; response spectrum; SRSS

混合結構的地震響應計算方法與單一材料結構采用的常用方法不同。由于不同材料的阻尼比不同，混合結構的阻尼矩陣并不與質量矩陣或剛度矩陣保持比例關系，故混合結構無法直接采用基于粘性阻尼模型的實模態(tài)疊加法，抗震計算時也無法直接采用SRSS法。

目前，常用的阻尼模型為粘性阻尼模型和復阻尼模型[1]。對于混合結構，選擇結構的重要振型，可得到相應的Rayleigh阻尼矩陣，進而實現模態(tài)疊加法[2-4]，但其結果不具有唯一性，高階振型地震作用效應亦偏小[5];采用等效阻尼比的思路，將阻尼矩陣等價為比例矩陣，同樣可實現模態(tài)疊加法[6-7]，但計算出局部結構的地震效應偏小或偏大[8];汪夢甫[9]、Neugebauer等[10-11]利用狀態(tài)空間法，提出了基于粘性阻尼模型的復模態(tài)疊加法，但矩陣維度增加了一倍，計算量較大;劉慶林等[12]基于復阻尼模型的復模態(tài)疊加法可直接適用于混合結構，但計算較為復雜，且直接剔除復阻尼運動方程通解中發(fā)散項的做法存在理論缺陷[13]。Wang[14]采用Rayleigh阻尼矩陣等效復阻尼矩陣，Reggio等[15]將Maxwell-Wiechert本構模型等效于復阻尼本構模型，進而解決了發(fā)散問題，但上述方法得到的計算結果與復阻尼理論存在一定的誤差。遲滯阻尼模型計算結果唯一，且克服了粘性阻尼模型中每周期耗散能量與外激勵頻率相關的缺點。同時，遲滯阻尼模型的時域計算結果穩(wěn)定收斂，避免了復阻尼模型的發(fā)散現象[16]。

遲滯阻尼模型克服了粘性阻尼模型和復阻尼模型的缺陷，但不滿足能量守恒的原則。本文利用結構每周期耗散能量與阻尼做功相等的特性，對遲滯阻尼模型的阻尼系數進行修正，得到改進遲滯阻尼模型，并提出了對應的混合結構模態(tài)疊加法。同時，分析改進遲滯阻尼模型的反應譜特點，計算最不利地震作用效應組合，并提出了基于改進遲滯阻尼模型的混合結構抗震分析SRSS法。

1 基于改進遲滯阻尼模型的SRSS法

1.1 改進遲滯阻尼模型

遲滯阻尼模型假定阻尼力與結構的位移大小成正比，與結構的速度方向相反[16-17]。單自由度體系的遲滯阻尼運動方程為

由式（15）可計算出對應的自振頻率和特征向量，采用改進遲滯阻尼模型時，Kf可視為結構的動態(tài)剛度矩陣，結構的振動過程將是剛度不斷變化的過程。由式（15）可知，方程的特征值和特征向量不僅與結構的質量、剛度和阻尼有關，還與運動過程中的符號矩陣有關。對于n自由度結構，符號矩陣的個數為2n，因此，需要計算對應的2n個特征值和特征向量。當符號矩陣確定時，利用每個符號矩陣對應的特征向量，可對式（11）進行解耦計算，進而實現改進遲滯阻尼模型的模態(tài)疊加法。

利用特征向量和加速度反應譜，可計算每個符號矩陣對應下的地震作用效應，但為了保證結構設計的安全性，用于設計的地震作用效應須選擇最不利情況，即多種地震作用效應中的最大值，從而可實現改進遲滯阻尼模型的抗震分析SRRS法。

2 基于改進遲滯阻尼模型的加速度反應譜

采用改進遲滯阻尼模型時，單自由度體系的時域積分計算仍按照時間步長Δt進行離散，任意時刻可表示為tk=kΔt（k=0，1，2…）。當Δt取值足夠小時，可利用tk時刻結構的動力響應確定tk+1時刻的符號矩陣，并得到tk+1時刻的運動方程，由此進行時域迭代計算。利用單自由度的加速度響應，可得到基于改進遲滯阻尼模型的加速度反應譜。

加速度反應譜是結構抗震分析SRSS法的重要依據。以兩條地震波為例，分析了改進遲滯阻尼模型加速度反應譜和粘性阻尼模型反應譜之間的差異。

在圖2和圖3中，計算所得的兩條地震波作用下改進遲滯阻尼模型（CZ）的加速度譜值和粘性阻尼模型（NZ）的加速度譜值表明：在短周期段，CZ計算的加速度譜值小于NZ計算結果;在中周期段和長周期段，CZ計算的加速度譜值大于NZ計算結果。

由于本文的重點在于研究遲滯阻尼理論的改進以及混合結構SRSS法的實現方法，對于遲滯阻尼理論反應譜的深入分析以及兩種阻尼模型對應的設計反應譜的差異和對地震響應計算結果的影響，有待后續(xù)研究進一步澄清。與此同時，考慮到基于粘性阻尼理論的規(guī)范設計反應譜存在人為抬高拉平的做法[18]，故在后續(xù)算例中采用改進遲滯阻尼模型進行SRSS計算時，仍暫時借用規(guī)范的設計反應譜。

3 算例分析

如圖4所示，以不同材料組成的3層剪切型框架結構為例，其具體參數為：抗震設防烈度為8度，設計基本地震加速度為0.20g，設計地震分組為第2組，場地類別為II類，鋼結構的阻尼比為0.02，鋼筋混凝土結構的阻尼比為0.05。

對于3層框架，其符號矩陣共有8個，因此，需要計算每個符號矩陣對應運動方程下的框架層間剪力。如上文所述，基于改進遲滯阻尼模型的SRSS法可暫時采用規(guī)范反應譜，據此可計算出8個符號矩陣對應下的層間剪力，其結果如表1和圖5所示。

在圖5中，繪制了層間剪力的散點圖，可得到最大剪力包絡線，即選擇最不利地震作用效應組合，從而計算出改進遲滯阻尼模型下的框架層間剪力。

混合結構基于粘性阻尼模型的SRSS法可采用等效阻尼比ξd計算層間剪力。選擇等效阻尼比時，若直接取阻尼比的平均值，可得等效阻尼比為0.035;若采用復頻率法[6]確定混合結構振型阻尼比，將質量參與系數作為權系數進行加權平均，得到的混合結構等效阻尼比為0.038 9。分別采用粘性阻尼模型的SRSS法（NZ_SRSS）和改進遲滯阻尼模型的SRSS法（CZ_SRSS）計算框架的層間剪力，其結果見表2。

CZ_SRSS的計算結果較NZ_SRSS（ξd=0.035）計算的層間剪力增加幅度為14.81%～16.75%，表明CZ_SRSS的計算結果普遍較大，且頂層剪力增幅最大，達到16.75%;CZ_SRSS的計算結果較NZ_SRSS（ξd=0.038 9）計算的層間剪力增加幅度為16.92%～19.14%，其增幅規(guī)律類似，且頂層剪力增幅達到19.14%。

對于N層剪切型框架，第k層結構在符號矩陣中對應的對角元素為1，其他層結構對應的對角元素為-1，即

τ 此時，相比粘性阻尼模型下原結構的剛度分布情況，第k層結構對應的剛度增大，其他層對應的剛度減小，進而導致計算結果中第k層結構對應的地震作用效應增大。

基于改進遲滯阻尼模型的SRSS法選擇每個符號矩陣對應下的最大地震作用效應，進而形成最不利地震作用效應組合，使得改進遲滯阻尼模型計算的地震作用效應大于粘性阻尼模型。

4 結論

經理論推導和算例分析，得以下結論：

1）改進遲滯阻尼模型具有每周期耗散能量與外激勵頻率無關的優(yōu)點，但改進遲滯阻尼系統不再保持線性特征，具有局部的非線性特征。

2）與粘性阻尼模型相比，基于改進遲滯阻尼模型的模態(tài)疊加法可適用于非比例阻尼矩陣的混合結構。

3）基于改進遲滯阻尼模型的混合結構抗震分析SRSS法需要計算每種符號矩陣下的地震作用效應，計算量相對較大;若選擇最不利地震作用效應組合，改進遲滯阻尼模型計算的地震作用效應比粘性阻尼模型計算的地震作用效應更大，其增大幅度約為15%～20%。

參考文獻：

[1] 梁超鋒， 劉鐵軍， 鄒篤建，等. 材料黏滯系數與損耗因子的頻率相關性研究[J]. 力學學報， 2012， 44（5）： 933-937.

LIANG C F， LIU T J， ZOU D J， et al. The frequency-dependent study on viscosity coefficient and loss tangent of viscoelastic materials [J]. Chinese Journal of Theoretical and Applied Mechanics， 2012， 44（5）： 933-937. （in Chinese）

[2] 潘旦光，高莉莉. Rayleigh阻尼系數解法比較及對結構地震反應影響[J]. 工程力學，2015，32（6）： 192-199.

PAN D G， GAO L L. Comparison of determination methods for Rayleigh damping coefficients and effects on seismic responses of structures [J]. Engineering Mechanics， 2015， 32（6）： 192-199. （in Chinese）

[3] 黃維，錢江，周知. 基于Rayleigh阻尼模型的豎向混合結構設計阻尼比研究[J].工程力學，2015， 32（10）： 60-67.

HUANG W， QIAN J， ZHOU Z. Research on equivalent damping ratio of vertically mixed structures based on Rayleigh damping model [J]. Engineering Mechanics， 2015， 32（10）： 60-67.（in Chinese）

[4] HUANG W， QIAN J， ZHOU Z， et al. An approach to equivalent damping ratio of vertically mixed structures based on response error minimization [J]. Soil Dynamics and Earthquake Engineering， 2015， 72： 119-128.

[5] 劉慶林，傅學怡. 基于復阻尼假定的不同材料阻尼特性混合結構抗震分析反應譜CCQC法[J]. 土木工程學報，2011， 44（3）： 61-71.

LIU Q L， FU X Y. A response spectrum CCQC method for seismic analysis of structures of multiple material damping characteristics based on complex damping assumption [J]. China Civil Engineering Journal， 2011， 44（3）： 61-71. （in Chinese）

[6] 黃本才. 組合結構振動的等效阻尼比[J]. 上海力學，1998， 19（2）：141-145.

HUANG B C. Equivalent modal damping ratios of vibration in composite structures [J]. Shanghai Journal of Mechanics， 1998， 19（2）：141-145. （in Chinese）

[7] 薛彥濤，韋承基，孫仁范，等. 采用不同材料加層時結構阻尼比計算方法（應變能法）[J]. 工程抗震與加固改造， 2008， 30（2）： 91-95.

XUE Y T， WEI C J， SUN R F， et al. Calculation method for damping ration of different story added structures （strain energy method） [J]. Earthquake Resistant Engineering and Retrofitting， 2008， 30（2）： 91-95. （in Chinese）

[8] 梁仁杰，吳京，孟少平. 能力譜法中等效阻尼比的誤差分析及假設檢驗[J]. 土木工程學報，2010， 43（Sup1）： 33-36.

LIANG R J， WU J， MENG S P. Error analysis and assumption test on the equivalent damping ratio in capacity spectrum method [J]. China Civil Engineering Journal， 2010， 43（Sup1）： 33-36. （in Chinese）

[9] 汪夢甫. 非比例阻尼線性體系地震反應計算的振型分解反應譜法[J]. 地震工程與工程振動，2007， 27（1）： 31-37.

WANG M F. On seismic response analysis methods of non-proportional damped linear MDOF systems [J]. Journal of Earthquake Engineering and Engineering Vibration， 2007， 27（1）： 31-37. （in Chinese）

[10] NEUGEBAUER R， SCHEFFLER C， WABNER M， et al. Advanced state space modeling of non-proportional damped machine tool mechanics [J]. CIRP Journal of Manufacturing Science and Technology， 2010， 3（1）： 8-13.

[11] NEUGEBAUER R， SCHEFFLER C， WABNER M， et al. State space modeling of non-proportional passive damping in machine tools [J]. The International Journal of Advanced Manufacturing Technology， 2011， 53： 945-952.

[12] 劉慶林，傅學怡，孫占琦. 基于復阻尼假定的不同材料阻尼特性混合結構抗震分析復模態(tài)疊加法[J]. 建筑結構學報，2011， 32（9）： 27-33.

LIU Q L， FU X Y， SUN Z Q. A complex mode superposition method for seismic analysis of structures of multiple material damping characteristics based on complex damping assumption [J]. Journal of Building Structures， 2011， 32（9）：27-33. （in Chinese）

[13] 何鐘怡，廖振鵬，王小華. 關于復阻尼理論的幾點注記[J]. 地震工程與工程振動，2002，22（1）： 1-6.

HE Z Y， LIAO Z P， WANG X H. Some notes on theory of complex damping [J]. Earthquake Engineering and Engineering Vibration， 2002， 22（1）： 1-6. （in Chinese）

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

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

[16] 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.

[17] 朱鏡清. 頻率相關粘性阻尼理論及有關問題的解[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）

[18] 朱鏡清，朱敏. 復阻尼地震反應譜的計算方法及其它[J]. 地震工程與工程振動，2000， 20（2）： 19-23.

ZHU J Q， ZHU M. Calculation of complex damping response spectra from earthquake records [J]. Earthquake Engineering and Engineering Vibration， 2000， 20（2）： 19-23. （in Chinese）

（編輯 王秀玲）