abaqus调用fortran子程序读取txt文件报错forrtl: severe (24)

wzjdse · 发表于 2020-8-19 22:10:55

各位大佬，我想请教一下，我在用abaqus调用fortran子程序做分析，提示错误：forrtl: severe (24): end-of-file during read, unit 11, file D:\temp\xdata.txt。在网上查网友说是对同一个txt多次读取，后来我就改了一下程序，设置只读一次。如下为代码：但改完依然报同样的错误

[Fortran] syntaxhighlighter_viewsource syntaxhighlighter_copycode

if (stepTime==0) then
          
          !读入轴承力Fx数据
          open(11,file='D:\temp\xdata.txt',status='old')
          do ix=1,10001
              read(11,*)Fxm(ix,1),Fxm(ix,2)    !循环写入
          end do
      
          !读入轴承力Fy数据
          open(11,file='D:\temp\ydata.txt',status='old')
          do iy=1,10001
              read(11,*)Fym(iy,1),Fym(iy,2)    !循环写入
          end do
          
      end if

根据本论坛的一个贴又做了如下修改：

[Fortran] syntaxhighlighter_viewsource syntaxhighlighter_copycode

!只读入一次
      if (stepTime.eq.0) then
          
          !读入轴承力Fx数据
          do ix=1,10001
              open(11,file='D:\temp\xdata.txt',ioStat=krl,
     &        action="readwrite",status="old")
              if(krl==0) then
                  read(11,*)Fxm(ix,1),Fxm(ix,2)    !循环写入
                  close(11)
              end if
          end do
      
          !读入轴承力Fy数据
          do iy=1,10001
              open(11,file='D:\temp\ydata.txt',ioStat=krl,
     &        action="readwrite",status="old")
              if(krl==0) then
                  read(11,*)Fym(iy,1),Fym(iy,2)    !循环写入
                  close(11)
              end if
          end do
          
      end if

报错信息改变为：forrtl: severe (29): file not found, unit 11, file C:\Users\WANGZI~1\AppData\Local\Temp\WangZijia_Housing_15704\fort.11

请教各位大神，我原本以为读取txt很简单，没成想在这卡了好久了

wzjdse · 发表于 2020-8-19 22:12:31

我的思路其实很简答，仿真最开始，读入一个已经存在的txt，之后就不读了（之前方法比较笨，想着是每次abaqus调用子程序时都要读一遍这个txt，但是发现会报错forrtl: severe (24): end-of-file during read，所以正好优化一下只读一次）。

风平老涡 · 发表于 2020-8-20 01:39:05

在你原先的俩循环之间加一个close(11),

风平老涡 · 发表于 2020-8-20 01:44:58

还有就是你的txt文件是否有20002个数据。

wzjdse · 发表于 2020-8-20 08:42:54

风平老涡发表于 2020-8-20 01:44
还有就是你的txt文件是否有20002个数据。

我的txt就是10001行，两列，如下：
0.000000 19476.300000
0.000100 19475.066748
0.000200 19473.864231
0.000300 19472.702313
0.000400 19471.590858
0.000500 19470.539728
0.000600 19469.558788
0.000700 19468.657901
0.000800 19467.846931
0.000900 19467.135741
0.001000 19466.534195
0.001100 19466.049638
0.001200 19465.659818
0.001300 19465.343557
0.001400 19465.083727

我有两个txt ，分别是xdata和ydata，分别读他们两个并循环给数组赋值，我是这么个思路

wzjdse · 发表于 2020-8-20 08:46:13

风平老涡发表于 2020-8-20 01:44
还有就是你的txt文件是否有20002个数据。

如果是10001行，2列，那确实是20002个数据

chiangtp · 发表于 2020-8-20 12:35:10

0. 程式裡沒有其他地方對 'D:\temp\xdata.txt'檔案做READ動作

1. 'D:\temp\xdata.txt'檔案确实是10001行

2. 可能有若干行不小心的 "沒有2个数据"

3. 一個小程式驗證一下

[Fortran] syntaxhighlighter_viewsource syntaxhighlighter_copycode

