大家好，目前試著在用VUSDFLD實(shí)現(xiàn)單元?jiǎng)h除。刪除準(zhǔn)則很簡(jiǎn)單，即x方向應(yīng)變大于0.001時(shí)，單元失效刪除。根據(jù)論壇上已有的例子，試著編寫了一段VUSDFLD子程序，但是運(yùn)行后，發(fā)現(xiàn)結(jié)果與不加入子程序并無任何差別，子程序似乎沒有起作用，請(qǐng)教大家解答。下面附上我的子程序，由于是第一次編寫，估計(jì)有錯(cuò)誤，請(qǐng)大家不吝賜教，感謝
c
c User subroutine VUSDFLD for user-defined fields
c
      subroutine vusdfld(
c Read only -
     *   nblock, nstatev, nfieldv, nprops, ndir, nshr, 
     *   jElemUid, kIntPt, kLayer, kSecPt, 
     *   stepTime, totalTime, dt, cmname, 
     *   coordMp, direct, T, charLength, props, 
     *   stateOld, 
c Write only -
     *   stateNew, field )
c
      include 'vaba_param.inc'
c
      dimension props(nprops),
     *          jElemUid(nblock), coordMp(nblock, *), 
     *          direct(nblock, 3, 3), T(nblock,3,3), 
     *          charLength(nblock),
     *          stateOld(nblock, nstatev), 
     *          stateNew(nblock, nstatev),
     *          field(nblock, nfieldv)
      character*80 cmname

c
      character*3 cData(maxblk*6)
      dimension jData(maxblk*6)
      dimension stress(maxblk*6)

   
c Get srain from previous increment
      jStatus = 1
      call vgetvrm( 'LE', strain, jData, cData, jStatus)
c
      call LE11strain(nblock, nstatev,
     *     nfieldv, ndir, nshr, 
     *     strain, 
     *     stateOld, 
     *     stateNew, field)
c
      return 
      end
c
      subroutine  LE11strain(nblock, nstatev,
     *     nfieldv, ndir, nshr, 
     *     strain, 
     *     stateOld, 
     *     stateNew, field)
c
      include 'vaba_param.inc'
c
      dimension strain(nblock,ndir+nshr),
     *     stateOld(nblock,nstatev), 
     *     stateNew(nblock,nstatev),
     *     field(nblock,nfieldv)
c
      parameter ( zero = 0.d0 )     

c Read properties from .inp
      dc = 0.001      
c Get peeq from previous increment
      jStatus = 1
      call vgetvrm( 'LE', strain, jData, cData, jStatus )      
      do k = 1, nblock
         if(strain(k,1).gt.dc) then 
c Element Deletion
           stateNew(k,1) = zero
       
             
         end if
      end do
c
      return 
      end