Télécharger la plaquette Cast3M


fr en

Télécharger blmast.eso

Retour à la liste

Numérotation des lignes : 

blmast

  1. C BLMAST    SOURCE    CHAT      05/01/12    21:40:59     5004
  2.       SUBROUTINE BLMAST(IGAU,MFR,NBNO,LRE,IFOU,NDDL,NN,XEL,
  3.      1                     SHPTOT,SHP,BLX,BLY,BLZ,BLT,DJAC)
  4. C=======================================================================
  5. C
  6. C    CALCULE LES MATRICES BL POUR LES ELEMENTS LIQUIDES
  7. C
  8. C          ROUTINE FORTRAN PUR
  9. C=======================================================================
  10. C  INPUT
  11. C     IGAU=NUMERO DU POINT DE GAUSS
  12. C     MFR =NUMERO DE LA FORMULATION
  13. C     NBNO=NOMBRE DE NOEUDS
  14. C     LRE =NOMBRE DE COLONNES DE LA MATRICE B
  15. C     IFOU=IFOUR DE CCOPTIO
  16. C     NDDL=NOMBRE DE DEGRES DE LIBERTE
  17. C     NN  =NUMERO DU MODE DE FOURIER
  18. C     XEL =COORDONNEES LOCALES DE L ELEMENT
  19. C     SHPTOT(6,NBNO,NBGAU)=FONCTIONS DE FORMES ET DERIVEES
  20. C  ZONE DE TRAVAIL
  21. C     SHP(6,NBNO)=TABLEAU DE TRAVAIL
  22. C  OUTPUT
  23. C     DJAC=JACOBIEN
  24. C     BLX = MATRICE CONTENANT LES DERIVEES DES FONCTIONS DE FORME / X
  25. C     BLY = MATRICE CONTENANT LES DERIVEES DES FONCTIONS DE FORME / Y
  26. C     BLZ = MATRICE CONTENANT LES DERIVEES DES FONCTIONS DE FORME / Z
  27. C     BLT = MATRICE CONTENANT LES DERIVEES DES FONCTIONS DE FORME / TETA
  28. C=======================================================================
  29.       IMPLICIT INTEGER(I-N)
  30.       IMPLICIT REAL*8(A-H,O-Z)
  31.       DIMENSION XEL(3,1),SHP(6,1),SHPTOT(6,NBNO,1)
  32.       DIMENSION BLX(NDDL,1),BLY(NDDL,1),BLZ(NDDL,1),BLT(NDDL,1)
  33.       CALL ZERO(BLX,NDDL,LRE)
  34.       CALL ZERO(BLY,NDDL,LRE)
  35.       CALL ZERO(BLZ,NDDL,LRE)
  36.       CALL ZERO(BLT,NDDL,LRE)
  37.       IFR=IFOU+4
  38.       DIM3=1.D0
  39.       RR=0.D0
  40.       GOTO (666,10,10,20,30,40) ,IFR
  41. C
  42. C     ELEMENTS LIQUIDES BIDIM   CONTRAINTES PLANES OU DEFRMTNS PLANES
  43. C
  44.   10  CONTINUE
  45.       DO 101 NP=1,NBNO
  46.       SHP(2,NP)=SHPTOT(2,NP,IGAU)
  47.       SHP(3,NP)=SHPTOT(3,NP,IGAU)
  48.   101 CONTINUE
  49.       CALL DEVOLU(XEL,SHP,MFR,NBNO,IFOU,NN,2,DIM3,RR,DJAC)
  50.       K=1
  51.       DO 102 NP=1,NBNO
  52.       BLX(2,K+1)=SHP(2,NP)
  53.       BLY(2,K+1)=SHP(3,NP)
  54.   102 K=K+2
  55.       GOTO  666
  56. C
  57. C     ELEMENTS LIQUIDES BIDIM AXISYMETRIQUE
  58. C
  59.   20  CONTINUE
  60.       DO 201 NP=1,NBNO
  61.       SHP(1,NP)=SHPTOT(1,NP,IGAU)
  62.       SHP(2,NP)=SHPTOT(2,NP,IGAU)
  63.       SHP(3,NP)=SHPTOT(3,NP,IGAU)
  64.   201 CONTINUE
  65.       CALL DEVOLU(XEL,SHP,MFR,NBNO,IFOU,NN,2,DIM3,RR,DJAC)
  66.       K=1
  67.       DO 202 NP=1,NBNO
  68.       BLX(2,K+1)=SHP(2,NP)
  69.       BLY(2,K+1)=SHP(3,NP)
  70.   202 K=K+2
  71.       GOTO  666
  72. C
  73. C     ELEMENTS LIQUIDES BIDIM FOURIER
  74. C
  75.   30  CONTINUE
  76.       DO 301 NP=1,NBNO
  77.       SHP(1,NP)=SHPTOT(1,NP,IGAU)
  78.       SHP(2,NP)=SHPTOT(2,NP,IGAU)
  79.       SHP(3,NP)=SHPTOT(3,NP,IGAU)
  80.   301 CONTINUE
  81.       CALL DEVOLU(XEL,SHP,MFR,NBNO,IFOU,NN,2,DIM3,RR,DJAC)
  82.       K=1
  83.       DO 302 NP=1,NBNO
  84.       BLX(2,K+1)= SHP(2,NP)
  85.       BLY(2,K+1)= SHP(3,NP)
  86.       BLT(2,K+1)= -SHP(1,NP)*NN/RR
  87.   302 K=K+2
  88.       GOTO  666
  89. C
  90. C     ELEMENTS LIQUIDES TRIDIM
  91. C
  92.   40  CONTINUE
  93.       DO 401 NP=1,NBNO
  94.       SHP(2,NP)=SHPTOT(2,NP,IGAU)
  95.       SHP(3,NP)=SHPTOT(3,NP,IGAU)
  96.       SHP(4,NP)=SHPTOT(4,NP,IGAU)
  97.   401 CONTINUE
  98.       CALL DEVOLU(XEL,SHP,MFR,NBNO,IFOU,NN,3,DIM3,RR,DJAC)
  99.       K=1
  100.       DO 402 NP=1,NBNO
  101.       BLX(2,K+1)=SHP(2,NP)
  102.       BLY(2,K+1)=SHP(3,NP)
  103.       BLZ(2,K+1)=SHP(4,NP)
  104.   402 K=K+2
  105.       GOTO 666
  106. C
  107.  666  CONTINUE
  108.       RETURN
  109.       END
  110.  
  111.  
  112.  
  113.  

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