shell63-fluid80直接法流固耦合模态分析。(附命令流)
浏览：134289 回答：5
用直接法做流固耦合模态的公开命令流成熟的真的不太多，看来看去都是那个fluid30的例子，用直接法fluid80做流固耦合模态的，目前任何书本或文献都没有提供一个简单有效的命令流，但是看到文献有这样用，我就相信应该是可以做得出来，以下是我一直到刚才所做的命令流，有几点不确定的地方。命令流下面有我的参考数据，是引用别人所述。
1.Fluid80的Keyoption(2)=1一定要用吗? 或是用0也可以? 因为我用1就出错误，用0就可以跑。但我看过有人用1。
2..耦合的地方我用CPINTF，共两小段，特别是这里我很不确定。主要是有参考一篇文章。
3.主自由度选取的地方有问题吗?
有做过的大大给点意见吧! 什么意见都好呀~
-------------------------------------------------------------------------------------------------
/PREP7
!定义壳材料与性质
!壳元素与材料
ET,1,shell63 $MP,EX,1,201E9 $MP,prxy,1,0.26$MP,dens,1,7.85E3 $r,1,0.006
!流体元素与材料
ET,2,FLUID80 $MP,EX,2,1.5e9 $MP,DENS,2,0.84e3$mp,visc,2,1.0e-10
!以下这个keyoption怎么用? 如过用1，就会显示[Element 877 may not have a positive Z coordinateIF KEYOPT(2) = 1.]，显示这个错误代表要做什么修正吗?所以我暂时用KEYOPT(2) = 0就可以跑。
KEYOPT,2,2,0
!建立壳关键点
K,1,10,0,0 $K,2,10,0,12
!建立中心线关键点
k,3,0,0,0 $k,4,0,0,20
!定义壳壁线
L,1,2 $L,1,3
!以关键点3,4为中心线旋转360度生成壳体
AROTAT,all,,,,,,3,4,360
!划分壳体网格
AATT,1,1,1 $esize,2 $mshape,0,3D $mshkey,2 $amesh,all$alls
!延伸出水位体积
VEXT,2,8,2,0,0,10,0,0,0 $vglue,all
csys,1
!划分水位网格
type,2 $mat,2 $esize,2 $mshape,0,3D $mshkey,1$vmesh,all
alls
!以上建模应该没太大问题
!以下是耦合，我在流固界面上的网格是重合节点，特别是下面这两段落我很不确定该怎么设定，感觉问题就出在这边了!这里解决了应该就可以。要怎么改?或是用CP? 或是NUMMRG? 重点是流体和固体要一起动，通常设定不好就流体自己动，或是流体都跑到壳体外面去了，流体跟壳不应该穿越，而是一起有行为。
csys,1 !将工作平面定义为柱坐标。
nsel,s,loc,x,10
nrotate,all !旋转节点坐标系。
CPINTF,UX,0.0001, !将径向约束（即X方向）加到节点上。
nsel,s,loc,z,0
nrotate,all
CPINTF,UZ,0.0001,
!边界条件，将底部固定，并给予Z方向加速度。
NSEL,S,LOC,Z,0 $D,ALL,ALL $acel,,,9.8
fini
/solu
antype,modal
modopt,reduc,10,,
mxpand,10,
csys,1 !Z上柱坐标系
!定义主自由度，由图显示感觉是没问题，但我也不太确定。
Esel,s,type,,1 !选择壳
Nsle,s,all !所有点
Nsel,u,loc,z,0 !排除边界条件
m,all,ux !(径向)x方向的主自由度
Esel,s,type,,2 !选择液体
Nsel,s,loc,z,10 !再选择液面表面
m,all,uz !(竖向)z方向的主自由度
alls
solve
fini
-------------------------------------------------------------------------------------------------
以下是我参考的数据
主自由度选取
1.做slosh分析，那么全部固体自由度都必须是主自由度。可以不要流体内部自由度，接口上耦合他们径向的位移，主自由度就是固体自由度，相应的流体自由度已经耦合了。流体的切向自由度，你可以都消掉，这些流体的接口切向运动无关紧要。必须把流体表面的竖向自由度留下来，其他自由度消掉。
2.选取了液面交界处x方向的主自由度和液面z方向的主自由度，算出来的结果和线性波理论的结果比较吻合。
耦合
1.在给定节点处生成并修改耦合自由度集
命令：CP
GUI: Main Menu>Preprocessor>Coupling /Ceqn>Couple DOFs
在生成一个耦合节点集之后，通过执行一个另外的耦合操作（保证用相同的参考编号集）将更多节点加到耦合集中来。也可用选择逻辑来耦合所选节点的相应自由度。用CP命令输入负的节点号来删除耦合集中的节点。要修改一耦合自由度集（即增、删节点或改变自由度标记）可用CPNGEN命令。（不能由GUI直接得到CPNBGEN命令）。
2.耦合重合节点。
CPINTF命令通过在每对重合节点上定义自由度标记生成一耦合集而实现对模型中重合节点的耦合。此操作对“扣紧”几对节点（诸如一条缝处）尤为有用。
命令：CPINTF
GUI: Main Menu>Preprocessor>Coupling /Ceqn>Coincident Nodes
3.除耦合重复节点外，还可用下列替换方法迫使节点有相同的表现方式：
如果对重复节点所有自由度都要进行耦合，常用NUMMRG命令（GUI：Main Menu>Preprocessor>NumberingCtrls>Merge Items）合并节点。
可用EINTF命令（GUI：Main Menu> Preprocessor>Create> Elements>At Coincid Nd）通在重复节点对之间生成2节点单元来连接它们。
用CEINTF命令（GUI：Main Menu>Preprocessor> Coupling/Ceqn>Adjacent Regions）将两个有不相似网格模式的区域连接起来。这项操作使一个区域的选定节点与另一个区域的选定单元连接起来生成约束方程。
1.Fluid80的Keyoption(2)=1一定要用吗? 或是用0也可以? 因为我用1就出错误，用0就可以跑。但我看过有人用1。
2..耦合的地方我用CPINTF，共两小段，特别是这里我很不确定。主要是有参考一篇文章。
3.主自由度选取的地方有问题吗?
有做过的大大给点意见吧! 什么意见都好呀~
-------------------------------------------------------------------------------------------------
/PREP7
!定义壳材料与性质
!壳元素与材料
ET,1,shell63 $MP,EX,1,201E9 $MP,prxy,1,0.26$MP,dens,1,7.85E3 $r,1,0.006
!流体元素与材料
ET,2,FLUID80 $MP,EX,2,1.5e9 $MP,DENS,2,0.84e3$mp,visc,2,1.0e-10
!以下这个keyoption怎么用? 如过用1，就会显示[Element 877 may not have a positive Z coordinateIF KEYOPT(2) = 1.]，显示这个错误代表要做什么修正吗?所以我暂时用KEYOPT(2) = 0就可以跑。
KEYOPT,2,2,0
!建立壳关键点
K,1,10,0,0 $K,2,10,0,12
!建立中心线关键点
k,3,0,0,0 $k,4,0,0,20
!定义壳壁线
L,1,2 $L,1,3
!以关键点3,4为中心线旋转360度生成壳体
AROTAT,all,,,,,,3,4,360
!划分壳体网格
AATT,1,1,1 $esize,2 $mshape,0,3D $mshkey,2 $amesh,all$alls
!延伸出水位体积
VEXT,2,8,2,0,0,10,0,0,0 $vglue,all
csys,1
!划分水位网格
type,2 $mat,2 $esize,2 $mshape,0,3D $mshkey,1$vmesh,all
alls
!以上建模应该没太大问题
!以下是耦合，我在流固界面上的网格是重合节点，特别是下面这两段落我很不确定该怎么设定，感觉问题就出在这边了!这里解决了应该就可以。要怎么改?或是用CP? 或是NUMMRG? 重点是流体和固体要一起动，通常设定不好就流体自己动，或是流体都跑到壳体外面去了，流体跟壳不应该穿越，而是一起有行为。
csys,1 !将工作平面定义为柱坐标。
nsel,s,loc,x,10
nrotate,all !旋转节点坐标系。
CPINTF,UX,0.0001, !将径向约束（即X方向）加到节点上。
nsel,s,loc,z,0
nrotate,all
CPINTF,UZ,0.0001,
!边界条件，将底部固定，并给予Z方向加速度。
NSEL,S,LOC,Z,0 $D,ALL,ALL $acel,,,9.8
fini
/solu
antype,modal
modopt,reduc,10,,
mxpand,10,
csys,1 !Z上柱坐标系
!定义主自由度，由图显示感觉是没问题，但我也不太确定。
Esel,s,type,,1 !选择壳
Nsle,s,all !所有点
Nsel,u,loc,z,0 !排除边界条件
m,all,ux !(径向)x方向的主自由度
Esel,s,type,,2 !选择液体
Nsel,s,loc,z,10 !再选择液面表面
m,all,uz !(竖向)z方向的主自由度
alls
solve
fini
-------------------------------------------------------------------------------------------------
以下是我参考的数据
主自由度选取
1.做slosh分析，那么全部固体自由度都必须是主自由度。可以不要流体内部自由度，接口上耦合他们径向的位移，主自由度就是固体自由度，相应的流体自由度已经耦合了。流体的切向自由度，你可以都消掉，这些流体的接口切向运动无关紧要。必须把流体表面的竖向自由度留下来，其他自由度消掉。
2.选取了液面交界处x方向的主自由度和液面z方向的主自由度，算出来的结果和线性波理论的结果比较吻合。
耦合
1.在给定节点处生成并修改耦合自由度集
命令：CP
GUI: Main Menu>Preprocessor>Coupling /Ceqn>Couple DOFs
在生成一个耦合节点集之后，通过执行一个另外的耦合操作（保证用相同的参考编号集）将更多节点加到耦合集中来。也可用选择逻辑来耦合所选节点的相应自由度。用CP命令输入负的节点号来删除耦合集中的节点。要修改一耦合自由度集（即增、删节点或改变自由度标记）可用CPNGEN命令。（不能由GUI直接得到CPNBGEN命令）。
2.耦合重合节点。
CPINTF命令通过在每对重合节点上定义自由度标记生成一耦合集而实现对模型中重合节点的耦合。此操作对“扣紧”几对节点（诸如一条缝处）尤为有用。
命令：CPINTF
GUI: Main Menu>Preprocessor>Coupling /Ceqn>Coincident Nodes
3.除耦合重复节点外，还可用下列替换方法迫使节点有相同的表现方式：
如果对重复节点所有自由度都要进行耦合，常用NUMMRG命令（GUI：Main Menu>Preprocessor>NumberingCtrls>Merge Items）合并节点。
可用EINTF命令（GUI：Main Menu> Preprocessor>Create> Elements>At Coincid Nd）通在重复节点对之间生成2节点单元来连接它们。
用CEINTF命令（GUI：Main Menu>Preprocessor> Coupling/Ceqn>Adjacent Regions）将两个有不相似网格模式的区域连接起来。这项操作使一个区域的选定节点与另一个区域的选定单元连接起来生成约束方程。
KEYOPT(2)
Location of gravity springs:
0 --
Place gravity springs on all sides of all elements
1 --
Place gravity springs only on face of elements located on Z = 0.0 plane (elements must not have positive Z coordinates)