C BAREP3 SOURCE CHAT 05/01/12 21:33:47 5004 C BAREP3 SOURCE DJER 89/05/26 22:54:40 SUBROUTINE BAREP3(XE,DEPLA,DEFO,QSIGAU,POIGAU,NBPGAU,IB) C----------------------------------------------------------------------- C ROUTINE DE CALCUL DES DEFORMATIONS POUR LA BARRE C C ENTREE C DEPLA DEPLACEMENTS C XE XEONNEES NOEUDS C C SORTIE C DEFO DEFORMATION AUX POINTS D INTEGRATION C----------------------------------------------------------------------- IMPLICIT INTEGER(I-N) IMPLICIT REAL*8(A-H,O-Z) -INC PPARAM -INC CCOPTIO DIMENSION DEFO(*),XE(3,*),V(3),B3(3),DEPLA(*),U(3) DIMENSION QSIGAU(*),POIGAU(*) C KERRE=0 V(3)=0.D0 DO 1 I=1,IDIM 1 V(I)=XE(I,3)-XE(I,1) XL=(V(1)*V(1)+V(2)*V(2)+V(3)*V(3)) XL=SQRT(XL) IF(XL.EQ.0.) GO TO 999 * DO 2 I=1,3 U(I)=0.D0 DO 3 J=1,IDIM IG=(I-1)*IDIM+J U(I)=U(I)+(DEPLA(IG)*V(J))/XL 3 CONTINUE 2 CONTINUE * DO 5 IGAU=1,NBPGAU X1=XL/2.D0*(QSIGAU(IGAU)+1.D0) B3(1)=(4.D0*X1/XL-3.D0)/XL B3(2)=(4.D0-8.D0*X1/XL)/XL B3(3)=(4.D0*X1/XL-1.D0)/XL DEFO(IGAU)=B3(1)*U(1)+B3(2)*U(2)+B3(3)*U(3) 5 CONTINUE C-------------------------------------------- * NRE=LRE * WRITE(*,*) '** DEPLA GLOBAL BAR3 (dans barep3)** ELEME:',IB * WRITE(*,1991) (DEPLA(IC),IC=1,NRE) * WRITE(*,*) '** U LOCAL BAR3 (dans barep3)** ' * WRITE(*,1991) (U(IC),IC=1,3) * WRITE(*,*) '** QSIGAU NBPGAU=',NBPGAU,'XL=',XL * WRITE(*,1991) (QSIGAU(IC),IC=1,NBPGAU) * WRITE(*,*) '** POIGAU ' * WRITE(*,1991) (POIGAU(IC),IC=1,NBPGAU) * WRITE(*,*) '** B BAR3 (dans barep3)** ' * WRITE(*,1991) (B3(IC),IC=1,3) * WRITE(*,*) '** DEFO BAR3 (dans barep3)** ' * WRITE(*,1991) (DEFO(IC),IC=1,3) *1991 FORMAT(18(1X,E12.5)) C-------------------------------------------- 100 CONTINUE RETURN 999 KERRE=1 RETURN END