陳 雋， 周戈杰， 王鵬程

(1.同濟大學 土木工程學院，上海 200092；2.同濟大學 土木工程防災國家重點實驗室，上海 200092)

單人、人群完成步行、跳躍、跑步、踏步、原地屈伸、跳舞、瞬時起立或坐下等動作時通過雙足對支撐結(jié)構(gòu)施加的動力作用稱為人致荷載，由此導致的工程結(jié)構(gòu)振動稱為人致結(jié)構(gòu)振動。振動幅度過大時會造成結(jié)構(gòu)上的使用者對振動有感、頭暈、心慌等振動舒適度問題，或者影響對振動敏感的高精密設備的正常使用，嚴重時則可能造成結(jié)構(gòu)局部構(gòu)件破壞甚至整體破壞、人員傷亡等安全問題[1-6]。

跑步運動是極為常見的一種運動形式，尤其隨著近年來的全民健身熱，參加5 km、馬拉松之類的體育活動的人數(shù)迅速增加。據(jù)報道，2014年我國馬拉松賽事運營的收入就達到20億元，市場規(guī)模超百億元，2015年又增了近50%，2019年更是達到了1 828場次，如表1所示。此外，一些多層體育場館室內(nèi)、高層建筑的內(nèi)部以及建筑屋頂?shù)纫查_始設計有跑道，方便建筑使用者就近開展跑步活動，如圖1所示。

在跑步運動中，當人的跑步頻率與支撐結(jié)構(gòu)的自振頻率接近時，支撐結(jié)構(gòu)會產(chǎn)生較大的振動，可能造成振動舒適度問題，嚴重時甚至可能誘發(fā)公共安全事件。尤其伴隨著建筑材料輕質(zhì)化、建筑功能多樣化的需求，以及開敞、輕盈建筑美學的普及，建筑樓蓋、人行橋等結(jié)構(gòu)設計向著大跨度、輕柔方向發(fā)展。此類結(jié)構(gòu)自振頻率與人的跑步頻率范圍重合，容易出現(xiàn)人致結(jié)構(gòu)振動問題。例如：美國舊金山金門大橋曾由于大量人群跑步過橋造成大幅振動，引起橋上行人的不適；成都某理工大學校園內(nèi)一跨度不大的混凝土人行橋，在學生早操跑步通過時劇烈振動，造成了人群的恐慌[7]。上述例子表明，在大跨度結(jié)構(gòu)設計過程中有必要進行人致結(jié)構(gòu)振動分析。

表1 近年來我國馬拉松賽事舉辦場次與參賽人次

(a)(b)

正確的荷載取值是分析跑步荷載所引起的結(jié)構(gòu)振動的基本前提，然而，我國目前關(guān)于跑步荷載模型的研究工作甚少，沒有可用的跑步荷載模型，相關(guān)的基礎(chǔ)性試驗研究工作也非常缺乏。國外有學者采用類似跳躍荷載的“半正弦波”模型建立了跑步荷載模型[8]，但跑步荷載在頻率、幅值、空間移動特性上與跳躍荷載有明顯的區(qū)別，不宜直接使用。工程設計中常常用“有節(jié)奏荷載”的模型代表跑步荷載，然而“有節(jié)奏運動荷載”一般指健美操、舞蹈、節(jié)奏性跳躍或Bounce等，其對應的試驗數(shù)據(jù)基礎(chǔ)也與跑步荷載不同。

因此，本文通過試驗收集了一批單人跑步荷載時程，建立了跑步荷載的傅里葉級數(shù)模型，給出了模型參數(shù)取值方法。進而開展了實際連廊結(jié)構(gòu)跑步作用下的振動實測，對荷載模型的應用效果進行了檢驗。

1 跑步荷載試驗

大量跑步荷載實測時程是建立荷載模型的前提。本研究采用Novel公司的Pedar高精度無線測力鞋墊技術(shù)開展了跑步荷載試驗,如圖2所示。測力鞋墊用于測量人在跑步時的足底壓力，通過鞋墊內(nèi)的分布式壓力傳感器實現(xiàn)，可以連續(xù)實時記錄所受壓力，并通過積分獲得整體壓力。測力鞋墊技術(shù)是一種新型試驗手段，解決了傳統(tǒng)固定測力板方式無法獲得連續(xù)跑步荷載、也無法實現(xiàn)長距離測試的問題。

試驗共有25名健康成年測試者參與，均為同濟大學學生。試驗前所有測試者均按照試驗倫理要求簽署了知情同意書。每位測試者共完成7組跑步工況的試驗，跑步頻率范圍為2.4～3.6 Hz，間隔為0.2 Hz，對應的跑步速度為2.8～6.6 m/s。測力鞋墊的采樣頻率為100 Hz，每個工況重復兩次。對于每個工況，試驗志愿者都需在節(jié)拍器的引導下勻速連續(xù)跑步20 s以上，結(jié)束后需休息充分至無疲憊感才可進行下一工況的試驗。

試驗總共收集了345條有效連續(xù)跑步荷載時程，各工況有效荷載數(shù)據(jù)數(shù)量的統(tǒng)計信息如表2所示，各個工況數(shù)據(jù)量的均衡為荷載模型的準確性提供了保障。

