请教UMAT出现的问题

forcamingseu

IF(KSTEP.EQ.2)THEN
       DO K1=1,NTENS
        DO K2=1,NTENS
          DDSDDE(K1,K2)=ZERO
        END DO
       END DO
       DO K1=1,NDI
        DO K2=1,NDI
          DDSDDE(K2,K1)=ELAM
        END DO
        DDSDDE(K1,K1)=EG2+ELAM
       END DO
       DO K1=NDI+1,NTENS
         DDSDDE(K1,K1)=EG
       END DO
       SMISES1=(STRESS(1)-STRESS(2))**2+(STRESS(2)-STRESS(3))**2
     1         +(STRESS(3)-STRESS(1))**2
       DO K1=NDI+1,NTENS
         SMISES1=SMISES1+SIX*STRESS(K1)**2
       END DO
       SMISES1=SQRT(SMISES1/TWO)
  EQVSER1=SMISES1*EXP(-1.0*PROPS(6)*TIME(1))/PROPS(5)
       EQVEER1=SMISES1*EXP(-1.0*PROPS(3)/PROPS(4)*TIME(1))/PROPS(4)
       EQCRER1=EQVSER1+EQVEER1
  DO K1=1,NTENS
   STREST(K1)=STRESS(K1)
  END DO
  DO K1=1,NTENS
    DECRT(K1)=0.0
  END DO
C
       LOOP1O KEWTON=1,NEWTON
C    CALCULATE EQUIVALENT VON MISES STRESS
C
        SMISES2=(STREST(1)-STREST(2))**2+(STREST(2)-STREST(3))**2
     1          +(STREST(3)-STREST(1))**2
        DO K1=NDI+1,NTENS
          SMISES2=SMISES2+SIX*STREST(K1)**2
        END DO
        SMISES2=SQRT(SMISES2/TWO)
        SMISES3=0.2*SMISES1+0.8*SMISES2
C
C    CALCULATE EQUIVALENT VISCO,VISCOELASTIC AND CREEP STRAIN
C    INCREMENTS USING MODIFIED BURGERS MODEL
C
        EQVSER2=SMISES2*EXP(-1.0*PROPS(6)*(TIME(1)+DTIME))/PROPS(5)
        EQVEER2=SMISES2*EXP(-1.0*PROPS(3)/PROPS(4)*(TIME(1)+DTIME))
     1          /PROPS(4)
        EQCRER2=EQVSER2+EQVEER2
        EQCRER3=0.2*EQCRER1+0.8*EQCRER2
        TERM1=(THREE/TWO*EQCRER3)/SMISES3
C
C    CALCULATE CREEP STRAIN INCREMENTS
C
        DO K1=1,NTENS
          STREST2(K1)=0.2*STREST(K1)+0.8*STRESS(K1)
        END DO
C    PRESSURE
        PRESS=0.0
        DO K1=1,NDI
          PRESS=PRESS+STREST2(K1)/THREE
        END DO
C    DEVIATORIC STRESS
        DO K1=1,NDI
          DVSTRESS(K1)=STREST2(K1)-PRESS
        END DO
        DO K2=NDI+1,NTENS
          DVSTRESS(K2)=STREST2(K2)
        END DO
C    CREEP STRAIN INCREMENTS
        DO K1=1,NTENS
          DECRE(K1)=TERM1*DTIME*DVSTRESS(K1)
          DEELA(K1)=DSTRAN(K1)-DECRE(K1)
        END DO
        DEQCRE=EQCRER3*DTIME
        DEQVSE=(0.2*EQVSER1+0.8*EQVSER2)*DTIME
        DEQVEE=(0.2*EQVEER1+0.8*EQVEER2)*DTIME
        IF(NOEL.EQ.3.AND.NPT.EQ.1)THEN
          WRITE(*,*) 'STREST,DECRE,DSTRAN,DEELA',KEWTON
          DO K1=1,NTENS
            WRITE(*,*) STREST(K1),DECRE(K1),DSTRAN(K1),DEELA(K1)
          END DO
        END IF
