第 19 卷 第 2 期
2010 年 6 月
计 算 机 辅 助 工 程
Computer Aided Engineering
Vol. 19 No. 2
Jun. 2010
STAR-CCM + 使用技巧
1 如 何 充 分 利 用 STAR-CCM + 的 一 体 化 集 成
优势?
STAR-CCM + 的操作是流程化的,整个 CAE 分
析流程都集成在 1 个界面中,用户可以完全摆脱学
习和掌握专业 CAD 造型软件,其他网格生成、表面
处理等前处理软件以及结果处理、动画制作等后处
理软件的漫长和痛苦的过程. 其中,3D-CAD 模块的
加入,更加强化一体化带来的高效便捷的优势. 如在
旋风式分离器的设计中,出口深度通常是个需要改
变的 参 数,以 寻 求 更 优 的 分 离 效 率. 用 户 可 以 在
STAR-CCM + 中绘制分离器草图时,勾选出口深度
的 Expose parameter? 选 项,即 将 其 指 定 为 Design
Parameter. 通过对其他相关边进行位置约束,即可实
现改变 Design Parameter 来改变整体几何外形的目
的( 见图 1) ,用户无须设置网格模型、尺寸及边界条
件,就可以直接生成网格并实施计算,大大缩短优化
改进设计的分析周期.
图 1 通过 Design Parameter 改变几何外形
2 多面体网格有哪些优势?
STAR-CCM + 中的多面体网格技术非常先进成
熟. 多面体网格具有六面体网格的精确度兼具四面
体网格的易生成性,在 STAR-CCM + 中是最常用的
网格类型. 多面体具有比四面体网格更好的收敛性
和更小的网格依赖性,大大降低用户的硬件资源要
求和计算时间. 用某赛车外流分析实例说明选择多
面体网格的优势,见图 2. 在该例中,若采用四面体
网格,则需要 210 万个网格才能消除网格依赖性,占
用内存 1. 3 GB;若采用多面体网格,则仅需 35 万个
网格就可消除网格依赖性,占用内存 900 MB. 图 3
和 4 分别为四面体网格和多面体网格在相同计算条
件下监控得到的阻力因数和升力因数曲线的收敛情
收稿日期: 2010-04-28
况,可以看出,后者收敛速度远快于前者.
图 2 某赛车外流场分析
图 3 四面体网格阻力因数和升力因数收敛曲线
图 4 多面体网格阻力因数和升力因数收敛曲线
3 STAR-CCM + 中如何局部加密体网格?
STAR-CCM + 中 可 以 对 局 部 区 域 内的 表 面 参
数、体网格参数等进行单独控制,常用于对空间网格
进行局部加密,见图 5.
图 5 体网格局部加密
第 2 期
STAR-CCM + 使用技巧
99
加密过程如下:首先,在管理树 Tools > Volume
Shapes 上点击右键,在 New Shape 下选择要加密区
域的 形 状,如 长 方 体 ( Brick ) 、锥 体 ( Cone ) 、柱 体
( Cylinder) 和球体( Sphere) 等. 进入编辑状态后,用
鼠标拖动或坐标输入的方式确定加密区域的大小范
围,并点击 Create 建立,在 Volume Shapes 节点下会
生成 1 个子节点( 如 Brick1) ,即新建的区域.
然后,右键点击 Continua > Mesh1 > Volumetric
Controls 选 择 New 新 建,出 现 子 节 点 Volumetric
Control1,在其属性窗口中将 Shapes 项选入前面新
建的加密区域 Brick1. 在 Volumetric Control1 > Mesh
Conditions 中 选 择 与 体 网 格 相 关 的 Mesher ( 如
Polyhedral Mesher,Trimmer 等) 并在其属性窗口中
将 Customize… 项勾 选. 在 Volumetric Control1 节 点
下会生成 Mesh Values 子节点,修改其参数可以单独
控制加密区域的体网格尺寸,实现局部加密的目的.
4 STAR-CCM + 中如何处理无厚度表面生成双面
边界层网格?
在生成体网格时,如果遇到空间中的无厚度表
面,按下述方法处理可解决拓扑封闭问题,并可在表
面两侧生成质量很高的边界层网格( 图 6) .
图 6 无厚度表面边界层网格
右键点击 Regions > Region 1 > Boundaries 节点
下的空间面名称,选择 Convert to Interface( s) ,将此
无厚度面转换成 Interface. 修改 Interfaces 节点下此
无厚度面 interface 的属性,将 Type 选择为 Baffle. 在
该 Interface 节 点 下 的 Mesh Conditions 下 Interface
Prism Layer Option 属性打勾选中.
从 Interface 节点向双面产生边界层,默认的边
界层参数是 在 Mesh Continua 里 设 定 的 全 局 参 数.
如果需要 定 义 向 双 面生 长 的 不 同 参 数 ,需 要 修 改
Regions > Region 1 > Boundaries 里,对 应 于 生 成
Interface 节 点 的 2 个 Boundary 节 点 的 Mesh
Conditions 里的 Customize Prism Mesh 参数来实现.
( 待续)
( 本文由西迪阿特信息科技( 上海) 有限公司技术部供稿. 读者
)
若对 STAR-CCM + 产品感兴趣,可以联系 support@ cdaj-china. com.
檿檿檿檿檿檿檿檿檿檿檿檿檿檿檿檿檿檿檿檿檿檿檿檿檿檿檿檿檿檿檿檿檿檿檿檿檿檿檿檿檿檿檿檿檿檿
( 上接第 97 页)
5 在 MSC Patran 中 如 何 撤 销 已 选 择 的 单 元 或
节点?
当模型较大且选择较多单元或节点时,如果使
用 Backspace 键撤销选择会浪费较多时间,但采用
以下 2 种方法可以快速撤销已选单元:(1) 双击选
项框;(2) 充分利用快捷键 Tab,同时按下 Tab 键和
Shift 键. 通过这 2 种方法就可以 1 次选中要删除的
所有单元或节点,然后 1 次性删除即可.
6 在 LS-DYNA 中 简 化 积 分 时 怎 样 避 免 或 减 小
沙漏?
在显式动力分析中采用简化积分可以极大节省
数据存储量和运算次数,并且在大变形分析中更加
适用. 但是,简化积分会出现沙漏( 零能模式) ,因此
需要有效控制分析中可能出现的沙漏变形,控制沙
漏的方法有:(1) 尽可能使用均匀的网格划分;(2)
尽量避免单点载荷;(3) 由于全积分单元不会出现
沙漏,用全积分单元定义模型的部分或全部以减小
沙漏;(4) 全局增加模型的体积黏性.
7 在 HyperMesh 中如何改变壳单元的方向?
前处理有限元软件 HyperMesh 划分好面网格
后,可能会存在同一个面上单元方向不一致的情况
( 颜色光亮度有一定差异) ,在后续施加面载荷或定
义单面接触时需要改变壳单元的方向,使同一面上
壳单元的方向一致. 修改方法为:(1) 在主菜单区选
择 Tools -> normals;(2) 通过组件或直接选择同一
面上的单元,并采用 color display normals 显示方式,
点击 display normals,由此直观地通过蓝红 2 种颜色
的 单 元 将 不 同 方向 的 单 元 区 分 开 来 ; (3) 在
orientation 中 选 择 面 上 的 1 个 单 元,点 击 adjust
normals,则所在面上壳单元的方向都改为 orientation
中所选单元的方向.
( 摘自同济大学郑百林教授《CAE 操作技能与实践》课程讲义. )
第 19 卷 第 3 期
2010 年 9 月
计 算 机 辅 助 工 程
Computer Aided Engineering
Vol. 19 No. 3
Sept. 2010
STAR-CCM + 使用技巧( 续 1)
1 如何通过 Sim 模板缩短从导入模型到执行计算
过程的工作时间?
对于几何模型不同但边界条件、计算条件甚至
网格尺寸范围相同的算例,STAR-CCM + 可将已计
算过的 sim 文件作为新模型的计算模板,大大节省
流程中网格尺寸调试、边界条件和计算条件等设定
所花费的工时.
图 1 为模板应用方法实例. 计算模型为长方体
Box,边界名称分别为壁面、入口和出口,有 1 套网格
模型物理模型及参数,同时有设定好的计算条件和
数据报告等. 新的几何模型为某分离器,几何形状差
异较大,但基本边界相同,可分为壁面、入口和出口.
网格尺度、物理模型、计算参数以及条件都可沿用
Box 设定.
Faces 上点右键选择 Split By Patch,用 Patch 方法对
其几何表面作边界分割. 用户可在编辑界面中选择
要分割的 表 面 或 在 Patch 列 表 中 选 择 表 面 对 应 的
Patch 编号,点击 Create 进行分割( 如 PIPING 几何体
表面分为 Faces 1,Faces 2 和 Faces 3 等) . (4) 将已
分割好的几何体表面与区域中的边界进行关联,即
在 Regions 下对应的区域( 如区域名称为 Body 1) 属
性中 Parts 项下选择对应的几何体( 即 PIPING) ,在
区域内各 Boundary 的属性中 Part Surfaces 项下选择
对应的几何体表面,对应关系见图 2. (5) 保留其他
设定( 可根据需要改动) . (6) 重新生成网格并进行
计算.
图 1 模板应用方法实例
以 Box 的设定作为分离器计算的模板,具体操
作步骤如下:(1) 将已经分析过的 sim 文件除几何模
型外的所有设定作为模板. 右键删除 Representations
下除 Geometry 外的所有已存网格( Import,Remeshed
Surface 和 Volume Mesh 等) ,保留包括边界在内的
其他所有设置,如壁面、入口、出口及其他各项参数
等. (2) 点击 Import Surface Mesh 按钮,导入新分析
对象的 CAD 模型,在跳出的对话框中 Import Mode
选项下选择 Create New Part,创建新的 Part,点击 OK
确认. (3) 在 Geometry →Part 下找到对应的几何体
( 如名为 PIPING 的几何体) ,在其子节点 Surfaces→
收稿日期: 2010-08-24
图 2 几何表面与物理边界关联关系
2 如何生成合理的薄壁网格?
在很多问题中常会遇到薄壁几何的情况,图 3
为流体区域中存在固体薄壁的情况示意.
图 3 流体区域中存在固体薄壁的情况示意
对于存在薄壁的区域,为满足薄壁厚度方向上
第 3 期
STAR-CCM + 使用技巧( 续 1)
901
流动或传热的计算要求,必须指定一定层数的网格
连续的情况下生成 3 层网格的形态.
以反映该方向上速度或温度等其他物理量的梯度变
化,同时在非厚度方向的网格尺寸可根据具体问题
设定为较大的值,保证整体网格数量不至于过多. 由
图 3 可知,厚度方向指定 2 层网格,非厚度方向的网
格尺寸远大于厚度,既保证厚度方向传热计算的要
求,又保证整体网格数量很少. 用 STAR-CCM + 自带
的 Thin Mesher 工具可方便地达到上述目的.
在网格连续体中选择 Thin Mesher 工具,用户可
指定相关参数控制薄壁厚度及厚度方向的网格层
数. 以某螺旋桨模型为例,叶片为薄壁,对该模型的
流固区域同时生成网格,方法如下: ( 1 ) 导入 CAD
文件,以 Split by Surface Topology 方法自动将流固分
为 2 个区域,在网格连续体 Mesh 1 的属性中勾选
Per-Region Meshing;(2) 在 Mesh 1 →Models 中右键
选择体网格模型 Thin Mesher;(3) 在 Thin Mesher 的
属性中 勾 选 Customize Thickness Threshold,用 户 就
可在 Reference Values 中定义薄壁厚度和网格层数;
(4) 在 Reference Values→Thin Mesher Layer 属性中
指定厚度方向上体网格的层数 ( 如 3 层) ,在 Thin
Solid Thickness 节点下指定薄壁厚度范围,凡小于该
值的区域厚度都被作为薄壁.
如果希望 2 区域网格连续且保证厚度方向有一
定数量层数的网格,可不做层状网格,只需在选择网
格模 型 时 选 择 Polyhedral / Tetrahedral + Embedded
Thin Mesher 即可实现. 图 4 为使用薄壁网格工具生
成的体网格形态,图 4( a) 为通过上述方法得到的层
状网格,叶片厚度方向为3 层网格;图4 ( b) 为网格
( a) 层状网格
图 4 使用薄壁网格工具生成的体网格形态
( b) 网格连续的薄壁网格
对于单区域的模型,可省略步骤(1) 进行处理.
图 5 为单区域薄壁网格示意图. 流体通道一端入口
径向厚度很薄,沿轴向其壁厚逐渐增大,较薄的一段
被划为 5 层网格以符合计算要求.
图 5 单区域薄壁网格,mm
( 本文由西迪阿特信息科技( 上海) 有限公司技术部供稿. 读者若对
STAR-CCM + 产品感兴趣,可以联系 support@ cdaj-china. com.
)
( 待续)
檿檿檿檿檿檿檿檿檿檿檿檿檿檿檿檿檿檿檿檿檿檿檿檿檿檿檿檿檿檿檿檿檿檿檿檿檿檿檿檿檿檿檿檿檿檿
西迪阿特信息科技( 上海) 有限公司招聘技术工程师
职 位:技术工程师( 发布时间:2010 年 9 月 2 日)
岗位描述:电磁场分析软件 Jmag 售前 / 售后技术支持工作
任职要求:
工作城市:上海
1 电机相关专业,本科或以上学历,硕士优先;
2 熟悉电机结构,掌握电磁场分析软件如 Ansoft,JMAG,FLUX 和 MAGNET 者优先,有相关工作经验者
优先;
3 较强的沟通、协调和人际交往能力,一定的演讲能力;
4 良好的英语阅读能力,熟练使用办公软件;
5 能够适应出差工作.
如您对以上职位感兴趣,请通过以下方式联系西迪阿特信息科技( 上海) 有限公司
地 址:上海东方路 971 号钱江大厦 25 楼 C 座
E-mail:info@ cdaj-china. com
电 话:021-50588290 /50588291
邮 编:200122
传 真:021-50588292
第 19 卷 第 4 期
2010 年 12 月
计 算 机 辅 助 工 程
Computer Aided Engineering
Vol. 19 No. 4
Dec. 2010
STAR-CCM + 使用技巧( 续 2)
1 如何使用 STAR-CCM + 的 3D-CAD 功能进行
参数化建模?
STAR-CCM + 从 V 5. 02 开始集成完全参数化
的 CAD 几何建模功能 3D-CAD. 3D-CAD 模块搭建
起 CAD 与 CAE 之间沟通的桥梁,用户无须担心在
仿真过程中使用到的不同软件之间的数据交换问
题. 3D-CAD 完全集成在 STAR-CCM + 的界面和工
作流程中,实现 STAR-CCM + 从参数化 CAD 建模到
表面准备、体网格生成,再到计算求解和后处理的一
体化工作流程. 3D-CAD 的完全参数化功能使用户
可方便地修改几何模型尺寸,快速投入到新的案例
计算中,同时也为与 modeFRONTIER 等优化软件连
接提供更大便利. 利用 3D-CAD 进行参数化建模的
基本流程为:
(1) 二维草图绘制. 通过软件提供的点、线、四
边形、圆弧和样条曲线等工具在基准平面上绘制二
维的草图.
(2) 草图特征约束. 所绘制的草图可添加各种
约束条件以确保草图的约束关系,如固定、共线、平
行、垂直和共圆心等.
(3) 三 维 实 体 成 型. 将 二 维 草 图 通 过 拉 伸、旋
转、放样和扫掠操作变为三维实体.
(4) 特征操作. 对通过倒角、倒圆和阵列等操作
编辑和优化三维实体模型.
(5) 参数设计. 在草图绘制、实体成型和特征操
作中将长度、角度和一般参数定义为设计参数. 具体
做法是勾选参数的“Expose Parameter?”选项,即将
其指定为“Design Parameter”. 通过对其他相关边进
行位置约束,可实现改变“Design Parameter”来改变
整体几何外形的目的,用户无须设置网格模型、尺寸
及边界条件就可直接生成网格并实施计算,大大缩
短优化改进设计的分析周期.
以图 1 为例说明 3D-CAD 建模的具体操作. 打
开 STAR-CCM + 后,左击树状菜单中 Geometry 旁边
的 + 号,拉出 3D-CAD Models,右击 3D-CAD Models,
在菜单中选择 New. 此时自动进入 3D-CAD 模式. 右
击 Features 下 xy,选择 Create Sketch,自动进入草图
收稿日期: 2010-10-29
绘制界面,点击利用圆心和半径做圆图标( 第 1 排
最右 1 个) ,在右边带网格线的平面上左击原点作
为圆心,用鼠标延格子边长方向移动 5 个格子后左
击确定半径,此时圆自动作出. 点击左边窗口中的
OK 按钮,结束草图绘制. 此时自动返回 3D-CAD 界
面. 仿照以上操作,进入 yz 平面内的草图绘制界面,
此时右边带网格线的平面为 yz 平面,点击利用 3 个
顶点作矩形按钮( 第 1 排左起第 3 个) ,第 1 点点击
在圆点,延 y 轴方向偏移 2 个网格点击第 2 点,再延
z 轴偏移 3 个网格点击第 3 点,然后点击左边窗口的
OK 按钮. 这样就建立了 1 个 yz 平面上的矩形. 自动
返回 3D-CAD 界面 后,可 看 到 树 状 菜 单 的 Features
下自动多出 Sketch 1 和 Sketch 2,分别对应所作出的
圆和矩形. 右击 Sketch 1,点击 Create Extrude,进入
拉伸界面,在 Distance 中填入 0. 1 m,此时右边将显
示拉伸后的形状,点击左边窗口的 OK 按钮,确认拉
伸. 同样,将 Sketch 2 向 x 轴正向拉伸 0. 1 m. 此时
CAD 模型就完成了,见图 1 右边窗口.
图 1 某简单几何模型
2 如何使用 STAR-CCM + 的流固耦合应力形变
功能?
STAR-CCM + 从 V 4. 04 开始增加流固耦合应
力形变功能,用户不需要其他第三方软件,只需在其
统一界面下即可完成流固耦合应力形变的模拟. 这
是 STAR-CCM + 作为多物理场仿真平台的重要功能
之一. 流固耦合应力形变所用到的模型为计算固体
应力的 Solid Stress 模型、计算流体变形的 Morphing
模型、计算固体变形的 Solid Displacement 模型. 以 V
5. 04 版本为例说明模型设定的要点.
021
计 算 机 辅 助 工 程
2010 年
( 1 ) 流 体 区 域 和 固 体 区 域 的 交 界 面 设 置 为
In-place Interface.
(2) 设置流体模型和固体模型. 流体模型无须
特殊设置,固体模型需选择 Solid Stress 模型,同时需
要设 置 固 体 物 性 的 泊 松 数 ( Poisson) 和 杨 氏 模 量
( Youngs Modulus) 值.
(3) 在树形管理窗口中的 Tools→Motions 节点
下添加 Morphing 模型和 Solid Displacement 模型.
(4) 分别将 Morphing 模型和 Solid Displacement
模型对应 到 流 体 计 算域 和 固 体 计 算 域 中 . 如 流 体
Regions → 计 算 域 Fluid → Physics Values → Motion
Specification,在 属 性 窗 口 中 将 Motion 设 置 为
Morphing.
(5) 将 流 体 域 中 交 界 面 Physics Conditions →
Morpher 的属性值设置为 Solid Stress.
(6 ) 将 固 体 域 中 不 能 变 形 的 边 界 Physics
Conditions→Stress / Displacement 设置为 Fixed.
图 2 为 STAR-CCM + 的流固耦合应力形变功能
示例.
图 2 STAR-CCM + 的流固耦合应力形变功能示例
3 如何使用 STAR-CCM + 的噪声计算功能?
STAR-CCM + 能预测气动噪声源,在 V 4. 06 及
后续的版本中有具体的各种气动噪声源计算模型可
选,分别适用于计算偶极子( 壁面上的压力脉动造
成) 、四极子( 流体空间内的湍流作用造成) 噪声源
等. 不同模型还有其不同的适用性,有只适用于稳态
计算的,也有适用于非稳态计算的. 使用稳态 MRF
方法,配合通常的 k-ε 或者 k-Omega 湍流模型,也可
预测流场噪声源分布情况. 但是,1 次计算只能选择
计算偶极子或四极子噪声源中的 1 种,若要得到流
场中某点上的噪声( 该点须在计算域以内) ,可使用
四极子噪声源计算方法得到在该点上由于湍流作用
造成的噪声源. 需要注意的是,以上所说的都是噪声
源的计算. 在最新的 V 5. 06 中可计算噪声传播. 另
一种方法是使用 DES 或 LES 湍流模型进行非稳态
计算,得更加细致、精确的流场压力脉动,配合监测
所需要的 点 的 压 力 值随 时 间 的 脉 动 数 据 ,再 通 过
STAR-CCM + 内 的 快 速 傅 里 叶 变 换 ( Fast Fourier
Transform,FFT) 功 能 ( V 4. 04 以 后 版 本 才 有 此 功
能) 转换到频域下的数据,作为近场噪声的参考值.
该方法需使用比较细密的网格,再加上 DES 或 LES
的非稳态计算,计算量非常大. 以下对瞬态模拟中
FFT 的设定和使用步骤进行简介.
(1) 设定压力监测点
( 见图 3) . 首先确定压力
监 测 点 的 位 置,然 后 在
Reports 目录中建立此监
测 点 的 压 力 report,在
Reports 目 录 中 生 成
Moni_P子目录. 在 Moni_P 子目录上单击右键,选择
Create Monitor from Report,则在 Monitors 目录下生
成 Moni_P Monitor 子目录. 设定 Moni_P Monitor 子
目录属性中的 Trigger 为 Time Step,这样压力监测点
的数据会按时间步进行更新.
图 3 计算模型及监控点
(2) 设定 FFT. 在 Tools→Data Set Functions 上单
击右键,选择 New Fast Fourier Transform,则建立名
为 Fast Fourier Transform 1 的 FFT 的子目录,将该子
目录 的 Amplitude Function 设 置 为 Sound Pressure
Level,其他选项保持默认. 关于 FFT 子目录的其他
属性设置参考 STAR-CCM + 的帮助文件.
( 3 ) 定 义 监 测 点 的 FFT. 在 Fast Fourier
Transform 1 → Derived 上 单 击 右 键 并 选 择 New
derived data from monitor,建立 derived monitor 子目
录,然后设置此子目录的 input data 属性为 Moni_P
Monitor.
(4) 建立 Monitor Plot 图. 在 Plots 目录上点击右
键并选择 New Plot→Monitor Plot,生成 Monitor Plot 1
子目录,将 Monitor Plot 1→Derived 中的 Derived data
属 性 设 置 为 derived
monitor.
(5 ) 打开 Monitor
Plot 1 图,执 行 计 算
时,Monitor Plot 1 的
曲线图( 见图 4) 会被
实时更新.
图 4 监控点 FFT 结果
( 待续)
( 本文由西迪阿特信息科技( 上海) 有限公司技术部供稿,读者若对
STAR-CCM + 产品感兴趣,可联系 support@ cdaj-china. com. )
第20卷第1期
2011年3月
计算机辅助工程
Computer Aided Engineering
V01.20 No.1
Mar.2011
STAR—CCM+使用技巧(续完)
1什么是Field Function?
集,使用起来非常简单.用户在编写Field Function
时,除可以使用各种表达式外,还可引用其他Field
Field Function是在STAR—CCM+中存取求解器
Function.
内单元和边界数据的机制.在STAR.CCM+中应用
Field Function的引用有其独特的方法.Scalar
十分广泛,而其他CFD软件绝大多数需使用用户子
程序才能实现有关功能;无须安装特别的编译器,具
有功能强大、简单易用的特点.
2
Field Function有什么用处?
Field Function常用于显示计算结果,如要显示
某个壁面上的压力分布云图,此处“压力”
(Pressure)就为一个Field Function.
Field Function还可用于定义边界和区域上的值
或定义初始条件.此时,用户一般需事先编写好一个
Field Function,然后将其应用到边界、区域或初始条
件的设定中.如要将一个速度入口边界的速度大小
值赋成随Y坐标呈抛物线规律变化,就需事先编写
一个满足该规律的Field Function,然后应用到速度
入口边界的设定中.
3
Field Function的分类
从Field Function的生成方式看,可分为System
和User这2类.System Field Function由软件根据所
选择的物理模型自动生成,User Field Function为用
户白定义的Field Function,定义时可引用System
Field Function或其他User Field Function.
从Field Function的数据类型看,可分为Scalar,
Vector,Array和Position等4种.Scalar Field Function
Field Function的引用是在函数名前加一个“$”符
号,如
$Temperature
Vector,Array和Position Field Function的引用则是
在函数名前加2个“$”符号,如
$$Velocity
其分量的引用方式为
8$Velocity[0]——i方向分量
$$Velocity[1]_f方向分量
$$Velocity[2]——l|}方向分量
Vector和Position Field Function在引用时还可指定
坐标系,如
$$Centroid(”Coordinate System 1”)
Field Function在定义时可灵活使用表1—5所
示的常用运算符、函数和表达式.其中,矢量函数适
用于Vector,Array和Position Field Function.
表1算术运算符
Tab.1
Arithmetic operators
表2逻辑运算符
Tab.2
Logical operators
为只有1个分量的标量场函数,如Temperature,
运算符
I=
●
>
(
≥
≤
&& ||
Pressure等.Vector,Array和Position Field Function
都具有3个分量,如Velocity属于Vector Field
Function,具有i,j,k等3个方向分量,另外还有合值
(Magnitude)这个属性.不过,Vector和Position Field
Function可在不同的坐标系互相转换,而Array Field
Function不可以.
4
Field Function的语法
Field Function用到的语法为C语言的一个子
功能
等于 不等于 大于 小于
大于 小于 逻辑 逻辑
等于 等于
与
或
表3三角函数
Tab.3
Trigonometric functions
函数名 mrcc08($x) aresin(1 x) arctan($x) arctan2($Y。$z) cos($x)
功能
反余弦
反正弦
反正切 y/x的反正切 余弦
函数名 cosh($z)
sin($x)
sinh($z)
tanh(s£)
tan($z)
功能 双曲余弦
正弦
双曲正弦
双曲正切
正切
收稿日期:2011-01-14
万方数据
http://www.ehinacae.C/1
计算机辅助工程
2011正
(8 nme≥0 01)’l 00。:300+70呻0+
$Time
例3在同一个边界不同位置定义不同的热{j|c
密度Y)0的部分热流密度是1 W/m 2.Y‘0的部分
W/m‘(注:$$Ceatroid和
热流密度是0 1
s$Positloa具有类似功能)
【S 8 P∞ition[I]>0 0)?1:01
倒4对不I可空问位置定望不l司的初始遵度
场在z‘4的空间,速度初场为[2.0.0]m/s;在
*≥4的部分.速度初场为[10,0,01一s速度韧场
为一个矢量场,需定义一个Vector Field Function
[($$Centmid『0j<4)'2:10,0,0]
倒5基于三维矢量场的分量定义一雏或二维
矢量场仅包古f方向速度的一维矢量场
f 3$Vel“:itv『0],0.0]
包含i』方向速度的二维矢鼠场
l 8$Vel叫ity【0j,}sVel∞ity[1],0】
倒6显示单元表面法线矢量函数$$^m是
单元表面面积的矢量,其矢量方向为沿单元表面的
法线方向,除以单兀表面面积的标量值函数
$8^瑚mug(),即得到单位矢量方向新建一个
Vector User Field Function
$$Ar蚰/$$Area mag()
在先量场(Vector
Field)里选择该新建的
User Field Function,显
示在单元表面上即可,
见剧2
例7定义衙c场衍
生参数
目2采面法线*量
o—。
2
5…rfa
m6
速度矢量场的散度甲.y:“E
div(S S Vel∞ib,)
密度与速度乘积的散度甲(p”)
div($Density+$$VelociW)
眶力场的梯度V·P
grad($P…)
表5#他常月自敲
Tatf 5^Ii∞eⅡane咄f帅cdon5
Field/:'unction中的条件表达式甘写方式为
(Ⅱ>b)’c:d
(1)
其含义为:如果n>6成立.则表达式值为‘;如果
n>6不成立,则表达式值为d
条件表达式可以嵌套,如
($$Cenlmid『1]>2)7 30:
(($$Centroid[1](1)’30:0)
对照式(I)可知.d
在式(2)中为条什表达
(2)
式整个表达J£的青义为
三段式分段函数.见
‘…
哺1
7 5
目1分段自救
5
Field IFmacfion应用举例Fill 1
s‘叩funell眦
例1缩放已有的Field Function已知x方向
的质心坐标晒数为$$Centmid[01,其单位为m欲
新建一个以英寸为单位的*方向质心坐标函数.需
对原始甬数进行缩放.除以0 025 4即可:
$$Centroid『0I/0 025 4
倒2定义瞬态计算中随时问变化的温度边界
速度场的旋度V”
条件已知温度边界条件的变化规律是在0叭s前
从300 K线性变化到I 000 K;在0 0l s后保持
1 000 K
curl(8$Vel们ity)
(#i自日4H##b#&(±≈)#m4日*$**#4}*”
$TAR_f2CM+,*$**,TE}joppn@dm小啪∞m)
(续完)
万方数据