Télécharger la plaquette Cast3M


fr en

Télécharger mkadim.eso

Retour à la liste

Numérotation des lignes : 

mkadim

  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