马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。
您需要 登录 才可以下载或查看,没有账号?我要加入
x
焊接热应力分析从热分析转入结构分析时温度载荷读不进去,部分程序如下:
!开始结构分析
/prep7
etchg,tts
et,3,combin14
r,1,1000000
mptemp
mptemp,1,0,300,600,1000,1300,1400,2000
mpdata,ex,1,1,2e11,1.86e11,1.35e11,2e10,2e7,2e7,2e7
mpdata,ex,2,1,2e11,1.86e11,1.35e11,2e10,2e7,2e7,2e7
mpdata,ex,3,1,2e11,1.86e11,1.35e11,2e10,2e7,2e7,2e7
mpdata,prxy,1,1,0.33,0.33,0.35,0.36,0.4,0.4
mpdata,prxy,2,1,0.33,0.33,0.35,0.36,0.4,0.4
mpdata,prxy,3,1,0.33,0.33,0.35,0.36,0.4,0.4
mpdata,alpx,1,1,12e-6,12.6e-6,13e-6,13.4e-6,0,0,0
mpdata,alpx,2,1,12e-6,12.6e-6,13e-6,13.4e-6,0,0,0
mpdata,alpx,3,1,12e-6,12.6e-6,13e-6,13.4e-6,0,0,0
!mp,reft,1,20
!mp,reft,2,20
!mp,reft,3,20
tb,bkin,1,5
tbtemp,15
tbdata,1,280e6,2e10
tbtemp,800
tbdata,1,10e6,1e7
tbtemp,2500
tbdata,1,10e6,1e6
tbcopy,bkin,1,2
k,50,1.1,0.3
k,51,1.1,0.3,length
l,50,9
l,51,21
type,3
real,1
esize,,1
lmesh,42,43,1
finish
/solu
dk,1,all
da,9,ux
da,13,ux
antype,trans
outpr,all,all
outres,all,all
!kbc,0
kbc,1
autots,on
time,0.02
esel,s,mat,,3
esel,a,mat,,2
ekill,all
esel,all
nsubst,1
solve
time=0.02
*do,kk,1,2
*if,kk,eq,1,then
v_num=1
time_inc=length/L_number/v
*endif
*if,kk,eq,2,then
v_num=2
time_inc=length/L_number/1.5/v
*endif
*do,j,1,L_number
vsel,s,,,v_num
eslv,s,1
nsel,r,loc,z,(j-1)*length/L_number+0.001,j*length/L_number-0.001
esln,r
ealive,all
alls
time=time+0.0001
time,time
solve
time=time+time_inc-0.0001
bfedele,all,temp
ldread,temp,,,time,,hanjie,rth
nsubst,10,50,3
lnsrch,on
time,time
!nlgeom,on
solve
*enddo
time=time+3600
bfedele,all,temp
ldread,temp,,,time,,hanjie,rth
time,time
solve
!save
*enddo
finish
|