Esope source files

   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 :