(a) 現(xiàn)場情況

表2 各工況有效荷載數(shù)據(jù)數(shù)量統(tǒng)計

圖3為某測試者在低頻(2.4 Hz)、中頻(3.0 Hz)、高頻(3.6 Hz)工況時的跑步荷載名義力(豎向雙足合力/體重)時程以及對應的傅里葉幅值譜。由圖3可見，跑步荷載時程具有顯著的周期性，中頻率時頻譜上前4～5階峰值明顯，分別對應于設定的跑步頻率及其倍頻。同時可以發(fā)現(xiàn)，在低頻和高頻工況，荷載的頻譜有一定的“雜化”現(xiàn)象，各階波峰呈現(xiàn)一定的寬度和“毛刺”，其對應的試驗現(xiàn)象和背景是：受人體生理機能的限制，較低或較高頻率，跑步者很難長時間準確按照節(jié)拍器的節(jié)奏跑動。

2 跑步荷載模型

2.1 傅里葉級數(shù)模型

基于第1章的觀察，可將人的連續(xù)跑步過程視作為周期性過程，并采用傅里葉級數(shù)模型，即

(1)

式中：F(t)為豎向跑步荷載時程；G為人的靜止重力；N為模型截斷階數(shù)；t為時間；rn為第n階傅里葉系數(shù)，也稱為動載因子(dynamic loading factor，DLF)；fp為跑步頻率；φn為第n階相位角。

可見，跑步荷載傅里葉級數(shù)模型是正弦函數(shù)的疊加，其用于結(jié)構(gòu)動力分析是十分方便的。此外，現(xiàn)有很多設計規(guī)范中所推薦的人致荷載模型均為傅里葉級數(shù)模型。綜上，本文采用傅里葉級數(shù)建立跑步荷載的模型。

2.2 模型階數(shù)

合理的模型截斷階數(shù)N，即需要保證包含荷載的絕大部分能量，也應避免因階數(shù)過高造成的模型參數(shù)取值困難和使用不便。基于上述原則，經(jīng)綜合分析全部的實測數(shù)據(jù)，建議取N=5，即考慮荷載的前5階的貢獻。這一取值與已有跳躍荷載(N=3)[9]和步行荷載(N=5)[10]，具有可比性。

(a) 跑步頻率2.4 Hz

2.3 動載因子計算原則

建立荷載模型的關(guān)鍵在于動載因子的取值，計算動載因子最常用的方法是峰值法和等效能量法(簡稱能量法)。峰值法是取跑步荷載傅里葉幅值譜上各階諧波峰值處的幅值作為該階動載因子FDL,n。由于測試者的連續(xù)跑步過程并非嚴格意義上的周期性過程，其傅里葉幅值譜上峰值附近范圍內(nèi)存在能量泄漏現(xiàn)象。峰值法忽略了泄漏的能量，僅僅考慮各階諧波處的最大峰值，低估了跑步荷載值；而能量法將各階諧波峰值附近所有能量集中在諧波頻率作為動載因子，又高估計跑步荷載。因此，有學者[11]提出介于這兩種方法之間的窄帶分析方法，即將各階諧波附近能量限制在一定帶寬范圍內(nèi)計入模型計算，不考慮距離較遠的能量，相對合理。因此，本文采用限定帶寬的分析方法，并將各階諧波附近能量擴散的范圍取為(1.00±0.05)fp[12]。則第n階主頻附近計算范圍內(nèi)的能量為

(2)

式中：En為跑步荷載幅值譜的第n階諧波附近計算帶寬內(nèi)的能量；fp為跑步荷載的頻率；An(f)為跑步荷載傅里葉幅值譜的第n階諧波附近計算帶寬內(nèi)的傅里葉幅值。則第n階FDL為

(3)

2.4 動載因子計算結(jié)果

依據(jù)2.3節(jié)原則，首先將跑步荷載時程的長度截取為跑步周期的整數(shù)倍再求其傅里葉幅值譜，利用上述限制帶寬法計算所有實測跑步荷載時程的前5階動載因子如圖4所示。

圖4顯示DLF隨頻率的變化比較離散。從與步行、跳躍等現(xiàn)有荷載模型一致的角度，計算7個頻率工況下DLF的均值，并對其隨頻率的變化進行線性擬合，結(jié)果見圖4。表3為前5階動載因子在各個頻率下的變化范圍(括號內(nèi)為均值標準差)，所得線性回歸方程如表4所示。結(jié)果顯示第1、第4階DLF隨頻率增加而增加，第2、第3、第5階DLF隨頻率增加而減小，與陳雋等研究中跳躍荷載DLF的變化趨勢類似。

由于DLF的離散性較大，設計時可選用表4中的DLF均值回歸方程進行初步分析，或表3中的DLF變化范圍、均值和標準差進行精細化分析。

2.5 相位角取值

