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    AM        09/12/08    21:15:12     6582
  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 5 I=1,NSTRS
  32.         DO 5 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. C
  37. C CONSTRUCTION DE LA MATRICE D ELASTICITE CORRIGE PAR LE FLUAGE
  38. C
  39. C  EQUATION : (2.4-1)
  40. C
  41.         DO 10 I=1,NSTRS
  42.           TEMP(I,I) = 1.D0 + TEMP(I,I)
  43.   10    CONTINUE
  44. C
  45.         CALL INVMA1(TEMP,NDI,NSTRS,ISING)
  46.         IF (ISING.EQ.1) THEN
  47.           WRITE(*,*) ' '
  48.           WRITE(*,*) 'ERREUR DANS LE CALCUL DE LA MATRICE DE HOOKE'
  49.           WRITE(*,*) 'CORRIGE PAR LE FLUAGE DANS ELASF1'
  50.           WRITE(*,*) 'MATRICE NON INVERSIBLE'
  51.           WRITE(*,*) ' '
  52.           KERRE=99
  53.           RETURN
  54.         ENDIF
  55. C
  56. C CALCUL DE (1 + CN)-1 * E0
  57. C
  58.         DO 20 I=1,NSTRS
  59.           DO 20 J=1,NSTRS
  60.             DO 20 K=1,NSTRS
  61.               EFLU(I,J) = EFLU(I,J) + TEMP(I,K)*DEP(K,J)
  62.   20    CONTINUE
  63. C
  64. C        DO 30 I=1,NSTRS
  65. C          DO 30 J=1,NSTRS
  66. C            DEP(I,J) = EFLU(I,J)
  67. C  30    CONTINUE
  68. C
  69. C CALCUL DES DEFORMATIONS DE FLUAGE PUIS ELASTIQUES
  70. C
  71.       DO 40 I=1,NSTRS
  72.         DEFLUN(I) = DEFLUN(I) + AN(I)
  73.         DO 40 J=1,NSTRS
  74.           DEFLUN(I) = DEFLUN(I) + (BN(I,J)+CN(I,J)) * SIGEF(J)
  75.   40  CONTINUE
  76. C
  77.       DO 50 I=1,NSTRS
  78.         DEFELA(I) = DEPST(I) - DEFLUN(I)
  79.   50  CONTINUE
  80. C
  81. C CALCUL DE L INCREMENT DE CONTRAINTES ELASTIQUES
  82. C
  83. C  EQUATION (2.4-2)
  84. C
  85.       DO 60 I=1,NSTRS
  86.         DO 60 J=1,NSTRS
  87.           DSIGT(I) = DSIGT(I) + EFLU(I,J) * DEFELA(J)
  88.   60  CONTINUE
  89. C
  90.       END
  91.  
  92.  
  93.  

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