1 : C PODHL3 SOURCE CHAT 05/01/13 02:16:14 5004 2 : SUBROUTINE PODHL3(RP,XL,VISC,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 VISC MODULE DE VISCOSITE 19 : C CISAIL MODULE DE VISCOSITE EN 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,12 32 : DO 120 J=1,12 33 : 120 RP(I,J)=0.D0 34 : C 35 : RL=1.D00/XL 36 : RL2=RL*RL 37 : X2=0.D00 38 : SRED = 0.5D0*(SREDY + SREDZ) 39 : XINYZ = 0.5D0*(XINY + XINZ) 40 : IF((SREDY.NE.0).AND.(SREDZ.NE.0)) 41 : . X2= VISC/CISAIL*XINYZ*RL2/SRED 42 : X12= 1./(1.+12.*X2) 43 : C 44 : RL3=RL2*RL 45 : EIYZ=VISC*XINYZ 46 : C 47 : C REMPLISSAGE DE LA MATRICE 48 : C 49 : C 50 : RP(2,3)= 12.*EIYZ*RL3*X12 51 : RP(3,2)= -12.*EIYZ*RL3*X12 52 : RP(2,5)= -6.*EIYZ*RL2*X12 53 : RP(5,2)= 6.*EIYZ*RL2*X12 54 : RP(2,9)= -12. *EIYZ*RL3*X12 55 : RP(9,2)= 12. *EIYZ*RL3*X12 56 : RP(2,11)= -6.*EIYZ*RL2*X12 57 : RP(11,2)= 6.*EIYZ*RL2*X12 58 : C 59 : RP(3,6)= -6.*EIYZ*RL2*X12 60 : RP(3,8)= 12.*EIYZ*RL3*X12 61 : RP(3,12)= -6.*EIYZ*RL2*X12 62 : RP(6,3)= 6.*EIYZ*RL2*X12 63 : RP(8,3)= -12.*EIYZ*RL3*X12 64 : RP(12,3)= 6.*EIYZ*RL2*X12 65 : C 66 : RP(5,6)= 4.*EIYZ*RL 67 : . *(1.+3.*X2)*X12 68 : RP(5,8)= -6.*EIYZ*RL2*X12 69 : RP(5,12)=2.*EIYZ*RL 70 : .*(1.-6.*X2)*X12 71 : RP(6,5)= -4.*EIYZ*RL 72 : . *(1.+3.*X2)*X12 73 : RP(8,5)= 6.*EIYZ*RL2*X12 74 : RP(12,5)=-2.*EIYZ*RL 75 : .*(1.-6.*X2)*X12 76 : C 77 : RP(6,9)=-6.*EIYZ*RL2*X12 78 : RP(6,11)=-2.*EIYZ*RL 79 : . *(1.-6.*X2)*X12 80 : RP(9,6)=6.*EIYZ*RL2*X12 81 : RP(11,6)=2.*EIYZ*RL 82 : . *(1.-6.*X2)*X12 83 : C 84 : RP(8,9)= 12.*EIYZ*RL3*X12 85 : RP(8,11)= 6.*EIYZ*RL2*X12 86 : RP(9,8)= -12.*EIYZ*RL3*X12 87 : RP(11,8)= -6.*EIYZ*RL2*X12 88 : C 89 : RP(9,12)= 6.*EIYZ*RL2*X12 90 : RP(12,9)= -6.*EIYZ*RL2*X12 91 : C 92 : RP(11,12)= 4.*EIYZ*RL 93 : . *(1.D0+3.*X2)*X12 94 : C 95 : RP(12,11)= -4.*EIYZ*RL 96 : . *(1.D0+3.*X2)*X12 97 : C 98 : RETURN 99 : END 100 : 101 : 102 :
© Cast3M 2003 - All rights reserved.
Disclaimer