Télécharger la plaquette Cast3M


fr en

Télécharger savcr1.eso

Retour à la liste

Numérotation des lignes : 

savcr1

  1. C SAVCR1    SOURCE    FD218221  26/02/13    21:15:29     12462          
  2.       SUBROUTINE SAVCR1(VARF,NVARI,CMAT,NMAT,SIGF,
  3.      $           SIGEF,NSTRS,IFOUR,NFLU,EPSFLU,AN,BN,CN)
  4. C
  5. C_______________________________________________________________________
  6. C
  7. C  SOUS PROGRAMME APPELE PAR SAVCR QUI CALCUL LES MATRICES DE 
  8. C  DEFORMATION DE FLUAGE PROPRE D APRES LE MODELE DE KELVIN
  9. C  A LA FIN DU PAS DE TEMPS
  10. C_______________________________________________________________________
  11. C
  12.       IMPLICIT INTEGER(I-N)
  13.       IMPLICIT REAL*8(A-H,O-Z)
  14.       DIMENSION VARF(NVARI),CMAT(NMAT),SIGF(NSTRS),SIGEF(NSTRS)
  15.       DIMENSION EPSFLU(6),FLUO(6),VEPSM(6),VEPSO(6),VEPSN(6),SIGFI(10)
  16.       DIMENSION AN(NSTRS),BN(NSTRS,NSTRS),CN(NSTRS,NSTRS)
  17. C
  18. C  TEST DE DIMENSION
  19. C
  20.       IF(NSTRS.GT.10) THEN
  21.           PRINT *, 'SAUFLD - DIMENSION INSUFFISANTE'
  22.           RETURN 
  23.       ENDIF
  24. C
  25. C RECUPERATION DES VALEURS DES PARAMETRES MATERIAU
  26. C
  27.       TDT  = CMAT(28)
  28. C
  29. C RECUPERATION DES VARIABLES INTERNES INITIALES
  30. C
  31. C   FLUAGE PROPRE
  32. C
  33.       EPSFLU(1) = VARF(9)
  34.       EPSFLU(2) = VARF(10)
  35.       EPSFLU(3) = VARF(11)
  36.       EPSFLU(4) = VARF(12)
  37.       EPSFLU(5) = VARF(13)
  38.       EPSFLU(6) = VARF(14)
  39. C
  40. C   VITESSE N
  41. C
  42.       VEPSO(1) = VARF(27)
  43.       VEPSO(2) = VARF(28)
  44.       VEPSO(3) = VARF(29)
  45.       VEPSO(4) = VARF(30)
  46.       VEPSO(5) = VARF(31)
  47.       VEPSO(6) = VARF(32)
  48. C
  49. C   FLUAGE PROPRE  N
  50. C
  51.       FLUO(1) = VARF(9)
  52.       FLUO(2) = VARF(10)
  53.       FLUO(3) = VARF(11)
  54.       FLUO(4) = VARF(12)
  55.       FLUO(5) = VARF(13)
  56.       FLUO(6) = VARF(14)
  57. C
  58. C  TEST SI LE TEMPS EST NUL
  59. C
  60.       IF (TDT.EQ.0.D0) THEN
  61.         GOTO 30
  62.       ENDIF
  63. C
  64. C INITIALISATION DES MATRICES
  65. C
  66.       CALL ZERO(AN,NSTRS,1)
  67.       CALL ZERO(BN,NSTRS,NSTRS)
  68.       CALL ZERO(CN,NSTRS,NSTRS)
  69.       CALL ZERO(SIGFI,NSTRS,1)
  70. C
  71.       DO 10 I=1,NSTRS
  72.           SIGFI(I) = SIGF(I)
  73.   10  CONTINUE
  74. C
  75. C_______________________________________________________________________
  76. C
  77. C  FLUAGE PROPRE KELVIN VOIGT
  78. C
  79.       NDI = NSTRS
  80.       CALL DEFCR1(VARF,NVARI,CMAT,NMAT,IFOUR,NSTRS,
  81.      $                 AN,BN,CN,NFLU,NDI)
  82. C
  83.       DO 20 I=1,NSTRS
  84.         EPSFLU(I) = EPSFLU(I) + AN(I)
  85.         DO 21 J=1,NSTRS
  86.           EPSFLU(I) = EPSFLU(I) + BN(I,J)*SIGEF(J) + CN(I,J)*SIGFI(J)
  87.   21    CONTINUE
  88.   20  CONTINUE
  89. C
  90. C
  91. C  VITESSE MOYENNE DE DEFORMATION
  92. C
  93.       VEPSM(1) = (EPSFLU(1) -  FLUO(1))/TDT
  94.       VEPSM(2) = (EPSFLU(2) -  FLUO(2))/TDT
  95.       VEPSM(3) = (EPSFLU(3) -  FLUO(3))/TDT
  96.       VEPSM(4) = (EPSFLU(4) -  FLUO(4))/TDT
  97.       VEPSM(5) = (EPSFLU(5) -  FLUO(5))/TDT
  98.       VEPSM(6) = (EPSFLU(6) -  FLUO(6))/TDT
  99. C
  100. C  VITESSE DE DEFORMATION (N+1)
  101. C
  102.       VEPSN(1) = (2.D0 * VEPSM(1)) - VEPSO(1)
  103.       VEPSN(2) = (2.D0 * VEPSM(2)) - VEPSO(2)
  104.       VEPSN(3) = (2.D0 * VEPSM(3)) - VEPSO(3)
  105.       VEPSN(4) = (2.D0 * VEPSM(4)) - VEPSO(4)
  106.       VEPSN(5) = (2.D0 * VEPSM(5)) - VEPSO(5)
  107.       VEPSN(6) = (2.D0 * VEPSM(6)) - VEPSO(6)
  108.  
  109. C
  110.   30  CONTINUE
  111. C
  112. C  SAUVEGARDE DES DEFORMATIONS ET VITESSE (N+1)
  113. C
  114. C    FLUAGE PROPRE
  115. C
  116.       VARF(9)  = EPSFLU(1)
  117.       VARF(10) = EPSFLU(2)
  118.       VARF(11) = EPSFLU(3)
  119.       VARF(12) = EPSFLU(4)
  120.       VARF(13) = EPSFLU(5)
  121.       VARF(14) = EPSFLU(6)
  122. C
  123.       VARF(27) = VEPSN(1)
  124.       VARF(28) = VEPSN(2)
  125.       VARF(29) = VEPSN(3)
  126.       VARF(30) = VEPSN(4)
  127.       VARF(31) = VEPSN(5)
  128.       VARF(32) = VEPSN(6)
  129. C
  130.       END
  131.  
  132.  

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