Télécharger la plaquette Cast3M


fr en

Télécharger pomal2.eso

Retour à la liste

Numérotation des lignes : 

pomal2

  1. C POMAL2    SOURCE    CB215821  26/08/24    21:17:46     12622          
  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 121 J=1,6
  30.       DO 120 I=1,6
  31.       RP(I,J)=0.D00
  32.   120 CONTINUE
  33.   121 CONTINUE
  34.       RL=1.D00/XL
  35.       XL2=XL*XL
  36.       RHOSL=RHO*SD*XL
  37.       RSP1= RHOSL*0.333333333333333333
  38.       RSP2= RHOSL*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.       U18=0.13333333333333333*XINZ*RHO*XL
  52.       U19=0.1*XINZ*RHO
  53.       U20=0.03333333333333333*XINZ*XL*RHO
  54. C
  55.       RP(1,1)=RSP1
  56.       RP(1,4)=RSP2
  57.       RP(2,2)=U1+U11
  58.       RP(2,3)=U2+U12
  59.       RP(2,5)=U3-U11
  60.       RP(2,6)=U4+U12
  61.       RP(3,3)=U5+U18
  62.       RP(3,5)=U6-U19
  63.       RP(3,6)=U7-U20
  64.       RP(4,4)=RSP1
  65.       RP(5,5)=U8+U11
  66.       RP(5,6)=U9-U12
  67.       RP(6,6)=U10+U18
  68.       DO 122 I=1,6
  69.       DO 9 J=I,6
  70.       RP(J,I)= RP(I,J)
  71.     9 CONTINUE
  72.   122 CONTINUE
  73.       RETURN
  74.       END
  75.  
  76.  
  77.  
  78.  
  79.  

© Cast3M 2003 - Tous droits réservés.
Mentions légales