ANSYS流固耦合分析实例命令流

1、ANSYS流固耦合分析实例命令流 DareDoc分享知识 传播快乐分享知识 传播快乐达尔文档 DareDocANSYSS流固耦耦合分析析实例命命令流本资料来源源于网络络，仅供供学习交交流2015年年10月 达尔文文档|DarreDooc 整理目 录TOC o 1-2 h z u HYPERLINK l _Toc433704654ANSYSS流固耦耦合例子子命令流流 PAGEREF _Toc433704654 h 3 HYPERLINK l _Toc433704655ANSYSS流固耦耦合的方方式 PAGEREF _Toc433704655 h 3 HYPERLINK l _Toc4337046

2、56一个流固耦耦合模态态分析的的例子11 PAGEREF _Toc433704656 h 3 HYPERLINK l _Toc433704657一个流固耦耦合模态态分析的的例子22 PAGEREF _Toc433704657 h 4 HYPERLINK l _Toc433704658一个流固耦耦合建模模的例子子 PAGEREF _Toc433704658 h 7 HYPERLINK l _Toc433704659一加筋板在在水中的的模态分分析 PAGEREF _Toc433704659 h 8 HYPERLINK l _Toc433704660一圆环在水水中的模模态分析析 PAGEREF _T

3、oc433704660 h 10 HYPERLINK l _Toc433704661接触分析实实例包含含初始间间隙 PAGEREF _Toc433704661 h 14 HYPERLINK l _Toc433704662耦合小程序序 PAGEREF _Toc433704662 h 19 HYPERLINK l _Toc433704663流固耦合练练习 PAGEREF _Toc433704663 h 21 HYPERLINK l _Toc433704664一个流固耦耦合的例例子 PAGEREF _Toc433704664 h 22 HYPERLINK l _Toc433704665使用物理环环境

4、法进进行流固固耦合的的实例及及讲解 PAGEREF _Toc433704665 h 23 HYPERLINK l _Toc433704666针对液面晃晃动问题题，ANNSYSS/LSS-DYYNA提提供三种种方法 PAGEREF _Toc433704666 h 30 HYPERLINK l _Toc4337046671、流固耦耦合 PAGEREF _Toc433704667 h 30 HYPERLINK l _Toc4337046682、 SPPH算法法 PAGEREF _Toc433704668 h 35 HYPERLINK l _Toc4337046693、 ALLE（接接触算法法） PA

5、GEREF _Toc433704669 h 38 HYPERLINK l _Toc433704670脱硫塔于浆浆液耦合合的分析析 PAGEREF _Toc433704670 h 42 HYPERLINK l _Toc433704671ANSYSS坝-库水流流固耦合合自振特特性的例例子 PAGEREF _Toc433704671 h 43 HYPERLINK l _Toc433704672空库时的IINP文文件 PAGEREF _Toc433704672 h 43 HYPERLINK l _Toc433704673满库时的IINP文文件 PAGEREF _Toc433704673 h 45 HY

6、PERLINK l _Toc433704674计算结果 PAGEREF _Toc433704674 h 48ANSYSS流固耦耦合的方方式一般说来，ANSSYS的流固固耦合主主要有4种方式式：1，seqquenntiaal这需要用户户进行APDDL编程进进行流固固耦合sequeentiia指的是是顺序耦耦合以采用MppCCII为例,你可以以利用ANSSYS和一个个第三方方CFD产品执执行流固固耦合分分析。在在这个方方法中，基基于网格格的平行行代码耦耦合界面面(MppCCII) 将ANSSYS和CFD程序耦耦合起来来。即使使网格上上存在差差别，MpCCCI也能够够实现流流固界面面的数据据转换。A

7、NSSYS CD中包含含有MpCCCI库和一一个相关关实例。关关于该方方法的详详细信息息，参见见ANSSYS Couupleed-FFielld AAnallysiis GGuidde中的Seqquenntiaal CCoupplinn2，FSII soolveer流固耦合的的设置过过程非常常简单，推推荐你使使用这种种方式3，mullti-fieeld sollverr这是FSII soolveer的扩展展，你可可以使用用它实现现流体，结结构，热热，电磁磁等的耦耦合4，直接采采用特殊殊的单元元进行直直接耦合合，耦合合计算直直接发生生在单元元刚度矩矩阵一个流固耦耦合模态态分析的的例子11这是一个

