Télécharger la plaquette Cast3M


fr en

Télécharger flud.eso

Retour à la liste

Numérotation des lignes : 

flud

  1. C FLUD      SOURCE    CB215821  16/04/21    21:16:54     8920
  2.       SUBROUTINE FLUD(VARF,NVARI,CMAT,NMAT,IFOUR,NSTRS,
  3.      $                       AN,BN,CN,NFLU,NDI)
  4. C
  5. C_______________________________________________________________________
  6. C
  7. C SOUS PROGRAMME QUI CALCUL LES MATRICES DE DEFORMATION
  8. C  DE FLUAGE PROPRE
  9. C   D APRES LE MODELE DE KELVIN
  10. C_______________________________________________________________________
  11. C
  12.       IMPLICIT INTEGER(I-N)
  13.       IMPLICIT REAL*8(A-H,O-Z)
  14.       DIMENSION VARF(NVARI),CMAT(NMAT),EPSFLU(6)
  15.       DIMENSION AN(*),BN(NDI,*),CN(NDI,*)
  16. C
  17.       CALL ZERO(AN,NSTRS,1)
  18.       CALL ZERO(BN,NDI,NSTRS)
  19.       CALL ZERO(CN,NDI,NSTRS)
  20. C
  21. C RECUPERATION DES VALEURS DES PARAMETRES MATERIAU
  22. C
  23.       RRDV = CMAT(8 + (2 * NFLU) )
  24.       TDT  = CMAT(NMAT)
  25.       TDEV = CMAT(9 + (2 * NFLU))
  26.       TDEXP = EXP(-TDT/TDEV)
  27.       XNU  = CMAT(24)
  28.  
  29. C
  30. C RECUPERATION DES VARIABLES INTERNES INITIALES
  31. C
  32. C    REMPLISSAGE DES VECTEURS DE DEFORMATION DE FLUAGE
  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
  42. C CONSTRUCTION DE LA MATRICE
  43. C
  44.       ADR = TDEXP - 1.D0
  45.       BDR = 1.D0/RRDV*( TDEXP*(-(2*TDEV/TDT+1.D0) + TDEV/TDT)
  46.      $      + (2.D0*(TDEV-TDT)/TDT+1.D0) - (TDEV-TDT)/TDT )
  47.       CDR = 1.D0/(TDT*RRDV)*( TDEV*TDEXP - (TDEV - TDT) )
  48.  
  49. C
  50.       DO 10 I=1,NSTRS
  51.           AN(I) = ADR * EPSFLU(I)
  52.   10  CONTINUE
  53. C
  54.       IF(IFOUR.EQ.-2) THEN
  55. C
  56. C  - CONTRAINTES PLANES
  57. C
  58.          BN(1,1)=1.D0
  59.          BN(1,2)=-XNU
  60.          BN(2,1)=BN(1,2)
  61.          BN(2,2)=BN(1,1)
  62.          BN(3,3)=1.D0+XNU
  63.       ELSE IF((IFOUR.EQ.-1).OR.(IFOUR.EQ.0)) THEN
  64. C
  65. C  - DEFORMATION PLANE OU AXISYMETRIQUE
  66. C
  67.          BN(1,1)=1.D0
  68.          BN(1,2)=-XNU
  69.          BN(1,3)=BN(1,2)
  70.          BN(2,1)=BN(1,2)
  71.          BN(2,2)=BN(1,1)
  72.          BN(2,3)=BN(1,2)
  73.          BN(3,1)=BN(1,2)
  74.          BN(3,2)=BN(1,2)
  75.          BN(3,3)=BN(1,1)
  76.          BN(4,4)=1.D0+XNU
  77.       ELSE
  78. C
  79. C  - CALCUL TRIDIMENSIONEL
  80. C
  81.          BN(1,1)=1.D0
  82.          BN(1,2)=-XNU
  83.          BN(1,3)=BN(1,2)
  84.          BN(2,1)=BN(1,2)
  85.          BN(2,2)=BN(1,1)
  86.          BN(2,3)=BN(1,2)
  87.          BN(3,1)=BN(1,2)
  88.          BN(3,2)=BN(1,2)
  89.          BN(3,3)=BN(1,1)
  90.          BN(4,4)=1.D0+XNU
  91.          BN(5,5)=BN(4,4)
  92.          BN(6,6)=BN(4,4)
  93.       ENDIF
  94. C
  95.       DO 20 I=1,NSTRS
  96.         DO 20 J=1,NSTRS
  97.           CN(I,J) = CDR * BN(I,J)
  98.           BN(I,J) = BDR * BN(I,J)
  99.   20  CONTINUE
  100. C
  101.       END
  102.  
  103.  
  104.  
  105.  
  106.  

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