Learning Forum

[Riccardo Petrelli]

hello, I'm using an ekill script that kills the elements that exceed a certain strain and stress threshold in particular 0.035 and 1900 Mpa in a compression test for lattice structures in Ti6Al4-V, I want to reproduce the decreasing part of the obtained stress strain curve in the laboratory. Below I am attaching the graph of the force and you can see the decreasing curves, which are due to the killing of the elements. Can anyone tell me how to increase the stiffness of the killed elements? because by default it is set to 0, I would like to reduce it perhaps by half. Thank you. I am attaching my script for Ekill

 

SCRIPT EKILL :

total_steps=20

subst_int=3

subst_max=50

subst_min=3

neqit, 1500!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

 

    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

      ekill,kill_stress%k%

      ekill,kill_strain%k%

!      *status

    *enddo              

*enddo

finish

 

 This is the mesh:

[Erik Kostson]

Hi

 

Perhaps seacrh the internet for convergence and ansys, e.g.,:

 

Hope this helps

 

Erik