postdf2d.asv

２Ｄ力学模型.有限单元法.巴黎六大实验教材.

%*****************************************************************************
%  PostDF2D.m




function []=PostDF2D(Tresult,Nx,Ny,nx,ny,L,H);

%1)Affichage de la solution calculee
%=======================================================================

figure
plot(Nx,Ny,Tresult,'ko');

%2) Affichage de la solution exacte
%=========================================================================

for i=1:nx
    Uex(i)=(-f*Nx(i)^2)/(2*E*A)+((P+f*L)/(E*A))*Nx(i);
    
end

hold on
plot(Nx,Uex,'k-');

legend('u^h','u^(ex)')

xlabel('x','Fontsize',12);
ylabel('u^h(x)','Fontsize',12);


%3=Calcul de l'erreur en norme L^(inf)
%=========================================================================

Err=Uex-Uresult';
err=max(abs(Err))