辛大波，姜海新，張洪福

(東北林業(yè)大學 土木工程學院，哈爾濱 150040)

0 引 言

風是使單體樹木和森林系統(tǒng)受災的主要非人為因素[1-4]。巨大的強風能量對樹木造成危害，不僅影響生態(tài)系統(tǒng)的穩(wěn)定性[4]，還會給林業(yè)生產(chǎn)帶來嚴重的經(jīng)濟損失[5-8]。保障樹木風致安全的首要問題是掌握強風下樹木-風相互作用特性。特別是，樹木的風致響應。

針對樹木-風相互作用特性這一科學問題，通常采用的研究手段主要包括理論分析[9-10]、數(shù)值模擬[3,11-12]和現(xiàn)場實測。陶嗣巍利用計算機輔助設計軟 件Pro/E（ Pro/Engineer） 中 的 簇 表 與UDF（ User Defined Function）功能，基于參數(shù)化建模思路和分形理論，提出了一種快速建立樹木結(jié)構(gòu)幾何模型的方法并分析了白毛楊的風振動特性[13]。宋曉鶴基于非線性假設，將東北魚鱗云杉模擬成彈性變截面圓錐桿，建立了東北云杉樹風倒非線性力學模型，并對其進行了靜力學分析[14]。王琳將東北云杉樹干近似為橫截面隨著高度成冪函數(shù)變化的彈性桿，一端固定，自由端連有質(zhì)量團，建立了力學模型，分析了云杉的受迫振動響應[15]。崔云靜通過Pro/E 三維繪圖軟件建立了新疆楊的實體模型，并運用有限元軟件Ansys workbench 對模型分別進行靜力和動力分析，得到新疆楊在模擬風載荷作用下的應力、應變和位移變化情況[16]。相比較其他兩種研究手段，現(xiàn)場實測具有真實反映樹木及風場的實際狀況、數(shù)據(jù)可靠的優(yōu)點，是開展相關研究較好的方式。在強風天氣下獲得樹木風致振動實測結(jié)果具有一定難度，到目前為止，大多數(shù)的研究主要是基于低風速條件下的現(xiàn)場實測[17-23]，且針葉樹的風致振動行為是目前現(xiàn)場實測研究關注的主要對象[17-20,23-25]。例如Dupont 等人基于現(xiàn)場實測結(jié)果發(fā)現(xiàn)了邊緣流中湍流風場頻率峰值與海岸松（Pinus pinasterAit）基本振動模式相關性較大[23]。

落葉闊葉樹相對針葉樹而言結(jié)構(gòu)要復雜的多，關于其在強風條件下的風致運動實測研究目前較少。但事實上，落葉闊葉樹與針葉樹一樣都具有很高的風敏感性，其覆蓋面積在各樹種中也占有較高的比重。在歐洲的許多森林中，針葉樹正逐漸被落葉闊葉樹所取代[5]；且大多數(shù)歐洲城市的行道樹均為落葉闊葉樹種[26]。因此，開展強風天氣下落葉闊葉樹的振動響應分析十分必要。

綜上，開展強風天氣下落葉闊葉林的風場特性及其風致運動特性研究具有重要意義。本文以黑龍江帽兒山森林生態(tài)系統(tǒng)國家野外科學觀測研究站的落葉闊葉林為研究對象，實測研究了三次強風時段的林間風場特性，并研究了某一白樺樹的風致動力行為，采用譜分析方法分析了脈動風特性、樹木脈動風荷載及風致響應。本文包括引言在內(nèi)共分為4 個部分：第1 節(jié)描述了樣本樹、現(xiàn)場實測的具體位置及監(jiān)測系統(tǒng)，還簡要介紹了功率譜分析的方法；第2 節(jié)介紹并討論了湍流風場特性、白樺樹的風致響應和風樹相互作用的結(jié)果；第3 節(jié)總結(jié)了我們的研究。

1 試驗概況

1.1 現(xiàn)場實測地點和樣本樹

觀測地點選為黑龍江帽兒山森林生態(tài)系統(tǒng)國家野外科學觀測研究站老爺嶺主站點（東經(jīng)127°40′，北緯45°24′）附近。觀測點位于東南-西北的山谷中，周圍的平均坡度約為9°[27]。其周圍的主要樹種包括榆樹、水曲柳、白樺、山楊和胡桃楸等[27]。觀測對象選為該站點的一棵白樺樹（Betula platyphyllaSuk），如圖1 所示。其中，高50 m 的通量塔位于樣本樹西南方向10 m 處，塔與樣本樹之間無其他樹木遮擋。

