1 : C BARIG3 SOURCE CHAT 05/01/12 21:34:03 5004 2 : C BARIG3 SOURCE CHAT 89/05/26 22:52:13 3 : SUBROUTINE BARIG3(REL,LRE,ES,XE,TGL,KERRE,IB) 4 : C----------------------------------------------------------------------- 5 : C ROUTINE DE CALCUL DE LA MATRICE DE RIGIDITE ELEMENTAIRE 6 : C DE L ELEMENT DE BARRE A 3 NOEUDS 7 : C----------------------------------------------------------------------- 8 : C ENTREE 9 : C LRE DIMENSION DU TABLEAU REL 10 : C ES YOUN * SECT 11 : C COORD COORDONNEES DES 2 NOEUDS 12 : C 13 : C SORTIE 14 : C REL MATRICE DE RIGIDITE ELEMENTAIRE 15 : C KERRE 0 SI PAS DE PB 16 : C 1 SI LONGUEUR NULLE 17 : C----------------------------------------------------------------------- 18 : IMPLICIT INTEGER(I-N) 19 : IMPLICIT REAL*8(A-H,O-Z) 20 : -INC CCOPTIO 21 : DIMENSION REL(LRE,LRE),TGL(3,LRE),RE(3,3),C(3),B(3),XE(3,3) 22 : C 23 : DO 13 I=1,3 24 : DO 13 J=1,LRE 25 : TGL(I,J)=0.D0 26 : 13 CONTINUE 27 : C 28 : KERRE=0 29 : B(3)=0.D0 30 : DO I=1,IDIM 31 : B(I)=XE(I,3)-XE(I,1) 32 : END DO 33 : XL=(B(1)*B(1)+B(2)*B(2)+B(3)*B(3)) 34 : XL=DSQRT(XL) 35 : IF(XL.LE.0.D0) GO TO 999 36 : C 37 : FAC=ES/XL/3.D0 38 : RE(1,1)= 7.D0*FAC 39 : RE(1,2)=-8.D0*FAC 40 : RE(1,3)= 1.D0*FAC 41 : RE(2,1)=-8.D0*FAC 42 : RE(2,2)=16.D0*FAC 43 : RE(2,3)=-8.D0*FAC 44 : RE(3,1)= 1.D0*FAC 45 : RE(3,2)=-8.D0*FAC 46 : RE(3,3)= 7.D0*FAC 47 : C 48 : DO 20 I=1,3 49 : DO 10 J=1,IDIM 50 : L1=(I-1)*IDIM+J 51 : TGL(I,L1)=B(J)/XL 52 : 10 CONTINUE 53 : 20 CONTINUE 54 : C 55 : DO 5 I=1,LRE 56 : DO 6 J1=1,3 57 : B(J1)=TGL(J1,I) 58 : 6 CONTINUE 59 : DO 1 J=1,3 60 : S=0.D0 61 : DO 2 K=1,3 62 : 2 S=S+B(K)*RE(K,J) 63 : 1 C(J)=S 64 : DO 3 J=1,LRE 65 : S=0.D0 66 : DO 4 K=1,3 67 : 4 S=S+C(K)*TGL(K,J) 68 : 3 REL(I,J)=S 69 : 5 CONTINUE 70 : * 71 : C-------------------------------------------- 72 : * WRITE(*,*) '** MATRICE KLOC BAR3 (dans barig3_1)** ELEME =',IB 73 : * DO 111 IA=1,3 74 : * WRITE(*,1991) (RE(IA,IC),IC=1,3) 75 : * 111 CONTINUE 76 : * WRITE(*,*) '** KGLOB BAR3 (dans barig3_1)** ATDH=',ATDH 77 : * DO 114 IA=1,LRE 78 : * WRITE(*,1991) (REL(IA,IC),IC=1,LRE) 79 : * 114 CONTINUE 80 : *1991 FORMAT(18(1X,E12.5)) 81 : C-------------------------------------------- 82 : RETURN 83 : 999 KERRE=1 84 : RETURN 85 : END 86 : 87 :
© Cast3M 2003 - All rights reserved.
Disclaimer