Fortran Coder

查看: 1530|回复: 2
打印 上一主题 下一主题

[文件读写] abaqus计算,文件读取,fortran语法save

[复制链接]

7

帖子

3

主题

0

精华

入门

F 币
39 元
贡献
17 点
跳转到指定楼层
楼主
发表于 2023-7-14 12:55:49 | 只看该作者 回帖奖励 |倒序浏览 |阅读模式
以下 是我abaqus dflux子程序的一段,想要实现对文件的一次读写,使用了data ii/1/, 并进行保存,然后再最后一句使用print打ii输出结果,在log文件里发现输出的i全是以下结果:
Run standard.exe
         241  1076423557
         241  1076423557
         241  1076423557
         241  1076423557
         241  1076423557
         241  1076423557
         ....
         241 1076423557
而且我的分析步只进行了一次增量步,但是输出的ii有三十万行。这个abaqus的计算步骤是怎么进行的啊。   
子程序语言:   
      data ii /1/
      save ii
      save TT
      save WENDU

      if(ii.eq.1) then
      !读入温度数据
      ii =  2
      call ReadDate(TT,WENDU,ntime)
      end if
      ! 对应时间步长
      do 100 j = 1,ntime
      if (TIME(2).ge.TT(j).and.TIME(2).lt.TT(j+1)) then
          SINK = WENDU(j)+(TIME(2)-TT(j))/(TT(j+1)-TT(j))*(WENDU(j+1)-
     1WENDU(j))
          PRINT*,ntime,ii


分享到:  微信微信
收藏收藏 点赞点赞 点踩点踩

7

帖子

3

主题

0

精华

入门

F 币
39 元
贡献
17 点
沙发
 楼主| 发表于 2023-7-14 12:58:08 | 只看该作者
SUBROUTINE ReadDate(TT,WENDU,ntime)
C
      INCLUDE 'ABA_PARAM.INC'
C
      DIMENSION TT(ntime),WENDU(ntime)
      
      open(12,file='D:\abaqusworks\abaqus\temperature.txt')
      do 10 i=1,ntime
          read(12,*) TT(i),WENDU(i)
10        continue
          close(12)
      RETURN
      END

7

帖子

3

主题

0

精华

入门

F 币
39 元
贡献
17 点
板凳
 楼主| 发表于 2023-7-14 12:58:45 | 只看该作者
parameter(ntime=241)
您需要登录后才可以回帖 登录 | 极速注册

本版积分规则

捐赠本站|Archiver|关于我们 About Us|小黑屋|Fcode ( 京ICP备18005632-2号 )

GMT+8, 2024-11-23 17:13

Powered by Tencent X3.4

© 2013-2024 Tencent

快速回复 返回顶部 返回列表