基于参数化设计的WORKBENCH CFX的动态流场计算方法
本文摘要（由AI生成）：
本文介绍了ANSYS CFX软件在CFD分析中的优势，包括稳健的算法、智能化选择参数、内嵌式开发语言等。同时，文章还讨论了CFX在动网格处理方面的局限性，并提出了一种基于WORKBENCH的解决方案，通过参数化驱动模块实现相关参数的几何驱动，从而解决了ANSYS CFX动态流场计算中几何面的识别问题。这种方法为动态流场计算提供了一种新的思路，并具有一定的借鉴意义。
0.前言
ANSYS CFX作为是全球第一个通过ISO9001质量认证的大型商业CFD软件。其基于有限元的有限体积法,在保证了有限体积法的守恒特性的基础上,吸收了有限元法的数值精确性,发展和使用全隐式多网格耦合求解技术的商业化软件,这种革命性的求解技术克服了FLUENT等软件中“假设压力项－求解－修正压力项”的反复迭代过程,而同时求解动量方程和连续性方程;因而可以看出,CFX采用了一种更加稳健的算法,对于一般的CFD分析问题,不需要进行参数的过多设置,而软件自身可智能化选择,因而对应用者的使用经验要求较低,大大降低了其在工程上推广的难度,这是FLUENT软件无法拥有的优越性[2-3]。
相比于FLUENT采用外部程序植入方法(UDF),CFX所具有的内嵌式开发语言如CFX命令语言(CFX Command Language,CCL)和CFX表达式语言(CFX Expression Language,CEL),不需要专业软件连接接口,且编程语法简单,实用性强。
由于CFX自身包含电磁场计算模型,因而在进行相关耦合场计算时,无需添加额外程序。
然而,由于自身软件的问题,CFX的动网格功能,只能够直接处理具有网格变形的小变形问题的研究,而不能像FLUENT那样实现网格较高级别的变形。针对具有大变形的CFX的动网格技术,目前存在几种网格重构方法,但由于编程、环境变量设置等种种原因增加了应用难度,其普遍应用还无法展开。
其对网格的变形只是基于简单的压缩和拉伸过程。目前,解决该问题的一种方法,是将一系列具有的大变形分解为只有变形的小变形,并将上一步迭代的初始值作为下一步计算的初值,这里具有微小差别的几何体,需要在外部CAD文件中处理后逐个输入到CAE软件中。
但是,由于不同软件间对几何面识别的差异,在CAD软件中绘出的几何模型,无法在CAE软件中被完全识别,因而可能造成相关面在边界定义时出现失败。具体来说,上一步定义好的相关面,在输入下一步发生微小变形后的几何体后,由于CAE软件自身无法实现相关面的自动识别,经常在网格划分时还需要重新定义相关面,因而导致计算流程程序化无法实施,而不得不采用手动方式完成,这将大大降低计算效率。
针对上述ANSYS CFX动态流场计算中存在的这一潜在问题,本文提出了一种解决方法。即将各个时间步骤下的几何模型导入,变为只导入同一几何模型,然后利用WORKBENCH中DM模块的几何处理能力,通过将涉及变动的相关位移通过parameter set模块实现参数化。这里,实施动态流场计算,只需导入最初始的几何体而不必输入批量输入,因而有效避免了将具有微小差别的批量几何体导入CAE软件时存在的相关面的识别问题；整个计算过程的实现是通过scripting程序化完成。该方法有效地解决了ANSYS CFX动态流场计算中几何面的识别问题,因而避免了上述潜在的流场计算失败；通过参数化驱动模块实现相关参数的几何驱动,来实现动态流场计算,也为动态流场计算提供了一种新的思路。另外,对于一般意义下的CAE优化分析问题,特别是模型来自实际工程图,重新绘制CAE分析所需模型十分繁琐时,该方法的解决思路也具有一定的借鉴意义。
因而,在workbench平台下,各模块的连接图如下图1所示。
图1 workbench连接流程图
1. 参数化计算原理
图2 计算程序流程图
本程序充分利用了workbench中DM模块的几何处理功能,通过将其中几何处理涉及到的变量参数化,即引入parameter set模块,从而在workbench软件内部便可实现模型的几何改变,故而解决了CAD和CAE软件不兼容问题对动态流场计算造成的可能问题。其计算流程图如图2所示。
在该研究方法中,首先应用了中DM的几何处理功能。DM作为workbench的专业几何处理软件,虽然没有专业CAD软件强大的绘图和几何处理能力,但是其在三维几何生成、几何修改和实体操作方面具有一定的功能,可以实现几何体的镜像、移动、复 制和布尔运算等,本文主要应用了其阵列(pattern)命令。在该命令中,通过对运动相关部件进行线性阵列从而达到移动的目的。而且在该命令中,可以将移动的距离参数化并传输parameter set中。对于初始输入的几何体,采用几何抑制(suppress)功能可使其不会出现在后续的流场计算流程中。
在workbench平台中,parameter set模块本身是主要用于模型参数优化设计,它能够将其他模块中相关变量参数化。在本文中,笔者将DM中相关运动部件的运动距离参数化,并通过parameter set模块实现对该参数的改变,具体操作过程则是通过scripting程序来实现的。
Workbench的scripting是一种采用PYTHON语言编写的程序,其生成文件后缀为.wbjn。PYTHON语言是一种开放源码软件,其代表简单主义思想的语言[15]. Workbench 通过采用这种程序,记录操作者通过GUI界面进行的相关操作,并进行简单编辑,实现命令流的功能。
可将,笔者参数化动态流程计算方法的实现,主要集成以上三个模块,即DM几何处理、parameter set参数化驱动及 scripting编程。
2. 一般步骤
本文中,基于参数化的动态流场计算,在其实施过程中,主要有几个步骤。
(1)数据预处理
需要根据原始文件中给出的等时间间隔的位移求出等位移的时间间隔,同时求出不同时间下的时间步长。本方法采用Ansys cfx进行动态流场计算的原理,就是基于每次计算的运动的距离相等,这即保证了在参数不变的情况下,网格质量不会出现问题,也实现了前一步计算的结果能够作为后步计算的初场使用。
(2)UG处理及文件导出
对于来自实际工程的3D模型,通过对其进行必要的简化处理,模型整合,可以将模型分为两部分,即运动几何部分和静止几何部分。根据涉及动网格的动态流场计算的原理,在UG中采用几何分域的方法,以提高网格划分效率。
因此,在最终的几何模型中,应将整个几何模型至少分为三个部分,即:
a.外部网格不变的静止域;
b.内部包含运动件的静止域;
c.运动部件;
通过后两者相减,可以得到流场计算中包含动边界的流体区域。
在将生成的几何导出时,首先将模型a导出至一个几何,而将模型b和c导出至一个几何中。
(3) Ansys CFX平台设置
本部分的主要工作是,根据上文所述的参数化计算原理,搭建诸如图1所示的计算流程图,并依据流体计算的基本步骤,实施几何导入DM中处理、meshing网格划分、求解设置、求解和结果输出等工作。
(4) workbench流程控制及计算
采用WORKBENCH的脚本语言功能,录制包括修改parameter set中参数在内的相关命令；采用PYTHON语言,对生成的scripting程序进行修改,加入程序循环实施的功能等；运行修改好后的程序,即可由计算机完成整个流场计算过程。
3.4 workbench流程控制及计算
采用WORKBENCH的命令流功能,录制包括修改parameter set中参数在内的相关命令；采用PYTHON语言,对生成的scripting程序进行修改,加入程序循环实施的功能等；运行修改好后的程序,即可由计算机完成整个流场计算过程。
