论文部分内容阅读
薄膜屋盖结构与风场之间的流固耦合作用是影响结构风振响应的重要因素。数值模拟是流固耦合研究的重要手段,准确模拟流体域的时变特性是掌握流固耦合效应产生机理和发展规律的前提。风与薄膜结构流固耦合运动中的流体域是含动边界的钝体绕流近地风场,在采用计算流体动力学(Computational Fluid Dynamics,CFD)技术对其进行大涡模拟(Large-eddy Simulation, LES)数值近似计算时,仍存在一些亟待解决的问题。建立适于大型工程应用的亚格子模型是大涡模拟方法理论研究中的热点问题之一。在求解近地风场大涡模拟控制方程时,投影法是一种高效的压力速度解耦算法。目前大多数投影格式对压力的计算在时间方向上比速度的计算至少低一阶精度,而压力场的时间精度将直接影响流固耦合效应的计算精度。因此,投影法的压力时问精度还有待提高。此外,动边界绕流模拟适宜在任意拉格朗日欧拉(Arbitrary Lagrangian Eulerian, ALE)动态网格框架下进行。模拟过程中网格更新频繁,其几何品质和更新前后的一致性是影响模拟精度的重要因素。目前大部分动网格更新算法仅关注更新后网格的几何品质,对更新前后网格一致性的测试不够;而动网格更新算法的计算耗时也是需要考虑的一个课题。针对上述问题,本文对亚格子模型、投影法及其在ALE动态网格下的应用以及动网格更新技术等内容进行了理论和数值研究,开展的研究工作和提出的解决方案主要包括:1、亚格子模型的测试和对比。在阐述近地风场的大涡模拟无量纲控制方程以及多种亚格子模型的构造过程后,在笛卡尔静态网格下编制了方程求解程序。以该程序为计算平台,对多种亚格子模型的适用性、准确性和计算效率进行了测试。结果表明,钝体绕流风场的大涡模拟工程应用中可以采用阻尼修正的标准S模型。2、动态网格下数值求解方法的建立。以ALE描述为基础,在动态网格下建立了大涡模拟控制方程(下文简称ALE-LES方程),并对方程中的网格运动参数和流动参数进行解耦,建立了交错求解方程的数值计算方法。期间,为提高数值方法的时间计算精度,构造了压力和速度能达到同一高阶时间精度的全精度连续投影方法,并将其与求解网格运动的二阶预测校正格式相结合,得到了ALE-LES方程的全二阶精度投影格式。将建立的数值求解方法进行坐标转换,在贴体坐标系下编制了ALE-LES方程求解程序的核心模块。3、动网格更新算法的对比研究。归纳整理出网格质量评价模型,并据此在两种畸变网格算例中,为适用于结构网格的多种动网格更新算法对网格几何品质的改善力度进行了对比;同时,对不同算法在更新前后网格的一致性和计算效率方面进行了测试,指出九点网格重构法综合性能最为优异。将九点网格重构法扩展到三维研究领域并更名为距离加权法,据此编制了网格动态更新子程序模块。4、数值模拟程序的开发和验证。将亚格子模型子程序模块和网格动态更新子程序模块与求解ALE-LES方程的核心程序相结合,开发了适用于模拟含运动边界钝体绕流风场的CFD计算程序。分别以二维方腔驱动流和Taylor涡列问题验证了程序的稳定性和时间计算精度。5、振动屋盖绕流特性的数值分析。以下部封闭式大跨度平屋盖结构为研究对象,应用本文开发程序,对屋盖以二阶谐波模态振动时的结构绕流场进行了数值模拟,通过与相同尺寸刚性屋盖结构绕流模拟结果的对比分析,研究了屋盖的振动对模型表面风压系数分布和绕流场流动结构的影响。结果表明,屋盖振动是影响建筑结构表面平均、脉动风压系数的重要因素;屋盖振动导致绕流场旋涡结构复杂化,特征湍流度增加;屋盖的振动对绕流场瞬时风压分布的影响明显,二阶谐波模态的振动形式容易引发脉动风与屋盖的共振,不利于屋盖的抗风稳定。自主开发程序的成功应用表明了该程序可作为薄膜结构流固耦合效应数值模拟计算平台的流体域求解器。