Learning Forum

[Riccardo Petrelli]

hello everyone, it's the umpteenth post I'm making on this topic because I can't find a solution, I should simulate this compression test performed in the laboratory (on three specimens, yellow, gray and orange curve, the blue curve is  to reproduce numerically). To be able to make the curve descend, I thought of using a script with the EKILL command which brings the stiffness to 0 but then using the ESTIF command I bring it back to 90% (I'm trying). The script works well but as the central rods are eliminated the program is no longer able to converge. Could you show me a new way in as much detail as possible, or how to improve mine. I am working on my thesis work. The script mine will

write in first comment

[Riccardo Petrelli]

 

I tried to improve the mesh and increase substeps and steps, but it doesn't solve much. Thank you all. SCRIPT:

 

total_steps=20

subst_int=10

subst_max=18

subst_min=10

neqit, 500 !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,s,live

!    *get,live_count%i%,elem,0,count

!    esel,inve,live

!    *get,dead_count%i%,elem,0,count

             

    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

 

 

    

    allsel

    cmwrite,kill_list,cm

    parsav,all

  finish

  

  /solu

    antype,,restart

    parres

!    *status   

      /input,kill_list,cm

    *do,k,1,i

      ekill,kill_stress%k%

      estif, 0.9

    

!      *status

    *enddo              

*enddo

finish

 

 

[Riccardo Petrelli]

I'm using static structural