| 1 | 1/1 | 返回列表 |
| 查看: 416 | 回復(fù): 0 | ||
微言1021木蟲 (小有名氣)
|
[求助]
尿素合成塔的流固耦合
|
|
本人正在學(xué)習(xí)ansys中的流固耦合,尿素合成塔的塔壁是用單元shell63模擬,塔內(nèi)的流體(水)用單元fluid80模擬。想用瞬態(tài)動力學(xué)分析的full方法求在地震作用下塔內(nèi)液體的晃動情況以及塔體隨液體的高度變化的變形情況。我采用了節(jié)點耦合命令cpintf,也參考了其他人的方法,但是求出的結(jié)果是:流體單元底部變形很大,甚至超出了塔體,而塔的變形很小。請各位高手幫忙解答一下,謝謝了。 下面是我的命令流,我直接黏貼上去,請大家?guī)臀倚薷男薷模?br /> /filname,ten-fuild-solid coupling6 /title,unconsiderable6 *set,d,1.8 *set,m_t,.07 *set,b_t,.16 *set,u_t,.6 *set,q_t,.06 *set,t_h1,.9 *set,t_h2,9.1 *set,t_h3,1.5 /prep7 et,1,shell63 et,2,fluid80 mp,ex,1,2.05e11 mp,ex,2,2e9 mp,prxy,1,.3 mp,dens,1,7.83e3 mp,dens,2,1000 mp,visc,1,1.0e-6 keyopt,2,2,1 r,1,m_t,m_t,m_t,m_t r,2,b_t,b_t,m_t,m_t r,3,m_t,m_t,u_t,u_t r,4,q_t,q_t,m_t,m_t /pnum,kp,1 /pnum,line,1 /pnum,area,1 k,1,,,-t_h2 k,2,t_h1,,-t_h2+t_h1 k,3,t_h1,,t_h1 k,4,,,t_h1 k,5,t_h1,,-t_h2-t_h3 k,6,t_h1,, k,7, larc,1,2,4,d/2 l,2,6 l,6,3 l,3,4 l,4,7 l,7,1 l,2,5 l,6,7 al,1,2,8,6 al,3,4,5,8 vrotat,1,2,,,,,1,4,, arotat,7,,,,,,1,4,, vglue,all nummrg,all numcmp,all asel,u,loc,z,0 asel,u,loc,x,0 asel,u,loc,y,0 cm,wall,area vdele,all csys,1 asel,s,loc,z,-t_h2+t_h1 asel,s,loc,z,t_h1 asel,a,loc,x,d/2 cm,wall,area lsel,s,loc,z,-t_h2,-t_h2-t_h3+0.1 lesize,all,,,20 lsel,s,loc,z,-t_h2+t_h1 lsel,a,loc,z,t_h1 lsel,a,loc,z,-t_h2-t_h3 lsel,a,loc,z,0 lesize,all,,,20 lsel,s,length,,t_h2 lsel,a,length,,-t_h1+t_h2 lesize,all,,,20 lsel,s,loc,z,0.1,t_h1-0.1 lesize,all,,,5 lsel,s,loc,z,-t_h2+0.1,-t_h2+t_h1-0.1 lesize,all,,,20 aplot,all asel,s,loc,z,-t_h2,-t_h2-t_h3 real,4 mshkey,0 mshape,0,3d amesh,all asel,s,loc,z,t_h1,-t_h2+t_h1 asel,r,loc,x,d/2 real,1 mshkey,0 mshape,0,3d amesh,all asel,s,loc,z,-t_h2,-t_h2+t_h1 real,2 mshkey,0 mshape,0,3d amesh,all asel,s,loc,z,t_h1 real,3 mshkey,0 mshape,0,3d amesh,all !建液體 lplot,all lsel,all lgen,2,1,2,1 l,7,22 l,7,20 al,43,44,45,46 vrotat,33,,,,,,1,4,, vglue,all cm,fluid,volu lsel,s,,,45,49,4 lsel,a,,,54,59,5 lsel,a,,,50,51,1 lsel,a,,,55,56,1 lsel,a,,,60,63,1 lesize,all,,,20 lsel,s,,,46 lesize,46,,,20 vplot type,2 mat,2 mshape,0,3d mshkey,1 vmesh,all 節(jié)點耦合 nplot,all nrotate,all cpintf,ux,0.0001 nsel,s,loc,z,-t_h2,-t_h2+t_h1 nrotate,all cpintf,uz,0.0001 施加約束 nsel,s,loc,z,-t_h2-t_h3 d,all,all !求解 *dim,tjx,array,2,190,0,,, *dim,tjy,array,2,190,0,,, *create,ansuitmp *vread,tjx(1,1),'ACELX','TXT','',190,,,,,, (e9.3,e11.3) *end /input,ansuitmp *creat,ansuitmp *vread,tjy(1,1),'ACELY','txt','',190,,,,,, (e9.3,e11.3) *end /input,ansuitmp fini /solu antype,4 trnopt,full timint,off outres,basic,all kbc,1 time,1e-10 nsub,4,8 acel,,,9.8 alls solve timint,on *do,t,1,6,1 time,0.1*t kbc,0 nsub,1 alphad,0.05 betad,0.01 acel,tjx(2,t),tjy(2,t) alls solve *enddo save [ Last edited by 微言1021 on 2012-8-1 at 11:17 ] |