圖1 樣本樹Fig. 1 Sample tree

在2016 年到2018 年期間，老爺嶺主站點樹木冠層的平均葉面積指數(shù)（mean leaf area index，LAI）和森林平均林分高度分別約為6.95±0.65 和18～20 m[27-28]，平均斷面積為24.16 m2/hm2[29]。所選的白樺樹觀測對象高23 m，胸徑（即，地面以上1.3 m 處樹干的直徑）0.28 m，樹冠中心高度16.8 m。

1.2 現(xiàn)場觀測系統(tǒng)

觀測系統(tǒng)包括：風速儀、加速度傳感器以及數(shù)據(jù)采集系統(tǒng)。風速儀為英國Gill 公司生產(chǎn)的Gill Wind Master 三維超聲波風速儀，風速測量范圍為0～45 m/s，分辨率為0.01 m/s，風速小于12 m/s 時精度小于1.5%；風向測量范圍為0°～359°，分辨率為0.1°，風速小于12 m/s 時精度為2°。本次試驗采樣頻率設置為10 Hz，長期離線采集。Gill Wind Master 三維超聲波風速儀如圖2 所示。在通量塔塔身的5、10、20、35、50 m 處各安裝一個Gill Wind Master 三維超聲波風速儀，通量塔如圖3 所示。

圖2 Wind Master 三維超聲波風速儀Fig. 2 Wind Master ultrasonic anemometers

圖3 50 m 通量塔Fig. 3 50-meter high flux tower

采用揚州晶明科技有限公司生產(chǎn)的三向加速度傳感器（型號：AIT2500）及JM5981A 數(shù)據(jù)采集系統(tǒng)測量并記錄樹木的搖擺，如圖4 和圖5 所示。加速度傳感器量程為±2.5g，最大沖擊極限±50g，橫向靈敏度小于5%，外形尺寸45 mm×45 mm×50 mm，重量200 g。本試驗在執(zhí)行中，加速度傳感器布設在白樺樹樹干上，距地面高度為7.97 m。試驗時，采集頻率設置為50 Hz，長期離線采集。試驗現(xiàn)場整體測試儀器布置如圖6 所示。

圖4 三向加速度傳感器Fig. 4 Three-axis accelerometer

圖5 JM5981A 多功能監(jiān)測分析系統(tǒng)Fig. 5 JM5981A multifunctional monitoring and analysis system

圖6 試驗現(xiàn)場整體布置示意圖Fig. 6 Layout of the measurement site

1.3 三次強風事件

選取帽兒山實測地點2019 年11 月2 日、11 月5 日和11 月8 日這三個強風天氣中所測的風特性和樹木響應觀測數(shù)據(jù)進行后續(xù)的研究與分析。在不到十天的時間里實測地發(fā)生了三次強風，每個強風天氣的一般氣象條件見表1。三次強風中11 月5 日的強風最強，其次為11 月8 日的。塔頂風速儀（z/h= 2.17，z為測點高度，h為樹高）記錄的11 月5 日強風最大平均風速為11.42 m/s，最大陣風為19.82 m/s。

表1 氣象情況Table 1 Meteorological conditions

1.4 風速及風向的處理計算方法

本試驗所用的Gill 超聲波風速儀可以獲得笛卡爾坐標系下三個方向的風速矢量數(shù)據(jù)，Ux、Uy、Uz。其中，Ux正方向定義為與方位角系統(tǒng)中正北方向一致的方向，Uy正方向定義為從正北方向逆時針旋轉(zhuǎn)90°的方向，Uz正方向為平行于風速儀安裝軸垂直向上。

依照中國規(guī)范，在計算分析中設置基本時距為10 min。則在基本時距內(nèi)的三維平均風速為：

式中N為基本計算時距內(nèi)的風速樣本數(shù)。

風向角定義為風向與正北方的夾角，沿順時針增加。則瞬時風向角及10 min 基本計算時距內(nèi)的風向角，即瞬時風向及平均風向可用下式表示：

基于前述10 min 基本計算時距內(nèi)的風向角，可計算得到順風向風速、橫風向風速及豎向風速分別為：

則順風向脈動風速、橫風向脈動風速及豎向脈動風速可分別表示為：

1.5 譜分析方法

