Télécharger la plaquette Cast3M


fr en

Télécharger porim2.eso

Retour à la liste

Numérotation des lignes : 

porim2

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

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