1 : C PORILO SOURCE AM 08/05/19 21:15:08 6108 2 : SUBROUTINE PORILO(RP,XL,YOUN,CISAIL,SD,SREDY,SREDZ, 3 : . XINX,XINY,XINZ) 4 : C----------------------------------------------------------------------- 5 : C 6 : C ROUTINE DE CALCUL DE LA MATRICE DE RIGIDITE ELEMENTAIRE D UNE 7 : C POUTRE DROITE EN AXES LOCAUX 8 : C D APRES TEDEL 9 : C 10 : C ENTREE 11 : C XL LONGUEUR 12 : C SD SECTION DROITE 13 : C SREDY SECTION REDUITE A L EFFORT TRANCHANT TY 14 : C SREDZ SECTION REDUITE A L EFFORT TRANCHANT TZ 15 : C XINX INERTIE IX 16 : C XINY INERTIE IY 17 : C XINZ INERTIE IZ 18 : C YOUN MODULE D YOUNG 19 : C CISAIL MODULE DE CISAILLEMENT 20 : C SORTIE 21 : C RP MATRICE DE RIGIDITE ELEMENTAIRE 12X12 22 : C 23 : C VIENT DE TEDEL AOUT 84 24 : C----------------------------------------------------------------------- 25 : IMPLICIT INTEGER(I-N) 26 : IMPLICIT REAL*8(A-H,O-Z) 27 : DIMENSION RP(12,12) 28 : C 29 : C INITIALISATIONS 30 : C 31 : DO 120 I=1,144 32 : 120 RP(I,1)=0.D00 33 : RL=1.D00/XL 34 : RL2=RL*RL 35 : RL3=RL2*RL 36 : X2=0.D00 37 : X3=0.D00 38 : 39 : IF( SREDZ. NE.0.) X2= YOUN/CISAIL*XINY*RL2/ SREDZ 40 : IF( SREDY. NE.0.) X3= YOUN/CISAIL*XINZ*RL2/ SREDY 41 : ESL=YOUN*SD*RL 42 : EIZ=YOUN*XINZ 43 : EIY=YOUN*XINY 44 : GJL=CISAIL*XINX*RL 45 : 46 : X12= 1./(1.+12.*X2) 47 : X13= 1./(1.+12.*X3) 48 : 49 : C 50 : C REMPLISSAGE DE LA MATRICE 51 : C 52 : RP(1,1)=ESL 53 : RP(1,7)=-ESL 54 : RP(2,2)= 12.*EIZ*RL3*X13 55 : RP(2,6)=6.*EIZ*RL2*X13 56 : RP(2,8)= -12. *EIZ*RL3*X13 57 : RP(2,12)= 6.*EIZ*RL2*X13 58 : RP(3,3)=12.*EIY*RL3*X12 59 : RP(3,5)= -6.*EIY*RL2*X12 60 : RP(3,9)= -12.*EIY*RL3*X12 61 : RP(3,11)= -6.*EIY*RL2*X12 62 : RP(4,4)= GJL 63 : RP(4,10)=-GJL 64 : RP(5,5)= 4.*EIY*RL 65 : . *(1.+3.*X2)*X12 66 : RP(5,9)= 6.*EIY*RL2*X12 67 : RP(5,11)=2.*EIY*RL 68 : .*(1.-6.*X2)*X12 69 : RP(6,6)= 4.*EIZ*RL 70 : . *(1.+3.*X3)*X13 71 : RP(6,8)=-6.*EIZ*RL2*X13 72 : RP(6,12)=2.*EIZ*RL 73 : . *(1.-6.*X3)*X13 74 : RP(7,7)=ESL 75 : RP(8,8)= 12.*EIZ*RL3*X13 76 : RP(8,12)=-6.*EIZ*RL2*X13 77 : RP(9,9)= 12.*EIY*RL3*X12 78 : RP(9,11)= 6.*EIY*RL2*X12 79 : RP(10,10)= GJL 80 : RP(11,11)= 4.*EIY*RL 81 : . *(1.+3.*X2)*X12 82 : RP(12,12)= 4.*EIZ*RL 83 : . *(1.+3.*X3)*X13 84 : DO 9 I=1,12 85 : DO 9 J=I,12 86 : 9 RP(J,I)=RP(I,J) 87 : RETURN 88 : END 89 : 90 : 91 :
© Cast3M 2003 - All rights reserved.
Disclaimer