Riccardo Petrelli
Subscriber

 

The sudden increase in value is due to the elimination of the element that exceeds a certain stress and strain threshold (actually with the command estif 0.5 I am lowering its stiffness). I am attaching my script.

total_steps=20
subst_int=10
subst_max=18
subst_min=10
neqit, 120 !numero massimo di iterazioni di default è impostato su 15
 
*do,i,1,total_steps
  /gopr
  time_step=i/total_steps
  time,time_step
  nsubst,subst_int,subst_max,subst_min
  outres,all,all
  esel,all
  solve
  finish
 
  /post1
    set,last
 
    esel,all
    etable,stress,s,eqv
    esel,s,etab,stress,ARG1
!    *get,stress_ecount%i%,elem,0,count
    cm,kill_stress%i%,elem
    cmsel,s,kill_stress%i%,elem
!    *get,stress_ecount%i%,elem,0,count
 
    esel,all
    etable,strain,epto,eqv
    esel,s,etab,strain,ARG2
!    *get,strain_ecount%i%,elem,0,count
    cm,kill_strain%i%,elem
    cmsel,s,kill_strain%i%,elem
!    *get,strain_ecount%i%,elem,0,count
    
    allsel
    cmwrite,kill_list,cm
    parsav,all
  finish
  
  /solu
    antype,,restart
    parres
!    *status   
      /input,kill_list,cm
    *do,k,1,i
estif, 0.5 !moltiplica la matrice di rigidezza di 0.5
      ekill,kill_stress%k%
      ekill,kill_strain%k%
!      *status
    *enddo              
*enddo
finish
 
my value for stress eq is 1750 and strain 0,032
 
note that the number of steps is not the one written in the script, because I don’t remember it for this case.