求教版主iomega和高手们一个瞬态热分析不收敛的问题

wnqiangjingjiu

2.822575e-9

solution not converged at time 1.0583333e-04(load step 1 substep 1),run termimated

求教版主一个瞬态热分析不收敛的问题
模型如下：三明治结构；BLOCK1相当于为移动热源，Block2、Block3为焊料层和结构件，它们之间存在接触热阻为2e4。单位制采用的是国际单位制，当体热源加到1e12 J/m^3时，出现不收敛的问题
“solution not converged at time 1.0583333e-04(load step 1 substep 1),run termimated”
请问怎么解决这个问题啊？
附上命令流：
/FILNAME,EX5-4,1                             ! set filename
/TITLE,The Thermal Analysis of the Welding   ! set title  
/PREP7                                       !  
!******************************
!   
!******************************
*SET,LENGTH1,12.7e-3                         !
*SET,LENGTH2,12.7e-3
*SET,LENGTH3,25.4e-3
*SET,WIDTH,6.35e-3  
*SET,HEIGHT1,35e-6  
*SET,HEIGHT2,25e-6  
*SET,HEIGHT3,5e-4
*SET,STEP,12.7e-4
*SET,COUNT,LENGTH1/STEP    ! count=20
*SET,V,6                    !stept=1e-4
!******************************
!   
!******************************
ET,1,SOLID70
!******************************
MPTEMP,,,,,,,,                       
MPTEMP,1,0                        
MPDATA,KXX,1,,160              
!定义比热容1
MPTEMP,,,,,,,,                 
MPTEMP,1,0                  
MPDATA,C,1,,830
!定义密度
MPTEMP,,,,,,,,                    
MPTEMP,1,0                        
MPDATA,DENS,1,,5550              

!!*************************************************************  ! PbSn  conductivity of PbSn
mptemp                                                           
mptemp,1,          290,    300,    310,    320,    330,   340  
mptemp,7,          350,    360,    370,   380
mpdata,kxx,2,1,    69, 70.71, 72.39, 74.05,  75.69, 77.32
mpdata,kxx,2,7, 78.93, 80.51, 82.09, 83.64
!!***********************                  !density of PbSn
!!mptemp,        1,   293,   300,   310,   320,  340,  350
!!mptemp,        7,   360,   370,  380,  390, 400, 410
!!mptemp,       13,  420,  430,  440,  450, 460, 470
!!mptemp,       19,  480,  490,  500
!!mpdata,dens,2,  1,   1.47e4,  1.4692e4,  1.4689e4, 1.4682e4, 1.4676e4, 1.4662e4
!!mpdata,dens,2,  7, 1.4654e4,  1.4647e4,  1.4639e4, 1.4632e4, 1.4624e4, 1.4615e4
!!mpdata,dens,2, 13, 1.4607e4,  1.4598e4,  1.4590e4, 1.4583e4, 1.4576e4, 1.4565e4
!!mpdata,dens,2, 19, 1.4554e4,  1.4545e4,  1.4534e4
!!*************************************************                   !material 2-specific heat of PbSn
mptemp
mptemp,  1,   298,   300,   350,   400,  450,  500   
mpdata,c,2, 1, 1.47e2,  1.47e2,  1.51e2, 1.53e2, 1.55e2, 1.57e2

