Sources Esope | Cast3M

mkredi
C MKREDI    SOURCE    AM        11/05/26    21:16:01     6982
      subroutine MKREDI(da,XMKel)
C     Matrice de rigidité élastique redimensionnée
 
        IMPLICIT INTEGER(I-N)
        IMPLICIT REAL*8(A-H,O-Z)
        DIMENSION da(20), XMKel(5,5), XMKelinv(5,5)
*
        XKelz=da(4)
        XKelh=da(5)
        XKelry=da(6)
        XKelrx=da(7)
*
        do i=1,5
           do j=1,5
              XMKel(i,j)=0.
              XMKelinv(i,j)=0.
           enddo
        enddo
*
        XMKel(1,1)=XKelz
        XMKel(2,2)=XKelh
        XMKel(3,3)=XKelry
        XMKel(4,4)=XKelh
        XMKel(5,5)=XKelrx
        XMKelinv(1,1)=1./(XKelz)
        XMKelinv(2,2)=1./(XKelh)
        XMKelinv(3,3)=1./(XKelry)
        XMKelinv(4,4)=1./(XKelh)
        XMKelinv(5,5)=1./(XKelrx)
*
        return
        end