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    CB215821  26/08/24    21:15:16     12622          
  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 667 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.   667 CONTINUE
  41.       GOTO  666
  42. C
  43. C     ELEMENTS LIQUIDES BIDIM FOURIER
  44. C
  45.   30  CONTINUE
  46.       DO 668 IX=1,LRE
  47.       DO 300 IY=1,IX
  48.       REL(IX,IY)=REL(IX,IY)+DJAC*VML22*(BLX(2,IY)*BLX(2,IX)+
  49.      1                        BLY(2,IY)*BLY(2,IX)+
  50.      2                        BLT(2,IY)*BLT(2,IX))
  51.       REL(IY,IX)=REL(IX,IY)
  52. 300   CONTINUE
  53.   668 CONTINUE
  54.       GOTO  666
  55. C
  56. C     ELEMENTS LIQUIDES TRIDIM
  57. C
  58.   40  CONTINUE
  59.       DO 669 IX=1,LRE
  60.       DO 400 IY=1,IX
  61.       REL(IX,IY)=REL(IX,IY)+DJAC*VML22*(BLX(2,IY)*BLX(2,IX)+
  62.      1                        BLY(2,IY)*BLY(2,IX)+
  63.      2                        BLZ(2,IY)*BLZ(2,IX))
  64.       REL(IY,IX)=REL(IX,IY)
  65. 400   CONTINUE
  66.   669 CONTINUE
  67.       GOTO  666
  68. C
  69.  666  CONTINUE
  70.       RETURN
  71.       END
  72.  
  73.  
  74.  

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