shell63fluid80直接法流固耦合模態(tài)分析

1、用直接法做流固耦合模態(tài)的公開(kāi)命令流成熟的真的不太多，看來(lái)看去都是那個(gè)fluid30的例子，用直接法fluid80做流固耦合模態(tài)的，目前任何書(shū)本或文獻(xiàn)都沒(méi)有提供一個(gè)簡(jiǎn)單有效的命令流，但是看到文獻(xiàn)有這樣用，我就相信應(yīng)該是可以做得出來(lái)，以下是我一直到剛才所做的命令流，有幾點(diǎn)不確定的地方。命令流下面有我的參考數(shù)據(jù)，是引用別人所述。l.Fluid80的Keyoption(2)=l定要用嗎？或是用0也可以？因?yàn)槲矣?就出錯(cuò)誤，用0就可以跑。但我看過(guò)有人用1。2.耦合的地方我用CPINTF,共兩小段，特別是這里我很不確定。主要是有參考一篇文章。3.主自由度選取的地方有問(wèn)題嗎？/PREP7!定義殼材料與性質(zhì)!

2、殼元素與材料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!以下這個(gè)keyoption怎么用？如過(guò)用1就會(huì)顯示Element877maynothaveapositiveZcoordinateIFKEYOPT(2)=1.，顯示這個(gè)錯(cuò)誤代表要做什么修正嗎?所以我暫時(shí)用KEYOPT(2)=0就可以跑。KEYOPT,2,2,0!建立殼關(guān)鍵點(diǎn)K,1,10,0,0$K,2,10,0

3、,12!建立中心線(xiàn)關(guān)鍵點(diǎn)k,3,0,0,0$k,4,0,0,20!定義殼壁線(xiàn)L,1,2$L,1,3!以關(guān)鍵點(diǎn)3,4為中心線(xiàn)旋轉(zhuǎn)360度生成殼體AROTAT,all,3,4,360!劃分殼體網(wǎng)格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!劃分水位網(wǎng)格type,2$mat,2$esize,2$mshape,0,3D$mshkey,1$vmesh,allalls!以上建模應(yīng)該沒(méi)太大問(wèn)題!以下是耦合，我在流固界面上的網(wǎng)格是重合節(jié)點(diǎn)，特別是下面這

4、兩段落我很不確定該怎么設(shè)定，感覺(jué)問(wèn)題就出在這邊了!這里解決了應(yīng)該就可以。要怎么改?或是用CP?或是NUMMRG?重點(diǎn)是流體和固體要一起動(dòng)，通常設(shè)定不好就流體自己動(dòng)，或是流體都跑到殼體外面去了，流體跟殼不應(yīng)該穿越，而是一起有行為。csys,1!將工作平面定義為柱坐標(biāo)。nsel,s,loc,x,10nrotate,allCPINTF,UX,0.0001,nsel,s,loc,z,0nrotate,allCPINTF,UZ,0.0001,!旋轉(zhuǎn)節(jié)點(diǎn)坐標(biāo)系。!將徑向約束（即X方向）加到節(jié)點(diǎn)上。!邊界條件，將底部固定，并給予Z方向加速度。NSEL,S,LOC,Z,0$D,ALL,ALL$acel,9.8

5、fini/soluantype,modalmodopt,reduc,10,mxpand,10,csys,1!Z上柱坐標(biāo)系!定義主自由度，由圖顯示感覺(jué)是沒(méi)問(wèn)題，但我也不太確定。Esel,s,type,1Nsle,s,allNsel,u,loc,z,0m,all,uxEsel,s,type,2Nsel,s,loc,z,10m,all,uzallssolvefini!選擇殼!所有點(diǎn)!排除邊界條件!（徑向）x方向的主自由度!選擇液體!再選擇液面表面!（豎向）z方向的主自由度主自由度選取做slosh分析，那么全部固體自由度都必須是主自由度。可以不要流體內(nèi)部自由度，接口上耦合他們徑向的位移，主自由度就是固

6、體自由度，相應(yīng)的流體自由度已經(jīng)耦合了。流體的切向自由度，你可以都消掉，這些流體的接口切向運(yùn)動(dòng)無(wú)關(guān)緊要。必須把流體表面的豎向自由度留下來(lái)，其他自由度消掉。選取了液面交界處x方向的主自由度和液面z方向的主自由度，算出來(lái)的結(jié)果和線(xiàn)性波理論的結(jié)果比較吻合。耦合在給定節(jié)點(diǎn)處生成并修改耦合自由度集命令：CP在生成一個(gè)耦合節(jié)點(diǎn)集之后，通過(guò)執(zhí)行一個(gè)另外的耦合操作（保證用相同的參考編號(hào)集）將更多節(jié)點(diǎn)加到耦合集中來(lái)。也可用選擇邏輯來(lái)耦合所選節(jié)點(diǎn)的相應(yīng)自由度。用CP命令輸入負(fù)的節(jié)點(diǎn)號(hào)來(lái)刪除耦合集中的節(jié)點(diǎn)。要修改一耦合自由度集（即增、刪節(jié)點(diǎn)或改變自由度標(biāo)記）可用CPNGEN命令。（不能由GUI直接得到CPNBGEN命令）。耦合重合節(jié)點(diǎn)。CPINTF命令通過(guò)在每對(duì)重合節(jié)點(diǎn)上定義自由度標(biāo)記生成一耦合集而實(shí)現(xiàn)對(duì)模型中重合節(jié)點(diǎn)的耦合。此操作對(duì)“扣緊”幾對(duì)節(jié)點(diǎn)（諸如一條縫處）尤為有用。命令：CPINTF除耦合重復(fù)節(jié)點(diǎn)外，還可用下列替換方法迫使節(jié)點(diǎn)有相同的表現(xiàn)方式：如果對(duì)重復(fù)節(jié)點(diǎn)所有自由度都要進(jìn)行耦合，常用NUMMRG命令合并節(jié)點(diǎn)。可用EINTF命令通在重復(fù)節(jié)點(diǎn)對(duì)之間生成2節(jié)點(diǎn)單元來(lái)連接它們。用CEINTF命令將兩個(gè)有不相似網(wǎng)格模式的區(qū)域連接起來(lái)。這項(xiàng)操作使一個(gè)區(qū)域的選定節(jié)點(diǎn)與另一個(gè)區(qū)域的選定單元連接起來(lái)生成約束方程。NODALSOLirilONSIEP=1SUB=10FRE(=,：331礎(chǔ)5U