1 : C PORIMO SOURCE AM 08/05/19 21:15:11 6108 2 : SUBROUTINE PORIMO(RP,XL,YOUN,CISAIL,SD,SREDY,SREDZ, 3 : . XINX,XINY,XINZ,IMAT,LHOOK,DDHOOK) 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 : DIMENSION DDHOOK(LHOOK,*) 29 : C 30 : C INITIALISATIONS 31 : C 32 : DO 120 I=1,144 33 : 120 RP(I,1)=0.D00 34 : RL=1.D00/XL 35 : RL2=RL*RL 36 : RL3=RL2*RL 37 : X2=0.D00 38 : X3=0.D00 39 : 40 : IF(IMAT.EQ.1) THEN 41 : IF( SREDZ. NE.0.) X2= YOUN/CISAIL*XINY*RL2/ SREDZ 42 : IF( SREDY. NE.0.) X3= YOUN/CISAIL*XINZ*RL2/ SREDY 43 : ESL=YOUN*SD*RL 44 : EIZ=YOUN*XINZ 45 : EIY=YOUN*XINY 46 : GJL=CISAIL*XINX*RL 47 : ELSE IF(IMAT.EQ.2) THEN 48 : IF(DDHOOK(3,3).NE.0.) X2= DDHOOK(5,5)*RL2/ DDHOOK(3,3) 49 : IF(DDHOOK(2,2).NE.0.) X3= DDHOOK(6,6)*RL2/ DDHOOK(2,2) 50 : ESL=DDHOOK(1,1)*RL 51 : EIZ=DDHOOK(6,6) 52 : EIY=DDHOOK(5,5) 53 : GJL=DDHOOK(4,4)*RL 54 : ENDIF 55 : 56 : X12= 1./(1.+12.*X2) 57 : X13= 1./(1.+12.*X3) 58 : 59 : C 60 : C REMPLISSAGE DE LA MATRICE 61 : C 62 : RP(1,1)=ESL 63 : RP(1,7)=-ESL 64 : RP(2,2)= 12.*EIZ*RL3*X13 65 : RP(2,6)=6.*EIZ*RL2*X13 66 : RP(2,8)= -12. *EIZ*RL3*X13 67 : RP(2,12)= 6.*EIZ*RL2*X13 68 : RP(3,3)=12.*EIY*RL3*X12 69 : RP(3,5)= -6.*EIY*RL2*X12 70 : RP(3,9)= -12.*EIY*RL3*X12 71 : RP(3,11)= -6.*EIY*RL2*X12 72 : RP(4,4)= GJL 73 : RP(4,10)=-GJL 74 : RP(5,5)= 4.*EIY*RL 75 : . *(1.+3.*X2)*X12 76 : RP(5,9)= 6.*EIY*RL2*X12 77 : RP(5,11)=2.*EIY*RL 78 : .*(1.-6.*X2)*X12 79 : RP(6,6)= 4.*EIZ*RL 80 : . *(1.+3.*X3)*X13 81 : RP(6,8)=-6.*EIZ*RL2*X13 82 : RP(6,12)=2.*EIZ*RL 83 : . *(1.-6.*X3)*X13 84 : RP(7,7)=ESL 85 : RP(8,8)= 12.*EIZ*RL3*X13 86 : RP(8,12)=-6.*EIZ*RL2*X13 87 : RP(9,9)= 12.*EIY*RL3*X12 88 : RP(9,11)= 6.*EIY*RL2*X12 89 : RP(10,10)= GJL 90 : RP(11,11)= 4.*EIY*RL 91 : . *(1.+3.*X2)*X12 92 : RP(12,12)= 4.*EIZ*RL 93 : . *(1.+3.*X3)*X13 94 : DO 9 I=1,12 95 : DO 9 J=I,12 96 : 9 RP(J,I)=RP(I,J) 97 : RETURN 98 : END 99 : 100 : 101 :
© Cast3M 2003 - All rights reserved.
Disclaimer