!!*************************************************                       ! material 3-heat conductivity of Si
MPTEMP,,,,,,,,                      !   
MPTEMP,1,293                         !   
MPTEMP,2,373
MPTEMP,3,423
MPTEMP,4,473
MPTEMP,5,523
MPTEMP,6,573
MPTEMP,7,623
MPTEMP,8,673
MPTEMP,9,773
MPTEMP,10,873
MPTEMP,11,973
MPTEMP,12,1073
MPTEMP,13,1173
MPTEMP,14,1273
MPDATA,KXX,3,,929                !
MPDATA,KXX,3,,929  
MPDATA,KXX,3,,422  
MPDATA,KXX,3,,264  
MPDATA,KXX,3,,191
MPDATA,KXX,3,,148  
MPDATA,KXX,3,,119  
MPDATA,KXX,3,,98.9  
MPDATA,KXX,3,,76.2  
MPDATA,KXX,3,,61.9  
MPDATA,KXX,3,,50.8  
MPDATA,KXX,3,,42.2  
MPDATA,KXX,3,,35.9
MPDATA,KXX,3,,31.2
!!*************************************************                     
MPTEMP,,,,,,,,                     !
MPTEMP,1,293                       !
MPTEMP,2,303
MPTEMP,3,313
MPTEMP,4,323
MPTEMP,5,333
MPTEMP,6,343   
MPTEMP,7,353   
MPTEMP,8,363   
MPTEMP,9,373   
MPTEMP,10,398  
MPTEMP,11,423
MPTEMP,12,448
MPTEMP,13,473  
MPTEMP,14,498
MPTEMP,15,523
MPTEMP,16,546  
MPTEMP,17,571
MPTEMP,18,623
MPTEMP,19,673  
MPTEMP,20,773
MPTEMP,21,873
MPTEMP,22,973
MPTEMP,23,1073
MPTEMP,24,1173
MPTEMP,25,1273
MPDATA,C,3,,3.37                  
MPDATA,C,3,,17.1
MPDATA,C,3,,44
MPDATA,C,3,,78.5
MPDATA,C,3,,115
MPDATA,C,3,,151
MPDATA,C,3,,188
MPDATA,C,3,,224
MPDATA,C,3,,260
MPDATA,C,3,,346
MPDATA,C,3,,425
MPDATA,C,3,,496
MPDATA,C,3,,557
MPDATA,C,3,,607
MPDATA,C,3,,649
MPDATA,C,3,,681
MPDATA,C,3,,712
MPDATA,C,3,,759
MPDATA,C,3,,793
MPDATA,C,3,,840
MPDATA,C,3,,869
MPDATA,C,3,,892
MPDATA,C,3,,912
MPDATA,C,3,,930
MPDATA,C,3,,946
!!*************************************************                     
MPTEMP,,,,,,,,  
MPTEMP,1,0  
MPDATA,DENS,3,,2330
!******************************
!   
!******************************
WPSTYLE,,,,,,,,0                     !  
BLOCK,0,LENGTH1,0,WIDTH,0,HEIGHT1,   !
BLOCK,0,LENGTH2,0,WIDTH,0,-HEIGHT2, !
BLOCK,0,LENGTH3,0,WIDTH,-HEIGHT2,-HEIGHT3-HEIGHT2, !
/VIEW,1,1,1,1   
/ANG,1  
/REP,FAST   
/VIEW,1,,-1
/ANG,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   
/DIST,1,0.729,1
/REP,FAST   
/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   
/DIST,1,0.729,1
/REP,FAST   
/DIST,1,0.729,1
/REP,FAST   
FLST,5,1,4,ORDE,1   
FITEM,5,1   
CM,_Y,LINE  
LSEL, , , ,P51X
CM,_Y1,LINE
CMSEL,,_Y   
!*  
FLST,5,3,4,ORDE,3   
FITEM,5,1   
FITEM,5,13  
FITEM,5,25  
CM,_Y,LINE  
LSEL, , , ,P51X
CM,_Y1,LINE
CMSEL,,_Y   
!*  
LESIZE,_Y1,1.27e-4, , , , , , ,1  
!*  
FLST,5,3,4,ORDE,3   
FITEM,5,4   
FITEM,5,16  
FITEM,5,28  
CM,_Y,LINE  
LSEL, , , ,P51X
CM,_Y1,LINE
CMSEL,,_Y   
!*  
LESIZE,_Y1,1.27e-4, , , , , , ,1  
!*  

