Télécharger masba3.eso

Retour à la liste

Numérotation des lignes :

masba3
  1. C MASBA3 SOURCE CHAT 05/01/13 01:30:56 5004
  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 103 i=1,3
  47. do 103 j=1,3
  48. 103 xma(i,j)=0.d0
  49. do 102 ipg=1,3
  50. xjpg=xj(ipg)*pga(ipg)
  51. do 102 i=1,3
  52. do 102 j=1,3
  53. xma(i,j)=xma(i,j)+ xjpg*pn(ipg,i)*pn(ipg,j)
  54. 102 continue
  55. C
  56. KERRE=0
  57. C
  58. do 60 I=1,3
  59. i1=(i-1)*3
  60. do 60 j=1,3
  61. j1=(j-1)*3
  62. do 60 k=1,3
  63. ii= i1 +k
  64. jj = j1 + k
  65. rel(ii,jj)=xma(i,j)
  66. 60 continue
  67. RETURN
  68. 999 KERRE=1
  69. RETURN
  70. END
  71.  
  72.  
  73.  

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