24小时热门版块排行榜

返回列表

微言1021

木虫 (小有名气)

应助: 2 (幼儿园)
金币: 1887.7
红花: 1
帖子: 113
在线: 60.2小时
虫号: 1919240
注册: 2012-08-01
性别: MM
专业: 力学中的基本问题和方法

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

本人正在学习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 ]

回复此楼