Télécharger masba3.eso

Retour à la liste

Numérotation des lignes :

  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. -INC CCOPTIO
  19. DIMENSION REL(LRE,LRE),TGL(3,LRE),XE(3,3)
  20. dimension xma(3,3),xga(3),pga(3),pn(3,3),dn(3,3),xj(3)
  21. pga(1)=0.5555555555555555D0 *ros
  22. pga(3)=pga(1)
  23. pga(2)=0.8888888888888888d0*ros
  24. xga(1)=-0.774596669241483D0
  25. xga(3)=-xga(1)
  26. xga(2)=0.d0
  27. do 100 i=1,3
  28. xx = xga(i)
  29. pn(i,1)= 0.5d0 * xx * ( xx -1.d0)
  30. pn(i,2)= -(xx -1.d0)*(xx+1.d0)
  31. pn(i,3)= 0.5d0 *xx*(xx+1.d0)
  32. dn(i,1)=xx-0.5d0
  33. dn(i,2)= -2.d0* xx
  34. dn(i,3)=xx +0.5d0
  35. 100 continue
  36. do 101 i=1,3
  37. dx = dn(i,1)*xe(1,1)+dn(i,2)*xe(1,2)+dn(i,3)*xe(1,3)
  38. dy = dn(i,1)*xe(2,1)+dn(i,2)*xe(2,2)+dn(i,3)*xe(2,3)
  39. dz = dn(i,1)*xe(3,1)+dn(i,2)*xe(3,2)+dn(i,3)*xe(3,3)
  40. ds2 = dx*dx + dy * dy + dz*dz
  41. ds = sqrt(ds2)
  42. xj(i)=ds
  43. 101 continue
  44. do 103 i=1,3
  45. do 103 j=1,3
  46. 103 xma(i,j)=0.d0
  47. do 102 ipg=1,3
  48. xjpg=xj(ipg)*pga(ipg)
  49. do 102 i=1,3
  50. do 102 j=1,3
  51. xma(i,j)=xma(i,j)+ xjpg*pn(ipg,i)*pn(ipg,j)
  52. 102 continue
  53. C
  54. KERRE=0
  55. C
  56. do 60 I=1,3
  57. i1=(i-1)*3
  58. do 60 j=1,3
  59. j1=(j-1)*3
  60. do 60 k=1,3
  61. ii= i1 +k
  62. jj = j1 + k
  63. rel(ii,jj)=xma(i,j)
  64. 60 continue
  65. RETURN
  66. 999 KERRE=1
  67. RETURN
  68. END
  69.  
  70.  
  71.  

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