多层焊接模拟分析热源加载求助

tto

多层焊接模拟分析热源加载求助


正在进行多层焊接模拟，采用单元生死，单元按生热率来加载，每个单元焊接作用时间为5s，计算后，时

间历程（温度和时间关系）图片如下：（和书上的曲线有很大差别）希望大家能一起讨论。命令流如下。

/config,nres,3000
/units,si
/prep7
et,1,PLANE55

MPTEMP,,,,,,,,  
MPTEMP,1,20  
MPTEMP,2,100
MPTEMP,3,200
MPTEMP,4,300
MPTEMP,5,400
MPTEMP,6,500
MPTEMP,7,750
MPTEMP,8,1000
MPTEMP,9,1500
MPTEMP,10,1700   
MPTEMP,11,2500  

MPDATA,KXX,1,,51.4
MPDATA,KXX,1,,47.2   
MPDATA,KXX,1,,36.8   
MPDATA,KXX,1,,30.25   
MPDATA,KXX,1,,25.31   
MPDATA,KXX,1,,20.2   
MPDATA,KXX,1,,33   
MPDATA,KXX,1,,35.2   
MPDATA,KXX,1,,36.5   
MPDATA,KXX,1,,38.9   
MPDATA,KXX,1,,40.3   


MPTEMP,,,,,,,,  
MPTEMP,1,20  
MPTEMP,2,100
MPTEMP,3,200
MPTEMP,4,300
MPTEMP,5,400
MPTEMP,6,500
MPTEMP,7,750
MPTEMP,8,1000
MPTEMP,9,1500
MPTEMP,10,1700   
MPTEMP,11,2500  

MPDATA,C,1,,480   
MPDATA,C,1,,485   
MPDATA,C,1,,498   
MPDATA,C,1,,506   
MPDATA,C,1,,526   
MPDATA,C,1,,539   
MPDATA,C,1,,675   
MPDATA,C,1,,670   
MPDATA,C,1,,660   
MPDATA,C,1,,780   
MPDATA,C,1,,820   

MPTEMP,,,,,,,,  
MPTEMP,1,20  
MPTEMP,2,100
MPTEMP,3,200
MPTEMP,4,300
MPTEMP,5,400
MPTEMP,6,500
MPTEMP,7,750
MPTEMP,8,1000
MPTEMP,9,1500
MPTEMP,10,1700   
MPTEMP,11,2500  

MPDATA,DENS,1,,7820
MPDATA,DENS,1,,7800
MPDATA,DENS,1,,7800
MPDATA,DENS,1,,7800
MPDATA,DENS,1,,7800
MPDATA,DENS,1,,7800
MPDATA,DENS,1,,7800
MPDATA,DENS,1,,7800
MPDATA,DENS,1,,7700
MPDATA,DENS,1,,7600
MPDATA,DENS,1,,7550

MPTEMP,,,,,,,,  
MPTEMP,1,20  
MPTEMP,2,100
MPTEMP,3,200
MPTEMP,4,300
MPTEMP,5,400
MPTEMP,6,500
MPTEMP,7,750
MPTEMP,8,1000
MPTEMP,9,1500
MPTEMP,10,1700   
MPTEMP,11,2500  

MPDATA,HF,1,,33.5
MPDATA,HF,1,,250
MPDATA,HF,1,,350
MPDATA,HF,1,,500
MPDATA,HF,1,,600
MPDATA,HF,1,,800
MPDATA,HF,1,,820
MPDATA,HF,1,,850
MPDATA,HF,1,,878
MPDATA,HF,1,,890
MPDATA,HF,1,,1000

!
MPTEMP,,,,,,,,  
MPTEMP,1,20  
MPTEMP,2,100
MPTEMP,3,200
MPTEMP,4,300
MPTEMP,5,400
MPTEMP,6,500
MPTEMP,7,750
MPTEMP,8,1000
MPTEMP,9,1500
MPTEMP,10,1700   
MPTEMP,11,2500  

MPDATA,KXX,2,,51.4
MPDATA,KXX,2,,47.2   
MPDATA,KXX,2,,36.8   
MPDATA,KXX,2,,30.25   
MPDATA,KXX,2,,25.31   
MPDATA,KXX,2,,20.2   
MPDATA,KXX,2,,33   
MPDATA,KXX,2,,35.2   
MPDATA,KXX,2,,36.5   
MPDATA,KXX,2,,38.9   
MPDATA,KXX,2,,40.3   


MPTEMP,,,,,,,,  
MPTEMP,1,20  
MPTEMP,2,100
MPTEMP,3,200
MPTEMP,4,300
MPTEMP,5,400
MPTEMP,6,500
MPTEMP,7,750
MPTEMP,8,1000
MPTEMP,9,1500
MPTEMP,10,1700   
MPTEMP,11,2500  

