Télécharger podefo.eso

Retour à la liste

Numérotation des lignes :

  1. C PODEFO SOURCE CHAT 11/03/22 21:15:13 6910
  2. SUBROUTINE PODEFO(V,DEFO,XL,IDEUX,carac,you,xnu)
  3. C-----------------------------------------------------------------------
  4. C
  5. C ROUTINE DE CALCUL DES DEFORMATIONS DE LA POUTRE DROITE
  6. C
  7. C ENTREE
  8. C V VECTEUR DES DEPLACEMENTS LOCAUX AUX 2 NOEUDS DE LA POUTRE
  9. C XL LONGUEUR
  10. C IDEUX INDICE POUR DEFORMATIONS DU 2-EME ORDRE
  11. C
  12. C SORTIE
  13. C DEFO DEFORMATIONS AUX 2 NOEUDS 6 PAR NOEUD
  14. C ( PAS D'EFFET D'EFFORT TRANCHANT || )
  15. C
  16. C-----------------------------------------------------------------------
  17. IMPLICIT INTEGER(I-N)
  18. IMPLICIT REAL*8(A-H,O-Z)
  19. DIMENSION DEFO(*),V(*),carac(*)
  20. C
  21. RL=1.D00/XL
  22. RL2=RL*RL
  23. X2=0.
  24. X3=0.
  25. xiny=carac(2)
  26. xinz=carac(3)
  27. sredy=carac(5)
  28. sredz=carac(6)
  29. cisail= you/(2.*(1+xnu))
  30. * of
  31. * if( sredy.ne.0.) X2=You/cisail*rl2*xiny/sredy
  32. * if( sredy.ne.0.) X2=You/cisail*rl2*xinz/sredz
  33. if( sredy.ne.0.) X3=You/cisail*rl2*xinz/sredy
  34. if( sredz.ne.0.) X2=You/cisail*rl2*xiny/sredz
  35. * of
  36. X13=1./(1.+12.*X3)
  37. X33=(1.+3.*X3)*X13
  38. X63=(1.-6.*X3)*X13
  39. X12=1./(1.+12.*X2)
  40. X32=(1.+3.*X2)*X12
  41. X62=(1.-6.*X2)*X12
  42. C
  43. C DEFORMATIONS GENERALISEES
  44. C
  45. EPS =(V(7)-V(1))*RL
  46. XSIX =(V(10)-V(4))*RL
  47. XSIY1=
  48. .RL2*( 6.*V(3)*X12-4.*XL*X32*V(5)-6.*V(9)*X12-2.*XL*X62*V(11))
  49. XSIY2=
  50. .RL2*(-6.*V(3)*X12+2.*XL*X62*V(5)+6.*V(9)*X12+4.*XL*X32*V(11))
  51. XSIZ1=
  52. .RL2*(-6.*V(2)*X13-4.*XL*X33*V(6)+6.*V(8)*X13-2.*XL*X63*V(12))
  53. XSIZ2=
  54. .RL2*( 6.*V(2)*X13+2.*XL*X63*V(6)-6.*V(8)*X13+4.*XL*X33*V(12))
  55. C
  56. C REMPLISSAGE DU TABLEAU DEFO
  57. C
  58. DEFO(1)=EPS
  59. DEFO(2)=0.D0
  60. DEFO(3)=0.D0
  61. DEFO(4)=XSIX
  62. DEFO(5)=XSIY1
  63. DEFO(6)=XSIZ1
  64. DEFO(7)=EPS
  65. DEFO(8)=0.D0
  66. DEFO(9)=0.D0
  67. DEFO(10)=XSIX
  68. DEFO(11)=XSIY2
  69. DEFO(12)=XSIZ2
  70. IF(IDEUX.EQ.0) RETURN
  71. C
  72. C TERMES DU SECOND ORDRE
  73. C
  74. AUX1=0.5D0*(EPS**2+V(5)**2+V(6)**2)
  75. AUX2=0.5D0*(EPS**2+V(11)**2+V(12)**2)
  76. DEFO(1)=DEFO(1)+0.5D0*(AUX1+AUX2)
  77. DEFO(7)=DEFO(7)+0.5D0*(AUX1+AUX2)
  78. C*OF? DEFO(1)=DEFO(1)+AUX1
  79. C*OF? DEFO(7)=DEFO(7)+AUX2
  80. RETURN
  81. END
  82.  
  83.  
  84.  
  85.  

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