求教:command里如何实现将计算结果输出到一个文本文件？

ygl0325

需要用command进行explicit计算，并把一些节点上的结果直接写到一个文本文件，该用哪些命令呢？试过*node print， 但是这条命令在explicit里不支持。恳切向各位求教，谢谢了先！！！

ygl0325

哪位热心人肯帮帮俺啊?

ygl0325

怎么没有人回帖呢?

ygl0325

自己顶一下。

ygl0325

自己顶,好可怜.

cafe0519

焦急的心情可以理解,严肃的求学或治学态度容不得半点浮躁。
论坛要泡,过尤不及。看看论坛规则,在等待的同时多看看书本。

你说的两个命令(*file output和*node file )分别输出结果到ODB和FIL文件,建议你使用ABAQUS viewer来打开ODB文件后使用report输出你想要的结果数据的ASCII格式的文本文件。如果你对第3方处理软件熟悉,可以对FIL文件进行同样的操作也可以得到你的想要。

ygl0325

首先非常感谢cafe对小弟的关心与教导.你的回帖对于无助的我在精神上实属莫大的鼓励.
在ODB里用report进行输出处理是很好的方法.但是由于我需要对输出结果进行非常大量的运算和处理(一次次重复手工操作几乎不可能),现在正试图使用自编程序对输出节点位移的ASCII文本文件进行处理,这样可以省掉庞大的重复操作工作量.因此,就需要一个直接含有输出结果的文件.但是苦于不得其法.
那么,到底ABAQUS能否在其inp文件里实现上述功能呢?
另外,在我对运算产生的.sel文件进行convert,成为fil文件,进而转成.fin文件后,发现里面是一些奇怪的数据,好象并不是节点位移(附示例文件),不知为何?
希望cafe和其他高人再次不吝赐教.谢谢!

ygl0325

这里是转后的.fin文件.

cafe0519

你说说你的处理的复杂程度我看看,一般的需要对输出结果进行大量的运算和处理(一次次重复手工操作几乎不可能),可以在report中实现，非常复杂的可以通过python来搞定,至于试图使用自编程序对输出节点位移的ASCII文本文件进行处理,这也是可行的,我很久以前也是干过的,几个循环就搞定啦。 要想一劳永逸,建议些个脚本,这样只要模型没有变,就可以坐享,只是你得花些工夫,不是在论坛能吹出来的。

cathayee

把你的数据用el file or node file 输出到.fil文件里，然后用hksmain子程序读出你要的数据并编成作所需要的运算。建议参看abaqus user analysis manual -〉 output -〉file output format 章节。

SUBROUTINE HKSMAIN
C
C     Calculate the maximum von Mises stress and its location
C
      INCLUDE 'aba_param.inc'
      CHARACTER*80 FNAME
      DIMENSION ARRAY(513),JRRAY(NPRECD,513),LRUNIT(2,1)
      EQUIVALENCE (ARRAY(1),JRRAY(1,1))
C
C     File initialization
C
      FNAME='TEST'
      NRU=1
      LRUNIT(1,1)=8
      LRUNIT(2,1)=2
      LOUTF=0
      CALL INITPF(FNAME,NRU,LRUNIT,LOUTF)
      JUNIT=8
      CALL DBRNU(JUNIT)
C
C     Loop on all records in results file
C
      STRESS=0.
      DO 100 K1=1,99999
C
        CALL DBFILE(0,ARRAY,JRCD)
        IF(JRCD.NE.0)GO TO 110
        KEY=JRRAY(1,2)
C
        IF(KEY.EQ.1) THEN
C
C          Element header record:
C          extract element, sec pt, int pt numbers
C
           JEL=JRRAY(1,3)
           JPNT=JRRAY(1,4)
           JSPNT=JRRAY(1,5)
C
C     Stress invariant record for ABAQUS/Standard
        ELSE IF(KEY.EQ.12)THEN
C     Stress invariant record for ABAQUS/Explicit
        ELSE IF(KEY.EQ.75)THEN
C
C          Extract von Mises stress
C
           IF(ARRAY(3).GT.STRESS)THEN
              STRESS=ARRAY(3)
              KEL=JEL
              KPNT=JPNT
              KSPNT=JSPNT
           END IF
C
        END IF
C
100  CONTINUE
110  CONTINUE
C
      WRITE(6,120) KEL,KPNT,KSPNT,STRESS
120  FORMAT(5X,'ELEMENT',I5,5X,'POINT',I4,5X,'SECTION POINT',
     1 I4,5X,'STRESS',1PG12.3)
      STOP
      END

ygl0325

cafe0519 wrote:
你说说你的处理的复杂程度我看看,一般的需要对输出结果进行大量的运算和处理(一次次重复手工操作几乎不可能),可以在report中实现，非常复杂的可以通过python来搞定,至于试图使用自编程序对输出节点位移的ASCII文本文件进行处理,这也是可行的,我很久以前也是干过的,几个循环就搞定啦。 要想一劳永逸,建议些个脚本,这样只要模型没有变,就可以坐享,只是你得花些工夫,不是在论坛能吹出来的。

麻烦之处就是需要对一个完整模型的温度,外载,约束做些许改动,产生类似的几十上百个模型,而所要求输出的数据完全相同.所以这种情况下用report就不大可行.
即使几个循环能够搞定,可是我现在的困难是不知道能否产生所需要的ASCII文本文件,以及如何产生.

ygl0325

多谢您的建议!

ygl0325

大家接着讨论啊。

ygl0325

哪位有处理过类似问题经验的,希望能发表一下自己的看法.谢谢!

lxuejin

现在我的情况相同，explicit不支持print， 希望自己编写个程序读fil文件，同时又不用abaqus子程序。

给自己点建议，先用abaqus子程序读取fil文件