1 : C BAREP3 SOURCE CHAT 05/01/12 21:33:47 5004 2 : C BAREP3 SOURCE DJER 89/05/26 22:54:40 3 : SUBROUTINE BAREP3(XE,DEPLA,DEFO,QSIGAU,POIGAU,NBPGAU,IB) 4 : C----------------------------------------------------------------------- 5 : C ROUTINE DE CALCUL DES DEFORMATIONS POUR LA BARRE 6 : C 7 : C ENTREE 8 : C DEPLA DEPLACEMENTS 9 : C XE XEONNEES NOEUDS 10 : C 11 : C SORTIE 12 : C DEFO DEFORMATION AUX POINTS D INTEGRATION 13 : C----------------------------------------------------------------------- 14 : IMPLICIT INTEGER(I-N) 15 : IMPLICIT REAL*8(A-H,O-Z) 16 : -INC CCOPTIO 17 : DIMENSION DEFO(*),XE(3,*),V(3),B3(3),DEPLA(*),U(3) 18 : DIMENSION QSIGAU(*),POIGAU(*) 19 : C 20 : KERRE=0 21 : V(3)=0.D0 22 : DO 1 I=1,IDIM 23 : 1 V(I)=XE(I,3)-XE(I,1) 24 : XL=(V(1)*V(1)+V(2)*V(2)+V(3)*V(3)) 25 : XL=SQRT(XL) 26 : IF(XL.EQ.0.) GO TO 999 27 : * 28 : DO 2 I=1,3 29 : U(I)=0.D0 30 : DO 3 J=1,IDIM 31 : IG=(I-1)*IDIM+J 32 : U(I)=U(I)+(DEPLA(IG)*V(J))/XL 33 : 3 CONTINUE 34 : 2 CONTINUE 35 : * 36 : DO 5 IGAU=1,NBPGAU 37 : X1=XL/2.D0*(QSIGAU(IGAU)+1.D0) 38 : B3(1)=(4.D0*X1/XL-3.D0)/XL 39 : B3(2)=(4.D0-8.D0*X1/XL)/XL 40 : B3(3)=(4.D0*X1/XL-1.D0)/XL 41 : DEFO(IGAU)=B3(1)*U(1)+B3(2)*U(2)+B3(3)*U(3) 42 : 5 CONTINUE 43 : C-------------------------------------------- 44 : * NRE=LRE 45 : * WRITE(*,*) '** DEPLA GLOBAL BAR3 (dans barep3)** ELEME:',IB 46 : * WRITE(*,1991) (DEPLA(IC),IC=1,NRE) 47 : * WRITE(*,*) '** U LOCAL BAR3 (dans barep3)** ' 48 : * WRITE(*,1991) (U(IC),IC=1,3) 49 : * WRITE(*,*) '** QSIGAU NBPGAU=',NBPGAU,'XL=',XL 50 : * WRITE(*,1991) (QSIGAU(IC),IC=1,NBPGAU) 51 : * WRITE(*,*) '** POIGAU ' 52 : * WRITE(*,1991) (POIGAU(IC),IC=1,NBPGAU) 53 : * WRITE(*,*) '** B BAR3 (dans barep3)** ' 54 : * WRITE(*,1991) (B3(IC),IC=1,3) 55 : * WRITE(*,*) '** DEFO BAR3 (dans barep3)** ' 56 : * WRITE(*,1991) (DEFO(IC),IC=1,3) 57 : *1991 FORMAT(18(1X,E12.5)) 58 : C-------------------------------------------- 59 : 100 CONTINUE 60 : RETURN 61 : 999 KERRE=1 62 : RETURN 63 : END 64 : 65 :
© Cast3M 2003 - All rights reserved.
Disclaimer