C
C       CONVERGE JUDGEMENT
C
        TERM2=0.0
        TERM3=0.0
        DO K1=1,NTENS
          TERM2=TERM2+(DECRE(K1)-DECRT(K1))**2
          TERM3=TERM3+DECRE(K1)**2
        END DO
        TERM4=SQRT(TERM2)
        TERM5=SQRT(TERM3)
        TERM6=TERM4/TERM5
        IF(TERM6.LE.TOLER)THEN
          EXIT LOOP1
        END IF      
C
C    CALCULATE UPDATED STRESS
C
        DO K1=1,NTENS
          DO K2=1,NTENS
            STREST(K2)=STRESS(K2)+DDSDDE(K2,K1)*DEELA(K1)
          END DO
        END DO
C
   DO K1=1,NTENS
     DECRT(K1)=DECRE(K1)
   END DO
       END DO LOOP1
C
       DO K1=1,NTENS
    STRESS(K1)=STREST(K1)
  END DO
C
C    UPDATE CREEP,ELASTIC STRAINS AND EQUIVALENT STRAINS
C
       DO K1=1,NTENS
         ECREE(K1)=ECREE(K1)+DECRE(K1)
         EELAS(K1)=EELAS(K1)+DEELA(K1)
       END DO
       EQCRE=EQCRE+DEQCRE
       EQVSE=EQVSE+DEQVSE
       EQVEE=EQVEE+DEQVEE
C
      END IF

各位大侠，这是我的UMAT子程序的一部分，其中  STREST(K2)=STRESS(K2)+DDSDDE(K2,K1)*DEELA(K1)是用来更新应力的，通过输出STREST发现，直接应力部分不更新，而剪切应力部分更新了，不知道是什么原因，那位大侠知道的话，还请赐教！
下面是输出的一部分结果，最后一列是总应变增量DSTRAN
STRESS,STRAIN           2           0  1.000000000000000E-009
  -312255.914057127       2.271389122026695E-004  0.000000000000000E+000
  -729161.657581303      -2.219916538917496E-003  0.000000000000000E+000
  -312260.810600813       2.271101716201623E-004  0.000000000000000E+000
  -752.508288874525      -8.833792956353115E-006  0.000000000000000E+000
STREST,DECRE           1
  -312255.914057127       5.085472553290782E-013  0.000000000000000E+000
  -729161.657581303      -1.017076592258290E-012  0.000000000000000E+000
  -312260.810600813       5.085293369292123E-013  0.000000000000000E+000
  -752.508288874525      -2.753726972695362E-015  0.000000000000000E+000
STREST,DECRE           2
  -312255.914057127       5.085472553290782E-013  0.000000000000000E+000
  -729161.657581303      -1.017076592258290E-012  0.000000000000000E+000
  -312260.810600813       5.085293369292123E-013  0.000000000000000E+000
  -752.508288639948      -2.753726972523680E-015  0.000000000000000E+000
STRESS,STRAIN           2           1  1.000000000000000E-003
  -312255.914057127       2.271389122026695E-004  0.000000000000000E+000
  -729161.657581303      -2.219916538917496E-003  0.000000000000000E+000
  -312260.810600813       2.271101716201623E-004  0.000000000000000E+000
  -752.508288874525      -8.833792956353115E-006  0.000000000000000E+000
STREST,DECRE           1
  -312255.914057127       5.085040822022346E-007  0.000000000000000E+000
  -729161.657581303      -1.016990247525786E-006  0.000000000000000E+000
  -312260.810600813       5.084861653235515E-007  0.000000000000000E+000
  -752.508288874525      -2.753493195002849E-009  0.000000000000000E+000
STREST,DECRE           2
  -312255.914057127       5.085040822022610E-007  0.000000000000000E+000
  -729161.657581303      -1.016990247525839E-006  0.000000000000000E+000
  -312260.810600813       5.084861653235779E-007  0.000000000000000E+000
  -752.273732046802      -2.753321542241547E-009  0.000000000000000E+000
STRESS,STRAIN           2           1  1.000000000000000E-003
  -312255.914057127       2.271389122026695E-004  2.403003333928053E-008
  -729161.657581303      -2.219916538917496E-003 -1.645759024605332E-008
  -312260.810600813       2.271101716201623E-004  2.402784255814391E-008
  -752.508288874525      -8.833792956353115E-006 -1.628170969281058E-009
