Télécharger la plaquette Cast3M


fr en

Télécharger masba3.eso

Retour à la liste

Numérotation des lignes : 

masba3

  1. C MASBA3    SOURCE    CB215821  26/08/24    21:17:15     12622          
  2. C MASBA3    SOURCE    DJER      93/03/02    22:53:00
  3.       SUBROUTINE MASBA3(REL,LRE,ROS,XE,TGL,KERRE)
  4. C-----------------------------------------------------------------------
  5. C  ROUTINE DE CALCUL DE LA MATRICE DE MASSE ELEMENTAIRE DE L ELEMENT DE
  6. C  BARRE A 3 NOEUDS
  7. C
  8. C  ENTREE
  9. C          LRE     DIMENSION DU TABLEAU REL
  10. C          ROS     RHO * SECT
  11. C          XE      COORDONNEES DES 3 NOEUDS
  12. C
  13. C  SORTIE
  14. C          REL     MATRICE DE MASSE ELEMENTAIRE
  15. C-----------------------------------------------------------------------
  16.       IMPLICIT INTEGER(I-N)
  17.       IMPLICIT REAL*8(A-H,O-Z)
  18.  
  19. -INC PPARAM
  20. -INC CCOPTIO
  21.       DIMENSION REL(LRE,LRE),TGL(3,LRE),XE(3,3)
  22.       dimension xma(3,3),xga(3),pga(3),pn(3,3),dn(3,3),xj(3)
  23.       pga(1)=0.5555555555555555D0 *ros
  24.       pga(3)=pga(1)
  25.       pga(2)=0.8888888888888888d0*ros
  26.       xga(1)=-0.774596669241483D0
  27.       xga(3)=-xga(1)
  28.       xga(2)=0.d0
  29.       do 100 i=1,3
  30.       xx = xga(i)
  31.       pn(i,1)= 0.5d0 * xx * ( xx -1.d0)
  32.       pn(i,2)= -(xx -1.d0)*(xx+1.d0)
  33.       pn(i,3)= 0.5d0 *xx*(xx+1.d0)
  34.       dn(i,1)=xx-0.5d0
  35.       dn(i,2)= -2.d0* xx
  36.       dn(i,3)=xx +0.5d0
  37.   100 continue
  38.       do 101 i=1,3
  39.       dx = dn(i,1)*xe(1,1)+dn(i,2)*xe(1,2)+dn(i,3)*xe(1,3)
  40.       dy = dn(i,1)*xe(2,1)+dn(i,2)*xe(2,2)+dn(i,3)*xe(2,3)
  41.       dz = dn(i,1)*xe(3,1)+dn(i,2)*xe(3,2)+dn(i,3)*xe(3,3)
  42.       ds2 = dx*dx + dy * dy + dz*dz
  43.       ds = sqrt(ds2)
  44.       xj(i)=ds
  45.   101 continue
  46.       do 1000 i=1,3
  47.       do 103 j=1,3
  48.       xma(i,j)=0.d0
  49.   103 CONTINUE
  50.  1000 CONTINUE
  51.       do 1002 ipg=1,3
  52.           xjpg=xj(ipg)*pga(ipg)
  53.         do 1001 i=1,3
  54.          do 102 j=1,3
  55.        xma(i,j)=xma(i,j)+ xjpg*pn(ipg,i)*pn(ipg,j)
  56.  102  continue
  57.  1001   CONTINUE
  58.  1002 CONTINUE
  59. C
  60.       KERRE=0
  61. C
  62.       do 1004 I=1,3
  63.        i1=(i-1)*3
  64.       do 1003 j=1,3
  65.          j1=(j-1)*3
  66.        do 60 k=1,3
  67.        ii= i1 +k
  68.        jj = j1 + k
  69.        rel(ii,jj)=xma(i,j)
  70.    60  continue
  71.  1003 CONTINUE
  72.  1004 CONTINUE
  73.       RETURN
  74.  999  KERRE=1
  75.       RETURN
  76.       END
  77.  
  78.  
  79.  
  80.  

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