1 : C PORIM2 SOURCE AM 08/05/19 21:15:10 6108 2 : SUBROUTINE PORIM2(RP,XL,YOUN,CISAIL,SD,SREDY,XINZ, 3 : & 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 EN 2D 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 XINZ INERTIE IZ 15 : C YOUN MODULE D YOUNG 16 : C CISAIL MODULE DE CISAILLEMENT 17 : C SORTIE 18 : C RP MATRICE DE RIGIDITE ELEMENTAIRE 12X12 19 : C 20 : C VIENT DE TEDEL AOUT 84 21 : C----------------------------------------------------------------------- 22 : IMPLICIT INTEGER(I-N) 23 : IMPLICIT REAL*8(A-H,O-Z) 24 : DIMENSION RP(6,6) 25 : DIMENSION DDHOOK(LHOOK,*) 26 : C 27 : C INITIALISATIONS 28 : C 29 : DO 120 I=1,36 30 : 120 RP(I,1)=0.D00 31 : RL=1.D00/XL 32 : RL2=RL*RL 33 : RL3=RL2*RL 34 : X2=0.D00 35 : X3=0.D00 36 : 37 : IF(IMAT.EQ.1) THEN 38 : IF( SREDY. NE.0.) X3= YOUN/CISAIL*XINZ*RL2/ SREDY 39 : ESL=YOUN*SD*RL 40 : EIZ=YOUN*XINZ 41 : ELSE IF(IMAT.EQ.2) THEN 42 : IF(DDHOOK(2,2).NE.0.) X3= DDHOOK(3,3)*RL2/ DDHOOK(2,2) 43 : ESL=DDHOOK(1,1)*RL 44 : EIZ=DDHOOK(3,3) 45 : ENDIF 46 : 47 : X12= 1./(1.+12.*X2) 48 : X13= 1./(1.+12.*X3) 49 : 50 : C 51 : C REMPLISSAGE DE LA MATRICE 52 : C 53 : RP(1,1)=ESL 54 : RP(1,4)=-ESL 55 : RP(2,2)= 12.*EIZ*RL3*X13 56 : RP(2,3)=6.*EIZ*RL2*X13 57 : RP(2,5)= -12. *EIZ*RL3*X13 58 : RP(2,6)= 6.*EIZ*RL2*X13 59 : RP(3,3)= 4.*EIZ*RL 60 : . *(1.+3.*X3)*X13 61 : RP(3,5)=-6.*EIZ*RL2*X13 62 : RP(3,6)=2.*EIZ*RL 63 : . *(1.-6.*X3)*X13 64 : RP(4,4)=ESL 65 : RP(5,5)= 12.*EIZ*RL3*X13 66 : RP(5,6)=-6.*EIZ*RL2*X13 67 : RP(6,6)= 4.*EIZ*RL 68 : . *(1.+3.*X3)*X13 69 : DO 9 I=1,6 70 : DO 9 J=I,6 71 : 9 RP(J,I)=RP(I,J) 72 : RETURN 73 : END 74 : 75 : 76 : 77 :
© Cast3M 2003 - All rights reserved.
Disclaimer