Télécharger mkadim.eso

Retour à la liste

Numérotation des lignes :

  1. C MKADIM SOURCE KK2000 13/11/08 21:15:50 7860
  2. subroutine MKADIM(da,XMKel)
  3. C Matrice de rigidité élastique adimensionnée
  4.  
  5. IMPLICIT INTEGER(I-N)
  6. IMPLICIT REAL*8(A-H,O-Z)
  7. -INC CCREEL
  8. DIMENSION da(20), XMKel(5,5)
  9. C DIMENSION XMKelinv(5,5)
  10. *
  11. Diam=da(1)
  12. XLx=da(2)
  13. XLy=da(3)
  14. XKelz=da(4)
  15. XKelh=da(5)
  16. XKelry=da(6)
  17. XKelrx=da(7)
  18. qmax=da(14)
  19. a9=da(15)
  20.  
  21. do i=1,5
  22. do j=1,5
  23. XMKel(i,j)=0.
  24. C XMKelinv(i,j)=0.
  25. enddo
  26. enddo
  27.  
  28.  
  29. if (a9.eq.1) then
  30. C semelle filante
  31. XMKel(1,1)=XKelz/qmax
  32. XMKel(2,2)=XKelh/qmax
  33. XMKel(3,3)=XKelry/(Diam**2*qmax)
  34. XMKel(4,4)=XKelh/qmax
  35. XMKel(5,5)=XKelrx/(Diam**2*qmax)
  36. C XMKelinv(1,1)=1./(XKelz/qmax)
  37. C XMKelinv(2,2)=1./(XKelh/qmax)
  38. C XMKelinv(3,3)=1./(XKelry/(Diam**2*qmax))
  39. C XMKelinv(4,4)=1./(XKelh/qmax)
  40. C XMKelinv(5,5)=1./(XKelry/(Diam**2*qmax))
  41.  
  42. else
  43. if (a9.eq.2) then
  44. C semelle rectangulaire
  45. XMKel(1,1)=XKelz/(qmax*sqrt(XLx**2+XLy**2))
  46. XMKel(2,2)=XKelh/(qmax*sqrt(XLx**2+XLy**2))
  47. XMKel(3,3)=XKelry/(qmax*XLx**2*sqrt(XLx**2+XLy**2))
  48. XMKel(4,4)=XKelh/(qmax*sqrt(XLx**2+XLy**2))
  49. XMKel(5,5)=XKelrx/(qmax*XLy**2*sqrt(XLx**2+XLy**2))
  50. C XMKelinv(1,1)=1./(XKelz/(qmax*sqrt(XLx**2+XLy**2)))
  51. C XMKelinv(2,2)=1./(XKelh/(qmax*sqrt(XLx**2+XLy**2)))
  52. C XMKelinv(3,3)=1./(XKelry/(qmax*XLx**2*sqrt(XLx**2+XLy**2)))
  53. C XMKelinv(4,4)=1./(XKelh/(qmax*sqrt(XLx**2+XLy**2)))
  54. C XMKelinv(5,5)=1./(XKelrx/(qmax*XLy**2*sqrt(XLx**2+XLy**2)))
  55.  
  56. else
  57. C semelle circulaire
  58. XMKel(1,1)=XKelz*Diam/ (qmax*XPI*Diam**2/4)
  59. XMKel(2,2)=XKelh*Diam/ (qmax*XPI*Diam**2/4)
  60. XMKel(3,3)=XKelry/(Diam* qmax*XPI*Diam**2/4)
  61. XMKel(4,4)=XKelh*Diam/ (qmax*XPI*Diam**2/4)
  62. XMKel(5,5)=XKelrx/(Diam* qmax*XPI*Diam**2/4)
  63. C XMKelinv(1,1)=1./(XKelz*Diam/ (qmax*XPI*Diam**2/4))
  64. C XMKelinv(2,2)=1./(XKelh*Diam/ (qmax*XPI*Diam**2/4))
  65. C XMKelinv(3,3)=1./(XKelry/(Diam*qmax*XPI*Diam**2/4))
  66. C XMKelinv(4,4)=1./(XKelh*Diam/ (qmax*XPI*Diam**2/4))
  67. C XMKelinv(5,5)=1./(XKelrx/(Diam*qmax*XPI*Diam**2/4))
  68. endif
  69. endif
  70.  
  71. return
  72. end
  73.  
  74.  
  75.  

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