模拟焊接生死单元时出了问题,大家帮我看看是怎么回事(文章内是我的程序)

zxw19851015

我模拟焊接生死单元时,双椭球体不能行走,请大家帮我看看
/prep7
width=0.05
wwidth=0.0046
wdepth=0.008
wheight=0.001
length=0.05
K,1,wwidth,0,0,
K,2,width,0,0,   
K,3,width,0,-0.01,   
K,4,-width,0,-0.01,  
K,5,-width,0,0,  
K,6,-wwidth,0,0,
K,7,0,0,-wdepth,
K,8,0,0,0.001,  
LSTR,1,2  
LSTR,2,3  
LSTR,3,4  
LSTR,4,5  
LSTR,6,5  
LSTR,6,7  
LSTR,7,1  
LARC,6,1,8
AL,5,4,3,2,1,7,6
AL,6,7,8
VOFFST,1,-length,
VOFFST,2,-length,
wprot,0,0,90
VSBW,1
NUMCMP,ALL
wpoff,0,0,0.015
VSBW,3
wpoff,0,0,-0.03
VSBW,2
wpoff,0,0,-0.015
VSBW,3
wpoff,0,0,0.06
VSBW,4
NUMCMP,ALL
ET,1,SOLID70
MPTEMP,,,,,,,,  
MPTEMP,1,298
MPTEMP,2,398
MPTEMP,3,498
MPTEMP,4,598
MPTEMP,5,698
MPTEMP,6,798
MPTEMP,7,898
MPTEMP,8,998
MPTEMP,9,1098   
MPDATA,KXX,1,,31
MPDATA,KXX,1,,32.5  
MPDATA,KXX,1,,33.5  
MPDATA,KXX,1,,34
MPDATA,KXX,1,,33
MPDATA,KXX,1,,33.5  
MPDATA,KXX,1,,33.8  
MPDATA,KXX,1,,34
MPDATA,KXX,1,,30
MPTEMP,,,,,,,,  
MPTEMP,1,298
MPTEMP,2,398
MPTEMP,3,498
MPTEMP,4,598
MPTEMP,5,698
MPTEMP,6,798
MPTEMP,7,898
MPTEMP,8,998
MPTEMP,9,1098   
MPDATA,C,1,,410
MPDATA,C,1,,500
MPDATA,C,1,,530
MPDATA,C,1,,570
MPDATA,C,1,,600
MPDATA,C,1,,720
MPDATA,C,1,,810
MPDATA,C,1,,1080
MPDATA,C,1,,800
MPTEMP,,,,,,,,                    
MPTEMP,1,0                        
MPDATA,DENS,1,,7890               
MPTEMP,,,,,,,,  
MPTEMP,1,298
MPTEMP,2,398
MPTEMP,3,498
MPTEMP,4,598
MPTEMP,5,698
MPTEMP,6,798
MPTEMP,7,898
MPTEMP,8,998
MPTEMP,9,1098
MPDATA,HF,1,,100                  
MPDATA,HF,1,,350
MPDATA,HF,1,,480
MPDATA,HF,1,,520
MPDATA,HF,1,,1000   
MPDATA,HF,1,,1050   
MPDATA,HF,1,,1500   
MPDATA,HF,1,,2000   
MPDATA,HF,1,,3000   
MPTEMP,,,,,,,,  
MPTEMP,1,298
MPTEMP,2,398
MPTEMP,3,498
MPTEMP,4,598
MPTEMP,5,698
MPTEMP,6,798
MPTEMP,7,898
MPTEMP,8,998
MPTEMP,9,1098   
MPDATA,KXX,2,,31
MPDATA,KXX,2,,32.5  
MPDATA,KXX,2,,33.5  
MPDATA,KXX,2,,34
MPDATA,KXX,2,,33
MPDATA,KXX,2,,33.5  
MPDATA,KXX,2,,33.8  
MPDATA,KXX,2,,34
MPDATA,KXX,2,,30
MPTEMP,,,,,,,,  
MPTEMP,1,298
MPTEMP,2,398
MPTEMP,3,498
MPTEMP,4,598
MPTEMP,5,698
MPTEMP,6,798
MPTEMP,7,898
MPTEMP,8,998
MPTEMP,9,1098   
MPDATA,C,2,,410
MPDATA,C,2,,500
MPDATA,C,2,,530
MPDATA,C,2,,570
MPDATA,C,2,,600
MPDATA,C,2,,720
MPDATA,C,2,,810
MPDATA,C,2,,1080
MPDATA,C,2,,800
MPTEMP,,,,,,,,                    
MPTEMP,1,0                        
MPDATA,DENS,2,,7890               
MPTEMP,,,,,,,,  
MPTEMP,1,298
MPTEMP,2,398
MPTEMP,3,498
MPTEMP,4,598
MPTEMP,5,698
MPTEMP,6,798
MPTEMP,7,898
MPTEMP,8,998
MPTEMP,9,1098
MPDATA,HF,2,,100                  
MPDATA,HF,2,,350
MPDATA,HF,2,,480
MPDATA,HF,2,,520
MPDATA,HF,2,,1000   
MPDATA,HF,2,,1050   
MPDATA,HF,2,,1500   
MPDATA,HF,2,,2000   
MPDATA,HF,2,,3000
SAVE
MAT,1
LESIZE,3,,,10,
LESIZE,27,,,10,
LESIZE,33,,,10,
LESIZE,2,,,10,
LESIZE,8,,,6,
LESIZE,4,,,6,
LESIZE,11,,,6,
LESIZE,1,,,6,
EXTOPT,ACLEAR,1
EXTOPT,VSWE,TETS,1  
EXTOPT,VSWE,AUTO,0
EXTOPT,ESIZE,24,0
VSWEEP,2,2,12
VSWEEP,3,5,16
SAVE
MAT,2
LESIZE,5,,,10,
LESIZE,6,,,10,
LESIZE,7,,,10,
EXTOPT,ACLEAR,1
EXTOPT,VSWE,TETS,1  
EXTOPT,VSWE,AUTO,0
EXTOPT,ESIZE,40,0
VSWEEP,1,1,8
SAVE
VSEL,S, , ,4
MAT,1
LESIZE,38,,,10,
LESIZE,34,,,10,
LESIZE,40,,,10,
LESIZE,30,,,4,
EXTOPT,ESIZE,40,0
VSWEEP,4,22,23
VSEL,S, , ,5
MAT,1
LESIZE,46,,,10,
LESIZE,32,,,10,
LESIZE,48,,,10,
EXTOPT,ESIZE,40,0
VSWEEP,5,27,28
VSEL,S, , ,6
MAT,1
LESIZE,54,,,6,
LESIZE,55,,,6,
LESIZE,56,,,6,
LESIZE,57,,,6,
mshkey,1
mshape,0,3d
vmesh,6
VSEL,S, , ,7
MAT,1
LESIZE,62,,,6,
LESIZE,63,,,6,
LESIZE,64,,,6,
LESIZE,65,,,6,
mshkey,1
mshape,0,3d
vmesh,7
save
WPSTYLE,,,,,,,,0
D=0.0092           
ah=D/2
chf=3*ah
chb=8*ah
bh=0.009
rf=2*chf/(chf+chb)
rb=2*chb/(chf+chb)
N=0.27                    
W0=2200                  
W=W0*N                     
ALLSEL,ALL  
VSEL,S,,,1      
ESLV,R  
cm,welddeath,elem
ekill,welddeath
ALLSEL,ALL

