Télécharger la plaquette Cast3M


fr en

Télécharger poksl2.eso

Retour à la liste

Numérotation des lignes : 

poksl2

  1. C POKSL2    SOURCE    CB215821  26/08/24    21:17:43     12622          
  2. C POKSL2    SOURCE    CASTEM    85/06/13    12:33:17
  3.       SUBROUTINE POKSL2(RP,XL,SIGM)
  4. C-----------------------------------------------------------------------
  5. C
  6. C    ROUTINE DE CALCUL DE LA MATRICE DE RIGIDITE GEOMETRIQUE
  7. C    ELEMENTAIRE D UNE POUTRE DROITE EN AXES LOCAUX
  8. C    D APRES TEDEL
  9. C
  10. C  ENTREE
  11. C    XL      LONGUEUR
  12. C    SIGM    EFFORT NORMAL MOYEN
  13. C  SORTIE
  14. C    RP      MATRICE DE RIGIDITE GEOMETRIQUE ELEMENTAIRE 12X12
  15. C
  16. C     MAI 85 MILLARD
  17. C-----------------------------------------------------------------------
  18.       IMPLICIT INTEGER(I-N)
  19.       IMPLICIT REAL*8(A-H,O-Z)
  20.       DIMENSION RP(6,6)
  21. C
  22. C  INITIALISATIONS
  23. C
  24.       DO 121 J=1,6
  25.       DO 120 I=1,6
  26.       RP(I,J)=0.D00
  27.   120 CONTINUE
  28.   121 CONTINUE
  29.       RL=1.D00/XL
  30.       RN=SIGM*RL
  31.       U1=1.2D00*RN
  32.       U2=0.1D00*SIGM
  33.       U3=0.13333333333333333D00*SIGM*XL
  34.       U4=0.03333333333333333D00*SIGM*XL
  35.       RP(1 ,1)=RN
  36.       RP(1,4)=-RN
  37.       RP(4,4)=RN
  38.       RP(2,2)=U1
  39.       RP(2,3)=U2
  40.       RP(2,5)=-U1
  41.       RP(2,6)=U2
  42.       RP(3,3)=U3
  43.       RP(3,5)=-U2
  44.       RP(3,6)=-U4
  45.       RP(5,5)=U1
  46.       RP(5,6)=-U2
  47.       RP(6,6)=U3
  48. C
  49. C  REMPLISSAGE DE LA MATRICE
  50. C
  51.       DO 122 I=1,6
  52.       DO 9 J=I,6
  53.       RP(J,I)=RP(I,J)
  54.     9 CONTINUE
  55.   122 CONTINUE
  56.       RETURN
  57.       END
  58.  
  59.  
  60.  
  61.  
  62.  

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