脈動風荷載作用下的樹木風致響應分析方法包括時程法和功率譜密度法。時程法通常需要知道樹木氣動力參數(shù)，如：樹木阻力系數(shù)CD和樹冠垂直于氣流的迎風面積A等，這些參數(shù)的獲得較為困難。功率譜密度法可直接進行樹木結(jié)構(gòu)實測響應分析，是常采用的分析方法[17-18,20,24,30-33]。功率譜密度法的一般過程如圖7 所示[24,31]。

圖7 功率譜密度法示意圖Fig. 7 Schematic of the power spectral method

對于一階線性系統(tǒng)，圖7 所示的傳遞關系可通過方程式[24]表達為：

式中：Pl(f)為 風荷載的功率譜密度； ρ為空氣密度，kg/m3；U為順風向平均風速，m/s；Ta(f)為氣動導納函數(shù)；Pu(f)為風速的功率譜密度。

式中：Pt（f）為樹響應（位移或者加速度）的功率譜密度； |Tm(if)|2為機械傳遞函數(shù)，其中i 為虛數(shù)單位。

在實際中由于CD和A并不容易獲得[34]，且風荷載的作用也并不僅僅是一維的，所以風荷載的功率譜并不能由式（13）得到，在研究樹木在風載下的響應的問題通常采用雷諾應力 τ作為風荷載的度量。雷諾應力是湍流脈動速度引起的動量傳遞，由法向和切向壓力構(gòu)成[18,24]。平均雷諾應力 τˉ可表達為：

式中u、w為順風向和豎向的脈動風速，m/s；同理，下文中的v為橫風向脈動風速，m/s。

由圖7 可知，順流向脈動風速功率譜和樹響應功率譜之間的關系可通過機械傳遞函數(shù) |Tm(if)|2表示，|Tm(if)|2提供了樹木對不同頻率湍流的響應的信息，同時 |Tm(if)|2也是對湍流風場到樹木的能量傳遞效率的一種測量[17-18,20]。通過式（15）的平均雷諾應力表達式以及功率譜估計可得到雷諾應力功率譜Pτˉ(f)。將Pτˉ(f)代 替Pl(f)代 入 式（14），則 機 械 傳 遞 函 數(shù)|Tm(if)|2可通過下式得到：

2 結(jié)果與分析

選取3 次強風過程中z= 2.17h高度處風速傳感器所測得的平均風速大于5 m/s 的時間段進行后續(xù)的所有分析。每個強風天氣內(nèi)所選取的時間段約為5 h。

2.1 平均風速與風向

如圖8 所示，為在三次強風天氣里測量并處理得到的順風向平均風速廓線，其中高度（z）采用z=0.87h（近冠頂）處的高度作為參考高度（zr）進行歸一化，風速（u）以對應的z= 0.87h處的風速作為參考風速（ur）進行歸一化。圖中虛線代表樣本樹冠頂高度。可以發(fā)現(xiàn)在三次強風天氣內(nèi)有著幾乎相同的順風向平均風速廓線。在冠層頂部以上平均風速沿測點高度近似呈現(xiàn)對數(shù)分布。

繪制了3 個強風天氣的平均風速和風向的時程曲線，如圖9 所示，圖中陰影區(qū)域為用于本文分析的強風時段。

從圖9 可以看出，三個強風展示了相似的風特性。從z= 1.25h到0.87h的高度，即到達冠層時，平均風速急劇下降（圖9（a）），這可能是由于冠層樹葉和樹枝等從氣流中提取動量導致的。為了說明該結(jié)論，我們選取2021 年5 月7 日研究區(qū)域附近某一空曠區(qū)域（無樹木影響區(qū)域）的風場數(shù)據(jù)與本文數(shù)據(jù)進行對比分析。2021 年5 月7 日研究區(qū)域附近某一空曠區(qū)域的20 m、35 m 和50 m 的平均風速時程曲線如圖10所示。

通過對比圖10 和圖9（a）可以發(fā)現(xiàn)，沒有樹木影響的區(qū)域下，隨著測點高度的降低平均風速稍有減小，幅值很??；從z= 1.25h到0.87h，風速并沒有急劇減小的現(xiàn)象。而在有樹木的影響下，從z= 1.25h到0.87h，即到達冠層時，平均風速急劇下降（圖9（a）），可知這是由于冠層樹葉和樹枝等從氣流中提取動量導致的。

