Télécharger barig3.eso

Retour à la liste

Numérotation des lignes :

barig3
  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.  
  21. -INC PPARAM
  22. -INC CCOPTIO
  23. DIMENSION REL(LRE,LRE),TGL(3,LRE),RE(3,3),C(3),B(3),XE(3,3)
  24. C
  25. DO 13 I=1,3
  26. DO 13 J=1,LRE
  27. TGL(I,J)=0.D0
  28. 13 CONTINUE
  29. C
  30. KERRE=0
  31. B(3)=0.D0
  32. DO I=1,IDIM
  33. B(I)=XE(I,3)-XE(I,1)
  34. END DO
  35. XL=(B(1)*B(1)+B(2)*B(2)+B(3)*B(3))
  36. XL=SQRT(XL)
  37. IF(XL.LE.0.D0) GO TO 999
  38. C
  39. FAC=ES/XL/3.D0
  40. RE(1,1)= 7.D0*FAC
  41. RE(1,2)=-8.D0*FAC
  42. RE(1,3)= 1.D0*FAC
  43. RE(2,1)=-8.D0*FAC
  44. RE(2,2)=16.D0*FAC
  45. RE(2,3)=-8.D0*FAC
  46. RE(3,1)= 1.D0*FAC
  47. RE(3,2)=-8.D0*FAC
  48. RE(3,3)= 7.D0*FAC
  49. C
  50. DO 20 I=1,3
  51. DO 10 J=1,IDIM
  52. L1=(I-1)*IDIM+J
  53. TGL(I,L1)=B(J)/XL
  54. 10 CONTINUE
  55. 20 CONTINUE
  56. C
  57. DO 5 I=1,LRE
  58. DO 6 J1=1,3
  59. B(J1)=TGL(J1,I)
  60. 6 CONTINUE
  61. DO 1 J=1,3
  62. S=0.D0
  63. DO 2 K=1,3
  64. 2 S=S+B(K)*RE(K,J)
  65. 1 C(J)=S
  66. DO 3 J=1,LRE
  67. S=0.D0
  68. DO 4 K=1,3
  69. 4 S=S+C(K)*TGL(K,J)
  70. 3 REL(I,J)=S
  71. 5 CONTINUE
  72. *
  73. C--------------------------------------------
  74. * WRITE(*,*) '** MATRICE KLOC BAR3 (dans barig3_1)** ELEME =',IB
  75. * DO 111 IA=1,3
  76. * WRITE(*,1991) (RE(IA,IC),IC=1,3)
  77. * 111 CONTINUE
  78. * WRITE(*,*) '** KGLOB BAR3 (dans barig3_1)** ATDH=',ATDH
  79. * DO 114 IA=1,LRE
  80. * WRITE(*,1991) (REL(IA,IC),IC=1,LRE)
  81. * 114 CONTINUE
  82. *1991 FORMAT(18(1X,E12.5))
  83. C--------------------------------------------
  84. RETURN
  85. 999 KERRE=1
  86. RETURN
  87. END
  88.  
  89.  
  90.  
  91.  
  92.  

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