Esope source files

   1 : C SIGLIA    SOURCE    CHAT      05/01/13    03:16:39     5004           
   2 : C SIGLIA    SOURCE    DJER     23/01/96
   3 :       SUBROUTINE SIGLIA(BPSS,XDDL,WORK,LRE,NBNN,XSTRS)
   4 : C====================================================================           
   5 : C     CALCULE DES CONTRAINTES  POUR LIA2 :
   6 : C      element de liaison a 2 noeuds (6 ddl par noeuds)                                                        
   7 : C     (ON CALCUL en fait DES EFFORTs NOEUDAUX ) 
   8 : C====================================================================           
   9 :       IMPLICIT INTEGER(I-N)
  10 :       IMPLICIT REAL*8(A-H,O-Z)                                                  
  11 :       DIMENSION BPSS(3,3),XDDL(LRE),WORK(9),XSTRS(*),DEPL(12)
  12 : *-------------------------------------------------------------
  12 :       CALL ZERO(DEPL,12,12)
  13 : C
  14 :       DO 3 L=1,4
  15 :          L1 = 3*(L-1)
  16 :       DO 1 I=1,2
  17 :       DO 2 K=1,3
  18 :     2   DEPL(I+L1)=DEPL(I+L1)+BPSS(I,K)*XDDL(K+L1)
  19 :     1 CONTINUE
  20 :     3 CONTINUE
  21 : C      
  22 :       DO I =1,6
  23 :         XSTRS(I)  =WORK(I)*DEPL(I)
  24 :         XSTRS(I+6)=WORK(I)*DEPL(I+6)
  25 :       END DO
  26 : *
  27 : *-------------------------------------------------------------
  28 : *      WRITE(*,*) '**  EFFORT  ( dans siglia) ** '
  29 : *      WRITE(*,1992) (XSTRS(IC),IC=1,12)
  30 : *1992  FORMAT(6(6(1X,E12.5),/))
  31 : *-------------------------------------------------------------
  32 :       RETURN                                                                    
  33 :       END                                                                       
  34 :  
  35 :