MPDATA,C,2,,480   
MPDATA,C,2,,485   
MPDATA,C,2,,498   
MPDATA,C,2,,506   
MPDATA,C,2,,526   
MPDATA,C,2,,539   
MPDATA,C,2,,675   
MPDATA,C,2,,670   
MPDATA,C,2,,660   
MPDATA,C,2,,780   
MPDATA,C,2,,820   

MPTEMP,,,,,,,,  
MPTEMP,1,20  
MPTEMP,2,100
MPTEMP,3,200
MPTEMP,4,300
MPTEMP,5,400
MPTEMP,6,500
MPTEMP,7,750
MPTEMP,8,1000
MPTEMP,9,1500
MPTEMP,10,1700   
MPTEMP,11,2500  

MPDATA,DENS,2,,7820
MPDATA,DENS,2,,7800
MPDATA,DENS,2,,7800
MPDATA,DENS,2,,7800
MPDATA,DENS,2,,7800
MPDATA,DENS,2,,7800
MPDATA,DENS,2,,7800
MPDATA,DENS,2,,7800
MPDATA,DENS,2,,7700
MPDATA,DENS,2,,7600
MPDATA,DENS,2,,7550

MPTEMP,,,,,,,,  
MPTEMP,1,20  
MPTEMP,2,100
MPTEMP,3,200
MPTEMP,4,300
MPTEMP,5,400
MPTEMP,6,500
MPTEMP,7,750
MPTEMP,8,1000
MPTEMP,9,1500
MPTEMP,10,1700   
MPTEMP,11,2500  

MPDATA,HF,2,,33.5
MPDATA,HF,2,,250
MPDATA,HF,2,,350
MPDATA,HF,2,,500
MPDATA,HF,2,,600
MPDATA,HF,2,,800
MPDATA,HF,2,,820
MPDATA,HF,2,,850
MPDATA,HF,2,,878
MPDATA,HF,2,,890
MPDATA,HF,2,,1000

MPTEMP,,,,,,,,  
MPTEMP,1,20  
MPTEMP,2,100
MPTEMP,3,200
MPTEMP,4,300
MPTEMP,5,400
MPTEMP,6,500
MPTEMP,7,750
MPTEMP,8,1000
MPTEMP,9,1500
MPTEMP,10,1700   
MPTEMP,11,2500  

MPDATA,ENTH,2,,345
MPDATA,ENTH,2,,737
MPDATA,ENTH,2,,1177
MPDATA,ENTH,2,,1664
MPDATA,ENTH,2,,2199
MPDATA,ENTH,2,,2779
MPDATA,ENTH,2,,3406
MPDATA,ENTH,2,,4079
MPDATA,ENTH,2,,4797
MPDATA,ENTH,2,,5561
MPDATA,ENTH,2,,6369

arc=0.023094  !0.025
theta1=60
theta2=90
x1=0
y1=0
x2=0.25
y2=0.02
lsize=0.004
la=0.05  !half of oringin

pi=3.1415926
I=135
U=15
effcient=0.85
r=0.010  
Qm=3*I*U*effcient/(pi*r**2)
kk=0.85
Qmax=Qm*kk
v=0.004 !速度
tm=lx/v !焊接总时间

PCIRC,arc, ,theta1,theta2,  
RECTNG,x1,x2,y1,y2,
FLST,2,2,5,ORDE,2   
FITEM,2,1   
FITEM,2,-2  
AOVLAP,P51X
/DIST,1,0.729,1
/REP,FAST   
/DIST,1,0.729,1
/REP,FAST   
/DIST,1,0.729,1
/REP,FAST   
/FOC,1,-0.3,,,1
/REP,FAST   
/FOC,1,-0.3,,,1
/REP,FAST   
/FOC,1,-0.3,,,1
/REP,FAST   
/FOC,1,-0.3,,,1
/REP,FAST   
/FOC,1,-0.3,,,1
/REP,FAST   
FLST,2,2,5,ORDE,2   
FITEM,2,3   
FITEM,2,-4  
AADD,P51X   
VOFFST,1,la, ,
APLOT   
!*  
VOFFST,5,la, ,
/VIEW,1,1,2,3
/ANG,1  
/REP,FAST     
/FOC,1,-0.3,,,1
/REP,FAST
/VIEW,  1, -0.246014646680    ,  0.221983578159    ,  0.943504151897
/ANG,   1, -0.584250030771  
/REPLO   
!*  
ET,2,SOLID70
keyopt,2,2,1
TYPE,   2   
MAT,       1
REAL,   
ESYS,       0   
SECNUM,
!*  
ESIZE,lsize,0,   
FLST,5,3,6,ORDE,2   
FITEM,5,1   
FITEM,5,-3  
CM,_Y,VOLU  
VSEL, , , ,P51X
CM,_Y1,VOLU
CHKMSH,'VOLU'   
CMSEL,S,_Y  
!*  
VSWEEP,_Y1  
!*  
CMDELE,_Y   
CMDELE,_Y1  
CMDELE,_Y2  
!*  
NUMMRG,ALL,,,,LOW
NUMCMP,ALL

