稳态低速风压分布计算.pdf免费

第十届全围现代结构T程学术研讨会 稳态低速风压分布计算 韩向科1钱若军1袁行飞2王春江3林智斌4 (1同济大学土木工程学院，上海，200092；2浙江大学土木T程学院，杭州，310027 3上海交通大学空问结构研究中心。上海． 200030：4威斯康星州大学美国) 提要： 利用计算流体动力学求解结构工程中建筑物风压分布。以虚拟时间求解技术求解稳态问题，推导了低速不压流体 基于特征线离散的有限元方程。以典型的矩形建筑物风场模拟作为算例，给出了风压沿建筑物高度的变化规律， 结果与实际情况较为符合。 关键词： 不可压，有限元，风场模拟 一、稳态低速风作用流体控制方程 对于工程中常见的风而言，其风速都V100m／s，故密度变化较小，所以都可以将其简化为不可压流 动，因此研究不可压风流动对于实际工程具有重要的意义。 对于不可压流动，能量方程不再同连续和动量方程耦合，可以单独求解。因此对于不可压流动问题主 要是求解连续方程和动量方程。 当建筑环境相对简单且只考虑风对结构的静态效应时，可将问题局限为稳态层流粘性不可压流体的流 动问题，此时NS方程组 。 娶+要+孥+Q：D af 毗 毗 (式中符号意义同前)便退化为稳态层流粘性不可压方程，为此略去方程中所有关于时间的项，并令密度为 常数，经过整理最终可得： 连续方程 塑=o 啦 式中， 匕i-pv,，其中，V为流体速度，p为流体密度。 动量方程 (2) 孤；、q∥…孤；毗广川。 一瓦a(嘶)鲁音嗍lf=o 式中，K为对流速度；g。为重力加速度；互，为偏应力分量，对于不可压流体： 基金项目：国家自然科学基金资助，浙江省自然科学基金(R1080363) 工业建筑2010增刊 935 第十届全国现代结构工程学术研讨会 例睁鼍] (3) 其中，∥为动力粘度系数；当f=J时磊=l，i≠，时磊=o。 上述方程便为稳态层流粘性不可压流体的控制方程。实际上稳态流动问题可以看作为瞬态流动问题在 某一时刻点上的稳态问题，即瞬态流动问题可以分解为一系列瞬时的稳态流动问题，因此研究稳态流动问 题是首要的。同时，如前可知湍流问题只需在原有的NS方程组添加相应的物理量平均项将问题简化为修 正的NS方程组获得最终求解，因此研究稳态层流粘性不可压流体在工程领域内具有普遍意义。 二、基于虚拟时间稳态问题求解方法 正如前所述，稳态流动问题可以看作为瞬态流动问题在某一时刻点上的稳态流动。这一方法可以通过 引入虚拟时间f来实现。这一方法首先由Jameson川针对瞬态问题提出，利用双时间步法(dual—timestepping method)在计算效率方面具有强大的优势，大量的研究工作充分证实了这一点。 对于方程 警=肿) (4) 若采用显式差分格式则势必时间步要严格控制，若采用隐式差分格式虽然是无条件稳定的，但需要求 解庞大的矩阵方程，特别是求解三维问题，也是相当耗时的。针对上述离散格式的缺陷，Jameson提出了 一种隐式时间离散方法，大大提高了瞬态计算的效