8、流流固耦合合模态分分析的典典型事例例，采用用ANSSYS/MECCHANNICAAL可以完完成。处处理过程程中需要要注意以以下几个个方面的的问题：1、单元的的选择；2、流体材材料模式式；3、流固耦耦合关系系的定义义；4、模态提提取方法法。lengtth=22widthh=3heighht=22/prepp7et,1,63et,2,30 !选用FLUUID330单元，用用于流固固耦合问问题r,1,00.011mp,exx,1,2e111mp,nuuxy,1,00.3mp,deens,1,778000mp,deens,2,110000 !定义Acooustticss材料来来描述流流体材料料-水mp

9、,soonc,2,114000mp,muu,0,!blockk,llenggth,wiidthh,hheigghtesizee,0.5mshkeey,11!type,1mat,11real,1asel,u,lloc,y,wwidtthameshh,alllalls!type,2mat,22vmeshh,alllfini/soluuantyppe,22modoppt,uunsyym,110 !非对称称模态提提取方法法处理流流固耦合合问题eqslvv,frronttmxpannd,110,1nsel,s,lloc,x,nsel,a,lloc,x,llenggthnsel,r,lloc,yd,all

10、l,ux,uy,uz,nsel,s,lloc,y,wwidtth,d,alll,prres,0allsasel,u,lloc,y,wwidtth,sfa,aall,fssi !定义流流固耦合合界面allssolvfini/postt1set,ffirsstplnsool,uu,suum,22,1fini一个流固耦耦合模态态分析的的例子22一实例，水水箱采用用SHEELL663单元，水水箱中的的水采用用FLUUID330单元，以以下即为为整个流流固耦合合模态计计算的命命令流文文件：lengtth=11widthh=0.6heighht=00.8/prepp7et,1,63et,2,30 !选用F

11、LUUID330单元，用用于流固固耦合问问题r,1,00.011mp,exx,1,2e111mp,nuuxy,1,00.3mp,deens,1,778000mp,deens,2,110000 !定义Acooustticss材料来来描述流流体材料料-水mp,soonc,2,114000mp,muu,2,!blockk,llenggth,wiidthh,hheigghtesizee,0.1mshkeey,11!type,1mat,11real,1asel,u,lloc,y,wwidtthameshh,alllalls!type,2mat,22vmeshh,alllfini/soluuantyppe

12、,22modoppt,uunsyym,110 !非对称称模态提提取方法法处理流流固耦合合问题eqslvv,frronttmxpannd,110,1nsel,s,lloc,x,nsel,a,lloc,x,llenggthnsel,r,lloc,yd,alll,ux,uy,uz,nsel,s,lloc,y,wwidtth,d,alll,prres,0allsasel,u,lloc,y,wwidtth,sfa,aall,fssi !定义流流固耦合合界面allssolvfini/postt1set,ffirsstplnsool,u,ssum,2,11fini我得出的结结果是：SET TTIMEE/FR

13、REQ LOAAD SSTEPP SUUBSTTEP CUMMULAATIVVE1 29.1955 1 1 112 0.000000 1 1 113 33.2022 1 2 224 0.000000 1 2 225 37.5988 1 3 336 0.000000 1 3 337 44.5922 1 4 448 0.000000 1 4 449 47.9177 1 5 5510 0.00000 11 5 511 888.0115 11 6 612 0.00000 11 6 613 888.5665 11 7 714 0.00000 11 7 715 977.1333 11 8 816 0.00

14、000 11 8 817 1009.449 11 9 918 0.00000 11 9 919 1009.881 11 100 10020 0.00000 11 100 100为了比较我我还计算算了水箱箱不装水水的情况况，（水水箱体模模型的建建立和约约束与上上面的相相同，只只是把流流体部分分去掉），计计算结果果是：SET TTIMEE/FRREQ LOAAD SSTEPP SUUBSTTEP CUMMULAATIVVE1 45.1088 1 1 112 52.5911 1 2 223 68.9222 1 3 334 76.3733 1 4 445 91.3988 1 5 556 1266.49

15、9 1 6 667 1400.433 1 7 778 1522.155 1 8 889 1566.777 1 9 9910 1660.779 11 100 100我想请问：1。水箱装装水时，怎怎么计算算结果中中有0频率啊啊？其中中不等于于0的频率率是水箱箱体的频频率么？2。水箱不不装水时时的结果果和装水水时的相相差怎么么这么大大啊？还还是我理理解错了了么？企盼指教一个流固耦耦合建模模的例子子/prepp7 !进行预预处理模模块et,1,30, !定义1号单元元为Fluuid330 流固耦耦合单元元et,2,29 !定义2号单元元为Fluuid229平面流流体单元元et,3,30 , ,1 !定

