shell63-fluid80直接法流固耦合模态分析

1、用直接法做流固耦合模态的公开命令流成熟的真的不太多，看来看去都是那个 fluid30 的例 子，用直接法 fluid80 做流固耦合模态的，目前任何书本或文献都没有提供一个简单有效的 命令流， 但是看到文献有这样用， 我就相信应该是可以做得出来， 以下是我一直到刚才所做 的命令流，有几点不确定的地方。命令流下面有我的参考数据，是引用别人所述。1. Fluid80 的 Keyoption(2)=1 一定要用吗 ? 或是用 0也可以? 因为我用 1就出错误，用 0就可 以跑。但我看过有人用 1。2.耦合的地方我用 CPINTF ，共两小段，特别是这里我很不确定。主要是有参考一篇文章。 3.主自由度

2、选取的地方有问题吗 ?/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 coordinate IF KEYOPT(2) = 1. ，显示这个错误代表要做什么修正吗 ?

3、所以我暂时用 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,allcsys,1!划分水位网格type,2 $mat,2 $esi

4、ze,2 $mshape,0,3D $mshkey,1 $vmesh,allalls!以上建模应该没太大问题 !以下是耦合，我在流固界面上的网格是重合节点，特别是下面这两段落我很不确定该怎么 设定，感觉问题就出在这边了 ! 这里解决了应该就可以。 要怎么改 ?或是用 CP? 或是 NUMMRG? 重点是流体和固体要一起动，通常设定不好就流体自己动，或是流体都跑到壳体外面去了， 流体跟壳不应该穿越，而是一起有行为。csys,1! 将工作平面定义为柱坐标。nsel,s,loc,x,10! 旋转节点坐标系。! 将径向约束（即 X 方向）加到节点上。nrotate,allCPINTF,UX,0.000

5、1, nsel,s,loc,z,0 nrotate,allCPINTF,UZ,0.0001,!边界条件，将底部固定，并给予Z 方向加速度。NSEL,S,LOC,Z,0 $D,ALL,ALL $acel,9.8 fini /solu antype,modalmodopt,reduc,10, mxpand,10,csys,1 !Z 上柱坐标系!定义主自由度，由图显示感觉是没问题，但我也不太确定。Esel,s,type,1Nsle,s,allNsel,u,loc,z,0 m,all,uxEsel,s,type,2 Nsel,s,loc,z,10 m,all,uz alls solve fini 主自

6、由度选取! 选择壳! 所有点! 排除边界条件!（径向 ）x 方向的主自由度 ! 选择液体 !再选择液面表面!（竖向）Z方向的主自由度1. 做 slosh 分析，那么全部固体自由度都必须是主自由度。可以不要流体内部自由度，接口 上耦合他们径向的位移， 主自由度就是固体自由度， 相应的流体自由度已经耦合了。 流体的 切向自由度， 你可以都消掉， 这些流体的接口切向运动无关紧要。 必须把流体表面的竖向自 由度留下来，其他自由度消掉。2. 选取了液面交界处 x方向的主自由度和液面 z方向的主自由度，算出来的结果和线性波理 论的结果比较吻合。耦合1. 在给定节点处生成并修改耦合自由度集命令： CP在生成

7、一个耦合节点集之后， 通过执行一个另外的耦合操作 （保证用相同的参考编号集） 将 更多节点加到耦合集中来。也可用选择逻辑来耦合所选节点的相应自由度。用CP命令输入负的节点号来删除耦合集中的节点。 要修改一耦合自由度集 （即增、删节点或改变自由度标 记）可用CPNGEN命令。（不能由 GUI直接得到CPNBGEN命令）。2. 耦合重合节点。CPINTF 命令通过在每对重合节点上定义自由度标记生成一耦合集而实现对模型中重合节点 的耦合。此操作对 “扣紧 ”几对节点（诸如一条缝处）尤为有用。命令：CPINTF3. 除耦合重复节点外，还可用下列替换方法迫使节点有相同的表现方式：如果对重复节点所有自由度都要进行耦合，常用NUMMRG命令合并节点。可用EINTF命令通在重复节点对之间生成 2节点单元来连接它们。用CEINTF命令将两个有不相似网格模式的区域连接起来。这项操作使一个区域的选定节点 与另一个区域的选定单元连接起来生成约束方程。NODAL SOLUTIONSTEF=iSUB =10 FRE=.331fi35 USUMR