圖9 三次強風期間的風特性Fig. 9 Wind dynamics during the three wind storms

圖10 2021 年5 月7 日平均風速時程Fig. 10 Time course of average wind speed on May 7, 2021

在冠層下面的兩個高度（即z= 0.43h和0.22h）的風速幾乎相同。通過表1 和圖9 可以發(fā)現(xiàn)對于同一測點高度，在風速較大的強風期段，平均風向變動范圍較小，基本保持穩(wěn)定；而在z= 2.17h～1.52h高度處平均風速近似小于2.0 m/s（z= 0.87h高度處平均風速近似小于1.0 m/s，z= 0.43h～0.22h高度處平均風速近似小于0.5 m/s）的低風速時段內(nèi)，風向變化范圍很大。還可以發(fā)現(xiàn)平均風速方向隨高度方向發(fā)生偏轉(zhuǎn)，冠層高度范圍以內(nèi)，偏轉(zhuǎn)梯度大，冠層高度以外，偏轉(zhuǎn)梯度較??；最大風偏轉(zhuǎn)角約為67°，從z= 2.17h處的西北風轉(zhuǎn)變?yōu)閦= 0.22h處的西南風，這是由于樹木的擾動作用對風特性產(chǎn)生了影響。三次強風時段的平均風向廓線如圖11 所示。

圖11 平均風向廓線Fig. 11 Average wind direction profiles

2.2 脈動風速功率譜的分析

脈動風速功率譜反映脈動風能量隨頻率變化的分布情況，是非常重要的脈動風特性之一。當脈動風頻率接近或等于高聳結(jié)構(gòu)的固有頻率時，湍流脈動分量的大小將會顯著影響結(jié)構(gòu)的響應。脈動風速功率譜包括水平縱向（順風向）脈動風速功率譜、水平橫向（橫風向）脈動風速功率譜及豎向脈動風速功率譜。在這項研究中不同測點高度處順風向、橫風向和豎向的歸一化時間平均脈動風速功率譜如圖12所示。

三次強風事件中的三個風速分量（u、v和w）在z= 1.52h和2.17h高度處的時間平均功率譜展現(xiàn)出與常見的大氣表層譜相似的形狀，即功率譜含能區(qū)的斜率近似以f+1的 函數(shù)上升，在慣性子區(qū)近似以f-2/3的函數(shù)下降（圖12）。三次強風事件中，相比其他測點的u-譜，z= 0.22h測點的u-譜在耗散區(qū)的斜率最陡。這是因為冠層以下植被背后的尾跡結(jié)構(gòu)發(fā)展，導致湍流耗散加速[23]。

由圖12（a）可知，順風向脈動風速功率譜（Su）的歸一化峰值頻率（fpu，即最大功率對應的頻率）隨風暴強度的增加而減??；對于較強的風暴（11 月5 日和11 月8 日）來說，順風向脈動風速功率譜（Su）的最大功率隨測點高度的減小而增加，結(jié)合圖12（a）的順風向湍流強度廓線圖可以發(fā)現(xiàn)這是由于樹木的擾動作用增加了冠層內(nèi)風場的湍流強度造成的。李正農(nóng)等人通過風洞試驗也提出樹木的背風面存在風影響區(qū)，此區(qū)域的流場特性會發(fā)生顯著改變，體現(xiàn)為風速降低和湍流強度增大[35-36]。由圖12（b）可知，橫風向脈動風速功率譜（Sv）的歸一化峰值頻率（fpv）隨風暴強度的變化規(guī)律與Su的相同；Sv的最大功率，對于較強的風暴（11 月5 日和11 月8 日），除了測點z= 0.22h，隨測點高度的減小而增加；對于較弱的風暴（11 月2 日），各測點的最大功率較為接近，無明顯規(guī)律。各個風暴事件的豎向脈動風速功率譜（Sw）（圖12（c））的最大功率在各測點之間較為接近，無明顯規(guī)律。

圖12 歸一化時間平均風速功率譜及順風向湍流強度廓線圖Fig. 12 Normalized time-averaged power spectra of the wind-speed components and turbulence intensity profile of stream-wise direction

