基于自激勵人體模型的人-結(jié)構(gòu)豎向相互作用研究

謝偉平， 馮金鵬， 何 衛(wèi)

(1.武漢理工大學 土木工程與建筑學院，武漢 430070; 2.武漢理工大學 理學院，武漢 430070)

基于自激勵人體模型的人-結(jié)構(gòu)豎向相互作用研究

謝偉平1， 馮金鵬1， 何 衛(wèi)2

(1.武漢理工大學 土木工程與建筑學院，武漢 430070; 2.武漢理工大學 理學院，武漢 430070)

基于自激勵人體模型，建立了考慮人-結(jié)構(gòu)豎向相互作用下行人過橋運動方程，提出了考慮步頻主諧和亞諧分量的單人步行力模型，并通過單人過橋試驗對該模型進行了驗證；分析了人-橋豎向相互作用的變化規(guī)律及影響因素。結(jié)果表明，人體步頻亞諧可能與結(jié)構(gòu)的基頻相近，從而引起較大的結(jié)構(gòu)響應(yīng)；結(jié)構(gòu)固有頻率和人/結(jié)構(gòu)質(zhì)量比對人-結(jié)構(gòu)相互作用影響較大，當結(jié)構(gòu)頻率與步頻及其主諧波分量頻率接近時，耦合效應(yīng)最顯著，其對人體足底力影響最大。隨著人/結(jié)構(gòu)質(zhì)量比的增加，人-結(jié)構(gòu)相互作用逐漸增強。

人致振動；豎向人體動力學模型；人-結(jié)構(gòu)相互作用；足底力

隨著經(jīng)濟的快速發(fā)展及高性能建筑材料的出現(xiàn)，造型獨特、結(jié)構(gòu)新穎的大跨度公用建筑大量涌現(xiàn)，如大跨人行天橋、大型車站、體育場館等。這些大跨度結(jié)構(gòu)大都具有自振頻率低，人群密度大等特點，所以其人致振動響應(yīng)問題引起了工程界和學術(shù)界的極大關(guān)注。這類結(jié)構(gòu)的人致振動問題并不是簡單意義上的強迫振動，而是由強迫振動引發(fā)的、包含人-結(jié)構(gòu)動力相互作用的一種復(fù)雜的結(jié)構(gòu)振動形式。然而，由于缺乏對人-結(jié)構(gòu)相互作用效應(yīng)深層次的認識，近年來學者們在分析此類結(jié)構(gòu)的人致振動響應(yīng)時多采用基于剛性測力平臺實測的足底反力提出的傅里葉級數(shù)模型[1-3]進行求解，而該模型并未考慮人與結(jié)構(gòu)間的相互作用對人體行走時足底反力的影響[4]，從而導(dǎo)致計算結(jié)果與實測值相偏離。因此，研究人-結(jié)構(gòu)相互作用對分析結(jié)構(gòu)的人致振動問題顯得尤為重要，這一問題的解決將有助于提高人致振動計算理論的準確性和可靠性。

有研究表明當行人在剛度較小、基頻較低的結(jié)構(gòu)上行走時，產(chǎn)生的足底反力與在剛性測力平臺上產(chǎn)生的足底反力相比有所變化[5]。Ohlsson[6]測得大跨輕型樓板上的豎向步行力不同于在剛性地面；Pimentel[7]研究發(fā)現(xiàn)實際人行橋測得的人行荷載動載因子不同于剛性地面的測試結(jié)果；Yao等[8]在研究人體跑步、跳躍荷載時發(fā)現(xiàn)在柔性結(jié)構(gòu)上測得的人行荷載小于在剛性地面的數(shù)據(jù)。在國內(nèi)，秦敬偉等[9]基于雙足模型對考慮人-橋豎向耦合作用下的過橋過程進行了模擬，發(fā)現(xiàn)人體于柔性結(jié)構(gòu)上行走時，產(chǎn)生的足底反力小于剛性地面時的結(jié)果。謝偉平等[10]就人-結(jié)構(gòu)相互作用對人體行走時足底反力的影響做了初步的探索。綜上所述，目前學界對于人-結(jié)構(gòu)相互作用會影響人體行走時的足底反力這一結(jié)論已形成共識，但對于考慮耦合作用下足底反力的影響因素及變化規(guī)律的研究還有待開展。