! divide into two
*SET,nna,1
ESEL,ALL
*GET,emax,ELEM,,NUM,MAX
VSEL,S, , ,       nna
ESLV
*GET,nse,ELEM,,COUNT
*DIM,ne,,nse
*DIM,n1,,nse
*DIM,n2,,nse
*DIM,n3,,nse
*DIM,nex,,nse
*DIM,ney,,nse
*DIM,neorder,,nse
ii=0

*DO,i,1,emax
*IF,ESEL(i),EQ,1,THEN
ii=ii+1
ne(ii)=i
*ENDIF
*ENDDO

! nse belong to welding material.
*DO,i,1,nse
*GET,ney(i),ELEM,ne(i),CENT,Y
*GET,nex(i),ELEM,ne(i),CENT,X
*ENDDO

! making a up welding
d=0
*DO,i,1,nse
*IF,ney(i),GE,y2*2/3,THEN
d=d+1
n2(d)=ne(i)
*ENDIF
*ENDDO

! making a down
c=0
*DO,i,1,nse
*IF,ney(i),LT,y2/3,THEN
c=c+1
n1(c)=ne(i)
*ENDIF
*ENDDO

! making a middle
e=0
*DO,i,1,nse
*IF,ney(i),GE,y2/3,AND,ney(i),LT,y2*2/3,THEN
e=e+1
n3(e)=ne(i)
*ENDIF
*ENDDO

ESEL,NONE
*DO,i,1,d
*IF,n2(i),NE,0,THEN
ESEL,A,ELEM,,n2(i)
*ENDIF
*ENDDO
EPLOT

mine=0
*DIM,ne2,,d
*DIM,nex2,,d
*DIM,ney2,,d
*DO,i1,1,d
ESEL,U,ELEM,,mine
*GET,nse2,ELEM,,COUNT
ii=0
*DO,i,1,emax
*IF,ESEL(i),EQ,1,THEN
ii=ii+1
ne2(ii)=i
*ENDIF
*ENDDO
*DO,i,1,nse2
*GET,ney2(i),ELEM,ne2(i),CENT,Y
*GET,nex2(i),ELEM,ne2(i),CENT,X
*ENDDO

miny=1e20
minx=1e20

*DO,i,1,nse2
*IF,ney2(i),LT,miny,THEN
miny=ney2(i)
minx=nex2(i)
mine=ne2(i)
*ELSE
*IF,ney2(i),EQ,miny,THEN
*IF,nex2(i),LT,minx,THEN
miny=ney2(i)
minx=nex2(i)
mine=ne2(i)
*ENDIF
*ENDIF
*ENDIF
*ENDDO
neorder(i1)=mine
*ENDDO

maxe=0
ESEL,NONE
*DO,i,1,c
*IF,n1(i),NE,0,THEN
ESEL,A,ELEM,,n1(i)
*ENDIF
*ENDDO
EPLOT

*DIM,ne1,,c
*DIM,nex1,,c
*DIM,ney1,,c

*DO,i1,1,c
ESEL,U,ELEM,,maxe
*GET,nse3,ELEM,,COUNT
ii=0
*DO,i,1,emax
*IF,ESEL(i),EQ,1,THEN
ii=ii+1
ne1(ii)=i
*ENDIF
*ENDDO

*DO,i,1,nse3
*GET,ney1(i),ELEM,ne1(i),CENT,Y
*GET,nex1(i),ELEM,ne1(i),CENT,X
*ENDDO

maxy=-1e10
maxx=-1e10
*DO,i,1,nse3
*IF,ney1(i),GT,maxy,THEN
maxy=ney1(i)
maxx=nex1(i)
maxe=ne1(i)
*ELSE
*IF,ney1(i),EQ,maxy,THEN
*IF,nex1(i),LT,maxx,THEN
maxy=ney1(i)
maxx=nex1(i)
maxe=ne1(i)
*ENDIF
*ENDIF
*ENDIF
*ENDDO
neorder(i1+d)=maxe
*ENDDO

midd=0
ESEL,NONE
*DO,i,1,e
*IF,n3(i),NE,0,THEN
ESEL,A,ELEM,,n3(i)
*ENDIF
*ENDDO
EPLOT

