查看: 5466|回复: 2
Fluent 经典问题答疑[6]
61 FLUENT help和GAMBIT help能教会我们（特别是刚入门的新手）学习什么基本知识？
可以了解其基本原理和基本的操作。不过我觉得对于新手熟悉软件最好的还是tutorial guide
& && && &个人认为：软件的帮助文档还是很有用的，无论对不初学者还是已经相对熟练的“老手”都具有参考意义。
63 FLUENT模拟飞行器外部流场，最高MA多少时就不准确了？MA达到一定的程度做模拟需注意哪些问题？
不准确的标准是什么？没有判断标准就没办法判断。一般来说fluent计算马赫数大于3-5之后就不是很理想了（不过相信版本越新结果越好）。计算的时候应该从低马赫数慢慢往上算。比如说如果计算马赫数是5的话，就在马赫数4的计算结果上算，另外，求解器需选择耦合和显式的。（对于6.3来说，选择基于密度的求解器）
68 做飞机设计时，经常计算一些翼型，可是经常出现计算出来的阻力是负值，出现负值究竟是什么原因，是网格的问题还是计算参数设置的问题？
如果这个问题对于某个人经常出现的话，那就比较奇怪了，阻力是负值，难道就是传说中的前缘吸力现身？呵呵，只是开个玩笑：），估计肯定是计算错了或者是设置错了。在飞机翼型气动里面，阻力主要有两种成份：压差阻力和摩擦阻力。应该是正值的。
排除是计算过程的其他问题，我觉得在使用Fluent进行这方面的计算时，需要注意两个方面：
1)参考值的设置，也就是Report-&Reference Values...
这些参考值，是用来计算Re，以及升力，阻力，力矩系数所要用到的。如果设置不当，即使计算过程是对的，所得到的升阻力等系数也是不对的。对于2D翼型仿真计算，比较容易出错的就是里面的Area该写什么，单位是平方米，这里应该填写翼型的弦长（Chord Length），The area here is actually area per unit depth；就是每单位展长的面积。
2)在监视力的时候，关于力的矢量方向设置，Solve-&Monitor-&Force...
这个矢量方向千万不要小看，不能填错，填错了就可能出现阻力是负值的错误，Fluent之前的版本所附带的例子，关于NACA0012翼型的计算中，这里的矢量就设置错了，受错误例子的影响，韩占忠那本书中三角形翼型的那个例子也设置错误，在书的第112页的第6步的第（7）小步就设置错误，升力系数的力方向矢量，应该是X=-0.087155，Y=0.996195；前面他也写到要注意：要确保阻力和升力分别与来流平行和垂直，那么这两个力矢量肯定是垂直的了，那么这两个矢量的点乘就肯定等于零了；所幸的是，在Fluent6.3版本的例子中，这个错误已经改正过来了。
74 大概需要划分100万个左右的单元，且只计算稳态流动，请问这样的问题PC机上算的了吗？如果能算至少需要怎样的计算机配置呢？
一般来说，按照1000个节点对1MB内存这样预估就差不多了，只计算稳态流动，pc机应该差不多了，不过因为一般的pc机可能在连续计算5、6天之后就出现浮点运算错误，所以如果计算不是很复杂，采用的求解器和湍流模型不是太好计算资源，应该还是可以的。
如果使用pc机计算，建议至少采用2GB内存，主板最好固态电容，不易爆浆，电源最好功率大典，应该差不多了，现在流行四核cpu的，可以考虑使用四核的，这样的配置下来也不比服务器差多少。
76 GAMBIT划分三维网格后，怎样知道结点数？如何知道总生成多少网格（整个模型）？
个人一般是将网格读入fluent后，通过grid-&info-&size来看：）
77 在FLUENT的后处理中可以显示一个管道的。某个标量的。圆截面平均值沿管道轴线（中心线）的变化曲线吗？何显示空间某一点的数值呀（比如某一点温度）？
先创建一条ling（中心线），然后在xyplot中生成曲线
80 如何在gambit中输入cad和Pro/e的图形？如何将FLUNET的结果EXPORT成ANSYS的文件？
autocad需要将图形转化为sat格式，pro/e可以将文件转化为igse或者stp格式。在fluent的flie/export 中可以选择导出ansys格式的文件
81 入口和出口处的k和epsilon值怎么设置？
可参考陶文铨教授的《数值传热学》、《计算传热学的近代进展》，里面有非常具体的介绍，且容易掌握和实施。
87&&courant数：在模拟高压的流场的时候,迭代的时候总是自动减小其数值，这是什么原因造成的，为什么？怎么修改？
这是流场的压力梯度较大，Fluent自身逐步降低时间步长，防止计算发散。一般的处理办法是：先将边界条件上的压力设置较低点，使得压力梯度较小一点，等到收敛的感觉差不多，在这个基础上，逐渐把压力增大，这样就不容易发散。
94 把带网格的几个volume，copy到另一处，但原来split的界面，现在都变成了wall，怎么才能把wall变成内部流体呢？
直接边界面定义为interior即可
98 Gambit的网格相连问题：如果物体是由两个相连的模型所结合，一个的网格划分比较密、另一个比较稀疏，用Gambit有办法将两个网格密度不同的物体，相连在一起吗？
请参考第16题答案。将两种网格交界的地方设置成一对interface即可。
100 在FLUENT里定义流体的密度时，定义为不可压理想流体是用在什么地方呀，讲义上说是用于可变密度的不可压流动，不知如何理解？
define/matirial 中定义。可变密度的不可压缩流动，就是说在该流动下，流体介质的密度可以认为不变。比如说空气在流速在0.3马赫的情况下都可以认为是密度不变的
101 已经建好的模型，想修改一些尺寸，但不知道顶点的座标，请问如何在gambit中显示点的座标？
在gambit中的geometr-〉vortex-&summarize vortices即可显示点的坐标。
102 在FLUENT模拟以后用display下的操作都无法显示，不过刚开始用的是好的，然后就不行了，为什么？
DirectX 控制面板中的“加速”功能禁用即可
103 能否同时设置进口和出口都为压力的边界条件？在这样的边界条件设置情况下发现没有收敛，研究的物理模型只是知道进口和出口的压力，不知道怎么修改才能使其收敛？
当然可以同时设置进口和出口都为压力的边界条件。如果没有收敛，需要首先看看求解器、湍流模型、气体性质和边界条件时有没有出现warning；其次，还是我上边的帖子所说的，对于可压流动，采用压力边界条件，不能一下把压力和温度加到所需值，应该首先设置较低的压力或温度，然后逐渐增大，最后达到自己所需的值。
104 在FLUENT计算时，有时候计算时间会特别长，为了避免断电或其它情况影响计算，应设置自动保存功能，如何设置自动保存功能？在非定常计算中读入自动保存文件时如下出现问题：
Writing &F:\propane\16\160575.cas&...
Error: sopenoutputfile&: unable to open file for output
Error Object: &F:\propane\16\160575.cas&
Error: Error writing &F:\propane\16\160575.cas&.
Error Object: #f
非定常的，算了一段之后停下来，改天继续算的时候，自动保存的时候出现问题，请问如何解决？
File-&write-&Autosave就可以实现自动保存，自动保存的是date文件阿，你的怎么是CASE文件?
105 gambit划分时运动部分与静止部分交接面：一个系统的两块，运动部分与静止部分交接部分近似认为没有空隙（无限小，虽然实际上是不可能的），假设考虑做成一个实体，那么似乎要一起运动或静止；假设分开做成两个实体，那么交接处的两个不完全重合的面要设为WALL还是什么呢，设成WALL不就不能过流了 吗？&&
将这一对接触面设置成Interface就行了，具体请参考第47题的解答。
108 想把gambit的图形保存成图片，可是底色怎么做成白色
首先点开GAMBIT的EDIT菜单,其次点GRAPHICS,在下拉列表中点到 WINDOWS BACKGROUND COLOR BLACK 一项& &在下面VALUE 中填写WHITE,再点左面的MODIFY,就可以了.
110 在分析一个转轮时，想求得转轮的转矩，不知道fluent中有什么方法可以提供该数据。本来想到用叶片上面的压力乘半径，然后做积分运算，但是由于叶片正反壁面统一定义的，即全部定义为wall-rn1，所以分不出方向来了。
report/force/moment 定义需要计算的面和旋转中心就ok了
111 如何在gambit中实现坐标轴的变换：有一个三维的网格，想在柱坐标中实现，可是gambit中一直显示直角坐标？
对于这个问题，你可以尝试一下：Operation---&Tools---&Coordinate System---&Activate Coordinate System.
113 利用vof非稳态求解，结果明显没有收敛的情况下，为什么就开始提示收敛,虽然可以不管它，继续算下去达到收敛。但是求解怎么会提前收敛？
好像非稳态不存在收敛这个概念吧。（除非是双时间推进中）
116 在Gambit中如何将两个dbs文件到入：把炉膛分成了三个dbs文件，现在想导入两个dbs文件，在Gambit中进行操作，但好象使用open命令就只能open一个dbs文件，请问这要怎么处理？
将其中一个导出成iges或者别的格式，然后就能和dbs一起导入了
120 scale是把你所画模型中的单位转化为Fluent默认的m，而unite是根据你自己的需要转化单位，也就是把Fluent中默认的m转画为其他的单位，两中方法对计算没有什么影响吗？ scale是对几何进行比例缩放，而unit只是改变单位，不改变几何外形的大小。比如，一个是1m的几何外形，通过scale将m变为mm，那么几何外形就变成了1mm。如果通过unite将m改为mm，那么几何外形不变，还是1000mm，只是表示的单位变成mm了
27 什么叫PDF方法？FLUENT中模拟煤粉燃烧的方法有哪些？
概率密度函数输运输运方程方法 (PDF方法)是近年来逐步建立起来的描述湍流两相流动的新模型方法。所谓的概率密度函数(Probability Density Function,简称PDF)方法是基于湍流场随机性和概率统计描述，将流场的速度、温度和组分浓度等特征量作为随机变量，研究其概率密度函数在相空间的传递行为的研究方法。PDF模型介于统观模拟和细观模拟之间，是从随机运动的分子动力论和两相湍流的基本守恒定律出发，探讨两相湍流的规律，因此可作为发展双流体模型框架内两相湍流模型的理论基础。它实质上是沟通E-L模型和E-E模型的桥梁，可以用颗粒运动的拉氏分析通过统计理论，即PDF方程的积分 建立封闭的E-E两相湍流模型。
非预混湍流燃烧过程的正确模拟要求同时模拟混合和化学反应过程。FLUENT 提供了四种反应模拟方法：即有限率反应法、混合分数PDF 法、不平衡（火焰微元）法和预混燃烧法。火焰微元法是混合分数PDF 方法的一种特例。该方法是基于不平衡反应的，混合分数PDF 法不能模拟的不平衡现象如火焰的悬举和熄灭，NOx 的形成等都可用该方法模拟。但由于该方法还未完善，在FLUENT 只能适用于绝热模型。
对许多燃烧系统，辐射式主要的能量传输方式，因此在模拟燃烧系统时，对辐射能量的传输的模拟也是非常重要的。在FLUENT 中，对于模拟该过程的模型也是非常全面的。包括DTRM、P-1、Rosseland、DO 辐射模型，还有用WSGG 模型来模拟吸收系数。
30&&FLUENT运行过程中，出现残差曲线震荡是怎么回事？如何解决残差震荡的问题？残差震荡对计算收敛性和计算结果有什么影响？在进行稳态计算时候，开始残差线是一直下降的，可是到后来各种残差线都显示为波形波动，是不是不收敛？
1)残差波动的主要原因：高精度格式；网格太粗；网格质量差；流场本身边界复杂，流动复杂；模型的不恰当使用。
有些复杂或流动环境恶劣情形下确实很难收敛。计算的精度（2 阶），网格太疏，网格质量太差，等都会使残差波动。经常遇到，一开始下降，然后出现波动，可以降低松弛系数，我的问题就能收敛，但如果网格质量不好，是很难的。通常，计算非结构网格，如果问题比较复杂，会出现这种情况，建议作网格时多下些功夫。理论上说，残差的震荡是数值迭代在计算域内传递遭遇障碍物反射 形成周期震荡导致的结果，与网格亚尺度雷诺数有关。例如，通常压力边界是主要的反射源，换成OUTFLOW 边界会好些。这主要根据经验判断。所以说网格和边界条件是主要因素。
2)(1)网格问题：比如流场内部存在尖点等突变，导致网格在局部质量存在问题，影响收敛。
(2)可以调整一下courant number，courant number实际上是指时间步长和空间步长的相对关系，系统自动减小courant数，这种情况一般出现在存在尖锐外形的计算域，当局部的流速过大或者压差过大时出错，把局部的网格加密再试一下。
& & 在fluent中，用courant number来调节计算的稳定性与收敛性。一般来说，随着courant number的从小到大的变化，收敛速度逐渐加快，但是稳定性逐渐降低。所以具体的问题，在计算的过程中，最好是把courant number从小开始设置，看看迭代残差的收敛情况，如果收敛速度较慢而且比较稳定的话，可以适当的增加courant number的大小，根据自己具体的问题，找出一个比较合适的courant number，让收敛速度能够足够的快，而且能够保持它的稳定性。
31&&数值模拟过程中，什么情况下出现伪扩散的情况？以及对于伪扩散在数值模拟过程中如何避免？
假扩散（false diffusion）的含义：
基本含义：由于对流—扩散方程中一阶导数项的离散格式的截断误差小于二阶而引起较大数值计算误差的现象。有的文献中将人工粘性（artificial viscosity）或数值粘性（numerical viscosity）视为它的同义词。
拓宽含义：现在通常把以下三种原因引起的数值计算误差都归在假扩散的名称下
1)非稳态项或对流项采用一阶截差的格式;
2)流动方向与网格线呈倾斜交叉(多维问题);
3)建立差分格式时没有考虑到非常数的源项的影响。
为克服或减轻数值计算中的假扩散(包括流向扩散及交叉扩散)误差，应当：
1) 采用截差阶数较高的格式；
2) 减轻流线与网格线之间的倾斜交叉现象或在构造格式时考虑到来流方向的影响。
3) 至于非常数源项的问题，目前文献中，还没有为克服这种影响而专门构造的格式，但是高阶格式显然对减轻其影响是有利的。
32 FLUENT轮廓（contour）显示过程中，有时候标准轮廓线显示通常不能精确地显示其细节，特别是对于封闭的3D物体（如柱体），其原因是什么？如何解决？
FLUENT等高线（contour）显示过程中，可以通过调节显示的水平等级来调节其显示细节，Levels...最大值允许设置为100.对于 封闭的3D物体，可以通过建立Surface，监视Surface上的量来显示计算结果。或者计算之后将结果导入到Tecplot中，作切片图显示
33 如果采用非稳态计算完毕后，如何才能更形象地显示出动态的效果图？
对于非定常计算，可以通过创建动画来形象地显示出动态的效果图。
Solve-&Animate-&Define...，具体操作请参考Fluent用户手册。
34&&在FLUENT的学习过程中，通常会涉及几个压力的概念，比如压力是相对值还是绝对值？参考压力有何作用？如何设置和利用它？
GAUGE PRESSURE 就是静压。
GAUGE total PRESSURE 是总压。
这里需要强调一下 Gauge为名义值，什么意思呢？如果， INITIAL Gauge PRESSURE ＝0 ,那么 GAUGE PRESSURE 就是实际的静压Pinf。
& &GAUGE total PRESSURE 是实际的总压Pt。如果INITIAL Gauge PRESSURE 不等于零
GAUGE PRESSURE ＝ Pinf - INITIAL Gauge PRESSURE GAUGE total PRESSURE ＝ Pt - INITIAL Gauge PRESSURE
35 在FLUENT结果的后处理过程中，如何将美观漂亮的定性分析的效果图和定量分析示意图插入到论文中来说明问题？
1)在Fluent中显示你想得到的效果图的窗口，可以直接在任务栏中右键该窗口将其复制到剪贴板，保存；或者打印到文件，保存。
2)在Fluent中，在你想要保存相关窗口的效果图时，首先激活效果图监视窗口，就是用鼠标左键监视窗口，然后在Fluent中操作，Fluent-&File-&Hardcopy...，选择好你想要的图片格式，然后就可以保存了。
3)将计算结果或者相关数据导入到Tecplot中，然后作出你想要的效果图，这种方法得出的图片，个人感觉比Fluent得到的图片美观简洁大方。
Powered by帐号:密码:下次自动登录{url:/nForum/slist.json?uid=guest&root=list-section}{url:/nForum/nlist.json?uid=guest&root=list-section}
贴数:1&分页:水印发信人: mildwater (水印), 信区: NumComp
标&&题: 请教：Fluent气液两相流中的升力系数
发信站: 水木社区 (Fri Mar 20 15:48:27 2009), 站内 &&&&&& 气液两相流，使用Eulerian 模型，用udf自己加曳力和升力。 &&&& 在user's guide 里有关于曳力系数的说明，是FD=KD(V(p)-V(q))这样的形式，KD为曳力系数，曳力已经成功加入。关于升力系数user's guide里没有这种说明，但在udf的说明里又说可以用和曳力系数一样的那个宏来设置升力系数，我拿曳力的表达式类比了一下弄出一个FL=KL(V(p)-V(q))的样子，KL为升力系数，但是很快就导致发散了。 &&&& 不知道哪位碰到过类似的问题，升力系数到底应该是什么样子的，不胜感激。
-- && ※ 修改:·mildwater 于 Mar 20 15:51:14 2009 修改本文·[FROM: 166.111.35.*]
※ 来源:·水木社区 ·[FROM: 166.111.35.*]
文章数:1&分页:Fluent使用中的参数设定经验（不定时更新） - 博客频道 - CSDN.NET
Spark_Alice的博客
分类：Fluent
在受壁面限制的流动中，因为壁面附件流场变量的梯度较大，所以壁面对湍流计算的影响很大，因此在壁面附件要进行特殊处理。一种办法是用半经验公式将自由流中的湍流与壁面附件的流动连接起来，这种方法被称为壁面函数法。另一种方法是通过在壁面附近加密网格，同时调整湍流模型以包含壁面附件低雷诺数流动的影响，这种方法称为近壁模型法。。
当在数据库中的物理属性找不到我们所需要的物质的物理参数时，可以根据相关资料直接在图1-47的Name文本框中输入相应物质的名称，并在Properties对应的选项中填入已知参数值，添加一个新物质。
因为使用了湍流模型，所以还要求设置入口的湍流参数，从湍流参数的指定方式（Specification Method）下拉列表中选择湍流强度和水力直径。关于湍流入口参数的设置一般都是经验性质的，这个参数的设置一般影响一定范围内的流动参数分布。
Outflow其物理意义是充分发展，Under-Relaxation
Factor选项表示所求解的控制方程和一些变量的松弛因子。
一般来说，初始解对求解的影响比较大，所以给出的初始解要尽量接近真实解。一般是用进口的量对全场进行初始化。残差的值越小表示计算的精度要求越高。选择Residual，在Options下面的Print和Plot复选框都选中
湍流模型中的大涡模拟(LES)和离散涡模拟（DES）只能用于非定常流动的模拟。由流体粘性所带来的边界层是尾迹形成的原因
PISO: pressureimplicit split operator, SIMPLE系列算法计算过程大致分为两步：预估步和校正步，显式修正步。SIMPLE和SIMPLEC算法常被用于定常流动的计算，PISO算法虽也可以用来计算定常流动，但更常用于计算非定常流动以及网格偏斜较为严重的情况。
密度基求解方法是针对可压缩流动而设计的，
虽然Density Based求解器更适合求解可压缩流动问题，但是Pressure
Based耦合求解器也很适合求解带有激波的高速空气动力问题，在GradientOption下的Green-Gauss
Node Based，基于节点的高斯克林函数求梯度的方法比基于控制体中心的精度要高，特别适合分结构化网格，可以更为精确地计算升阻力
Density选择ideal-gas，在Viscosity右侧选择sutherland，选择三系数方法，用Sutherland定律求粘性非常适合可压缩流动。在Specification
Method右侧下拉框中选择Turbulent ViscosityRatio为10的默认设置，对于外部绕流，粘度比通常为1~10之间。
对高阶差分格式推荐较低的松弛因子，在explicit relaxation factors下的Momentum和Pressure为0.5.
执行Full Multigrid（FMG）初始化，FMG初始化可以使初场更为合理，容易在较大的库郎数下启动，可以在较小的迭代次数内获得收敛
在设置阻力、升力收敛曲线监视器之前先进行少量次数的迭代求解，可以壁面监视器初始阶段过大的振荡变化，缩小了纵坐标方向的数值范围，更容易判断收敛情况。
设置翼型阻力系数监视器，plot前面打勾是为了动态显式阻力系数的变化，write前面打勾表示阻力系数随迭代过程变化的情况被写入一个文件，
设置用于计算升力系数、阻力系数和力矩系数的参考值，这里参考值是用于对作用于翼上的力进行无量纲化，得到的无量纲的力和力矩系数分别为升力、阻力、力矩系数。阻力和升力系数的定义一般为阻力、升力与来流动压的比值，具体的表达式为：
排名：千里之外
（2）（0）（0）（1）（2）（0）（3）第3页/共3页
70机械设计与制造
在拍动频率为40Hz，时间步长设为0.001s，计算200个时间
0.080.060.040.02CI
0.00-0.02-0.04-0.06-0.08
No.5May.2011
升力系数变化曲线，如图7所示。步长，迭代了4000次后，
0.080.060.040.02CI
0.00-0.02-0.04-0.06-0.08
00.0050.01
图10纵坐标为升力系数，横纵坐标拍动时间
图7纵坐标为升力系数，横纵坐标为拍动周期数
计算流体力学分析软件Fluent提供的非结构化生成程序把计算相对复杂的几何结构问题变得容易和轻松。应用时可以根据网计算过程中网格的变化实时调整网格以满足计算要求。同时，格自适应和调整只是在加密的流动区域里实施，而不是整个流动场，因此可以节省计算时间。应用动网格技术可以有效地解决数值求解流场的N-S方程所存在的“动边界”问题。结果表明应用Fluent进行微扑翼飞行器的动力学仿真其结果是基本可靠的。
研究表明翅翼以不同的模式拍打运动将会产生不同的升力和升力系数。以铰链四杆机构为翅翼机构所得到的升力系数稍大于翅翼匀速拍动和以正弦规律运动所得到的升力系数。
研究的是二维流场中翅翼运动产生的升力和升力系数。对于更加复杂的翅翼的运动规律，则问题将会变得更加复杂，还需要进行更加深入的研究。
翅翼在下拍到最低点附近时流场的速度及压力分布，如图8、图9所示。如图8所示，可以看出，当翅翼在下拍到最低点时，上表面的空气的流速大于下表面空气流速，并且可以看到绕流在翅翼上表面形成了漩涡，说明出现了非定常的空气流动；由于上下表面空气流速的差别，造成如图9所示的上表面压力小于下表面压力，上下表面压力差便形成了升力。
4.85e+014.65e+014.35e+014.06e+013.83e+013.55e+013.32e+013.07e+012.81e+012.56e+012.35e+012.05e+011.75e+011.53e+011.28e+011.02e+017.67e+005.12e+002.54e+003.92e+00
5.11e+014.85e+014.65e+014.35e+014.06e+013.83e+013.55e+013.32e+013.07e+012.81e+012.56e+012.35e+012.05e+011.75e+011.53e+011.28e+011.02e+017.67e+005.12e+002.54e+003.92e+00
（a）速度分布云图
图8流场速度
1.02e+051.02e+051.02e+051.02e+051.02e+051.02e+051.02e+051.02e+051.01e+051.01e+051.01e+051.01e+051.01e+051.01e+051.01e+051.00e+051.00e+051.00e+051.00e+059.99e+04
1.03e+051.02e+051.02e+051.02e+051.02e+051.02e+051.02e+051.02e+051.02e+051.01e+051.01e+051.01e+051.01e+051.01e+051.01e+051.01e+051.00e+051.00e+051.00e+051.00e+059.99e+04
（b）速度等值线图
［1］孙茂.昆虫飞行的高升力机理［J］.力学进展，2002（32）：425-434.
［2］孙茂，吴江浩.昆虫飞行的高升力机理和能耗［J］.北京航空航天大学
）：970-977.学报，
［3］刘岚，方宗德，侯宇，傅卫平，吴立言.微型扑翼飞行器的气动建模分析
与试验［J］.航空动力学报，）：22-28.
［4］刘岚.微型扑翼飞行器的仿生翼设计技术研究［D］：［博士学位论文］.
西安：西北工业大学，2007.［5］宗光华，贾明，毕树生，徐一村.扑翼式微型飞行器的升力测量与分析
［J］.机械工程学报，）：120-124.［6］江敬强，谢进，陈永.一种仿生微扑翼机构的动力学分析和模拟［J］.机
械设计，）:24-26.
［7］张孝松.基于蜻蜓翅翼的仿生微扑翼飞行器机翼的有限元分析［D］：
［硕士学位论文］.哈尔滨：哈尔滨工业大学，2006.
［8］张正娟，屠恒章，沈怀荣.基于Fluent的低雷诺数下翼型数值模拟
［A］.2007系统仿真技术及其应用学术会议论文集［C］，2007.［9］WootonRJ.Themechanicaldesignofinsectwings，ScientificAmerican，
：114-120.
［10］周建华.微型拍翅式飞行机器人翅运动及控制系统研究［D］：［博士
学位论文］.南京：东南大学，2005.［11］侯宇，方宗德，刘岚，傅卫平.仿生微扑翼飞行器机构动态分析与工程
设计方法［J］.航空学报，）：173-178.［12］金晓怡，颜景平，周建华.翅变形观点分析昆虫飞行时的运动和力［J］.
机械设计，）：40-44.
［13］C.Ellington.Theaerodynamicsofhoveringinsectflight:I.thequasi-steady
analysis.PhilosophicalTransactionsoftheRoyalSocietyofLondonB.
305：1-15.1984，
［14］王大燕.考虑弹簧和运动副间隙的微扑翼飞行器翅翼机构动力学研究
［D］：［硕士学位论文］.成都：西南交通大学，2008.［15］AnsariSA，ZbikowskiR，et.alAerodynamicmodelofinseat-likeflapping
feightformicroairvehicles.ProgressinAerospaceSciences，2006（42）：129-172.［16］韩占忠，王敬，兰小平.FLUENT流体工程仿真计算实例与应用［M］.
北京:北京理工大学出版社，2004.［17］李仁宪.有限体积法基础［M］.北京:国防工业出版社，2005.
（a）压力分布云图（b）压力等值线图
图9流场压力分布
4.3翅翼以其他拍打模式运动升力系数的比较
翅翼拍打过程中升力系数对于微扑翼飞行器翅翼机构的设将上述铰链四杆机构的翅翼机构翅翼计和分析具有重要的意义。
运动的升力系数与翅翼匀速拍打以及按正弦规律拍打的升力系数进行比较。
为了便于比较，将图1所示机构中构件3的摆动范围85°作为翅翼的拍动幅度。当翅翼匀速拍打时，其运动规律为：
≤≤≤≤≤≤≤≤≤
85π×2ft+54π，（0≤t≤0.0125）
-85π×2ft-116π，（0.0125≤t≤0.025）（7）
当翅翼按正弦规律拍打时，其运动规律为：π
准（t）=1×85πsin2πft--23π，（0≤t≤0.025）
式中：f—频率均取40Hz。
采用与上述相同的过程，可以得到升力系数的变化曲线，如图10所示。图中：曲线A—以铰链四杆机构为翅翼机构时的升力系数；曲线B—翅翼匀速拍打的升力系数；曲线C—翅翼按正弦规律运动时的升力系数。从图中可以看出，以铰链四杆机构为翅翼机构获得的升力系数稍大。
第3页/共3页
寻找更多 ""