请教流固耦合分析方面的几个问题（帮我看看命令流哪有问题）

临时抱佛脚

（1）流固耦合分析中流体单元选用问题
在一篇论文中看到有人处理流固耦合问题时，流体单元选用的是三维容器流体单元FLUID80，但网上或者帮助里的例子,大多用的却是三维声学流体单元FLUID30。不知道这两种单元的主要区别是什么?假如要模拟一个水电站引水管道上的调压塔(相当于一个装着部分水的圆柱形开口容器,不忽略其塔体的厚度),塔内的水体该选择什么单元比较合适?塔内的水体在考虑有速度或者没有速度的情况下，该怎么选单元？
（2）固体和液体交界面的处理问题。
以上面提到的调压塔的问题为例，假如塔体采用SOLID45单元，水体采用FLUID30单元，建模时，容器和液体有各自不同的表面，即在调压塔与水体的交界面上，同一位置有两个面，那么在进行耦合设置的时候，要怎么处理这两个面的耦合问题？设置耦合时有SF和SFA两个命令，我看了些例子，好象SF命令是用来耦合结点或者单元的，SFA是直接耦合面的，不知道是不是这样？那如果要耦合塔和水的交界面，是不是把交界面上的所有NODE选中后，用命令SF，ALL，FSI；或者把交界面上的固体面和液体面上的面都选定后，用SFA，ALL，，FSI这两个方式都能耦合？还有就是看到有例子把流体和固体交界位置的一层流体单元设置成Morphing Fluid的（帮助里有个这样的例子，好象是用物理环境求解耦合问题）。
（3）表面效应单元SURF154的使用问题
还是以上面的调压塔为例，如果要模拟其在地震作用下的动力响应，按照水工结构抗震规范，可以采用韦斯特伽特附加质量公式将水体以附加质量方式加在塔体上模拟水体，水深z处一定面积上的附加质量为：Mp=7/8*ρ*A*SQRT（H0*z），其中H0为塔内的总水深。要实现这一目标，采用SURF154单元时，怎么能设置出这种质量以Z方向的坐标变化的效果？


关于（1）和（2）中的问题，我根据网上的一个例子的命令流（一个开口容器中装满水的模态计算，容器壁采用的是壳单元），将容器的单元从SHELL63该为SOLID45后进行计算，但得不到结果，不知道问题出在那里。不知道有没有这方面的高手给我解答一下，万分感谢！
命令流如下：
/prep7
et,1,45
et,2,30
!选用FLUID30单元，用于流固耦合问题
mp,ex,1,2e10
mp,nuxy,1,0.167
mp,dens,1,2600
mp,dens,2,1000
!定义Acoustics材料来描述流体材料-水
mp,sonc,2,1400
mp,mu,2,
!
block,-3,3,-3,3,0,3
block,-2,2,-2,2,1,3
VSBV,1,2,SEPO,DELETE,KEEP

/PNUM,KP,1

/PNUM,LINE,1
/PNUM,AREA,1
/PNUM,VOLU,1
/PNUM,NODE,0
/PNUM,TABN,0
/PNUM,SVAL,0
/NUMBER,0

KWPAVE,
11
VSBW,
3
wpro,,,90.000000
VSBW,ALL
wpro,,90.000000,
VSBW,ALL

KWPAVE,
10
VSBW,ALL
wpro,,90.000000,
VSBW,ALL

VSEL,S, , ,
2

CM,V-WATER,VOLU
!建立水体的体集合
ALLS
VSEL,U, , ,
2
CM,V-CONTAINER,VOLU
!建立容器的体集合
alls
!划分网格
MSHAPE,0,3D

MSHKEY,1
ESIZE,0.5
cmsel,s,V-CONTAINER
VATT,1,,1
VMESH,ALL

ALLS
CMSEL,S,V-WATER
VATT,2,,2
VMESH,ALL

ASEL,S, , ,
8
NSLA,S,1
CM,N-SURF,NODE
ALLS
ASEL,S,,,9,12,1
ASEL,A,,,7
CM,W-OUHE,AREA
ALLS
ASEL,S,,,18,22,1
CM,C-OUHE,AREA
ALLS
CMGRP,OUHE,W-OUHE,C-OUHE

/solu
antype,2
modopt,unsym,10
!非对称模态提取方法处理流固耦合问题
eqslv,front
mxpand,10,,,1
nsel,s,loc,z,
d,all,,,,,,ux,uy,uz,
CMSEL,S,N-SURF
d,all,pres,0
alls

CMSEL,OUHE
sfa,all,,fsi
!定义流固耦合界面
alls
solv

fini

临时抱佛脚

没有人发表一下见解呀

ruiboy

对命令流一窍不通的新手

lhj西部野狼

不是能看明白楼主最终要解决的是什么问题？

fawwi

你的流体单元应该设置为type1

hahaxucong

费解！期待高手

jog

我觉得流体单元和固体单元的交界面公用节点，模型中节点没有nummrg

一阶模态计算结果如下图：



也不知道对不对

crruose

flotran不能用surf154

liuge711

需要将壁面与液面贴合，就是重合的面GLUE一下，但是好像分析带速度的流体不好实现

kalen001

我也在考虑一个与楼主在物理性质上类似的问题，关于流体和固体界面我在想到底是用接触问题分析还是用流固耦合问题分析，哪个比较好。有高手提提建议吗？