Télécharger la plaquette Cast3M


fr en

Télécharger saufld.eso

Retour à la liste

Numérotation des lignes : 

saufld

  1. C SAUFLD    SOURCE    AM        09/12/08    21:15:27     6582
  2.       SUBROUTINE SAUFLD(VARF,NVARI,CMAT,NMAT,SIGF,
  3.      $           SIGEF,NSTRS,IFOUR,NFLU,EPSFLU,AN,BN,CN)
  4. C
  5. C_______________________________________________________________________
  6. C
  7. C SOUS PROGRAMME QUI CALCUL LA DEFORMATION
  8. C DE FLUAGE PROPRE A LA FIN DU PAS DE TEMPS
  9. C_______________________________________________________________________
  10. C
  11.       IMPLICIT INTEGER(I-N)
  12.       IMPLICIT REAL*8(A-H,O-Z)
  13.       DIMENSION VARF(NVARI),CMAT(NMAT),SIGF(NSTRS),SIGEF(NSTRS)
  14.       DIMENSION EPSFLU(6),SIGFI(10)
  15.       DIMENSION AN(NSTRS),BN(NSTRS,NSTRS),CN(NSTRS,NSTRS)
  16.  
  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.  
  25. C
  26. C RECUPERATION DES VALEURS DES PARAMETRES MATERIAU
  27. C
  28.       TDT  = CMAT(NMAT)
  29. C
  30. C RECUPERATION DES VARIABLES INTERNES INITIALES
  31. C
  32. C   FLUAGE PROPRE
  33. C
  34.       EPSFLU(1) = VARF(3 + (6 * NFLU) )
  35.       EPSFLU(2) = VARF(4 + (6 * NFLU) )
  36.       EPSFLU(3) = VARF(5 + (6 * NFLU) )
  37.       EPSFLU(4) = VARF(6 + (6 * NFLU) )
  38.       EPSFLU(5) = VARF(7 + (6 * NFLU) )
  39.       EPSFLU(6) = VARF(8 + (6 * NFLU) )
  40. C
  41. C  TEST SI LE TEMPS EST NUL
  42. C
  43.       IF (TDT.EQ.0.D0) THEN
  44.         GOTO 30
  45.       ENDIF
  46. C
  47. C INITIALISATION DES MATRICES
  48. C
  49.       CALL ZERO(AN,NSTRS,1)
  50.       CALL ZERO(BN,NSTRS,NSTRS)
  51.       CALL ZERO(CN,NSTRS,NSTRS)
  52.       CALL ZERO(SIGFI,NSTRS,1)
  53. C
  54.       DO 10 I=1,NSTRS
  55.           SIGFI(I) = SIGF(I)
  56.   10  CONTINUE
  57. C
  58. C_______________________________________________________________________
  59. C
  60. C  FLUAGE PROPRE KELVIN VOIGT
  61. C
  62.       NDI = NSTRS
  63.       CALL FLUD(VARF,NVARI,CMAT,NMAT,IFOUR,NSTRS,
  64.      $                 AN,BN,CN,NFLU,NDI)
  65. C
  66.       DO 20 I=1,NSTRS
  67.         EPSFLU(I) = EPSFLU(I) + AN(I)
  68.         DO 20 J=1,NSTRS
  69.           EPSFLU(I) = EPSFLU(I) + BN(I,J)*SIGEF(J) + CN(I,J)*SIGFI(J)
  70.   20  CONTINUE
  71. C
  72.   30  CONTINUE
  73. C
  74. C  SAUVEGARDE DES DEFORMATIONS
  75. C
  76. C    FLUAGE PROPRE
  77. C
  78.       VARF(3 + (6 * NFLU) ) = EPSFLU(1)
  79.       VARF(4 + (6 * NFLU) ) = EPSFLU(2)
  80.       VARF(5 + (6 * NFLU) ) = EPSFLU(3)
  81.       VARF(6 + (6 * NFLU) ) = EPSFLU(4)
  82.       VARF(7 + (6 * NFLU) ) = EPSFLU(5)
  83.       VARF(8 + (6 * NFLU) ) = EPSFLU(6)
  84. C
  85.       END
  86.  
  87.  
  88.  

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