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