16、义3号单元元为Fluuid330流体介介质单元元et,4,63 !定义4号单元元为Sheell663壳体单单元et,5,1888 !定义5号单元元为Beaam1888梁单元元r,4,00.0002 !定义4号单元元的厚度度为2mp,deens,4,778000 !定义4号物理理属性包包括有密密度mp,exx,4,2.11e111 !杨氏模模量、mp,nuuxy,4,00.3 !泊松比比mp,soonc,1,114600 !设置水水中声速速mp,deens,1,110000 !设置流流体密度度sectyype, 1, beeam, T, ! 选取T型梁secofffseet,oriig !设置梁

17、梁的方向向secdaata,0.004,00.055,0.0022,0.02,0,00,0,0,00,0 !所建立的的圆柱壳壳体的参参数：圆圆柱长为为 500 ,半径为 255 ,壳体的的壁厚为为 2 , cyl4,0,00,0.25,5 !形成圆圆面k,9,00,0,0 ! 定义原原点k,10,0,00,0.5lstr,9,110 !通过原原点作直直线adragg,5,6,77,8,99 !通过放放样形成成圆柱wpofff,0,0,00.1asel,s,2,5asbw,alll, !移动工工作平面面与选取取的侧面面相切!重复上面面操作,形成四四个环肋肋面wpofff,0,0,-0.44 !工

18、作平平面回到到原点位位置上k,31,0.22,0,0.11 !定义环环肋的方方向点lsel,s,200 !选择要要划分为为环肋的的线段latt,4,55,5,311,400,1 !定义线线段物理理属性lesizze,220,6 !划分数数目secnuum,11lmeshh,200 !划分线线段type，2 ! 选取第第二种单单元lsel，s，1，4 ！选取取线段lesizze，all，10 ！线段段划分数数目lesizze，all，6ameshh，1 ！将通通过工作作平面的的面 1 进行划划分esizee，0.1，0 ！选取取沿面放放样的网网格大小小mat，11 ！定义义放样的的形成单单元的物

19、物理属性性vdragg，1，9 ！进行放样形成流体介质 一加筋板在在水中的的模态分分析再给大家一一个实例例！考虑结构在在水中的的自振频频率：例例子是一一加筋板板在水中中的模态态分析。命令流如下下：FINISSH/CLEAAR/FILEENAMME,pplanne/UNITTS,SSI/TITLLE,pplanne/PREPP7!*EELEMMENTT DEEFINNE*ET,633,633ET,4,beaam4et,300,flluidd30!*MATTERIIAL DEFFINEE*MP,EXX,1,2.110E111MP,DEENS,1,778500MP,NUUXY,1,00.3mp,de

20、ens,30,10225mp,soonc,30,15000mp,muu,300,0.5!*REAAL CCONSSTANNT*r,30,1e-06r,50,0.005r,75,0.3375ee-022,0.781125ee-066,0.000001664066k,1k,4,11kfilll,1,4,22,11kgen,4,11,4,1,1/33,110a,1,22,122,111*do,ii,0,2*do,jj,0,2*110,110a,1+ii+j,2+ii+j,12+i+jj,111+i+j*endddo*endddo!*ffluiid eelemmentt*k,1000,-114.55,

21、-114.55k,1011,-114.55,155.5k,1022,155.5,15.5k,1033,155.5,-144.5k,1400,-114.55,-114.55,300k,1411,-114.55,155.5,30k,1422,155.5,15.5,330k,1433,155.5,-144.5,30a,1000,1001,1102,1033,4,14,24,34,33,32,31,21,11,1a,1,22,3,4,1103,1000a,1400,1441,1142,1433a,1000,1001,1141,1400a,1011,1002,1142,1411a,1422,1443,1

22、103,1022a,1400,1443,1103,1000a,14,24,34,33,32,31,21,11,1,22,3,4asel,u,1,FLST,2,88,5,ORDDE,22FITEMM,2,10FITEMM,2,-177VA,P551Xnummrrg,aallallsMSHKEEY,00MSHAPPE,00esizee,1lsel,s,lloc,y,11/3lsel,r,lloc,x,00,1lsel,r,lloc,z,00latt,1,775,44lmeshh,allllsel,s,lloc,y,22/3lsel,r,lloc,x,00,1lsel,r,lloc,z,00latt

