复杂帧式起落架的三维风压绕流特性.docxVIP

复杂帧式起落架的三维风压绕流特性 cd-rom照明系统工作在远离高日照明的户外环境中。该钢结架具有较大的风速和较差的透风性。风载荷直接影响机架刚度、强度以及驱动系统的设计, 且大风作用导致的聚光器镜面破裂是较常见的损坏形式, 同时随碟式系统功率级别的增加其迎风面积将更大, 风载荷已成为机架设计的主要控制载荷。然而碟式光热发电技术在国内仍处于起步阶段, 对机架风载荷的研究文献报道甚少。虽机架承风主体与抛物雷达天线为同体型结构, 但由于抛物焦半径及结构存在的差异必然会引起风载荷数值及分布的不同, 已有抛物天线风洞试验仅供参考, 所以开展碟式钢构机架的风载荷研究尤显迫切。 实际工程中结构风载荷的研究方法主要有现场实测、缩尺模型风洞试验、数值风洞模拟等, 其中前两者是公认的能较准确获得结构风载荷的方法, 但都存在着试验经费大, 周期长等问题, 并且现场实测无法在结构建成前进行测试, 还受到实测场地因素的制约。同时, 碟式机架空间尺寸大且带有复杂桁架结构, 风洞试验时其缩尺模型制作也存在一定困难。然而, 数值风洞可进行结构全尺模型风载荷和绕流特性的模拟, 且有成本低、周期短等优点, 随着计算机软硬件的发展和计算流体力学技术的完善, 数值风洞模拟已成为一种有效的结构风载荷研究手段。Naeeni等采用风洞试验与数值模拟的方法对250 k W功率的太阳能发电站槽式聚光镜进行研究, 得到不同工况下聚光镜表面的风载荷。郑德乾等采用Realizable k-ε湍流模型对复杂世博轴膜面结构平均风压进行模拟, 将结果与风洞试验对比表明了模拟的有效性, 并进一步通过绕流特性探讨结构载荷的产生机理。刘若斐等也选用上述湍流模型对大型冷却塔风载荷进行模拟, 模拟结果与载荷规范比较接近。许多学者的研究均表明, 相比风洞试验方法, 数值风洞模拟可以对结构风载荷进行全方位多层次的分析, 并能给出满足工程应用精度的载荷数据, 这对结构抗风预研有着重要的意义。 文中将在Fluent 6.3软件平台, 选用基于Reynolds时均的Realizable k-ε湍流模型开展碟式机架多工况的风载荷研究, 获得不同工况的总风载荷、风压分布及绕流特性等结果。分析风载荷随不同工况的变化规律, 并将载荷曲线同抛物雷达天线风洞试验结果进行定性的对比分析, 验证模拟结果的合理性。探讨了聚光器镜面风载荷占机架总风载荷的比值, 给出了典型工况聚光器体型系数分布及各分区载荷极值, 通过分析流场绕流特性揭示了载荷产生的机制。所作分析能够为碟式光热太阳能整机的抗风设计提供参考。 1 数值场的基本理论 1.1 基于湍流模型的求解 工程中钝体结构绕流问题一般都是复杂的湍流运动, 流态表现为三维非定常且伴有强烈脉动的漩涡, 目前还不具备直接模拟计算的条件。然而, 有学者通过大量试验研究表明, 湍流运动的物理参数如速度、压力等的统计平均值存在着确定性规律且可重复再现, 于是工程中广泛采用了统计平均的方法来研究复杂的湍流运动。依据雷诺的观点将描述湍流物理量中的湍流瞬时量分解成平均和脉动两部分, 并且将流体基本控制方程 (N-S方程) 逐项平均就可得到基于雷诺应力的时均N-S方程: 由式 (1) 、 (2) 构成的方程组不封闭, 必须引入湍流模型才能对方程求解, 目前的方法主要有雷诺应力模型和湍流黏性系数法, 由于雷诺应力法求解的计算量非常大, 文中选择了工程应用较为广泛的后者。黏性系数法中标准k-ε模型能在一定程度上反映湍流的特征, 从而在工程中有着一定的应用, 但标准k-ε模型对雷诺应力的模拟是采用推广的Boussinsq各向同性的涡黏性假设, 导致该模型可能会出现不正确的计算雷诺正应力。有研究就表明, 标准k-ε模型应用于时均应变率特别大的情况时, 有可能导致负的正应力。鉴于标准k-ε模型对于一些各向异性较强的流动如有分离、漩涡的流动将产生较大的误差, 不能恰当的描述流动的发展, 特别是二次流动问题。为了使流动符合湍流运动的规律, 有学者提出了Realizable k-ε湍流模型, 其主要改进在于将湍流黏度计算式中的系数Cμ看成与应变率相联系的变量。杨伟等的研究也表明Realizable k-ε模型较标准k-ε模型的模拟值更接近于风洞试验值。因此文中选用Realizable k-ε湍流模型, 该模型关于湍动能k和湍流耗散率ε的输送方程及其参数定义详见文献。 1.2 近壁面湍流状态的非平衡壁面函数法处理u 近壁面的湍流流动受到分子黏性的显著影响, 这种情况的处理很大程度上决定了整个数值模拟结果的准确性。然而, Realizable k-ε模型一般只适用于湍流核心区域, 针对近壁面低Re数的流动需要进行特殊处理-壁面函数法, 文中采用非平衡壁面函数法来处理近壁面的湍流状态。大量的试验表明, 流场的近壁面区