!*  
FLST,5,3,4,ORDE,1  
FITEM,5,33  
CM,_Y,LINE  
LSEL, , , ,P51X
CM,_Y1,LINE
CMSEL,,_Y   
!*  
LESIZE,_Y1,1e-4, , , , , , ,1
!*  
!*  
FLST,5,3,4,ORDE,2   
FITEM,5,9   
FITEM,5,21   
CM,_Y,LINE  
LSEL, , , ,P51X
CM,_Y1,LINE
CMSEL,,_Y   
!*  
LESIZE,_Y1,5e-6, , , , , , ,1
CM,_Y,VOLU  
VSEL, , , ,       1
CM,_Y1,VOLU
CMSEL,S,_Y  
!*  
CMSEL,S,_Y1
VATT,       1, ,   1,       0   
CMSEL,S,_Y  
CMDELE,_Y   
CMDELE,_Y1  
!*  
CM,_Y,VOLU  
VSEL, , , ,       1
CM,_Y1,VOLU
CHKMSH,'VOLU'   
CMSEL,S,_Y  
!*  
MSHAPE,0,3d
MSHKEY,1
VMESH,_Y1   
MSHKEY,0
!*  
CMDELE,_Y   
CMDELE,_Y1  
CMDELE,_Y2  
!*  
CM,_Y,VOLU  
VSEL, , , ,       2
CM,_Y1,VOLU
CMSEL,S,_Y  
!*  
CMSEL,S,_Y1
VATT,       2, ,   1,       0   
CMSEL,S,_Y  
CMDELE,_Y   
CMDELE,_Y1  
!*  
CM,_Y,VOLU  
VSEL, , , ,       2
CM,_Y1,VOLU
CHKMSH,'VOLU'   
CMSEL,S,_Y  
!*  
MSHAPE,0,3d
MSHKEY,1
VMESH,_Y1   
MSHKEY,0
!*  
CMDELE,_Y   
CMDELE,_Y1  
CMDELE,_Y2  
!*  
CM,_Y,VOLU  
VSEL, , , ,       3
CM,_Y1,VOLU
CMSEL,S,_Y  
!*  
CMSEL,S,_Y1
VATT,       3, ,   1,       0   
CMSEL,S,_Y  
CMDELE,_Y   
CMDELE,_Y1  
!*  
CM,_Y,VOLU  
VSEL, , , ,       3
CM,_Y1,VOLU
CHKMSH,'VOLU'   
CMSEL,S,_Y  
!*  
MSHAPE,0,3d
MSHKEY,1
VMESH,_Y1   
MSHKEY,0
!*  
CMDELE,_Y   
CMDELE,_Y1  
CMDELE,_Y2  
!*  
/REPLOT,RESIZE  
!*  
!*  
/COM, CONTACT PAIR CREATION - START
CM,_NODECM,NODE
CM,_ELEMCM,ELEM
CM,_KPCM,KP
CM,_LINECM,LINE
CM,_AREACM,AREA
CM,_VOLUCM,VOLU
/GSAV,cwz,gsav,,temp
MP,MU,1,
MAT,1   
R,3
REAL,3  
ET,2,170
ET,3,174
KEYOPT,3,9,0
KEYOPT,3,10,2   
R,3,
RMORE,  
RMORE,,2e4  
RMORE,0
KEYOPT,3,1,2
! Generate the target surface   
ASEL,S,,,1  
CM,_TARGET,AREA
TYPE,2  
NSLA,S,1
ESLN,S,0
ESLL,U  
ESEL,U,ENAME,,188,189   
ESURF   
CMSEL,S,_ELEMCM
! Generate the contact surface  
ASEL,S,,,8  
CM,_CONTACT,AREA
TYPE,3  
NSLA,S,1
ESLN,S,0
ESURF   
ALLSEL  
ESEL,ALL
ESEL,S,TYPE,,2  
ESEL,A,TYPE,,3  
ESEL,R,REAL,,3  
/PSYMB,ESYS,1   
/PNUM,TYPE,1
/NUM,1  
EPLOT   
ESEL,ALL
ESEL,S,TYPE,,2  
ESEL,A,TYPE,,3  
ESEL,R,REAL,,3  
CMSEL,A,_NODECM
CMDEL,_NODECM   
CMSEL,A,_ELEMCM
CMDEL,_ELEMCM   
CMSEL,S,_KPCM   
CMDEL,_KPCM
CMSEL,S,_LINECM
CMDEL,_LINECM   
CMSEL,S,_AREACM
CMDEL,_AREACM   
CMSEL,S,_VOLUCM
CMDEL,_VOLUCM   
/GRES,cwz,gsav  
CMDEL,_TARGET   
CMDEL,_CONTACT  
/COM, CONTACT PAIR CREATION - END   
!*  
!*  
/COM, CONTACT PAIR CREATION - START
CM,_NODECM,NODE
CM,_ELEMCM,ELEM
CM,_KPCM,KP
CM,_LINECM,LINE
CM,_AREACM,AREA
CM,_VOLUCM,VOLU
/GSAV,cwz,gsav,,temp
MP,MU,1,0   
MAT,1   
R,4
REAL,4  
ET,4,170
ET,5,175
KEYOPT,5,9,0
KEYOPT,5,10,2   
R,4,
RMORE,  
RMORE,,2e4  
RMORE,0
KEYOPT,5,1,2
! Generate the target surface   
ASEL,S,,,7  
CM,_TARGET,AREA
TYPE,4  
NSLA,S,1
ESLN,S,0
ESLL,U  
ESEL,U,ENAME,,188,189   
ESURF   
CMSEL,S,_ELEMCM
! Generate the contact surface  
NSEL,S,LOC,X,0,12.7e-3  
NSEL,R,LOC,Y,0,6.35e-3  
NSEL,R,LOC,Z,-25e-6


