Télécharger inva.eso

Retour à la liste

Numérotation des lignes :

  1. C INVA SOURCE CB215821 16/04/21 21:17:13 8920
  2. C INVA SOURCE INSL 24/10/96
  3. SUBROUTINE INVA(S1,DJ2,DI1,CO3T,DJ3,NSTRS)
  4. C ****************************************************
  5. C * CALCUL DES IVARIANT : I1, J2, J3 et COS(3TETA) *
  6. C * ( MODEL D'OTTOSEN ) *
  7. C ****************************************************
  8. IMPLICIT INTEGER(I-N)
  9. IMPLICIT REAL*8(A-H,O-Z)
  10. DIMENSION S1(NSTRS)
  11. C
  12. CALL PARINV(S1,SX,SY,SZ,SXY,SXZ,SYZ,NSTRS)
  13. C
  14. DI1=SX+SY+SZ
  15. SSX=SX-DI1/3.D0
  16. SSY=SY-DI1/3.D0
  17. SSZ=SZ-DI1/3.D0
  18. C
  19. DJ2=(SSX*SSX+SSY*SSY+SSZ*SSZ)/2.D0+SXY*SXY+SXZ*SXZ+SYZ*SYZ
  20. DJ=SSX*SSY*SSZ+2.D0*SXY*SXZ*SYZ-SSX*SYZ*SYZ-SSY*SXZ*SXZ
  21. DJ3=DJ-SSZ*SXY*SXY
  22. CO3T=1.5D0*SQRT(3.D0)*DJ3/SQRT(DJ2*DJ2*DJ2)
  23. C
  24. IF(CO3T .LT. -1.D0) CO3T=-1.D0
  25. IF(CO3T .GT. 1.D0) CO3T=1.D0
  26. C---------------------------------------------------------------------
  27. IECR=0
  28. IF(IECR.EQ.3) THEN
  29. WRITE(*,100) SX,SY,SZ,SXY,SXZ,SYZ,DI1,DJ2,DJ3,CO3T
  30. 100 FORMAT('SX=',E9.3,1X,'SY=',E9.3,1X,'SZ=',E9.3,1X,'SXY=',E9.3,1X
  31. *,'SXZ=',E9.3,1X,'SYZ=',E9.3,2X,/,'DI1=',E9.3,2X,'DJ2=',E9.3,2X,
  32. * 'DJ3=',E9.3,1X,'CO3T=',E10.4)
  33. ENDIF
  34. C---------------------------------------------------------------------
  35. RETURN
  36. END
  37.  
  38.  
  39.  
  40.  
  41.  

© Cast3M 2003 - Tous droits réservés.
Mentions légales