Télécharger cmaxob.eso

Retour à la liste

Numérotation des lignes :

  1. C CMAXOB SOURCE BP208322 17/03/01 21:15:42 9325
  2. SUBROUTINE CMAXOB(wrk52,wrk53,wrk54,WR12,NWA,NCHAIN,EPSFLU)
  3. *
  4. * MODELE MAXOTT : 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. -INC DECHE
  15.  
  16. SEGMENT WR12
  17. REAL*8 EM0(2,NWA(1)),EM1(2,NWA(2)),EM2(2,NWA(3))
  18. REAL*8 EM3(2,NWA(4)),EM4(2,NWA(5)),EM5(2,NWA(6))
  19. REAL*8 EM6(2,NWA(7)),EM7(2,NWA(8)),EM8(2,NWA(9))
  20. REAL*8 SM0(NSTRS),SM1(NSTRS),SM2(NSTRS),SM3(NSTRS)
  21. REAL*8 SM4(NSTRS),SM5(NSTRS),SM6(NSTRS),SM7(NSTRS)
  22. REAL*8 SM8(NSTRS)
  23. ENDSEGMENT
  24.  
  25. DIMENSION NWA(9)
  26. DIMENSION EPSFLU(8)
  27.  
  28. *
  29. * ON CREE LES PROGS ASSOCIES AUX VARIABLES INTERNES
  30. *
  31. * Formulation tridimensionnelle massive
  32. *
  33. IF(IFOUR.EQ.2.AND.MFR.EQ.1) THEN
  34. MXVREE = 22
  35. *
  36. * Formulation deformations planes / axisymetrique
  37. *
  38. ELSE IF(IFOUR.EQ.-1.OR.IFOUR.EQ.-3
  39. & .OR.IFOUR.EQ.0.OR.IFOUR.EQ.1) THEN
  40. MXVREE = 17
  41. *
  42. * Formulation contraintes planes / 3D coques minces
  43. *
  44. ELSE
  45. MXVREE = 14
  46. ENDIF
  47. DO JC=0,NCHAIN
  48. JG=NSTRS
  49. SEGINI MLREEL
  50. VARF(JC+MXVREE)=DBLE(MLREEL)
  51. IF(JC.EQ.0) THEN
  52. DO JD=1,NSTRS
  53. PROG(JD)=EPSFLU(JD)
  54. END DO
  55. ELSE IF(JC.EQ.1) THEN
  56. DO JD=1,NSTRS
  57. PROG(JD)=SM0(JD)
  58. END DO
  59. ELSE IF(JC.EQ.2) THEN
  60. DO JD=1,NSTRS
  61. PROG(JD)=SM1(JD)
  62. END DO
  63. ELSE IF(JC.EQ.3) THEN
  64. DO JD=1,NSTRS
  65. PROG(JD)=SM2(JD)
  66. END DO
  67. ELSE IF(JC.EQ.4) THEN
  68. DO JD=1,NSTRS
  69. PROG(JD)=SM3(JD)
  70. END DO
  71. ELSE IF(JC.EQ.5) THEN
  72. DO JD=1,NSTRS
  73. PROG(JD)=SM4(JD)
  74. END DO
  75. ELSE IF(JC.EQ.6) THEN
  76. DO JD=1,NSTRS
  77. PROG(JD)=SM5(JD)
  78. END DO
  79. ELSE IF(JC.EQ.7) THEN
  80. DO JD=1,NSTRS
  81. PROG(JD)=SM6(JD)
  82. END DO
  83. ELSE IF(JC.EQ.8) THEN
  84. DO JD=1,NSTRS
  85. PROG(JD)=SM7(JD)
  86. END DO
  87. ELSE IF(JC.EQ.9) THEN
  88. DO JD=1,NSTRS
  89. PROG(JD)=SM8(JD)
  90. END DO
  91. ENDIF
  92. SEGDES MLREEL
  93.  
  94. END DO
  95. DO JC=NCHAIN+1,9
  96. JG=0
  97. SEGINI MLREEL
  98. VARF(JC+MXVREE)=DBLE(MLREEL)
  99. SEGDES MLREEL
  100. END DO
  101. *
  102. RETURN
  103. END
  104.  
  105.  
  106.  
  107.  
  108.  
  109.  
  110.  
  111.  
  112.  
  113.  
  114.  
  115.  
  116.  
  117.  
  118.  
  119.  
  120.  
  121.  
  122.  
  123.  
  124.  
  125.  
  126.  
  127.  
  128.  
  129.  
  130.  
  131.  
  132.  
  133.  
  134.  
  135.  
  136.  
  137.  
  138.  

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