STREST,DECRE           1
  -312255.914057127       5.085040822022346E-007  2.403003333928053E-008
  -729161.657581303      -1.016990247525786E-006 -1.645759024605332E-008
  -312260.810600813       5.084861653235515E-007  2.402784255814391E-008
  -752.508288874525      -2.753493195002849E-009 -1.628170969281058E-009
STREST,DECRE           2
  -312255.914057127       5.085040822022453E-007  2.403003333928053E-008
  -729161.657581303      -1.016990247525807E-006 -1.645759024605332E-008
  -312260.810600813       5.084861653235622E-007  2.402784255814391E-008
  -752.412428092333      -2.753423042417126E-009 -1.628170969281058E-009
STRESS,STRAIN           2           1  1.000000000000000E-003
  -312255.914057127       2.271389122026695E-004  2.862047561066369E-008
  -729161.657581303      -2.219916538917496E-003 -1.621784271014096E-008
  -312260.810600813       2.271101716201623E-004  2.861777658792718E-008
  -752.508288874525      -8.833792956353115E-006 -2.148888764286767E-009
STREST,DECRE           1
  -312255.914057127       5.085040822022346E-007  2.862047561066369E-008
  -729161.657581303      -1.016990247525786E-006 -1.621784271014096E-008
  -312260.810600813       5.084861653235515E-007  2.861777658792718E-008
  -752.508288874525      -2.753493195002849E-009 -2.148888764286767E-009
STREST,DECRE           2
  -312255.914057127       5.085040822022403E-007  2.862047561066369E-008
  -729161.657581303      -1.016990247525797E-006 -1.621784271014096E-008
  -312260.810600813       5.084861653235572E-007  2.861777658792718E-008
  -752.456785534130      -2.753455503963620E-009 -2.148888764286767E-009
STRESS,STRAIN           2           1  1.000000000000000E-003
  -312255.914057127       2.271389122026695E-004  3.014559188844961E-008
  -729161.657581303      -2.219916538917496E-003 -1.473590995920184E-008
  -312260.810600813       2.271101716201623E-004  3.014271516321107E-008
  -752.508288874525      -8.833792956353115E-006 -2.389669615661425E-009
STREST,DECRE           1
  -312255.914057127       5.085040822022346E-007  3.014559188844961E-008
  -729161.657581303      -1.016990247525786E-006 -1.473590995920184E-008
  -312260.810600813       5.084861653235515E-007  3.014271516321107E-008
  -752.508288874525      -2.753493195002849E-009 -2.389669615661425E-009
STREST,DECRE           2
  -312255.914057127       5.085040822022380E-007  3.014559188844961E-008
  -729161.657581303      -1.016990247525793E-006 -1.473590995920184E-008
  -312260.810600813       5.084861653235549E-007  3.014271516321107E-008
  -752.477296495544      -2.753470514241503E-009 -2.389669615661425E-009
STRESS,STRAIN           2           1  1.000000000000000E-003
  -312255.914057127       2.271389122026695E-004  3.078165416848692E-008
  -729161.657581303      -2.219916538917496E-003 -1.347363301190560E-008
  -312260.810600813       2.271101716201623E-004  3.077865554133651E-008
  -752.508288874525      -8.833792956353115E-006 -2.525881364883103E-009
STREST,DECRE           1
  -312255.914057127       5.085040822022346E-007  3.078165416848692E-008
  -729161.657581303      -1.016990247525786E-006 -1.347363301190560E-008
  -312260.810600813       5.084861653235515E-007  3.077865554133651E-008
  -752.508288874525      -2.753493195002849E-009 -2.525881364883103E-009
STREST,DECRE           2
  -312255.914057127       5.085040822022366E-007  3.078165416848692E-008
  -729161.657581303      -1.016990247525790E-006 -1.347363301190560E-008
  -312260.810600813       5.084861653235535E-007  3.077865554133651E-008
  -752.488899718625      -2.753479005681728E-009 -2.525881364883103E-009
STRESS,STRAIN           2           1  1.000000000000000E-003
  -312255.914057127       2.271389122026695E-004  3.107618653158698E-008
  -729161.657581303      -2.219916538917496E-003 -1.257287149643815E-008
  -312260.810600813       2.271101716201623E-004  3.107308397226353E-008
  -752.508288874525      -8.833792956353115E-006 -2.609931072131535E-009