CM,_CONTACT,NODE
TYPE,5  
ESLN,S,0
ESURF   
ALLSEL  
ESEL,ALL
ESEL,S,TYPE,,4  
ESEL,A,TYPE,,5  
ESEL,R,REAL,,4  
/PSYMB,ESYS,1   
/PNUM,TYPE,1
/NUM,1  
EPLOT   
ESEL,ALL
ESEL,S,TYPE,,4  
ESEL,A,TYPE,,5  
ESEL,R,REAL,,4  
CMSEL,A,_NODECM
CMDEL,_NODECM   
CMSEL,A,_ELEMCM
CMDEL,_ELEMCM   
CMSEL,S,_KPCM   
CMDEL,_KPCM
CMSEL,S,_LINECM
CMDEL,_LINECM   
CMSEL,S,_AREACM
CMDEL,_AREACM   
CMSEL,S,_VOLUCM
CMDEL,_VOLUCM   
/GRES,cwz,gsav  
CMDEL,_TARGET   
CMDEL,_CONTACT  
/COM, CONTACT PAIR CREATION - END   
ALLSEL,ALL  
/MREP,EPLOT
FINISH                              !   
!******************************
!   
!******************************
!   
*DO,I,0,COUNT                        !
*SET,DISY,I*STEP                           !
*SET,T,(DISY+STEP)/V                       !
/SOLU  
                              !
ANTYPE,4                              !
TRNOPT,FULL                           !     set full transient analysis
TUNIF,293,                            !      
OUTRES,BASIC,LAST                     !     output last-substep basic result
TIME,T                                !
NROPT,FULL,,ON                        !      set N-R option
DELTIM,2E-6,1E-6,2E-6,1                             !      open auto time step
NSUBST,10,,,1                          !      set substep=10
KBC,1                                 !
LNSRCH,1                       !             open linear search
TSRES,ERASE                           !       clear time-step option
TINTP,,,,1,,,,,
TIMINT,1,THERM                        !        set time intergetion
LNSRCH,AUTO
SOLCONTROL,ON

!*******************************
!   
!*******************************
NSEL,S,EXT                     !
NSEL,U,LOC,Z,35E-6                 !  
SF,ALL,CONV,35,293                   !
ALLSEL,ALL                            !
!*******************************
   
!*******************************
LOCAL,11,0,DISY,0,0,       !
NSEL,S,LOC,X,0,12.7e-4                   !  
NSEL,R,LOC,Y,0,WIDTH                      !
NSEL,R,LOC,Z,0,HEIGHT1                      !                       !  
BF,ALL,HGEN,1.468886598e12               !  
CMDELE,NODEF                        !
ALLSEL,ALL                          !   
SOLVE                                   !   
BFDELE,ALL,HGEN                     !   
*ENDDO                                  !   
SAVE                                    !   
!***************************
!   
!***************************
CSYS,0
NSEL,S,EXT
NSEL,U,LOC,Z,35E-6   
SF,ALL,CONV,35,293
ALLSEL,ALL
TIME,LENGTH1/V+0.2                                  !  
NSUBST,10,,,1   
DELTIM,1E-2,5E-3,1E-2,1                         !
SOLVE
CSYS,0
ALLSEL,ALL
SAVE                                   !  

FINISH  
/POST26
nsol,9,node(6e-3,0,-125e-6),temp,,
nsol,2,1775,temp,,                       !
nsol,3,1278,temp,,                       !  
nsol,4,1279,temp,,                       !
nsol,5,1374,temp,,                       !
nsol,6,1375,temp,,                       !  
nsol,7,2204,temp,,
nsol,8,2022,temp,,                       !  
plvar,2,3,4,5,6,7,8,9                    !

wnqiangjingjiu

希望高手们能帮我看一下，怎样解决不收敛的问题

iomega

Simple way to check your model:
(1) run static simulation case.
(2) change boundary conditions, e.g., isothermal, or set much large thermal contact conductance...
BTW, the 2e4 value is the thermal contact conductance value, not  thermal contact resistance value

tonnyw

Your model is too complicated. Why not remove the contact element and run a simple 2D linear model first? If everything goes well, add nonlinearity to the 2D model.

wnqiangjingjiu

3# iomega

to iomega
sorry! i have made a mistake!
the 2e4 value indeed is the themal contact conductance,not a themal contact resistance value.

and thank you ,i will try it

wnqiangjingjiu

4# tonnyw


thank you ! let me try first,