CFD课程设计报告说明书.docx
《CFD课程设计报告说明书.docx》由会员分享，可在线阅读，更多相关《CFD课程设计报告说明书.docx（22页珍藏版）》请在冰豆网上搜索。
CFD课程设计报告说明书
CFD课程设计
翼身组合体流场分析
院系：
航空航天工程学部
专业：
飞行器设计与工程
班级：
24030301
学号：
2012040303023
姓名：
摘要
此次课程设计是利用ANSYS软件中的ICEM和Fluent求解器计算不同迎角下，翼身组合体的升力系数，阻力系数，力矩系数以及各个状态下的流场分布情况，机身为方截面机身，机翼为三角上单翼，翼型选择NACA4412，计算结束后，利用tecplot软件绘制Cy-α，Cy-Cx，Mz-Cy曲线，得出Cy0，最大升阻比等气动力特征参数。
关键词ICEMFluent翼身组合体tecplot
第一章绪论
1.1ANSYS软件介绍
ANSYS软件是融结构、流体、电场、磁场、声场分析于一体的大型通用有限元分析软件，是一个多用途的有限元法计算机设计程序，可以用来求解结构、流体、电力、电磁场及碰撞等问题。
由世界上最大的有限元分析软件公司之一的美国ANSYS开发，它能与多数CAD软件接口，实现数据的共享和交换，如Pro/Engineer，NASTRAN，Alogor，I－DEAS， AutoCAD等，是现代产品设计中的高级CAD工具之一。
在此次的课题中，主要用到其中的ICEM及Fluent部分。
1.2主要内容
本次课程设计的主要内容就是通过CATIA建立机身和机翼的组合体模型，通过fluent解算器进行有限元分析，从而得到该组合体的一些相关的气动数据。
此次课程设计的重点在于模型的建立，通过CATIA建立基础的模型，然后导入到ANSYS.ICEM中进行模型的处理以及网格包括壳网格、体网格及附面层网格的划分。
完成之后导入到fluent解算器设置属性，相关参数等，然后进行计算不同迎角下的翼身组合体的相关气动参数及压力云图分布情况。
第二章模型的建立
2.1CATIA建立模型及导出
通过CATIA建立模型，机身为方截面机身，机翼为三角上单翼，展弦比为3，翼型选取NACA4412，翼身组合体及流畅区如图所示。
其中机身2m，弦长1m，前场8m，后场16m。
图2.1翼身组合体模型
图2.2流场区域
选择文件-另存为，在保存类型里选择model，然后选择保存即可。
图2.3导出模型
第三章ANSYS.ICEM处理
3.1导入模型
打开ANSYS.ICEM，设置后工作目录，然后选择file-ImportGeometry，选择CATIAV4，选择保存的model文件打开，单击Apply。
然后选择Geometry下的RepairGeometry修复模型。
均为红线，没有问题。
划分Part，定义对称面、入口、出口、机翼、机身及流场边界。
图3.1导入模型
3.2划分网格
（1）壳网格划分
在Mesh选项卡下，选择GlobalMeshSetup，设置全局网格大小为800，点击Apply；
图3.2.1全局网格尺寸
定义壳网格参数，点击Apply，如图；
图3.2.2壳网格参数
定义Part的网格尺寸，其中对称面，流场边界，以及入口出口的maxsize设置为1000，机身设置为20，机翼设置为10；选择Mesh下的ComputeMesh，选择生成壳网格点击Compute；
图3.2.3生成壳网格
查看网格质量；发现网格生成的质量还可以，机身附近有明显的加密情况。
图3.2.4壳网格
（2）体网格及附面层网格划分
设置体网格和棱柱网格参数，点击Apply；
图3.2.5体网格参数
图3.2.6棱柱网格参数
在机翼后缘向后创建网格加密区；然后设置要生成附面层的Part，勾选机翼和机身的part，选择生成体网格，勾选CreatePrismLayers，点击Compute生成网格；
图3.2.7生成体网格
检查对称面处棱柱网格生成情况，发现有良好的附面层生成，如图。
图3.2.8附面层网格
3.3导出网格
选择EditMesh，检查网格质量如图；网格质量良好，保存网格。
图3.3.1网格质量
选择Output，选择OutputSolver为FluentV6，点击Apply；选择WriteInput，选择刚保存的网格打开，命名导出的文件名及文件路径，点击Done。
图3.3.2导出网格
第四章Fluent计算
4.1设置参数计算
（1）定义网格
打开Fluent，选择File-Read-case，选择保存的msh网格文件，打开。
在General下，选择Scale，在MeshWasCreatedIn下拉列表中选择mm，点击scale，然后关闭。
选择Check，检查网格，MinimumVolume应大于1。
Solver框里的VelocityFormulation中选择Relative。
（2）定义求解模型
选择Models-Viscous，双击，选择Spalart-allmaras（1eqn）模型。
图4.1.1求解器模型
选择Materials，定义材料，默认为空气，在编辑菜单中的Density中选择Idel-gas，Viscosity栏中选择sutherland，在弹出的菜单中选择OK，点击Change/Create，然后点击Close关闭。
图4.1.2流体材料
（3）定义边界条件
定义流场域材料，在zone中选择airplane，在type栏中选择fluid，及之前定义的air的fluid材料。
定义壁面，在机翼和机身的type类型中选择wall，弹出对话框点击ok默认。
定义对称面，在对称面的part的type下选择symmetry，点击OK默认。
定义远场，在入口，出口以及流场边界三个zone的type栏中均选择pressure-far-field，在弹出的对话框中，设置machnumber为0.5，输入来流的方向向量的三维坐标值（改变迎角），Temperature输入300，点击OK。
图4.1.3定义远场
（4）初始化计算
选择ReferenceValues，在ComputeFrom下拉栏中选择入口的part，在Area栏中输入参考面积（0.376mm2）。
选择SolutionControls，定义松弛系数，均为默认值的一半，
选择Monitors，定义监视器。
显示残差曲线（默认显示），设定各个参数的收敛残差值为1e-3，点击OK；显示升力系数变化曲线，点击Create，选择Lift，勾选plot，在zone列表中选择机身和机翼，输入坐标（0，0，1），点击OK；显示阻力系数变化曲线，点击Create，选择Drag，勾选plot，同样在zone列表中选择机身和机翼，坐标输入（-1，0，0），点击OK；显示力矩系数变化曲线，点击Create，选择Moment，勾选plot，zone选择机身和机翼，MomentCenter输入（-0.67，0.32，0.032），MomentAxis输入（0，1，0），点击OK。
图4.1.4定义监视器
选择SolutionInitialization，点击Initialize初始化流场。
选择RunCalculation，在NumberofIterations栏中输入迭代次数，这里输入600，点击Calculate，开始计算，在大约330步左右达到收敛要求，计算结束，改变来流方向，重新计算。
记录不同来流方向下的计算结果及压力分布云图。
4.2计算结果
此次计算状态为翼身组合体在理想气体中进行计算，其中，流体速度为0.5Ma，迎角有-4-18度变化范围，共十次计算，得到不同迎角下的Cy，CX，Mz以及压力云图如下：
表4.2.1气动参数表
a
-4
0
4
6
8
10
12
14
16
18
Cl
0.144
0.451
0.704
0.833
0.951
1.069
1.186
1.295
1.393
1.477
Cx
0.003
0.004
0.017
0.031
0.046
0.062
0.078
0.093
0.106
0.116
Cm
-0.037
-0.042
-0.068
-0.07
-0.08
-0.095
-0.103
-0.106
-0.107
-0.108
图4.2.1-4度压力云图
图4.2.20度压力云图
图4.2.34度压力云图
图4.2.46度压力云图
图4.2.58度压力云图
图4.2.610度压力云图
图4.2.712度压力云图
图4.2.814度压力云图
图4.2.916度压力云图
图4.2.1018度压力云图
第五章数据处理分析
5.1气动参数曲线
通过计算结果，使用tecplot软件画出Cy-α，Cy-Cx以及Mz-Cy曲线，如图：
图5.1.1Cy-α曲线
图5.1.2Cy-Cx曲线
图5.1.3Mz-Cy曲线
通过对以上曲线分析得到参数如下表：
表5.2.1气动力特征参数表
气动参数
Cya
Cy0
Kmax
Cx0
Mzcy
Mz0
数值
0.0583
-3.7
90.138
0.011
-0.037
-0.042
参考文献
[1]陈再新等.《空气动力学》.航空工业出版社2014
[2]ANSYS公司.Fluent软件用户帮助.航空工业出版社2000
[3]ANSYS公司.ANSYSICEM用户帮助.航空工业出版社1997
升级会员 