纳闷了好几天了,求达人指导(附图和数据)

evalin

我计算的是一个三维轴对称问题,选择二维rng k-e模型计算,计算模型如下图(下图为位移0和位移10两个模型的叠加)
水流方向向右

                               
登录/注册后可看大图

位移0如图,中间部件每向右移1mm建一个模型,利用fluent积分求解部件前后面的作用力,和压强
数据如下

位移	前作用力/N	后作用力/N	前压力/pa	后压力/pa	备注
0	2.064	-0.784	4846	2391
					未求
2	2.460	-0.863	5776	2921
3	2.714	-0.894	6372	3219
4	2.970	-0.898	6974	3462
5	3.332	-0.923	7823	3821
6	3.802	-0.958	8926	4272
7	4.302	-0.967	10102	4665
8	4.774	-0.929	11211	4889
9	5.207	-0.844	12228	4953
10	-0.103	4.918	-242	-8487

11	-0.106	5.329	-248	-9556
12	-0.111	5.252	-260	-9728

观察这数据,很容易发现位移10的时候,数值突然变负了,后面位移13,14也都是这样
可以保证所以设计,网格,计算过程都是跟前面一样的,而且出错后,还整体重先算了一遍,
结果依旧这样
这个问题困扰了我好几天了
忘哪位达人能指点一下

evalin

忘了说:边界条件是,1m/s速度入口,自由出流

若还需提供其他信息,随时奉上

xrs333

应该是流动固有的特性。严格的说。这种情况阀芯下游的流动是非定常的（即使上游条件是定常的），开度小的时候更明显，所以应该做非稳态计算。这恐怕也是通常阀门开度小的时候产生振动的根本原因。在建立模型之前应该对系统的行为（behavior）作一定的分析。
你愿意把cas和dat文件发到我邮箱吗（9和10两个位置的）：xrs333@yahoo.com.cn，文件压缩后应该不太大。

evalin

给您发过去了
thank you!

evalin

都没人帮忙解答嘛?

周华

哪边是入口，哪边是出口？

xrs333

1. 问题出在参考压力（operating pressure）。
2. 你的模型是没有压力边界的不可压缩流动，也就是说模型里没有压力条件，也不能通过状态方程确定压力，这样压力场是浮动的。
3. 为了避免浮动的压力场，需要指定参考压力位置，即指定该点的压力等于参考压力。缺省参考压力为101325Pa，参考压力位置x=0，y=0。
4. 你#10及其后的模型阀芯移动，使点（0，0）位于计算域外，因此这些模型实际没有参考压力。
5. 解决办法：定义合适的参考压力位置。如有实验数据，将参考压力位置定义在有数据的点，便于比较。
6. 实际参考压力位置与所指定位置会稍有不同，可使用文本命令报告实际参考压力位置：defin / operating-conditions / used - ref - pressure - location。
7. 我在前面说的问题你仍然需要考虑。

注：这里所称的参考压力（operating pressure）与后处理时所用的参考压力（reference pressure）不同，尽管后者翻译符合英文原文，但我仍愿意用“参考压力”一词，更符合其所实际被使用的含义，将“operating pressure”翻译成“运行压力”没有反应出其物理意义。大家意见如何？

[ 本帖最后由 xrs333 于 2009-8-26 17:55 编辑 ]