1 : C POMALO SOURCE CHAT 05/01/13 02:17:28 5004 2 : SUBROUTINE POMALO(RP,XL,RHO,SD,XINX,XINY,XINZ) 3 : C----------------------------------------------------------------------- 4 : C 5 : C ROUTINE DE CALCUL DE LA MATRICE DE MASSE ELEMENTAIRE D UNE 6 : C POUTRE DROITE EN AXES LOCAUX 7 : C D APRES TEDEL 8 : C DANS CETTE VERSION,ON PREND EN COMPTE LES INERTIES DE ROTATION 9 : C MAIS PAS LA DEFORMABILITE A L EFFORT TRANCHANT 10 : C 11 : C ENTREE 12 : C XL LONGUEUR 13 : C SD SECTION DROITE 14 : C XINX INERTIE IX 15 : C XINY INERTIE IY 16 : C XINZ INERTIE IZ 17 : C RHO MASSE VOLUMIQUE 18 : C SORTIE 19 : C RP MATRICE DE MASSE ELEMENTAIRE 12X12 20 : C 21 : C VIENT DE TEDEL AOUT 84 22 : C 23 : C----------------------------------------------------------------------- 24 : IMPLICIT INTEGER(I-N) 25 : IMPLICIT REAL*8(A-H,O-Z) 26 : DIMENSION RP(12,12) 27 : C 28 : C INITIALISATIONS 29 : C 30 : DO 120 I=1,144 31 : 120 RP(I,1)=0.D00 32 : RL=1.D00/XL 33 : XL2=XL*XL 34 : RHOSL=RHO*SD*XL 35 : RSP1= RHOSL*0.333333333333333333 36 : RSP2= RHOSL*0.166666666666666666 37 : TOR1= RHO*XINX*XL*0.333333333333333333 38 : TOR2= RHO*XINX*XL*0.166666666666666666 39 : U 1= 0.3714285714285714D+00 * RHOSL 40 : U 2= 0.5238095238095238D-01 * RHOSL *XL 41 : U 3= 0.1285714285714286D+00 * RHOSL 42 : U 4= -0.3095238095238095D-01 * RHOSL *XL 43 : U 5= 0.9523809523809524D-02 * RHOSL *XL2 44 : U 6= 0.3095238095238095D-01 * RHOSL *XL 45 : U 7= -0.7142857142857143D-02 * RHOSL *XL2 46 : U 8= 0.3714285714285714D+00 * RHOSL 47 : U 9= -0.5238095238095238D-01 * RHOSL *XL 48 : U10= 0.9523809523809524D-02 * RHOSL *XL2 49 : U11=1.2*RHO*RL*XINZ 50 : U12=0.1*RHO*XINZ 51 : U13=1.2*RHO*RL*XINY 52 : U14=0.1*RHO*XINY 53 : U15=0.13333333333333333*XINY*RHO*XL 54 : U16=0.1*XINY*RHO 55 : U17=0.03333333333333333*XINY*XL*RHO 56 : U18=0.13333333333333333*XINZ*RHO*XL 57 : U19=0.1*XINZ*RHO 58 : U20=0.03333333333333333*XINZ*XL*RHO 59 : C 60 : RP(1,1)=RSP1 61 : RP(1,7)=RSP2 62 : RP(2,2)=U1+U11 63 : RP(2,6)=U2+U12 64 : RP(2,8)=U3-U11 65 : RP(2,12)=U4+U12 66 : RP(3,3)=U1+U13 67 : RP(3,5)=-U2-U14 68 : RP(3,9)=U3-U13 69 : RP(3,11)=-U4-U14 70 : RP(4,4)= TOR1 71 : RP(4,10)= TOR2 72 : RP(5,5)=U5+U15 73 : RP(5,9)=-U6+U16 74 : RP(5,11)=U7-U17 75 : RP(6,6)=U5+U18 76 : RP(6,8)=U6-U19 77 : RP(6,12)=U7-U20 78 : RP(7,7)=RSP1 79 : RP(8,8)=U8+U11 80 : RP(8,12)=U9-U12 81 : RP(9,9)=U8+U13 82 : RP(9,11)=-U9+U14 83 : RP(10,10)=TOR1 84 : RP(11,11)=U10+U15 85 : RP(12,12)=U10+U18 86 : DO 9 I=1,12 87 : DO 9 J=I,12 88 : 9 RP(J,I)= RP(I,J) 89 : RETURN 90 : END 91 :
© Cast3M 2003 - All rights reserved.
Disclaimer