Télécharger la plaquette Cast3M


fr en

Télécharger modmax.eso

Retour à la liste

Numérotation des lignes : 

modmax

  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