Télécharger la plaquette Cast3M


fr en

Télécharger defcr3.eso

Retour à la liste

Numérotation des lignes : 

defcr3

  1. C DEFCR3    SOURCE    FD218221  26/02/13    21:15:11     12462          
  2.       SUBROUTINE DEFCR3(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 DE DESSICATION D APRES LE MODELE 
  9. C  ADOPTE PAR F. BENBOUDJEMA (similaire au modele de Bazant) 
  10. C_______________________________________________________________________
  11. C
  12.       IMPLICIT INTEGER(I-N)
  13.       IMPLICIT REAL*8(A-H,O-Z)
  14.       DIMENSION VARF(NVARI),CMAT(NMAT)
  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.       TDT  = CMAT(28)
  24.       XFDE = CMAT(15)
  25.       XNU  = CMAT(14)
  26. C
  27. C RECUPERATION DU INCREMENT PRESSION CAPILLAIRE (DPC) DANS LE PAS
  28. C
  29.       DPC = CMAT(17)
  30. C
  31. C CONSTRUCTION DE LA MATRICE
  32. C
  33.       ADR = 0.D0
  34.       BDR = (XFDE * DPC * TDT)/2.D0
  35.       CDR = BDR
  36. C
  37.       DO 10 I=1,NSTRS
  38.           AN(I) = 0.D0
  39.   10  CONTINUE
  40. C
  41.       IF(IFOUR.EQ.-2) THEN
  42. C
  43. C  - CONTRAINTES PLANES 
  44. C
  45.          BN(1,1)=1.D0
  46.          BN(1,2)=-XNU
  47.          BN(2,1)=BN(1,2)
  48.          BN(2,2)=BN(1,1)
  49.          BN(3,3)=2.D0*(1.D0+XNU)
  50.       ELSE IF((IFOUR.EQ.-1).OR.(IFOUR.EQ.0)) THEN
  51. C
  52. C  - DEFORMATION PLANE OU AXISYMETRIQUE 
  53. C
  54.          BN(1,1)=1.D0
  55.          BN(1,2)=-XNU
  56.          BN(1,3)=BN(1,2)
  57.          BN(2,1)=BN(1,2)
  58.          BN(2,2)=BN(1,1)
  59.          BN(2,3)=BN(1,2)
  60.          BN(3,1)=BN(1,2)
  61.          BN(3,2)=BN(1,2)
  62.          BN(3,3)=BN(1,1)
  63.          BN(4,4)=2.D0*(1.D0+XNU)
  64.       ELSE
  65. C
  66. C  - CALCUL TRIDIMENSIONEL 
  67. C
  68.          BN(1,1)=1.D0
  69.          BN(1,2)=-XNU
  70.          BN(1,3)=BN(1,2)
  71.          BN(2,1)=BN(1,2)
  72.          BN(2,2)=BN(1,1)
  73.          BN(2,3)=BN(1,2)
  74.          BN(3,1)=BN(1,2)
  75.          BN(3,2)=BN(1,2)
  76.          BN(3,3)=BN(1,1)
  77.          BN(4,4)=2.D0*(1.D0+XNU)
  78.          BN(5,5)=BN(4,4)
  79.          BN(6,6)=BN(4,4)
  80.       ENDIF
  81. C
  82.       DO 20 I=1,NSTRS
  83.         DO 21 J=1,NSTRS
  84.           CN(I,J) = CDR * BN(I,J)
  85.           BN(I,J) = BDR * BN(I,J)
  86.   21    CONTINUE
  87.   20  CONTINUE
  88. C
  89.       END
  90.  
  91.  

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