subroutine uhyper(bi1,bi2,aj,u,ui1,ui2,ui3,temp,noel,
     .                 cmname,incmpflag,numstatev,statev,
     .                 numfieldv,fieldv,fieldvinc,numprops,props)

      include 'aba_param.inc'
c
      character*8 cmname
      dimension ui1(3),ui2(6),ui3(6),statev(*),fieldv(*),
     .          fieldvinc(*),props(*)
c
      if (noel.eq.1) then
      end if
c
      c10 = props(7)
      c01 = props(8)
c
      if(fieldv(1) .le. props(1)) then
        d1 = props(2)
      else if(fieldv(1) .ge. props(5)) then
        d1 = props(6)
      else if(fieldv(1).gt.props(1) .and.
     .        fieldv(1).le.props(3)) then
        d1 = props(2)+(props(4)-props(2))/
     .  (props(3)-props(1))*(fieldv(1)-props(1))
      else
        d1 = props(4)+(props(6)-props(4))/
     .  (props(5)-props(3))*(fieldv(1)-props(3))
      end if
c
      u = c10*(bi1-3.)+c01*(bi2-3.)+((aj-1.)**2)/d1
      ui1(1) = c10
      ui1(2) = c01
      ui1(3) = 2./d1*(aj-1.)
      ui2(1) = 0.
      ui2(2) = 0.
      ui2(3) = 2./d1
      ui2(4) = 0.
      ui2(5) = 0.
      ui2(6) = 0.
      ui3(1) = 0.
      ui3(2) = 0.
      ui3(3) = 0.
      ui3(4) = 0.
      ui3(5) = 0.
      ui3(6) = 0.

      return
      end