基于双向流固耦合的林木风致响应及防风效果研究.pdf

摘要 摘要 本文以山东省东营市军马五分场刺槐人工防护林作为实验样地，在刺槐林场进行 风速、风向的实测实验并确定标准木。基于分形理论的L-系统，根据标准木的树高、 枝下高、胸径和冠幅等每木检尺数据，利用Matlab 和Ansys 建立刺槐的三维模型。 设置刺槐的9个弹性常数并命名为tree，作为刺槐主干和分支的材料属性。利用线性 滤波法中的AR 法，基于Matlab 编程，以实测风速为初始风速，以Davenport 谱作为 目标谱，采用Shiotami 的空间相关性表达式，获得脉动风速时程，并添加到脉动风场 的入口处。在Ansysworkbench 平台实现刺槐与脉动风场的双向流固耦合数值模拟。 研究结果表明，防护林风速的实测值与模拟值具有较高的一致性，拥有标准k-ε 湍流模型的数值模拟可以为林木抗风提供良好的预测。刺槐会对风产生阻碍作用，在 刺槐的背风面形成一个低风速的椭圆形的庇护区。刺槐模型对风产生绕流，从而在刺 槐的周围产生两个风速增大的区域。刺槐受风会进行摆动，其中刺槐的主干与一级分 枝的位移具有相似性，主干顶点的位移小于一级分枝顶点的位移。刺槐在脉动风场的 作用下将会产生风致响应，其中YZ 方向的剪应力最大，Z 方向的压应力最大。刺槐 Z 方向的拉应力最大，Y 方向与X 方向的拉应力大小相近。树高与刺槐的防风效果成 正相关，增大树高可以显著增加刺槐的防风距离，疏透度与刺槐的防风效果呈负相关， 减小刺槐的疏透度，可以降低刺槐背风处的风速，刺槐的防风效果越好。随着刺槐行 数增多，产生的低风速庇护区的范围越大，对风速的减弱效果越明显。单行刺槐对于 风的减弱效率是最高的，虽然增加行数可以有效地增大刺槐的防风距离，但是也会使 其防风效率极大地降低。缩减刺槐的行间距可以获得更好的防风效果，3米行间距的 两行刺槐可以提供一个整体性的低风速的庇护区。当行间距增大到15米以上时，不 需要将两行刺槐当作一个整体，可以近似的看成两行分离的刺槐。线性排列的分布方 式可以获得较宽的防护区域，交错排列的分布方式防护区域则较长，交错排列的刺槐 林的防风效果优于矩形排列的防风效果。 关键词：分形理论，脉动风场，林木抗风，数值模拟，防护林优化 I 目录 目 录 摘要I ABSTRACTII 目 录IV 1 引言1 1.1 研究背景与意义1 1.2 国内外研究现状2 1.2.1 林木抗风的研究2 1.2.2 树木力学模型的研究4 1.2.3 风场的研究5 1.3 研究内容与创新点7 1.3.1 研究内容7 1.3.2 创新点8 1.4 技术路线8 2 材料与方法9 2.1 研究区概况9 2.2 风场实测的实验仪器与方法9 2.2.1 实验仪器9 2.2.2 实验方法 11 2.3 流固耦合的基础理论 12 2.3.1 流固耦合的基本概念12 2.3.2 流固耦合的分类12 2.4 刺槐模型建立的基础理论 13 2.4.1 分形理论的基本概念13 2.4.2 分形理论的应用14 2.5 脉动风场实现的基础理论 14 2.5.1 风的基本概念14 2.5.2 风速与风压15 2.5.3 基本风速与基本风压15 2.5.4 平均风与脉动风15 2.5.5 脉动风速功率谱16 2.5.6 脉动风的空间相关性17 2.5.7 线性滤波法的计算过程18 IV 目录 3 双向流固耦合数值模拟19 3.1 刺槐三维力学模型的建立 19 3.1.1 刺槐标准木的选择19 3.1.2 利用Matlab 建立刺槐的三维模型 19 3.1.3 刺槐力学模型的建立20 3.2 脉动风场的建立22 3.2