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

1、ANSYS流固耦合分析實例命令流 DareDoc分享知識 傳播快樂分享知識 傳播快樂達爾文檔 DareDocANSYSS流固耦耦合分析析實例命命令流本資料來源源于網(wǎng)絡(luò)絡(luò)，僅供供學(xué)習(xí)交交流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一個流固耦耦合模態(tài)態(tài)分析的的例子11 PAGEREF _Toc433704656 h 3 HYPERLINK l _Toc433704657一個流固耦耦合模態(tài)態(tài)分析的的例子22 PAGEREF _Toc433704657 h 4 HYPERLINK l _Toc433704658一個流固耦耦合建模模的例子子 PAGEREF _Toc433704658 h 7 HYPERLINK l _Toc433704659一加筋板在在水中的的模態(tài)分分析 PAGEREF _Toc433704659 h 8 HYPERLINK l _Toc433704660一圓環(huán)在水水中的模模態(tài)分析析 PAGEREF _T

3、oc433704660 h 10 HYPERLINK l _Toc433704661接觸分析實實例包含含初始間間隙 PAGEREF _Toc433704661 h 14 HYPERLINK l _Toc433704662耦合小程序序 PAGEREF _Toc433704662 h 19 HYPERLINK l _Toc433704663流固耦合練練習(xí) PAGEREF _Toc433704663 h 21 HYPERLINK l _Toc433704664一個流固耦耦合的例例子 PAGEREF _Toc433704664 h 22 HYPERLINK l _Toc433704665使用物理環(huán)環(huán)境

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計算結(jié)果 PAGEREF _Toc433704674 h 48ANSYSS流固耦耦合的方方式一般說來，ANSSYS的流固固耦合主主要有4種方式式：1，seqquenntiaal這需要用戶戶進行APDDL編程進進行流固固耦合sequeentiia指的是是順序耦耦合以采用MppCCII為例,你可以以利用ANSSYS和一個個第三方方CFD產(chǎn)品執(zhí)執(zhí)行流固固耦合分分析。在在這個方方法中，基基于網(wǎng)格格的平行行代碼耦耦合界面面(MppCCII) 將ANSSYS和CFD程序耦耦合起來來。即使使網(wǎng)格上上存在差差別，MpCCCI也能夠夠?qū)崿F(xiàn)流流固界面面的數(shù)據(jù)據(jù)轉(zhuǎn)換。A

7、NSSYS CD中包含含有MpCCCI庫和一一個相關(guān)關(guān)實例。關(guān)關(guān)于該方方法的詳詳細信息息，參見見ANSSYS Couupleed-FFielld AAnallysiis GGuidde中的Seqquenntiaal CCoupplinn2，F(xiàn)SII soolveer流固耦合的的設(shè)置過過程非常常簡單，推推薦你使使用這種種方式3，mullti-fieeld sollverr這是FSII soolveer的擴展展，你可可以使用用它實現(xiàn)現(xiàn)流體，結(jié)結(jié)構(gòu)，熱熱，電磁磁等的耦耦合4，直接采采用特殊殊的單元元進行直直接耦合合，耦合合計算直直接發(fā)生生在單元元剛度矩矩陣一個流固耦耦合模態(tài)態(tài)分析的的例子11這是一個

8、流流固耦合合模態(tài)分分析的典典型事例例，采用用ANSSYS/MECCHANNICAAL可以完完成。處處理過程程中需要要注意以以下幾個個方面的的問題：1、單元的的選擇；2、流體材材料模式式；3、流固耦耦合關(guān)系系的定義義；4、模態(tài)提提取方法法。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 !非對稱稱模態(tài)提提取方法法處理流流固耦合合問題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一個流固耦耦合模態(tài)態(tài)分析的的例子22一實例，水水箱采用用SHEELL663單元，水水箱中的的水采用用FLUUID330單元，以以下即為為整個流流固耦合合模態(tài)計計算的命命令流文文件：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 !非對稱稱模態(tài)提提取方法法處理流流固耦合合問題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我得出的結(jié)結(jié)果是：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為了比較我我還計算算了水箱箱不裝水水的情況況，（水水箱體模模型的建建立和約約束與上上面的相相同，只只是把流流體部分分去掉），計計算結(jié)果果是：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。水箱裝裝水時，怎怎么計算算結(jié)果中中有0頻率啊啊？其中中不等于于0的頻率率是水箱箱體的頻頻率么？2。水箱不不裝水時時的結(jié)果果和裝水水時的相相差怎么么這么大大啊？還還是我理理解錯了了么？企盼指教一個流固耦耦合建模模的例子子/prepp7 !進行預(yù)預(yù)處理模模塊et,1,30, !定義1號單元元為Fluuid330 流固耦耦合單元元et,2,29 !定義2號單元元為Fluuid229平面流流體單元元et,3,30 , ,1 !定

16、義3號單元元為Fluuid330流體介介質(zhì)單元元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 !設(shè)置水水中聲速速mp,deens,1,110000 !設(shè)置流流體密度度sectyype, 1, beeam, T, ! 選取T型梁secofffseet,oriig !設(shè)置梁

