Shell63-Fluid80水槽,直接法流固耦合.docxVIP

本人最近在做流固耦合，一个三维水槽，槽体用shell63，水用fluid80，想做模态分析，其行为要真的流固耦合，也就是说水槽不是刚体，水槽会因为水的作用而变形，而水也会因为水槽的变形而变形，即为流固耦合。? ? 首先，我在国际期刊及国内论文发表上都有看到shell63-fluid80做模态和瞬态的分析，瞬态大多是以地震力为之，但在中文网络上最常见到是shell63-fluid30的例子，用声学元素fluid30来仿真水，或是利用ANSYS-FLOTRAN做间接耦合，以就是fluid142，但是间接法不管是ANSYS-FLOTRAN还是ANSYS-CFX在我理解范围中是解稳态和瞬态，并无法解模态，所以我想以直接法进行模态分析，而fluid30和fluid80的选择上，由于国际期刊及国内学者的发表都是用fluid80，以及元素性质的判断，用fluid80会比用fluid30来的更仿真。? ???在“Finite Elements in Anslysus and Design 44 (2008)”看到一篇加拿到作者做的论文“Three-dimensional modeling of curved structures containing and/or submerged in fluid”，在5. Numerical results and discussions中对于壳体和流体的几何尺寸、材料有详细的说明，以及阐述了在ANSYS的应用，是很好的验证例子，因为总不能ANSYS做出来图很漂亮就是答案了，数据验证的过程更为重要，此篇文献的矩形水槽有空槽和满水位两种情况，分别求出前五个频率，我在空槽的情况做出来和文献是无误差的，但是在满水位的情况下，一直做不出来，在此我就要请教各位前辈了，如何做出shell63-fluid80的模态分析呢?? ?? ? 以下附件是文献的论文全文和一张数据表，这是我的验证对象，empty是空槽的情况，做出来和文献无误差，fill-water是满槽水位的，一直做不出来，我不知道要弄接触对呢?还是弄标记接口?实在花了时间还是搞不定，refer other people是我参考的一个命令流，当中似乎是用标记接口。? ? 不知能否在此有前辈做过这方面的仿真，如果有的话还请不吝赐教，我每天都用信箱，亦可用信箱跟我联系 我看我直接附上命令流，這樣大家要看就不用載了。這是空槽的情況，也就是empty的命令流/prep7et,1,shell63MP,EX,1,2E11? ?? ??MP,nuxy,1,0.3MP,dens,1,7970r,1,0.005k,1,0,0,0k,2,0.5,0,0k,3,0.5,0.7,0k,4,0,0.7,0k,5,0,0,0.4k,6,0.5,0,0.4k,7,0.5,0.7,0.4k,8,0,0.7,0.4A,1,2,3,4A,5,1,4,8A,6,2,3,7A,5,6,2,1A,8,7,3,4AATT,1,1,1mshape,0mshkey,0esize,0.1amesh,allNSEL,S,LOC,z,0? ??D,ALL,ALLALLSEPLOT/VIEW,1,1,1,1/soluAntype,modalmodopt,lanb,5,mxpand,5,,,0savesolvefinish 以下是我做滿水位的情況，不知諸位大大是否做過，有做過的給個意見。因為一直做不出跟文獻一樣的值，讓我很困惑，我做出來的值都很小。或是不知道的錯誤。懇請了~/prep7et,1,shell63MP,EX,1,2E11 $MP,nuxy,1,0.3 $MP,dens,1,7970 $r,1,0.005ET,2,FLUID80?MP,dens,2,1000 $MP,EX,2,2110E5keyopt,2,1? ?? ?? ?? ?? ?? ???!不知道這對不對!建模與分網k,1,0,0,0 $k,2,0.5,0,0 $k,3,0.5,0.7,0 $k,4,0,0.7,0k,5,0,0,0.4 $k,6,0.5,0,0.4 $k,7,0.5,0.7,0.4 $k,8,0,0.7,0.4A,1,2,3,4 $A,5,1,4,8 $A,6,2,3,7 $A,5,6,2,1 $A,8,7,3,4AATT,1,1,1 $mshape,0,3D $mshkey,2 $esize,0.1 $amesh,allVEXT,1,,,0,0,0.4,0,0,0type,2 $mat,2 $mshape,0,3D $mshkey,1 $esize,0.1 $vmesh,1!---------