C SIGLIA SOURCE CHAT 05/01/13 03:16:39 5004 C SIGLIA SOURCE DJER 23/01/96 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,12) C DO 3 L=1,4 L1 = 3*(L-1) DO 1 I=1,2 DO 2 K=1,3 2 DEPL(I+L1)=DEPL(I+L1)+BPSS(I,K)*XDDL(K+L1) 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