17、梁的方向向secdaata,0.004,00.055,0.0022,0.02,0,00,0,0,00,0 !所建立的的圓柱殼殼體的參參數(shù)：圓圓柱長為為 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, !移動工工作平面面與選取取的側(cè)面面相切!重復(fù)上面面操作,形成四四個環(huán)肋肋面wpofff,0,0,-0.44 !工

18、作平平面回到到原點位位置上k,31,0.22,0,0.11 !定義環(huán)環(huán)肋的方方向點lsel,s,200 !選擇要要劃分為為環(huán)肋的的線段latt,4,55,5,311,400,1 !定義線線段物理理屬性lesizze,220,6 !劃分?jǐn)?shù)數(shù)目secnuum,11lmeshh,200 !劃分線線段type，2 ! 選取第第二種單單元lsel，s，1，4 ！選取取線段lesizze，all，10 ！線段段劃分?jǐn)?shù)數(shù)目lesizze，all，6ameshh，1 ！將通通過工作作平面的的面 1 進行劃劃分esizee，0.1，0 ！選取取沿面放放樣的網(wǎng)網(wǎng)格大小小mat，11 ！定義義放樣的的形成單單元的物

19、物理屬性性vdragg，1，9 ！進行放樣形成流體介質(zhì) 一加筋板在在水中的的模態(tài)分分析再給大家一一個實例例！考慮結(jié)構(gòu)在在水中的的自振頻頻率：例例子是一一加筋板板在水中中的模態(tài)態(tài)分析。命令流如下下：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總是出現(xiàn)eerroor 說矩陣陣不對稱稱，不可可以用lannb計算。總結(jié)：流體體單元不不能用對對稱的解解法應(yīng)該采用非非對稱解解法。一圓環(huán)在水水中的模模態(tài)分析析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! 創(chuàng)建四四分之一一模型CYL4,0,00,0.2544,0,0.2260335,990CYL4,0,00,0.260035,0,00.3112422,900! 選擇屬

26、屬性,網(wǎng)格劃劃分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! 關(guān)于YY軸鏡像像

27、nsym,x,110000,alll ! offfseet nnodee nuumbeer bby 110000esym,10000,alll! 關(guān)于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! 標(biāo)識流流固交接接面nsel,s,lloc,x,00.2660355esel,s,ttypee,22sf,alll,ffsi,1nsel,alllesel,alllFINISSH/soluuantyppe,mmodaalmodoppt,ddampp,100mxpannd,110,yeessolveefinissh為了便于對對比，也也對圓環(huán)環(huán)在空氣氣中做了了模態(tài)分分析finissh/cleaar/PREPP7!定義單元元類型ET,1,PLAANE442 !

29、sttruccturral eleemennt!材料屬性性MP,EXX,1,2.0068ee11MP,DEENS,1,779299MP,NUUXY,1,00! 創(chuàng)建四四分之一一模型CYL4,0,00,0.2544,0,0.2260335,990! 選擇屬屬性,網(wǎng)格劃劃分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! 關(guān)于YY軸鏡像像nsym,x,1100

30、00,alll ! offfseet nnodee nuumbeer bby 110000esym,10000,alll! 關(guān)于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）考慮流流固耦合合，做模模態(tài)分析析時流體體單元是是否只能能用fluuid229（2d）和fluuid330（3d),對于fluuid1129和fluuid1130在耦合合中具體體起到什什么作用用，能不不能不設(shè)設(shè)，而用用邊界約約束條件件代替？2）流體范范圍怎樣樣確定，如如本例中中（CYLL4,00,0,0.2260335,00,

34、0.312242,90），外外半徑為為0.3312442。如果果不是環(huán)環(huán)形的，如如一塊當(dāng)當(dāng)水板，該該怎樣考考慮？3）如果不不考慮流流體的壓壓縮性，把把聲速設(shè)設(shè)的很大大，MP,SONNC,22,1ee20，就可可以了。4）從自振振頻率可可以看出出，在水水中和在在空氣中中，圓環(huán)環(huán)的自振振頻率差差別特別別大，且且振型也也大相徑徑庭，為為什么？在水中中時，模模態(tài)提取取方法用用dammp（為什什么不能能用unssym），特特征值的的虛部代代表角頻頻率，為為什么第第一階為為正，第第二階為為負，而而第三階階和第四四階都為為0，第六六階、八八階、十十階都為為負。應(yīng)應(yīng)該是從從小到大大才對？5）在空氣氣中時，模模

35、態(tài)提取取方法用用lannb，為什什么第一一階第二二階的頻頻率都為為0。請高手指點點迷津，急急盼中對以上問題題的解答答：頻率為零，一一般是發(fā)發(fā)生了剛剛體位移移，估計計你是把把水抽走走，而沒沒有限制制圓環(huán)。1。圓環(huán)在在水中振振動必然然導(dǎo)致波波動(其實就就是聲波波)在水中中傳播，當(dāng)當(dāng)聲波到到達水的的另一個個界面時時就會發(fā)發(fā)生反射射(除非水水和另一一個相鄰鄰體的聲聲阻抗是是匹配的的)。水和和金屬中中的聲速速相差不不大，即即可壓縮縮性相差差不大。兩兩種可壓壓縮性相相差不大大的物質(zhì)質(zhì)的相互互作用對對兩者影影響都很很大。圓圓環(huán)在水水中振動動，水對對圓環(huán)的的反作用用是由于于反射波波引起的的，流固固耦合中中采用

