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