同樣的，可采用2.4節(jié)研究各階相位角隨跑步頻率的變化規(guī)律，圖5所示為前5階相位角隨跑步頻率的變化趨勢。數(shù)據(jù)顯示，各階相位角的試驗結(jié)果相當離散。由于人致結(jié)構(gòu)振動通常是單模態(tài)主導，DLF值對響應最大值起決定性作用，相位角只影響最大值出現(xiàn)(或相對于外激勵滯后)的時間。因此，本文采用目前常用的直接對各階相位角取均值的方式，得第1、第3、第5階相位角取π/12，第2、第4階相位角取0°。

3 結(jié)構(gòu)實測對比

3.1 連廊結(jié)構(gòu)

某連廊結(jié)構(gòu)如圖6(a)所示，在該結(jié)構(gòu)上進行單人跑步試驗，將實測結(jié)構(gòu)響應與采用第2章跑步荷載模型的計算響應進行比較，來檢驗荷載模型用于結(jié)構(gòu)設計的可行性。該連廊為鋼結(jié)構(gòu)，材料強度等級為Q345，彈性模量Es=2.06×105N/mm2，平面尺寸為24.0 m×4.5 m，高度300 mm，連廊兩端簡支。經(jīng)測試，結(jié)構(gòu)自振頻率為6.6 Hz，阻尼比為0.025。

(c)

(e)

表3 前5階動載因子在各個頻率下的變化范圍和均值、標準差

3.2 試驗過程

測試者穿戴好測力鞋墊后，按照節(jié)拍器引導頻率(如2 Hz，則為120 beat/min)在連廊上沿環(huán)形連續(xù)跑步30 s以上，以激發(fā)結(jié)構(gòu)的穩(wěn)態(tài)反應。同時，采用5個高靈敏度加速度傳感器(LC0132T)記錄連廊跨中等位置處的響應。測力鞋墊采樣頻率為100 Hz，加速度傳感器采樣頻率為200 Hz[13-14]。試驗現(xiàn)場和跑步路徑見圖6。

(a) 1階相位角

(d) 4階相位角

表4 各階動載因子回歸方程

(a) 試驗現(xiàn)場

(b) 跑步路徑及檢驗點

某測試者的體重為50 kg。按上述試驗要求完成了5種跑步頻率共12組的試驗。圖7為兩次跑步頻率為3.3 Hz時的跨中加速度響應時程，各工況結(jié)構(gòu)跨中加速度峰值見表5?？梢钥闯觯敎y試者的跑步頻率為3.3 Hz時，荷載時程的第2階諧波分量引發(fā)結(jié)構(gòu)共振，此時的最大加速度是所有工況下的最大值，并且跑步頻率越接近3.3 Hz，結(jié)構(gòu)跨中加速度響應峰值越大。

(a) 第一次

(b) 第二次

表5 不同頻率跑步激勵下連廊結(jié)構(gòu)的實測最大加速度響應

3.3 數(shù)值模擬與試驗結(jié)果對比

為驗證所建荷載模型，首先對比模擬荷載時程與實測荷載時程。如圖8所示，可見荷載時程的形狀相似，幅值相差也較小。為進一步驗證，建立連廊的有限元模型并計算結(jié)構(gòu)在模擬荷載作用下的加速度響應，圖9為連廊有限元模型及跑步路線示意圖。模型基頻為6.65 Hz，與實測基頻6.6 Hz相近。計算不同跑步頻率下DLF分別取均值、均值±標準差時連廊在模擬荷載作用下的跨中加速度響應峰值[15-18]，跑步頻率增量取0.1 Hz，計算結(jié)果如圖10所示，結(jié)構(gòu)的實測加速度響應峰值與DLF取均值時的計算結(jié)果大體接近，且?guī)缀醵悸湓诎j線(DLF取均值±標準差)內(nèi)。

(a) 模擬3.0 Hz跑步荷載時程

(c) 實測3.0 Hz跑步荷載時程

圖9 連廊有限元模型和荷載布置

圖10 模型峰值加速度響應包絡線與實測結(jié)構(gòu)響應對比

此外，如果采用美國鋼結(jié)構(gòu)協(xié)會AISC規(guī)范Design Guide 11-1997所推薦的節(jié)奏性荷載，則對于2.75 Hz(該規(guī)范給出的荷載頻率的上限值)，計算結(jié)果見圖10中星號，比實測值大很多。對本例而言采用節(jié)奏性荷載代替跑步荷載是不合適的。由上述分析可得，本文提出的跑步荷載模型與實測荷載較為吻合，所引起的結(jié)構(gòu)響應大體一致，可用于相關(guān)工程的設計與分析。

4 結(jié) 論

跑步荷載是室內(nèi)、屋頂運動場以及人行橋等結(jié)構(gòu)振動舒適度設計需要考慮的重要荷載類型。我國目前尚無跑步荷載模型。本研究利用高精度無線測力鞋墊技術(shù)開展了單人跑步荷載試驗，并提出了傅里葉級數(shù)的荷載模型，由試驗數(shù)據(jù)給出了連續(xù)跑步荷載的前5階動載因子和相位角的建議值。通過將數(shù)值模擬與實際結(jié)構(gòu)的響應對比表明，本文所建議的跑步荷載模型可供相關(guān)人致結(jié)構(gòu)振動的工程設計與分析使用。