df1d poutre.m

１Ｄ力学模型

%******************************************************%
%          Probleme 1D: Elasticite 1D                  %
%          METHODE des DIFFERENCES FINIES              %
%             PROBLEME ELLIPTIQUE                      %
%******************************************************%


clear all
close all
format long

%0) Constantes du probleme
%======================================================================
L=1;
E=1;
f=1;
A=1;
P=1;

%1) Creation du maillage
%=======================================================================

nx=input('entrez le nombre de noeuds suivant x');

[Nx,h]=Maillage1D(L,nx);

%2) Construction de la matrice de Raideur K et du seconde membre F %
%======================================================================

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

%4)Resolution du probleme
%=====================================================================

Uresult=K\F;

%5) PostTraitement
%====================================================================
PostDF1D(Uresult,Nx,f,E,nx,L,A,P);