STREST,DECRE           1
  -312255.914057127       5.085040822022346E-007  3.107618653158698E-008
  -729161.657581303      -1.016990247525786E-006 -1.257287149643815E-008
  -312260.810600813       5.084861653235515E-007  3.107308397226353E-008
  -752.508288874525      -2.753493195002849E-009 -2.609931072131535E-009
STREST,DECRE           2
  -312255.914057127       5.085040822022359E-007  3.107618653158698E-008
  -729161.657581303      -1.016990247525788E-006 -1.257287149643815E-008
  -312260.810600813       5.084861653235528E-007  3.107308397226353E-008
  -752.496059508502      -2.753484245340341E-009 -2.609931072131535E-009
STRESS,STRAIN           2           1  1.000000000000000E-003
  -312255.914057127       2.271389122026695E-004  3.122116969646961E-008
  -729161.657581303      -2.219916538917496E-003 -1.197670947896611E-008
  -312260.810600813       2.271101716201623E-004  3.121797995252829E-008
  -752.508288874525      -8.833792956353115E-006 -2.663487493391817E-009
STREST,DECRE           1
  -312255.914057127       5.085040822022346E-007  3.122116969646961E-008
  -729161.657581303      -1.016990247525786E-006 -1.197670947896611E-008
  -312260.810600813       5.084861653235515E-007  3.121797995252829E-008
  -752.508288874525      -2.753493195002849E-009 -2.663487493391817E-009
STREST,DECRE           2
  -312255.914057127       5.085040822022354E-007  3.122116969646961E-008
  -729161.657581303      -1.016990247525787E-006 -1.197670947896611E-008
  -312260.810600813       5.084861653235524E-007  3.121797995252829E-008
  -752.500621722165      -2.753487584047556E-009 -2.663487493391817E-009
STRESS,STRAIN           2           1  1.000000000000000E-003
  -312255.914057127       2.271389122026695E-004  3.129569095778052E-008
  -729161.657581303      -2.219916538917496E-003 -1.160245478507218E-008
  -312260.810600813       2.271101716201623E-004  3.129243197750543E-008
  -752.508288874525      -8.833792956353115E-006 -2.697715299903093E-009
STREST,DECRE           1
  -312255.914057127       5.085040822022346E-007  3.129569095778052E-008
  -729161.657581303      -1.016990247525786E-006 -1.160245478507218E-008
  -312260.810600813       5.084861653235515E-007  3.129243197750543E-008
  -752.508288874525      -2.753493195002849E-009 -2.697715299903093E-009
STREST,DECRE           2
  -312255.914057127       5.085040822022351E-007  3.129569095778052E-008
  -729161.657581303      -1.016990247525787E-006 -1.160245478507218E-008
  -312260.810600813       5.084861653235521E-007  3.129243197750543E-008
  -752.503537424201      -2.753489717808958E-009 -2.697715299903093E-009
STRESS,STRAIN           2           1  1.000000000000000E-003
  -312255.914057127       2.271389122026695E-004  3.133537614940193E-008
  -729161.657581303      -2.219916538917496E-003 -1.137871747212581E-008
  -312260.810600813       2.271101716201623E-004  3.133206487334493E-008
  -752.508288874525      -8.833792956353115E-006 -2.719335562086548E-009
STREST,DECRE           1
  -312255.914057127       5.085040822022346E-007  3.133537614940193E-008
  -729161.657581303      -1.016990247525786E-006 -1.137871747212581E-008
  -312260.810600813       5.084861653235515E-007  3.133206487334493E-008
  -752.508288874525      -2.753493195002849E-009 -2.719335562086548E-009
STREST,DECRE           2
  -312255.914057127       5.085040822022348E-007  3.133537614940193E-008
  -729161.657581303      -1.016990247525786E-006 -1.137871747212581E-008
  -312260.810600813       5.084861653235517E-007  3.133206487334493E-008
  -752.505379150239      -2.753491065616067E-009 -2.719335562086548E-009
STRESS,STRAIN           2           1  1.000000000000000E-003
  -312255.914057127       2.271389122026695E-004  3.135721918636361E-008
  -729161.657581303      -2.219916538917496E-003 -1.125213451271390E-008
  -312260.810600813       2.271101716201623E-004  3.135387001768270E-008
  -752.508288874525      -8.833792956353115E-006 -2.732732516687487E-009
