1 : C POKSL2 SOURCE CHAT 05/01/13 02:16:58 5004 2 : C POKSL2 SOURCE CASTEM 85/06/13 12:33:17 3 : SUBROUTINE POKSL2(RP,XL,SIGM) 4 : C----------------------------------------------------------------------- 5 : C 6 : C ROUTINE DE CALCUL DE LA MATRICE DE RIGIDITE GEOMETRIQUE 7 : C ELEMENTAIRE D UNE POUTRE DROITE EN AXES LOCAUX 8 : C D APRES TEDEL 9 : C 10 : C ENTREE 11 : C XL LONGUEUR 12 : C SIGM EFFORT NORMAL MOYEN 13 : C SORTIE 14 : C RP MATRICE DE RIGIDITE GEOMETRIQUE ELEMENTAIRE 12X12 15 : C 16 : C MAI 85 MILLARD 17 : C----------------------------------------------------------------------- 18 : IMPLICIT INTEGER(I-N) 19 : IMPLICIT REAL*8(A-H,O-Z) 20 : DIMENSION RP(6,6) 21 : C 22 : C INITIALISATIONS 23 : C 24 : DO 120 I=1,36 25 : 120 RP(I,1)=0.D00 26 : RL=1.D00/XL 27 : RN=SIGM*RL 28 : U1=1.2D00*RN 29 : U2=0.1D00*SIGM 30 : U3=0.13333333333333333D00*SIGM*XL 31 : U4=0.03333333333333333D00*SIGM*XL 32 : RP(1 ,1)=RN 33 : RP(1,4)=-RN 34 : RP(4,4)=RN 35 : RP(2,2)=U1 36 : RP(2,3)=U2 37 : RP(2,5)=-U1 38 : RP(2,6)=U2 39 : RP(3,3)=U3 40 : RP(3,5)=-U2 41 : RP(3,6)=-U4 42 : RP(5,5)=U1 43 : RP(5,6)=-U2 44 : RP(6,6)=U3 45 : C 46 : C REMPLISSAGE DE LA MATRICE 47 : C 48 : DO 9 I=1,6 49 : DO 9 J=I,6 50 : 9 RP(J,I)=RP(I,J) 51 : RETURN 52 : END 53 : 54 :
© Cast3M 2003 - All rights reserved.
Disclaimer