Sources Esope | Cast3M

flutra
C FLUTRA    SOURCE    CB215821  16/04/21    21:16:55     8920
      SUBROUTINE FLUTRA(SIG0,NSTRS,DEPST,VAR0,NVARI,XMAT,NMATT,
     *            IFOUR,TDT,IB,IGAU,TETA1,TETA2,ITHER,NMAT,
     *            SIGF,VARF,CMAT,DEPSI,SIGI,DEP,EPSILO,EPSEQU,
     *            EPSFLU,DSIGT,SIGEF,AN,BN,CN,KERRE)
C
C___________________________________________________________________
C
C             MODELE N. REVIRON
C
C DEFINITION DES DIFFERENTES VARIABLES ENTREES, SORTIES,
C PARAMETRES
C
C___________________________________________________________________
C
C     NSTRS       = TAILLE DU VECTEUR DES CONTRAINTES
C
C ENTREES 3D
C
C     SIG0(NSTRS) = VECTEUR DES CONTRAINTES INITIALES APPARENTES
C                          1 : XX
C                          2 : YY
C                          3 : ZZ
C                          4 : XY
C                          5 : ZX
C                          6 : YZ
C
C     SIGI(NSTRS) = VECTEUR DES CONTRAINTES INITIALES EFFECTIVES
C     (SE11,SE22,SE12)
C                          idem
C
C     DEPST(NSTRS)= VECTEUR DES INCREMENTS DE DEFORMATIONS TOTALES
C                          idem
C
C ENTREES 2D CONTRAINTES PLANES
C
C     SIG0(NSTRS) = VECTEUR DES CONTRAINTES INITIALES APPARENTES
C                          1 : XX
C                          2 : YY
C                         (3 : ZZ)
C                          4 : XY
C
C     SIGI(NSTRS) = VECTEUR DES CONTRAINTES INITIALES EFFECTIVES
C     (SE11,SE22,SE12)
C                          idem
C
C     DEPST(NSTRS)= VECTEUR DES INCREMENTS DE DEFORMATIONS TOTALES
C                          1 : XX
C                          2 : YY
C                         (3 : ZZ)
C                          4 : XY
C
C     DEPSI (UTILSER DANS LES SUBROUTINES)
C                          1 : XX
C                          2 : YY
C                          3 : XY
C                         (4 : ZZ)
C
C ENTREES 2D DEFORMATIONS PLANES
C
C     SIG0(NSTRS) = VECTEUR DES CONTRAINTES INITIALES APPARENTES
C                          1 : XX
C                          2 : YY
C                          3 : ZZ
C                          4 : XY
C
C     SIGI(NSTRS) = VECTEUR DES CONTRAINTES INITIALES EFFECTIVES
C     (SE11,SE22,SE12)
C                          idem
C
C     DEPST(NSTRS)= VECTEUR DES INCREMENTS DE DEFORMATIONS TOTALES
C                          1 : XX
C                          2 : YY
C                         (3 : ZZ = 0)
C                          4 : XY
C
C ENTREES AXISYMETRIE
C
C     SIG0(NSTRS) = VECTEUR DES CONTRAINTES INITIALES APPARENTES
C                          1 : RR
C                          2 : ZZ
C                          3 : TT
C                          4 : RZ
C
C     SIGI(NSTRS) = VECTEUR DES CONTRAINTES INITIALES EFFECTIVES
C                          idem
C
C     DEPST(NSTRS)= VECTEUR DES INCREMENTS DE DEFORMATIONS TOTALES
C                          idem
C___________________________________________________________
C
C     VAR0(NVARI) = VECTEUR DES VARIABLES INTERNES INITIALES
C___________________________________________________________
C
C       VARIABLES INTERNES MODELE MECANIQUE ISOTROPE
C
C     VAR0(1)    = EPSTIM
C     VAR0(2)    = DINI
C     VAR0(3)    = EF11
C     VAR0(4)    = EF22
C     VAR0(5)    = EF33
C     VAR0(6)    = EF12
C     VAR0(7)    = EF31
C     VAR0(8)    = EF23
C     VAR0(9)    = SE11
C     VAR0(10)   = SE22
C     VAR0(11)   = SE12
C__________________________________________________________
C
C     NVARI       = TAILLE DU VECTEUR DE VARIABLES INTERNES
C__________________________________________________________
C
C     XMAT(NVARI) = VECTEUR DES PARAMETRES DU MATERIAU
C
C       PARAMETRES STANDARD
C     XMAT(1)     = YOUN : MODULE D YOUNG
C     XMAT(2)     = XNU  : COEFFICIENT DE POISSON ELASTIQUE
C     XMAT(3)     = RHO  : MASSE VOLUMIQUE
C     XMAT(4)     = ALPH : COEFFICIENT DE DILATATION THERMIQUE
C
C       PARAMETRES MODELE MECANIQUE
C     XMAT(5)     = KTR0 : KTR0 Mazars
C     XMAT(6)    = ATRA : ATRA Mazars
C     XMAT(7)    = BTRA : BTRA Mazars
C     XMAT(11)    = A : COEF De Schutter Module d YOUNG
C     XMAT(12)    = B : COEF De Schutter fc : non utilise
C     XMAT(13)    = C : COEF De Schutter ft
C     XMAT(14)    = D : COEF De Schutter Gft
C___________________________________________________
C
C     NMATT       = TAILLE DU VECTEUR DES PARAMETRES
C     MELE        = NUMERO DE L ELEMENT FINI UTILISE
C     XE          = VECTEUR DES COORDONNEES DE L ELEMENT
C     NBNN        = TAILLE DU VECTEUR XE
C     IFOUR       = INDICATEUR DU TYPE DE CALCUL
C                     -3 : DEFORMATION PLANE GENERALISEE
C                     -2 : CONTRAINTES PLANES
C                     -1 : DEFORMATION PLANE
C                      0 : AXISYMETRIQUE
C                      1 : SERIE DE FOURIER
C                      2 : TRIDIMENSIONNEL
C     MFR         = NUMERO DE FORMULATION
C                      1 : ELEMENTS MASSIFS
C                   AUTRE: VOIR nummfr.eso
C     ITTHER      = INDICATEUR DE CHARGEMENT THERMIQUE
C                      0 : PAS DE CHARGEMENT THERMIQUE
C                      1 : PRESENCE DE CHARGEMENT THERMIQUE
C     TINI        = TEMPERATURE EN DEBUT DE PAS
C     TFIN        = TEMPERATURE EN FIN DE PAS
C     TDT         = PAS DE TEMPS
C     IB          = NUMERO DE L ELEMENT CALCULE
C     IGAU        = NUMERO DU POINT DE GAUSS
C     DEP(6,6)    = MATRICE ELASTO-PLASTIQUE TANGENTE
C     DINI        = ENDOMMAGEMENT EN DEBUT DE PAS
C     D           = ENDOMMAGEMENT EN FIN DE PAS
C     SIG0        = VECTEUR DES CONTRAINTES INITIALES APPARENTES
C     SIGF        = VECTEUR DES CONTRAINTES FINALES APPARENTES
C     SIGI        = VECTEUR DES CONTRAINTES INITIALES EFFECTIVES
C     EPSILO      = VECTEUR DES DEFORMATIONS ELASTIQUES FIN PAS
C____________________________________________________
C
C  VARIABLES EN SORTIES
C____________________________________________________
C
C     SIGF(NSTRS) = VECTEUR DES CONTRAINTES FINALES
C REMARQUE :
C SIGF est dans un premier temps les contraintes effectives finales
C puis les contraintes apparentes finales par multiplication par D
C
C     DEPSFT      = VECTEUR INCREMENT DES DEFORMATIONS DE FLUAGE
C                   (EF11,EF22,EF33,EF12,EF31,EF23)
C     VARF(NVARI) = VECTEUR DES VARIABLES INTERNES FINALES
C     KERRE       = INDICATEUR D ERREUR
C                      0 : PAS D ERREUR
C
C____________________________________________________
C____________________________________________________
C____________________________________________________
C____________________________________________________
C
C
      IMPLICIT INTEGER(I-N)
      IMPLICIT REAL*8(A-H,O-Z)
      DIMENSION XMAT(NMATT),CMAT(NMAT)
      DIMENSION VAR0(NVARI),VARF(NVARI)
      DIMENSION DEPST(NSTRS),DEPSI(NSTRS)
      DIMENSION SIG0(NSTRS),SIGF(NSTRS),SIGI(NSTRS)
      DIMENSION DEP(NSTRS,NSTRS)
      DIMENSION EPSILO(NSTRS)
      DIMENSION EPSEQU(NSTRS),EPSFLU(NSTRS)
      DIMENSION DSIGT(NSTRS),SIGEF(NSTRS)
      DIMENSION AN(NSTRS),BN(NSTRS,NSTRS),CN(NSTRS,NSTRS)
