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    CB215821  26/08/24    21:17:42     12622          
  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 121 I=1,12
  32.       DO 120 J=1,12
  33.       RP(I,J)=0.D0
  34.   120 CONTINUE
  35.   121 CONTINUE
  36. C
  37.       RL=1.D00/XL
  38.       RL2=RL*RL
  39.       X2=0.D00
  40.       SRED = 0.5D0*(SREDY + SREDZ)
  41.       XINYZ = 0.5D0*(XINY + XINZ)
  42.       IF((SREDY.NE.0).AND.(SREDZ.NE.0))
  43.      .     X2= VISC/CISAIL*XINYZ*RL2/SRED
  44.       X12= 1./(1.+12.*X2)
  45. C
  46.       RL3=RL2*RL
  47.       EIYZ=VISC*XINYZ
  48. C
  49. C  REMPLISSAGE DE LA MATRICE
  50. C
  51. C
  52.        RP(2,3)= 12.*EIYZ*RL3*X12
  53.        RP(3,2)= -12.*EIYZ*RL3*X12
  54.        RP(2,5)=  -6.*EIYZ*RL2*X12
  55.        RP(5,2)=  6.*EIYZ*RL2*X12
  56.        RP(2,9)=  -12. *EIYZ*RL3*X12
  57.        RP(9,2)=  12. *EIYZ*RL3*X12
  58.        RP(2,11)= -6.*EIYZ*RL2*X12
  59.        RP(11,2)= 6.*EIYZ*RL2*X12
  60. C
  61.        RP(3,6)= -6.*EIYZ*RL2*X12
  62.        RP(3,8)=  12.*EIYZ*RL3*X12
  63.        RP(3,12)=  -6.*EIYZ*RL2*X12
  64.        RP(6,3)= 6.*EIYZ*RL2*X12
  65.        RP(8,3)=  -12.*EIYZ*RL3*X12
  66.        RP(12,3)=  6.*EIYZ*RL2*X12
  67. C
  68.        RP(5,6)= 4.*EIYZ*RL
  69.      . *(1.+3.*X2)*X12
  70.        RP(5,8)= -6.*EIYZ*RL2*X12
  71.        RP(5,12)=2.*EIYZ*RL
  72.      .*(1.-6.*X2)*X12
  73.        RP(6,5)= -4.*EIYZ*RL
  74.      . *(1.+3.*X2)*X12
  75.        RP(8,5)= 6.*EIYZ*RL2*X12
  76.        RP(12,5)=-2.*EIYZ*RL
  77.      .*(1.-6.*X2)*X12
  78. C
  79.       RP(6,9)=-6.*EIYZ*RL2*X12
  80.       RP(6,11)=-2.*EIYZ*RL
  81.      .      *(1.-6.*X2)*X12
  82.       RP(9,6)=6.*EIYZ*RL2*X12
  83.       RP(11,6)=2.*EIYZ*RL
  84.      .      *(1.-6.*X2)*X12
  85. C
  86.        RP(8,9)= 12.*EIYZ*RL3*X12
  87.        RP(8,11)= 6.*EIYZ*RL2*X12
  88.        RP(9,8)= -12.*EIYZ*RL3*X12
  89.        RP(11,8)= -6.*EIYZ*RL2*X12
  90. C
  91.        RP(9,12)= 6.*EIYZ*RL2*X12
  92.        RP(12,9)= -6.*EIYZ*RL2*X12
  93. C
  94.        RP(11,12)= 4.*EIYZ*RL
  95.      .  *(1.D0+3.*X2)*X12
  96. C
  97.        RP(12,11)= -4.*EIYZ*RL
  98.      .   *(1.D0+3.*X2)*X12
  99. C
  100.       RETURN
  101.       END
  102.  
  103.  
  104.  
  105.  
  106.  

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