rbconvinit_stressfree

 init, different BC, lam and lammax compared to rbconvinit

0001 function p=rbconvinit_stressfree(p) 
0002 % init, different BC, lam and lammax compared to rbconvinit
0003 pre=sprintf('%s',inputname(1)); p=setfn(p,pre);
0004 p=stanparam(p); p.pstyle=2; 
0005 p.neq=3; p.f=@rbconvf; p.jac=@rbconvjac; 
0006 p.bcf=@rbconvbc_stressfree; p.outfu=@stanbra; 
0007 lx=2; lz=0.5; p.geo=rec(lx,lz); p=stanmesh(p,100,25); p=setbmesh(p);
0008 p.ds=0.1; p.dsmin=0.00001; p.dlammax=10; p.nsteps=5; p.neig=20; 
0009 p.jsw=0; p.nsw=0; p.bisecmax=5; p.amod=0; p.pcomp=2; p.tol=1e-10;
0010 p.lam=664; p.u=0*ones(p.neq*p.np,1); p.lammax=900;
0011 p.xi=1/p.np; p.tau=zeros(p.neq*p.np+1,1); p.tau(1)=1; 
0012