C
C TEST SUR L OPTION DE CALCUL
C  LE CALCUL N EST POSSIBLE QUE DANS LE CAS
C   CONTRAINTES PLANES : IFOUR = -2
C   DEFORMATION PLANE  : IFOUR = -1
C   AXISYMETRIQUE      : IFOUR =  0
C   TRIDIMENSIONEL     : IFOUR =  2
C
      IF ((IFOUR.NE.-2).AND.(IFOUR.NE.-1).AND.(IFOUR.NE.0)
     $    .AND.(IFOUR.NE.2)) THEN
        WRITE(*,*) 'L OPTION DE CALCUL N EST PAS IMPLEMENTE'
        WRITE(*,*) 'IL FAUT CHOISIR PARMIS LES OPTIONS SUIVANTES :'
        WRITE(*,*) '  CONTRAINTES PLANES'
        WRITE(*,*) '  DEFORMATIONS PLANES'
        WRITE(*,*) '  AXISYMETRIQUE'
        WRITE(*,*) '  TRIDIMENSIONEL'
        KERRE = 99
        RETURN
      ENDIF
C
C INITIALISATION DES MATRICES
C
      XPRECIS = 1.D-6
      CALL ZERO(CMAT,NMAT,1)
      CALL ZERO(SIGF,NSTRS,1)
