I am writing Abaqus USDFLD to change my Young's Modulus after strain crosses the value of 0.035 and start deleting elements when strain value crosses 0.05 in a heterogeneous model (composed of two different sections).

I am using linear degradation; (1-D)*E where D goes from 0 to 1 as strain goes from 0.035 to 0.05.

See the image below, and I want that degradation to follow the linear curve.

So, my problem is how do I infer the young's modulus (E) and which section I am at when I use the GETVRM subroutine to get values at the integration point? I have checked the Abaqus documentation but couldn't find the method to get E at the integration points. Is there any alternative way to accomplish this?

Is there any way to use different subroutines for different material properties, as this may solve my integration point?

EDIT:

After reading some comments, I have written this code but getting compilation errors; please help:

I have defined field(1) in Young's Modulus table (as suggested in the comments), see images below, and my element deletion state variable is 2.

​​​​​​​

      subroutine usdfld(field,statev,pnewdt,direct,t,celent,
     1 time,dtime,cmname,orname,nfield,nstatv,noel,npt,layer,
     2 kspt,kstep,kinc,ndi,nshr,coord,jmac,jmatyp,matlayo,laccfla)
c
      include 'aba_param.inc'
c
      character*80 cmname,orname
      character*3  flgray(15)
      dimension field(nfield),statev(nstatv),direct(3,3),
     1 t(3,3),time(2)
      dimension array(15),jarray(15),jmac(*),jmatyp(*),coord(*)

      real :: start_strain, end_strain, D, strain_max
      start_strain = 0.035
      end_strain = 0.05
    
    
        
      CALL GETVRM('EP', ARRAY, JARRAY, FLGRAY, JRCD,
     1 JMAC, JMATYP, MATLAYO, LACCFLA)
       strain_max = max(abs(array(1)), abs(array(3)))
    
       if (strain_max .gt. 0.035 .and. strain_max .lt. 0.05) then
          D = (strain_max - start_strain)/(end_strain - start_strain)
          field(1) = (1 - D)*field(1)
          statev(1) = field(1)
          statev(2) = 1
       end if
            
       if (strain_max .ge. 0.05) then
          statev(2) = 0
       else
          statev(2) = 1
       end if
    
        
      return
      end