Télécharger la plaquette Cast3M


fr en

Télécharger defcr1.eso

Retour à la liste

Numérotation des lignes : 

defcr1

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

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