|
|
楼主 |
发表于 2007-1-10 14:26
|
显示全部楼层
感谢guxinchun的回复
:hug: :victory: :call: :handshake :kiss: :time: :@) :lol
感谢guxinchun的回复，这里能发附件吗？如果能发，我就把data文件发过来
前两天又去买内存了，现在内存是1.5个G，模拟复杂的结构还是不行。我想把截图和文件传上来大家讨论讨论，却不知道怎么发，谁知道请跟个贴说说，谢谢了。
这两天还在用fluent，昨天弄到半夜，模拟一个栅栏结构。遇到了点问题，请高手帮忙分析一下是什么地方出问题了。
栅栏结构是由78个结构相同的薄片组成。每个薄片的尺寸x*y*z为0.5mm*50mm*210mm（薄片上还均匀分布着贯穿x,y方向,沿z向长度为17mm的沟，每个薄片被三个沟分为四份），78个薄片在空间上沿x轴均匀等距分布，每相邻两个中间缝隙为2mm，所以在空间上结构为78*4个相互独立且平行的实体。
因为gambit划分网格的前提条件就是把实体模型必须合并成1个。所以就想办法，在y方向上的两头分别建立了1mm厚的连接板(尺寸为190mm*1mm*210mm，其中190mm为78个2.5mm积累而成的)，所以加上连接板后，栅栏结构总体的模型变为1个实体，尺寸就变为190mm*52mm*210mm。outlet,inlet分别设定在y方向上两端的侧面。
这个结构比较复杂，为了便于理解，我把在gambit里面建立模型的步骤写一下。
1。solver=fluent5/6,建立单个薄片实体模型 L1=0.5mm*50mm*210mm，
2。切出3个贯穿x,y，z向17mm的沟，建立沟的实体模型L2。然后用布尔运算的减法切除。建立L2=0.5mm*50mm*17mm，在复制两个均匀分布在z向，然后用布尔运算L1-L2-L2copy1-L2copy2
3。沿x向，阵列薄片模型77份，共78份。间距2.5mm。
4。在y向两头加上190mm*1m*210m薄片，整个栅栏结构y方向长度为1+50+1=52mm
5。布尔加法，合并实体。整个变为一个实体
6。划分网格 Elements:Tet/Hybrid 。Type:TGrid。Interval size: 0.5
划分网格，到2个小时左右，提示out of memory。
-----------------------------------------------------
为了减少运算量，而且模型是沿x,y,z都是对称的，所以我沿x,z轴对切，得到原来模型的四分之一研究，然后再mesh完毕后，设定边界的时候，把x向，z向的切除面都设定为对称面就行了，估计应该计算结果和原来一样。
7。简化完毕模型后，划分网格，还是，提示out of memory。
？？？？
问题1：我现在就搞不明白，这个问题该怎么简化才会能够计算，还是我网格单元类型的选择有问题，请高手指教。
然后，我取独立的一个小薄片单元体研究，尺寸0.5*52*210，中间均匀分布三个贯穿x,y,且z向17毫米的矩形通孔。步骤如下：
8。gambit中，划分网格，类型和前面一样，划分成功。
9。设置边界条件，设定y向两端面为inlet和outlet，inlet为velocity inlet，outlet为 pressure outlet，其他面默认为wall
10。导出mesh文件
11。fluent3d打开,
设定入口inlet条件，入口速度=1.7m/s,入口温度293k，湍动能强度10%，水利直径5mm
设定出口outlet条件，出口表压=99671pa,出口温度293k，湍动能强度10%，水利直径5mm
其他设定，如单位mm，couple,implicit,energy equation，k-$2模型，等等。
12。经过初始化，计算等步骤后，到时结果收敛，还可以，但是计算结果中有个问题，如下。
？？？？
问题2：我现在的已知条件是流速和入口压力，求整个留场的压力和流速分布情况。但是现在fluent里面设定的边界条件却是入口速度和出口的压力。而且求出的压力分布图是出口各点压力相等，入口的压力各点压力不等。和我期待的结果正好相反，现在就是在这个怎么设定边界条件上面卡住了，请大家帮忙。 |
|
发表于 2006-12-21 22:22
提升卡
变色卡
显身卡
发表于 2007-1-1 11:20
