Télécharger la plaquette Cast3M


fr en

Télécharger pocolo.eso

Retour à la liste

Numérotation des lignes : 

pocolo

  1. C POCOLO    SOURCE    CHAT      05/01/13    02:15:54     5004
  2.       SUBROUTINE POCOLO(V,CONT,XL,YOUN,CISAIL,SD,SREDY,SREDZ,
  3.      .                  XINX,XINY,XINZ,IDEUX)
  4. C-----------------------------------------------------------------------
  5. C
  6. C    ROUTINE DE CALCUL DES CONTRAINTES DE LA POUTRE DROITE
  7. C    D APRES TEDEL
  8. C
  9. C  ENTREE
  10. C    V       VECTEUR DES DEPLACEMENTS LOCAUX AUX 2 NOEUDS DE LA POUTRE
  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    IDEUX   INDICE POUR DEFORMATIONS DU SECOND ORDRE
  21. C  SORTIE
  22. C    CONT    CONTRAINTES AUX 2 NOEUDS    6 PAR NOEUD
  23. C            SEULS LES MOMENTS DE FLEXION VARIENT LINEAIREMENT
  24. C
  25. C    AOUT 84
  26. C-----------------------------------------------------------------------
  27.       IMPLICIT INTEGER(I-N)
  28.       IMPLICIT REAL*8(A-H,O-Z)
  29.       DIMENSION CONT(*),V(*)
  30. C
  31.       RL=1.D00/XL
  32.       RL2=RL*RL
  33.       RL3= RL2*RL
  34.           X2=0.
  35.           X3=0.
  36.           IF( SREDZ. NE.0.) X2= YOUN/CISAIL*RL2*XINY/SREDZ
  37.           IF( SREDY. NE.0.) X3= YOUN/CISAIL*RL2*XINZ/SREDY
  38.              X13=1./(1.+12.*X3)
  39.              X33=(1.+3.*X3)*X13
  40.              X63=(1.-6.*X3)*X13
  41.              X12=1./(1.+12.*X2)
  42.              X32=(1.+3.*X2)*X12
  43.              X62=(1.-6.*X2)*X12
  44. C
  45. C  DEFORMATIONS GENERALISEES
  46. C
  47.       EPS  =(V(7)-V(1))*RL
  48.       XSIX =(V(10)-V(4))*RL
  49.       XSIY1=
  50.      .RL2*( 6.*V(3)*X12-4.*XL*X32*V(5)-6.*V(9)*X12-2.*XL*X62*V(11))
  51.       XSIY2=
  52.      .RL2*(-6.*V(3)*X12+2.*XL*X62*V(5)+6.*V(9)*X12+4.*XL*X32*V(11))
  53.       XSIZ1=
  54.      .RL2*(-6.*V(2)*X13-4.*XL*X33*V(6)+6.*V(8)*X13-2.*XL*X63*V(12))
  55.       XSIZ2=
  56.      .RL2*( 6.*V(2)*X13+2.*XL*X63*V(6)-6.*V(8)*X13+4.*XL*X33*V(12))
  57. C
  58. C  TERMES DU SECOND ORDRE
  59. C
  60.       IF(IDEUX.EQ.0) THEN
  61.       EPS1=EPS
  62.       EPS2=EPS
  63.       ELSE IF(IDEUX.NE.0) THEN
  64.       EPS1=EPS+0.5D0*(EPS**2+V(5)**2+V(6)**2)
  65.       EPS2=EPS+0.5D0*(EPS**2+V(11)**2+V(12)**2)
  66.       ENDIF
  67. C
  68. C  CONTRAINTES GENERALISEES
  69. C
  70.       SIG1  =YOUN*SD*(EPS1+EPS2)*0.5D0
  71.       SIG2  =SIG1
  72.       TOR   =CISAIL*XINX*XSIX
  73.       XMOMZ1=YOUN*XINZ*XSIZ1
  74.       XMOMZ2=YOUN*XINZ*XSIZ2
  75.       XMOMY1=YOUN*XINY*XSIY1
  76.       XMOMY2=YOUN*XINY*XSIY2
  77.       TRANY =(XMOMZ1-XMOMZ2)*RL
  78.       TRANZ =(XMOMY2-XMOMY1)*RL
  79. C
  80. C  REMPLISSAGE DU TABLEAU CONT
  81. C
  82.       CONT(1)=SIG1
  83.       CONT(2)=TRANY
  84.       CONT(3)=TRANZ
  85.       CONT(4)=TOR
  86.       CONT(5)=XMOMY1
  87.       CONT(6)=XMOMZ1
  88.       CONT(7)=SIG2
  89.       CONT(8)=TRANY
  90.       CONT(9)=TRANZ
  91.       CONT(10)=TOR
  92.       CONT(11)=XMOMY2
  93.       CONT(12)=XMOMZ2
  94.       RETURN
  95.       END
  96.  
  97.  

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