Télécharger la plaquette Cast3M


fr en

Télécharger porilo.eso

Retour à la liste

Numérotation des lignes : 

porilo

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

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