Télécharger pocolo.eso

Retour à la liste

Numérotation des lignes :

  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