24小时热门版块排行榜    

查看: 501  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 为什么资助数各大高校都创新高,自己申请怎么就这么难 +8 Kittylucky 2026-08-27 8/400 2026-08-27 18:33 by raiy123
[教师之家] 导师吐槽:我怎么摊上了这么个极品研究生! +7 苏东坡二世 2026-08-23 7/350 2026-08-27 18:14 by 瞬息宇宙
[基金申请] 青基如何破局 +6 lyz123lyz 2026-08-27 7/350 2026-08-27 17:31 by 宁静致远sy
[基金申请] 投票:  有多少人是今天查系统知道结果的? +15 爱看书的可乐 2026-08-26 17/850 2026-08-27 17:25 by xxniao123
[基金申请] 看板上这么多中的,有点像50人群里49个人都是骗子的那种感觉…… +5 a089 2026-08-26 6/300 2026-08-27 14:05 by jonewore
[基金申请] 面上合作单位盖章 +3 ssyjh 2026-08-27 3/150 2026-08-27 09:59 by 阿拉贡aragon
[基金申请] 范进中举一文的中心思想 +9 炎黄贵胄 2026-08-22 10/500 2026-08-26 15:40 by semaglutide
[基金申请] 2026国自然函评费到账 +22 羊腰板 2026-08-21 25/1250 2026-08-26 15:00 by zuocuiping
[基金申请] 能否退出参与的面上项目解除限项 +23 koalala 2026-08-24 26/1300 2026-08-26 14:29 by 宝贝虫子
[基金申请] 出来了 +9 trojank 2026-08-26 9/450 2026-08-26 14:25 by 宝贝虫子
[基金申请] 为什么国自然不能直接公布 +4 bjdxyxy 2026-08-26 4/200 2026-08-26 13:12 by qingmu1201
[基金申请] 国合里面能看到了 +7 一怀馨秋 2026-08-26 7/350 2026-08-26 11:23 by zhaosm1982
[基金申请] 明天应该可查了!? +6 chengyan1220 2026-08-23 6/300 2026-08-25 19:45 by zfd97
[基金申请] 某些机构,以效率低为荣,以效率低作为存在感 +9 yuleib84 2026-08-25 10/500 2026-08-25 17:14 by alexon
[基金申请] 建议基金发布提前给出明确的时间点 +13 kulium 2026-08-21 16/800 2026-08-24 16:27 by superceng
[教师之家] 跳槽后在研项目怎么办? +5 简单化xn 2026-08-22 10/500 2026-08-23 12:38 by 简单化xn
[基金申请] 时间戳今天,20号变了 +5 archvillain 2026-08-20 5/250 2026-08-22 06:12 by hui_daxiao
[基金申请] 看来今天不会放榜了? +8 chengyan1220 2026-08-21 11/550 2026-08-21 17:52 by dcqxinyang
[基金申请] 应该是下周三26日公布了吧? +4 哈哈蛤? 2026-08-21 4/200 2026-08-21 10:58 by Vivilian
[基金申请] 基金啊基金 +4 longfie172 2026-08-20 4/200 2026-08-21 08:58 by mark mao
信息提示
请填处理意见