*DIM,ne3,,c
*DIM,nex3,,c
*DIM,ney3,,c

*DO,i1,1,e
ESEL,U,ELEM,,maxe
*GET,nse1,ELEM,,COUNT
ii=0
*DO,i,1,emax
*IF,ESEL(i),EQ,1,THEN
ii=ii+1
ne3(ii)=i
*ENDIF
*ENDDO

*DO,i,1,nse1
*GET,ney3(i),ELEM,ne3(i),CENT,Y
*GET,nex3(i),ELEM,ne3(i),CENT,X
*ENDDO

midy=-1e20
midx=1e20
*DO,i,1,nse1
*IF,ney3(i),GT,midy,THEN
midy=ney3(i)
midx=nex3(i)
midd=ne3(i)
*ELSE
*IF,ney3(i),EQ,midy,THEN
*IF,nex3(i),LT,midx,THEN
midy=ney3(i)
midx=nex3(i)
midd=ne3(i)
*ENDIF
*ENDIF
*ENDIF
*ENDDO
neorder(i1+d+c)=midd
*ENDDO
ALLSEL,ALL

ESEL,ALL
EPLOT
/AUTO,1
/REPLOT
/DIST,1,0.729,1
/REP,FAST   
/DIST,1,0.729,1
/REP,FAST   
/FOC,1,-0.3,,,1
/REP,FAST   
/FOC,1,-0.3,,,1
/REP,FAST   
/FOC,1,-0.3,,,1
/REP,FAST  
/DIST,1,0.729,1
/REP,FAST   
/FOC,1,-0.3,,,1
/REP,FAST   
*DO,i,1,nse
EKILL,neorder(i)
ESEL,S,LIVE
EPLOT
*ENDDO
ALLSEL,ALL

/solu
u=d+1 !compute number
p=c
h=2000
dt1=1e-3
dt=5
t=0
ANTYPE,4
NSUBST,1
AUTOTS,1
LNSRCH,1
NLGEOM,1
NROPT,FULL,,ON !turn on large sharp


OUTRES,ALL,ALL
IC,ALL,TEMP,20
KBC,1
TIMINT,0,STRUCT
TIMINT,1,THERM
TIMINT,0,MAG
TINTP,0.005,,,1,0.5,0.2

nsub1=2
nsub2=20
*DO,i,nse,u,-1 !do down weld
EALIVE,neorder(i)
ESEL,S,LIVE
EPLOT
ESEL,ALL
t=t+dt1
TIME,t

*DO,j,1,8
D,NELEM(neorder(i),j),TEMP,h
*ENDDO

MPCHG,2,neorder(i) !change material number
SOLVE
t=t+dt1
TIME,t
SOLVE

*DO,j,1,8
DDELE,NELEM(neorder(i),j),TEMP
*ENDDO

t=t+dt-2*dt1
TIME,t
NSUBST,nsub1
SOLVE
*ENDDO

nsub1=2
nsub2=20
*DO,i,nse,u,-1 !do down weld
EALIVE,neorder(i)
ESEL,S,LIVE
EPLOT
ESEL,ALL
t=t+dt1
TIME,t

BFE,ALL,HGEN,1,h

MPCHG,2,neorder(i) !change material number
SOLVE
t=t+dt1
TIME,t
SOLVE

BFEDELE,ALL,HGEN

t=t+dt-2*dt1
TIME,t
NSUBST,nsub1
SOLVE
*ENDDO

FINISH

!
/config,nres,3000
/solu
o=d !compute number
h=2000
dt1=1e-3
dt=5

ANTYPE,4
NSUBST,1
AUTOTS,1
LNSRCH,1
NLGEOM,1

OUTRES,ALL,ALL
IC,ALL,TEMP,20
KBC,1
TIMINT,0,STRUCT
TIMINT,1,THERM
TIMINT,0,MAG
TINTP,0.005,,,1,0.5,0.2

nsub1=2
nsub2=20

*DO,i,1,o  !do up weld
EALIVE,neorder(i)
ESEL,S,LIVE
EPLOT
ESEL,ALL
t=t+dt1
TIME,t
*DO,j,1,8
D,NELEM(neorder(i),j),TEMP,h
*ENDDO
MPCHG,2,neorder(i)
SOLVE
t=t+dt1
TIME,t
SOLVE

*DO,j,1,8
DDELE,NELEM(neorder(i),j),TEMP
*ENDDO

t=t+dt-2*dt1
TIME,t
NSUBST,nsub1
SOLVE
*ENDDO

FINISH

tto

顶下 希望有人来讨论啊
qq：418456682

tto

顶下 希望有人来讨论啊
qq：418456682