總的來說，觀察三個強風事件的u-譜和v-譜發(fā)現(xiàn)，相鄰兩個測點之間的最大功率的差值隨風暴強度的增加而增加。同時，對于同一強風事件的同一速度分量在不同高度處的功率譜歸一化峰值頻率幾乎相同。11 月2 日、5 日和8 日強風的u、v和w譜的歸一化峰值頻率分別為f=nz/u= 0.33、0.35±0.01 和0.60±0.09、f=nz/u= 0.28±0.02、0.28±0.02 和0.67±0.10 和f=nz/u= 0.30±0.02、0.29 和0.57±0.01?；诿}動風速譜及歸一化頻率公式整理可得各測點高度的脈動風峰值頻率，發(fā)現(xiàn)三次強風期間各測點的順風向和橫風向脈動風速（u和v）的峰值頻率約為0.03 Hz，豎向脈動風速（w）的峰值頻率在0.05～0.07 Hz 之間，不同測點高度之間有細微差別。

對同一強風事件同一測點高度下三個風速分量（u、v和w）的功率譜密度圖進行了繪制，如圖13所示。

圖13 三個風速分量的歸一化時間平均功率譜Fig. 13 Normalized time-averaged power spectra of the three wind-speed components (u, v and w)

對于11 月2 日的強風事件，5 個測點高度的功率譜具有相似的變化規(guī)律，由u→v→w，峰值頻率越來越大，最大功率越來越小。對于11 月5 日和8 日的強風時段，所有測點高度的峰值頻率變化規(guī)律具有相似變化規(guī)律，即fpu≈fpv＜fpw。11 月5 日的最大功率，除z= 0.87h測點外，由u→v→w，越來越小。11 月8 日的最大功率，對于測點z= 0.22h和0.43h，由u→v→w，越來越?。粚τ跍y點z= 0.87h、1.52h和2.17h，v-譜最大，u-譜次之。三個風速分量的不同的功率譜密度表明了林分樹冠上部、頂部及冠層內(nèi)湍流的各向異性特征?？傮w看來，隨風暴強度的增加，湍流各向異性增加；對于同一風暴，湍流各向異性隨測點高度的增加而減小。

2.3 樹響應及風荷載的分析

用三次強風期間樹干加速度分量（東北-西南和東南-西北方向）來表征風荷載作用下樣本樹的運動，如圖14 所示為三次強風期間樣本樹加速度響應的時程曲線。通過圖14 可以發(fā)現(xiàn)，在三次強風期間，白樺樹均呈現(xiàn)出不規(guī)則的復雜的搖擺運動。比較三個強風天氣里樣本樹的運動大小可以發(fā)現(xiàn)，隨著風強度的增加（“11 月2 日”＜“11 月8 日”＜“11 月5 日”），白樺樹在兩個方向的運動變大。且在強風的影響下，白樺在東北-西南和東南-西北兩個方向的運動大小相近，可知樣本樹在強風期間的運動軌跡近似為圓形。

圖14 加速度響應時程Fig. 14 Time series of the tree acceleration response

z= 2.17h、1.52h、0.87h、0.43h和0.22h的5 個測點處的雷諾應力功率譜、基于z= 0.87h測點處的雷諾應力功率譜和樹木加速度響應功率譜得到的機械傳遞函數(shù)以及z= 0.35h高度處的樹木水平縱向和橫向加速度響應功率譜分別如圖15～圖17 所示。

圖17 機械傳遞函數(shù)Fig. 17 Mechanical transfer function

通過圖15 可以發(fā)現(xiàn)，總體而言，在三次強風期間各測點低頻區(qū)雷諾應力功率譜密度增加，到達功率譜密度的最大值后，又逐漸減小，但隨著風暴強度的增加，減小的斜率越來越平緩。但在11 月5 日的強風暴期間，在z= 0.87h測點處（近冠層頂部），當功率譜密度值達到最大后，在一定頻率范圍內(nèi)，呈現(xiàn)出一段功率譜密度幾乎在相同水平的頻率區(qū)間。11 月2 日、5 日和8 日的τ-譜的歸一化峰值頻率分別為f=nz/u= 0.58 ± 0.08、0.84 ± 0.30 和0.55 ± 0.09。

圖15 雷諾應力的歸一化時間平均功率譜Fig. 15 Normalized time-averaged power spectra of the Reynolds stress

