Télécharger la plaquette Cast3M


fr en

Télécharger poril2.eso

Retour à la liste

Numérotation des lignes : 

poril2

  1. C PORIL2    SOURCE    CB215821  26/08/24    21:17:48     12622          
  2.       SUBROUTINE PORIL2(RP,XL,YOUN,CISAIL,SD,SREDY,XINZ)
  3. C-----------------------------------------------------------------------
  4. C
  5. C    ROUTINE DE CALCUL DE LA MATRICE DE RIGIDITE ELEMENTAIRE D UNE
  6. C    POUTRE DROITE EN AXES LOCAUX EN 2D
  7. C    D APRES TEDEL
  8. C
  9. C  ENTREE
  10. C    XL      LONGUEUR
  11. C    SD      SECTION DROITE
  12. C    SREDY   SECTION REDUITE A L EFFORT TRANCHANT TY
  13. C    XINZ    INERTIE IZ
  14. C    YOUN    MODULE D YOUNG
  15. C    CISAIL  MODULE DE CISAILLEMENT
  16. C  SORTIE
  17. C    RP      MATRICE DE RIGIDITE ELEMENTAIRE 12X12
  18. C
  19. C     VIENT DE TEDEL AOUT 84
  20. C-----------------------------------------------------------------------
  21.       IMPLICIT INTEGER(I-N)
  22.       IMPLICIT REAL*8(A-H,O-Z)
  23.       DIMENSION RP(6,6)
  24. C
  25. C  INITIALISATIONS
  26. C
  27.       DO 121 J=1,6
  28.       DO 120 I=1,6
  29.       RP(I,J)=0.D00
  30.   120 CONTINUE
  31.   121 CONTINUE
  32.       RL=1.D00/XL
  33.       RL2=RL*RL
  34.       RL3=RL2*RL
  35.       X2=0.D00
  36.       X3=0.D00
  37.  
  38.            IF( SREDY. NE.0.) X3= YOUN/CISAIL*XINZ*RL2/  SREDY
  39.            ESL=YOUN*SD*RL
  40.            EIZ=YOUN*XINZ
  41.  
  42.         X12= 1./(1.+12.*X2)
  43.         X13= 1./(1.+12.*X3)
  44.  
  45. C
  46. C  REMPLISSAGE DE LA MATRICE
  47. C
  48.       RP(1,1)=ESL
  49.       RP(1,4)=-ESL
  50.        RP(2,2)= 12.*EIZ*RL3*X13
  51.       RP(2,3)=6.*EIZ*RL2*X13
  52.         RP(2,5)=  -12. *EIZ*RL3*X13
  53.       RP(2,6)= 6.*EIZ*RL2*X13
  54.       RP(3,3)= 4.*EIZ*RL
  55.      . *(1.+3.*X3)*X13
  56.       RP(3,5)=-6.*EIZ*RL2*X13
  57.       RP(3,6)=2.*EIZ*RL
  58.      .      *(1.-6.*X3)*X13
  59.       RP(4,4)=ESL
  60.        RP(5,5)= 12.*EIZ*RL3*X13
  61.       RP(5,6)=-6.*EIZ*RL2*X13
  62.       RP(6,6)= 4.*EIZ*RL
  63.      . *(1.+3.*X3)*X13
  64.       DO 122 I=1,6
  65.       DO 9 J=I,6
  66.       RP(J,I)=RP(I,J)
  67.     9 CONTINUE
  68.   122 CONTINUE
  69.       RETURN
  70.       END
  71.  
  72.  
  73.  
  74.  
  75.  
  76.  
  77.  

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