STEP=0.00125               
N=length/STEP
V=0.004                  
*DO,I,0,N
YY=I*STEP
T=(YY+STEP)/V
/SOLU
ANTYPE,4
TRNOPT,FULL
TUNIF,298,
OUTRES,NSOL,all
TIME,T
NROPT,FULL,,ON
AUTOTS,1
NSUBST,5, , ,1   !!!
KBC,1
lnsrch,1     !!!
TSRES,ERASE
TINTP,,,,1,,,,,!!!
TIMINT,0,struct
TIMINT,1,THERM
TIMINT,0,MAG
sfa,14,1,conv,200,298
sfa,32,1,conv,200,298
sfa,25,1,conv,200,298
sfa,20,1,conv,200,298
sfa,18,1,conv,200,298
sfa,38,1,conv,200,298
sfa,2,1,conv,200,298
sfa,30,1,conv,200,298
sfa,27,1,conv,200,298
sfa,1,1,conv,200,298
sfa,22,1,conv,200,298
sfa,35,1,conv,200,298
sfa,5,1,conv,200,298
sfa,12,1,conv,200,298
sfa,31,1,conv,200,298
sfa,28,1,conv,200,298
sfa,8,1,conv,200,298
sfa,23,1,conv,200,298
sfa,36,1,conv,200,298
sfa,16,1,conv,200,298
sfa,17,1,conv,200,298
sfa,37,1,conv,200,298
sfa,21,1,conv,200,298
sfa,26,1,conv,200,298
sfa,33,1,conv,200,298
sfa,15,1,conv,200,298
sfa,3,1,conv,200,298
sfa,4,1,conv,200,298
LOCAL,15,2,0,YY,wheight, , , ,chf/ah,bh/ah,  
NSEL,S,LOC,X,-ah,ah
NSEL,R,LOC,Y,,180
NSEL,R,LOC,Z,-90,0
ESLN,S
CM,headlive,ELEM
EALIVE,headlive
ALLSEL,ALL
LOCAL,16,2,0,YY,wheight, , , ,chb/ah,bh/ah,   
NSEL,S,LOC,X,-ah,ah   
NSEL,R,LOC,Y,180,360
NSEL,R,LOC,Z,-90,0
ESLN,S
CM,behindlive,ELEM
EALIVE,behindlive
ALLSEL,ALL
CMDELE,headlive
CMDELE,behindlive
ALLSEL,ALL
!!!!!!添加各边对流换热条件
NSEL,S,EXT              
NSEL,U,LOC,X,0           
SF,ALL,CONV,200,298                 
ALLSEL,ALL
!!!双椭球热源的加载
!!!前半球节点的选取
LOCAL,11,2,0,YY,wheight, , , ,chf/ah,bh/ah,  
NSEL,S,LOC,X,-ah,ah
NSEL,R,LOC,Y,,180
NSEL,R,LOC,Z,-90,0
CM,FRONT,NODE
ALLSEL,ALL
!!!!!!!!!!!!!!!!!!!!!!!!!!!
LOCAL,13,0,0,YY,wheight, , , ,1,1,
!!!前半球函数、表格定义
*DEL,_FNCNAME   
*DEL,_FNCMTID   
*DEL,_FNC_C1
*DEL,_FNC_C2
*DEL,_FNC_C3
*DEL,_FNC_C4
*DEL,_FNC_C5
*DEL,_FNCCSYS
*DEL,'TFRONT'    !!!!!!!!
*SET,_FNCNAME,'TFRONT'  
*DIM,_FNC_C1,,1
*DIM,_FNC_C2,,1
*DIM,_FNC_C3,,1
*DIM,_FNC_C4,,1
*DIM,_FNC_C5,,1  
*SET,_FNC_C1(1),rf           !!!rf
*SET,_FNC_C2(1),w            !!!Q
*SET,_FNC_C3(1),chf          !!!chf
*SET,_FNC_C4(1),ah           !!!ah
*SET,_FNC_C5(1),bh           !!!bh
*SET,_FNCCSYS,13            
*DIM,%_FNCNAME%,TABLE,7,38,1,,,,%_FNCCSYS%  
*SET,%_FNCNAME%(0,0,1), 0.0, -999   
*SET,%_FNCNAME%(2,0,1), 0.0
*SET,%_FNCNAME%(3,0,1), %_FNC_C1(1)%
*SET,%_FNCNAME%(4,0,1), %_FNC_C2(1)%
*SET,%_FNCNAME%(5,0,1), %_FNC_C3(1)%
*SET,%_FNCNAME%(6,0,1), %_FNC_C4(1)%
*SET,%_FNCNAME%(7,0,1), %_FNC_C5(1)%
*SET,%_FNCNAME%(0,1,1), 1.0, -1, 16, 3, 0, 0, 0
*SET,%_FNCNAME%(0,2,1), 0.0, -2, 0, 6, 0, 0, -1
*SET,%_FNCNAME%(0,3,1),   0, -3, 0, 1, -2, 3, -1
*SET,%_FNCNAME%(0,4,1), 0.0, -1, 0, 1, -3, 3, 17
*SET,%_FNCNAME%(0,5,1), 0.0, -2, 0, 1, -1, 3, 18
*SET,%_FNCNAME%(0,6,1), 0.0, -1, 0, 0, 0, 0, 0  
*SET,%_FNCNAME%(0,7,1), 0.0, -3, 0, 1, 0, 0, -1
*SET,%_FNCNAME%(0,8,1), 0.0, -4, 0, 1, -1, 2, -3
*SET,%_FNCNAME%(0,9,1), 0.0, -1, 0, 2, 0, 0, 2  
*SET,%_FNCNAME%(0,10,1), 0.0, -3, 0, 1, 2, 17, -1   
*SET,%_FNCNAME%(0,11,1), 0.0, -1, 0, 2, 0, 0, 19
*SET,%_FNCNAME%(0,12,1), 0.0, -5, 0, 1, 19, 17, -1  
*SET,%_FNCNAME%(0,13,1), 0.0, -1, 0, 1, -3, 4, -5   
*SET,%_FNCNAME%(0,14,1), 0.0, -3, 0, 2, 0, 0, 3
*SET,%_FNCNAME%(0,15,1), 0.0, -5, 0, 1, 3, 17, -3   
*SET,%_FNCNAME%(0,16,1), 0.0, -3, 0, 2, 0, 0, 20
*SET,%_FNCNAME%(0,17,1), 0.0, -6, 0, 1, 20, 17, -3  
*SET,%_FNCNAME%(0,18,1), 0.0, -3, 0, 1, -5, 4, -6   
*SET,%_FNCNAME%(0,19,1), 0.0, -5, 0, 1, -1, 1, -3   
*SET,%_FNCNAME%(0,20,1), 0.0, -1, 0, 2, 0, 0, 4
*SET,%_FNCNAME%(0,21,1), 0.0, -3, 0, 1, 4, 17, -1   
*SET,%_FNCNAME%(0,22,1), 0.0, -1, 0, 2, 0, 0, 21
*SET,%_FNCNAME%(0,23,1), 0.0, -6, 0, 1, 21, 17, -1   
*SET,%_FNCNAME%(0,24,1), 0.0, -1, 0, 1, -3, 4, -6   
*SET,%_FNCNAME%(0,25,1), 0.0, -3, 0, 1, -5, 1, -1   
*SET,%_FNCNAME%(0,26,1), 0.0, -1, 0, 3, 0, 0, -3
*SET,%_FNCNAME%(0,27,1), 0.0, -5, 0, 1, -1, 3, -3   
*SET,%_FNCNAME%(0,28,1), 0.0, -1, 0, 1, -4, 3, -5   
*SET,%_FNCNAME%(0,29,1), 0.0, -1, 7, 1, -1, 0, 0
*SET,%_FNCNAME%(0,30,1), 0.0, -3, 0, 1, -2, 3, -1   
*SET,%_FNCNAME%(0,31,1), 0.0, -1, 0, 1, 20, 3, 21   
*SET,%_FNCNAME%(0,32,1), 0.0, -2, 0, 1, -1, 3, 19   
*SET,%_FNCNAME%(0,33,1), 0.0, -1, 0, 3.14159265358979310, 0, 0, -2  
*SET,%_FNCNAME%(0,34,1), 0.0, -4, 0, 1, -2, 3, -1   
*SET,%_FNCNAME%(0,35,1), 0.0, -1, 16, 3.14159265358979310, 0, 0, 0  
*SET,%_FNCNAME%(0,36,1), 0.0, -2, 0, 1, -4, 3, -1   
*SET,%_FNCNAME%(0,37,1), 0.0, -1, 0, 1, -3, 4, -2   
*SET,%_FNCNAME%(0,38,1), 0.0, 99, 0, 1, -1, 0, 0
!!!!!!!!!!!!!!!!!!
BF,FRONT,HGEN, %TFRONT%     
CMDELE,FRONT   
ALLSEL,ALL
!!!后半球节点的选取
LOCAL,12,2,0,YY,wheight, , , ,chb/ah,bh/ah,   
NSEL,S,LOC,X,-ah,ah   
NSEL,R,LOC,Y,180,360
NSEL,R,LOC,Z,-90,0
CM,BEHIND,NODE
!!!后半球加载
!!!!!!!!!!!!!!!!!!!
LOCAL,14,0,0,YY,wheight, , , ,1,1,
!!!后半球函数、表格定义
*DEL,_FNCNAME   
*DEL,_FNCMTID   
*DEL,_FNC_C1
*DEL,_FNC_C2
*DEL,_FNC_C3
*DEL,_FNC_C4
*DEL,_FNC_C5
*DEL,_FNCCSYS
*DEL,'TBEHIND'  !!!!!!!!!
*SET,_FNCNAME,'TBEHIND'
*DIM,_FNC_C1,,1
*DIM,_FNC_C2,,1
*DIM,_FNC_C3,,1
*DIM,_FNC_C4,,1
*DIM,_FNC_C5,,1
*SET,_FNC_C1(1),rb                   !!!rb
*SET,_FNC_C2(1),w                   !!!Q
*SET,_FNC_C3(1),chb                !!!chb
*SET,_FNC_C4(1),ah                 !!!ah
*SET,_FNC_C5(1),bh                 !!!bh
*SET,_FNCCSYS,14                        !!!   0
*DIM,%_FNCNAME%,TABLE,7,38,1,,,,%_FNCCSYS%  
*SET,%_FNCNAME%(0,0,1), 0.0, -999   
*SET,%_FNCNAME%(2,0,1), 0.0
*SET,%_FNCNAME%(3,0,1), %_FNC_C1(1)%
*SET,%_FNCNAME%(4,0,1), %_FNC_C2(1)%
*SET,%_FNCNAME%(5,0,1), %_FNC_C3(1)%
*SET,%_FNCNAME%(6,0,1), %_FNC_C4(1)%
*SET,%_FNCNAME%(7,0,1), %_FNC_C5(1)%
*SET,%_FNCNAME%(0,1,1), 1.0, -1, 16, 3, 0, 0, 0
*SET,%_FNCNAME%(0,2,1), 0.0, -2, 0, 6, 0, 0, -1
*SET,%_FNCNAME%(0,3,1),   0, -3, 0, 1, -2, 3, -1
*SET,%_FNCNAME%(0,4,1), 0.0, -1, 0, 1, -3, 3, 17
*SET,%_FNCNAME%(0,5,1), 0.0, -2, 0, 1, -1, 3, 18
*SET,%_FNCNAME%(0,6,1), 0.0, -1, 0, 0, 0, 0, 0  
*SET,%_FNCNAME%(0,7,1), 0.0, -3, 0, 1, 0, 0, -1
*SET,%_FNCNAME%(0,8,1), 0.0, -4, 0, 1, -1, 2, -3
*SET,%_FNCNAME%(0,9,1), 0.0, -1, 0, 2, 0, 0, 2  
*SET,%_FNCNAME%(0,10,1), 0.0, -3, 0, 1, 2, 17, -1   
*SET,%_FNCNAME%(0,11,1), 0.0, -1, 0, 2, 0, 0, 19
*SET,%_FNCNAME%(0,12,1), 0.0, -5, 0, 1, 19, 17, -1  
*SET,%_FNCNAME%(0,13,1), 0.0, -1, 0, 1, -3, 4, -5   
*SET,%_FNCNAME%(0,14,1), 0.0, -3, 0, 2, 0, 0, 3
*SET,%_FNCNAME%(0,15,1), 0.0, -5, 0, 1, 3, 17, -3   
*SET,%_FNCNAME%(0,16,1), 0.0, -3, 0, 2, 0, 0, 20
*SET,%_FNCNAME%(0,17,1), 0.0, -6, 0, 1, 20, 17, -3  
*SET,%_FNCNAME%(0,18,1), 0.0, -3, 0, 1, -5, 4, -6   
*SET,%_FNCNAME%(0,19,1), 0.0, -5, 0, 1, -1, 1, -3   
*SET,%_FNCNAME%(0,20,1), 0.0, -1, 0, 2, 0, 0, 4
*SET,%_FNCNAME%(0,21,1), 0.0, -3, 0, 1, 4, 17, -1   
*SET,%_FNCNAME%(0,22,1), 0.0, -1, 0, 2, 0, 0, 21
*SET,%_FNCNAME%(0,23,1), 0.0, -6, 0, 1, 21, 17, -1  
*SET,%_FNCNAME%(0,24,1), 0.0, -1, 0, 1, -3, 4, -6   
*SET,%_FNCNAME%(0,25,1), 0.0, -3, 0, 1, -5, 1, -1   
*SET,%_FNCNAME%(0,26,1), 0.0, -1, 0, 3, 0, 0, -3
*SET,%_FNCNAME%(0,27,1), 0.0, -5, 0, 1, -1, 3, -3   
*SET,%_FNCNAME%(0,28,1), 0.0, -1, 0, 1, -4, 3, -5   
*SET,%_FNCNAME%(0,29,1), 0.0, -1, 7, 1, -1, 0, 0
*SET,%_FNCNAME%(0,30,1), 0.0, -3, 0, 1, -2, 3, -1   
*SET,%_FNCNAME%(0,31,1), 0.0, -1, 0, 1, 20, 3, 21   
*SET,%_FNCNAME%(0,32,1), 0.0, -2, 0, 1, -1, 3, 19   
*SET,%_FNCNAME%(0,33,1), 0.0, -1, 0, 3.14159265358979310, 0, 0, -2  
*SET,%_FNCNAME%(0,34,1), 0.0, -4, 0, 1, -2, 3, -1   
*SET,%_FNCNAME%(0,35,1), 0.0, -1, 16, 3.14159265358979310, 0, 0, 0  
*SET,%_FNCNAME%(0,36,1), 0.0, -2, 0, 1, -4, 3, -1   
*SET,%_FNCNAME%(0,37,1), 0.0, -1, 0, 1, -3, 4, -2   
*SET,%_FNCNAME%(0,38,1), 0.0, 99, 0, 1, -1, 0, 0
BF,BEHIND,HGEN, %TBEHIND%   
CMDELE,BEHIND
ALLSEL,ALL
SOLVE
BFDELE,all,HGEN
*ENDDO
SAVE

zxw19851015

请大家帮我解决这个问题，捆饶了一段时间了，我的QQ1003460493

半个苹果

期待 热心人出现   遇到同样问题 QQ 80976658

bijiaoshuai

期待学习中

spiderman216

期待高手出现

lingke125

热源移动可以采用建立局部坐标系的方法，局部坐标移动就可以实现啦。
你这个是直接建立载荷表么？

crruose

期待解决  呼唤高手

crruose

期待解决  呼唤高手

sd9366

期待高手解决

chutianyao

期待高手出现

wxpand2004

期待高手出现解决问题！

sd9366

期待高手

wxpand2004

谁要解决了期望快点贴出来！