matriceraideurdf1d.m

１Ｄ力学模型

%**********************************************************
% MatriceRaideurDF1D.m
%
%Creation de la matrice de raideur par discretisation
%par differences finies. Les conditions aux limites
%pour le probleme 1 sont egalement imposees.
%**********************************************************


function [K,F]=MatriceRaideurDF1D(h,nx,f,E,Nx,A,P)

%1) remplissage de la matrice avec des zeros.

K=sparse([]);
K(nx,nx)=0;

for i=1:nx
    if i==1;
        K(i,i)=1
    elseif i==nx
        K(i,i-1)=-1;
        K(i,i)=1;
    else
        K(i,i-1)=1;
        K(i,i)=-2;
        K(i,i+1)=1;
    end
end

%2) Creation du vecteur second membre

F=(-f*h^2/(E*A))*ones(nx,1);
F(1)=0;
F(nx)=P*h/(E*A);