1 : C POMAL2 SOURCE CHAT 05/01/13 02:17:24 5004 2 : C POMAL2 SOURCE CASTEM 85/06/13 12:33:23 3 : SUBROUTINE POMAL2(RP,XL,RHO,SD,XINZ) 4 : C----------------------------------------------------------------------- 5 : C 6 : C ROUTINE DE CALCUL DE LA MATRICE DE MASSE ELEMENTAIRE D UNE 7 : C POUTRE DROITE EN AXES LOCAUX 8 : C D APRES TEDEL 9 : C DANS CETTE VERSION,ON PREND EN COMPTE LES INERTIES DE ROTATION 10 : C MAIS PAS LA DEFORMABILITE A L EFFORT TRANCHANT 11 : C 12 : C ENTREE 13 : C XL LONGUEUR 14 : C SD SECTION DROITE 15 : C XINZ INERTIE IZ 16 : C RHO MASSE VOLUMIQUE 17 : C SORTIE 18 : C RP MATRICE DE MASSE ELEMENTAIRE 12X12 19 : C 20 : C VIENT DE TEDEL AOUT 84 21 : C 22 : C----------------------------------------------------------------------- 23 : IMPLICIT INTEGER(I-N) 24 : IMPLICIT REAL*8(A-H,O-Z) 25 : DIMENSION RP(6,6) 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 : XL2=XL*XL 33 : RHOSL=RHO*SD*XL 34 : RSP1= RHOSL*0.333333333333333333 35 : RSP2= RHOSL*0.166666666666666666 36 : U 1= 0.3714285714285714D+00 * RHOSL 37 : U 2= 0.5238095238095238D-01 * RHOSL *XL 38 : U 3= 0.1285714285714286D+00 * RHOSL 39 : U 4= -0.3095238095238095D-01 * RHOSL *XL 40 : U 5= 0.9523809523809524D-02 * RHOSL *XL2 41 : U 6= 0.3095238095238095D-01 * RHOSL *XL 42 : U 7= -0.7142857142857143D-02 * RHOSL *XL2 43 : U 8= 0.3714285714285714D+00 * RHOSL 44 : U 9= -0.5238095238095238D-01 * RHOSL *XL 45 : U10= 0.9523809523809524D-02 * RHOSL *XL2 46 : U11=1.2*RHO*RL*XINZ 47 : U12=0.1*RHO*XINZ 48 : U18=0.13333333333333333*XINZ*RHO*XL 49 : U19=0.1*XINZ*RHO 50 : U20=0.03333333333333333*XINZ*XL*RHO 51 : C 52 : RP(1,1)=RSP1 53 : RP(1,4)=RSP2 54 : RP(2,2)=U1+U11 55 : RP(2,3)=U2+U12 56 : RP(2,5)=U3-U11 57 : RP(2,6)=U4+U12 58 : RP(3,3)=U5+U18 59 : RP(3,5)=U6-U19 60 : RP(3,6)=U7-U20 61 : RP(4,4)=RSP1 62 : RP(5,5)=U8+U11 63 : RP(5,6)=U9-U12 64 : RP(6,6)=U10+U18 65 : DO 9 I=1,6 66 : DO 9 J=I,6 67 : 9 RP(J,I)= RP(I,J) 68 : RETURN 69 : END 70 : 71 :
© Cast3M 2003 - All rights reserved.
Disclaimer