如圖16 所示，在11 月2 日、5 日和8 日的強風期間a1和a2譜的歸一化峰值振動頻率分別為f=nz/u=3.44 和3.28、2.67 和2.55 以及3.34 和3.19。同一強風期間，白樺樹兩個加速度分量的功率譜差異很小，這與Peltola 等對黑松進行測量后得到的結(jié)果相同[19]，也與Mayer 得到的結(jié)果一致[24,31]。Mayer 曾指出這種較小的差異性可能是由樹木樹冠結(jié)構(gòu)部分的不對稱性分布、平均風向以及根莖的發(fā)展導致的[37]?；诩铀俣软憫V及歸一化頻率公式整理可得到白樺樹在三次強風中a1和a2譜的峰值振動頻率分別為0.30 和0.28、0.27 和0.26 以及0.31 和0.29 Hz，不同強風天氣之間有細微差別。

圖16 加速度響應的歸一化時間平均功率譜Fig. 16 Normalized time-averaged power spectra of the tree acceleration

在三個強風的加速度響應譜中還發(fā)現(xiàn)了兩加速度分量基本一致的高于峰值振動頻率的其他兩個顯著的峰值，這可能對應于白樺樹的二階和三階諧波頻率，分別約為0.76 和1.72 Hz、0.69 和1.52 Hz 以及0.65 和1.65 Hz。Peltola 等人提出相比于間伐過的林分中的黑松，未間伐的林分中的黑松因具有相對大的靈活性的小直徑的樹木元素而使入射能量被劃分為更高的共振模式[19]。Scannell 也提出這種更高階的振動模式可能是由于枝葉等樹木元素造成[38]。

為了得到樹木自由振動頻率，還進行了樹木拉力試驗，并對相應的加速度響應數(shù)據(jù)進行傅里葉變換分析后，得到白樺樹自由振動下兩個加速度分量的峰值頻率，分別為0.28 Hz 和0.26 Hz，與強風下白樺樹的峰值振動頻率一致性較好。Mayer、Amtmann 等對針葉樹進行研究也得到了相同的結(jié)果[17-18,31-32,39]。

如圖17 所示，三次強風下白樺樹的機械傳遞函數(shù)均在白樺樹加速度響應功率譜的峰值頻率處達到峰值，然后在更高頻率處下降。這與Gardiner、Stacey 等對針葉樹在風下進行測試得到的結(jié)果相同；Gardiner、Stacey 等表明，通過將針葉樹視為具有單個集總質(zhì)量的簡單阻尼諧振子并使用冠層質(zhì)量、莖質(zhì)量、莖彈性和阻尼比的測量值對風致樹木響應進行模擬后，得到的模擬機械傳遞函數(shù)與通過試驗得到的樹響應的機械傳遞函數(shù)非常吻合，也就是說，對針葉樹在風下進行測試得到風致樹木響應的機械傳遞函數(shù)幾乎與受迫阻尼諧振子的傳遞函數(shù)相同[17-18,30,32,39]。總的來說，森林可以被看作是一組有效地吸收樹木共振頻率處的湍流譜能量的諧波振蕩器[17]。此外，和加速度響應譜一樣，兩個加速度分量的機械傳遞函數(shù)之間只有微小的差異。

3 結(jié) 論

以黑龍江帽兒山森林生態(tài)系統(tǒng)國家野外科學觀測研究站的落葉闊葉林為研究對象，實測研究了三次強風時段的林間平均風與脈動風特性，分析了林中某一白樺樹的脈動風荷載及風致響應。主要得到以下幾點結(jié)論：

1）樹木的擾動作用導致平均風隨高度發(fā)生偏轉(zhuǎn)，在冠層高度范圍以內(nèi)偏轉(zhuǎn)梯度大，冠層高度以外偏轉(zhuǎn)梯度較小；最大風偏轉(zhuǎn)角度近似為67°。較大的風偏轉(zhuǎn)角可能對林木或者林中結(jié)構(gòu)產(chǎn)生扭曲作用，在開展林木抗風分析、防風減災能力評估以及林中柔性結(jié)構(gòu)物抗風設計等方面應引起足夠的重視。

2）樹木的擾動作用增加了冠層內(nèi)風場的湍流強度，這一影響隨著來流風強度的增加而增加。三次強風作用下，樹木在脈動風作用下呈現(xiàn)出復雜的隨機振動，振動響應的峰值頻率與樹木自由振動的峰值頻率一致，與來流風的脈動風特性無關。

3）闊葉樹的機械傳遞函數(shù)的峰值頻率與響應功率譜的峰值頻率基本一致，樹木風致響應的機械傳遞函數(shù)幾乎與受迫阻尼諧振子的傳遞函數(shù)相同，這一點與針葉樹具有相似的特征。