C
C RECUPERATION DES VALEURS DES PARAMETRES MATERIAU
C
C
C REMPLISSAGE D UN NOUVEAU VECTEUR DE MATERIAU
C
      DO 300 I =1, NMATT
        CMAT(I)= XMAT(I)
  300 CONTINUE
C
C   AJOUT DU TEMPS A LA FIN
C
      CMAT(NMAT) = TDT
C
C   GESTION DES TEMPERATURES
C
      CMAT(29) = TETA1
      CMAT(30) = TETA2
C
C RECUPERATION DES VARIABLES INTERNES INITIALES
C
C REMPLISSAGE DES VARIABLES INTERNES FINALES
C
      DO 200 I =1, NVARI
        VARF(I)= VAR0(I)
  200 CONTINUE
C
C_______________________________________________________________________
C
C REDIMENSIONNEMENT DES VECTEURS DANS LE CAS DE CONTRAINTES PLANES
C_______________________________________________________________________
C
C
C
C CHANGEMENT DE NOM POUR TRAITEMENT
C
      DO 100 I =1, NSTRS
        DEPSI(I)= DEPST(I)
  100 CONTINUE
C
      IPLAN = 0
      IF ((NSTRS.EQ.4).AND.(IFOUR.EQ.-2)) THEN
         IPLAN = 1
         DEPSI(3) = DEPST(4)
         DEPSI(4) = DEPST(3)
         IF (ABS(SIG0(3)).GE.1.D-3) THEN
           WRITE(*,*)'PROBLEME AVEC LES CONTRAINTES PLANES'
           KERRE=99
           RETURN
         ENDIF
      ENDIF
C_______________________________________________________________________
C
C MODELE MECANIQUE
C_______________________________________________________________________
C
       CALL CAMECA(NSTRS,DEPSI,VAR0,NVARI,CMAT,NMAT,
     *            IFOUR,IB,IGAU,SIGF,VARF,SIGI,DEP,
     *            EPSILO,EPSEQU,EPSFLU,DSIGT,SIGEF,
     *            AN,BN,CN,KERRE)
C
C
C_______________________________________________________________________
C
C REDIMENSIONNEMENT DES VECTEURS DANS LE CAS DE CONTRAINTES PLANES
C_______________________________________________________________________
C
      IF (IPLAN.EQ.1) THEN
         temp = SIGF(4)
         SIGF(4)=SIGF(3)
         SIGF(3)=temp
      ENDIF
C
      END