abaqus调用fortran子程序读写txt时出现steptime重复，为什么？ - Fortran Coder

[Fortran] 纯文本查看 复制代码

C

C User subroutine VDLOAD

            subroutine vdload (

C Read only -

     *     nblock, ndim, stepTime, totalTime, amplitude,

     *     curCoords, velocity, dircos, jltyp, sname,

C Write only -

     *     value )

C

      include 'vaba_param.inc'

      parameter ( const = 1.d5 ) 

C

      dimension curCoords(nblock,ndim), velocity(nblock,ndim),

     *     dircos(nblock,ndim,ndim), value(nblock)

*

      character*80 sname

*

      integer ::i

c      character(len=20)::filename="duxie.txt"

      REAL::ta

      real ::ta1

      real::tz



      open(93,FILE="D:\abaqus2016\SIMULIA\duxie.txt",status="old")      

      do

      read(93,*,iostat = ios) ta,ta1

      if ( ios == -1 ) exit 

       end do

       backspace(93)

       backspace(93)



              ta=ta+2*stepTime

              ta1=ta1+5*stepTime

              tz=0.5*ta

              tx=0.5*ta1

      write(93,*) ta,ta1

      close(93) 

      open(89,FILE="D:\abaqus2016\SIMULIA\duxiet.txt",position="append")

      write(89,*) ta,ta1,stepTime, tz

      close(89) 

      

      if( jltyp .eq. 2 ) then

        do k = 1, nblock        

        value(k) =tz

        end do

      end if

        if( jltyp .eq. 1 ) then

        do k = 1, nblock        

        value(k) =tx

        end do

        end if

    

      return

      end

[Fortran] 纯文本查看 复制代码

            subroutine vdload (

C Read only -

     *     nblock, ndim, stepTime, totalTime, amplitude,

     *     curCoords, velocity, dircos, jltyp, sname,

C Write only -

     *     value )

C

      include 'vaba_param.inc'

      parameter ( const = 1.d5 ) 

C

      dimension curCoords(nblock,ndim), velocity(nblock,ndim),

     *     dircos(nblock,ndim,ndim), value(nblock)

*

      character*80 sname

*

      REAL::ta

      open(93,FILE="D:\abaqus2016\SIMULIA\duxie.txt",status="old")      

       read(93,*) ta

       rewind(93)

       

       ta=ta+4*stepTime

       

      write(93,*) ta

      close(93) 

      open(89,FILE="D:\abaqus2016\SIMULIA\duxiet.txt",position="append")

      write(89,*) ta,stepTime

      close(89) 

      

      if( jltyp .eq. 2 ) then

        do k = 1, nblock        

        value(k) =-ta/4    ！体积是4， value(k)表示体力，体力X体积等于施加的反力（ value(k)是给到abaqus里的体力，abaqus能识别，cae里面后处理得到的反力等于 value(k)X4）

        end do

      end if    

      return

      end