用ANSYS分析过整个桥梁施工过程
Q: I must build the whole model and kill the elements that don't take part in the
analysis of certain erection stage, so i think the multiple steps is not
a useful way to solve the problem
A: You may try to utilize multiple steps and set one load case a step. Meanwhile,large deformation effects option is recommended.
用单元死活模拟浇铸过程中的温度分布
要想对浇铸过程的温度场进行分析,必须熟悉下面两个方面的内容:1.瞬态温度场的分析,2.单元死活的应用。
瞬态温度场分析:在进行瞬态温度分析时,我们常遇到的一个问题是温度结果明显不合理,:计算得到的温度高于给定的最高温度或低于给定的最低温度。造成这种结果的原因有两个:1、单元不合理,网格太大。2、时间积分方法不恰当。针对上述两个原因,我们采取的措施是:1、细划网格。2、选择向后积分法(THETA=1)。
单元死活:应用单元死活进行分析时,需要用到多个载荷步。在不同的载荷步之间单元的死活状态不同,ANSYS的通用动画生成工具生成的动画无法表现单元状态随时间的变化,因此必须使用ANSYS的DISPLAY Utility来生成能够表现单元死活状态的动画。
下面是一个简单的模拟浇铸过程的例子,环境温度为70,浇铸温度为600,计算结果应是:最低温度不低于70,最高温度不高于600。如果我们划分的网格不适当,则会出现最低温度低于70或最高温度高于600的情况,因此必须细划网格并选择THETA=1。计算完成后使用ANSYS的DISPLAY Utility来生成动画文件。
fini
/cle
/prep7
et,1,55
et,2,55
KXX , 1,1.99/3600
DENS,1,0.00073
C ,1,0.100
KXX,2,7.38E-6
DENS,2,9.35E-5
c,2,1
k,1,-.5,-.5
k,4,6,-.5
k,12
k,13,5.5
k,21,-.5,6
k,22,,6
k,23,5.5,6
k,24,6,6
a,12,13,23,22
a,1,12,22,21
a,1,4,13,12
a,4,24,23,13
lsel,s,,,5,6
lsel,a,,,9,11,2
lesize,all,,,8
lsel,all
lesize,all,,,20
asel,s,,,1
aatt,2,2,2
asel,all
amesh,all
fini
/solu
antype,tran
!solc,on
timint,off
lsel,s,,,7,8
lsel,a,,,10
nsll,s
sf,all,conv,100,70
alls
esel,s,type,,2
nropt,full
ekill,all
alls
tunif,100
time,1
outres,all,all
solve
dtop=0
tim=1
delt=.05
timint,on
tintp,,,,1
*do,ii,1,20
ddele,all,temp
esel,s,type,,2
nsle,s
tim=tim+delt
time,tim
dtop=dtop+.3
nsel,r,loc,y,0,dtop
esln,s,1
ealive,all
nsel,r,loc,y,dtop
esln,s
nsle,s,1
d,all,temp,600
esel,s,live
nall
eall
solve
*enddo
fini
/show,tmpplt,f33
/post1
dtop=0
*do,ii,2,21
set,ii
dtop=dtop+.3
nsel,r,loc,y,0,dtop
esln,s,1
esel,a,type,,1
nall
plns,temp
*enddo
在ANSYS5.6中如何施加函数变化的表面载荷
ANSYS5.7版本具有函数加载功能,可以很方便地在模型表面施加函数变化的各种载荷,在ANSYS5.6版本中,也可以通过变通的方式来实现此功能,其思路是:
首先选定所要施加函数变化表面载荷的表面上的节点,利用ANSYS的参数数组和嵌入函数知识写一简单的命令流,定义好相应节点位置的面载荷值,然后通过在节点上施加面载荷来完成。
下面以在一圆柱表面施加函数变化载荷为例:
/prep7
et,1,45
cyl4,,,0.5,,,,3
vsweep,all
asel,s,loc,y,0.01,1
nsla
!
*get,nmax,node,,num,max,
*get,nmin,node,,num,min,
*afun,deg
*dim,t1,array,nmax,1,1,
csys,1
*do,k,nmin,nmax
*if,nsel,eq,1,then
t1=1000*sin(ny)
*else
t1=0
*endif
*enddo
!
sffun,pres,t1(1)
sf,all,pres,0 |