Sources Esope | Cast3M

siglia
C SIGLIA    SOURCE    PV        22/04/15    13:20:13     11344          
      SUBROUTINE SIGLIA(BPSS,XDDL,WORK,LRE,NBNN,XSTRS)
C====================================================================
C     CALCULE DES CONTRAINTES  POUR LIA2 :
C      element de liaison a 2 noeuds (6 ddl par noeuds)
C     (ON CALCUL en fait DES EFFORTs NOEUDAUX )
C====================================================================
      IMPLICIT INTEGER(I-N)
      IMPLICIT REAL*8(A-H,O-Z)
      DIMENSION BPSS(3,3),XDDL(LRE),WORK(9),XSTRS(*),DEPL(12)
*-------------------------------------------------------------
      CALL ZERO(DEPL,12,1)
C
      DO 3 L=1,4
         L1 = 3*(L-1)
      DO 1 I=1,2
      DO  K=1,3
        DEPL(I+L1)=DEPL(I+L1)+BPSS(I,K)*XDDL(K+L1)
      enddo
    1 CONTINUE
    3 CONTINUE
C
      DO I =1,6
        XSTRS(I)  =WORK(I)*DEPL(I)
        XSTRS(I+6)=WORK(I)*DEPL(I+6)
      END DO
*
*-------------------------------------------------------------
*      WRITE(*,*) '**  EFFORT  ( dans siglia) ** '
*      WRITE(*,1992) (XSTRS(IC),IC=1,12)
*1992  FORMAT(6(6(1X,E12.5),/))
*-------------------------------------------------------------
      RETURN
      END