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