Télécharger la plaquette Cast3M


fr en

Télécharger defcr2.eso

Retour à la liste

Numérotation des lignes : 

defcr2

  1. C DEFCR2    SOURCE    FD218221  26/02/13    21:15:10     12462          
  2.       SUBROUTINE DEFCR2(VARF,NVARI,CMAT,NMAT,IFOUR,NSTRS,
  3.      $                       AN,BN,CN,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 AMORTISSEUR
  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)
  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.       VIDV = CMAT(13)
  23.       TDT  = CMAT(28)
  24.       XNU  = CMAT(14)
  25. C
  26. C RECUPERATION DES VARIABLES INTERNES INITIALES
  27. C
  28. C    REMPLISSAGE DES VECTEURS DE DEFORMATION DE FLUAGE
  29. C
  30.       EPSFLU(1) = VARF(15)
  31.       EPSFLU(2) = VARF(16)
  32.       EPSFLU(3) = VARF(17)
  33.       EPSFLU(4) = VARF(18)
  34.       EPSFLU(5) = VARF(19)
  35.       EPSFLU(6) = VARF(20)
  36. C
  37. C
  38. C CONSTRUCTION DE LA MATRICE 
  39. C
  40.       ADR = 0.D0
  41.       BDR = TDT/(2.D0*VIDV)
  42.       CDR = TDT/(2.D0*VIDV)
  43. C
  44.       DO 10 I=1,NSTRS
  45.           AN(I) = ADR * EPSFLU(I)
  46.   10  CONTINUE
  47. C
  48.       IF(IFOUR.EQ.-2) THEN
  49. C
  50. C  - CONTRAINTES PLANES 
  51. C
  52.          BN(1,1)=1.D0
  53.          BN(1,2)=-XNU
  54.          BN(2,1)=BN(1,2)
  55.          BN(2,2)=BN(1,1)
  56.          BN(3,3)=2.D0*(1.D0+XNU)
  57.       ELSE IF((IFOUR.EQ.-1).OR.(IFOUR.EQ.0)) THEN
  58. C
  59. C  - DEFORMATION PLANE OU AXISYMETRIQUE 
  60. C
  61.          BN(1,1)=1.D0
  62.          BN(1,2)=-XNU
  63.          BN(1,3)=BN(1,2)
  64.          BN(2,1)=BN(1,2)
  65.          BN(2,2)=BN(1,1)
  66.          BN(2,3)=BN(1,2)
  67.          BN(3,1)=BN(1,2)
  68.          BN(3,2)=BN(1,2)
  69.          BN(3,3)=BN(1,1)
  70.          BN(4,4)=2.D0*(1.D0+XNU)
  71.       ELSE
  72. C
  73. C  - CALCUL TRIDIMENSIONEL 
  74. C
  75.          BN(1,1)=1.D0
  76.          BN(1,2)=-XNU
  77.          BN(1,3)=BN(1,2)
  78.          BN(2,1)=BN(1,2)
  79.          BN(2,2)=BN(1,1)
  80.          BN(2,3)=BN(1,2)
  81.          BN(3,1)=BN(1,2)
  82.          BN(3,2)=BN(1,2)
  83.          BN(3,3)=BN(1,1)
  84.          BN(4,4)=2.D0*(1.D0+XNU)
  85.          BN(5,5)=BN(4,4)
  86.          BN(6,6)=BN(4,4)
  87.       ENDIF
  88. C
  89.       DO 20 I=1,NSTRS
  90.         DO 21 J=1,NSTRS
  91.           CN(I,J) = CDR * BN(I,J)
  92.           BN(I,J) = BDR * BN(I,J)
  93.   21    CONTINUE
  94.   20  CONTINUE
  95. C
  96.       END
  97.  
  98.  

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