23、,1,775,44lmeshh,allllsel,s,lloc,x,11/3lsel,r,lloc,y,00,1lsel,r,lloc,z,00latt,1,775,44lmeshh,allllsel,s,lloc,x,22/3lsel,r,lloc,y,00,1lsel,r,lloc,z,00latt,1,775,44lmeshh,alllasel,s,1,9aatt,1,550,663ameshh,alllallsMSHAPPE,11,3ddesizee,3vsel,s,1type,30 $maat,330 $reaal,330vmeshh,alllallsFINISSH/soluuall

24、s!* 求解 *!*ANTYPPE,MMODAAL!MODOOPT,lannb,225,00MODOPPT,UUNSYYMM,25,0SOLVEEFINISSH总是出现eerroor 说矩阵阵不对称称，不可可以用lannb计算。总结：流体体单元不不能用对对称的解解法应该采用非非对称解解法。一圆环在水水中的模模态分析析finissh/cleaar/PREPP7!定义单元元类型ET,1,PLAANE442 ! sttruccturral eleemenntET,2,FLUUID229 ! accoussticc flluidd ellemeent witth uux & uyyET,3,1299

25、! acoousttic inffiniite linne eelemmenttr,3,00.3112422,0,0ET,4,FLUUID229,1,00 ! acoousttic fluuid eleemennt wwithhoutt uxx & uy!材料属性性MP,EXX,1,2.0068ee11MP,DEENS,1,779299MP,NUUXY,1,00MP,DEENS,2,110300MP,SOONC,2,114600! 创建四四分之一一模型CYL4,0,00,0.2544,0,0.2260335,990CYL4,0,00,0.260035,0,00.3112422,900! 选择属

26、属性,网格划划分ASEL,S,AAREAA,11AATT,1,11,1,0LESIZZE,11,16,1LESIZZE,33,16,1LESIZZE,22,1,11LESIZZE,44,1,11MSHKEEY,11MSHAPPE,00,2DD ! mapppedd quuad messhAMESHH,1ASEL,S,AAREAA,22AATT,2,11,2,0LESIZZE,55,16,1LESIZZE,77,16,1LESIZZE,66,5LESIZZE,88,5MSHKEEY,00MSHAPPE,00,2DD ! mapppedd quuad messhAMESHH,2! 关于YY轴镜像像

27、nsym,x,110000,alll ! offfseet nnodee nuumbeer bby 110000esym,10000,alll! 关于yy轴镜像像nsym,y,220000,alll ! offfseet nnodee nuumbeer bby 220000esym,20000,alllNUMMRRG,AALL ! mmergge aall quaantiitieesesel,s,ttypee,11nsle,sesln,s,00nsle,sesel,invvensle,semodiif,aall,typpe,44esel,alllnsel,alll! 指定无无限吸收收边界csy

28、s,1nsel,s,lloc,x,00.3112422type,3real,3mat,22esurffesel,alllnsel,alll! 标识流流固交接接面nsel,s,lloc,x,00.2660355esel,s,ttypee,22sf,alll,ffsi,1nsel,alllesel,alllFINISSH/soluuantyppe,mmodaalmodoppt,ddampp,100mxpannd,110,yeessolveefinissh为了便于对对比，也也对圆环环在空气气中做了了模态分分析finissh/cleaar/PREPP7!定义单元元类型ET,1,PLAANE442 !

29、sttruccturral eleemennt!材料属性性MP,EXX,1,2.0068ee11MP,DEENS,1,779299MP,NUUXY,1,00! 创建四四分之一一模型CYL4,0,00,0.2544,0,0.2260335,990! 选择属属性,网格划划分ASEL,S,AAREAA,11AATT,1,11,1,0LESIZZE,11,16,1LESIZZE,33,16,1LESIZZE,22,1,11LESIZZE,44,1,11MSHKEEY,11MSHAPPE,00,2DD ! mapppedd quuad messhAMESHH,1! 关于YY轴镜像像nsym,x,1100

30、00,alll ! offfseet nnodee nuumbeer bby 110000esym,10000,alll! 关于yy轴镜像像nsym,y,220000,alll ! offfseet nnodee nuumbeer bby 220000esym,20000,alllNUMMRRG,AALL/soluuantyppe,mmodaalmodoppt,llanbb,100mxpannd,110,yeessolveefinissh在水中的自自振频率率为 SEET TTIMEE/FRREQ LOAAD SSTEPP SUUBSTTEP CUMMULAATIVVE 1-00.1995444