本文采用單自由度人體模型，將人體簡化為帶有生物自激勵的彈簧-質(zhì)量-阻尼系統(tǒng)，建立考慮人-結(jié)構(gòu)豎向耦合作用下行人過橋時的運動方程，提出了考慮步頻主諧和亞諧分量的單人步行力模型，并通過在人行鋼橋上進行試驗實測，與模擬結(jié)果進行對比，從而驗證模型的合理性與有效性。同時還進一步分析人/結(jié)構(gòu)質(zhì)量比和結(jié)構(gòu)固有頻率等因素對行人足底反力的影響，揭示人-橋相互作用的影響因素及變化規(guī)律。

1 人-結(jié)構(gòu)豎向相互作用模型

1.1運動方程的建立

假設(shè)結(jié)構(gòu)為等截面，即抗彎剛度EI為常數(shù)，線密度為mb，阻尼系數(shù)為c，梁上的外加激勵是由人體行走而引起的，如圖1所示。

圖1 人-結(jié)構(gòu)相互作用模型

行人激勵下，梁的運動方程表達式為[11]

(1)

式中：w為梁的撓度;Fp(x,t)為人對梁的作用力，表達式為

Fp(x,t)=fp(t)·δ(x-xp(t))

(2)

fp(t)為足底作用力時程，利用人體質(zhì)心加速度表示如下

fp(t)=-mpg-mpac(t)

(3)

式中：mp為人體質(zhì)量;ac(t)為人體質(zhì)心加速度。梁振動時，人體的質(zhì)心加速度分為兩個部分：相對于結(jié)構(gòu)運動產(chǎn)生的相對加速度acr(t)和隨結(jié)構(gòu)共同運動產(chǎn)生的牽連加速度ace(t)，即

ac(t)=acr(t)+ace(t)

(4)

ace(t)指人所處位置梁的振動加速度，可表示為

(5)

不考慮結(jié)構(gòu)振動對人體行走步態(tài)的影響，則有

(6)

(7)

可見式(7)實質(zhì)上將行人過橋時的作用力等效為剛性地面足底力與橋面振動引起的慣性力之和。利用振型疊加法求解式(7)，簡支歐拉-伯努利梁振型函數(shù)為

(8)

則有

φi(x)

(9)

將式(9)代入式(7)則有

(10)

將式(10)中的每一項乘以第j個振型函數(shù)φj(x)，并沿梁的全長進行積分，右端的結(jié)果為

(11)

則廣義力為

(12)

將式(12)代入式(10)整理可得

(13)

則有

(14)

這是一個無窮多自由度的聯(lián)立方程組。如果位移級數(shù)取N項，則結(jié)構(gòu)的廣義自由度為N，系統(tǒng)的N階矩陣表達式為

(15)

式中,廣義位移向量U=[q1,q2,…,qN]T，廣義的質(zhì)量矩陣、阻尼矩陣、剛度矩陣以及力向量分別如式(16)～(19)所示。

(16)

(17)

(18)

F=[ρFφ1,ρFφ2,…ρFφN]T

(19)

人體在結(jié)構(gòu)上移動時，方程中的系數(shù)會不斷變化，因此式(15)為時變系統(tǒng)的多階微分方程組，應(yīng)采用非線性的動力分析方法。本文基于Newmark-β法對上述方程進行求解。

1.2單人步行力模型的建立

(20)

2 荷載模型的驗證

本文通過在一輕柔單跨人行鋼橋上進行單人行走實驗來驗證模型的有效性。該人行橋全長10.5 m，有效長度10.3 m，兩邊簡支，線密度為73.7 kg/m，經(jīng)實測，該橋一階豎彎頻率為2.734 Hz。如圖2所示。

圖2 人行鋼橋

開展單人過橋試驗時令體重為85 kg的行人在節(jié)拍器的提示下以不同步頻通過人行橋，圖3為步頻為1.83 Hz時實測跨中加速度時程曲線和文中模擬結(jié)果對比，圖4為兩者的頻譜對比。

(a) 實測結(jié)果

(b) 本文模擬結(jié)果

由圖3可以看出，模擬結(jié)果與實測響應(yīng)波形走勢基本吻合。從頻域上來說，兩者的頻率成分基本吻合，都可以明顯的看出結(jié)構(gòu)的一階頻率和人體的前三階步頻。實測頻率成分相對更為豐富，包含了結(jié)構(gòu)的高階模態(tài)和人體的高階步頻成分。實測頻譜在結(jié)構(gòu)一階頻率處達到最大峰值，原因是人體行走步頻的二階亞諧(1.5fs)與該階頻率相近，從而導(dǎo)致此處響應(yīng)較大。模擬頻譜中步頻二階亞諧處也出現(xiàn)了一個峰值，原因同上。此外，實測結(jié)果的幅值略小于模擬結(jié)果，這是由于實測的頻率成分更加豐富，能量分布更加分散。

