Télécharger maxtrb.eso

Retour à la liste

Numérotation des lignes :

  1. C MAXTRB SOURCE BP208322 17/03/01 21:17:54 9325
  2. SUBROUTINE MAXTRB(WTRAV,WRK1,WRK5,WR12,NWA,NSTRS,NCHAIN,CMATE)
  3. *
  4. * MODELE DE MAXWELL : RECUPERATION DES VARIABLES INTERNES
  5. * APRES ECOULEMENT
  6. *
  7. IMPLICIT INTEGER(I-N)
  8. IMPLICIT REAL*8(A-H,O-Z)
  9. *
  10. -INC CCOPTIO
  11. -INC SMEVOLL
  12. -INC SMLREEL
  13. -INC CCHAMP
  14. *
  15. SEGMENT WRK1
  16. REAL*8 DDHOOK(LHOOK,LHOOK),SIG0(NSTRS),DEPST(NSTRS)
  17. REAL*8 SIGF(NSTRS),VAR0(NVARI),VARF(NVARI)
  18. REAL*8 DEFP(NSTRS),XCAR(ICARA)
  19. ENDSEGMENT
  20. *
  21. SEGMENT WRK5
  22. REAL*8 EPIN0(NSTRS),EPINF(NSTRS),EPST0(NSTRS)
  23. ENDSEGMENT
  24. *
  25. SEGMENT WR12
  26. REAL*8 EM0(2,NWA(1)),EM1(2,NWA(2)),EM2(2,NWA(3))
  27. REAL*8 EM3(2,NWA(4)),EM4(2,NWA(5)),EM5(2,NWA(6))
  28. REAL*8 EM6(2,NWA(7)),EM7(2,NWA(8)),EM8(2,NWA(9))
  29. REAL*8 SM0(NSTRS),SM1(NSTRS),SM2(NSTRS),SM3(NSTRS)
  30. REAL*8 SM4(NSTRS),SM5(NSTRS),SM6(NSTRS),SM7(NSTRS)
  31. REAL*8 SM8(NSTRS)
  32. ENDSEGMENT
  33. *
  34.  
  35.  
  36.  
  37. SEGMENT WTRAV
  38. REAL*8 DDAUX(LHOOK,LHOOK),VALMAT(NUMAT)
  39. REAL*8 VALCAR(NUCAR),DSIGT(NSTRS)
  40. REAL*8 TXR(IDIM,IDIM),DDHOMU(LHOOK,LHOOK)
  41. REAL*8 XLOC(3,3),XGLOB(3,3)
  42. REAL*8 D1HOOK(LHOOK,LHOOK),ROTHOO(LHOOK,LHOOK)
  43. ENDSEGMENT
  44. *
  45.  
  46. CHARACTER*8 CMATE
  47.  
  48.  
  49.  
  50. DIMENSION NWA(9)
  51.  
  52. *
  53. * ON SE REMET DANS LES BONS AXES SI BESOIN
  54. *
  55. *
  56. ******* CAS D'UNE FORMULATION UNIDIRECTIONNELLE
  57. *
  58. IF(CMATE.EQ.'UNIDIREC') THEN
  59. ICAS=2
  60. CALL MAXRO2(ICAS,WTRAV,WRK1,WRK5,WR12,NCHAIN)
  61. IF(IERR.NE.0) RETURN
  62. ENDIF
  63.  
  64. *
  65. * ON CREE LES PROGS ASSOCIES AUX VARIABLES INTERNES
  66. *
  67. DO JC=1,NCHAIN
  68. * zzzzzzzz il faudra modifier le test sur var0
  69.  
  70.  
  71. JG=NSTRS
  72. SEGINI MLREEL
  73. VARF(JC+1)=DBLE(MLREEL)
  74. IF(JC.EQ.1) THEN
  75. DO JD=1,NSTRS
  76. PROG(JD)=SM0(JD)
  77. END DO
  78. ELSE IF(JC.EQ.2) THEN
  79. DO JD=1,NSTRS
  80. PROG(JD)=SM1(JD)
  81. END DO
  82. ELSE IF(JC.EQ.3) THEN
  83. DO JD=1,NSTRS
  84. PROG(JD)=SM2(JD)
  85. END DO
  86. ELSE IF(JC.EQ.4) THEN
  87. DO JD=1,NSTRS
  88. PROG(JD)=SM3(JD)
  89. END DO
  90. ELSE IF(JC.EQ.5) THEN
  91. DO JD=1,NSTRS
  92. PROG(JD)=SM4(JD)
  93. END DO
  94. ELSE IF(JC.EQ.6) THEN
  95. DO JD=1,NSTRS
  96. PROG(JD)=SM5(JD)
  97. END DO
  98. ELSE IF(JC.EQ.7) THEN
  99. DO JD=1,NSTRS
  100. PROG(JD)=SM6(JD)
  101. END DO
  102. ELSE IF(JC.EQ.8) THEN
  103. DO JD=1,NSTRS
  104. PROG(JD)=SM7(JD)
  105. END DO
  106. ELSE IF(JC.EQ.9) THEN
  107. DO JD=1,NSTRS
  108. PROG(JD)=SM8(JD)
  109. END DO
  110. ENDIF
  111. SEGDES MLREEL
  112.  
  113. END DO
  114. DO JC=NCHAIN+1,9
  115. JG=0
  116. SEGINI MLREEL
  117. VARF(JC+1)=DBLE(MLREEL)
  118. SEGDES MLREEL
  119. END DO
  120. *
  121. RETURN
  122. END
  123.  
  124.  
  125.  
  126.  
  127.  
  128.  
  129.  
  130.  
  131.  
  132.  
  133.  
  134.  
  135.  
  136.  
  137.  
  138.  
  139.  
  140.  
  141.  
  142.  
  143.  
  144.  
  145.  
  146.  

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