Télécharger la plaquette Cast3M


fr en

Télécharger flude1.eso

Retour à la liste

Numérotation des lignes : 

flude1

  1. C FLUDE1    SOURCE    CB215821  16/04/21    21:16:54     8920
  2.       SUBROUTINE FLUDE1(VARF,NVARI,CMAT,NMAT,IFOUR,NSTRS,
  3.      $                       AN,BN,CN,NDI)
  4. C
  5. C_______________________________________________________________________
  6. C
  7. C SOUS PROGRAMME QUI CALCUL LES MATRICES DE DEFORMATION
  8. C  DE FLUAGE DESSICCATION
  9. C   D APRES LE MODELE BAZANT
  10. C
  11. C
  12. C_______________________________________________________________________
  13. C
  14.       IMPLICIT INTEGER(I-N)
  15.       IMPLICIT REAL*8(A-H,O-Z)
  16.       DIMENSION VARF(NVARI),CMAT(NMAT)
  17.       DIMENSION AN(*),BN(NDI,*),CN(NDI,*)
  18. C
  19.       CALL ZERO(AN,NSTRS,1)
  20.       CALL ZERO(BN,NDI,NSTRS)
  21.       CALL ZERO(CN,NDI,NSTRS)
  22. C
  23. C RECUPERATION DES VALEURS DES PARAMETRES MATERIAU
  24. C
  25.       XFDE = CMAT(25)
  26.       XNU  = CMAT(24)
  27. C
  28. C RECUPERATION DES TEMPERATURES
  29. C
  30.       TINI = CMAT(29)
  31.       TFIN = CMAT(30)
  32.       HINI = EXP(0.018*TINI/(1000.*8.314*293.))
  33.       HFIN = EXP(0.018*TFIN/(1000.*8.314*293.))
  34. C
  35. C CONSTRUCTION DE LA MATRICE
  36. C
  37.       ADR = 0.D0
  38.       BDR = XFDE*( ABS(HFIN - HINI) )/2.D0
  39.       CDR = BDR
  40. C
  41.       DO 10 I=1,NSTRS
  42.           AN(I) = 0.D0
  43.   10  CONTINUE
  44. C
  45.       IF(IFOUR.EQ.-2) THEN
  46. C
  47. C  - CONTRAINTES PLANES
  48. C
  49.          BN(1,1)=1.D0
  50.          BN(1,2)=-XNU
  51.          BN(2,1)=BN(1,2)
  52.          BN(2,2)=BN(1,1)
  53.          BN(3,3)=1.D0+XNU
  54.       ELSE IF((IFOUR.EQ.-1).OR.(IFOUR.EQ.0)) THEN
  55. C
  56. C  - DEFORMATION PLANE OU AXISYMETRIQUE
  57. C
  58.          BN(1,1)=1.D0
  59.          BN(1,2)=-XNU
  60.          BN(1,3)=BN(1,2)
  61.          BN(2,1)=BN(1,2)
  62.          BN(2,2)=BN(1,1)
  63.          BN(2,3)=BN(1,2)
  64.          BN(3,1)=BN(1,2)
  65.          BN(3,2)=BN(1,2)
  66.          BN(3,3)=BN(1,1)
  67.          BN(4,4)=1.D0+XNU
  68.       ELSE
  69. C
  70. C  - CALCUL TRIDIMENSIONEL
  71. C
  72.          BN(1,1)=1.D0
  73.          BN(1,2)=-XNU
  74.          BN(1,3)=BN(1,2)
  75.          BN(2,1)=BN(1,2)
  76.          BN(2,2)=BN(1,1)
  77.          BN(2,3)=BN(1,2)
  78.          BN(3,1)=BN(1,2)
  79.          BN(3,2)=BN(1,2)
  80.          BN(3,3)=BN(1,1)
  81.          BN(4,4)=1.D0+XNU
  82.          BN(5,5)=BN(4,4)
  83.          BN(6,6)=BN(4,4)
  84.       ENDIF
  85. C
  86.       DO 20 I=1,NSTRS
  87.         DO 20 J=1,NSTRS
  88.           CN(I,J) = CDR * BN(I,J)
  89.           BN(I,J) = BDR * BN(I,J)
  90.   20  CONTINUE
  91. C
  92.       END
  93.  
  94.  
  95.  
  96.  
  97.  

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