留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

稳定风速和阵风交错作用下的林木风倒分析

贾杰

贾杰. 稳定风速和阵风交错作用下的林木风倒分析[J]. 浙江农林大学学报, 2014, 31(4): 604-610. doi: 10.11833/j.issn.2095-0756.2014.04.017
引用本文: 贾杰. 稳定风速和阵风交错作用下的林木风倒分析[J]. 浙江农林大学学报, 2014, 31(4): 604-610. doi: 10.11833/j.issn.2095-0756.2014.04.017
JIA Jie. Tree windthrow with stable winds and gusts[J]. Journal of Zhejiang A&F University, 2014, 31(4): 604-610. doi: 10.11833/j.issn.2095-0756.2014.04.017
Citation: JIA Jie. Tree windthrow with stable winds and gusts[J]. Journal of Zhejiang A&F University, 2014, 31(4): 604-610. doi: 10.11833/j.issn.2095-0756.2014.04.017

稳定风速和阵风交错作用下的林木风倒分析

doi: 10.11833/j.issn.2095-0756.2014.04.017
基金项目: 

黑龙江省教育厅科学技术研究项目 12513033

详细信息
    作者简介: 贾杰, 讲师, 从事非线性动力学的分析与控制研究。E-mail:jiajiecontrol@126.com
  • 中图分类号: S718.5;S727.2

