Télécharger podef2.eso

Retour à la liste

Numérotation des lignes :

  1. C PODEF2 SOURCE PV 08/09/11 21:16:16 6150
  2. C PODEF2 SOURCE TC1 92/03/16 21:03:26 501
  3. SUBROUTINE PODEF2(V,DEFO,XL,IDEUX,xinz,sredy,you,xnu)
  4. C-----------------------------------------------------------------------
  5. C
  6. C ROUTINE DE CALCUL DES DEFORMATIONS DE LA POUTRE DROITE EN 2D
  7. C
  8. C ENTREE
  9. C V VECTEUR DES DEPLACEMENTS LOCAUX AUX 2 NOEUDS DE LA POUTRE
  10. C XL LONGUEUR
  11. C IDEUX INDICE POUR DEFORMATIONS DU 2-EME ORDRE
  12. C
  13. C SORTIE
  14. C DEFO DEFORMATIONS AUX 2 NOEUDS 6 PAR NOEUD
  15. C ( PAS D'EFFET D'EFFORT TRANCHANT || )
  16. C
  17. C-----------------------------------------------------------------------
  18. IMPLICIT INTEGER(I-N)
  19. IMPLICIT REAL*8(A-H,O-Z)
  20. DIMENSION DEFO(*),V(*)
  21. C
  22. RL=1.D00/XL
  23. RL2=RL*RL
  24. X2=0.
  25. X3=0.
  26. YOUN=YOU
  27. CISAIL=YOU/(2.D0*(1.D0+XNU))
  28.  
  29. IF( SREDY. NE.0.) X3= YOUN/CISAIL*RL2*XINZ/SREDY
  30. X13=1./(1.+12.*X3)
  31. X33=(1.+3.*X3)*X13
  32. X63=(1.-6.*X3)*X13
  33. X12=1./(1.+12.*X2)
  34. X32=(1.+3.*X2)*X12
  35. X62=(1.-6.*X2)*X12
  36. C
  37. C DEFORMATIONS GENERALISEES
  38. C
  39. EPS =(V(4)-V(1))*RL
  40. XSIZ1=
  41. .RL2*(-6.*V(2)*X13-4.*XL*X33*V(3)+6.*V(5)*X13-2.*XL*X63*V(6))
  42. XSIZ2=
  43. .RL2*( 6.*V(2)*X13+2.*XL*X63*V(3)-6.*V(5)*X13+4.*XL*X33*V(6))
  44. C
  45. C REMPLISSAGE DU TABLEAU DEFO
  46. C
  47. DEFO(1)=EPS
  48. DEFO(2)=0.D0
  49. DEFO(3)=XSIZ1
  50. DEFO(4)=EPS
  51. DEFO(5)=0.D0
  52. DEFO(6)=XSIZ2
  53. IF(IDEUX.EQ.0) RETURN
  54. C
  55. C TERMES DU SECOND ORDRE
  56. C
  57. AUX1=0.5D0*(EPS**2+V(3)**2)
  58. AUX2=0.5D0*(EPS**2+V(6)**2)
  59. DEFO(1)=DEFO(1)+0.5D0*(AUX1+AUX2)
  60. DEFO(4)=DEFO(4)+0.5D0*(AUX1+AUX2)
  61. RETURN
  62. END
  63.  
  64.  
  65.  
  66.  
  67.  

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