PROGRAM test                                                                    
  IMPLICIT NONE                                                                 
                                                                                
  INTEGER :: i, ier                                                             
  REAL    :: a, b                                                               
                                                                                
  !-----------------                                                            
                                                                                
  OPEN(UNIT=11, FILE='D:\temp\xdata.txt', status='old', IOSTAT=ier)             
  IF( ier /= 0 ) STOP 'Open Error'                                              
                                                                                
  OPEN(UNIT=22, FILE='check.txt')                                               
                                                                                
  !-------                                                                      
                                                                                
  DO i = 0, 10000                                                               
    READ(11, *, IOSTAT=ier) a, b                                                
                                                                                
    IF( ier /= 0 ) THEN                                                         
      WRITE(*,*) i                                                              
      STOP 'Reading Error'                                                      
    END IF                                                                      
                                                                                
   !IF( NINT(a*10000.0) /= i ) THEN                                             
   !  WRITE(*,*) i, a, b                                                        
   !  STOP 'Data Error?'                                                        
   !END IF                                                                      
                                                                                
    WRITE(22,*) i, a, b                                                         
  END DO                                                                        
                                                                                
  !-------                                                                      
                                                                                
  CLOSE(UNIT=11)                                                                
  CLOSE(UNIT=22)                                                                
                                                                                
END PROGRAM test

wzjdse · 发表于 2020-8-20 14:54:04

chiangtp 发表于 2020-8-20 12:35
0. 程式裡沒有其他地方對 'D:\temp\xdata.txt'檔案做READ動作

1. 'D:\temp\xdata.txt'檔案确实是10001 ...

您好，我只读一次的问题已经，解决，我换了一个判据。但是我想把数组Fxm和Fym只读一次就保存下来，下次ABAQUS再调用子程序的时候就能用，所以用了save。但是一加save后，又出现了上述错误：forrtl: severe (24): end-of-file during read, unit 55, file D:\temp\xdata.txt. 下面是我的代码：

[Fortran] syntaxhighlighter_viewsource syntaxhighlighter_copycode

dimension Fxm(10001,2)
      dimension Fym(10001,2)
      double precision Fxm,Fym,Fx,Fy,F

      data iread /1/
      save iread
      save Fxm
      save Fym
      
      !只读入一次
      if (iread.eq.1) then
          
          iread=2
          !读入轴承力Fx数据
          open(55,file='D:\temp\xdata.txt',status='old')
          do ix=1,10001
              read(55,*)Fxm(ix,1),Fxm(ix,2)    !循环写入
          end do
          close(55)
      
          !读入轴承力Fy数据
          open(22,file='D:\temp\ydata.txt',status='old')
          do iy=1,10001
              read(22,*)Fym(iy,1),Fym(iy,2)    !循环写入
          end do
          close(22)
          
      end if

想请教您这是为何

chiangtp · 发表于 2020-8-20 15:49:01

1. 用abaqus调用fortran子程序, "abaqus" 是 Fortran library?
2. 上述程序是子程序, 如果是被某一Fortran程序調用, SAVE: 應該只會讀一次, 数组Fxm和Fym會保存下来
3. 試試

[Fortran] syntaxhighlighter_viewsource syntaxhighlighter_copycode

      if (iread.eq.1) then                                           
          open(55,file='D:\temp\xdata.txt',status='old')             
          do ix=1,10001                                              
              read(55,*) Fxm(ix,1), Fxm(ix,2)                        
              IF( NINT(Fxm(ix,1)*10000.0) /= ix ) THEN               
                WRITE(*,*) iread, ix, Fxm(ix,1), Fxm(ix,2)           
                STOP 'Data Error?'                                   
              END IF                                                 
          end do                                                     
          close(55)                                                  
          iread = 2 !---> move to here

wzjdse · 发表于 2020-8-20 16:01:42

chiangtp 发表于 2020-8-20 15:49
1. 用abaqus调用fortran子程序, "abaqus" 是 Fortran library?
2. 上述程序是子程序, 如果是被某一Fortran ...

如果save的数组太大，会不会出现我这个错误，像我的这个数组10001*2，这么大的

[文件读写] abaqus调用fortran子程序读取txt文件报错forrtl: severe (24)

规矩勋章

贡献勋章

管理勋章

帅哥勋章

元老勋章

星光勋章