31、E-110 11 1 1 2 00.2996400E-003 11 1 1 3-00.2116633E-110 11 2 2 4-00.2996400E-003 11 2 2 5 00.3008700E-003 11 3 3 6 00.00000 1 33 3 7-00.3008700E-003 11 4 4 8 00.00000 1 44 4 9-00.5337266E-003 11 5 5 110 00.5775222E-111 11 5 5 111 00.5337266E-003 11 6 6 112-00.8990577E-111 11 6 6 113 00.9880599E-001

32、11 7 7 114 335.2232 1 77 7 115 00.9880599E-001 11 8 8 116 -35.2322 1 8 88 117 00.9880611E-001 11 9 9 118 335.2233 1 99 9 119 00.9880611E-001 11 100 100 220 -35.2333 1 10 10在空气中的的自振频频率为 SETT TIIME/FREEQ LLOADD STTEP SUBBSTEEP CCUMUULATTIVEE 1 00.00000 1 11 1 2 00.00000 1 22 2 3 00.7336099E-003 11 3 3

33、 4 660.8805 1 44 4 5 660.8805 1 55 5 6 1172.97 1 66 6 7 1172.97 1 77 7 8 3334.40 1 88 8 9 3334.40 1 99 9 110 5546.59 1 110 110主要有以下下疑问：1）考虑流流固耦合合，做模模态分析析时流体体单元是是否只能能用fluuid229（2d）和fluuid330（3d),对于fluuid1129和fluuid1130在耦合合中具体体起到什什么作用用，能不不能不设设，而用用边界约约束条件件代替？2）流体范范围怎样样确定，如如本例中中（CYLL4,00,0,0.2260335,00,

34、0.312242,90），外外半径为为0.3312442。如果果不是环环形的，如如一块当当水板，该该怎样考考虑？3）如果不不考虑流流体的压压缩性，把把声速设设的很大大，MP,SONNC,22,1ee20，就可可以了。4）从自振振频率可可以看出出，在水水中和在在空气中中，圆环环的自振振频率差差别特别别大，且且振型也也大相径径庭，为为什么？在水中中时，模模态提取取方法用用dammp（为什什么不能能用unssym），特特征值的的虚部代代表角频频率，为为什么第第一阶为为正，第第二阶为为负，而而第三阶阶和第四四阶都为为0，第六六阶、八八阶、十十阶都为为负。应应该是从从小到大大才对？5）在空气气中时，模模

35、态提取取方法用用lannb，为什什么第一一阶第二二阶的频频率都为为0。请高手指点点迷津，急急盼中对以上问题题的解答答：频率为零，一一般是发发生了刚刚体位移移，估计计你是把把水抽走走，而没没有限制制圆环。1。圆环在在水中振振动必然然导致波波动(其实就就是声波波)在水中中传播，当当声波到到达水的的另一个个界面时时就会发发生反射射(除非水水和另一一个相邻邻体的声声阻抗是是匹配的的)。水和和金属中中的声速速相差不不大，即即可压缩缩性相差差不大。两两种可压压缩性相相差不大大的物质质的相互互作用对对两者影影响都很很大。圆圆环在水水中振动动，水对对圆环的的反作用用是由于于反射波波引起的的，流固固耦合中中采用

36、fluuid1129和130就是最最大程度度的减弱弱反射波波。2。声波从从圆环开开始传播播，随着着传播距距离的增增加，波波阵面不不断增大大，振幅幅不断减减小。同同时由于于水的衰衰减，声声波也不不断减弱弱。如果果水的空空间越大大，则反反射波返返回圆环环的路径径越长，衰衰减也就就越多，影影响也就就越小。fluuid1129和130对反射射波的衰衰减(通过很很小的反反射实现现)有限，因因此还需需要水要要有足够够的空间间。fluuid1129和130离结构构应该大大于0.2(=c/f,cc为水中中声速)。以上上的做法法在误差差允许的的情况下下等效于于水在无无限大水水空间中中的情况况。如果果是挡水水板，

37、水水就是有有限空间间了，情情况也不不一样。3。声速加加大情况况也不一一样，就就是不知知是不是是你所要要的情况况？4。空气作作为介质质，由于于其声速速比金属属小很多多，可压压缩性大大很多，影影响可以以忽略不不计。而而水的影影响就不不同了。这这可能就就是频率率和振型型不同的的原因吧吧？我试试了你的的例子，各各种提取取方法都都可以。5。空气的的影响忽忽略不计计，因此此需要对对圆环进进行约束束。你没没有约束束，那么么就会发发生静态态位移即即频率为为零。圆圆环有两两个对称称轴，因因此会发发生频率率成对出出现的情情况。也也就是说说，两个个方向上上有同样样的振型型。接触分析实实例包含初初始间隙隙fini/c