STREST,DECRE           1
  -312255.914057127       5.085040822022346E-007  3.135721918636361E-008
  -729161.657581303      -1.016990247525786E-006 -1.125213451271390E-008
  -312260.810600813       5.084861653235515E-007  3.135387001768270E-008
  -752.508288874525      -2.753493195002849E-009 -2.732732516687487E-009
STREST,DECRE           2
  -312255.914057127       5.085040822022348E-007  3.135721918636361E-008
  -729161.657581303      -1.016990247525786E-006 -1.125213451271390E-008
  -312260.810600813       5.084861653235517E-007  3.135387001768270E-008
  -752.506520372298      -2.753491900782196E-009 -2.732732516687487E-009
STRESS,STRAIN           2           1  1.000000000000000E-003
  -312255.914057127       2.271389122026695E-004  3.136966614618904E-008
  -729161.657581303      -2.219916538917496E-003 -1.118566440455488E-008
  -312260.810600813       2.271101716201623E-004  3.136629049109822E-008
  -752.508288874525      -8.833792956353115E-006 -2.740835983028529E-009
STREST,DECRE           1
  -312255.914057127       5.085040822022346E-007  3.136966614618904E-008
  -729161.657581303      -1.016990247525786E-006 -1.118566440455488E-008
  -312260.810600813       5.084861653235515E-007  3.136629049109822E-008
  -752.508288874525      -2.753493195002849E-009 -2.740835983028529E-009
STREST,DECRE           2
  -312255.914057127       5.085040822022346E-007  3.136966614618904E-008
  -729161.657581303      -1.016990247525786E-006 -1.118566440455488E-008
  -312260.810600813       5.084861653235515E-007  3.136629049109822E-008
  -752.507210667579      -2.753492405952274E-009 -2.740835983028529E-009
STRESS,STRAIN           2           2  1.000000000000000E-003
  -312255.914057127       2.271702818688157E-004  3.136966614618904E-008
  -729161.657581303      -2.219927724581901E-003 -1.118566440455488E-008
  -312260.810600813       2.271415379106535E-004  3.136629049109822E-008
  -752.507210667579      -8.836533792336143E-006 -2.740835983028529E-009
STREST,DECRE           1
  -312255.914057127       5.084501204551296E-007  3.136966614618904E-008
  -729161.657581303      -1.016882325932890E-006 -1.118566440455488E-008
  -312260.810600813       5.084322054777609E-007  3.136629049109822E-008
  -752.507210667579      -2.753197053294231E-009 -2.740835983028529E-009
STREST,DECRE           2
  -312255.914057127       5.084501204551298E-007  3.136966614618904E-008
  -729161.657581303      -1.016882325932891E-006 -1.118566440455488E-008
  -312260.810600813       5.084322054777611E-007  3.136629049109822E-008
  -752.506157687519      -2.753196282786905E-009 -2.740835983028529E-009
STRESS,STRAIN           2           2  1.000000000000000E-003
  -312255.914057127       2.271702818688157E-004  3.138187698329411E-008
  -729161.657581303      -2.219927724581901E-003 -1.112228519160210E-008
  -312260.810600813       2.271415379106535E-004  3.137846577800797E-008
  -752.507210667579      -8.836533792336143E-006 -2.750179755075248E-009
STREST,DECRE           1
  -312255.914057127       5.084501204551296E-007  3.138187698329411E-008
  -729161.657581303      -1.016882325932890E-006 -1.112228519160210E-008
  -312260.810600813       5.084322054777609E-007  3.137846577800797E-008
  -752.507210667579      -2.753197053294231E-009 -2.750179755075248E-009
STREST,DECRE           2
  -312255.914057127       5.084501204551297E-007  3.138187698329411E-008
  -729161.657581303      -1.016882325932890E-006 -1.112228519160210E-008
  -312260.810600813       5.084322054777610E-007  3.137846577800797E-008
  -752.506953638471      -2.753196865215826E-009 -2.750179755075248E-009