Tree windthrow with stable winds and gusts

  • 摘要: 风倒是形成林隙的主要原因, 是发生频率较高、波及范围较广、造成损失较严重的一种自然灾害。因此, 研究风倒具有重要的理论和实践意义。针对相关理论模型的不足, 利用动力学理论建立了林木风倒的非线性动态模型, 分别从稳定风速作用和阵风作用2个方面对林木的风倒机理进行了研究, 通过严格的数学推导获得了林木发生风倒的判断条件。结果表明:稳定风速和阵风作用下林木发生风倒的条件不同, 且互为补充。基于理论分析, 以西加云杉Picea sitchensis为例进行数值计算, 得到了风倒的标准风速值, 研究结果可以应用到其他树种的风倒问题中。同时, 稳定风速和阵风交错作用下发生风倒的标准风速值与还与风荷载的作用时间有关。研究结果为后续的理论研究和工程应用提供了指导。
  • 图  1  树木基本受力图

    Figure  1  Basic force diagram of trees

    图  2  总体抵抗力矩和倾角θ的关系

    Figure  2  Relationship between whole resisting moment and angle θ

    图  3  抵抗力矩mR(θ)的模型图

    Figure  3  Mode diagram of resisting moment mR(θ)

    图  4  式(9)的函数图像[a表示在条件(10)满足的时候,b表示条件(10)不满足的情况]

    Figure  4  Function diagram of the formula (a show the condition (10) is satisfied, b show the condition (10) is dissatisfied)

    图  5  引起风倒的阵风风速和连续风速的函数关系(在τ=0.2,0.5,1,2,5,10s的情况下)

    Figure  5  Function relationship between gust wind and continuous wind about windthrow (in the case of τ=0.2, 0.5, 1, 2, 5, 10s)

    图  6  紊流强度为0.2的时候破坏的平均, 风速图

    Figure  6  Diagram of average wind speed when turbulence intensity is 0.2

    图  7  紊流强度为0.2,τ=1,2,5s破坏的平均风速图

    Figure  7  Average wind speed diagram of damage when turbulence intensity is 0.2, τ=1, 2, 5s

  • [1] 于振良, 赵士洞.林隙(gap)模型研究进展[J].生态学杂志, 1996, 16(2):42-46.

    YU Zhenliang, ZHAO Shidong. Research progress on gap model[J]. Chin J Ecol, 1996, 16(2):42-46.
    [2] 李国旗, 安树青, 张纪林, 等.海岸带防护林4种树木的风压应力分析[J].南京林业大学学报, 1999, 23(4):76-80.

    LI Guoqi, AN Shuqing, ZHANG Jilin, et al. The bending stress analysis of 4 species of woods caused by wind pressure in coastal shelter forest[J]. J Nanjing For Univ, 1999, 23(4):76-80.
    [3] 陈少雄, 王观明, 罗建中.桉树幼林不同株行距配置抗台风效果的研究[J].林业科学研究, 1995, 8(5):582-585.

    CHEN Shaoxiong, WANG Guanming, LUO Jianzhong. Effect of different spacing on typhoon resistance of young Eucalyptus plantation[J]. For Res, 1995, 8(5):582-585.
    [4] 陈少雄, 杨民胜, 王理平.尾叶桉造林密度与蓄积量、抗风和材性关系研究[J].林业科学研究, 1998, 11(4):435-438.

    CHEN Shaoxiong, YANG Minsheng, WANG Liping. Effect of spacing on volume, storm-resistance and wood quality of Eucalyptus urophylla[J]. For Res, 1998, 11(4):435-438.
    [5] 陈士银.木麻黄不同株行距配置抗台风效果[J].广东林业科技, 1997, 13(2):44-46.

    CHEN Shiyin. Effects of spacing on the resisitance to typhoon in Casuarina equisetifolia forest belt[J]. J Guangdong For Sci Technol, 1997, 13(3):44-46.
    [6] 代力民, 徐振邦, 杨丽韫.红松阔叶林倒木储量动态的研究[J].应用生态学报, 1999, 10(5):513-517.

    DAI Limin, XU Zhenbang, YANG Linyun, et al. Storage dynamics of fallen trees in Korean pine broad-leaved forest[J]. Chin J Appl Ecol, 1999, 10(5):513-517.
    [7] 宋晓鹤.云杉风倒静力学模型的建立及其分析[D].哈尔滨:哈尔滨工业大学, 2006.

    SONG Xiaohe. Builing and Analysis of Mechanics Model for Windthrow of Spruce[D]. Harbin:Harbin Institute of Technology, 2006.
    [8] 王琳.云杉风倒动力学模型的建立与分析[D].哈尔滨:哈尔滨工业大学, 2006.

    WANG Lin. Building and Analysis of Dynamic Model for Windthrow of Spruce[D]. Harbin:Harbin Institute of Technology, 2006.
    [9] 赖秋明.云杉风倒动力学问题的研究[D].哈尔滨:哈尔滨工业大学, 2008. LAI Qiuming.

    Dynamics Research for Windthrow of Spruce[D]. Harbin:Harbin Institute of Technology, 2008.
    [10] 吴志华, 李天会, 张华林, 等.沿海防护林树种木麻黄和相思生长和抗风性状比较研究[J].草业学报, 2010, 19(4):166-175.

    WU Zhihua, LI Tianhui, ZHANG Hualin, et al. Studies on growth and wind-resistance traits of Casuarina and Aacia stands from coastal protection forest[J]. Acta Pratac Sin, 2010, 19(4):166-175.
    [11] 段溪.亚热带泥质海岸68个防护林树种抗风性能比较分析[D].南京:南京林业大学, 2012.

    DUAN Xi. A Comparative Analysis of Wind Resistance for 68 Shelter Tree Species in Subtropical Muddy coast[D]. Nanjing:Nanjing Forestry University, 2012.
    [12] RUEL J C, PIN D, SPACE L, et al. The estimation of wind exposure for windthrow hazard rating:comparison between Strong blow, MC2, topex and a wind tunnel study[J]. Forestry, 1997, 70(3):254-266.
    [13] ENGLAND A H, BAKER D J, SAUNDERSON S E T. A dynamic analysis of windthrow of trees[J]. Forestry, 2000, 73(3):225-237.
    [14] BENOIT R, DESGAGNÉM, PELLERIN P, et al. The Canadian MC2:a semi-lagrangian, semi-implicit wideband atmospheric model suited for fine-scale process studies and simulation[J]. Mon Weather Rev, 1997, 125(10):2382-2415.
    [15] ACHIM A, RUEL J C, GARDINER B A, et al. Modeling the vulnerability of balsam fir forests to wind damage[J]. For Ecol Manage, 2005, 204:35-50.
    [16] 贾杰, 李静辉.林木风倒动态模型的建立与分析[J].森林工程, 2013, 29(4):63-67.

    JIA Jie, LI Jinghui. Establishment and analysis of dynamical model for windthrow of trees[J]. For Eng, 2013, 29(4):63-67.
    [17] 贾杰, 李洪峰, 李静辉.基于Von-Karman谱的独立林木风倒机理的频率域分析[J].东北林业大学学报, 2013, 41(8):165-169.

    JIA Jie, LI Hongfeng, LI Jinghui. Frequency domain analysis of independent trees windthrow based on Von-Karman spectrum[J]. J Northeast For Univ, 2013, 41(8):165-169.
    [18] MERRENS E J, PEART D R. Effects of hurricane damage on individual growth and stand structure in a hardwood structure in a hardwood forest in New Hampshire, USA[J]. J Ecol, 1992, 80(4):787-795.
  • [1] 陈涵, 张超, 余树全.  基于Stacking模型集成算法的莲都区南方红豆杉潜在分布区 . 浙江农林大学学报, 2019, 36(3): 494-500. doi: 10.11833/j.issn.2095-0756.2019.03.009
    [2] 李军, 陆云峰, 杨安娜, 柏志靓, 王洋, 楼炉焕, 童再康.  紫楠天然群落物种多样性对不同干扰强度的响应 . 浙江农林大学学报, 2019, 36(2): 279-288. doi: 10.11833/j.issn.2095-0756.2019.02.009
    [3] 邱问心, 张勇, 俞佳骏, 张超, 郑超超, 余树全.  InVEST模型水源涵养模块实地应用的可行性验证 . 浙江农林大学学报, 2018, 35(5): 810-817. doi: 10.11833/j.issn.2095-0756.2018.05.004
    [4] 朱宇颐, 解潍嘉, 黄华国.  基于三维模型ENVI-met对黑河森林和北方森林的潜热及显热通量模拟 . 浙江农林大学学报, 2018, 35(3): 440-452. doi: 10.11833/j.issn.2095-0756.2018.03.007
    [5] 陈丽萍, 李平衡, 莫路锋, 周国模, 李金荣.  基于通量源区模型的雷竹林生态系统碳通量信息提取 . 浙江农林大学学报, 2016, 33(1): 1-10. doi: 10.11833/j.issn.2095-0756.2016.01.001
    [6] 何涛, 孙玉军.  基于InVEST模型的森林碳储量动态监测 . 浙江农林大学学报, 2016, 33(3): 377-383. doi: 10.11833/j.issn.2095-0756.2016.03.002
    [7] 杨爱芳, 韩有志, 杨秀清, 刘泽民, 杨莎, 王官.  不同密度云杉林下草本植物多样性 . 浙江农林大学学报, 2014, 31(5): 676-682. doi: 10.11833/j.issn.2095-0756.2014.05.003
    [8] 李妍, 李登武, 韩东辰.  贺兰山东坡青海云杉种群动态 . 浙江农林大学学报, 2014, 31(1): 50-56. doi: 10.11833/j.issn.2095-0756.2014.01.008
    [9] 李菲菲, 饶良懿, 吕琨珑, 李会杰, 宋丹丹.  Priestley-Taylor模型参数修正及在蒸散发估算中的应用 . 浙江农林大学学报, 2013, 30(5): 748-754. doi: 10.11833/j.issn.2095-0756.2013.05.018
    [10] 魏书精, 孙龙, 魏书威, 胡海清.  森林生态系统粗木质残体研究进展 . 浙江农林大学学报, 2013, 30(4): 585-598. doi: 10.11833/j.issn.2095-0756.2013.04.019
    [11] 王金亮, 程鹏飞, 徐申, 王小花, 程峰.  基于遥感信息模型的香格里拉森林生物量估算 . 浙江农林大学学报, 2013, 30(3): 325-329. doi: 10.11833/j.issn.2095-0756.2013.03.003
    [12] 俞静芳, 余树全, 张超, 李土生.  应用CASA模型估算浙江省植被净初级生产力 . 浙江农林大学学报, 2012, 29(4): 473-481. doi: 10.11833/j.issn.2095-0756.2012.04.001
    [13] 杨前宇, 谢锦忠, 张玮, 林振清.  椽竹各器官生物量模型 . 浙江农林大学学报, 2011, 28(3): 519-526. doi: 10.11833/j.issn.2095-0756.2011.03.027
    [14] 张利阳, 温国胜, 王圣杰, 刘兆玲.  毛竹光响应模型适用性分析 . 浙江农林大学学报, 2011, 28(2): 188-193. doi: 10.11833/j.issn.2095-0756.2011.02.003
    [15] 刘烁, 周国模, 白尚斌.  基于光照强度变化的毛竹扩张对杉木影响的探讨 . 浙江农林大学学报, 2011, 28(4): 550-554. doi: 10.11833/j.issn.2095-0756.2011.04.005
    [16] 高铭聪, 蒋文伟, 金竹秀, 郭慧慧, 梅艳霞.  西径山森林公园夏季空气负离子日变化 . 浙江农林大学学报, 2011, 28(4): 667-673. doi: 10.11833/j.issn.2095-0756.2011.04.023
    [17] 接程月, 辛赞红, 信晓颖, 江洪, 魏晓华.  FORECAST模型的原理、方法和应用 . 浙江农林大学学报, 2009, 26(6): 909-915.
    [18] 许慧敏, 张劲峰, 向红梅, 郭华.  不同强度人为干扰对香格里拉亚高山森林树种结构的影响 . 浙江农林大学学报, 2008, 25(5): 591-596.
    [19] 哀建国, 俞琳, 章丽英, 钱柳钦, 张腾超.  浙江雪胆群落学特征研究 . 浙江农林大学学报, 2007, 24(6): 706-710.
    [20] 何莹, 韦新良, 蔡霞, 李可追, 王珍.  生态景观林群落结构定量分析 . 浙江农林大学学报, 2007, 24(6): 711-718.
  • 加载中
  • 链接本文:

    https://zlxb.zafu.edu.cn/article/doi/10.11833/j.issn.2095-0756.2014.04.017

    https://zlxb.zafu.edu.cn/article/zjnldxxb/2014/4/604

