Numérotation des lignes :

C ELASF1    SOURCE    AM        09/12/08    21:15:12     6582      SUBROUTINE ELASF1(DEP,AN,BN,CN,DEPST,SIGEF,     \$                   NSTRS,DSIGT,KERRE)CC ROUTINE APPELE DANS CAMECAC_______________________________________________________________________CC SOUS PROGRAMME QUI CALCUL LES CONTRAINTES DE TEST FINALESC  ELASTIQUES CORRIGE PAR LE FLUAGE TOTALC_______________________________________________________________________C      IMPLICIT INTEGER(I-N)      IMPLICIT REAL*8(A-H,O-Z)      DIMENSION AN(NSTRS),BN(NSTRS,NSTRS),CN(NSTRS,NSTRS)      DIMENSION DEP(NSTRS,NSTRS),SIGEF(NSTRS),DSIGT(NSTRS)      DIMENSION DEPST(NSTRS)*      DIMENSION EFLU(10,10),TEMP(10,10)      DIMENSION DEFLUN(10),DEFELA(10) CC INITIALISATION DES VARIABLESC      NDI = 10      CALL ZERO(DSIGT,NSTRS,1)      CALL ZERO(EFLU,NDI,NSTRS)      CALL ZERO(TEMP,NDI,NSTRS)      CALL ZERO(DEFELA,NSTRS,1)      CALL ZERO(DEFLUN,NSTRS,1)C      DO 5 I=1,NSTRS        DO 5 J=1,NSTRS          DO 5 K=1,NSTRS            TEMP(I,J) = TEMP(I,J) + DEP(I,K) * CN(K,J)  5   CONTINUECC CONSTRUCTION DE LA MATRICE D ELASTICITE CORRIGE PAR LE FLUAGECC  EQUATION : (2.4-1)C        DO 10 I=1,NSTRS          TEMP(I,I) = 1.D0 + TEMP(I,I)  10    CONTINUEC        CALL INVMA1(TEMP,NDI,NSTRS,ISING)        IF (ISING.EQ.1) THEN          WRITE(*,*) ' '          WRITE(*,*) 'ERREUR DANS LE CALCUL DE LA MATRICE DE HOOKE'          WRITE(*,*) 'CORRIGE PAR LE FLUAGE DANS ELASF1'          WRITE(*,*) 'MATRICE NON INVERSIBLE'          WRITE(*,*) ' '          KERRE=99          RETURN        ENDIFCC CALCUL DE (1 + CN)-1 * E0C        DO 20 I=1,NSTRS          DO 20 J=1,NSTRS            DO 20 K=1,NSTRS              EFLU(I,J) = EFLU(I,J) + TEMP(I,K)*DEP(K,J)  20    CONTINUECC        DO 30 I=1,NSTRSC          DO 30 J=1,NSTRSC            DEP(I,J) = EFLU(I,J)C  30    CONTINUECC CALCUL DES DEFORMATIONS DE FLUAGE PUIS ELASTIQUESC      DO 40 I=1,NSTRS        DEFLUN(I) = DEFLUN(I) + AN(I)        DO 40 J=1,NSTRS          DEFLUN(I) = DEFLUN(I) + (BN(I,J)+CN(I,J)) * SIGEF(J)  40  CONTINUEC      DO 50 I=1,NSTRS        DEFELA(I) = DEPST(I) - DEFLUN(I)  50  CONTINUECC CALCUL DE L INCREMENT DE CONTRAINTES ELASTIQUESCC  EQUATION (2.4-2)C      DO 60 I=1,NSTRS        DO 60 J=1,NSTRS          DSIGT(I) = DSIGT(I) + EFLU(I,J) * DEFELA(J)  60  CONTINUEC      END

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