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