Télécharger la plaquette Cast3M


fr en

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