(a) 實測結(jié)果

(b) 本文模擬結(jié)果

表1則給出了不同步頻下跨中均方根加速度、峰值加速度的實測結(jié)果與模擬結(jié)果對比，從表中可以看出，模擬結(jié)果最大誤差不超過10%，基本滿足精度要求。

表1 不同步頻下模擬與實測結(jié)果對比

通過上述分析可知，文中的人-橋豎向耦合系統(tǒng)的動力分析和求解程序基本正確，可用于研究人-橋豎向耦合系統(tǒng)的人致振動相關(guān)問題。

3 足底反力影響因素分析

目前在對人行橋進行振動分析時，選用的人體步行荷載多為傅里葉級數(shù)模型荷載或單步落足荷載，然而這些模型幾乎都是基于固定地面或剛性結(jié)構(gòu)上的測試結(jié)果提出來的，與實際存在差異。近年來隨著柔性結(jié)構(gòu)的逐漸增多，人們越來越認識到人與結(jié)構(gòu)之間存在相互作用影響。有學者研究發(fā)現(xiàn)當人體在柔性結(jié)構(gòu)上行走時的足底力與剛性地面所測得足底力有所不同，并且兩者的差異與相互作用的大小有關(guān)。

為分析人-結(jié)構(gòu)豎向相互作用對足底力的影響，將人體于柔性結(jié)構(gòu)上行走時的足底力與剛性測力平臺上所測足底力的差值進行無量綱化處理，定義相互作用力影響函數(shù)[13]

(21)

由式(10)可知，式(21)可表示為

(22)

即DI(t)為行人所在位置處結(jié)構(gòu)加速度和重力加速度的比值。

3.1人/結(jié)構(gòu)質(zhì)量比對足底反力的影響

算例梁長9 m，寬0.8 m，高0.12 m，阻尼比0.003，基頻為2 Hz，為研究人/結(jié)構(gòu)質(zhì)量比對足底反力的影響，令人/結(jié)構(gòu)質(zhì)量比為u，保持結(jié)構(gòu)基頻不變，令人體步頻與結(jié)構(gòu)基頻相同，并通過對梁體的材料參數(shù)(密度和彈性模量)進行改變以實現(xiàn)不同的人/結(jié)構(gòu)質(zhì)量比。

圖5所示為不同人/結(jié)構(gòu)質(zhì)量比下無量綱相互作用力影響函數(shù)曲線。

相互作用力影響函數(shù)DI曲線反映了人運動位置處結(jié)構(gòu)的加速度狀況，由圖5可以看出，當行人由梁端部向跨中行走時，DI逐漸增大，隨著行人逐漸遠離，DI逐漸減小直到行人離開結(jié)構(gòu)。DI達到最大值的位置并不在結(jié)構(gòu)的跨中，而是在其附近的某一位置，具體位置與人/結(jié)構(gòu)質(zhì)量比有關(guān)，原因是人體行走時人-橋耦合系統(tǒng)為一時變體系，結(jié)構(gòu)的振型跡線在不斷發(fā)生變化；當人/結(jié)構(gòu)質(zhì)量比較小時，人/結(jié)構(gòu)相互作用隨著質(zhì)量比的增加逐漸增強，相互作用力影響系數(shù)逐漸增大，但是當質(zhì)量比增大到一定程度時，人-結(jié)構(gòu)相互作用力反而減小，原因是較大的質(zhì)量比使此時結(jié)構(gòu)的自振頻率發(fā)生了改變，導(dǎo)致激振頻率遠離了結(jié)構(gòu)的自振頻率。

(a) u=0.01

(b) u=0.05

(c) u=0.1

Fig.5 The interaction force time-history curve of different human/structure mass ratio

3.2結(jié)構(gòu)頻率對足底反力的影響

算例保持3.1節(jié)算例中梁的截面參數(shù)及阻尼特性不變，線密度為230.4 kg/m，梁1楊氏模量E=30 GPa，基頻f=2.375 Hz；梁2楊氏模量E=14.5 GPa，基頻f=1.65 Hz。人體質(zhì)量80 kg，步頻取2 Hz。

