166
7 • Formulation des éléments finis
deformee_ref(idof) = vecforce(ntotdof-1)*...
(lb*abscisse_ref(idof)^2/2-abscisse_ref(idof)^3/6)/
(E*I);
end
plot(abscisse_ref,deformee_ref,abscisse,deformee)
legend(’exacte’,’approchee’,2)
xlabel(’abscisse (m)’)
ylabel(’deformee (m)’);
Comme attendu, la solution éléments finis aux nœuds de chaque élément
est exacte : la fonction d’interpolation est en effet du degré polynomial de la
déformée. La comparaison entre la solution exacte et la solution numérique
obtenue avec cinq éléments est représentée ci-dessous.
0
0.1
0.2
0.3
0.4
0.5
0.6
0.7
0.8
0.9
1
0
0.002
0.004
0.006
0.008
0.01
0.012
0.014
0.016
abscisse (m)
deformee (m)
exacte
approchee
Figure 7.1