24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

lhj769349654

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 6.5
帖子: 27
在线: 7.1小时
虫号: 3576901
注册: 2014-12-04
专业: 水工结构和材料及施工

[求助] ANSYS 焊接残余应力模拟已有1人参与

热分析时热源荷载是取得生热率和生死单元的方法，这是其中一个荷载步的温度场。主要命令流如下：
/CLEAR
/FILNAME, exercise10
/title,weldingstress10
/prep7
!**************定义参数及热源********************
L=0.2
LSIZE=0.002
MAX_X=L/LSIZE
A= 0.2827E-04
V=0.004
U=25
I=180
K=0.75
HGENUM=(K*U*I)/(A*V)
!**************定义单元及材料参数********************
et,1,mesh200
KEYOPT,1,1,6
ET,2,solid70
dof,temp
mp,dens,1,7850
mptemp,1,20,100,200,400,600,800
mptemp,7,1000,1200,1500,1700
mpdata,kxx,1,1,48.2,50.7,46.5,38.1,30,26
mpdata,kxx,1,7,26.8,29.2,35,14
mpdata,c,1,1,470,490,515,590,725,835
mpdata,c,1,7,650,654,660,780
mptemp,1,20,200,400,700,750,880
mptemp,7,1000,1200,1500,20000
mpdata,ex,1,1,20.6e10,19e10,17e10,9e10,2e10,1.5e10
mpdata,ex,1,7,1.3e10,1.0e10,0.5e10,0.1e10
mpdata,prxy,1,1,0.28,0.29,0.31,0.35,0.40,0.45
mpdata,prxy,1,7,0.46,0.47,0.48,0.48
tb,biso,1,6,2
tbtemp,20,1
tbdata,,345e6,
tbtemp,200,2,
tbdata,,250e6,
tbtemp,400,3,
tbdata,,100e6,
tbtemp,700,4,
tbdata,,50e6,
tbtemp,880,5,
tbdata,,30e6,
tbtemp,2500,6,
tbdata,,20e6,
mp,alpx,1,1.06e-5
mpcopy,,1,2
tb,biso,2,6,2
tbtemp,20,1
tbdata,,345e6,
tbtemp,200,2,
tbdata,,250e6,
tbtemp,400,3,
tbdata,,100e6,
tbtemp,700,4,
tbdata,,50e6,
tbtemp,880,5,
tbdata,,30e6,
tbtemp,2500,6,
tbdata,,20e6,
!**************建立模型********************
rectng,-0.002,0.002,-0.016,0
rectng,-0.0104853,-0.002,-0.016,0
rectng,0.002,0.0104853,-0.016,0
rectng,-0.030,-0.0104853,-0.016,0
rectng,0.0104853,0.030,-0.016,0
rectng,-0.120,-0.030,-0.016,0
rectng,0.030,0.120,-0.016,0
cyl4,-0.002,0,0.0084853,135,,180
cyl4,0.002,0,0.0084853,0,,45
aglue,all
k,34,-0.008,0.016
k,35,0.008,0.016
a,29,4,3,33,35,34,29
rectng,-0.008,0.008,0.016,0.036
rectng,-0.008,0.008,0.036,0.180
aglue,all
numcmp,all
  type,1
esize,0.002
amesh,1,2
amesh,5,8
  esize,0.005
amesh,4
amesh,11,12
  esize,0.002
amesh,3
amesh,9,10
EXTOPT,aclear,1,
EXTOPT,attr,0,
ESIZE,,100
type,2
mat,1
Vext,5,6,,,,0.2
EXTOPT,aclear,1,
EXTOPT,attr,0,
ESIZE,,100
type,2
mat,2
Vext,1,4,1,,,0.2
Vext,7,12,1,,,0.2
Nummrg,elem
Nummrg,node
Nummrg,kp
!**************非线性分析求解选项设置********************
/solu
antype,trans
trnopt,full
nropt,full,,on
AUTOTS,on
lnsrch,on
!**************杀死焊缝单元********************
allsel
/replot
esel,s,mat,,1
ekill,all
allsel
esel,s,live
eplot
!**************施加对流荷载********************
allsel
asel,s,,,1,4
asel,a,,,7,12
asel,a,,,21,22
asel,a,,,26
asel,a,,,29
asel,a,,,31,33
asel,a,,,35,41
asel,a,,,45,46
asel,a,,,50,52
asel,a,,,54,55
asel,a,,,57,65
esla,s
nsle,s
sf,all,conv,30,20
!**************施加初始温度条件********************
nsel,all
tunif,20
!**************开始热分析********************
timint,on,therm
*do,i,1,40
esel,s,mat,,1
nsle,s
nsel,r,loc,z,0.002*(i-1),0.002*i
esln,s,all
ealive,all
eplot
bfe,all,HGEN,,HGENUM
time,0.5*i
kbc,1
DELTIM,0.25,0.25,0.5
allsel,all
   OUTRES,ALL,last
   tsres,erase
SOLVE
bfeDELE,all,HGEN,,HGENUM
allsel,all
*ENDDO
!**************冷却过程********************
time,60
AUTOTS,1
NSUBST,20,20,20
KBC,0
TSRES,ERASE
solve
time,100
AUTOTS,1
NSUBST,8,8,8
KBC,0
TSRES,ERASE
solve
time,1100
AUTOTS,1
NSUBST,20,20,20
KBC,0
TSRES,ERASE
solve
time,2100
AUTOTS,1
NSUBST,10,10,10
KBC,0
TSRES,ERASE
solve
!**************进入应力分析********************
/prep7
et,3,solid45
dof,ux,uy,uz
allsel
eplot
etchg,tts
lsclear,solid
!**************施加位移边界条件********************
asel,s,,,38
nsla,s,1
dsym,symm,y
asel,s,,,64
nsla,s,1
d,all,ux,0
asel,s,,,61
nsla,s,1
d,all,ux,0,,,,uy
allsel
eplot
!**************非线性分析求解选项设置********************
/solu
antype,trans
trnopt,full
nropt,full,,on
nlgeom,on
AUTOTS,on
lnsrch,on
neqit,50
!**************杀死焊缝单元********************
allsel
/replot
esel,s,mat,,1
ekill,all
allsel
esel,s,live
eplot
!**************应力分析********************
timint,on
*do,i,1,100
esel,s,mat,,1
nsle,s
nsel,r,loc,z,0.002*(i-1),0.002*i
esln,s,all
ealive,all
eplot
time,0.5*i
kbc,1
DELTIM,0.25,0.25,0.5
ldread,temp,i,,,,,rth
allsel,all
   OUTRES,ALL,last
   tsres,erase
SOLVE
*ENDDO
allsel,all
eplot
time,60
AUTOTS,1
NSUBST,20,20,20
KBC,0
TSRES,ERASE
solve
time,100
AUTOTS,1
NSUBST,8,8,8
KBC,0
TSRES,ERASE
solve
time,1100
AUTOTS,1
NSUBST,20,20,20
KBC,0
TSRES,ERASE
solve
time,2100
AUTOTS,1
NSUBST,10,10,10
KBC,0
TSRES,ERASE
solve
可产生的温度为什么没有扩散，是模型网格还是单元的问题，希望能帮忙修改下，然后就是转应力场分析时有问题，收敛也存在问题。

360截图20160623094633177.jpg

回复此楼