Learning Forum

[Riccardo Petrelli]

 

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.