Télécharger la plaquette Cast3M


fr en

Télécharger bmlbst.eso

Retour à la liste

Numérotation des lignes : 

bmlbst

  1. C BMLBST    SOURCE    CHAT      05/01/12    21:41:28     5004
  2. C BMLBST    SOURCE    BROC      85/07/01    23:05:03
  3.       SUBROUTINE BMLBST(BLX,BLY,BLZ,BLT,DJAC,VML22,IFOU,LRE,NDDL,REL)
  4. C=======================================================================
  5. C
  6. C    INTEGRE LE PRODUIT BLT ML BL
  7. C       POUR LES ELEMENTS LIQUIDES
  8. C
  9. C          ROUTINE FORTRAN PUR
  10. C=======================================================================
  11. C  INPUT
  12. C     BLX =MATRICE CONTENANT LES DERIVEES DES FONCTIONS DE FORME /X
  13. C     BLY =MATRICE CONTENANT LES DERIVEES DES FONCTIONS DE FORME /Y
  14. C     BLZ =MATRICE CONTENANT LES DERIVEES DES FONCTIONS DE FORME /Z
  15. C     BLT =MATRICE CONTENANT LES DERIVEES DES FONCTIONS DE FORME /TETA
  16. C     DJAC=JACOBIEN
  17. C     VML22=-(COEFPI**2)/RHO
  18. C     IFOU=IFOUR DE CCOPTIO
  19. C     LRE =NOMBRE DE COLONNES DES MATRICES BL
  20. C     NDDL=NOMBRE DE DEGRES DE LIBERTE
  21. C  OUTPUT
  22. C     REL(LRE,LRE) =MATRICE MASSE
  23. C=======================================================================
  24.       IMPLICIT INTEGER(I-N)
  25.       IMPLICIT REAL*8(A-H,O-Z)
  26.       DIMENSION BLX(NDDL,1),BLY(NDDL,1),BLZ(NDDL,1),BLT(NDDL,1)
  27.       DIMENSION REL(LRE,1)
  28.       IFR=IFOU+4
  29.       GOTO (666,10,10,10,30,40) ,IFR
  30. C
  31. C     ELEMENTS LIQUIDES 2DIM   CONT PLANES OU DEFO PLANES OU AXIS
  32. C
  33.   10  CONTINUE
  34.       DO 100 IX=1,LRE
  35.       DO 100 IY=1,IX
  36.       REL(IX,IY)=REL(IX,IY)+DJAC*VML22*(BLX(2,IY)*BLX(2,IX)+
  37.      1                        BLY(2,IY)*BLY(2,IX))
  38.       REL(IY,IX)=REL(IX,IY)
  39. 100   CONTINUE
  40.       GOTO  666
  41. C
  42. C     ELEMENTS LIQUIDES BIDIM FOURIER
  43. C
  44.   30  CONTINUE
  45.       DO 300 IX=1,LRE
  46.       DO 300 IY=1,IX
  47.       REL(IX,IY)=REL(IX,IY)+DJAC*VML22*(BLX(2,IY)*BLX(2,IX)+
  48.      1                        BLY(2,IY)*BLY(2,IX)+
  49.      2                        BLT(2,IY)*BLT(2,IX))
  50.       REL(IY,IX)=REL(IX,IY)
  51. 300   CONTINUE
  52.       GOTO  666
  53. C
  54. C     ELEMENTS LIQUIDES TRIDIM
  55. C
  56.   40  CONTINUE
  57.       DO 400 IX=1,LRE
  58.       DO 400 IY=1,IX
  59.       REL(IX,IY)=REL(IX,IY)+DJAC*VML22*(BLX(2,IY)*BLX(2,IX)+
  60.      1                        BLY(2,IY)*BLY(2,IX)+
  61.      2                        BLZ(2,IY)*BLZ(2,IX))
  62.       REL(IY,IX)=REL(IX,IY)
  63. 400   CONTINUE
  64.       GOTO  666
  65. C
  66.  666  CONTINUE
  67.       RETURN
  68.       END
  69.  
  70.  

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