38、leaar,nnosttartt/prepp7et,1,82KEYOPPT,11,3,3r,1,00.5mp,exx,1,1e99mp,prrxy,1,00.3k,1,00,0k,2,110,00k,3,110,55k,4,66.2,5k,5,77.5,3.44k,6,22.5,3.44k,7,33.8,5k,8,00,5a,1,22,3,4,55,6,7,88LFILLLT,66,5,0.118, , LFILLLT,55,4,0.118, , FLST,2,33,4 FITEMM,2,9 FITEMM,2,11 FITEMM,2,10 AL,P551XFLST,2,33,4 FITEMM,

39、2,13 FITEMM,2,14 FITEMM,2,12 AL,P551XFLST,2,33,5,ORDDE,22 FITEMM,2,1 FITEMM,2,-3 AADD,P511Xrect,0,110,44.8,5ASBA, 4, 1 gap=00.022k,24,6.22-gaap,55k,25,7.55-gaap,33.4k,26,2.55+gaap,33.4k,27,3.88+gaap,55a,24,25,26,27LFILLLT,44,3,0.22, , LFILLLT,33,2,0.22, , FLST,2,33,4 FITEMM,2,7 FITEMM,2,10 FITEMM,2,

40、8 AL,P551XFLST,2,33,4 FITEMM,2,13 FITEMM,2,14 FITEMM,2,11 AL,P551XFLST,3,22,5,ORDDE,22 FITEMM,3,3 FITEMM,3,-4 ASBA, 1,P511X rect,3.88+gaap,66.2-gapp,5,10rect,3.88+gaap,33.8+gapp+8,10,12FLST,2,33,5,ORDDE,33 FITEMM,2,1 FITEMM,2,3 FITEMM,2,5 AADD,P511X rect,3.88+gaap+88,3.8+ggap+8+22,100,122FLST,2,22,5

41、,ORDDE,22 FITEMM,2,1 FITEMM,2,4 AGLUEE,P551X CYL4,2.00,1.8,00.6 CYL4,7.00,1.8,00.6FLST,2,33,5,ORDDE,33 FITEMM,2,2 FITEMM,2,4 FITEMM,2,-5 AOVLAAP,PP51XXesizee,0.2ameshh,alllFLST,5,1135,2,OORDEE,322FITEMM,5,4855FITEMM,5,5766FITEMM,5,-5777FITEMM,5,6211FITEMM,5,-6225FITEMM,5,7077FITEMM,5,-7111FITEMM,5,7

42、166FITEMM,5,7411FITEMM,5,-7445FITEMM,5,7500FITEMM,5,-7551FITEMM,5,7666FITEMM,5,7977FITEMM,5,-7998FITEMM,5,8544FITEMM,5,8888FITEMM,5,-9338FITEMM,5,11001FITEMM,5,11003FITEMM,5,14220FITEMM,5,16228FITEMM,5,16553FITEMM,5,16996FITEMM,5,16999FITEMM,5,-17702 FITEMM,5,17226FITEMM,5,-17728 FITEMM,5,18552FITEM

43、M,5,-18874 FITEMM,5,20444FITEMM,5,-20066 CM,_YY,ELLEM ESEL, , , ,P511XCM,_YY1,EELEMMCMSELL,S,_Y CMDELLE,_Y EREF,_Y11, , ,11,0,1,11CMDELLE,_Y1 ET,2,TARRGE1169 ET,3,CONNTA1172 R,3, , ,R,3,00,0,0.11, 110,00,0R,4, , ,R,4,00,0,0.11, 110,-0.002,00 lsel,s,9lsel,a,5lsel,a,122nsll,s,11type,3real,3esurff,alll

44、alls,lsel,s,199lsel,a,200nsll,s,11type,3real,4esurff,alllalls,lsel,s,7lsel,a,3lsel,a,111nsll,s,11type,2real,3esurff,alllalls,lsel,s,255lsel,a,266nsll,s,11type,2real,4esurff,alllalls,FLST,2,22,5,ORDDE,22 FITEMM,2,4 FITEMM,2,-5 DA,P551X,ALLL,FLST,2,11,4,ORDDE,11 FITEMM,2,6 SFL,PP51XX,PRRES,5000, /solu

