7.4 Application à l’étude d’une poutre sollicitée en flexion
165
c = ddlfixe(idof);
for jdof=1:size(kglobal)
kglobal(c,jdof) = 0. ;
end
kglobal(c,c) = 1. ;
vecforce(c) = valfixe(idof) ;
end
% Inversion du systeme
vecdisp = kglobal\vecforce ;
% Deplacement calcule
disp([’Deformee max calculee :
’,num2str(vecdisp(ntotdof-1))])
% Deplacement theorique
disp([’Deformee max theorique :
’,num2str(vecforce(ntotdof-1)*lb^3/(3.*E*I))])
% Construction de la solution approchee
abscisse = zeros(nel+1,1) ; deformee = zeros(nel+1,1) ;
i = 0 ;
for idof=1:ndof:ntotdof
i = i+1 ;
abscisse(i) = (i-1)*l ;
deformee(i) = vecdisp(idof) ;
end
% Construction de la solution exacte
npref = 100 ;
abscisse_ref = zeros(npref,1);
deformee_ref = zeros(npref,1);
for idof=1:npref
abscisse_ref(idof) = lb/(npref-1)*(idof-1);
© Dunod – La photocopie non autorisée est un délit