Télécharger la plaquette Cast3M


fr en

Télécharger podhl3.eso

Retour à la liste

Numérotation des lignes : 

podhl3

  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.  
  103.  

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