Télécharger la plaquette Cast3M


fr en

Télécharger pocol2.eso

Retour à la liste

Numérotation des lignes : 

pocol2

  1. C POCOL2    SOURCE    CHAT      05/01/13    02:15:51     5004
  2. C POCOL2    SOURCE    TC1       92/03/16    21:04:19     501
  3.       SUBROUTINE POCOL2(V,CONT,XL,YOUN,CISAIL,SD,SREDY,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    XINZ    INERTIE IZ
  15. C    YOUN    MODULE D YOUNG
  16. C    CISAIL  MODULE DE CISAILLEMENT
  17. C    IDEUX   INDICE POUR DEFORMATIONS DU SECOND ORDRE
  18. C  SORTIE
  19. C    CONT    CONTRAINTES AUX 2 NOEUDS    6 PAR NOEUD
  20. C            SEULS LES MOMENTS DE FLEXION VARIENT LINEAIREMENT
  21. C
  22. C    AOUT 84
  23. C-----------------------------------------------------------------------
  24.       IMPLICIT INTEGER(I-N)
  25.       IMPLICIT REAL*8(A-H,O-Z)
  26.       DIMENSION CONT(*),V(*)
  27. C
  28.       RL=1.D00/XL
  29.       RL2=RL*RL
  30.       RL3= RL2*RL
  31.           X2=0.
  32.           X3=0.
  33.           IF( SREDY. NE.0.) X3= YOUN/CISAIL*RL2*XINZ/SREDY
  34.              X13=1./(1.+12.*X3)
  35.              X33=(1.+3.*X3)*X13
  36.              X63=(1.-6.*X3)*X13
  37.              X12=1./(1.+12.*X2)
  38.              X32=(1.+3.*X2)*X12
  39.              X62=(1.-6.*X2)*X12
  40. C
  41. C  DEFORMATIONS GENERALISEES
  42. C
  43.       EPS  =(V(4)-V(1))*RL
  44.       XSIZ1=
  45.      .RL2*(-6.*V(2)*X13-4.*XL*X33*V(3)+6.*V(5)*X13-2.*XL*X63*V(6))
  46.       XSIZ2=
  47.      .RL2*( 6.*V(2)*X13+2.*XL*X63*V(3)-6.*V(5)*X13+4.*XL*X33*V(6))
  48. C
  49. C  TERMES DU SECOND ORDRE
  50. C
  51.       IF(IDEUX.EQ.0) THEN
  52.       EPS1=EPS
  53.       EPS2=EPS
  54.       ELSE IF(IDEUX.NE.0) THEN
  55.       EPS1=EPS+0.5D0*(EPS**2+V(3)**2)
  56.       EPS2=EPS+0.5D0*(EPS**2+V(6)**2)
  57.       ENDIF
  58. C
  59. C  CONTRAINTES GENERALISEES
  60. C
  61.       SIG1  =YOUN*SD*(EPS1+EPS2)*0.5D0
  62.       SIG2  =SIG1
  63.       XMOMZ1=YOUN*XINZ*XSIZ1
  64.       XMOMZ2=YOUN*XINZ*XSIZ2
  65.       TRANY =(XMOMZ1-XMOMZ2)*RL
  66. C
  67. C  REMPLISSAGE DU TABLEAU CONT
  68. C
  69.       CONT(1)=SIG1
  70.       CONT(2)=TRANY
  71.       CONT(3)=XMOMZ1
  72.       CONT(4)=SIG2
  73.       CONT(5)=TRANY
  74.       CONT(6)=XMOMZ2
  75.       RETURN
  76.       END
  77.  
  78.  
  79.  

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