% visualisation des contraintes de Von Mises

[K,M,F,Q,G,H,R]=assempde(b,p,e,t,c,a,f);
u = assempde(K,M,F,Q,G,H,R);
mis = pdesmech(p,t,c,u,'tensor','von Mises');
pdeplot(p,e,t,'xydata',mis,...
 'xystyle','flat','mesh','off');