圖6所示為人體分別在梁1和梁2上行走時的歸一化足底反力(足底力/人體體重)時程圖，通過圖中可以看出，人體在柔性結(jié)構(gòu)上行走時的足底反力與剛性地面的有所不同，隨著人體向跨中移動，足底反力增大或減小，具體變化規(guī)律與結(jié)構(gòu)頻率有關(guān)。考慮到人/結(jié)構(gòu)質(zhì)量比很小，為3.7%，因此可以認為：當結(jié)構(gòu)剛度較小，基頻較低時，人與結(jié)構(gòu)件的相互作用較為顯著，結(jié)構(gòu)振動對人行荷載產(chǎn)生的影響也較為明顯。

為進一步研究結(jié)構(gòu)頻率對足底反力的影響，定義DI的幅值為相互作用力影響系數(shù)，記為DF

DF=DI(t)max

(a) 梁1(f=2.375 Hz)

(b) 梁2(f=1.65 Hz)

(23)

圖7為不同人/結(jié)構(gòu)質(zhì)量比時相互作用力影響系數(shù)DF與結(jié)構(gòu)頻率的關(guān)系曲線。

圖7 相互作用力影響系數(shù)與結(jié)構(gòu)頻率的關(guān)系

Fig.7 The relationship between the influence coefficient of the interaction force and the frequency of the structure

由圖7可知，當結(jié)構(gòu)頻率與激振頻率(2 Hz)及其某階主諧波頻率相近時，人-結(jié)構(gòu)耦合效應(yīng)十分顯著，此時行人足底反力與剛性地面相比有顯著變化，而隨著結(jié)構(gòu)頻率的增大，相互作用力影響系數(shù)下降較快且趨于平緩。當結(jié)構(gòu)頻率大于8 Hz時，人-結(jié)構(gòu)相互作用已經(jīng)很弱，表明行人足底力與在剛性地面行走時差別很小，此時在人/結(jié)構(gòu)質(zhì)量比較小的情況下可以不考慮人-結(jié)構(gòu)相互作用的影響，而將行人動荷載作為外激勵；隨著人/結(jié)構(gòu)質(zhì)量比的增加，相互作用力影響系數(shù)逐漸增大，這主要是由于人體參振質(zhì)量增大從而使人-結(jié)構(gòu)相互作用逐漸增強。因此對于大跨度輕質(zhì)、柔性結(jié)構(gòu)(如大跨人行橋等)，當結(jié)構(gòu)上人員密度較大時，需考慮人-結(jié)構(gòu)相互作用的影響。

4 結(jié) 論

本文基于自激勵彈簧-質(zhì)量-阻尼模型，建立了人-結(jié)構(gòu)相互作用系統(tǒng)的運動方程，與實測對比驗證了模型有效性后，進行了數(shù)值分析，結(jié)果表明：

(1)人體于結(jié)構(gòu)上行走時，步頻亞諧(0.5fs，1.5fs，2.5fs)可能與結(jié)構(gòu)的基頻相近，從而引起較大的結(jié)構(gòu)響應(yīng)。

(2)結(jié)構(gòu)固有頻率和人/結(jié)構(gòu)質(zhì)量比對人-結(jié)構(gòu)相互作用有較大影響。當結(jié)構(gòu)頻率與步頻及其主諧波分量頻率接近時，耦合效應(yīng)最顯著，其對人體足底反力影響最大。結(jié)構(gòu)頻率越大，人-結(jié)構(gòu)相互作用就越弱，此時對于人/結(jié)構(gòu)質(zhì)量比較小的情況，可將結(jié)構(gòu)上的人行荷載視為外激勵。隨著人/結(jié)構(gòu)質(zhì)量比的增加，相互作用力系數(shù)逐漸增大，人-結(jié)構(gòu)相互作用逐漸增強。因此對于大跨度輕質(zhì)柔性結(jié)構(gòu)在人員密度較大時，需考慮人-結(jié)構(gòu)相互作用的影響。

[1] 陳雋, 王浩祺, 彭怡欣. 行走激勵的傅里葉級數(shù)模型及其參數(shù)的實驗研究[J]. 振動與沖擊, 2014, 33(8):11-15.

