24小时热门版块排行榜    

查看: 503  |  回复: 0

微言1021

木虫 (小有名气)

[求助] 尿素合成塔的流固耦合

本人正在学习ansys中的流固耦合,尿素合成塔的塔壁是用单元shell63模拟,塔内的流体(水)用单元fluid80模拟。想用瞬态动力学分析的full方法求在地震作用下塔内液体的晃动情况以及塔体随液体的高度变化的变形情况。我采用了节点耦合命令cpintf,也参考了其他人的方法,但是求出的结果是:流体单元底部变形很大,甚至超出了塔体,而塔的变形很小。请各位高手帮忙解答一下,谢谢了。
下面是我的命令流,我直接黏贴上去,请大家帮我修改修改:
/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
节点耦合
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 ]
回复此楼
言之命至,心随己愿
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

智能机器人

Robot (super robot)

我们都爱小木虫

相关版块跳转 我要订阅楼主 学员8JxXkC 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 面上意见出来了 +4 黄鸟于飞Chao 2026-08-29 8/400 2026-08-29 14:50 by 黄鸟于飞Chao
[博后之家] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +4 ASdOkHsho7FD 2026-08-28 5/250 2026-08-29 14:24 by jCd0dEvKHShX
[公派出国] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 ASdOkHsho7FD 2026-08-28 5/250 2026-08-29 14:23 by jCd0dEvKHShX
[论文投稿] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 ASdOkHsho7FD 2026-08-28 4/200 2026-08-29 14:13 by jCd0dEvKHShX
[基金申请] 麻烦专家们看看评委们的意见(F口面上) +5 gdd2018 2026-08-28 10/500 2026-08-29 14:10 by Jacob678
[公派出国] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 ASdOkHsho7FD 2026-08-28 4/200 2026-08-29 14:02 by jCd0dEvKHShX
[硕博家园] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 ASdOkHsho7FD 2026-08-28 9/450 2026-08-29 13:41 by jCd0dEvKHShX
[找工作] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 ASdOkHsho7FD 2026-08-28 6/300 2026-08-29 11:51 by jCd0dEvKHShX
[基金申请] 投票:  有多少人是今天查系统知道结果的? +16 爱看书的可乐 2026-08-26 18/900 2026-08-29 10:18 by winsaint
[教师之家] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +4 ASdOkHsho7FD 2026-08-28 7/350 2026-08-29 09:26 by G6APbkg8SA6w
[基金申请] 怎么查啊 +6 huang1991js 2026-08-26 6/300 2026-08-28 08:42 by winsaint
[基金申请] 申请删除本帖 +6 lyz123lyz 2026-08-27 7/350 2026-08-27 17:31 by 宁静致远sy
[文学芳草园] 梦想 +3 myrtle 2026-08-26 3/150 2026-08-27 10:01 by angelyueyi
[基金申请] 为什么 国际(地区)合作与交流项目 没有放榜? 10+3 majunge000 2026-08-26 11/550 2026-08-27 08:42 by 北京莱茵编辑
[基金申请] 2026年8月25日国自然放榜前突然收到列入评审专家邮件,有关系吗? +25 木水思豆 2026-08-25 28/1400 2026-08-26 14:53 by draco1987
[基金申请] 能否退出参与的面上项目解除限项 +23 koalala 2026-08-24 26/1300 2026-08-26 14:29 by 宝贝虫子
[基金申请] 为什么国自然不能直接公布 +4 bjdxyxy 2026-08-26 4/200 2026-08-26 13:12 by qingmu1201
[基金申请] 今天务委会开完了,明天出结果吗 +19 angus9576 2026-08-25 23/1150 2026-08-26 10:03 by zp519
[基金申请] 牛来!米来!面来! +8 beefly 2026-08-26 8/400 2026-08-26 08:37 by xuzhipiao
[基金申请] 没有任何消息-是不是就凉了 +9 图啦图啦 2026-08-24 10/500 2026-08-25 11:59 by 南海小哥
信息提示
请填处理意见