关于回弹警告和错误

byb

最近在学习回弹算法。从论坛上下载到一个实例，经过仔细查看，发现步骤应该没有问题。却出现不少警告和错误。

请问做过回弹的朋友，是不是显式过渡到隐式计算时候，会出现如同：
1，必须设置关闭单元形状检测，可是会出现警告： Deactivation of element shape checking is not recommended.   
2， 只选择板料节点，忽视刚体时候，提示警告： Out of the 535 defined elements, only 50 elements are selected.   
3，还有如同：

*** WARNING ***                        
The program chosen initial timestep/load-factor is arbitrary.  It is  necessary for the user to supply a suitable initial  timestep/load-factor through the NSUB or DELTIM command for  convergence and overall efficiency.      

*** WARNING ***                     
Element order list length 535 does not match number of selected   element 50     
  New order vector being created.   

*** WARNING ***                     
There are 1 small equation solver pivot terms.   

*** WARNING ***   SUPPRESSED MESSAGE   
The calculated reference MOMENT CONVERGENCE VALUE = 1.395454264E-04 is   less than 1E-2 a threshold zero, A value of 1E-2 or specified MINREF    is used.  Check results carefully.   

这些警告是不是无法避免的？

4，步骤没有错，竟然不收敛，是什么原因，难道只有通过设置收敛准则吗？
*** ERROR ***     SUPPRESSED MESSAGE   
Solution not converged at time 1.E-02 (load step 1 substep 12).     Run terminated.   

*** WARNING ***   SUPPRESSED MESSAGE   
The unconverged solution (identified as time 1.E-02 substep 999999) is   output for analysis debug purposes.  Results should not be used for   any other purpose.      

附：
得出上面问题的名流，请大侠们帮忙看下。

         
! stamp.inp is used in the Chapter 10 exercise
/title,Explicit-to-Implicit Sequential Solution (Springback Example)
/filnam,stamp
/prep7
/view,,1,2,3
/plopts,info,1
/pnum,mat,1
/num,1
et,1,SHELL163
keyopt,1,1,10   ! 10th element formulation recommended for springback
r,1,,3,0.05
edshell,,,1     ! membrane straining causes thickness change (default)
et,2,SOLID164
r,2
et,3,SOLID164
r,3
mp,ex,1,30e6
mp,nuxy,1,0.29
mp,dens,1,2.28e-3
tb,biso,1,1
tbdata,1,30e3,100e3
edmp,hgls,1,5
mp,ex,2,30e6
mp,nuxy,2,0.29
mp,dens,2,2.28e-3
edmp,rigid,2,6,7
mp,ex,3,30e6
mp,nuxy,3,0.29
mp,dens,3,2.28e-3
edmp,rigid,3,7,7
wpoffs,,-0.069
wprota,,-90
rectng,-0.99,-0.01,0.01,0.49
type,1
mat,1
real,1
lesize,1,,,10,1
lesize,2,,,5,1
lesize,3,,,10,1
lesize,4,,,5,1
amesh,1
wpcsys,-1,0
wpstyl
esize,0.1

k, 11, 0.00, 0.00, 0.00
k, 12, 0.00, 0.50, 0.00
k, 13,-1.00, 0.50, 0.00
k, 14,-1.00, 0.10, 0.00
k, 15, 0.50, 0.05, 0.00
k, 16, 0.00,-0.50, 0.00
k, 17, 0.00,-0.20, 0.00
k, 18,-1.00,-0.10, 0.00
k, 19,-1.00,-0.50, 0.00
k, 20, 0.50,-0.15, 0.00
l,11,12
l,12,13
l,13,14
larc,14,11,15,6
al,5,6,7,8
mat,2
amesh,2
l,16,17
larc,17,18,20,6
l,18,19
l,19,16
al,9,10,11,12
mat,3
amesh,3
esize,,5
mat,2
type,2
voffst,2,-0.5
aclear,2
mat,3
type,3
voffst,3,-0.5
aclear,3
eplot
edwrite,pcreate
edwrite,plist
edcgen,assc
fini
/solu

esel,s,mat,,1
cm,eblank,elem
nsle
cm,nblank,node
nsel,r,loc,x,-0.995,-0.985
nsel,u,loc,z,-0.485,-0.015
d,all,ux,0.0
nsel,r,loc,z,-0.495,-0.485
d,all,uz,0.0
esel,s,mat,,2
cm,epunch,elem
nsle
cm,npunch,node
esel,s,mat,,3
cm,edie,elem
nsle
cm,ndie,node
nsel,all
esel,all
! Punch displacement function = (amplitude*cos(omega*time))-amplitude
*dim,time1,,101,1,1           ! dimension array for time
*dim,disp1,,101,1,1           ! dimension array for punch displacement
curtime=0.0                   ! initialize current time to zero
tottime=0.0101                ! total load curve time, seconds
deltime=tottime/100.0         ! time increment, seconds
amp1=0.0825                   ! displacement amplitude, inches
omega1=180.0/tottime          ! product of this constant and time = wt
*afun,deg                     ! work in degrees
*do,i,1,101,1
   time1(i,1,1)=curtime
   disp1(i,1,1)=(amp1*cos(omega1*curtime))-amp1  ! maximum uy = -0.165"
   curtime=curtime+deltime
*enddo
edload,add,rbuy,,2,time1,disp1,0      ! punch = part number 2
time,0.01
edrst,50
edhtime,100
edhist,nblank
edhist,eblank
edenergy,1,1,1,1
edout,glstat
save
eplot

SAVE
/STAT,SOLU  
SOLVE   
ESEL,S,MAT,,1   
NSLE,S  
FINISH  
/POST1  
SET,LAST
!*  
/DSCALE,1,1.0   
/REPLOT
!*  
!*  
PLNS,S,EQV  
ANDATA,0.5, ,0,1,51,5,0,0   
!*  
FINISH  
/POST26
NUMVAR,30,  
FILE,stamp,his,
!*  
!*  
!*  
!*  
NSOL,2,2,U,Y,node2uy
!*  
STORE,MERGE, ,  
PLVAR,2, , , , , , , , , ,  
ALLSEL,ALL  
SAVE
!*  
!*  
FINISH  
/FILNAM,spring  
!*  
/PREP7  
ETCHG,ETI   
!*  
TBDELE,ALL,ALL, ,   
SHPP,OFF
!*  
UPGEOM,1,1,51,stamp,rst,
!*  
EPLOT   
ESEL,S,MAT,,1   
NSLE,S  
EPLOT   
FINISH  
/SOLU   

FLST,2,2,1,ORDE,2   
FITEM,2,1   
FITEM,2,17  
!*  
/GO
D,P51X, ,0.0, , , ,UY, , , , ,  
RIMPORT,DYNA,STRESS,ELEM,1,51,stamp,rst,
!*  
NLGEOM,1
NROPT,AUTO, ,   
LUMPM,0
EQSLV,SPAR, ,0,
PREC,0  
PIVCHECK,1  
SSTIF   
PSTRES  
TOFFST,0,   
!*  
SAVE
/STAT,SOLU  
SOLVE   
FINISH  
/POST1  
SET,LAST
!*  
/DSCALE,1,10.0
/REPLOT
!*  
*SET,_ZF,1  
AVPRIN,0,0,
!*  
/EFACE,_ZF  
PLNSOL,U,SUM,2,1
FINISH  
!*  
! /EXIT,NOSAV   


最后谢谢大家！

jubinskd

正在学习