STRESS,STRAIN           2           2  1.000000000000000E-003
  -312255.914057127       2.271702818688157E-004  3.139056924124746E-008
  -729161.657581303      -2.219927724581901E-003 -1.110090899237382E-008
  -312260.810600813       2.271415379106535E-004  3.138713477773862E-008
  -752.507210667579      -8.836533792336143E-006 -2.755508618739744E-009
STREST,DECRE           1
  -312255.914057127       5.084501204551296E-007  3.139056924124746E-008
  -729161.657581303      -1.016882325932890E-006 -1.110090899237382E-008
  -312260.810600813       5.084322054777609E-007  3.138713477773862E-008
  -752.507210667579      -2.753197053294231E-009 -2.755508618739744E-009
STREST,DECRE           2
  -312255.914057127       5.084501204551296E-007  3.139056924124746E-008
  -729161.657581303      -1.016882325932890E-006 -1.110090899237382E-008
  -312260.810600813       5.084322054777609E-007  3.138713477773862E-008
  -752.507407578709      -2.753197197381925E-009 -2.755508618739744E-009
STRESS,STRAIN           2           3  1.000000000000000E-003
  -312255.914057127       2.272016724380569E-004  3.139056924124746E-008
  -729161.657581303      -2.219938825490893E-003 -1.110090899237382E-008
  -312260.810600813       2.271729250454312E-004  3.138713477773862E-008
  -752.507407578709      -8.839289300954882E-006 -2.755508618739744E-009
STREST,DECRE           1
  -312255.914057127       5.083961646017836E-007  3.139056924124746E-008
  -729161.657581303      -1.016774416127305E-006 -1.110090899237382E-008
  -312260.810600813       5.083782515255214E-007  3.138713477773862E-008
  -752.507407578709      -2.752905609114510E-009 -2.755508618739744E-009
STREST,DECRE           2
  -312255.914057127       5.083961646017834E-007  3.139056924124746E-008
  -729161.657581303      -1.016774416127304E-006 -1.110090899237382E-008
  -312260.810600813       5.083782515255212E-007  3.138713477773862E-008
  -752.507629316566      -2.752905771351687E-009 -2.755508618739744E-009
STRESS,STRAIN           2           3  1.000000000000000E-003
  -312255.914057127       2.272016724380569E-004  3.139237265062094E-008
  -729161.657581303      -2.219938825490893E-003 -1.112991890173692E-008
  -312260.810600813       2.271729250454312E-004  3.138892725219766E-008
  -752.507407578709      -8.839289300954882E-006 -2.756309367996016E-009
STREST,DECRE           1
  -312255.914057127       5.083961646017836E-007  3.139237265062094E-008
  -729161.657581303      -1.016774416127305E-006 -1.112991890173692E-008
  -312260.810600813       5.083782515255214E-007  3.138892725219766E-008
  -752.507407578709      -2.752905609114510E-009 -2.756309367996016E-009
STREST,DECRE           2
  -312255.914057127       5.083961646017834E-007  3.139237265062094E-008
  -729161.657581303      -1.016774416127304E-006 -1.112991890173692E-008
  -312260.810600813       5.083782515255212E-007  3.138892725219766E-008
  -752.507697528540      -2.752905821259800E-009 -2.756309367996016E-009

[ 本帖最后由 forcamingseu 于 2007-5-24 10:52 编辑 ]

forcamingseu

为什么
C
C    CALCULATE UPDATED STRESS
C
        DO K1=1,NTENS
          DO K2=1,NTENS
            STREST(K2)=STRESS(K2)+DDSDDE(K2,K1)*DEELA(K1)
          END DO
        END DO
只会更新剪切应力，而直接应力不更新？到底什么原因，大侠们赐教！

aakkaa

STREST(K2)=STRESS(K2)+DDSDDE(K2,K1)*DEELA(K1)
---------------
STREST和STRESS是什么关系?

forcamingseu

已经找到原因了，我这样更新应力STREST(K2)=STRESS(K2)+DDSDDE(K2,K1)*DEELA(K1)是错误的！改正了以后就好了!

wwwaba

TIME(1)+DTIME

forcamingseu
请问你这里的TIME(1)+DTIME 是什么意思?

wwwaba

ding !!!!   
TIME(1)+DTIME

forcamingseu
请问你这里的TIME(1)+DTIME 是什么意思? 为何不用DTIME?????????????????????????