|
楼主 |
发表于 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里面设定的边界条件却是入口速度和出口的压力。而且求出的压力分布图是出口各点压力相等,入口的压力各点压力不等。和我期待的结果正好相反,现在就是在这个怎么设定边界条件上面卡住了,请大家帮忙。 |
|