三个流固耦合分析实例

length=2!定义体各种变量参数，长宽高width=3height=2/prep7et,1,63 !选用壳模型et,2,30 !选用FLUID30单元，用于流固耦合问题r,1,0.01 增加实常数，壳厚为0.01mp,ex,1,2e11mp,nuxy,1,0.3-水mp,dens,1,7800! -水mp,dens,2,1000! 定义Acoustics材料来描述流体材料mp,sonc,2,1400! 定义声单元声速mp,mu,0,!定义吸声系数!block,,length,,width,,height!建立长方体esize,0.5!选择壳单元mshkey,1!选择壳单元type,1mat,1real,1asel,u,loc,y,width!选择面amesh,all!划分面单元alls !选择所有项type,2!选择声单元mat,2vmesh,all!划分体单元fini/soluantype,2modopt,unsym,10!非对称模态提取方法处理流固耦合问题eqslv,frontmxpand,10,,,1nsel,s,loc,x,nsel,a,loc,x,lengthnsel,r,loc,yd,all,,,,,,ux,uy,uz,nsel,s,loc,y,width,d,all,pres,0!上面几步为定义边界条件和约束allsasel,u,loc,y,width,sfa,all,,fsi!定义流固耦合界面allssolv选择所有项求解fini/post1后处理!显示图形!显示图形/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，就会显示[Element877maynothaveapositiveZcoordinateIFKEYOPT(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,allcsys,1!划分水位网格type,2$mat,2$esize,2$mshape,0,3D$mshkey,1$vmesh,allalls!以上建模应该没太大问题!以下是耦合，我在流固界面上的网格是重合节点，特别是下面这两段落我很不确定该怎么设定，感觉问题就出在这边了!这里解决了应该就可以。要怎么改?或是用CP?或是NUMMRG?重点是流体和固体要一起动，通常设定不好就流体自己动，或是流体都跑到壳体外面去了，流体跟壳不应该穿越，而是一起有行为。csys,1 !将工作平面定义为柱坐标。nsel,s,loc,x,10nrotate,all!旋转节点坐标系。CPINTF,UX,0.0001, !将径向约束（即X方向）加到节点上。nsel,s,loc,z,0nrotate,allCPINTF,UZ,0.0001,!边界条件，将底部固定，并给予 Z方向加速度。NSEL,S,LOC,Z,0$D,ALL,ALL$acel,,,9.8fini/soluantype,modalmodopt,reduc,10,,mxpand,10,csys,1!Z上柱坐标系!定义主自由度，由图显示感觉是没问题，但我也不太确定。Esel,s,type,,1Nsle,s,allNsel,u,loc,z,0m,all,ux!选择壳!所有点!排除边界条件!（径向）X方向的主自由度Esel,s,type,,2Nsel,s,loc,z,10m,all,uz!选择液体!再选择液面表面!（竖向）z方向的主自由度allssolvefini水坝空库/BATCHKEYW,PR_SET,1KEYW,PR_STRUC,1KEYW,PR_THERM,0KEYW,PR_FLUID,0KEYW,PR_ELMAG,0KEYW,MAGNOD,0KEYW,MAGEDG,0KEYW,MAGHFE,0KEYW,MAGELC,0KEYW,PR_MULTI,0KEYW,PR_CFD,0/GO!*/prep7!*definematerialprotertiesmp,dens,1,2650.!mat1fordammp,ex,1,3.15e10mp,prxy,1,.167!*defineelementtypeet,1,PLANE42,,,2!*definegeometrylwater=618.k,1,0.0,0.0,0.0k,2,70.2,0.,0.k,3,0.,66.5,0.k,4,21.9875,66.5,0.k,5,0.,103.,0.k,6,14.8,103.,0.a,1,2,4,3a,3,4,6,5asel,s,loc,x,0.,100.aatt,1,,1

cm,adam,areaallsel,all!*meshgeometryESIZE,0,10lsel,s,loc,y,0.1,66.lesize,all,,,15mshape,0,2Dmshkey,1allsel,allamesh,allfinish/soluantype,modalMODOPT,LANB,30MXPAND,30,,,0esel,s,mat,,1nsle,snsel,r,loc,y,-1.0,1.0d,all,ux,0.d,all,uy,0.!坝下部剖分分数!!坝下部剖分分数!坝体约束/pnum,type,1/number,1gplotallsel,allsavesolve满库/BATCHKEYW,PR_SET,1KEYW,PR_STRUC,1KEYW,PR_THERM,0KEYW,PR_FLUID,0KEYW,PR_ELMAG,0KEYW,MAGNOD,0KEYW,MAGEDG,0KEYW,MAGHFE,0KEYW,MAGELC,0KEYW,PR_MULTI,0KEYW,PR_CFD,0/GO!*/prep7!*definematerialprotertiesmp,dens,1,2650.!mat1fordammp,ex,1,3.15e10mp,prxy,1,.167mp,dens,2,1000.!mat2forwatermp,sonc,2,1440!*defineelementtypeet,1,PLANE42,,,2et,2,29et,3,29,,1!*definegeometrylwater=618.k,1,0.0,0.0,0.0k,2,70.2,0.,0.k,3,0.,66.5,0.k,4,21.9875,66.5,0.k,5,0.,103.,0.k,6,14.8,103.,0.k,7,-1.*lwater,103.,0.k,8,-1.*lwater,66.5,0.k,9,-1.*lwater,0.,0.a,1,2,4,3a,3,4,6,5a,8,3,5,7a,9,1,3,8asel,s,loc,x,0.,100.aatt,1,,1cm,adam,areaasel,s,loc,x,-1*lwater,0.aatt,2,,3cm,awater,areaallsel,all!*meshgeometryESIZE,0,10lsel,s,loc,y,0.1,66. !坝下部剖分分数lesize,all,,,15lsel,s,loc,x,-0.1,-1*lwater-1. !水体长度方向剖分分数lsel,r,loc,y,-0.1,67.lesize,all,,,40,0.5lsel,s,loc,x,-0.1,-1*lwater-1. !水体长度方向剖分分数lsel,r,loc,y,100.,104.lesize,all,,,40,2.0mshape,0,2Dmshkey,1allsel,allamesh,all!更改与水体接触的单元类型esel,s,type,,3nsle,snsel,r,loc,x,-1.,1.esln,remodif,all,type,2allsel,allfinish/soluantype,modalMODOPT,UNSYM,30MXPAND,30,,,0

/pnum,type,1/number,1nsel,s,loc,x,-1.,1.esln,sesel,r,mat,,2sf,all,f