Télécharger barig3.eso

Retour à la liste

Numérotation des lignes :

  1. C BARIG3 SOURCE CB215821 16/04/21 21:15:13 8920
  2. C BARIG3 SOURCE CHAT 89/05/26 22:52:13
  3. SUBROUTINE BARIG3(REL,LRE,ES,XE,TGL,KERRE,IB)
  4. C-----------------------------------------------------------------------
  5. C ROUTINE DE CALCUL DE LA MATRICE DE RIGIDITE ELEMENTAIRE
  6. C DE L ELEMENT DE BARRE A 3 NOEUDS
  7. C-----------------------------------------------------------------------
  8. C ENTREE
  9. C LRE DIMENSION DU TABLEAU REL
  10. C ES YOUN * SECT
  11. C COORD COORDONNEES DES 2 NOEUDS
  12. C
  13. C SORTIE
  14. C REL MATRICE DE RIGIDITE ELEMENTAIRE
  15. C KERRE 0 SI PAS DE PB
  16. C 1 SI LONGUEUR NULLE
  17. C-----------------------------------------------------------------------
  18. IMPLICIT INTEGER(I-N)
  19. IMPLICIT REAL*8(A-H,O-Z)
  20. -INC CCOPTIO
  21. DIMENSION REL(LRE,LRE),TGL(3,LRE),RE(3,3),C(3),B(3),XE(3,3)
  22. C
  23. DO 13 I=1,3
  24. DO 13 J=1,LRE
  25. TGL(I,J)=0.D0
  26. 13 CONTINUE
  27. C
  28. KERRE=0
  29. B(3)=0.D0
  30. DO I=1,IDIM
  31. B(I)=XE(I,3)-XE(I,1)
  32. END DO
  33. XL=(B(1)*B(1)+B(2)*B(2)+B(3)*B(3))
  34. XL=SQRT(XL)
  35. IF(XL.LE.0.D0) GO TO 999
  36. C
  37. FAC=ES/XL/3.D0
  38. RE(1,1)= 7.D0*FAC
  39. RE(1,2)=-8.D0*FAC
  40. RE(1,3)= 1.D0*FAC
  41. RE(2,1)=-8.D0*FAC
  42. RE(2,2)=16.D0*FAC
  43. RE(2,3)=-8.D0*FAC
  44. RE(3,1)= 1.D0*FAC
  45. RE(3,2)=-8.D0*FAC
  46. RE(3,3)= 7.D0*FAC
  47. C
  48. DO 20 I=1,3
  49. DO 10 J=1,IDIM
  50. L1=(I-1)*IDIM+J
  51. TGL(I,L1)=B(J)/XL
  52. 10 CONTINUE
  53. 20 CONTINUE
  54. C
  55. DO 5 I=1,LRE
  56. DO 6 J1=1,3
  57. B(J1)=TGL(J1,I)
  58. 6 CONTINUE
  59. DO 1 J=1,3
  60. S=0.D0
  61. DO 2 K=1,3
  62. 2 S=S+B(K)*RE(K,J)
  63. 1 C(J)=S
  64. DO 3 J=1,LRE
  65. S=0.D0
  66. DO 4 K=1,3
  67. 4 S=S+C(K)*TGL(K,J)
  68. 3 REL(I,J)=S
  69. 5 CONTINUE
  70. *
  71. C--------------------------------------------
  72. * WRITE(*,*) '** MATRICE KLOC BAR3 (dans barig3_1)** ELEME =',IB
  73. * DO 111 IA=1,3
  74. * WRITE(*,1991) (RE(IA,IC),IC=1,3)
  75. * 111 CONTINUE
  76. * WRITE(*,*) '** KGLOB BAR3 (dans barig3_1)** ATDH=',ATDH
  77. * DO 114 IA=1,LRE
  78. * WRITE(*,1991) (REL(IA,IC),IC=1,LRE)
  79. * 114 CONTINUE
  80. *1991 FORMAT(18(1X,E12.5))
  81. C--------------------------------------------
  82. RETURN
  83. 999 KERRE=1
  84. RETURN
  85. END
  86.  
  87.  
  88.  
  89.  
  90.  

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