图(7)
计量
  • 文章访问数:  2626
  • HTML全文浏览量:  413
  • PDF下载量:  533
  • 被引次数: 0
出版历程
  • 收稿日期:  2013-11-29
  • 修回日期:  2014-02-27
  • 刊出日期:  2014-08-20

稳定风速和阵风交错作用下的林木风倒分析

doi: 10.11833/j.issn.2095-0756.2014.04.017
    基金项目:

    黑龙江省教育厅科学技术研究项目 12513033

    作者简介:

    贾杰, 讲师, 从事非线性动力学的分析与控制研究。E-mail:jiajiecontrol@126.com

  • 中图分类号: S718.5;S727.2

摘要: 风倒是形成林隙的主要原因, 是发生频率较高、波及范围较广、造成损失较严重的一种自然灾害。因此, 研究风倒具有重要的理论和实践意义。针对相关理论模型的不足, 利用动力学理论建立了林木风倒的非线性动态模型, 分别从稳定风速作用和阵风作用2个方面对林木的风倒机理进行了研究, 通过严格的数学推导获得了林木发生风倒的判断条件。结果表明:稳定风速和阵风作用下林木发生风倒的条件不同, 且互为补充。基于理论分析, 以西加云杉Picea sitchensis为例进行数值计算, 得到了风倒的标准风速值, 研究结果可以应用到其他树种的风倒问题中。同时, 稳定风速和阵风交错作用下发生风倒的标准风速值与还与风荷载的作用时间有关。研究结果为后续的理论研究和工程应用提供了指导。

English Abstract

