多层焊接模拟分析热源加载求助
正在进行多层焊接模拟，采用单元生死，单元按生热率来加载，每个单元焊接作用时间为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
BF,NELEM(neorder(i),j),HGEN,h
*ENDDO
MPCHG,2,neorder(i) !change material number
SOLVE
t=t+dt1
TIME,t
SOLVE
*DO,j,1,8
BFDELE,NELEM(neorder(i),j),HGEN
*ENDDO
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
BF,NELEM(neorder(i),j),HGEN,h
*ENDDO
MPCHG,2,neorder(i)
SOLVE
t=t+dt1
TIME,t
SOLVE
*DO,j,1,8
BFDELE,NELEM(neorder(i),j),HGEN
*ENDDO
t=t+dt-2*dt1
TIME,t
NSUBST,nsub1
SOLVE
*ENDDO
FINISH
[ 本帖最后由 tto 于 2006-7-3 15:54 编辑 ] 顶下 希望有人来讨论啊
qq：418456682 做焊接的人实在是太少了，最近在做点焊，属于热电分析，和你的不太一样，都卡了好几天了， 一直也没碰到作过的同仁，真是同病相连啊。帮你顶一下 单元按生热率?
生热率要用hgen阿！
束上的曲线贴出来看看 hyd1130你是做熔化焊温度场和应力场分析的吗?
我做的是点焊的温度场和应力场分析,不知道你了解不?可否指导一下,急啊. 生热率要用hgen 我是直接定义了数值 然后加载上去的
*DO,j,1,8
BF,NELEM(neorder(i),j),HGEN,h
*ENDDO
是按照节点加载的。帮我分析下 看哪里出问题了 谢谢啊
qq：418456682
不好意思 以前贴的命令流有问题 现在我修改了 希望大家运行看下 大家怎不来讨论下嘛
我每天都来自己在顶了 放假没人来了吗？
我在顶下 在顶下 希望大家来讨论 在顶 在顶 一定要等有人来讨论 在顶 在顶 一定要等有人来讨论 顶下 希望大家来看 在顶下 顶下 每天都顶下 每天都顶下 每天都顶下 BF和BFE是有区别的，一般说来用BF结果温度要高一些你可以试试BFE，另在精华帖中有BF和BFE区别的介绍你也可以搜搜 谢谢 bf和bfe我都试过了 还是不行啊 焊接温度和冷却曲线 就是和经典的不一致
请大家在看看
页:
[1]
2