Télécharger la plaquette Cast3M


fr en

Télécharger elasf1.eso

Retour à la liste

Numérotation des lignes : 

elasf1

  1. C ELASF1    SOURCE    CB215821  26/08/24    21:16:20     12622          
  2.       SUBROUTINE ELASF1(DEP,AN,BN,CN,DEPST,SIGEF,
  3.      $                   NSTRS,DSIGT,KERRE)
  4. C
  5. C ROUTINE APPELE DANS CAMECA
  6. C_______________________________________________________________________
  7. C
  8. C SOUS PROGRAMME QUI CALCUL LES CONTRAINTES DE TEST FINALES
  9. C  ELASTIQUES CORRIGE PAR LE FLUAGE TOTAL
  10. C_______________________________________________________________________
  11. C
  12.       IMPLICIT INTEGER(I-N)
  13.       IMPLICIT REAL*8(A-H,O-Z)
  14.       DIMENSION AN(NSTRS),BN(NSTRS,NSTRS),CN(NSTRS,NSTRS)
  15.       DIMENSION DEP(NSTRS,NSTRS),SIGEF(NSTRS),DSIGT(NSTRS)
  16.       DIMENSION DEPST(NSTRS)
  17. *
  18.       DIMENSION EFLU(10,10),TEMP(10,10)
  19.       DIMENSION DEFLUN(10),DEFELA(10)
  20.  
  21. C
  22. C INITIALISATION DES VARIABLES
  23. C
  24.       NDI = 10
  25.       CALL ZERO(DSIGT,NSTRS,1)
  26.       CALL ZERO(EFLU,NDI,NSTRS)
  27.       CALL ZERO(TEMP,NDI,NSTRS)
  28.       CALL ZERO(DEFELA,NSTRS,1)
  29.       CALL ZERO(DEFLUN,NSTRS,1)
  30. C
  31.       DO 62 I=1,NSTRS
  32.         DO 61 J=1,NSTRS
  33.           DO 5 K=1,NSTRS
  34.             TEMP(I,J) = TEMP(I,J) + DEP(I,K) * CN(K,J)
  35.   5   CONTINUE
  36.    61   CONTINUE
  37.    62 CONTINUE
  38. C
  39. C CONSTRUCTION DE LA MATRICE D ELASTICITE CORRIGE PAR LE FLUAGE
  40. C
  41. C  EQUATION : (2.4-1)
  42. C
  43.         DO 10 I=1,NSTRS
  44.           TEMP(I,I) = 1.D0 + TEMP(I,I)
  45.   10    CONTINUE
  46. C
  47.         CALL INVMA1(TEMP,NDI,NSTRS,ISING)
  48.         IF (ISING.EQ.1) THEN
  49.           WRITE(*,*) ' '
  50.           WRITE(*,*) 'ERREUR DANS LE CALCUL DE LA MATRICE DE HOOKE'
  51.           WRITE(*,*) 'CORRIGE PAR LE FLUAGE DANS ELASF1'
  52.           WRITE(*,*) 'MATRICE NON INVERSIBLE'
  53.           WRITE(*,*) ' '
  54.           KERRE=99
  55.           RETURN
  56.         ENDIF
  57. C
  58. C CALCUL DE (1 + CN)-1 * E0
  59. C
  60.         DO 64 I=1,NSTRS
  61.           DO 63 J=1,NSTRS
  62.             DO 20 K=1,NSTRS
  63.               EFLU(I,J) = EFLU(I,J) + TEMP(I,K)*DEP(K,J)
  64.   20    CONTINUE
  65.    63     CONTINUE
  66.    64   CONTINUE
  67. C
  68. C        DO 30 I=1,NSTRS
  69. C          DO 30 J=1,NSTRS
  70. C            DEP(I,J) = EFLU(I,J)
  71. C  30    CONTINUE
  72. C
  73. C CALCUL DES DEFORMATIONS DE FLUAGE PUIS ELASTIQUES
  74. C
  75.       DO 65 I=1,NSTRS
  76.         DEFLUN(I) = DEFLUN(I) + AN(I)
  77.         DO 40 J=1,NSTRS
  78.           DEFLUN(I) = DEFLUN(I) + (BN(I,J)+CN(I,J)) * SIGEF(J)
  79.   40  CONTINUE
  80.    65 CONTINUE
  81. C
  82.       DO 50 I=1,NSTRS
  83.         DEFELA(I) = DEPST(I) - DEFLUN(I)
  84.   50  CONTINUE
  85. C
  86. C CALCUL DE L INCREMENT DE CONTRAINTES ELASTIQUES
  87. C
  88. C  EQUATION (2.4-2)
  89. C
  90.       DO 66 I=1,NSTRS
  91.         DO 60 J=1,NSTRS
  92.           DSIGT(I) = DSIGT(I) + EFLU(I,J) * DEFELA(J)
  93.   60  CONTINUE
  94.    66 CONTINUE
  95. C
  96.       END
  97.  
  98.  
  99.  
  100.  

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