C MKADIM    SOURCE    KK2000    13/11/08    21:15:50     7860
      subroutine MKADIM(da,XMKel)
C     Matrice de rigidité élastique adimensionnée

        IMPLICIT INTEGER(I-N)
        IMPLICIT REAL*8(A-H,O-Z)
-INC CCREEL
        DIMENSION da(20), XMKel(5,5)
C        DIMENSION XMKelinv(5,5)
*
        Diam=da(1)
        XLx=da(2)
        XLy=da(3)
        XKelz=da(4)
        XKelh=da(5)
        XKelry=da(6)
        XKelrx=da(7)
        qmax=da(14)
        a9=da(15)

        do i=1,5
           do j=1,5
              XMKel(i,j)=0.
C              XMKelinv(i,j)=0.
           enddo
        enddo


        if (a9.eq.1) then
C       semelle filante
        XMKel(1,1)=XKelz/qmax
        XMKel(2,2)=XKelh/qmax
        XMKel(3,3)=XKelry/(Diam**2*qmax)
        XMKel(4,4)=XKelh/qmax
        XMKel(5,5)=XKelrx/(Diam**2*qmax)
C        XMKelinv(1,1)=1./(XKelz/qmax)
C        XMKelinv(2,2)=1./(XKelh/qmax)
C        XMKelinv(3,3)=1./(XKelry/(Diam**2*qmax))
C        XMKelinv(4,4)=1./(XKelh/qmax)
C        XMKelinv(5,5)=1./(XKelry/(Diam**2*qmax))

        else
        if (a9.eq.2) then
C       semelle rectangulaire
        XMKel(1,1)=XKelz/(qmax*sqrt(XLx**2+XLy**2))
        XMKel(2,2)=XKelh/(qmax*sqrt(XLx**2+XLy**2))
        XMKel(3,3)=XKelry/(qmax*XLx**2*sqrt(XLx**2+XLy**2))
        XMKel(4,4)=XKelh/(qmax*sqrt(XLx**2+XLy**2))
        XMKel(5,5)=XKelrx/(qmax*XLy**2*sqrt(XLx**2+XLy**2))
C        XMKelinv(1,1)=1./(XKelz/(qmax*sqrt(XLx**2+XLy**2)))
C        XMKelinv(2,2)=1./(XKelh/(qmax*sqrt(XLx**2+XLy**2)))
C        XMKelinv(3,3)=1./(XKelry/(qmax*XLx**2*sqrt(XLx**2+XLy**2)))
C        XMKelinv(4,4)=1./(XKelh/(qmax*sqrt(XLx**2+XLy**2)))
C        XMKelinv(5,5)=1./(XKelrx/(qmax*XLy**2*sqrt(XLx**2+XLy**2)))

        else
C       semelle circulaire
        XMKel(1,1)=XKelz*Diam/        (qmax*XPI*Diam**2/4)
        XMKel(2,2)=XKelh*Diam/        (qmax*XPI*Diam**2/4)
        XMKel(3,3)=XKelry/(Diam*       qmax*XPI*Diam**2/4)
        XMKel(4,4)=XKelh*Diam/        (qmax*XPI*Diam**2/4)
        XMKel(5,5)=XKelrx/(Diam*       qmax*XPI*Diam**2/4)
C        XMKelinv(1,1)=1./(XKelz*Diam/ (qmax*XPI*Diam**2/4))
C        XMKelinv(2,2)=1./(XKelh*Diam/ (qmax*XPI*Diam**2/4))
C        XMKelinv(3,3)=1./(XKelry/(Diam*qmax*XPI*Diam**2/4))
C        XMKelinv(4,4)=1./(XKelh*Diam/ (qmax*XPI*Diam**2/4))
C        XMKelinv(5,5)=1./(XKelrx/(Diam*qmax*XPI*Diam**2/4))
        endif
        endif

        return
        end