贾杰. 稳定风速和阵风交错作用下的林木风倒分析[J]. 浙江农林大学学报, 2014, 31(4): 604-610. doi: 10.11833/j.issn.2095-0756.2014.04.017
引用本文: 贾杰. 稳定风速和阵风交错作用下的林木风倒分析[J]. 浙江农林大学学报, 2014, 31(4): 604-610. doi: 10.11833/j.issn.2095-0756.2014.04.017
JIA Jie. Tree windthrow with stable winds and gusts[J]. Journal of Zhejiang A&F University, 2014, 31(4): 604-610. doi: 10.11833/j.issn.2095-0756.2014.04.017
Citation: JIA Jie. Tree windthrow with stable winds and gusts[J]. Journal of Zhejiang A&F University, 2014, 31(4): 604-610. doi: 10.11833/j.issn.2095-0756.2014.04.017
  • 在众多自然灾害中,风荷载作用所引起的林木风倒是最常见一种,同时风倒也是形成林隙的主要原因,逐渐成为森林生态的一个扰动因子。风倒现象主要表现为干折、根折、连根拔起等3种方式。全世界每年非灾害性风倒造成的木材减产在15%以上[1]。近年来,人类对自然环境的破坏日益严重,导致全球气候环境的急剧变化,风倒引发的森林破坏频繁发生。因此,展开树木风倒机制的相关研究日益受到科学研究人员的重视,研究结果对于森林生态系统和木材生产的理论研究和生产实践都有极其重要的意义。迄今为止,国内外很多专家学者开展了林木风倒机制方面的研究,并取得了丰富的研究成果。其中,李国旗等[2]研究了位于不同风压和高度处树木的应力分布规律。陈少雄等[3-4]、陈士银[5]基于统计和野外调查的方法分别研究了泡桐Paulownia tomentosa,桉树Eucalyptus robusta,木麻黄Casuarina equisetifolia等树种的抗风性与个体特征之间的关系。代力民等[6]从生态学角度出发研究林隙的形成原因。宋晓鹤[7],王琳[8],赖秋明[9]分别从静力学和动力学角度建立了风倒的模型,并展开了相应的分析。吴志华等[10]进行了抗风性状对抗风效果的影响规律的研究。段溪[11]运用层次分析法研究了杭州湾地区亚热带泥质海岸的68种防护树种的抗风性能。英国、加拿大、美国等一些学者对樟子松Pinus sylvestris var. mongolica等树种的风倒机制进行了深入广泛的研究,以动力学理论,结合野外试验中林木和土壤力学参数的测定,地形数字高程模型和实物模型风洞试验,获得了复杂地形下和不同林分密度时的风场分布,建立了一系列力学方程来描述林木的抗风性[12],对风倒的力学机制也进行了相应的阐述,形成若干关于林木风倒机制的模型,如MC2[13]和HWIND[14]等。尽管这些研究代表了风倒机制研究的国际先进水平,但主要是基于材料力学方法的线性力学模型[15]。虽然描述了风倒机制的许多方面,但受研究条件等的限制在系统设计时对关键环节进行了简化,造成模型与实际情况之间的误差。作者在已有研究成果的基础上,建立了林木发生风倒的非线性动力学模型[16],并且利用Von-Karman风谱模型从频率域内分析了林木风倒机制[17],验证了模型的有效性。相比其他模型,这个模型更接近实际物理系统且考虑了实际林木发生风倒的一些不确定性影响因素。但模型也存在一些缺陷,首先,并未考虑树冠的形状、树木个体差异、枝条和树叶的作用等对风倒的影响,其次,也未考虑风荷载作用下的“流固耦合”现象,这也是后续研究的方向和内容。本研究在前期的研究基础上,着重考虑孤立木在稳定风速和阵风作用下所引起的倾覆,通过严格的数学推导给出风倒发生的条件,并以西加云杉Picea sitchenrsis为例,进行了数值计算,给出树木发生风倒的标准风速值。研究结果对森林抗风防护建设和经营提供现实的指导意义。

    • 当风荷载作用于树木的时候,空气动力将会对树杆和树枝产生一个对于基础部分的倾覆力矩,由根部所产生的力矩与其抵抗(抵抗力矩)。我们将倾覆力矩定义为mT(θ),抵抗力矩定义为mR(θ),并且假设树为刚性体,我们可以得到如下的模型:

      $$ {I_0}\theta ''\left( t \right) = {m_{\rm{T}}}\left( \theta \right)-{m_{\rm{R}}}\left( \theta \right)。 $$ (1)

      式(1)中:I0表示树木相对根部的转动惯量,θ表示树杆在风力作用下相对于垂直位置的倾角。如图 1所示。以西加云杉为例,通过对其根部的抵抗力矩进行测量,发现抵抗力矩取决于根部铰接部分的抵抗矩、土壤的压力、迎风面根系的拉力、根系土壤的质量等等,如图 2所示。

      图  1  树木基本受力图

      Figure 1.  Basic force diagram of trees

      图  2  总体抵抗力矩和倾角θ的关系

      Figure 2.  Relationship between whole resisting moment and angle θ

      图 2可知:总体抵抗力矩是关于倾角θ的函数,是很多树木拉力试验中的典型曲线。并且当θ很小的时候,总体抵抗力矩急剧上升,在2°到5°左右,抵抗力矩达到最大值,在10°的时候,抵抗力矩降到其最大抵抗力矩的一半左右。基于本研究的目的,可以将抵抗力矩模拟成一个二次方程,在θ=θm的时候有最大抵抗力矩mR(θ)=Mm,在θ=0时,定义mR(θ)=M0。因此可以得到关于mR(θ)的方程:

      $$ {m_{\rm{R}}}\left( \theta \right) = {M_{\rm{m}}}-\frac{{\left( {{M_{\rm{m}}}-{M_{\rm{0}}}} \right){{\left( {{\theta _{\rm{m}}}-\theta } \right)}^2}}}{{{\theta ^2}_{\rm{m}}}}。 $$ (2)

      在0°到0.5°之间,由于其mR(θ)曲线陡然增加,近似用割线代替,从坐标轴上选取(0, M0)作为起点,得抵抗力矩模型曲线如图 3所示。

      图  3  抵抗力矩mR(θ)的模型图

      Figure 3.  Mode diagram of resisting moment mR(θ)

    • 图 1所示,空气动力对于树木一个微小的单元作用力为νn2dAνn即对单元作用的风速,dA表示树干上具有代表性的单元面积。假设风具有连续的水平速度u0,当树木垂直的时候,风力对树木根部产生的弯矩为Fm(u0),当树木在风力作用下发生倾斜,倾角为θ的时候,其弯矩变为Fm(u0)cos2(θ),同时重心也将发生偏移,自身重力将产生一个倾覆力矩,定义为Wmsin(θ)。Wm表示θ=90°时,树木自身重力对根部所产生的力矩。则方程变为:

      $$ {I_0}\theta ''\left( t \right) = {F_{\rm{m}}}\left( {{u_0}} \right){\cos ^2}\left( \theta \right) + {W_{\rm{m}}}\sin \left( \theta \right) - {M_{\rm{m}}} + \frac{{\left( {{M_{\rm{m}}} - {M_{\rm{0}}}} \right){{\left( {{\theta _{\rm{m}}} - \theta } \right)}^2}}}{{\theta _{\rm{m}}^2}}。$$ (3)

      当根部抵抗力矩较小时,可以近似地认为:sinθθ,cos(θ) ≈1-θ2/2,因此式变为:

      $$ {I_0}\theta ''\left( t \right) = {F_{\rm{m}}}\left( {{u_0}} \right)\left( {1 - {\theta ^2}} \right) + {W_{\rm{m}}}\theta - {M_{\rm{m}}} + \left( {{M_{\rm{m}}} - {M_{\rm{0}}}} \right){\left( {1 - \theta /{\theta _{\rm{m}}}} \right)^2}。 $$ (4)

      θ较小时,角加速度θ″(t)可用旋转角速度ω=θ′(t)表示为:θ′(t)=ω′(t)=ωdω/dθ,则式(4)变为:

      $$ {I_0}\omega \frac{{{\rm{d}}\omega }}{{{\rm{d}}\theta }} = {F_{\rm{m}}}\left( {{u_0}} \right) - {M_{\rm{0}}} + \left[ {{W_{\rm{m}}} - \frac{{2\left( {{M_{\rm{m}}} - {M_{\rm{0}}}} \right)}}{{{\theta _{\rm{m}}}}}} \right]\theta +\\ \left[ {\frac{{{M_{\rm{m}}} - {M_{\rm{0}}}}}{{{\theta ^2}_{\rm{m}}}} - {F_{\rm{m}}}\left( {{u_0}} \right)} \right]{\theta ^2}。 $$ (5)

      式(5)综合了树木旋转时的动能,同时角速度ω也是θ的函数,则有:

      $$ \frac{1}{2}{I_0}{\omega ^2} = \left[ {{F_{\rm{m}}}\left( {{u_0}} \right) - {M_{\rm{0}}}} \right]\theta + \left( {\frac{1}{2}{W_{\rm{m}}} - \frac{{{M_{\rm{m}}} - {M_{\rm{0}}}}}{{{\theta _{\rm{m}}}}}} \right){\theta ^2} + \\ \frac{1}{3}\left[ {\frac{{{M_{\rm{m}}} - {M_{\rm{0}}}}}{{\theta _{\rm{m}}^2}} - {F_{\rm{m}}}\left( {{u_0}} \right)} \right]{\theta ^3} + C。 $$ (6)

      由式(6)可知,当θ→0时,树木在风载作用下开始倾斜,此时ω=0,C=0。则式(6)可表示为:

      $$ \frac{1}{2}{I_0}{\omega ^2} = \\ \theta \left[ {\left( {{F_{\rm{m}}}\left( {{u_0}} \right) - {M_{\rm{0}}}} \right) + \left( {\frac{1}{2}{W_{\rm{m}}} - \frac{{{M_{\rm{m}}} - {M_{\rm{0}}}}}{{{\theta _{\rm{m}}}}}} \right)\theta + \frac{1}{3}\frac{{{M_{\rm{m}}} - {M_{\rm{0}}}}}{{\theta _{\rm{m}}^2}} - {F_{\rm{m}}}\left( {{u_0}} \right){\theta ^2}} \right]。 $$ (7)

      由以上分析可知,式(7)只有当θm<5°时才成立,此时θ可以用θm线性表示:

      $$ \theta = \varphi \theta 。 $$ (8)

      则式(7)可变为:

      $$ \frac{1}{2}{I_0}{\omega ^2} = \\ \phi {\theta _{\rm{m}}}\left[ {\left( {{F_{\rm{m}}}\left( {{u_0}} \right) - {M_{\rm{0}}}} \right) - {M_{\rm{0}}} - \left( {{M_{\rm{m}}} - {M_{\rm{0}}}} \right)\phi + \frac{1}{3}\left( {{M_{\rm{m}}} - {M_{\rm{0}}}} \right){\phi ^2}} \right]。 $$ (9)

      φ取正值的时候,式(9)右边如果也保持正值,那么角速度ω也将保持正值,这时随着θ的增大,树木将继续倾斜,根据式(9)可以得到如图 4的函数图像。从图 4中可知,当满足条件(10)时式(9)没有正根:

      $$ \begin{array}{l} {\left( {{M_{\rm{m}}} - {M_{\rm{0}}}} \right)^2} < \frac{4}{3}\left( {{F_{\rm{m}}}\left( {{u_0}} \right) - {M_{\rm{0}}}} \right)\left( {{M_{\rm{m}}} - {M_{\rm{0}}}} \right)\\ or{\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {F_{\rm{m}}}\left( {{u_0}} \right) > \frac{1}{4}\left( {3{M_{\rm{m}}} + {M_{\rm{0}}}} \right)。 \end{array} $$ (10)

      图  4  式(9)的函数图像[a表示在条件(10)满足的时候,b表示条件(10)不满足的情况]

      Figure 4.  Function diagram of the formula (a show the condition (10) is satisfied, b show the condition (10) is dissatisfied)

      u0足够大能满足不等式(10)时,树木风倒便会发生。

      针对图 1中的树木模型,将树干模拟成一个圆柱,其半径为r0,高度为h,则作用于单元dy上的荷载为:ρacDr0u02dy。其中:ρa表示空气密度,cD表示无量纲的阻力系数(与树干的半径有关),u0表示作用于树干的平均连续稳定风速。则风荷载相对于根部所产生的弯矩为:

      $$ {F_{\rm{m}}}\left( {{u_0}} \right) = \int {_0^h} {\rho _{\rm{a}}}{c_{\rm{D}}}{r_0}u_0^2y{\rm{d}}y = \frac{1}{2}{\rho _{\rm{a}}}{c_{\rm{D}}}{r_0}{h^2}u_0^2。 $$ (11)

      实际上,树木在风载作用下产生倾斜,树木重心位置会发生偏移,重力矩Wm将会对弯矩Fm(u0)产生一个附加的作用,这个附加作用可以通过对式(11)进行修正来体现。

      $$ {{\bar F}_{\rm{m}}}\left( {{u_0}} \right) = \frac{1}{2}{\rho _{\rm{a}}}{c_{\rm{D}}}{r_0}{h^2}u_0^2\left( {1 + \varepsilon } \right) = \frac{1}{2}{\rho _{\rm{a}}}{{c'}_{\rm{D}}}{r_0}{h^2}u_0^2。 $$ (12)

      式(12)中:ε=0.4ρgh3/(Er02),ρ表示树木的密度,E表示弹性模量。式(12)相当于式(11)中的阻力系数cD被提高了1+ε倍。

    • 基于2.1的理论分析,以西加云杉为例进行数值计算,选取相应参数为:h=19 m,r0=0.1 m,ρ=900 kg·m-3ρa=1.2 kg·m-3cD=0.6。可计算得到抵抗力矩:θm=2.5°=0.044 rad,M0=8 × 103 N·m,Mm=12.8 × 103 N·m。通过对根系强度的可变性进行分析,对上述大小的西加云杉树木,其根倒弯矩值可能在10×103 N·m到50 × 103 N·m,并且能给出图 1模型树木对根部的转动惯量:

      $$ {I_0} = \frac{1}{3}{\rm{\pi }}\rho r_0^2{h^3}。 $$ (13)

      对于连续稳定风速引起的风倒,满足条件式的最小风速为:

      $$ \frac{1}{2}{\rho _{\rm{a}}}{c_{\rm{D}}}{r_0}{h^2}u_0^2 = \frac{1}{4}\left( {3{M_{\rm{m}}} + {M_{\rm{0}}}} \right)。 $$ (14)

      可得u0=29.9 m·s-1,由于式(13)并未考虑根系和树枝惯性矩的影响,因此计算得到风倒发生时的风速可能是被低估了的风倒标准风速。在临界风速u0=29.9 m·s-1作用下,树木至少需要2 s的时间让倾斜达到θm,而需要4.5 s的时间让树木旋转角速度ω达到最小值。总之,树木如果受到连续稳定风速作用并产生弯矩Fm(u0),当满足条件式(10),那么树木将在连续的风速作用下发生倾覆。

      除西加云杉以外,还对稻类植物、悬铃木Platanus acerifolia也进行了相应的分析,基于各类参数计算得到稻类植物发生风倒的临界风速u0=11.6 m·s-1,悬铃木发生风倒的临界风速u0=22.8 m·s-1,而根据所查资料[18],3种林木实际发生风倒的数值分别:26.3, 8.8和19.1 m·s-1。注意到,理论计算值和实际数据比较吻合,但是计算值看起来要比实际数据大一些。究其原因是:实际发生风倒与林木个体参数、土壤条件、根系条件以及枝条、树叶与风荷载所形成的“流固耦合”效应有关。以典型的小麦Triticum aestivum来说,一般情况下都有2到6个分支,假设分支数目为n,如果每一个分支所承受的弯矩都有根系来支持,就降低了理论计算所得到的风倒临界风速。再比如“流固耦合”效应导致流体边界的变化,从而对流动产生反作用。不考虑“流固耦合”作用会使得计算结果是偏于安全的,也即“流固耦合”效应导致实际发生风倒的临界风速降低。根系在风倒过程中有至关重要的地位,长期的连续风荷载作用,会降低根系的支持力和抵抗弯矩的能力,以至于林木可能在较小的风速下发生风倒。因此,理论计算值与实验数据结果比较吻合,同时风的动态作用可以得到合理的解释。

    • 实际上树木更容易在强大的阵风作用下发生风倒,假设其作用时间为τ,在风速为u0连续风载作用下根部产生倾覆弯矩Fm(u0),则t=0时刻阵风所引起的冲量矩为最大值,树木的角动量为[I0ω]0-0+=(Mg-M0)τ,其中Mgτ表示阵风所产生的冲击荷载引起的冲量矩,M0τ表示阵风引起根部的抵抗力矩。则树木的角速度为:

      $$ {\omega _0} = \left( {{M_{\rm{g}}} - {M_{\rm{0}}}} \right)\tau /{I_0}。 $$ (15)

      式(7)变为:

      $$ \frac{1}{2}{I_0}{\omega ^2} = \\ \left[ {{F_{\rm{m}}}\left( {{u_0}} \right) - {M_{\rm{0}}}} \right]\theta + \left( {\frac{1}{2}{W_{\rm{m}}} - \frac{{{M_{\rm{m}}} - {M_{\rm{0}}}}}{{{\theta _{\rm{m}}}}}} \right){\theta ^2} +\\ \frac{1}{3}\left[ {\frac{{{M_{\rm{m}}} - {M_{\rm{0}}}}}{{\theta _{\rm{m}}^2}} - {F_{\rm{m}}}\left( {{u_0}} \right)} \right]{\theta ^3} + \frac{1}{2}{I_0}\omega _0^2。 $$ (16)

      t=0+时,积分常数C=I0ω02/2,阵风的作用,导致产生了力矩Mg,将Mg表达成连续稳定风速作用下弯矩Fm(u0)的函数:

      $$ {M_{\rm{g}}} = {F_{\rm{m}}}\left( {{u_0}} \right){G^2}_{\rm{v}}。 $$ (17)

      式(17)中,${G_{\rm{v}}} = \frac{{\hat u}}{{\bar u}} $为阵风系数,${\hat u} $觠指阵风风速极值,u指阵风风速的平均值。让θ=φθm,由式(16)可得:

      $$ \frac{1}{2}{I_0}{\omega ^2} = \\ \left[ {\left( {{F_{\rm{m}}}\left( {{u_0}} \right) - {M_{\rm{0}}}} \right) - \left( {{M_{\rm{m}}} - {M_{\rm{0}}}} \right)\phi + \frac{1}{3}\left( {{M_{\rm{m}}} - {M_{\rm{0}}}} \right){\phi ^2}} \right]\phi {\theta _{\rm{m}}} +\\ \frac{1}{{2{I_0}}}{\left[ {{F_{\rm{m}}}\left( {{u_0}} \right)G_{\rm{v}}^2 - {M_{\rm{0}}}} \right]^2}{\tau ^2}。 $$ (18)

      显然,式(18)右边必须保持正值,树木才会继续发生倾斜,在其他所有参数已知的情况下,式(18)是关于φ的一个三次函数,假设其梯度为0,则退化为关于φ的二次函数,在$\varphi = {\varphi _{\rm{r}}} = 1 + {\left( {\frac{{{M_{\rm{m}}}-{F_{\rm{m}}}\left( {{u_0}} \right)}}{{{M_{\rm{m}}}-{M_{\rm{0}}}}}} \right)^{\frac{1}{2}}} $处可以得到根的最大值。要保证式(18)右边是正值,则必须保证:

      $$ G_{\rm{v}}^2\tau > \frac{{{M_{\rm{0}}}\tau }}{{{F_{\rm{m}}}\left( {{u_0}} \right)}} + {\left[ {\frac{{2{I_0}\left( {{M_{\rm{m}}} - {M_{\rm{0}}}} \right){\theta _{\rm{m}}}}}{{3F_{_{\rm{m}}}^2\left( {{u_0}} \right)}}\left( {2{\varphi _r} - 3} \right){\varphi _r}^2} \right]^{\frac{1}{2}}}。 $$ (19)

      成立。由式可以得到相应的判断条件:

      $$ {F_{\rm{m}}}\left( {{u_0}} \right) < \frac{1}{4}\left( {3{M_{\rm{m}}} + {M_{\rm{0}}}} \right)。 $$ (20)

      可以看出:式(20)是对式(10)的很好补充。其物理含义是:式(10)指出当风速足够大的时候,在连续稳定风速的作用下,风倒便会发生。式(20)指出当连续风速不够大,以至于不能够引起风倒的情况下,连续的风载将和阵风共同作用下引起风倒[阵风系数GV足够大,并且作用时间τ足够长,能满足式(19)]。式(19)取决于如下的无量纲量:Fm=Fm(u0)/M0Mm=Mm/M0I0=I0θm/M0τ2。则条件表示的破坏阵风系数必须满足:

      $$ G_{\rm{v}}^2 > \frac{1}{{\overline {{F_{\rm{m}}}} }}\left\{ {1 + {{\left[ {\frac{2}{3}\overline {{I_0}} \left( {\overline {{M_{\rm{m}}}} - 1} \right)\left( {2{\phi _r} - 3} \right){\varphi _r}^2} \right]}^{\frac{1}{2}}}} \right\},\\ {\phi _r} = 1 + {\left( {\frac{{\overline {{M_{\rm{m}}}} - \overline {{F_{\rm{m}}}} }}{{\overline {{M_{\rm{m}}}} - 1}}} \right)^{\frac{1}{2}}}。 $$ (21)
    • 由无量纲化的风倒破坏式(21)可知,风倒条件取决于所定义的无量纲量,且都是风速u0和阵风作用时间τ的函数。以西加云杉为例进行数值计算,所选参数与2.2中的相同,则无量纲参数变为:

      $$ \overline {{F_{\rm{m}}}} = 1.62{u_0}^2{10^{-3}}, \overline {{M_{\rm{m}}}} = 1.6, \overline {{I_{\rm{0}}}} = 0.352{\tau ^{-2}}, {\phi _r} = 1 + {\left( {2.67-2.7{u_0}^2{{10}^{ - 3}}} \right)^{\frac{1}{2}}} $$

      当阵风风速最大值$ {\hat u}$觠和作用时间τ子满足式(21),风倒便会发生,可得到风倒的风速临界值:

      $$ \hat u = 24.81\left\{ {1 + 0.376\frac{{{\phi _r}}}{\tau }{{\left( {2{\phi _r} - 3} \right)}^{\frac{1}{2}}}} \right\}。 $$ (22)

      图 5表示出阵风风速${\hat u} $和连续风速u0(可以由阵风风速的平均值u代替)在不同的时间段τ内的关系。从图中可以看到,当u0接近于临界风速u0=29.9 m·s-1的时候,所有曲线都将终止。如果将1 h风速的作用时间按照τf时间段来划分,可以得到这段时间内发生的最大风速。

      图  5  引起风倒的阵风风速和连续风速的函数关系(在τ=0.2,0.5,1,2,5,10s的情况下)

      Figure 5.  Function relationship between gust wind and continuous wind about windthrow (in the case of τ=0.2, 0.5, 1, 2, 5, 10s)

      $$ \frac{{\hat u}}{{\bar u}} = 1 + 0.42\left( {\frac{\sigma }{{\bar u}}} \right)\ln \left( {\frac{{3600}}{{{\tau _f}}}} \right)。 $$ (23)

      式(23)中,σ/u表示紊流强度,对于独立树木,σ/u一般在0.2左右;而对于森林,σ/u会更高一些。依据式(23),给出不同的τ值,可以得到如图 6所示的一族直线,这些直线穿过图 5中的曲线,相应产生交点,将交点连在一起,如图 7所示,得到风倒发生的风速标准线。在这个直线下方,不可能出现阵风风速${\hat u} $。由图 7可知:风倒可能在较小风速17.7 m·s-1就可能发生。而对于阵风来说,在风速为27.5 m·s-1,作用时间5 s的情况下风倒就可能发生。尽管5 s时间给树木以冲击荷载看起来比较短,但于树木来说,从受风载作用开始倾斜到其抵抗弯矩最大值位置,大约需要2 s。所以,可以寻找作用时间更短的阵风,对于阵风风速为35.1 m·s-1,作用时间τ=1 s叠加在风速为21.0 m·s-1的风速上就会引起风倒。

      图  6  紊流强度为0.2的时候破坏的平均, 风速图

      Figure 6.  Diagram of average wind speed when turbulence intensity is 0.2

      图  7  紊流强度为0.2,τ=1,2,5s破坏的平均风速图

      Figure 7.  Average wind speed diagram of damage when turbulence intensity is 0.2, τ=1, 2, 5s

    • 林木的破坏模式有2种,一种是破坏值超过林木材料的允许值而发生的折断现象,另外一种破坏值超过根系的支持力而发生的连根拔起现象。本研究综合考虑在稳定风速和阵风作用下树木的倾覆问题。将树木风倒用一个简单的模型来模拟根系,从动力学角度建立了风倒的动态方程,得到在稳定风速下树木发生风倒的条件。同时将阵风作用模拟为一个冲量荷载作用在树木上部,建立阵风作用下树木风倒的条件。以西加云杉为例,计算了树木在稳定风速和阵风作用下破坏的标准值,将阵风风速的作用时间与平均风速作用时间比较,得到树木发生风倒的极限风速。但是,实际林木风倒临界风速值低于文章所计算出的理论值,这是由于林木树冠形状、个体参数差异、枝条、树叶与风荷载之间的“流固耦合”效应以及复杂的非线性动力学行为都导致实际发生风倒的风速值偏低。这也为今后的研究提供了方向和思路。总之,文章从理论方面的研究结果对森林的抗风防护以及木材的生产实践有很重要的现实指导意义。

参考文献 (18)

目录

    /

    返回文章
    返回