找到一些相關(guān)的精華帖子,希望有用哦~
| 1 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 317求調(diào)劑 +3 | 十閑wx 2026-03-24 | 3/150 |
|
|---|---|---|---|---|
|
[考研] 266求調(diào)劑 +5 | 陽陽哇塞 2026-03-27 | 5/250 |
|
|
[考研] 求調(diào)劑 +3 | 劉柯@ 2026-03-24 | 4/200 |
|
|
[考研] 085600,材料與化工321分調(diào)劑 +3 | 大饞小子 2026-03-27 | 4/200 |
|
|
[考研] 0856材料專碩353求調(diào)劑 +9 | NIFFFfff 2026-03-20 | 9/450 |
|
|
[考研] 329求調(diào)劑 +7 | 鈕恩雪 2026-03-25 | 7/350 |
|
|
[考研] 343求調(diào)劑 +4 | 贈我一本書 2026-03-23 | 4/200 |
|
|
[考研] 333求調(diào)劑 +6 | wfh030413@ 2026-03-23 | 6/300 |
|
|
[考研] 327求調(diào)劑 +7 | prayer13 2026-03-23 | 7/350 |
|
|
[考研] 生物學(xué)學(xué)碩,一志愿湖南大學(xué),初試成績338 +4 | YYYYYNNNNN 2026-03-26 | 4/200 |
|
|
[考研] 材料277求調(diào)劑 +5 | min3 2026-03-24 | 5/250 |
|
|
[考研] 315分求調(diào)劑 +5 | 26考研上岸版26 2026-03-26 | 5/250 |
|
|
[考研] 一志愿上海交大生物與醫(yī)藥專碩324分,求調(diào)劑 +6 | jiajunX 2026-03-22 | 6/300 |
|
|
[考研] 一志愿南航 335分 | 0856材料化工 | GPA 4.07 | 有科研經(jīng)歷 +6 | cccchenso 2026-03-23 | 6/300 |
|
|
[考研] 材料專碩 335 分求調(diào)劑 +4 | 拒絕冷暴力 2026-03-25 | 4/200 |
|
|
[考研] 一志愿吉林大學(xué)材料與化工303分求調(diào)劑 +4 | 為學(xué)666 2026-03-24 | 4/200 |
|
|
[考研]
|
黃粱一夢千年 2026-03-24 | 3/150 |
|
|
[考研] 一志愿河北工業(yè)大學(xué)0817化工278分求調(diào)劑 +7 | jhybd 2026-03-23 | 12/600 |
|
|
[考研] 材料 271求調(diào)劑 +5 | 展信悅_ 2026-03-21 | 5/250 |
|
|
[考研] 一志愿重慶大學(xué)085700資源與環(huán)境總分308求調(diào)劑 +7 | 墨墨漠 2026-03-20 | 7/350 |
|