Sources Esope | Cast3M

barig3
C BARIG3    SOURCE    CB215821  26/08/24    21:15:11     12622          
C BARIG3    SOURCE    CHAT      89/05/26    22:52:13
      SUBROUTINE BARIG3(REL,LRE,ES,XE,TGL,KERRE,IB)
C-----------------------------------------------------------------------
C       ROUTINE DE CALCUL DE LA MATRICE DE RIGIDITE ELEMENTAIRE
C                      DE L ELEMENT DE BARRE A 3 NOEUDS
C-----------------------------------------------------------------------
C  ENTREE
C          LRE     DIMENSION DU TABLEAU REL
C          ES      YOUN * SECT
C          COORD   COORDONNEES DES 2 NOEUDS
C
C  SORTIE
C          REL     MATRICE DE RIGIDITE ELEMENTAIRE
C          KERRE   0 SI PAS DE PB
C                  1 SI LONGUEUR NULLE
C-----------------------------------------------------------------------
      IMPLICIT INTEGER(I-N)
      IMPLICIT REAL*8(A-H,O-Z)
 
-INC PPARAM
-INC CCOPTIO
      DIMENSION REL(LRE,LRE),TGL(3,LRE),RE(3,3),C(3),B(3),XE(3,3)
C
      DO 1000 I=1,3
      DO 13 J=1,LRE
       TGL(I,J)=0.D0
   13 CONTINUE
 1000 CONTINUE
C
      KERRE=0
      B(3)=0.D0
      DO I=1,IDIM
        B(I)=XE(I,3)-XE(I,1)
      END DO
      XL=(B(1)*B(1)+B(2)*B(2)+B(3)*B(3))
      XL=SQRT(XL)
      IF(XL.LE.0.D0) GO TO 999
C
      FAC=ES/XL/3.D0
      RE(1,1)= 7.D0*FAC
      RE(1,2)=-8.D0*FAC
      RE(1,3)= 1.D0*FAC
      RE(2,1)=-8.D0*FAC
      RE(2,2)=16.D0*FAC
      RE(2,3)=-8.D0*FAC
      RE(3,1)= 1.D0*FAC
      RE(3,2)=-8.D0*FAC
      RE(3,3)= 7.D0*FAC
C
      DO 20 I=1,3
      DO 10 J=1,IDIM
        L1=(I-1)*IDIM+J
        TGL(I,L1)=B(J)/XL
   10 CONTINUE
   20 CONTINUE
C
      DO 5 I=1,LRE
      DO 6 J1=1,3
        B(J1)=TGL(J1,I)
    6 CONTINUE
      DO 1 J=1,3
      S=0.D0
      DO 2 K=1,3
      S=S+B(K)*RE(K,J)
    2 CONTINUE
      C(J)=S
    1 CONTINUE
      DO 3 J=1,LRE
      S=0.D0
      DO 4 K=1,3
      S=S+C(K)*TGL(K,J)
    4 CONTINUE
      REL(I,J)=S
    3 CONTINUE
    5 CONTINUE
*
C--------------------------------------------
*      WRITE(*,*) '** MATRICE  KLOC BAR3 (dans barig3_1)** ELEME =',IB
*      DO 111 IA=1,3
*       WRITE(*,1991) (RE(IA,IC),IC=1,3)
*  111 CONTINUE
*      WRITE(*,*) '**   KGLOB BAR3 (dans barig3_1)** ATDH=',ATDH
*      DO 114 IA=1,LRE
*       WRITE(*,1991) (REL(IA,IC),IC=1,LRE)
*  114 CONTINUE
*1991  FORMAT(18(1X,E12.5))
C--------------------------------------------
      RETURN
 999  KERRE=1
      RETURN
      END