vkdemo

 script for von Karman plate

0001 % script for von Karman plate
0002 close all; clear all; 
0003 mypause('Init p for trivial branch of von Karman plate and run cont'); 
0004 p=[]; p=vkinit(p); p=cont(p); % the trivial solution
0005 mypause('2 primary and one secondary bifurcations'); 
0006 q=swibra('p','bp1','q'); q.dlammax=1; q=cont(q);
0007 r=swibra('p','bp2','r'); r.dlammax=1; r=cont(r);
0008 w=swibra('r','bp1','w'); w.dlammax=1; w.smod=3;w=cont(w);
0009 mypause('Plot bifurcation diagram'); 
0010 clf(3); cmp=2;
0011 plotbra(q,3,cmp,'ms',10,'lwst',2,'lwun',2,'cl','k');
0012 plotbra(r,3,cmp,'ms',10,'lwst',2,'lwun',2,'cl','b');
0013 plotbra(w,3,cmp,'ms',10,'lab',[3,6],'lw',5,'cl','r');
0014 axis([5 8 0 8]);xlabel('\lambda'); ylabel('||u_1||_2');
0015 mypause('Plot some solutions'); 
0016 plotsolf('q','bp1',4,1,2);
0017 plotsolf('r','bp1',5,1,2);
0018 plotsolf('w','p3',6,1,2);
0019 plotsolf('w','p6',7,1,2);
0020 plotsolf('q','bp1',8,3,2);
0021 plotsolf('r','bp1',9,3,2);