Télécharger modmax.eso

Retour à la liste

Numérotation des lignes :

  1. C MODMAX SOURCE CB215821 17/11/30 21:16:53 9639
  2. SUBROUTINE MODMAX(TP,EU,M,CODU,CODL,COD,BETJEF,BETFLU)
  3. C
  4. IMPLICIT INTEGER(I-N)
  5. IMPLICIT REAL*8(A-H,O-Z)
  6. DIMENSION CODU(9,9),CODL(8,8),COD(8)
  7. C
  8. *
  9. SEGMENT BETJEF
  10. REAL*8 AA,BETA,FC,PALF,YOUN,XNU,GFC,GFT,CAR,ETA,TDEF,
  11. & TCON,DPSTF1,DPSTF2,TETA,PDT,TP00
  12. INTEGER ICT,ICC,IMOD,IVIS,ITR,
  13. & ISIM,IBB,IGAU,IZON
  14. ENDSEGMENT
  15. SEGMENT BETFLU
  16. REAL*8 DATCOU,DATCUR,DATSEC,E28,PGTZO,PGDUR,TAU1,TAU2,
  17. & TP0,TZER
  18. INTEGER ITYPE,IMD,NBRC,NCOE,NTZERO,NTPS,IFOR
  19. ENDSEGMENT
  20. C
  21. NC = NCOE+1
  22. C
  23. C----------------------------------------------------------------
  24. IF (IFOR.EQ.1) THEN
  25. C Modele visco elasto plastique solide de MAXWELL
  26. C----------------------------------------------------------------
  27. EU=0.D0
  28. C
  29. DO 20 J=1,M
  30. C
  31. IF (J.EQ.M) THEN
  32. C
  33. C
  34. DO 30 N=1,NC
  35. IF(N.EQ.1) THEN
  36. C
  37. EU=EU+CODU(J,N)
  38. C
  39. ELSE
  40. C
  41. EUD=0.D0
  42. EUD=CODU(J,N)*EXP((DATCUR-(TP/86400-DATCOU))
  43. * *COD(N-1))
  44. EU=EU+EUD
  45. C
  46. ENDIF
  47. 30 CONTINUE
  48. C
  49. C ELSE
  50. C WRITE(*,*)'PAS DE VALEURS TROUVE POUR LE MODULE'
  51. C
  52. ENDIF
  53. C
  54. 20 CONTINUE
  55. C
  56. C
  57. C MODIFICATION DU 19/01/00
  58. C********************************************************************
  59. C
  60. C----------------------------------------------------------
  61. ELSE IF (IFOR.EQ.2) THEN
  62. C Modele visco elasto plastique liquide de MAXWELL
  63. C----------------------------------------------------------
  64. C
  65. EU=0.D0
  66. C
  67. DO 40 J=1,M
  68. C
  69. IF (J.EQ.M) THEN
  70. C
  71. DO 50 N=1,NCOE
  72. EUD=0.D0
  73. EUD=CODL(J,N)*EXP((DATCUR-(TP/86400-DATCOU))
  74. * *COD(N))
  75. EU=EU+EUD
  76. 50 CONTINUE
  77. C
  78. C ELSE
  79. C WRITE(*,*)'PAS DE VALEURS TROUVE POUR LE MODULE'
  80. C
  81. ENDIF
  82. C
  83. 40 CONTINUE
  84. C
  85. C------------------------------------------------------------
  86. ENDIF
  87. C------------------------------------------------------------
  88. C
  89. RETURN
  90. END
  91.  
  92.  
  93.  
  94.  

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