有奖求助：生死单元两道焊温度场模拟

talentiankun

/units,si

/config,nres,10000
!定义参数值
pi=3.1415926
y_dist=0.04
!钢板宽度
x_dist=0.05
!钢板长度
z_dist=0.008
!钢板厚度
x_size=0.002
!x向单元网格长度
y_size=x_size
!y向单元网格长度
weld_v=0.01
!焊接速度
weld_tm=x_dist/weld_v
!焊接时间

tinc=x_size/weld_v
!单元子步时间
weld_I=180
!焊接电流
weld_U=20
!焊接电压
ef=0.75
!焊接热效率
weld_r=0.0071
!有效加热半径
Qmax =3*weld_I*weld_U*ef/(pi*weld_r**2)
!最大热流高斯密度
vdist1=0.008

vdist2=0.014
max_time=1+weld_tm/tinc
max_x=1+x_dist/x_size

max_y=1+2*vdist1/x_size



/prep7
et,1,solid70
keyopt,1,2,1
et,2,surf152
r,2
keyopt,2,4,1
keyopt,2,5,0
keyopt,2,8,1
mptemp,1,20,250,500,750,1000,1500,1700,2500

mpdata,kxx,1,1,50,47,40,27,30,35,140,142
mpdata,c,1,1,460,480,530,675,670,660,780,820
mpdata,hf,1,1,100,350,520,1000,1500,3000,3100,3500
mpdata,dens,1,1,7820,7700,7610,7550,7490,7350,7300,7090

!建模,划分网格
/PREP7

ET,1,SOLID70
ET,2,SURF152
BLOCK,0,0.05,-0.02,0.02,0,0.008,
ESIZE,0.002,0,

MSHAPE,0,3D
MSHKEY,1
CM,_Y,VOLU

VSEL, , , ,
1
CM,_Y1,VOLU
CHKMSH,'VOLU'

CMSEL,S,_Y

VMESH,_Y1

CMDELE,_Y

CMDELE,_Y1

CMDELE,_Y2

FINISH

Allsel
Nsel,s,loc,z,0,0.004
Nsel,r,loc,y,-0.002，0.002
Nsel,r,loc,x,0,0.05
Esln,,1
Cm,bdelem1, elem
Allsel
Nsel,s,loc,z,0.004,0.008
Nsel,r,loc,y,-0.002，0.002
Nsel,r,loc,x,0,0.05
Esln,,1
Cm,bdelem2, elem

finish


边界条件，求解
/solu

Ekill,bdelem2
*dim,weld_f,table,max_x,max_y,max_time,x,y,time

*do,k,1,max_time,1

*do,i,1,max_x,1

weld_f(i,0,k)=(i-1)*x_size

*enddo

*do,j,1,max_y,1


weld_f(0,j,k)=(j-1)*x_size-vdist1

*enddo

*do,i,1,max_x,1

*do,j,1,max_y,1

xcenter=weld_v*(k-1)*tinc

ycenter=0


  distance=sqrt(((i-1)*x_size-xcenter)**2+((j-1)*x_size-vdist1-ycenter)**2)

*if,distance,le,weld_r,then
weld_f(i,j,k)=Qmax*exp(-3*distance**2/(weld_r**2))!高斯移动热源公式

*else

weld_f(i,j,k)=0

*endif


*enddo

*enddo

weld_f(0,0,k)=(k-1)*tinc

*enddo

antype,trans
tunif,20
!初始温度20度
kbc,1
timint,on
autots,off
NROPT,FULL !设定牛顿－拉夫森选项

esel,s,type,,2
!选择表面效应单元
nsle,s
sf,all,conv,-1,20
!施加对流边界条件
allsel

tinc=tinc
!载荷子步时间
*do,t,1e-6,weld_tm,tinc
time,t

autots,0
deltim,,,,1
kbc,0
!渐变载荷
outres,all,all
tsres,erase
nsel, s,loc,z,0               有问题

sf, all,hflux,%weld_f%
solve
*enddo




sfdele,all,hflux
!删除热流密度
allsel

t=t+4
!冷却4秒
time,t
deltim,0.2,0.1,0.5
solve
t=t+10
!冷却14秒
time,t
deltim,0.5,0.5,1
solve
t=t+10
!冷却24秒
time,t
deltim,1,1,2
solve
t=t+50
!冷却74秒
time,t
deltim,5,5,10
t=t+40
!冷却114秒
time,t
deltim,10,10,20

T=T+400
TIME,T
DELTIM,50,50,100
solve
finish
第二次无法加载
/solu
Ealive,bdelem2
*dim,weld_f1,table,max_x,max_y,max_time,x,y,time

*do,m,1,max_time,1

*do,n,1,max_x,1

weld_f1(n,0,m)=(n-1)*x_size

*enddo

*do,g,1,max_y,1


weld_f1(0,g,m)=(g-1)*x_size-vdist1

*enddo

*do,n,1,max_x,1

*do,g,1,max_y,1

xcenter=weld_v*(m-1)*tinc

ycenter=0


  distance=sqrt(((n-1)*x_size-xcenter)**2+((g-1)*x_size-vdist1-ycenter)**2)

*if,distance,le,weld_r,then
weld_f1(n,g,m)=Qmax*exp(-3*distance**2/(weld_r**2))!高斯移动热源公式

*else

weld_ f1(n,g,m)=0

*endif


*enddo

*enddo

weld_f1(0,0,m)=(m-1)*tinc

*enddo

antype,trans
kbc,1
timint,on
autots,off


esel,s,type,,2
!选择表面效应单元
nsle,s
sf,all,conv,-1,20
!施加对流边界条件
allsel
Tsres,erase

tinc=tinc
!载荷子步时间
*do,t, 520+1e-6,weld_tm+520,tinc
time,t

autots,0
deltim,,,,1
kbc,0
!渐变载荷
outres,all,all
tsres,erase
nsel,s,loc,z,z_dist/2

sf,all,hflux,%weld_f1%
solve
*enddo


sfdele,all,hflux
!删除热流密度
allsel

t=t+4
!冷却4秒
time,t
deltim,0.2,0.1,0.5
solve
t=t+10
!冷却14秒
time,t
deltim,0.5,0.5,1
solve
t=t+10
!冷却24秒
time,t
deltim,1,1,2
solve
t=t+50
!冷却74秒
time,t
deltim,5,5,10
t=t+40
!冷却114秒
time,t
deltim,10,10,20

T=T+400
TIME,T
DELTIM,50,50,100
solve

finish


又知道那里错的高手，请把改正意见发到我的邮箱（tiankun789789141@sina.com），咨询费RMB100，决不食言。

xbzhang

都什么错误提示啊？

zxw19851015

错误比较多

ycs2011

学习一下，看来你研究很长时间了