图1 植物模型参数和计算域边界条件（A），以及柔性方柱阵列流固耦合（FSI，B）示意图

Fig.1 Schematic diagram of plant model parameters and boundary conditions of computational domain （A）， and numerical simulation scheme for FSI using a cylinder array as the plant （B）

图2 不同时刻的瞬态速度（A）、瞬态压力（B）与涡结构（C）云图

Fig.2 Transient velocity （A）， transient pressure （B）， and vortex structure （C）

图3 不同模型植物内部瞬态速度随时间变化曲线

Fig.3 Transient velocity variation curves within the canopy layer over time， with numerical calculation results incorporating flexible plant sway compared to previous experimental data

图4 VPM（A）与FSI（B）模拟下柔性植物速度侧视云图

Fig.4 Velocity contour of flexible plants by VPM （A） and FSI （B）

图5 3种网格对应的不同流向位置的风廓线

Fig.5 Three mesh refinement forms and wind profiles at different flow positions

图6 不同工况的柔韧性植物的阻力系数随时间变化曲线

Fig.6 Drag coefficient variation over time for flexible plants under different operating conditions

图7 不同工况对应的瞬态（点线）与平均（虚线）挡风系数对比

Fig.7 Comparison of transient （solid line） and time-averaged （dashed line） wind blockage coefficient

图8 不同植物柔韧性对应的平均速度（A）、瞬态速度（B）与瞬态压力（C）剖面云图

Fig.8 Time-average velocity （A）， transient velocity （B）， and transient pressure （C） contour plots corresponding to different plant flexibility levels

图9 植物冠层流动分区与监测点示意图

Fig.9 Schematic diagram of flow zoning and monitoring points within the plant canopy

图10 前端接触区的速度脉动与湍动能（A，点线为脉动速度均方，虚线为湍动能TKE）、速度功率谱密度（B）、瞬态速度（C，D，实线为瞬态流向速度，虚线为时间平均流速）

Fig.10 Fluctuating velocity （A， dot solid lines for RMS values and dashed lines for turbulent kinetic energy）， spectrogram of velocity energy density （B）， transient velocity （C， D， solid lines for instantaneous and dashed lines for time-averaged velocity） in the front contact zone

图11 冠层内部脉动速度与湍动能曲线（A）、速度功率谱密度（B）、瞬态流速（C， D）

Fig.11 Fluctuating velocity and turbulent kinetic energy （A）， velocity energy density spectrum （B） and transient velocity curves （C， D） inside of canopy

图12 近尾流区脉动速度与湍动能（A）、速度功率谱密度（B）与瞬态流速（C）

Fig.12 Fluctuating velocity， mean turbulent kinetic energy （A）， velocity power spectral density （B）， and instantaneous flow velocity （C） in near wake region

图13 远尾流区脉动速度与湍动能（A）、速度功率谱密度（B）与瞬态流速（C）

Fig.13 Fluctuating velocity， mean turbulent kinetic energy （A）， velocity power spectral density （B）， and instantaneous flow velocity （C） in far wake region

图14 工况4与工况5对应的瞬态速度与时间平均流速在近尾流区（A）与远尾流区（B）的对比

Fig.14 Comparison of time-frequency curves of flow velocity between operating case 4 and 5

图15 不同柔韧性的植物在不同流向位置处的风廓线，以及刚性孔隙植物模型与柔性植物模型对风廓线计算的差异

Fig.15 Wind profiles at different streamwise locations， and the differences in wind profile between rigid porous plant models and flexible plant models

图16 不同时刻下的柔性植物附近瞬态地表剪切应力分布注：白色的虚线为植物所在的区域，红色虚线为数据监测的位置

Fig.16 The transient distribution of wall shear stress near flexible plants， the white dashed line represents the area occupied by the plants

图17 不同工况对应的平均地表剪切应力沿流向距离变化曲线（A），以及不同时刻下中轴线位置处瞬态地表剪切应力沿流向距离变化曲线（B）

Fig.17 Variation of average wall shear stress along the flow direction （A）， and variation curves of transient wall shear stress along the flow direction （B）