SUN老师和各位高手帮帮忙呀(fsi例子)

movboy

我把您提供的第一个例子(壳单元耦合及实体单元声固耦合)作了一遍,但运行总不出结果,总显示The degree of freedom solution is not availabe.The PLNSOL command is ignored.,我想应该是Fsi处的问题,但不知道怎么改,希望您能解答.部分内容不会用命令流,只好用gui,谢谢!

/PREP7   
ET,1,FLUID30  
ET,2,FLUID30
ET,3,FLUID30
ET,4,FLUID30
ET,10,SHELL63   
ET,11,SHELL63
KEYOPT,1,2,1
KEYOPT,3,2,1
R,1,2e-5,   
R,2,2e-5,   
R,3,2e-5,   
R,4,2e-5,   
R,5,0.003376
R,6,0.003376
MP,EX,1,3.515e9
MP,PRXY,1,0.35
MP,DENS,1,1350
MP,DENS,2,1.21
MP,SONC,2,344  
MP,DENS,3,1.21  
MP,SONC,3,344  
MP,MU,3,1     

RECTNG,0,0.5,0,0.5,
RECTNG,0,0.5,0,0.5,  

k,9,0,0,5
k,10,0.5,0,5
k,12,0,0.5,5
k,11,0.5,0.5,5
k,13,0,0,-5
k,14,0.5,0,-5
k,16,0,0.5,-5
k,15,0.5,0.5,-5

l,9,10
l,10,11
l,11,12
l,12,9
l,10,2
l,11,3
l,12,4
l,9,1

l,13,14
l,14,15
l,15,16
l,16,13
l,6,14
l,7,15
l,8,16
l,5,13

al,9,10,11,12
al,10,13,2,14
al,11,14,3,15
al,12,15,4,16
al,9,13,1,16

al,17,18,19,20
al,6,21,18,22
al,7,22,19,23
al,8,23,20,24
al,5,21,17,24

va,1,3,4,5,6,7
va,2,8,9,10,11,12

lesize,1,,,5,1
lesize,2,,,5,1
lesize,3,,,5,1
lesize,4,,,5,1
lesize,5,,,5,1
lesize,6,,,5,1
lesize,7,,,5,1
lesize,8,,,5,1
lesize,9,,,5,1
lesize,10,,,5,1
lesize,11,,,5,1
lesize,12,,,5,1
lesize,13,,,10,1
lesize,14,,,10,1
lesize,15,,,10,1
lesize,16,,,10,1
lesize,17,,,5,1
lesize,18,,,5,1
lesize,19,,,5,1
lesize,20,,,5,1
lesize,21,,,10,1
lesize,22,,,10,1
lesize,23,,,10,1
lesize,24,,,10,1

asel,s,,,1
aatt,1,5,10
MSHKEY,1
amesh,1
alls

asel,s,,,2
aatt,1,6,11
MSHKEY,1
amesh,2
alls

!vsel,s,,,1
vatt,2,2,2
MSHKEY,1
vmesh,1
alls

!vsel,s,,,2
vatt,2,4,4
MSHKEY,1
vmesh,2
alls


esel,s,type,,10
nsle,s
esln,s,0
nsle,s
esel,inve
esel,r,type,,2
nsle,s
emodif,all,type,1   
emodif,all,real,1
alls

esel,s,type,,11
nsle,s
esln,s,0
nsle,s
esel,inve
esel,r,type,,4
nsle,s
emodif,all,type,3   
emodif,all,real,3   
alls                  !转化不与结构相连的流体单元


esel,s,type,,10
nsle
esel,s,type,,2
sf,all,fsi
alls

esel,s,type,,11
nsle
esel,s,type,,4
sf,all,fsi
alls                               !定义流固耦合面

nsel,s,loc,z,0
NPLOT     
FLST,5,80,1,ORDE,8  
FITEM,5,1   
FITEM,5,-20
FITEM,5,37  
FITEM,5,-56
FITEM,5,73  
FITEM,5,-92
FITEM,5,505
FITEM,5,-524
NSEL,R, , ,P51X
d,all,uz
alls                            !简支固定边界

esel,s,type,,10,11
nsle
CPINTF,UX,0.0001,
CPINTF,UY,0.0001,
CPINTF,UZ,0.0001,
CPINTF,ROTX,0.0001,
CPINTF,ROTY,0.0001,
CPINTF,ROTZ,0.0001,             !耦合两板
ALLS

nsel,s,loc,z,5
nsel,a,loc,z,-5
esln
emodif,all,mat,3
emodif,all,real,3
sf,all,impd,1            !消声端
alls

esel,s,type,,1
esel,a,type,,3  
nsle
d,all,ux
d,all,uy
d,all,uz                      !固定
alls

/SOLU

antype,harmic
hropt,full  
nsel,s,loc,z,2
d,all,pres,1
alls
HARF,100,100  
SOLVE
FINISH

[ 本帖最后由 movboy 于 2006-12-4 22:02 编辑 ]

北极熊甲

不好意思，不懂声啊
看来我得学习一下声场了

xiejin

可以运行啊. 其实这是网上一个例子
http://ansys.net/ansys/?mycat=se ... mp;mystring=shell63
最下面那个pdf 就是例子的详细介绍. pressure分布是这样

movboy

xiejin兄,是用我的命令流作的图吗,您觉得对吗

真奇怪,刚才我也运行成功了,和你贴的图一样,但再运行就不对了,又出现The degree of freedom solution is not available, The PLNSOL command is ignored,究竟是怎么回事?

[ 本帖最后由 movboy 于 2006-12-5 15:03 编辑 ]

xiejin

我是运行你的文件.
你的问题我不是很清楚,你把当前的working directory里的文件都删掉试试.

[ 本帖最后由 xiejin 于 2006-12-5 16:05 编辑 ]

movboy

我换了个工作目录，不行，然后重装了一下ansys，还是不行
有谁做过这个例子？帮帮忙呀

[ 本帖最后由 movboy 于 2006-12-6 15:43 编辑 ]

xiejin

建议你在第一行/prep7前加上
/CLEAR,START
....
再试试.

movboy

我现在能看到图了，但不是在plots results 里，是在results reviewer里,但我做的好像与例子结果不一样

[ 本帖最后由 movboy 于 2006-12-7 22:22 编辑 ]