CHEN Jun, WANG Haoqi, PENG Yixin. Experimental investigation on fourier-series model of walking load and its coefficients [J]. Journal of Vibration and Shock, 2014, 33(8): 11-15.

[2] 李愛群, 陳鑫, 張志強. 大跨樓蓋結(jié)構(gòu)減振設(shè)計與分析[J]. 建筑結(jié)構(gòu)學報, 2010, 31(6):160-170.

LI Aiqun, CHEN Xin, ZHANG Zhiqiang. Design and analysis on vibration control of long-span floor structures [J]. Journal of Building Structures, 2010, 31(6):160-170.

[3] 聶建國, 陳宇, 樊健生. 步行荷載作用下單跨人行橋振動的均方根加速度反應(yīng)譜法[J]. 土木工程學報, 2010, 43(9): 109-116.

NIE Jianguo, CHEN Yu, FAN Jiansheng. RMS acceleration response spectrum method for single-span footbridges under pedestrian load [J]. China Civil Engineering Journal, 2010, 43(9):109-116.

[4] 袁旭斌. 人行橋人致振動特性研究[D]. 上海: 同濟大學, 2006.

[5] 李紅利, 陳政清. 人-橋豎向動力相互作用響應(yīng)理論與試驗研究[J]. 土木工程學報, 2014, 47(6):78-87.

LI Hongli, CHEN Zhengqing. Analytical and experimental study on vertically dynamic interaction between human and bridge [J]. China Civil Engineering Journal, 2014, 47(6):78-87.

[6] OHISSON S. Floor vibration and human discomfort [D]. Gortborg: Chalmers University of Technology, 1982.

[7] PIMENTEL R L. Vibrational performance of pedestrian bridges due to human-induced loads [D]. Sheffield: University of Sheffield, 1997.

[8] YAO S, WRIGHT J R, PAVIC A, et al. Forces generated when bouncing or jumping on a flexible structure [C].Proceedings of the International Conference on Noise and Vibration, 2002, 2: 563-572.

[9] 秦敬偉, 楊慶山. 基于雙足步行模型和反饋機制的人體-結(jié)構(gòu)相互作用[J]. 建筑結(jié)構(gòu)學報, 2014, 35(增刊1):18-24.

QIN Jingwei, YANG Qingshan. Human-structure interaction based on bipedal walking model and feedback mechanism [J]. Journal of Building Structures, 2014, 35(Sup1):18-24.

[10] 謝偉平, 魯偉, 何衛(wèi). 人-板耦合系統(tǒng)動力響及相互作用研究[J]. 武漢理工大學學報, 2013, 37(1):56-62.

XIE Weiping, LU Wei, HE Wei. Research on dynamic response and interaction of human-floor coupled system [J]. Journal of Wuhan University of Technology, 2013, 37(1):56-62.

[11] 劉隆. 人群與梁的豎向相互作用及人行橋側(cè)向振動機理研究[D]. 武漢: 武漢理工大學, 2013.

[13] 何衛(wèi). 大跨度高鐵車站振動舒適度理論與試驗研究[D]. 武漢: 武漢理工大學, 2014.

Verticalhuman-structureinteractionbasedonaself-excitinghumanmodel

XIEWeiping1,FENGJinpeng1,HEWei2

(1. School of Civil Engineering and Architecture, Wuhan University of Technology, Wuhan 430070, China; 2. School of Science, Wuhan University of Technology, Wuhan 430070, China)

Based on a self-exciting human model, the motion equation of a footbridge was established in considering the vertical human-structure interaction and a walk force model was proposed in considering the main harmonic and subharmonic frequencies of walking. The model was verified by measuring the footbridge vibration responses caused by single person walking. The variation of vertical human-structure interaction and its influencing factors were analyzed. The results indicate that the walking subharmonic frequency may be close to the fundamental frequency of the structure, resulting in large structural responses. The natural frequency of the structure and the human/structure mass ratio have a great influence on the human-structure interaction. When the structural frequency is closer to the walking frequency or its main harmonic components, the coupling effect is significant and has the greatest impact on footfall forces. Moreover, with the increasing of human/structure mass ratio, the interaction will strengthen gradually.

human-induced vibration; vertical human dynamic model; human-structure interaction; footfall force

TU311.3

A

10.13465/j.cnki.jvs.2017.21.005

國家自然科學基金(51508431);中國博士后科學基金(2015M582288)

2016-02-02 修改稿收到日期：2016-06-15

謝偉平 男，博士，教授，1965年6月