请教大神，怎么样将应变输出，为什么我输出的都是0

淋雨后的晴朗

1. SUBROUTINE UMAT(STRESS,STATEV,DDSDDE,SSE,SPD,SCD,RPL,DDSDDT,
   
2.      1 DRPLDE,DRPLDT,STRAN,DSTRAN,TIME,DTIME,TEMP,DTEMP,PREDEF,DPRED,
   
3.      2 CMNAME,NDI,NSHR,NTENS,NSTATV,PROPS,NPROPS,COORDS,DROT,
   
4.      3 PNEWDT,CELENT,DFGRDO,DFGRDl,NOEL,NPT,LAYER,KSPT,KSTEP,KINC)
   
5.       
   
6.         include 'aba_param.inc'
   
7.          CHARACTER*8 CMNAME
   
8.           DIMENSION STRESS(NTENS),STATEV(NSTATV),DDSDDE(NTENS,NTENS),
   
9.      1 DDSDDT(NTENS),DRPLDE(NTENS),STRAN(NTENS),DSTRAN(NTENS),
   
10.      2 TIME(2),PREDEF(1),DPRED(1),PROPS(NPROPS),COORDS(3),DROT(3,3),
    
11.      3 DFGRD0(3,3),DFGRDl(3,3)
    
12. C UMAT FOR ISOTROPIC ELASTICITY
    
13. C CANNOT BE USED FOR PLANE STRESS
    
14. C -  ------ -
    
15. C PROPS(1) - E
    
16. C PROPS(2) - NU
    
17. CC
    
18.                
    
19.           IF (NDI.NE.3) THEN
    
20.              WRITE (6,*) 'THIS UMAT MAY ONLY BE USED FOR ELEMENTS
    
21.      1 WITH THREE DIRECT STRESS COMPONENTS'
    
22.             CALL XIT
    
23.           ENDIF
    
24. C
    
25. C ELASTIC PROPERTIES
    
26.         EMOD=PROPS(1)                     
    
27.         ENU=PROPS(2)                     
    
28.         EBULK3=EMOD/(1-2*ENU)
    
29.         EG2=EMOD/(1+ENU)
    
30.         EG=EG2/2
    
31.         EG3=3*EG
    
32.         ELAM=(EBULK3-EG2)/3
    
33. C
    
34. C ELASTIC STIFFNESS
    
35. C
    
36.             DO K1=1, NDI
    
37.                DO K2=1, NDI
    
38.             DDSDDE(K2, K1)=ELAM
    
39.                            END DO
    
40.             DDSDDE(K1, K1)=EG2+ELAM
    
41.             END DO
    
42.             DO K1=NDI+1, NTENS
    
43.              DDSDDE(K1 ,K1)=EG
    
44.             END DO
    
45. C
    
46. C CALCULATE STRESS
    
47. C
    
48.              DO K1=1, NTENS
    
49.              DO K2=1, NTENS
    
50.             STRESS(K2)=STRESS(K2)+DDSDDE(K2, K1)*DSTRAN(K1)
    
51.             END DO
    
52.           END DO
    
53. 
    
54.                
    
55. C                WRITE (6,*),'NSTATV',NSTATV
    
56. C                WRITE (6,*),'STRESS',STRESS
    
57. C                WRITE (6,*),'DDSDDE',DDSDDE
    
58.                 WRITE (6,*),'DSTRAN',DSTRAN
    
59.                 WRITE (6,*),'STRAN' ,STRAN
    
60. C                WRITE (6,*),'STATEV',STATEV
    
61.           RETURN
    
62.         END

复制代码

空中飞客

请问6代表什么，是某个打开的文本吗，输出内容会显示在CLI窗口吗

淋雨后的晴朗

空中飞客 发表于 2017-9-15 22:25
请问6代表什么，是某个打开的文本吗，输出内容会显示在CLI窗口吗

我查了在fortran语言中6代表的是屏幕，我运行完了可以在.dat文件中查看到我输出的应变值

cassiliyas

你的NDI， NTENS 怎么没有看见在哪里赋值了。

caezyf1990

我在调试子程序的时候一般都是open一个文件，然后将结果write到此文件中，然后看程序是否运行正确，还没用过直接print到屏幕上

caezyf1990

还有就是应变和应力在ABAQUS中输出的都是分量，阁下也没有指出要输出哪个分量，我懂的也不多，希望说的是对的吧