36、fluuid1129和130就是最最大程度度的減弱弱反射波波。2。聲波從從圓環(huán)開開始傳播播，隨著著傳播距距離的增增加，波波陣面不不斷增大大，振幅幅不斷減減小。同同時由于于水的衰衰減，聲聲波也不不斷減弱弱。如果果水的空空間越大大，則反反射波返返回圓環(huán)環(huán)的路徑徑越長，衰衰減也就就越多，影影響也就就越小。fluuid1129和130對反射射波的衰衰減(通過很很小的反反射實現(xiàn)現(xiàn))有限，因因此還需需要水要要有足夠夠的空間間。fluuid1129和130離結(jié)構(gòu)構(gòu)應(yīng)該大大于0.2(=c/f,cc為水中中聲速)。以上上的做法法在誤差差允許的的情況下下等效于于水在無無限大水水空間中中的情況況。如果果是擋水水板，

37、水水就是有有限空間間了，情情況也不不一樣。3。聲速加加大情況況也不一一樣，就就是不知知是不是是你所要要的情況況？4。空氣作作為介質(zhì)質(zhì)，由于于其聲速速比金屬屬小很多多，可壓壓縮性大大很多，影影響可以以忽略不不計。而而水的影影響就不不同了。這這可能就就是頻率率和振型型不同的的原因吧吧？我試試了你的的例子，各各種提取取方法都都可以。5?？諝獾牡挠绊懞龊雎圆挥嬘嫞虼舜诵枰獙A環(huán)進進行約束束。你沒沒有約束束，那么么就會發(fā)發(fā)生靜態(tài)態(tài)位移即即頻率為為零。圓圓環(huán)有兩兩個對稱稱軸，因因此會發(fā)發(fā)生頻率率成對出出現(xiàn)的情情況。也也就是說說，兩個個方向上上有同樣樣的振型型。接觸分析實實例包含初初始間隙隙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耦合小程序序最近用到耦耦合,寫了一一段小程程序,奉獻出出來,與大家家共享。如果有很多多節(jié)點，每每兩個節(jié)節(jié)點位置置相同，如如果將這這些雜亂亂無章的的節(jié)點耦合，是件件很麻煩煩的事，可可用這段段程序，輕輕松解決決。cpnumm=0cmsell,s,n-zzhonng !需要耦耦合的節(jié)節(jié)點*GET,n_nnum,NODDE,COUUNT, , , , !節(jié)點總總數(shù)*do,ii,1,n_nnumcmsell,s,n-zzhonng*GET,n_llowe

46、est,NODDE,NUMM,MIIN, , , , !號碼最最小的節(jié)節(jié)點*GET,n_xx,NOODE,n_lloweest,LOCC,X !該節(jié)點點坐標(biāo)*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 !尋找與與該節(jié)點點位置相相同的節(jié)節(jié)點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 !位置相相同的節(jié)節(jié)點形成成一個組組cmsell,

47、s,n-zzhonngcmsell,u,n_ccp_ccpcm,n-zhoong,nodde !取消這這些點后后剩余的的點形成成組*GET,n_nnum_1,NNODEE,CCOUNNT, , , , !節(jié)點總總數(shù)*if,nn_nuum_11,ltt,2,exiit !如果節(jié)節(jié)點數(shù)小小于二則則退出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因為你選用用的單元元有六個個自由度度,如果只只約束三三個,程序是是不會運運行的.另:三次循循環(huán)語句句的I相等,約束UY時,UX的耦合合就被刪刪掉了,最后只只剩UZ了這樣修改：流固耦合練練習(xí)/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一個流固耦耦合的例例子這個例子關(guān)關(guān)于裝有有水的水水杯旋轉(zhuǎn)轉(zhuǎn)，是軸軸對稱問問題，為為了簡化化，所以以選擇了了平面模模型。*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本人拋磚引引玉：使用物理環(huán)環(huán)境法進進行流固固耦合的的實例及及講解ansyss從9.0發(fā)展到10.0，一個個最大的的進步就就是流固固耦合計計算更加加規(guī)范，這這一點已已遠領(lǐng)先先于其他他同類軟軟件，實實現(xiàn)了單單向耦合合到即時時雙向耦耦合的飛飛躍，使使用戶對對于解決決流固耦耦合問題題又多了了一種選選擇，希希望大家家對多種種方法物理環(huán)環(huán)境轉(zhuǎn)換換，fsi，mfx等進行行討論，提提供一下下案例流道中有一一橡膠墊墊阻礙水水的流動動，入口口速度為為2m/s，其他他參數(shù)將將在命令令流中詳詳細給出出。求解解水通過過此流道道的壓力力降，以以及穩(wěn)態(tài)態(tài)條件下下

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