Télécharger maxro2.eso

Retour à la liste

Numérotation des lignes :

  1. C MAXRO2 SOURCE BP208322 17/03/01 21:17:51 9325
  2. SUBROUTINE MAXRO2(ICAS,WTRAV,WRK1,WRK5,WR12,NCHAIN)
  3. *
  4. IMPLICIT INTEGER(I-N)
  5. IMPLICIT REAL*8(A-H,O-Z)
  6. *
  7. -INC CCOPTIO
  8. -INC SMEVOLL
  9. -INC SMLREEL
  10. -INC CCHAMP
  11. *
  12. SEGMENT WRK1
  13. REAL*8 DDHOOK(LHOOK,LHOOK),SIG0(NSTRS),DEPST(NSTRS)
  14. REAL*8 SIGF(NSTRS),VAR0(NVARI),VARF(NVARI)
  15. REAL*8 DEFP(NSTRS),XCAR(ICARA)
  16. ENDSEGMENT
  17. *
  18. *
  19. SEGMENT WRK5
  20. REAL*8 EPIN0(NSTRS),EPINF(NSTRS),EPST0(NSTRS)
  21. ENDSEGMENT
  22. *
  23. SEGMENT WR12
  24. REAL*8 EM0(2,NWA(1)),EM1(2,NWA(2)),EM2(2,NWA(3))
  25. REAL*8 EM3(2,NWA(4)),EM4(2,NWA(5)),EM5(2,NWA(6))
  26. REAL*8 EM6(2,NWA(7)),EM7(2,NWA(8)),EM8(2,NWA(9))
  27. REAL*8 SM0(NSTRS),SM1(NSTRS),SM2(NSTRS),SM3(NSTRS)
  28. REAL*8 SM4(NSTRS),SM5(NSTRS),SM6(NSTRS),SM7(NSTRS)
  29. REAL*8 SM8(NSTRS)
  30. ENDSEGMENT
  31. *
  32. SEGMENT WTRAV
  33. REAL*8 DDAUX(LHOOK,LHOOK),VALMAT(NUMAT)
  34. REAL*8 VALCAR(NUCAR),DSIGT(NSTRS)
  35. REAL*8 TXR(IDIM,IDIM),DDHOMU(LHOOK,LHOOK)
  36. REAL*8 XLOC(3,3),XGLOB(3,3)
  37. REAL*8 D1HOOK(LHOOK,LHOOK),ROTHOO(LHOOK,LHOOK)
  38. ENDSEGMENT
  39. *
  40. SEGMENT WTRA2
  41. REAL*8 A(NDIM,NDIM)
  42. REAL*8 RR(NDIM,NDIM),RT(NDIM,NDIM),RTRA(NDIM,NDIM)
  43. ENDSEGMENT
  44. *
  45. IDIM=XGLOB(/1)
  46. NDIM=IDIM
  47. IF(IFOUR.EQ.1) NDIM=IDIM+1
  48.  
  49. NSTRS=SIG0(/1)
  50.  
  51. SEGINI WTRA2
  52.  
  53. * ICAS = 1 : RT A R
  54. * ICAS = 2 : R A RT
  55.  
  56. DO I=1,IDIM
  57. DO J=1,IDIM
  58. RR(I,J)=XGLOB(I,J)
  59. ENDDO
  60. ENDDO
  61. IF(IDIM.EQ.2.AND.IFOUR.EQ.1) RR(3,3)=1.D0
  62. CALL TRSPOD(RR,NDIM,NDIM,RT)
  63. *
  64. * ROTATION DES TENSEURS
  65. *
  66. IDEFO=0
  67. CALL MAXRO3(ICAS,IDEFO,SM0,A,RR,RT,RTRA,NDIM,IDIM,IFOUR)
  68. CALL MAXRO3(ICAS,IDEFO,SM1,A,RR,RT,RTRA,NDIM,IDIM,IFOUR)
  69. CALL MAXRO3(ICAS,IDEFO,SM2,A,RR,RT,RTRA,NDIM,IDIM,IFOUR)
  70. CALL MAXRO3(ICAS,IDEFO,SM3,A,RR,RT,RTRA,NDIM,IDIM,IFOUR)
  71. CALL MAXRO3(ICAS,IDEFO,SM4,A,RR,RT,RTRA,NDIM,IDIM,IFOUR)
  72. IF(NCHAIN.GE.6) THEN
  73. CALL MAXRO3(ICAS,IDEFO,SM5,A,RR,RT,RTRA,NDIM,IDIM,IFOUR)
  74. ENDIF
  75. IF(NCHAIN.GE.7) THEN
  76. CALL MAXRO3(ICAS,IDEFO,SM6,A,RR,RT,RTRA,NDIM,IDIM,IFOUR)
  77. ENDIF
  78. IF(NCHAIN.GE.8) THEN
  79. CALL MAXRO3(ICAS,IDEFO,SM7,A,RR,RT,RTRA,NDIM,IDIM,IFOUR)
  80. ENDIF
  81. IF(NCHAIN.GE.9) THEN
  82. CALL MAXRO3(ICAS,IDEFO,SM8,A,RR,RT,RTRA,NDIM,IDIM,IFOUR)
  83. ENDIF
  84.  
  85.  
  86. IF(ICAS.EQ.1) THEN
  87.  
  88. CALL MAXRO3(ICAS,IDEFO,SIG0,A,RR,RT,RTRA,NDIM,IDIM,IFOUR)
  89. IDEFO=1
  90. CALL MAXRO3(ICAS,IDEFO,DEPST,A,RR,RT,RTRA,NDIM,IDIM,IFOUR)
  91. CALL MAXRO3(ICAS,IDEFO,EPIN0,A,RR,RT,RTRA,NDIM,IDIM,IFOUR)
  92. *
  93.  
  94. ELSE IF(ICAS.EQ.2) THEN
  95.  
  96. CALL MAXRO3(ICAS,IDEFO,SIGF,A,RR,RT,RTRA,NDIM,IDIM,IFOUR)
  97. IDEFO=1
  98. CALL MAXRO3(ICAS,IDEFO,DEFP ,A,RR,RT,RTRA,NDIM,IDIM,IFOUR)
  99. CALL MAXRO3(ICAS,IDEFO,EPINF,A,RR,RT,RTRA,NDIM,IDIM,IFOUR)
  100.  
  101. ENDIF
  102.  
  103. SEGSUP WTRA2
  104.  
  105. RETURN
  106. END
  107.  
  108.  
  109.  
  110.  
  111.  
  112.  
  113.  
  114.  
  115.  
  116.  

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