45、uantyppe,00nlgeoom,oonoutrees,aall,alllnsubsst,2200,2000,2neqitt,10000solvee耦合小程序序最近用到耦耦合,写了一一段小程程序,奉献出出来,与大家家共享。如果有很多多节点，每每两个节节点位置置相同，如如果将这这些杂乱乱无章的的节点耦合，是件件很麻烦烦的事，可可用这段段程序，轻轻松解决决。cpnumm=0cmsell,s,n-zzhonng !需要耦耦合的节节点*GET,n_nnum,NODDE,COUUNT, , , , !节点总总数*do,ii,1,n_nnumcmsell,s,n-zzhonng*GET,n_llowe

46、est,NODDE,NUMM,MIIN, , , , !号码最最小的节节点*GET,n_xx,NOODE,n_lloweest,LOCC,X !该节点点坐标*GET,n_yy,NOODE,n_lloweest,LOCC,Y*GET,n_zz,NOODE,n_lloweest,LOCC,ZNSEL,s,LLOC,X,nn_x-0.33,n_x+00.3 !寻找与与该节点点位置相相同的节节点NSEL,R,LLOC,Y,nn_y-0.33,n_y+00.3NSEL,R,LLOC,z,nn_z-0.33,n_z+00.3cm,n_cp_cp,nodde !位置相相同的节节点形成成一个组组cmsell,

47、s,n-zzhonngcmsell,u,n_ccp_ccpcm,n-zhoong,nodde !取消这这些点后后剩余的的点形成成组*GET,n_nnum_1,NNODEE,CCOUNNT, , , , !节点总总数*if,nn_nuum_11,ltt,2,exiit !如果节节点数小小于二则则退出cmsell,s,n_ccp_ccp*GET,n_nnum,NODDE,COUUNT, , , ,*if,nn_nuum,ggt,11,thhenCP,cppnumm+1,ux,alllCP,cppnumm+2,uy,alllCP,cppnumm+3,uz,alllcpnumm=cppnumm+3*e

48、lsee*endiif*endddo该段程序可可用CPINTTF，UX，0.0001CPINTTF，UY，0.0001CPINTTF，UZ，0.0001代替*DO,II,2,2966,3CP,I,UX,I,II+2*ENDDDO*DO,II,2,2966,3CP,I,UY,I,II+2*ENDDDO*DO,II,2,2966,3CP,I,UZ,I,II+2*ENDDDODK,1, , , ,0,UUX,UUY,UUZ, , , ,以上几句改改为:*DO,II,2,2966,3CP,NEEXT,ALLL,I,I+22*ENDDDODK,1, , , ,0,AALL或CPINNTF,ALLL,0.

49、0011因为你选用用的单元元有六个个自由度度,如果只只约束三三个,程序是是不会运运行的.另:三次循循环语句句的I相等,约束UY时,UX的耦合合就被删删掉了,最后只只剩UZ了这样修改：流固耦合练练习/PREPP7K,1,00,0K,2,00.1,0L,1,22K,3000,0,-1000000LGEN,1000,1,00.1,22ET,1,BEAAM1888MP,EXX,1,2.11e111MP,PRRXY,1,00.3MP,DEENS,1,00.7883e44SECTYYPE, 1, BEEAM, T, , 0SECOFFFSEET, CENNTSECDAATA,0.006,00.033,0.

50、0033,0.0066,0,0,00,0,0,00LSEL,ALLLLATT,1,11,1,3000LESIZZE,AALL,11,11,1LMESHH,ALLLcpinttf,aallDK,1,ux,0, , ,UY,UZDK,2000, , , , ,UYY,UZZACEL,0,99.8,0,FINISSH一个流固耦耦合的例例子这个例子关关于装有有水的水水杯旋转转，是轴轴对称问问题，为为了简化化，所以以选择了了平面模模型。*SET,RADD,0.8*SET,h,11 *SET,g,99.8*SET,OMEEGARR,2 *SET,ROUU,10000 /PREPP7ET,1,FLUUID7

51、79KEYOPPT,11,3,1MP,EXX,1,2E99MP,DEENS,1,RROU K,1K,2,RRADK,3,RRAD,H K,4,HK,4,HA,1,22,3,4 LESIZZE,AALL,110AMESHH,ALLL FINISSH/SOLDL,2,UXXDL,1,UYYNSEL,S,LLOC,XDSYM,SYMMM,XXD,ALLL,UXXD,ALLL,UXXNSEL,ALLLACEL,GOMEGAA,OOMEGGAR SOLVEE FINISSH/POSTT1SET,LLASTTPLNSOOL,UU,X,0,11*SET,UCEENT,UY(22) *SET,UEDDGE,U

