stress free BC
0001 function bc=rbconvbc_stressfree(p,u,lam) 0002 % stress free BC 0003 stifffac=10^4; % stiff spring factor 0004 qd = stifffac*[[1;0;0] [0;1;0] [0;0;1]];% Dirichlet: 0005 gv = zeros(p.neq,1);% values 0006 qm = qd; qm(3,3)=0; % mixed: Dirichlet for psi and om, and T_x=0 0007 0008 % ordering: bottom, right, top, left 0009 bc = gnbc(p.neq,qd,gv,qm,gv,qd,gv,qm,gv);