52、Y(12) *SET,UELLEV,UEDDGE-UCEENT本人抛砖引引玉：使用物理环环境法进进行流固固耦合的的实例及及讲解ansyss从9.0发展到10.0，一个个最大的的进步就就是流固固耦合计计算更加加规范，这这一点已已远领先先于其他他同类软软件，实实现了单单向耦合合到即时时双向耦耦合的飞飞跃，使使用户对对于解决决流固耦耦合问题题又多了了一种选选择，希希望大家家对多种种方法物理环环境转换换，fsi，mfx等进行行讨论，提提供一下下案例流道中有一一橡胶垫垫阻碍水水的流动动，入口口速度为为2m/s，其他他参数将将在命令令流中详详细给出出。求解解水通过过此流道道的压力力降，以以及稳态态条件下下

53、橡胶垫垫的变形形。/prepp7/sho,gasskett,grrphshpp,offfET,1,1411 ! Flluidd - staaticc meeshET,2,56, ! Hyyperrelaastiic eelemmentt! Fluuid Strructturee Innterracttionn - Mulltipphyssicss! Deeforrmattionn off a gasskett inn a floow ffielld.! Eleemennt pplotts aare wriitteen tto tthe fille ggaskket.grpph.! - Waat

54、err fllowss inn a verrticcal pippe tthrooughh a connstrructtionn! fformmed by a rrubbber gasskett.! - Deeterrminne tthe equuiliibriium possitiion of thee gaaskeet aand! tthe ressulttingg fllow fieeld! | |! | |! | Boounddaryy off mmorpphinng ffluiid! | _|! | |_ gasskett! | |! | Boounddaryy off mmorpp

55、hinng ffluiid (sff)! | |! 1. Buuildd thhe mmodeel oof tthe enttiree doomaiin:! FFluiid rregiion - sstattic messh! ! Gasskett leeavees aa hoole in thee ceenteer oof tthe ducct! Morrphiing Fluuid reggionn iss a useer ddefiinedd reegioon aarouund! tthe gasskett. TThe fluuid messh hheree wiill defformm

56、annd bbe! uupdaatedd ass thhe ggaskket defformms.! Parrameeterrizee Diimennsioons in thee fllow dirrecttionn!yent = 00.0 ! Y ccoorrdinnatee off thhe eentrrancce tto tthe pippedyen = 11.0 ! Unddefoormeed ggeommetrry ffloww enntraancee leengtthysf1 = yyentt+dyyen ! YY cooorddinaate of enttrannce to th

57、ee moorphhingg flluidd reegioondsf1 = 00.5 ! Thiicknnesss off uppstrreammygas = yysf11+dssf1 ! YY cooorddinaate of thee boottoom oof tthe gasskettdg = 0.002 ! TThicckneess of thee gaaskeetdg2=ddg/22.ytg = yggas+dg ! Y cooordiinatte oof tthe iniitiaal ttop of thee gaaskeetdsf2 = 00.5 ! Thiicknnesss o

58、ff doownsstreeam reggionnysf2 = yytg + ddsf22! YY off Toop oof tthe dowwnsttreaam mmorpphinng ffluiids reggionndyex = 66.0 ! Exiit ffluiid llenggthx = 00. ! Loocattionn off thhe aaxissymmmetrric Cennterrlinnedgasrr =.20 ! Iniitiaal sspann off gaaskeetpiperr = 0.33 ! Raadiuus oof tthe pippexrgapp = p

59、ipper-dgaasr! rradiius of commpleetelly uunobbtruucteed ffloww paassaage! CCreaate geoomettry!rect,xrggap,pipper,ygaas,yytg ! A11:Gaaskeet (keyypoiintss 1-4)rect,x,ppipeer,yysf11,yssf2 ! A2: Moorphhingg flluidd reegioonrect,x,ppipeer,yyentt,yssf1 ! A3: Flluidd reegioon wwithh sttatiic mmeshhrect,x,p

60、pipeer,yysf22,yssf2+dyeex ! A4: Flluidd reegioon wwithh sttatiic mmeshhaovlaap,aallk,22,xrggap+dg22,yggas+dg22rarc = ddg2*1.11larc,1,44,222,raarcal,6,4adeleete,7al,6,3,222,77,8,5,221,11!Messh DDiviisioon iinfoormaatioonngap = 110 ! Nuumbeer eelemmentts aacrooss thee gaapngas = 110 ! Nuumbeer oof eel