�� sr !com.femlab.server.ModelFileHeader�D���% L tagst Lcom/femlab/util/i;L typesq ~ L vrsnt Lcom/femlab/util/FlVersion;xpw sr com.femlab.util.FlVersion��%�/B = I buildI majorL datet Ljava/lang/String;L extq ~ L nameq ~ L rcsq ~ L reactionExtq ~ L reactionNameq ~ L scriptExtq ~ L scriptNameq ~ xpw t COMSOL Script 1.1t t #COMSOL Reaction Engineering Lab 1.2q ~ t COMSOL 3.3q ~ w �t $Name: $t $Date: 2006/08/31 18:03:47 $xur [Ljava.lang.String;��V��{G xp t modelinfot xfemt guit fem0t g7t g8t fem1t fem1.0q ~ q ~ t fem44t fem44.0q ~ q ~ t mfileinfouq ~ q ~ t femstructt guistructq ~ t drawq ~ t geomt mesht solutiont xmeshq ~ q ~ q ~ q ~ !q ~ xsr com.femlab.api.client.ModelInfo�^���% L descrq ~ L docURLq ~ [ imaget [Bxpw pt pxuq ~ t�clear xfem clear vrsn vrsn.name = 'COMSOL 3.3'; vrsn.ext = ''; vrsn.major = 0; vrsn.build = 405; vrsn.rcs = '$Name: $'; vrsn.date = '$Date: 2006/08/31 18:03:47 $'; xfem.version = vrsn; xfem.id = 43; xfem.geomdata = 'geom'; xfem.eqvars = 'on'; xfem.cplbndeq = 'on'; xfem.cplbndsh = 'off'; xfem.drawvalid = 'on'; xfem.geomvalid = 'on'; xfem.solvalid = 'on'; xfem.meshtime = 'TIME'; clear appl appl.mode.class = 'SmeSolid3'; appl.mode.type = 'cartesian'; appl.dim = {'u','v','w','p'}; appl.sdim = {'X','Y','Z'}; appl.name = 'smsld'; appl.module = 'SME'; appl.shape = {'shlag(2,''u'')','shlag(2,''v'')','shlag(2,''w'')'}; appl.gporder = 4; appl.cporder = 2; appl.sshape = 2; appl.border = 'off'; appl.assignsuffix = '_smsld'; appl.assign = {'Ex','EX_smsld','Ey','EY_smsld','Ez','EZ_smsld','Fx','FX_smsld', ... 'FxAmp','FXAmp_smsld','FxPh','FXPh_smsld','Fxg','FXg_smsld','Fy','FY_smsld', ... 'FyAmp','FYAmp_smsld','FyPh','FYPh_smsld','Fyg','FYg_smsld','Fz','FZ_smsld', ... 'FzAmp','FZAmp_smsld','FzPh','FZPh_smsld','Fzg','FZg_smsld','Gxy','GXY_smsld', ... 'Gxz','GXZ_smsld','Gyz','GYZ_smsld','Tax','TaX_smsld','Tay','TaY_smsld', ... 'Taz','TaZ_smsld','alphax','alphaX_smsld','alphay','alphaY_smsld', ... 'alphaz','alphaZ_smsld','ex','eX_smsld','exi','eXi_smsld','exy','eXY_smsld', ... 'exyi','eXYi_smsld','exz','eXZ_smsld','exzi','eXZi_smsld','ey','eY_smsld', ... 'eyi','eYi_smsld','eyz','eYZ_smsld','eyzi','eYZi_smsld','ez','eZ_smsld', ... 'ezi','eZi_smsld','nuxy','nuXY_smsld','nuxz','nuXZ_smsld','nuyz','nuYZ_smsld', ... 'sx','sX_smsld','sxi','sXi_smsld','sxy','sXY_smsld','sxyi','sXYi_smsld', ... 'sxz','sXZ_smsld','sxzi','sXZi_smsld','sy','sY_smsld','syi','sYi_smsld', ... 'syz','sYZ_smsld','syzi','sYZi_smsld','sz','sZ_smsld','szi','sZi_smsld'}; clear prop prop.elemdefault='Lag2'; prop.analysis='static'; prop.largedef='on'; prop.impl='weak'; prop.eigtype='freq'; prop.createframe='off'; prop.deformframe='ref'; prop.frame='ref'; clear weakconstr weakconstr.value = 'off'; weakconstr.dim = {'lm1','lm2','lm3'}; prop.weakconstr = weakconstr; appl.prop = prop; clear pnt pnt.Fx = {'0'}; pnt.FxAmp = {'1'}; pnt.FxPh = {'0'}; pnt.Fy = {'0'}; pnt.FyAmp = {'1'}; pnt.FyPh = {'0'}; pnt.Fz = {'0'}; pnt.FzAmp = {'1'}; pnt.FzPh = {'0'}; pnt.loadcoord = {'global'}; pnt.constrcoord = {'global'}; pnt.constrtype = {'standard'}; pnt.H = {{'0','0','0';'0','0','0';'0','0','0'}}; pnt.R = {{'0';'0';'0'}}; pnt.Rx = {'0'}; pnt.Hx = {'0'}; pnt.Ry = {'0'}; pnt.Hy = {'0'}; pnt.Rz = {'0'}; pnt.Hz = {'0'}; pnt.style = {{{'0'},{'0','0','0'}}}; pnt.ind = [1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1]; appl.pnt = pnt; clear edg edg.Fx = {'0'}; edg.FxAmp = {'1'}; edg.FxPh = {'0'}; edg.Fy = {'0'}; edg.FyAmp = {'1'}; edg.FyPh = {'0'}; edg.Fz = {'0'}; edg.FzAmp = {'1'}; edg.FzPh = {'0'}; edg.loadcoord = {'global'}; edg.constrcoord = {'global'}; edg.constrtype = {'standard'}; edg.H = {{'0','0','0';'0','0','0';'0','0','0'}}; edg.R = {{'0';'0';'0'}}; edg.Rx = {'0'}; edg.Hx = {'0'}; edg.Ry = {'0'}; edg.Hy = {'0'}; edg.Rz = {'0'}; edg.Hz = {'0'}; edg.style = {{{'0'},{'0','0','0'},{'solid'}}}; edg.ind = [1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1, ... 1,1,1,1,1,1]; appl.edg = edg; clear bnd bnd.name = {'Fluid load','',''}; bnd.Fx = {'-lm7*dvol_ale/dvol','0','0'}; bnd.FxAmp = {'1','1','1'}; bnd.FxPh = {'0','0','0'}; bnd.Fy = {'-lm8*dvol_ale/dvol','0','0'}; bnd.FyAmp = {'1','1','1'}; bnd.FyPh = {'0','0','0'}; bnd.Fz = {'-lm9*dvol_ale/dvol','0','0'}; bnd.FzAmp = {'1','1','1'}; bnd.FzPh = {'0','0','0'}; bnd.loadcoord = {'global','global','global'}; bnd.constrcoord = {'global','global','global'}; bnd.constrtype = {'standard','standard','standard'}; bnd.constrcond = {'free','free','fixed'}; bnd.velx = {'0','0','0'}; bnd.Hvelx = {'0','0','0'}; bnd.accx = {'0','0','0'}; bnd.Haccx = {'0','0','0'}; bnd.vely = {'0','0','0'}; bnd.Hvely = {'0','0','0'}; bnd.accy = {'0','0','0'}; bnd.Haccy = {'0','0','0'}; bnd.velz = {'0','0','0'}; bnd.Hvelz = {'0','0','0'}; bnd.accz = {'0','0','0'}; bnd.Haccz = {'0','0','0'}; bnd.H = {{'0','0','0';'0','0','0';'0','0','0'},{'0','0','0';'0','0','0'; ... '0','0','0'},{'0','0','0';'0','0','0';'0','0','0'}}; bnd.R = {{'0';'0';'0'},{'0';'0';'0'},{'0';'0';'0'}}; bnd.Rx = {'0','0','0'}; bnd.Hx = {'0','0','0'}; bnd.Ry = {'0','0','0'}; bnd.Hy = {'0','0','0'}; bnd.Rz = {'0','0','0'}; bnd.Hz = {'0','0','0'}; bnd.offset = {'0','0','0'}; bnd.pn = {'E_smsld/h*min(1e-4*5^auglagiter,0.1)','E_smsld/h*min(1e-4*5^auglagiter,0.1)', ... 'E_smsld/h*min(1e-4*5^auglagiter,0.1)'}; bnd.pt = {'E_smsld/h*min(1e-4*5^auglagiter,0.1)','E_smsld/h*min(1e-4*5^auglagiter,0.1)', ... 'E_smsld/h*min(1e-4*5^auglagiter,0.1)'}; bnd.frictiontype = {'nofriction','nofriction','nofriction'}; bnd.searchmethod = {'fast','fast','fast'}; bnd.searchdist = {'auto','auto','auto'}; bnd.mandist = {'1e-2','1e-2','1e-2'}; bnd.contacttol = {'auto','auto','auto'}; bnd.mantol = {'1e-6','1e-6','1e-6'}; bnd.dynfric = {'0','0','0'}; bnd.mustat = {'0','0','0'}; bnd.mudyn = {'0','0','0'}; bnd.dcfric = {'0','0','0'}; bnd.cohe = {'0','0','0'}; bnd.Ttmax = {'Inf','Inf','Inf'}; bnd.Tni = {'1e6','1e6','1e6'}; bnd.contact_oldi = {'0','0','0'}; bnd.Ttxi = {'1e6','1e6','1e6'}; bnd.xm_oldi = {'0','0','0'}; bnd.Ttyi = {'1e6','1e6','1e6'}; bnd.ym_oldi = {'0','0','0'}; bnd.Ttzi = {'1e6','1e6','1e6'}; bnd.zm_oldi = {'0','0','0'}; bnd.style = {{{'0'},{'217','217','217'}},{{'0'},{'0','0','255'}},{{'0'},{'0', ... '255','0'}}}; bnd.ind = [2,2,2,2,2,2,2,2,2,2,2,1,1,3,3,1,1,2]; appl.bnd = bnd; clear equ equ.shape = {[1;2;3],[1;2;3]}; equ.gporder = {{1;1;1;1},{1;1;1;1}}; equ.cporder = {{1;1;1;1},{1;1;1;1}}; equ.init = {{'0';'0';'0';'0'},{'0';'0';'0';'0'}}; equ.usage = {0,1}; equ.name = {'Fluid domain','Solid domain'}; equ.E = {'2.0e11','31100'}; equ.rho = {'7850','7850'}; equ.dampingtype = {'Rayleigh','Rayleigh'}; equ.alphadM = {'1','1'}; equ.betadK = {'0.001','0.001'}; equ.eta = {'0','0'}; equ.dinit = {{'0';'0';'0';'0'},{'0';'0';'0';'0'}}; equ.materialcoord = {'global','global'}; equ.materialmodel = {'iso','iso'}; equ.mixedform = {'0','0'}; equ.hypertype = {'neo_hookean','neo_hookean'}; equ.mu = {'8e5','8e5'}; equ.C10 = {'2e5','2e5'}; equ.C01 = {'2e5','2e5'}; equ.kappa = {'1e10','1e10'}; equ.hardeningmodel = {'iso','iso'}; equ.yieldtype = {'mises','mises'}; equ.isodata = {'tangent','tangent'}; equ.ETkin = {'2.0e10','2.0e10'}; equ.ETiso = {'2.0e10','2.0e10'}; equ.Sys = {'2.0e8','2.0e8'}; equ.Syfunc = {'mises_smsld','mises_smsld'}; equ.Syfunc_kin = {'misesKin_smsld','misesKin_smsld'}; equ.Shard = {'2.0e10/(1-2.0e10/2.0e11)*epe_smsld','2.0e10/(1-2.0e10/2.0e11)*epe_smsld'}; equ.ini_stress = {'0','0'}; equ.ini_strain = {'0','0'}; equ.alpha = {'1.2e-5','1.2e-5'}; equ.nu = {'0.33','.49'}; equ.Ex = {'2.0e11','2.0e11'}; equ.nuxy = {'0.33','0.33'}; equ.Gxy = {'7.52e10','7.52e10'}; equ.alphax = {'1.2e-5','1.2e-5'}; equ.sxi = {'0','0'}; equ.sxyi = {'0','0'}; equ.exi = {'0','0'}; equ.exyi = {'0','0'}; equ.Ey = {'2.0e11','2.0e11'}; equ.nuyz = {'0.33','0.33'}; equ.Gyz = {'7.52e10','7.52e10'}; equ.alphay = {'1.2e-5','1.2e-5'}; equ.syi = {'0','0'}; equ.syzi = {'0','0'}; equ.eyi = {'0','0'}; equ.eyzi = {'0','0'}; equ.Ez = {'2.0e11','2.0e11'}; equ.nuxz = {'0.33','0.33'}; equ.Gxz = {'7.52e10','7.52e10'}; equ.alphaz = {'1.2e-5','1.2e-5'}; equ.szi = {'0','0'}; equ.sxzi = {'0','0'}; equ.ezi = {'0','0'}; equ.exzi = {'0','0'}; equ.D = {{'2.0e11/((1+0.33)*(1-2*0.33))*(1-0.33)','2.0e11/((1+0.33)*(1-2*0.33))*0.33', ... '2.0e11/((1+0.33)*(1-2*0.33))*0.33','0','0','0';'2.0e11/((1+0.33)*(1-2*0.33))*0.33', ... '2.0e11/((1+0.33)*(1-2*0.33))*(1-0.33)','2.0e11/((1+0.33)*(1-2*0.33))*0.33', ... '0','0','0';'2.0e11/((1+0.33)*(1-2*0.33))*0.33','2.0e11/((1+0.33)*(1-2*0.33))*0.33', ... '2.0e11/((1+0.33)*(1-2*0.33))*(1-0.33)','0','0','0';'0','0','0','2.0e11/((1+0.33)*2)', ... '0','0';'0','0','0','0','2.0e11/((1+0.33)*2)','0';'0','0','0','0','0','2.0e11/((1+0.33)*2)'}, ... {'2.0e11/((1+0.33)*(1-2*0.33))*(1-0.33)','2.0e11/((1+0.33)*(1-2*0.33))*0.33', ... '2.0e11/((1+0.33)*(1-2*0.33))*0.33','0','0','0';'2.0e11/((1+0.33)*(1-2*0.33))*0.33', ... '2.0e11/((1+0.33)*(1-2*0.33))*(1-0.33)','2.0e11/((1+0.33)*(1-2*0.33))*0.33', ... '0','0','0';'2.0e11/((1+0.33)*(1-2*0.33))*0.33','2.0e11/((1+0.33)*(1-2*0.33))*0.33', ... '2.0e11/((1+0.33)*(1-2*0.33))*(1-0.33)','0','0','0';'0','0','0','2.0e11/((1+0.33)*2)', ... '0','0';'0','0','0','0','2.0e11/((1+0.33)*2)','0';'0','0','0','0','0','2.0e11/((1+0.33)*2)'}}; equ.alphavector = {{'1.2e-5';'1.2e-5';'1.2e-5';'0';'0';'0'},{'1.2e-5';'1.2e-5';'1.2e-5'; ... '0';'0';'0'}}; equ.Fx = {'0','0'}; equ.FxAmp = {'1','1'}; equ.FxPh = {'0','0'}; equ.Fy = {'0','0'}; equ.FyAmp = {'1','1'}; equ.FyPh = {'0','0'}; equ.Fz = {'0','0'}; equ.FzAmp = {'1','1'}; equ.FzPh = {'0','0'}; equ.loadcoord = {'global','global'}; equ.Tflag = {'0','0'}; equ.Temp = {'0','0'}; equ.Tempref = {'0','0'}; equ.constrcoord = {'global','global'}; equ.constrtype = {'standard','standard'}; equ.constrcond = {'free','free'}; equ.velx = {'0','0'}; equ.Hvelx = {'0','0'}; equ.accx = {'0','0'}; equ.Haccx = {'0','0'}; equ.vely = {'0','0'}; equ.Hvely = {'0','0'}; equ.accy = {'0','0'}; equ.Haccy = {'0','0'}; equ.velz = {'0','0'}; equ.Hvelz = {'0','0'}; equ.accz = {'0','0'}; equ.Haccz = {'0','0'}; equ.H = {{'0','0','0';'0','0','0';'0','0','0'},{'0','0','0';'0','0','0'; ... '0','0','0'}}; equ.R = {{'0';'0';'0'},{'0';'0';'0'}}; equ.Rx = {'0','0'}; equ.Hx = {'0','0'}; equ.Ry = {'0','0'}; equ.Hy = {'0','0'}; equ.Rz = {'0','0'}; equ.Hz = {'0','0'}; equ.dx = {'1','1'}; equ.Lx = {'dx_smsld','dx_smsld'}; equ.dy = {'1','1'}; equ.Ly = {'dy_smsld','dy_smsld'}; equ.dz = {'1','1'}; equ.Lz = {'dz_smsld','dz_smsld'}; equ.dr = {'1','1'}; equ.Lr = {'dr_smsld','dr_smsld'}; equ.R0 = {'1','1'}; equ.userPML = {{'0';'0';'0'},{'0';'0';'0'}}; equ.PMLtype = {'none','none'}; equ.coordOnx = {'0','0'}; equ.coordOny = {'0','0'}; equ.coordOnz = {'0','0'}; equ.x0 = {'0','0'}; equ.y0 = {'0','0'}; equ.z0 = {'0','0'}; equ.srcaxis = {{'0';'0';'0'},{'0';'0';'0'}}; equ.coord = {'global','global'}; equ.rOn = {'0','0'}; equ.style = {{{'0'},{'217','217','217'}},{{'0'},{'0','0','255'}}}; equ.ind = [1,1,2]; appl.equ = equ; appl.var = {'t_old_ini','-1'}; xfem.appl{1} = appl; clear appl appl.mode.class = 'MovingMesh'; appl.mode.type = 'cartesian'; appl.dim = {}; appl.sdim = {'X','Y','Z'}; appl.name = 'ale'; appl.shape = {'shlag(2,''lm4'')','shlag(2,''lm5'')','shlag(2,''lm6'')','shlag(2,''x'')','shlag(2,''y'')','shlag(2,''z'')'}; appl.gporder = {10,4}; appl.cporder = 2; appl.sshape = 2; appl.border = 'off'; appl.assignsuffix = '_ale'; clear prop prop.elemdefault='Lag2'; prop.smoothing='winslow'; prop.analysis='static'; prop.allowremesh='off'; prop.definesframe='ale'; prop.origrefframe='ale'; prop.frame='ref'; clear weakconstr weakconstr.value = 'non-ideal'; weakconstr.dim = {'lm4','lm5','lm6'}; prop.weakconstr = weakconstr; appl.prop = prop; clear bnd bnd.name = {'Fixed','','Structural displacement'}; bnd.constrcoord = {'global','global','global'}; bnd.type = {'def','def','def'}; bnd.deform = {{'0';'0';'0'},{'0';'0';'0'},{'u';'v';'w'}}; bnd.veldeform = {{'0';'0';'0'},{'0';'0';'0'},{'0';'0';'0'}}; bnd.defflag = {{'1';'1';'1'},{'0';'0';'0'},{'1';'1';'1'}}; bnd.veldefflag = {{'0';'0';'0'},{'0';'0';'0'},{'0';'0';'0'}}; bnd.weakconstr = {'1','1','1'}; bnd.wcshape = {[1;2;3],[1;2;3],[1;2;3]}; bnd.wcgporder = {{1;1;1},{1;1;1},{1;1;1}}; bnd.wcinit = {{'0';'0';'0'},{'0';'0';'0'},{'0';'0';'0'}}; bnd.style = {{{'0'},{'217','217','217'}},{{'0'},{'0','255','0'}},{{'0'},{'217', ... '217','217'}}}; bnd.ind = [1,1,1,1,1,1,1,1,1,2,1,3,3,2,2,3,3,1]; appl.bnd = bnd; clear equ equ.shape = {[4;5;6],[4;5;6]}; equ.gporder = {{2;2;2},{2;2;2}}; equ.cporder = {{1;1;1},{1;1;1}}; equ.init = {{'xinit_ale';'yinit_ale';'zinit_ale'},{'xinit_ale';'yinit_ale'; ... 'zinit_ale'}}; equ.usage = {1,1}; equ.name = {'Fluid domain','Solid domain'}; equ.type = {'free','phys'}; equ.physexpr = {{'0';'0';'0'},{'u';'v';'w'}}; equ.presexpr = {{'0';'0';'0'},{'0';'0';'0'}}; equ.dinit = {{'0';'0';'0'},{'0';'0';'0'}}; equ.style = {{{'0'},{'217','217','217'}},{{'0'},{'217','217','217'}}}; equ.ind = [1,1,2]; appl.equ = equ; xfem.appl{2} = appl; clear appl appl.mode.class = 'NavierStokes'; appl.mode.type = 'cartesian'; appl.dim = {'u2','v2','w2','p2','logk','logd','nxw','nyw','nzw'}; appl.sdim = {'x','y','z'}; appl.name = 'chns'; appl.module = 'CHEM'; appl.shape = {'shlag(2,''lm7'')','shlag(2,''lm8'')','shlag(2,''lm9'')','shlag(2,''u2'')','shlag(2,''v2'')','shlag(2,''w2'')','shlag(1,''p2'')'}; appl.gporder = {10,4,2}; appl.cporder = {2,1}; appl.sshape = 2; appl.border = 'off'; appl.assignsuffix = '_chns'; clear prop prop.elemdefault='Lagp2p1'; prop.analysis='static'; prop.stensor='full'; prop.cornersmoothing='Off'; prop.nisot='Off'; prop.turbmodel='None'; prop.nonnewtonian='Off'; prop.brinkmandef='Off'; prop.GLSdef='false'; prop.frame='ale'; clear weakconstr weakconstr.value = 'non-ideal'; weakconstr.dim = {'lm7','lm8','lm9','lm10','lm11'}; prop.weakconstr = weakconstr; appl.prop = prop; clear pnt pnt.p0 = {'0'}; pnt.pnton = {'0'}; pnt.style = {{{'0'},{'0','0','0'}}}; pnt.ind = [1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1]; appl.pnt = pnt; clear edg edg.style = {{{'0'},{'0','0','0'},{'solid'}}}; edg.ind = [1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1, ... 1,1,1,1,1,1]; appl.edg = edg; clear bnd bnd.name = {'','','','Structural velocity'}; bnd.type = {'uv','neutral','noslip','uv'}; bnd.u0 = {'.5','0','0','ut'}; bnd.v0 = {'0','0','0','vt'}; bnd.w0 = {'0','0','0','wt'}; bnd.p0 = {'0','0','0','0'}; bnd.k0 = {'0','0','0','0'}; bnd.d0 = {'0','0','0','0'}; bnd.TinType = {'LTIT','LTIT','LTIT','LTIT'}; bnd.LT = {'0.07','0.07','0.07','0.07'}; bnd.IT = {'0.054','0.054','0.054','0.054'}; bnd.uw = {'0','0','0','0'}; bnd.vw = {'0','0','0','0'}; bnd.ww = {'0','0','0','0'}; bnd.dwType = {'dw','dw','dw','dw'}; bnd.dw = {'h','h','h','h'}; bnd.yplus = {'100','100','100','100'}; bnd.relExpr = {{},{},{},{}}; bnd.U0 = {'0','0','0','0'}; bnd.V0 = {'0','0','0','0'}; bnd.p0_inl = {'0','0','0','0'}; bnd.L0 = {'1','1','1','1'}; bnd.flowtype = {'velocity','velocity','velocity','velocity'}; bnd.constr0 = {'0','0','0','0'}; bnd.weakconstr = {'1','1','1','1'}; bnd.wcshape = {[1;2;3],[1;2;3],[1;2;3],[1;2;3]}; bnd.wcgporder = {{1;1;1},{1;1;1},{1;1;1},{1;1;1}}; bnd.wcinit = {{'0';'0';'0';'0';'0'},{'0';'0';'0';'0';'0'},{'0';'0';'0';'0'; ... '0'},{'0';'0';'0';'0';'0'}}; bnd.style = {{{'0'},{'0','255','255'}},{{'0'},{'0','255','0'}},{{'0'},{'0', ... '0','255'}},{{'0'},{'217','217','217'}}}; bnd.ind = [1,2,2,2,2,2,3,3,3,2,3,3,3,3,3,3,3,2]; appl.bnd = bnd; clear equ equ.shape = {[4;5;6;7],[4;5;6;7]}; equ.gporder = {{2;2;2;3},{2;2;2;3}}; equ.cporder = {{1;1;1;2},{1;1;1;2}}; equ.init = {{'0';'0';'0';'0';'0';'0';'0';'0';'0'},{'0';'0';'0';'0';'0';'0'; ... '0';'0';'0'}}; equ.usage = {1,0}; equ.name = {'Fluid domain','Solid domain'}; equ.rho = {'1000','1'}; equ.eta = {'.001','1'}; equ.F_x = {'0','0'}; equ.F_y = {'0','0'}; equ.F_z = {'0','0'}; equ.relExpr = {{},{}}; equ.eta0 = {'1','1'}; equ.m = {'0','0'}; equ.n = {'0','0'}; equ.lambda = {'0','0'}; equ.eta_inf = {'0','0'}; equ.type_visc = {'power','power'}; equ.kappadv = {'0','0'}; equ.k = {'1','1'}; equ.brinkmaneqns = {'0','0'}; equ.idon = {{'0'},{'0'}}; equ.delid = {{'0.5'},{'0.5'}}; equ.sdon = {{'0'},{'0'}}; equ.sdtype = {{'pgc'},{'pgc'}}; equ.delsd = {{'0.25'},{'0.25'}}; equ.cdon = {{'0'},{'0'}}; equ.cdtype = {{'sc'},{'sc'}}; equ.delcd = {{'0.35'},{'0.35'}}; equ.pson = {{'0'},{'0'}}; equ.delps = {{'1'},{'1'}}; equ.idTon = {'1','1'}; equ.delidT = {'0.5','0.5'}; equ.sdTon = {'1','1'}; equ.sdTtype = {'ad','ad'}; equ.delsdT = {'0.5','0.5'}; equ.style = {{{'0'},{'217','217','217'}},{{'0'},{'217','217','217'}}}; equ.ind = [1,1,2]; appl.equ = equ; xfem.appl{3} = appl; xfem.geom = flbinary('fem1','geom','tutorial5.mph'); xfem.mesh = flbinary('fem1.0','mesh','tutorial5.mph'); xfem.sdim = {{'X','Y','Z'},{'x','y','z'}}; xfem.frame = {'ref','ale'}; xfem.shape = {'shlag(''order'',2,''basename'',''x'',''frame'',''ale'')','shlag(''order'',2,''basename'',''y'',''frame'',''ale'')','shlag(''order'',2,''basename'',''z'',''frame'',''ale'')','shlag(''order'',2,''basename'',''u2'',''frame'',''ale'')','shlag(''order'',2,''basename'',''v2'',''frame'',''ale'')','shlag(''order'',2,''basename'',''w2'',''frame'',''ale'')','shlag(''order'',1,''basename'',''p2'',''frame'',''ale'')','shlag(2,''u'')','shlag(2,''v'')','shlag(2,''w'')','shlag(2,''lm4'')','shlag(2,''lm5'')','shlag(2,''lm6'')','shlag(2,''lm7'')','shlag(2,''lm8'')','shlag(2,''lm9'')'}; xfem.gporder = {4,2,10}; xfem.cporder = {2,1}; clear sshape sshape.dvolname = 'dvol'; sshape.frame = 'ref'; sshape.sorder = 2; sshape.type = 'fixed'; xfem.sshape{1} = sshape; clear sshape sshape.dvolname = 'dvol_ale'; sshape.frame = 'ale'; sshape.sorder = 2; sshape.type = 'moving_abs'; sshape.refframe = 'ref'; xfem.sshape{2} = sshape; clear sshape sshape.dvolname = 'dvol_ale'; sshape.frame = 'ale'; sshape.sorder = 2; sshape.type = 'moving_rel'; sshape.sdimdofs = {'u','v','w'}; sshape.refframe = 'ref'; xfem.sshape{3} = sshape; xfem.simplify = 'on'; xfem.border = 1; xfem.form = 'weak'; clear units; units.basesystem = 'SI'; xfem.units = units; clear equ equ.shape = {[1;2;3;4;5;6;7],[8;9;10]}; equ.gporder = {{1;1;1;1;1;1;1;1;1;1;2},{1;1;1;1;1;1;1;1;1;1;2}}; equ.cporder = {{1;1;1;1;1;1;1;1;1;1;2},{1;1;1;1;1;1;1;1;1;1;2}}; equ.init = {{'0';'0';'0';'0';'xinit_ale';'yinit_ale';'zinit_ale';'0';'0'; ... '0';'0'},{'0';'0';'0';'0';'0';'0';'0';'0';'0';'0';'0'}}; equ.dinit = {{'0';'0';'0';'0';'0';'0';'0';'0';'0';'0';'0'},{'0';'0';'0';'0'; ... '0';'0';'0';'0';'0';'0';'0'}}; equ.weak = {{'0';'0';'0';'0';'-Xx*test(xx)-Xy*test(xy)-Xz*test(xz)';'-Yx*test(yx)-Yy*test(yy)-Yz*test(yz)'; ... '-Zx*test(zx)-Zy*test(zy)-Zz*test(zz)';'(-2*eta_chns*u2x+p2)*test(u2x)-eta_chns*(u2y+v2x)*test(u2y)-eta_chns*(u2z+w2x)*test(u2z)-rho_chns*(u2*u2x+v2*u2y+w2*u2z)*test(u2)'; ... '-eta_chns*(v2x+u2y)*test(v2x)+(-2*eta_chns*v2y+p2)*test(v2y)-eta_chns*(v2z+w2y)*test(v2z)-rho_chns*(u2*v2x+v2*v2y+w2*v2z)*test(v2)'; ... '-eta_chns*(w2x+u2z)*test(w2x)-eta_chns*(w2y+v2z)*test(w2y)+(-2*eta_chns*w2z+p2)*test(w2z)-rho_chns*(u2*w2x+v2*w2y+w2*w2z)*test(w2)'; ... '-divU_chns*test(p2)'},{'-eX_smsld_test*SX_smsld-eY_smsld_test*SY_smsld-eZ_smsld_test*SZ_smsld-2*eXY_smsld_test*SXY_smsld-2*eYZ_smsld_test*SYZ_smsld-2*eXZ_smsld_test*SXZ_smsld'; ... '0';'0';'0';'0';'0';'0';'0';'0';'0';'0'}}; equ.dweak = {{'0';'0';'0';'0';'0';'0';'0';'rho_chns*u2t*test(u2)';'rho_chns*v2t*test(v2)'; ... 'rho_chns*w2t*test(w2)';'0'},{'0';'0';'0';'0';'0';'0';'0';'0';'0';'0';'0'}}; equ.constr = {{'0';'0';'0';'0';'0';'0';'0';'0';'0';'0';'0'},{'0';'0';'0';'0'; ... '0';'0';'0';'0';'0';'0';'0'}}; equ.sshape = {[1;2],[1;3]}; equ.sshapedim = {{1;1;1;1;2;2;2;2;2;2;2},{1;1;1;1;3;3;3;3;3;3;3}}; equ.ind = [1,1,2]; equ.dim = {'u','v','w','p','x','y','z','u2','v2','w2','p2'}; equ.var = {'FXg_smsld',{'','0'}, ... 'FYg_smsld',{'','0'}, ... 'FZg_smsld',{'','0'}, ... 'disp_smsld',{'','sqrt(real(u)^2+real(v)^2+real(w)^2)'}, ... 'sX_smsld',{'','(F11_smsld*(SX_smsld*F11_smsld+SXY_smsld*F12_smsld+SXZ_smsld*F13_smsld)+F12_smsld*(SXY_smsld*F11_smsld+SY_smsld*F12_smsld+SYZ_smsld*F13_smsld)+F13_smsld*(SXZ_smsld*F11_smsld+SYZ_smsld*F12_smsld+SZ_smsld*F13_smsld))/J_smsld'}, ... 'sY_smsld',{'','(F21_smsld*(SX_smsld*F21_smsld+SXY_smsld*F22_smsld+SXZ_smsld*F23_smsld)+F22_smsld*(SXY_smsld*F21_smsld+SY_smsld*F22_smsld+SYZ_smsld*F23_smsld)+F23_smsld*(SXZ_smsld*F21_smsld+SYZ_smsld*F22_smsld+SZ_smsld*F23_smsld))/J_smsld'}, ... 'sZ_smsld',{'','(F31_smsld*(SX_smsld*F31_smsld+SXY_smsld*F32_smsld+SXZ_smsld*F33_smsld)+F32_smsld*(SXY_smsld*F31_smsld+SY_smsld*F32_smsld+SYZ_smsld*F33_smsld)+F33_smsld*(SXZ_smsld*F31_smsld+SYZ_smsld*F32_smsld+SZ_smsld*F33_smsld))/J_smsld'}, ... 'sXY_smsld',{'','(F11_smsld*(SX_smsld*F21_smsld+SXY_smsld*F22_smsld+SXZ_smsld*F23_smsld)+F12_smsld*(SXY_smsld*F21_smsld+SY_smsld*F22_smsld+SYZ_smsld*F23_smsld)+F13_smsld*(SXZ_smsld*F21_smsld+SYZ_smsld*F22_smsld+SZ_smsld*F23_smsld))/J_smsld'}, ... 'sYZ_smsld',{'','(F21_smsld*(SX_smsld*F31_smsld+SXY_smsld*F32_smsld+SXZ_smsld*F33_smsld)+F22_smsld*(SXY_smsld*F31_smsld+SY_smsld*F32_smsld+SYZ_smsld*F33_smsld)+F23_smsld*(SXZ_smsld*F31_smsld+SYZ_smsld*F32_smsld+SZ_smsld*F33_smsld))/J_smsld'}, ... 'sXZ_smsld',{'','(F11_smsld*(SX_smsld*F31_smsld+SXY_smsld*F32_smsld+SXZ_smsld*F33_smsld)+F12_smsld*(SXY_smsld*F31_smsld+SY_smsld*F32_smsld+SYZ_smsld*F33_smsld)+F13_smsld*(SXZ_smsld*F31_smsld+SYZ_smsld*F32_smsld+SZ_smsld*F33_smsld))/J_smsld'}, ... 'eX_smsld',{'','uX+0.5*(uX^2+vX^2+wX^2)'}, ... 'eY_smsld',{'','vY+0.5*(uY^2+vY^2+wY^2)'}, ... 'eZ_smsld',{'','wZ+0.5*(uZ^2+vZ^2+wZ^2)'}, ... 'eXY_smsld',{'','0.5*(uY+vX+uX*uY+vX*vY+wX*wY)'}, ... 'eYZ_smsld',{'','0.5*(vZ+wY+uY*uZ+vY*vZ+wY*wZ)'}, ... 'eXZ_smsld',{'','0.5*(uZ+wX+uX*uZ+vX*vZ+wX*wZ)'}, ... 'SX_smsld',{'','E_smsld*(1-nu_smsld)*eX_smsld/((1+nu_smsld)*(1-2*nu_smsld))+E_smsld*nu_smsld*eY_smsld/((1+nu_smsld)*(1-2*nu_smsld))+E_smsld*nu_smsld*eZ_smsld/((1+nu_smsld)*(1-2*nu_smsld))'}, ... 'SY_smsld',{'','E_smsld*nu_smsld*eX_smsld/((1+nu_smsld)*(1-2*nu_smsld))+E_smsld*(1-nu_smsld)*eY_smsld/((1+nu_smsld)*(1-2*nu_smsld))+E_smsld*nu_smsld*eZ_smsld/((1+nu_smsld)*(1-2*nu_smsld))'}, ... 'SZ_smsld',{'','E_smsld*nu_smsld*eX_smsld/((1+nu_smsld)*(1-2*nu_smsld))+E_smsld*nu_smsld*eY_smsld/((1+nu_smsld)*(1-2*nu_smsld))+E_smsld*(1-nu_smsld)*eZ_smsld/((1+nu_smsld)*(1-2*nu_smsld))'}, ... 'SXY_smsld',{'','E_smsld*eXY_smsld/(1+nu_smsld)'}, ... 'SYZ_smsld',{'','E_smsld*eYZ_smsld/(1+nu_smsld)'}, ... 'SXZ_smsld',{'','E_smsld*eXZ_smsld/(1+nu_smsld)'}, ... 'K_smsld',{'','E_smsld/(3*(1-2*nu_smsld))'}, ... 'G_smsld',{'','0.5*E_smsld/(1+nu_smsld)'}, ... 'mises_smsld',{'','sqrt(sX_smsld^2+sY_smsld^2+sZ_smsld^2-sX_smsld*sY_smsld-sY_smsld*sZ_smsld-sX_smsld*sZ_smsld+3*sXY_smsld^2+3*sYZ_smsld^2+3*sXZ_smsld^2)'}, ... 'Ws_smsld',{'','0.5*(sX_smsld*eX_smsld+sY_smsld*eY_smsld+sZ_smsld*eZ_smsld+2*sXY_smsld*eXY_smsld+2*sYZ_smsld*eYZ_smsld+2*sXZ_smsld*eXZ_smsld)'}, ... 'evol_smsld',{'','-1+Jel_smsld'}, ... 'F11_smsld',{'','1+uX'}, ... 'F12_smsld',{'','uY'}, ... 'F13_smsld',{'','uZ'}, ... 'F21_smsld',{'','vX'}, ... 'F22_smsld',{'','1+vY'}, ... 'F23_smsld',{'','vZ'}, ... 'F31_smsld',{'','wX'}, ... 'F32_smsld',{'','wY'}, ... 'F33_smsld',{'','1+wZ'}, ... 'detF_smsld',{'','F11_smsld*F22_smsld*F33_smsld+F12_smsld*F23_smsld*F31_smsld+F13_smsld*F21_smsld*F32_smsld-F11_smsld*F23_smsld*F32_smsld-F12_smsld*F21_smsld*F33_smsld-F13_smsld*F22_smsld*F31_smsld'}, ... 'J_smsld',{'','detF_smsld'}, ... 'Jel_smsld',{'','J_smsld'}, ... 'invF11_smsld',{'','(F22_smsld*F33_smsld-F23_smsld*F32_smsld)/detF_smsld'}, ... 'invF12_smsld',{'','(F13_smsld*F32_smsld-F12_smsld*F33_smsld)/detF_smsld'}, ... 'invF13_smsld',{'','(F12_smsld*F23_smsld-F13_smsld*F22_smsld)/detF_smsld'}, ... 'invF21_smsld',{'','(F31_smsld*F23_smsld-F21_smsld*F33_smsld)/detF_smsld'}, ... 'invF22_smsld',{'','(F11_smsld*F33_smsld-F31_smsld*F13_smsld)/detF_smsld'}, ... 'invF23_smsld',{'','(F21_smsld*F13_smsld-F11_smsld*F23_smsld)/detF_smsld'}, ... 'invF31_smsld',{'','(F21_smsld*F32_smsld-F31_smsld*F22_smsld)/detF_smsld'}, ... 'invF32_smsld',{'','(F31_smsld*F12_smsld-F11_smsld*F32_smsld)/detF_smsld'}, ... 'invF33_smsld',{'','(F11_smsld*F22_smsld-F21_smsld*F12_smsld)/detF_smsld'}, ... 'tresca_smsld',{'','max(max(abs(s1_smsld-s2_smsld),abs(s2_smsld-s3_smsld)),abs(s1_smsld-s3_smsld))'}, ... 'xinit_ale',{'X','X'}, ... 'yinit_ale',{'Y','Y'}, ... 'zinit_ale',{'Z','Z'}, ... 'dx_ale',{'x-X','x-X'}, ... 'dy_ale',{'y-Y','y-Y'}, ... 'dz_ale',{'z-Z','z-Z'}, ... 'U_chns',{'sqrt(u2^2+v2^2+w2^2)',''}, ... 'Vx_chns',{'w2y-v2z',''}, ... 'Vy_chns',{'u2z-w2x',''}, ... 'Vz_chns',{'v2x-u2y',''}, ... 'V_chns',{'sqrt(Vx_chns^2+Vy_chns^2+Vz_chns^2)',''}, ... 'divU_chns',{'u2x+v2y+w2z',''}, ... 'cellRe_chns',{'rho_chns*U_chns*h_ale/eta_chns',''}, ... 'res_u2_chns',{'rho_chns*(u2*u2x+v2*u2y+w2*u2z)+p2x-F_x_chns-eta_chns*(2*u2xx+u2yy+v2xy+u2zz+w2xz)', ... ''}, ... 'res_tst_u2_chns',{'nojac(rho_chns)*(nojac(u2)*u2x+nojac(v2)*u2y+nojac(w2)*u2z)+p2x-nojac(eta_chns)*(2*u2xx+u2yy+v2xy+u2zz+w2xz)', ... ''}, ... 'res_sc_u2_chns',{'rho_chns*(u2*u2x+v2*u2y+w2*u2z)+p2x-F_x_chns', ... ''}, ... 'res_v2_chns',{'rho_chns*(u2*v2x+v2*v2y+w2*v2z)+p2y-F_y_chns-eta_chns*(v2xx+u2yx+2*v2yy+v2zz+w2yz)', ... ''}, ... 'res_tst_v2_chns',{'nojac(rho_chns)*(nojac(u2)*v2x+nojac(v2)*v2y+nojac(w2)*v2z)+p2y-nojac(eta_chns)*(v2xx+u2yx+2*v2yy+v2zz+w2yz)', ... ''}, ... 'res_sc_v2_chns',{'rho_chns*(u2*v2x+v2*v2y+w2*v2z)+p2y-F_y_chns', ... ''}, ... 'res_w2_chns',{'rho_chns*(u2*w2x+v2*w2y+w2*w2z)+p2z-F_z_chns-eta_chns*(w2xx+u2zx+w2yy+v2zy+2*w2zz)', ... ''}, ... 'res_tst_w2_chns',{'nojac(rho_chns)*(nojac(u2)*w2x+nojac(v2)*w2y+nojac(w2)*w2z)+p2z-nojac(eta_chns)*(w2xx+u2zx+w2yy+v2zy+2*w2zz)', ... ''}, ... 'res_sc_w2_chns',{'rho_chns*(u2*w2x+v2*w2y+w2*w2z)+p2z-F_z_chns', ... ''}, ... 'beta_x_chns',{'rho_chns*u2',''}, ... 'beta_y_chns',{'rho_chns*v2',''}, ... 'beta_z_chns',{'rho_chns*w2',''}, ... 'Dm_chns',{'eta_chns',''}, ... 'da_chns',{'rho_chns',''},'E_smsld',{'','31100'}, ... 'rho_smsld',{'','7850'}, ... 'alphadM_smsld',{'','1'}, ... 'betadK_smsld',{'','0.001'}, ... 'eta_smsld',{'','0'}, ... 'mixedform_smsld',{'','0'}, ... 'mu_smsld',{'','800000'}, ... 'C10_smsld',{'','200000'}, ... 'C01_smsld',{'','200000'}, ... 'kappa_smsld',{'','10000000000'}, ... 'ETkin_smsld',{'','20000000000'}, ... 'ETiso_smsld',{'','20000000000'}, ... 'Sys_smsld',{'','200000000'}, ... 'Syfunc_smsld',{'','mises_smsld'}, ... 'Syfunc_kin_smsld',{'','misesKin_smsld'}, ... 'Shard_smsld',{'','22222222222.22222*epe_smsld'}, ... 'alpha_smsld',{'','1.2e-005'}, ... 'nu_smsld',{'','0.49'}, ... 'EX_smsld',{'','200000000000'}, ... 'nuXY_smsld',{'','0.33'}, ... 'GXY_smsld',{'','75200000000'}, ... 'alphaX_smsld',{'','1.2e-005'}, ... 'sXi_smsld',{'','0'}, ... 'sXYi_smsld',{'','0'}, ... 'eXi_smsld',{'','0'}, ... 'eXYi_smsld',{'','0'}, ... 'EY_smsld',{'','200000000000'}, ... 'nuYZ_smsld',{'','0.33'}, ... 'GYZ_smsld',{'','75200000000'}, ... 'alphaY_smsld',{'','1.2e-005'}, ... 'sYi_smsld',{'','0'}, ... 'sYZi_smsld',{'','0'}, ... 'eYi_smsld',{'','0'}, ... 'eYZi_smsld',{'','0'}, ... 'EZ_smsld',{'','200000000000'}, ... 'nuXZ_smsld',{'','0.33'}, ... 'GXZ_smsld',{'','75200000000'}, ... 'alphaZ_smsld',{'','1.2e-005'}, ... 'sZi_smsld',{'','0'}, ... 'sXZi_smsld',{'','0'}, ... 'eZi_smsld',{'','0'}, ... 'eXZi_smsld',{'','0'}, ... 'FX_smsld',{'','0'}, ... 'FXAmp_smsld',{'','1'}, ... 'FXPh_smsld',{'','0'}, ... 'FY_smsld',{'','0'}, ... 'FYAmp_smsld',{'','1'}, ... 'FYPh_smsld',{'','0'}, ... 'FZ_smsld',{'','0'}, ... 'FZAmp_smsld',{'','1'}, ... 'FZPh_smsld',{'','0'}, ... 'Temp_smsld',{'','0'}, ... 'Tempref_smsld',{'','0'}, ... 'velx_smsld',{'','0'}, ... 'accx_smsld',{'','0'}, ... 'vely_smsld',{'','0'}, ... 'accy_smsld',{'','0'}, ... 'velz_smsld',{'','0'}, ... 'accz_smsld',{'','0'}, ... 'Rx_smsld',{'','0'}, ... 'Ry_smsld',{'','0'}, ... 'Rz_smsld',{'','0'}, ... 'dx_smsld',{'','1'}, ... 'Lx_smsld',{'','dx_smsld'}, ... 'dy_smsld',{'','1'}, ... 'Ly_smsld',{'','dy_smsld'}, ... 'dz_smsld',{'','1'}, ... 'Lz_smsld',{'','dz_smsld'}, ... 'dr_smsld',{'','1'}, ... 'Lr_smsld',{'','dr_smsld'}, ... 'R0_smsld',{'','1'}, ... 'x0_smsld',{'','0'}, ... 'y0_smsld',{'','0'}, ... 'z0_smsld',{'','0'}, ... 'H11_smsld',{'','0'}, ... 'H12_smsld',{'','0'}, ... 'H13_smsld',{'','0'}, ... 'H21_smsld',{'','0'}, ... 'H22_smsld',{'','0'}, ... 'H23_smsld',{'','0'}, ... 'H31_smsld',{'','0'}, ... 'H32_smsld',{'','0'}, ... 'H33_smsld',{'','0'}, ... 'rho_chns',{'1000',''}, ... 'eta_chns',{'0.001',''}, ... 'F_x_chns',{'0',''}, ... 'F_y_chns',{'0',''}, ... 'F_z_chns',{'0',''}, ... 'eta0_chns',{'1',''}, ... 'm_chns',{'0',''}, ... 'n_chns',{'0',''}, ... 'lambda_chns',{'0',''}, ... 'eta_inf_chns',{'0',''}, ... 'type_visc_chns',{'power',''}, ... 'kappadv_chns',{'0',''}, ... 'k_chns',{'1',''}, ... 'brinkmaneqns_chns',{'0',''}}; equ.expr = {}; equ.bnd.weak = {{'0';'0';'0';'0';'0';'0';'0';'0';'0';'0';'0'}}; equ.bnd.gporder = {{1;1;1;1;1;1;1;1;1;1;2}}; equ.bnd.ind = [1,1,1]; equ.bnd.expr = {}; equ.lock = [0,0,0]; equ.mlock = {[0,0,0]}; xfem.equ = equ; clear bnd bnd.weak = {{'0';'0';'0';'0';'0';'0';'0';'-lm4_test*(-x+X)+x_test*lm4';'-lm5_test*(-y+Y)+y_test*lm5'; ... '-lm6_test*(-z+Z)+z_test*lm6';'0';'0';'0';'0';'-lm7_test*(-u2+u0_chns)+u2_test*lm7'; ... 'lm8_test*v2+v2_test*lm8';'lm9_test*w2+w2_test*lm9'},{'0';'0';'0';'0';'0'; ... '0';'0';'-lm4_test*(-x+X)+x_test*lm4';'-lm5_test*(-y+Y)+y_test*lm5';'-lm6_test*(-z+Z)+z_test*lm6'; ... '0';'0';'0';'0';'0';'0';'0'},{'0';'0';'0';'0';'0';'0';'0';'-lm4_test*(-x+X)+x_test*lm4'; ... '-lm5_test*(-y+Y)+y_test*lm5';'-lm6_test*(-z+Z)+z_test*lm6';'0';'0';'0'; ... '0';'lm7_test*u2+u2_test*lm7';'lm8_test*v2+v2_test*lm8';'lm9_test*w2+w2_test*lm9'}, ... {'0';'0';'0';'0';'0';'0';'0';'0';'0';'0';'0';'0';'0';'0';'0';'0';'0'},{'FXg_smsld*u_test'; ... 'FYg_smsld*v_test';'FZg_smsld*w_test';'0';'0';'0';'0';'-lm4_test*(u-x+X)+x_test*lm4'; ... '-lm5_test*(v-y+Y)+y_test*lm5';'-lm6_test*(w-z+Z)+z_test*lm6';'0';'0';'0'; ... '0';'lm7_test*u2+u2_test*lm7';'lm8_test*v2+v2_test*lm8';'lm9_test*w2+w2_test*lm9'}, ... {'0';'0';'0';'0';'0';'0';'0';'0';'0';'0';'0';'0';'0';'0';'0';'0';'0'},{'FXg_smsld*u_test'; ... 'FYg_smsld*v_test';'FZg_smsld*w_test';'0';'0';'0';'0';'-lm4_test*(u-x+X)+x_test*lm4'; ... '-lm5_test*(v-y+Y)+y_test*lm5';'-lm6_test*(w-z+Z)+z_test*lm6';'0';'0';'0'; ... '0';'lm7_test*u2+u2_test*lm7';'lm8_test*v2+v2_test*lm8';'lm9_test*w2+w2_test*lm9'}}; bnd.dweak = {{'0';'0';'0';'0';'0';'0';'0';'0';'0';'0';'0';'0';'0';'0';'0'; ... '0';'0'},{'0';'0';'0';'0';'0';'0';'0';'0';'0';'0';'0';'0';'0';'0';'0';'0'; ... '0'},{'0';'0';'0';'0';'0';'0';'0';'0';'0';'0';'0';'0';'0';'0';'0';'0';'0'}, ... {'0';'0';'0';'0';'0';'0';'0';'0';'0';'0';'0';'0';'0';'0';'0';'0';'0'},{'0'; ... '0';'0';'0';'0';'0';'0';'0';'0';'0';'0';'0';'0';'0';'0';'0';'0'},{'0';'0'; ... '0';'0';'0';'0';'0';'0';'0';'0';'0';'0';'0';'0';'0';'0';'0'},{'0';'0';'0'; ... '0';'0';'0';'0';'0';'0';'0';'0';'0';'0';'0';'0';'0';'0'}}; bnd.constr = {{'0';'0';'0';'0';'0';'0';'0';'0';'0';'0';'0';'0';'0';'0';'0'; ... '0';'0'},{'0';'0';'0';'0';'0';'0';'0';'0';'0';'0';'0';'0';'0';'0';'0';'0'; ... '0'},{'0';'0';'0';'0';'0';'0';'0';'0';'0';'0';'0';'0';'0';'0';'0';'0';'0'}, ... {'0';'0';'0';'0';'0';'0';'0';'0';'0';'0';'0';'0';'0';'0';'0';'0';'0'},{'0'; ... '0';'0';'0';'0';'0';'0';'0';'0';'0';'0';'0';'0';'0';'0';'0';'0'},{'-u'; ... '-v';'-w';'0';'0';'0';'0';'0';'0';'0';'0';'0';'0';'0';'0';'0';'0'},{'0'; ... '0';'0';'0';'0';'0';'0';'0';'0';'0';'0';'0';'0';'0';'0';'0';'0'}}; bnd.shape = {[1;2;3;11;12;13;4;5;6;7;14;15;16],[1;2;3;11;12;13;4;5;6;7],[1; ... 2;3;11;12;13;4;5;6;7;14;15;16],[1;2;3;4;5;6;7],[8;9;10;1;2;3;11;12;13;4; ... 5;6;7;14;15;16],[8;9;10],[8;9;10;1;2;3;11;12;13;4;5;6;7;14;15;16]}; bnd.sshape = {[1;2],[1;2],[1;2],[1;2],[1;2],[1;3],[1;2]}; bnd.sshapedim = {{1;1;1;1;2;2;2;2;2;2;2;2;2;2;2;2;2},{1;1;1;1;2;2;2;2;2;2;2;2; ... 2;2;2;2;2},{1;1;1;1;2;2;2;2;2;2;2;2;2;2;2;2;2},{1;1;1;1;2;2;2;2;2;2;2;2; ... 2;2;2;2;2},{1;1;1;1;2;2;2;2;2;2;2;2;2;2;2;2;2},{1;1;1;1;3;3;3;3;3;3;3;3; ... 3;3;3;3;3},{1;1;1;1;2;2;2;2;2;2;2;2;2;2;2;2;2}}; bnd.gporder = {{1;1;1;1;1;1;1;3;3;3;1;1;1;2;3;3;3},{1;1;1;1;1;1;1;3;3;3;1;1; ... 1;2;3;3;3},{1;1;1;1;1;1;1;3;3;3;1;1;1;2;3;3;3},{1;1;1;1;1;1;1;3;3;3;1;1; ... 1;2;3;3;3},{1;1;1;1;1;1;1;3;3;3;1;1;1;2;3;3;3},{1;1;1;1;1;1;1;3;3;3;1;1; ... 1;2;3;3;3},{1;1;1;1;1;1;1;3;3;3;1;1;1;2;3;3;3}}; bnd.cporder = {{1;1;1;1;1;1;1;1;1;1;1;1;1;2;1;1;1},{1;1;1;1;1;1;1;1;1;1;1;1; ... 1;2;1;1;1},{1;1;1;1;1;1;1;1;1;1;1;1;1;2;1;1;1},{1;1;1;1;1;1;1;1;1;1;1;1; ... 1;2;1;1;1},{1;1;1;1;1;1;1;1;1;1;1;1;1;2;1;1;1},{1;1;1;1;1;1;1;1;1;1;1;1; ... 1;2;1;1;1},{1;1;1;1;1;1;1;1;1;1;1;1;1;2;1;1;1}}; bnd.init = {{'';'';'';'';'';'';'';'0';'0';'0';'';'';'';'';'0';'0';'0'},{''; ... '';'';'';'';'';'';'0';'0';'0';'';'';'';'';'0';'0';'0'},{'';'';'';'';''; ... '';'';'0';'0';'0';'';'';'';'';'0';'0';'0'},{'';'';'';'';'';'';'';'0';'0'; ... '0';'';'';'';'';'0';'0';'0'},{'';'';'';'';'';'';'';'0';'0';'0';'';'';''; ... '';'0';'0';'0'},{'';'';'';'';'';'';'';'0';'0';'0';'';'';'';'';'0';'0';'0'}, ... {'';'';'';'';'';'';'';'0';'0';'0';'';'';'';'';'0';'0';'0'}}; bnd.dinit = {{'';'';'';'';'';'';'';'';'';'';'';'';'';'';'';'';''},{'';'';''; ... '';'';'';'';'';'';'';'';'';'';'';'';'';''},{'';'';'';'';'';'';'';'';''; ... '';'';'';'';'';'';'';''},{'';'';'';'';'';'';'';'';'';'';'';'';'';'';''; ... '';''},{'';'';'';'';'';'';'';'';'';'';'';'';'';'';'';'';''},{'';'';'';''; ... '';'';'';'';'';'';'';'';'';'';'';'';''},{'';'';'';'';'';'';'';'';'';''; ... '';'';'';'';'';'';''}}; bnd.ind = [1,2,2,2,2,2,3,3,3,4,3,5,5,6,6,7,7,2]; bnd.dim = {'u','v','w','p','x','y','z','lm4','lm5','lm6','u2','v2','w2', ... 'p2','lm7','lm8','lm9'}; bnd.var = {'FXg_smsld',{'','','','','FX_smsld','0','FX_smsld'}, ... 'FYg_smsld',{'','','','','FY_smsld','0','FY_smsld'}, ... 'FZg_smsld',{'','','','','FZ_smsld','0','FZ_smsld'}, ... 'disp_smsld',{'','','','','sqrt(real(u)^2+real(v)^2+real(w)^2)','sqrt(real(u)^2+real(v)^2+real(w)^2)', ... 'sqrt(real(u)^2+real(v)^2+real(w)^2)'}, ... 'TaX_smsld',{'','','','','(F11_smsld*SX_smsld+F12_smsld*SXY_smsld+F13_smsld*SXZ_smsld)*nX_smsld+(F11_smsld*SXY_smsld+F12_smsld*SY_smsld+F13_smsld*SYZ_smsld)*nY_smsld+(F11_smsld*SXZ_smsld+F12_smsld*SYZ_smsld+F13_smsld*SZ_smsld)*nZ_smsld', ... '(F11_smsld*SX_smsld+F12_smsld*SXY_smsld+F13_smsld*SXZ_smsld)*nX_smsld+(F11_smsld*SXY_smsld+F12_smsld*SY_smsld+F13_smsld*SYZ_smsld)*nY_smsld+(F11_smsld*SXZ_smsld+F12_smsld*SYZ_smsld+F13_smsld*SZ_smsld)*nZ_smsld', ... '(F11_smsld*SX_smsld+F12_smsld*SXY_smsld+F13_smsld*SXZ_smsld)*nX_smsld+(F11_smsld*SXY_smsld+F12_smsld*SY_smsld+F13_smsld*SYZ_smsld)*nY_smsld+(F11_smsld*SXZ_smsld+F12_smsld*SYZ_smsld+F13_smsld*SZ_smsld)*nZ_smsld'}, ... 'TaY_smsld',{'','','','','(F21_smsld*SX_smsld+F22_smsld*SXY_smsld+F23_smsld*SXZ_smsld)*nX_smsld+(F21_smsld*SXY_smsld+F22_smsld*SY_smsld+F23_smsld*SYZ_smsld)*nY_smsld+(F21_smsld*SXZ_smsld+F22_smsld*SYZ_smsld+F23_smsld*SZ_smsld)*nZ_smsld', ... '(F21_smsld*SX_smsld+F22_smsld*SXY_smsld+F23_smsld*SXZ_smsld)*nX_smsld+(F21_smsld*SXY_smsld+F22_smsld*SY_smsld+F23_smsld*SYZ_smsld)*nY_smsld+(F21_smsld*SXZ_smsld+F22_smsld*SYZ_smsld+F23_smsld*SZ_smsld)*nZ_smsld', ... '(F21_smsld*SX_smsld+F22_smsld*SXY_smsld+F23_smsld*SXZ_smsld)*nX_smsld+(F21_smsld*SXY_smsld+F22_smsld*SY_smsld+F23_smsld*SYZ_smsld)*nY_smsld+(F21_smsld*SXZ_smsld+F22_smsld*SYZ_smsld+F23_smsld*SZ_smsld)*nZ_smsld'}, ... 'TaZ_smsld',{'','','','','(F31_smsld*SX_smsld+F32_smsld*SXY_smsld+F33_smsld*SXZ_smsld)*nX_smsld+(F31_smsld*SXY_smsld+F32_smsld*SY_smsld+F33_smsld*SYZ_smsld)*nY_smsld+(F31_smsld*SXZ_smsld+F32_smsld*SYZ_smsld+F33_smsld*SZ_smsld)*nZ_smsld', ... '(F31_smsld*SX_smsld+F32_smsld*SXY_smsld+F33_smsld*SXZ_smsld)*nX_smsld+(F31_smsld*SXY_smsld+F32_smsld*SY_smsld+F33_smsld*SYZ_smsld)*nY_smsld+(F31_smsld*SXZ_smsld+F32_smsld*SYZ_smsld+F33_smsld*SZ_smsld)*nZ_smsld', ... '(F31_smsld*SX_smsld+F32_smsld*SXY_smsld+F33_smsld*SXZ_smsld)*nX_smsld+(F31_smsld*SXY_smsld+F32_smsld*SY_smsld+F33_smsld*SYZ_smsld)*nY_smsld+(F31_smsld*SXZ_smsld+F32_smsld*SYZ_smsld+F33_smsld*SZ_smsld)*nZ_smsld'}, ... 'K_x_chns',{'2*nx_chns*eta_chns*u2x+ny_chns*eta_chns*(u2y+v2x)+nz_chns*eta_chns*(u2z+w2x)', ... '2*nx_chns*eta_chns*u2x+ny_chns*eta_chns*(u2y+v2x)+nz_chns*eta_chns*(u2z+w2x)', ... '2*nx_chns*eta_chns*u2x+ny_chns*eta_chns*(u2y+v2x)+nz_chns*eta_chns*(u2z+w2x)', ... '2*nx_chns*eta_chns*u2x+ny_chns*eta_chns*(u2y+v2x)+nz_chns*eta_chns*(u2z+w2x)', ... '2*nx_chns*eta_chns*u2x+ny_chns*eta_chns*(u2y+v2x)+nz_chns*eta_chns*(u2z+w2x)', ... '','2*nx_chns*eta_chns*u2x+ny_chns*eta_chns*(u2y+v2x)+nz_chns*eta_chns*(u2z+w2x)'}, ... 'T_x_chns',{'-nx_chns*p2+2*nx_chns*eta_chns*u2x+ny_chns*eta_chns*(u2y+v2x)+nz_chns*eta_chns*(u2z+w2x)', ... '-nx_chns*p2+2*nx_chns*eta_chns*u2x+ny_chns*eta_chns*(u2y+v2x)+nz_chns*eta_chns*(u2z+w2x)', ... '-nx_chns*p2+2*nx_chns*eta_chns*u2x+ny_chns*eta_chns*(u2y+v2x)+nz_chns*eta_chns*(u2z+w2x)', ... '-nx_chns*p2+2*nx_chns*eta_chns*u2x+ny_chns*eta_chns*(u2y+v2x)+nz_chns*eta_chns*(u2z+w2x)', ... '-nx_chns*p2+2*nx_chns*eta_chns*u2x+ny_chns*eta_chns*(u2y+v2x)+nz_chns*eta_chns*(u2z+w2x)', ... '','-nx_chns*p2+2*nx_chns*eta_chns*u2x+ny_chns*eta_chns*(u2y+v2x)+nz_chns*eta_chns*(u2z+w2x)'}, ... 'K_y_chns',{'nx_chns*eta_chns*(v2x+u2y)+2*ny_chns*eta_chns*v2y+nz_chns*eta_chns*(v2z+w2y)', ... 'nx_chns*eta_chns*(v2x+u2y)+2*ny_chns*eta_chns*v2y+nz_chns*eta_chns*(v2z+w2y)', ... 'nx_chns*eta_chns*(v2x+u2y)+2*ny_chns*eta_chns*v2y+nz_chns*eta_chns*(v2z+w2y)', ... 'nx_chns*eta_chns*(v2x+u2y)+2*ny_chns*eta_chns*v2y+nz_chns*eta_chns*(v2z+w2y)', ... 'nx_chns*eta_chns*(v2x+u2y)+2*ny_chns*eta_chns*v2y+nz_chns*eta_chns*(v2z+w2y)', ... '','nx_chns*eta_chns*(v2x+u2y)+2*ny_chns*eta_chns*v2y+nz_chns*eta_chns*(v2z+w2y)'}, ... 'T_y_chns',{'-ny_chns*p2+nx_chns*eta_chns*(v2x+u2y)+2*ny_chns*eta_chns*v2y+nz_chns*eta_chns*(v2z+w2y)', ... '-ny_chns*p2+nx_chns*eta_chns*(v2x+u2y)+2*ny_chns*eta_chns*v2y+nz_chns*eta_chns*(v2z+w2y)', ... '-ny_chns*p2+nx_chns*eta_chns*(v2x+u2y)+2*ny_chns*eta_chns*v2y+nz_chns*eta_chns*(v2z+w2y)', ... '-ny_chns*p2+nx_chns*eta_chns*(v2x+u2y)+2*ny_chns*eta_chns*v2y+nz_chns*eta_chns*(v2z+w2y)', ... '-ny_chns*p2+nx_chns*eta_chns*(v2x+u2y)+2*ny_chns*eta_chns*v2y+nz_chns*eta_chns*(v2z+w2y)', ... '','-ny_chns*p2+nx_chns*eta_chns*(v2x+u2y)+2*ny_chns*eta_chns*v2y+nz_chns*eta_chns*(v2z+w2y)'}, ... 'K_z_chns',{'nx_chns*eta_chns*(w2x+u2z)+ny_chns*eta_chns*(w2y+v2z)+2*nz_chns*eta_chns*w2z', ... 'nx_chns*eta_chns*(w2x+u2z)+ny_chns*eta_chns*(w2y+v2z)+2*nz_chns*eta_chns*w2z', ... 'nx_chns*eta_chns*(w2x+u2z)+ny_chns*eta_chns*(w2y+v2z)+2*nz_chns*eta_chns*w2z', ... 'nx_chns*eta_chns*(w2x+u2z)+ny_chns*eta_chns*(w2y+v2z)+2*nz_chns*eta_chns*w2z', ... 'nx_chns*eta_chns*(w2x+u2z)+ny_chns*eta_chns*(w2y+v2z)+2*nz_chns*eta_chns*w2z', ... '','nx_chns*eta_chns*(w2x+u2z)+ny_chns*eta_chns*(w2y+v2z)+2*nz_chns*eta_chns*w2z'}, ... 'T_z_chns',{'-nz_chns*p2+nx_chns*eta_chns*(w2x+u2z)+ny_chns*eta_chns*(w2y+v2z)+2*nz_chns*eta_chns*w2z', ... '-nz_chns*p2+nx_chns*eta_chns*(w2x+u2z)+ny_chns*eta_chns*(w2y+v2z)+2*nz_chns*eta_chns*w2z', ... '-nz_chns*p2+nx_chns*eta_chns*(w2x+u2z)+ny_chns*eta_chns*(w2y+v2z)+2*nz_chns*eta_chns*w2z', ... '-nz_chns*p2+nx_chns*eta_chns*(w2x+u2z)+ny_chns*eta_chns*(w2y+v2z)+2*nz_chns*eta_chns*w2z', ... '-nz_chns*p2+nx_chns*eta_chns*(w2x+u2z)+ny_chns*eta_chns*(w2y+v2z)+2*nz_chns*eta_chns*w2z', ... '','-nz_chns*p2+nx_chns*eta_chns*(w2x+u2z)+ny_chns*eta_chns*(w2y+v2z)+2*nz_chns*eta_chns*w2z'},'FX_smsld',{'','','','','-lm7*dvol_ale/dvol','0','-lm7*dvol_ale/dvol'}, ... 'FXAmp_smsld',{'','','','','1','1','1'}, ... 'FXPh_smsld',{'','','','','0','0','0'}, ... 'FY_smsld',{'','','','','-lm8*dvol_ale/dvol','0','-lm8*dvol_ale/dvol'}, ... 'FYAmp_smsld',{'','','','','1','1','1'}, ... 'FYPh_smsld',{'','','','','0','0','0'}, ... 'FZ_smsld',{'','','','','-lm9*dvol_ale/dvol','0','-lm9*dvol_ale/dvol'}, ... 'FZAmp_smsld',{'','','','','1','1','1'}, ... 'FZPh_smsld',{'','','','','0','0','0'}, ... 'velx_smsld',{'','','','','0','0','0'}, ... 'accx_smsld',{'','','','','0','0','0'}, ... 'vely_smsld',{'','','','','0','0','0'}, ... 'accy_smsld',{'','','','','0','0','0'}, ... 'velz_smsld',{'','','','','0','0','0'}, ... 'accz_smsld',{'','','','','0','0','0'}, ... 'Rx_smsld',{'','','','','0','0','0'}, ... 'Ry_smsld',{'','','','','0','0','0'}, ... 'Rz_smsld',{'','','','','0','0','0'}, ... 'nX_smsld',{'','','','','unX','nX','dnX'}, ... 'nY_smsld',{'','','','','unY','nY','dnY'}, ... 'nZ_smsld',{'','','','','unZ','nZ','dnZ'}, ... 'H11_smsld',{'','','','','0','1','0'}, ... 'H12_smsld',{'','','','','0','0','0'}, ... 'H13_smsld',{'','','','','0','0','0'}, ... 'H21_smsld',{'','','','','0','0','0'}, ... 'H22_smsld',{'','','','','0','1','0'}, ... 'H23_smsld',{'','','','','0','0','0'}, ... 'H31_smsld',{'','','','','0','0','0'}, ... 'H32_smsld',{'','','','','0','0','0'}, ... 'H33_smsld',{'','','','','0','1','0'}, ... 'nX_ale',{'nX','nX','nX','nX','dnX','','unX'}, ... 'nY_ale',{'nY','nY','nY','nY','dnY','','unY'}, ... 'nZ_ale',{'nZ','nZ','nZ','nZ','dnZ','','unZ'}, ... 'u0_chns',{'0.5','0','0','0','0','','0'}, ... 'v0_chns',{'0','0','0','0','0','','0'}, ... 'w0_chns',{'0','0','0','0','0','','0'}, ... 'p0_chns',{'0','0','0','0','0','','0'}, ... 'k0_chns',{'0','0','0','0','0','','0'}, ... 'd0_chns',{'0','0','0','0','0','','0'}, ... 'LT_chns',{'0.07000000000000001','0.07000000000000001','0.07000000000000001', ... '0.07000000000000001','0.07000000000000001','','0.07000000000000001'}, ... 'IT_chns',{'0.054','0.054','0.054','0.054','0.054','','0.054'}, ... 'uw_chns',{'0','0','0','0','0','','0'}, ... 'vw_chns',{'0','0','0','0','0','','0'}, ... 'ww_chns',{'0','0','0','0','0','','0'}, ... 'dw_chns',{'h','h','h','h','h','','h'}, ... 'yplus_chns',{'100','100','100','100','100','','100'}, ... 'U0_chns',{'0','0','0','0','0','','0'}, ... 'V0_chns',{'0','0','0','0','0','','0'}, ... 'p0_inl_chns',{'0','0','0','0','0','','0'}, ... 'L0_chns',{'1','1','1','1','1','','1'}, ... 'nx_chns',{'nx','nx','nx','nx','dnx','','unx'}, ... 'ny_chns',{'ny','ny','ny','ny','dny','','uny'}, ... 'nz_chns',{'nz','nz','nz','nz','dnz','','unz'}}; bnd.expr = {}; bnd.lock = [0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0]; bnd.mlock = {[0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0]}; xfem.bnd = bnd; clear edg edg.weak = {{'0';'0';'0';'0';'0';'0';'0';'0';'0';'0';'0'},{'0';'0';'0';'0'; ... '0';'0';'0';'0';'0';'0';'0'}}; edg.dweak = {{'0';'0';'0';'0';'0';'0';'0';'0';'0';'0';'0'},{'0';'0';'0';'0'; ... '0';'0';'0';'0';'0';'0';'0'}}; edg.constr = {{'0';'0';'0';'0';'0';'0';'0';'0';'0';'0';'0'},{'0';'0';'0';'0'; ... '0';'0';'0';'0';'0';'0';'0'}}; edg.shape = {[1;2;3;4;5;6;7],[8;9;10;1;2;3;4;5;6;7]}; edg.sshape = {[1;2],[1;2]}; edg.sshapedim = {{1;1;1;1;2;2;2;2;2;2;2},{1;1;1;1;2;2;2;2;2;2;2}}; edg.gporder = {{1;1;1;1;1;1;1;1;1;1;2},{1;1;1;1;1;1;1;1;1;1;2}}; edg.cporder = {{1;1;1;1;1;1;1;1;1;1;2},{1;1;1;1;1;1;1;1;1;1;2}}; edg.init = {{'';'';'';'';'';'';'';'';'';'';''},{'';'';'';'';'';'';'';'';''; ... '';''}}; edg.dinit = {{'';'';'';'';'';'';'';'';'';'';''},{'';'';'';'';'';'';'';'';''; ... '';''}}; edg.ind = [1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,2,2,2,2,2,2,2,2,2,2, ... 2,2,1,1,1,1]; edg.dim = {'u','v','w','p','x','y','z','u2','v2','w2','p2'}; edg.var = {'FXg_smsld',{'','0'}, ... 'FYg_smsld',{'','0'}, ... 'FZg_smsld',{'','0'}, ... 'disp_smsld',{'','sqrt(real(u)^2+real(v)^2+real(w)^2)'},'FX_smsld',{'','0'}, ... 'FXAmp_smsld',{'','1'}, ... 'FXPh_smsld',{'','0'}, ... 'FY_smsld',{'','0'}, ... 'FYAmp_smsld',{'','1'}, ... 'FYPh_smsld',{'','0'}, ... 'FZ_smsld',{'','0'}, ... 'FZAmp_smsld',{'','1'}, ... 'FZPh_smsld',{'','0'}, ... 'Rx_smsld',{'','0'}, ... 'Ry_smsld',{'','0'}, ... 'Rz_smsld',{'','0'}, ... 'H11_smsld',{'','0'}, ... 'H12_smsld',{'','0'}, ... 'H13_smsld',{'','0'}, ... 'H21_smsld',{'','0'}, ... 'H22_smsld',{'','0'}, ... 'H23_smsld',{'','0'}, ... 'H31_smsld',{'','0'}, ... 'H32_smsld',{'','0'}, ... 'H33_smsld',{'','0'}}; edg.expr = {}; edg.lock = [0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0]; edg.mlock = {[0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0]}; xfem.edg = edg; clear pnt pnt.weak = {{'0';'0';'0';'0';'0';'0';'0';'0';'0';'0';'0';'0'},{'0';'0';'0'; ... '0';'0';'0';'0';'0';'0';'0';'0';'0'}}; pnt.dweak = {{'0';'0';'0';'0';'0';'0';'0';'0';'0';'0';'0';'0'},{'0';'0';'0'; ... '0';'0';'0';'0';'0';'0';'0';'0';'0'}}; pnt.constr = {{'0';'0';'0';'0';'0';'0';'0';'0';'0';'0';'0';'0'},{'0';'0';'0'; ... '0';'0';'0';'0';'0';'0';'0';'0';'0'}}; pnt.shape = {[1;2;3;4;5;6;7],[8;9;10;1;2;3;4;5;6;7]}; pnt.sshape = {[1;2],[1;2]}; pnt.sshapedim = {{1;1;1;1;2;2;2;2;2;2;2;2},{1;1;1;1;2;2;2;2;2;2;2;2}}; pnt.init = {{'';'';'';'';'';'';'';'';'';'';'';''},{'';'';'';'';'';'';'';''; ... '';'';'';''}}; pnt.dinit = {{'';'';'';'';'';'';'';'';'';'';'';''},{'';'';'';'';'';'';'';''; ... '';'';'';''}}; pnt.ind = [1,1,1,1,1,1,1,1,1,1,1,1,2,2,2,2,2,2,2,2,1,1,1,1]; pnt.dim = {'u','v','w','p','x','y','z','u2','v2','w2','p2','Pinl_chns'}; pnt.var = {'FXg_smsld',{'','0'}, ... 'FYg_smsld',{'','0'}, ... 'FZg_smsld',{'','0'}, ... 'disp_smsld',{'','sqrt(real(u)^2+real(v)^2+real(w)^2)'},'FX_smsld',{'','0'}, ... 'FXAmp_smsld',{'','1'}, ... 'FXPh_smsld',{'','0'}, ... 'FY_smsld',{'','0'}, ... 'FYAmp_smsld',{'','1'}, ... 'FYPh_smsld',{'','0'}, ... 'FZ_smsld',{'','0'}, ... 'FZAmp_smsld',{'','1'}, ... 'FZPh_smsld',{'','0'}, ... 'Rx_smsld',{'','0'}, ... 'Ry_smsld',{'','0'}, ... 'Rz_smsld',{'','0'}, ... 'H11_smsld',{'','0'}, ... 'H12_smsld',{'','0'}, ... 'H13_smsld',{'','0'}, ... 'H21_smsld',{'','0'}, ... 'H22_smsld',{'','0'}, ... 'H23_smsld',{'','0'}, ... 'H31_smsld',{'','0'}, ... 'H32_smsld',{'','0'}, ... 'H33_smsld',{'','0'}, ... 'p0_chns',{'0','0'}}; pnt.expr = {}; pnt.lock = [0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0]; pnt.mlock = {[0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0]}; xfem.pnt = pnt; xfem.var = {'t_old_ini_smsld','-1','ut','uTIME','utt','uTIMETIME','uXt','uXTIME','uYt','uYTIME','uZt','uZTIME','vt','vTIME','vtt','vTIMETIME','vXt','vXTIME','vYt','vYTIME','vZt','vZTIME','wt','wTIME','wtt','wTIMETIME','wXt','wXTIME','wYt','wYTIME','wZt','wZTIME','xt','xTIME','xtt','xTIMETIME','xxt','xxTIME','xyt','xyTIME','xzt','xzTIME','yt','yTIME','ytt','yTIMETIME','yxt','yxTIME','yyt','yyTIME','yzt','yzTIME','zt','zTIME','ztt','zTIMETIME','zxt','zxTIME','zyt','zyTIME','zzt','zzTIME','u2t','u2TIME-xTIME*u2x-yTIME*u2y-zTIME*u2z','u2tt','u2TIMETIME','u2xt','u2xTIME-xTIME*u2xx-yTIME*u2xy-zTIME*u2xz','u2yt','u2yTIME-xTIME*u2yx-yTIME*u2yy-zTIME*u2yz','u2zt','u2zTIME-xTIME*u2zx-yTIME*u2zy-zTIME*u2zz','v2t','v2TIME-xTIME*v2x-yTIME*v2y-zTIME*v2z','v2tt','v2TIMETIME','v2xt','v2xTIME-xTIME*v2xx-yTIME*v2xy-zTIME*v2xz','v2yt','v2yTIME-xTIME*v2yx-yTIME*v2yy-zTIME*v2yz','v2zt','v2zTIME-xTIME*v2zx-yTIME*v2zy-zTIME*v2zz','w2t','w2TIME-xTIME*w2x-yTIME*w2y-zTIME*w2z','w2tt','w2TIMETIME','w2xt','w2xTIME-xTIME*w2xx-yTIME*w2xy-zTIME*w2xz','w2yt','w2yTIME-xTIME*w2yx-yTIME*w2yy-zTIME*w2yz','w2zt','w2zTIME-xTIME*w2zx-yTIME*w2zy-zTIME*w2zz','p2t','p2TIME-xTIME*p2x-yTIME*p2y-zTIME*p2z','p2tt','p2TIMETIME','p2xt','p2xTIME-xTIME*p2xx-yTIME*p2xy-zTIME*p2xz','p2yt','p2yTIME-xTIME*p2yx-yTIME*p2yy-zTIME*p2yz','p2zt','p2zTIME-xTIME*p2zx-yTIME*p2zy-zTIME*p2zz'}; xfem.expr = {}; clear elemmph clear elem elem.elem = 'elpric'; elem.g = {'1'}; clear equ equ.tensor = {{'real(sX_smsld)','real(sXY_smsld)','real(sY_smsld)','real(sXZ_smsld)', ... 'real(sYZ_smsld)','real(sZ_smsld)'}}; equ.ind = {{'3'}}; elem.geomdim = {{{},{},{},equ}}; elem.basename = 's'; elem.postname = 'smsld'; elem.sdim = {'X','Y','Z'}; elemmph{1} = elem; clear elem elem.elem = 'elpric'; elem.g = {'1'}; clear equ equ.tensor = {{'real(eX_smsld)','real(eXY_smsld)','real(eY_smsld)','real(eXZ_smsld)', ... 'real(eYZ_smsld)','real(eZ_smsld)'}}; equ.ind = {{'3'}}; elem.geomdim = {{{},{},{},equ}}; elem.basename = 'e'; elem.postname = 'smsld'; elem.sdim = {'X','Y','Z'}; elemmph{2} = elem; clear elem elem.elem = 'elcplscalar'; elem.g = {'1'}; src = cell(1,1); clear equ equ.expr = {{'-qual_ale'}}; equ.ipoints = {{'4'}}; equ.ind = {{'1','2','3'}}; src{1} = {{},{},{},equ}; elem.src = src; geomdim = cell(1,1); geomdim{1} = {}; elem.geomdim = geomdim; elem.var = {'minqual_tmp1_ale'}; elem.global = {'1'}; elem.maxvars = {'1'}; elemmph{3} = elem; clear elem elem.elem = 'elconst'; elem.var = {'minqual1_ale','-minqual_tmp1_ale'}; elemmph{4} = elem; clear elem elem.elem = 'elcplextr'; elem.g = {'1'}; src = cell(1,1); src{1} = {{}}; elem.src = src; geomdim = cell(1,1); geomdim{1} = {{},{},{},{}}; elem.geomdim = geomdim; elem.var = {}; map = cell(1,0); elem.map = map; elemmph{5} = elem; clear elem elem.elem = 'elcplscalar'; elem.g = {'1'}; src = cell(1,1); src{1} = {{},{},{},{}}; elem.src = src; geomdim = cell(1,1); geomdim{1} = {{},{},{},{}}; elem.geomdim = geomdim; elem.var = {}; elem.global = {}; elemmph{6} = elem; xfem.elemmph = elemmph; clear draw draw.p.objs = {}; draw.p.name = {}; draw.c.objs = {}; draw.c.name = {}; draw.f.objs = {}; draw.f.name = {}; draw.s.objs = {flbinary('g7','draw','tutorial5.mph'),flbinary('g8','draw','tutorial5.mph')}; draw.s.name = {'CO1','BLK1'}; xfem.draw = draw; xfem.const = {}; xfem.globalexpr = {}; clear fcns xfem.functions = {}; xfem.sol = flbinary('xfem','solution','tutorial5.mph'); xfem.xmcases = [0]; xfem.mcases = [0]; flbinary clear; xfem.rulingmode = 'smsld'; xfem.solform = 'weak'; uq ~ t��gui.solvemodel.toutcomp='off'; gui.solvemodel.currsolver='stationary'; gui.solvemodel.solveroption=''; gui.solvemodel.postsolver='stationary'; gui.solvemodel.nonlin='auto'; gui.solvemodel.ntol='1.0E-6'; gui.solvemodel.maxiter='25'; gui.solvemodel.manualdamp='off'; gui.solvemodel.hnlin='off'; gui.solvemodel.initstep='1.0'; gui.solvemodel.minstep='1.0E-4'; gui.solvemodel.rstep='10.0'; gui.solvemodel.useaugsolver='off'; gui.solvemodel.autoaugcomp='on'; gui.solvemodel.augcomp=''; gui.solvemodel.augtol='1.0E-6'; gui.solvemodel.augmaxiter='25'; gui.solvemodel.augsolver='lumped'; gui.solvemodel.atol='0.0010'; gui.solvemodel.rtol='0.01'; gui.solvemodel.tlist='0:0.1:1'; gui.solvemodel.tout='tlist'; gui.solvemodel.tsteps='free'; gui.solvemodel.manualreassem='off'; gui.solvemodel.emassconst='on'; gui.solvemodel.massconst='on'; gui.solvemodel.loadconst='on'; gui.solvemodel.constrconst='on'; gui.solvemodel.jacobianconst='on'; gui.solvemodel.constrjacobianconst='on'; gui.solvemodel.manualstep='off'; gui.solvemodel.initialstep='0.0010'; gui.solvemodel.maxorder='5'; gui.solvemodel.maxstep='1.0'; gui.solvemodel.timeusestopcond='off'; gui.solvemodel.paramusestopcond='off'; gui.solvemodel.masssingular='maybe'; gui.solvemodel.consistent='bweuler'; gui.solvemodel.estrat='0'; gui.solvemodel.complex='off'; gui.solvemodel.neigs='6'; gui.solvemodel.shift='0'; gui.solvemodel.maxeigit='300'; gui.solvemodel.etol='0.0'; gui.solvemodel.krylovdim='0'; gui.solvemodel.eigname='lambda'; gui.solvemodel.eigref='0'; gui.solvemodel.pname=''; gui.solvemodel.plist=''; gui.solvemodel.porder='1'; gui.solvemodel.manualparam='off'; gui.solvemodel.pinitstep='0.0'; gui.solvemodel.pminstep='0.0'; gui.solvemodel.pmaxstep='0.0'; gui.solvemodel.autooldcomp='on'; gui.solvemodel.oldcomp=''; gui.solvemodel.outform='auto'; gui.solvemodel.symmetric='auto'; gui.solvemodel.symmhermit='auto'; gui.solvemodel.method='eliminate'; gui.solvemodel.nullfun='auto'; gui.solvemodel.blocksize='5000'; gui.solvemodel.uscale='auto'; gui.solvemodel.manscale=''; gui.solvemodel.rowscale='on'; gui.solvemodel.conjugate='off'; gui.solvemodel.complexfun='off'; gui.solvemodel.matherr='on'; gui.solvemodel.solfile='off'; gui.solvemodel.adaptgeom='currgeom'; gui.solvemodel.maxt='10000000'; gui.solvemodel.rmethod='longest'; gui.solvemodel.resmethod='weak'; gui.solvemodel.resorder='0'; gui.solvemodel.l2scale='1'; gui.solvemodel.l2staborder='2'; gui.solvemodel.eigselect='1'; gui.solvemodel.tpfun='fltpft'; gui.solvemodel.ngen='1'; gui.solvemodel.tpmult='1.1'; gui.solvemodel.tpworst='0.8'; gui.solvemodel.tpfract='0.2'; gui.solvemodel.autosolver='on'; gui.solvemodel.solcomp='z,lm5,u2,w,lm9,p2,lm8,lm4,v,lm7,u,w2,y,v2,x,lm6'; gui.solvemodel.outcomp='z,lm5,u2,w,lm9,p2,lm8,lm4,v,lm7,u,w2,y,v2,x,lm6'; gui.solvemodel.inittype='init_expr_currsol_radio'; gui.solvemodel.initsolnum='Automatic'; gui.solvemodel.inittime='0'; gui.solvemodel.utype='u_init_radio'; gui.solvemodel.usolnum='Automatic'; gui.solvemodel.utime='0'; gui.solvemodel.scriptcommands=''; gui.solvemodel.usescript='off'; gui.solvemodel.autoscript='off'; gui.solvemodel.sameaxis='off'; gui.solvemodel.linsolvernode.currlinsolver='umfpack'; gui.solvemodel.linsolvernode.type='linsolver'; gui.solvemodel.linsolvernode.droptol='0.0'; gui.solvemodel.linsolvernode.thresh='0.1'; gui.solvemodel.linsolvernode.umfalloc='0.7'; gui.solvemodel.linsolvernode.preorder='mmd'; gui.solvemodel.linsolvernode.preroworder='on'; gui.solvemodel.linsolvernode.pivotperturb='1.0E-8'; gui.solvemodel.linsolvernode.errorchk='on'; gui.solvemodel.linsolvernode.iter='2'; gui.solvemodel.linsolvernode.itol='1.0E-6'; gui.solvemodel.linsolvernode.rhob='400.0'; gui.solvemodel.linsolvernode.maxlinit='10000'; gui.solvemodel.linsolvernode.itrestart='50'; gui.solvemodel.linsolvernode.seconditer='1'; gui.solvemodel.linsolvernode.relax='1.0'; gui.solvemodel.linsolvernode.amgauto='3'; gui.solvemodel.linsolvernode.mglevels='6'; gui.solvemodel.linsolvernode.mgcycle='v'; gui.solvemodel.linsolvernode.maxcoarsedof='5000'; gui.solvemodel.linsolvernode.oocmemory='128.0'; gui.solvemodel.linsolvernode.oocfilename=''; gui.solvemodel.linsolvernode.modified='off'; gui.solvemodel.linsolvernode.fillratio='2.0'; gui.solvemodel.linsolvernode.respectpattern='on'; gui.solvemodel.linsolvernode.droptype='droptol'; gui.solvemodel.linsolvernode.vankavars=''; gui.solvemodel.linsolvernode.vankasolv='gmres'; gui.solvemodel.linsolvernode.vankatol='0.02'; gui.solvemodel.linsolvernode.vankarestart='100'; gui.solvemodel.linsolvernode.vankarelax='0.8'; gui.solvemodel.linsolvernode.mgauto='shape'; gui.solvemodel.linsolvernode.rmethod='regular'; gui.solvemodel.linsolvernode.coarseassem='on'; gui.solvemodel.linsolvernode.meshscale='2'; gui.solvemodel.linsolvernode.mgautolevels='2'; gui.solvemodel.linsolvernode.mgkeep='off'; gui.solvemodel.linsolvernode.mggeom='Geom1'; gui.solvemodel.linsolvernode.mcase0='on'; gui.solvemodel.linsolvernode.mgassem0='on'; gui.solvemodel.defaults.prefun='gmg'; gui.solvemodel.defaults.linsolver='spooles'; gui.registry.general_currentmodel='Geom1'; gui.registry.general_currmeshcase='0'; gui.registry.general_savedonserver='off'; gui.registry.general_savedchanges='off'; gui.registry.general_rulingmode=''; gui.registry.general_incompletemfilehistory='off'; gui.registry.saved_license='1027370'; gui.registry.saved_version='COMSOL 3.3.0.405'; gui.registry.info_modelname=''; gui.registry.info_author=''; gui.registry.info_company=''; gui.registry.info_department=''; gui.registry.info_reference=''; gui.registry.info_url=''; gui.registry.info_saveddate='1278450746734'; gui.registry.info_creationdate='1278433404140'; gui.registry.info_modelresult=''; gui.reportregistry.report_contents=''; gui.reportregistry.report_outputformat='html'; gui.reportregistry.report_filename=''; gui.reportregistry.report_autoopen='off'; gui.reportregistry.report_paperformat='a4'; gui.reportregistry.report_includedefaults='off'; gui.reportregistry.report_template='full'; gui.reportregistry.report_showemptysections='off'; gui.flmodel{1}.modelname='Geom1'; gui.flmodel{1}.currmode='edge'; gui.flmodel{1}.currappl='2'; gui.flmodel{1}.axis.xmin='-9.999999747378752E-6'; gui.flmodel{1}.axis.xmax='1.9999999494757503E-4'; gui.flmodel{1}.axis.ymin='0.0'; gui.flmodel{1}.axis.ymax='1.9999999494757503E-4'; gui.flmodel{1}.axis.zmin='0.0'; gui.flmodel{1}.axis.zmax='1.9999999494757503E-4'; gui.flmodel{1}.axis.xspacing='4.9999999999999996E-5'; gui.flmodel{1}.axis.yspacing='4.9999999999999996E-5'; gui.flmodel{1}.axis.zspacing='4.9999999999999996E-5'; gui.flmodel{1}.axis.extrax=''; gui.flmodel{1}.axis.extray=''; gui.flmodel{1}.axis.extraz=''; gui.flmodel{1}.camera.xmin='-2.10234396334303'; gui.flmodel{1}.camera.xmax='2.10234396334303'; gui.flmodel{1}.camera.ymin='-1.8271196161067038'; gui.flmodel{1}.camera.ymax='1.8271196161067038'; gui.flmodel{1}.camera.camposx='-10.307033115682344'; gui.flmodel{1}.camera.camposy='13.410930427848111'; gui.flmodel{1}.camera.camposz='-4.431205746436035'; gui.flmodel{1}.camera.camtargetx='0.9499999760009814'; gui.flmodel{1}.camera.camtargety='0.9999999747378752'; gui.flmodel{1}.camera.camtargetz='0.9999999747378752'; gui.flmodel{1}.camera.camupx='0.22253243627166236'; gui.flmodel{1}.camera.camupy='-0.2143585613779249'; gui.flmodel{1}.camera.camupz='-0.9510676747587338'; gui.flmodel{1}.lightmodel.headlight.type='point'; gui.flmodel{1}.lightmodel.headlight.name='headlight'; gui.flmodel{1}.lightmodel.headlight.enable='on'; gui.flmodel{1}.lightmodel.headlight.colorr='255'; gui.flmodel{1}.lightmodel.headlight.colorg='255'; gui.flmodel{1}.lightmodel.headlight.colorb='255'; gui.flmodel{1}.lightmodel.headlight.xpos='-10.307034'; gui.flmodel{1}.lightmodel.headlight.ypos='13.410931'; gui.flmodel{1}.lightmodel.headlight.zpos='-4.4312057'; gui.flmodel{1}.lightmodel.scenelight{1}.type='directional'; gui.flmodel{1}.lightmodel.scenelight{1}.name='light 1'; gui.flmodel{1}.lightmodel.scenelight{1}.enable='on'; gui.flmodel{1}.lightmodel.scenelight{1}.colorr='255'; gui.flmodel{1}.lightmodel.scenelight{1}.colorg='255'; gui.flmodel{1}.lightmodel.scenelight{1}.colorb='255'; gui.flmodel{1}.lightmodel.scenelight{1}.xdir='1.0'; gui.flmodel{1}.lightmodel.scenelight{1}.ydir='-1.0'; gui.flmodel{1}.lightmodel.scenelight{1}.zdir='1.0'; gui.flmodel{1}.lightmodel.scenelight{2}.type='directional'; gui.flmodel{1}.lightmodel.scenelight{2}.name='light 2'; gui.flmodel{1}.lightmodel.scenelight{2}.enable='on'; gui.flmodel{1}.lightmodel.scenelight{2}.colorr='255'; gui.flmodel{1}.lightmodel.scenelight{2}.colorg='255'; gui.flmodel{1}.lightmodel.scenelight{2}.colorb='255'; gui.flmodel{1}.lightmodel.scenelight{2}.xdir='-1.0'; gui.flmodel{1}.lightmodel.scenelight{2}.ydir='1.0'; gui.flmodel{1}.lightmodel.scenelight{2}.zdir='-1.0'; gui.flmodel{1}.lightmodel.scenelight{3}.type='directional'; gui.flmodel{1}.lightmodel.scenelight{3}.name='light 3'; gui.flmodel{1}.lightmodel.scenelight{3}.enable='on'; gui.flmodel{1}.lightmodel.scenelight{3}.colorr='255'; gui.flmodel{1}.lightmodel.scenelight{3}.colorg='255'; gui.flmodel{1}.lightmodel.scenelight{3}.colorb='255'; gui.flmodel{1}.lightmodel.scenelight{3}.xdir='-1.0'; gui.flmodel{1}.lightmodel.scenelight{3}.ydir='-1.0'; gui.flmodel{1}.lightmodel.scenelight{3}.zdir='1.0'; gui.flmodel{1}.lightmodel.scenelight{4}.type='directional'; gui.flmodel{1}.lightmodel.scenelight{4}.name='light 4'; gui.flmodel{1}.lightmodel.scenelight{4}.enable='on'; gui.flmodel{1}.lightmodel.scenelight{4}.colorr='255'; gui.flmodel{1}.lightmodel.scenelight{4}.colorg='255'; gui.flmodel{1}.lightmodel.scenelight{4}.colorb='255'; gui.flmodel{1}.lightmodel.scenelight{4}.xdir='1.0'; gui.flmodel{1}.lightmodel.scenelight{4}.ydir='-1.0'; gui.flmodel{1}.lightmodel.scenelight{4}.zdir='-1.0'; gui.flmodel{1}.registry.axis_visible='on'; gui.flmodel{1}.registry.axis_auto='on'; gui.flmodel{1}.registry.axis_autoy='on'; gui.flmodel{1}.registry.axis_autoz='on'; gui.flmodel{1}.registry.axis_box='off'; gui.flmodel{1}.registry.axis_equal='on'; gui.flmodel{1}.registry.axis_csys='on'; gui.flmodel{1}.registry.grid_visible='on'; gui.flmodel{1}.registry.grid_auto='on'; gui.flmodel{1}.registry.grid_autoz='on'; gui.flmodel{1}.registry.grid_labels='on'; gui.flmodel{1}.registry.labels_object='on'; gui.flmodel{1}.registry.labels_vertex='off'; gui.flmodel{1}.registry.labels_edge='off'; gui.flmodel{1}.registry.labels_face='off'; gui.flmodel{1}.registry.labels_subdomain='off'; gui.flmodel{1}.registry.symbols_vertexscale='1.0'; gui.flmodel{1}.registry.symbols_edgescale='1.0'; gui.flmodel{1}.registry.symbols_facescale='1.0'; gui.flmodel{1}.registry.select_draw2d='size'; gui.flmodel{1}.registry.select_adj='cycle'; gui.flmodel{1}.registry.light_headlight='off'; gui.flmodel{1}.registry.light_scenelight='off'; gui.flmodel{1}.registry.light_shininess='0.5'; gui.flmodel{1}.registry.camera_mouse='orbit'; gui.flmodel{1}.registry.camera_camconstr='none'; gui.flmodel{1}.registry.camera_mouseconstr='none'; gui.flmodel{1}.registry.camera_perspective='off'; gui.flmodel{1}.registry.camera_moveasbox='off'; gui.flmodel{1}.registry.draw_assembly='off'; gui.flmodel{1}.registry.draw_dialog='off'; gui.flmodel{1}.registry.draw_keepborders='on'; gui.flmodel{1}.registry.draw_keepedges='off'; gui.flmodel{1}.registry.draw_multi='off'; gui.flmodel{1}.registry.draw_snap2grid='on'; gui.flmodel{1}.registry.draw_snap2vtx='on'; gui.flmodel{1}.registry.draw_solid='on'; gui.flmodel{1}.registry.draw_workplane_coordsys='on'; gui.flmodel{1}.registry.draw_workplane_showgeom='on'; gui.flmodel{1}.registry.draw_repair='off'; gui.flmodel{1}.registry.draw_repairtol=''; gui.flmodel{1}.registry.draw_projection='intersection'; gui.flmodel{1}.registry.transparency_value='1.0'; gui.flmodel{1}.registry.mesh_geomdetail='normal'; gui.flmodel{1}.registry.mesh_showquality='off'; gui.flmodel{1}.registry.post_cameraview='2'; gui.flmodel{1}.registry.graphics_scale='10000.0'; gui.flmodel{1}.registry.render_mesh='on'; gui.flmodel{1}.registry.render_bndarrow='on'; gui.flmodel{1}.registry.render_vertex='off'; gui.flmodel{1}.registry.render_edge='on'; gui.flmodel{1}.registry.render_face='off'; gui.flmodel{1}.registry.highlight_vertex='off'; gui.flmodel{1}.registry.highlight_edge='on'; gui.flmodel{1}.registry.highlight_face='on'; gui.flmodel{1}.meshparam.hauto='7'; gui.flmodel{1}.meshparam.usehauto='on'; gui.flmodel{1}.meshparam.hmax=''; gui.flmodel{1}.meshparam.hmaxfact='5'; gui.flmodel{1}.meshparam.hcurve='1'; gui.flmodel{1}.meshparam.hgrad='2'; gui.flmodel{1}.meshparam.hcutoff='0.07'; gui.flmodel{1}.meshparam.hnarrow='0.1'; gui.flmodel{1}.meshparam.hpnt='20'; gui.flmodel{1}.meshparam.xscale='1.0'; gui.flmodel{1}.meshparam.yscale='1.0'; gui.flmodel{1}.meshparam.zscale='1.0'; gui.flmodel{1}.meshparam.jiggle='on'; gui.flmodel{1}.meshparam.mcase='0'; gui.flmodel{1}.meshparam.rmethod='longest'; gui.flmodel{1}.meshparam.hmaxvtx={'','','','','','','','','','','','','','','','','','','','','','','',''}; gui.flmodel{1}.meshparam.hgradvtx={'','','','','','','','','','','','','','','','','','','','','','','',''}; gui.flmodel{1}.meshparam.hgradedg={'','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','',''}; gui.flmodel{1}.meshparam.hcurveedg={'','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','',''}; gui.flmodel{1}.meshparam.hmaxedg={'','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','',''}; gui.flmodel{1}.meshparam.hcutoffedg={'','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','',''}; gui.flmodel{1}.meshparam.hmaxfac={'','','','','','','','','','','','','','','','','',''}; gui.flmodel{1}.meshparam.hcutofffac={'','','','','','','','','','','','','','','','','',''}; gui.flmodel{1}.meshparam.hcurvefac={'','','','','','','','','','','','','','','','','',''}; gui.flmodel{1}.meshparam.hgradfac={'','','','','','','','','','','','','','','','','',''}; gui.flmodel{1}.meshparam.hgradsub={'','',''}; gui.flmodel{1}.meshparam.hmaxsub={'','',''}; gui.flmodel{1}.mapmeshparam.edgegroups={{'2','6','4','1'},{'5','14','3','1'},{'15','2','3','7'},{'16','4','5','8'},{'8','17','7','6'},{'10','12','9','18'},{'13','33','11','9'},{'','','',''},{'','','',''},{},{},{'22','26','24','21'},{'25','29','23','21'},{},{},{},{},{}}; gui.flmodel{1}.mapmeshparam.hauto=7gui.flmodel{1}.sweptmeshparam.nsubdomains=3; gui.flmodel{1}.sweptmeshparam.sourceface={1,'2',2,'8',3,'14'}; gui.flmodel{1}.sweptmeshparam.targetface={1,'5',2,'9',3,'15'}; gui.flmodel{1}.sweptmeshparam.hauto=-1gui.flmodel{1}.postmodel.postplot.sliceplot='on'; gui.flmodel{1}.postmodel.postplot.slicedata={'mises_smsld'}; gui.flmodel{1}.postmodel.postplot.slicerangeauto='on'; gui.flmodel{1}.postmodel.postplot.slicerangemin='0.044799589399362766'; gui.flmodel{1}.postmodel.postplot.slicerangemax='6.214998989201082'; gui.flmodel{1}.postmodel.postplot.slicecont='on'; gui.flmodel{1}.postmodel.postplot.sliceunit='Pa'; gui.flmodel{1}.postmodel.postplot.slicexspacing='5'; gui.flmodel{1}.postmodel.postplot.slicexvectorspacing=''; gui.flmodel{1}.postmodel.postplot.slicexisvector='off'; gui.flmodel{1}.postmodel.postplot.sliceyspacing='0'; gui.flmodel{1}.postmodel.postplot.sliceyvectorspacing=''; gui.flmodel{1}.postmodel.postplot.sliceyisvector='off'; gui.flmodel{1}.postmodel.postplot.slicezspacing='0'; gui.flmodel{1}.postmodel.postplot.slicezvectorspacing=''; gui.flmodel{1}.postmodel.postplot.slicezisvector='off'; gui.flmodel{1}.postmodel.postplot.slicemap='jet'; gui.flmodel{1}.postmodel.postplot.slicemapdepth='1024'; gui.flmodel{1}.postmodel.postplot.slicebar='on'; gui.flmodel{1}.postmodel.postplot.sliceusemap='on'; gui.flmodel{1}.postmodel.postplot.slicecolor='255,0,0'; gui.flmodel{1}.postmodel.postplot.slicecoloring='interp'; gui.flmodel{1}.postmodel.postplot.slicefill='fill'; gui.flmodel{1}.postmodel.postplot.isoplot='off'; gui.flmodel{1}.postmodel.postplot.isodata={'mises_smsld'}; gui.flmodel{1}.postmodel.postplot.isocont='on'; gui.flmodel{1}.postmodel.postplot.isounit='Pa'; gui.flmodel{1}.postmodel.postplot.isolevels='5'; gui.flmodel{1}.postmodel.postplot.isovectorlevels=''; gui.flmodel{1}.postmodel.postplot.isoisvector='off'; gui.flmodel{1}.postmodel.postplot.isocolordata={'mises_smsld'}; gui.flmodel{1}.postmodel.postplot.isocolorrangeauto='on'; gui.flmodel{1}.postmodel.postplot.isocolorrangemin=''; gui.flmodel{1}.postmodel.postplot.isocolorrangemax=''; gui.flmodel{1}.postmodel.postplot.isocolordatacheck='off'; gui.flmodel{1}.postmodel.postplot.isocolorunit='Pa'; gui.flmodel{1}.postmodel.postplot.isomap='jet'; gui.flmodel{1}.postmodel.postplot.isomapdepth='1024'; gui.flmodel{1}.postmodel.postplot.isobar='on'; gui.flmodel{1}.postmodel.postplot.isousemap='on'; gui.flmodel{1}.postmodel.postplot.isocolor='255,0,0'; gui.flmodel{1}.postmodel.postplot.isocoloring='interp'; gui.flmodel{1}.postmodel.postplot.isofill='fill'; gui.flmodel{1}.postmodel.postplot.tetplot='off'; gui.flmodel{1}.postmodel.postplot.tetkeepfraction='1'; gui.flmodel{1}.postmodel.postplot.tetkeeptype='random'; gui.flmodel{1}.postmodel.postplot.tetdata={'mises_smsld'}; gui.flmodel{1}.postmodel.postplot.tetrangeauto='on'; gui.flmodel{1}.postmodel.postplot.tetrangemin=''; gui.flmodel{1}.postmodel.postplot.tetrangemax=''; gui.flmodel{1}.postmodel.postplot.tetcont='on'; gui.flmodel{1}.postmodel.postplot.tetunit='Pa'; gui.flmodel{1}.postmodel.postplot.tetmap='jet'; gui.flmodel{1}.postmodel.postplot.tetmapdepth='1024'; gui.flmodel{1}.postmodel.postplot.tetbar='on'; gui.flmodel{1}.postmodel.postplot.tetusemap='on'; gui.flmodel{1}.postmodel.postplot.tetcolor='255,0,0'; gui.flmodel{1}.postmodel.postplot.tetcoloring='interp'; gui.flmodel{1}.postmodel.postplot.tetfill='fill'; gui.flmodel{1}.postmodel.postplot.triplot='on'; gui.flmodel{1}.postmodel.postplot.tridata={'disp_smsld'}; gui.flmodel{1}.postmodel.postplot.trirangeauto='on'; gui.flmodel{1}.postmodel.postplot.trirangemin='0.0'; gui.flmodel{1}.postmodel.postplot.trirangemax='6.03800714871271E-9'; gui.flmodel{1}.postmodel.postplot.tricont='on'; gui.flmodel{1}.postmodel.postplot.triunit='m'; gui.flmodel{1}.postmodel.postplot.trimap='jet'; gui.flmodel{1}.postmodel.postplot.trimapdepth='1024'; gui.flmodel{1}.postmodel.postplot.tribar='on'; gui.flmodel{1}.postmodel.postplot.triusemap='on'; gui.flmodel{1}.postmodel.postplot.tricolor='255,0,0'; gui.flmodel{1}.postmodel.postplot.tricoloring='interp'; gui.flmodel{1}.postmodel.postplot.trifill='fill'; gui.flmodel{1}.postmodel.postplot.linplot='off'; gui.flmodel{1}.postmodel.postplot.lindata={'disp_smsld'}; gui.flmodel{1}.postmodel.postplot.linrangeauto='on'; gui.flmodel{1}.postmodel.postplot.linrangemin=''; gui.flmodel{1}.postmodel.postplot.linrangemax=''; gui.flmodel{1}.postmodel.postplot.lincont='on'; gui.flmodel{1}.postmodel.postplot.linunit='m'; gui.flmodel{1}.postmodel.postplot.linmap='jet'; gui.flmodel{1}.postmodel.postplot.linmapdepth='1024'; gui.flmodel{1}.postmodel.postplot.linbar='on'; gui.flmodel{1}.postmodel.postplot.linusemap='on'; gui.flmodel{1}.postmodel.postplot.lincolor='255,0,0'; gui.flmodel{1}.postmodel.postplot.lincoloring='interp'; gui.flmodel{1}.postmodel.postplot.arrowplot='off'; gui.flmodel{1}.postmodel.postplot.arrowploton='sub'; gui.flmodel{1}.postmodel.postplot.arrowdata={'u','v','w'}; gui.flmodel{1}.postmodel.postplot.arrowunit='m'; gui.flmodel{1}.postmodel.postplot.arrowbnddata={'u','v','w'}; gui.flmodel{1}.postmodel.postplot.arrowbndunit='m'; gui.flmodel{1}.postmodel.postplot.arrowedgdata={'u','v','w'}; gui.flmodel{1}.postmodel.postplot.arrowedgunit='m'; gui.flmodel{1}.postmodel.postplot.arrowxspacing='7'; gui.flmodel{1}.postmodel.postplot.arrowxvectorspacing=''; gui.flmodel{1}.postmodel.postplot.arrowxisvector='off'; gui.flmodel{1}.postmodel.postplot.arrowyspacing='7'; gui.flmodel{1}.postmodel.postplot.arrowyvectorspacing=''; gui.flmodel{1}.postmodel.postplot.arrowyisvector='off'; gui.flmodel{1}.postmodel.postplot.arrowzspacing='7'; gui.flmodel{1}.postmodel.postplot.arrowzvectorspacing=''; gui.flmodel{1}.postmodel.postplot.arrowzisvector='off'; gui.flmodel{1}.postmodel.postplot.arrowtype='cone'; gui.flmodel{1}.postmodel.postplot.arrowlength='proportional'; gui.flmodel{1}.postmodel.postplot.arrowcolor='255,0,0'; gui.flmodel{1}.postmodel.postplot.arrowautoscale='on'; gui.flmodel{1}.postmodel.postplot.arrowscale='1'; gui.flmodel{1}.postmodel.postplot.princplot='off'; gui.flmodel{1}.postmodel.postplot.princploton='sub'; gui.flmodel{1}.postmodel.postplot.princdata={'s1_smsld','s1X_smsld','s1Y_smsld','s1Z_smsld','s2_smsld','s2X_smsld','s2Y_smsld','s2Z_smsld','s3_smsld','s3X_smsld','s3Y_smsld','s3Z_smsld'}; gui.flmodel{1}.postmodel.postplot.princunit='Pa'; gui.flmodel{1}.postmodel.postplot.princbnddata={'s1_smsld','s1X_smsld','s1Y_smsld','s1Z_smsld','s2_smsld','s2X_smsld','s2Y_smsld','s2Z_smsld','s3_smsld','s3X_smsld','s3Y_smsld','s3Z_smsld'}; gui.flmodel{1}.postmodel.postplot.princbndunit='Pa'; gui.flmodel{1}.postmodel.postplot.princxspacing='5'; gui.flmodel{1}.postmodel.postplot.princxvectorspacing=''; gui.flmodel{1}.postmodel.postplot.princxisvector='off'; gui.flmodel{1}.postmodel.postplot.princyspacing='5'; gui.flmodel{1}.postmodel.postplot.princyvectorspacing=''; gui.flmodel{1}.postmodel.postplot.princyisvector='off'; gui.flmodel{1}.postmodel.postplot.princzspacing='5'; gui.flmodel{1}.postmodel.postplot.princzvectorspacing=''; gui.flmodel{1}.postmodel.postplot.princzisvector='off'; gui.flmodel{1}.postmodel.postplot.princtype='cone'; gui.flmodel{1}.postmodel.postplot.princlength='proportional'; gui.flmodel{1}.postmodel.postplot.princcolor='0,153,0'; gui.flmodel{1}.postmodel.postplot.princautoscale='on'; gui.flmodel{1}.postmodel.postplot.princscale='1'; gui.flmodel{1}.postmodel.postplot.flowplot='off'; gui.flmodel{1}.postmodel.postplot.flowdata={'u2','v2','w2'}; gui.flmodel{1}.postmodel.postplot.flowunit='m/s'; gui.flmodel{1}.postmodel.postplot.flowuseexpression='off'; gui.flmodel{1}.postmodel.postplot.flowcolor='255,0,0'; gui.flmodel{1}.postmodel.postplot.flowcolordata={'mises_smsld'}; gui.flmodel{1}.postmodel.postplot.flowcolorunit='Pa'; gui.flmodel{1}.postmodel.postplot.flowmap='jet'; gui.flmodel{1}.postmodel.postplot.flowmapdepth='1024'; gui.flmodel{1}.postmodel.postplot.flowbar='on'; gui.flmodel{1}.postmodel.postplot.flowlinetype='line'; gui.flmodel{1}.postmodel.postplot.flowradiusdata={'mises_smsld'}; gui.flmodel{1}.postmodel.postplot.flowradiusdatacheck='off'; gui.flmodel{1}.postmodel.postplot.flowradiusunit='Pa'; gui.flmodel{1}.postmodel.postplot.flowradiusautoscale='on'; gui.flmodel{1}.postmodel.postplot.flowradiusscale='1'; gui.flmodel{1}.postmodel.postplot.flowlines='20'; gui.flmodel{1}.postmodel.postplot.flowstart='sub'; gui.flmodel{1}.postmodel.postplot.flowstartx=''; gui.flmodel{1}.postmodel.postplot.flowstarty=''; gui.flmodel{1}.postmodel.postplot.flowstartz=''; gui.flmodel{1}.postmodel.postplot.flowisstartvector='off'; gui.flmodel{1}.postmodel.postplot.flowtol='0.001'; gui.flmodel{1}.postmodel.postplot.flowstattol='0.01'; gui.flmodel{1}.postmodel.postplot.flowlooptol='0.01'; gui.flmodel{1}.postmodel.postplot.flowmaxtime='Inf'; gui.flmodel{1}.postmodel.postplot.flowmaxsteps='5000'; gui.flmodel{1}.postmodel.postplot.flowback='on'; gui.flmodel{1}.postmodel.postplot.flownormal='off'; gui.flmodel{1}.postmodel.postplot.flowdistuniform='0.15'; gui.flmodel{1}.postmodel.postplot.flowdistmin='0.05'; gui.flmodel{1}.postmodel.postplot.flowdistmax='0.15'; gui.flmodel{1}.postmodel.postplot.flowseedmanual='off'; gui.flmodel{1}.postmodel.postplot.flowseed1=''; gui.flmodel{1}.postmodel.postplot.flowseed2=''; gui.flmodel{1}.postmodel.postplot.flowseed3=''; gui.flmodel{1}.postmodel.postplot.flowinitref='1'; gui.flmodel{1}.postmodel.postplot.flowignoredist='0.5'; gui.flmodel{1}.postmodel.postplot.flowsat='1.3'; gui.flmodel{1}.postmodel.postplot.flowdistend='0.5'; gui.flmodel{1}.postmodel.postplot.flowdens='none'; gui.flmodel{1}.postmodel.postplot.partplot='off'; gui.flmodel{1}.postmodel.postplot.partmasstype='mass'; gui.flmodel{1}.postmodel.postplot.partplotas='lines'; gui.flmodel{1}.postmodel.postplot.predefforces='khanforce_chns'; gui.flmodel{1}.postmodel.postplot.partforceparam0='radius,partr,1e-4'; gui.flmodel{1}.postmodel.postplot.partmass='4*pi/3*1e-9'; gui.flmodel{1}.postmodel.postplot.partforce={'','',''}; gui.flmodel{1}.postmodel.postplot.part_massless_flowdata={'u','v','w'}; gui.flmodel{1}.postmodel.postplot.part_massless_flowunit='m'; gui.flmodel{1}.postmodel.postplot.parttstartauto='on'; gui.flmodel{1}.postmodel.postplot.parttstart=''; gui.flmodel{1}.postmodel.postplot.partvelstart={'0','0','0'}; gui.flmodel{1}.postmodel.postplot.explicitcoord={'0','0','0'}; gui.flmodel{1}.postmodel.postplot.partuseexpression='off'; gui.flmodel{1}.postmodel.postplot.partcolor='255,0,0'; gui.flmodel{1}.postmodel.postplot.partcolordata={'mises_smsld'}; gui.flmodel{1}.postmodel.postplot.partmap='jet'; gui.flmodel{1}.postmodel.postplot.partmapdepth='1024'; gui.flmodel{1}.postmodel.postplot.partbar='on'; gui.flmodel{1}.postmodel.postplot.partpointcolor='255,0,0'; gui.flmodel{1}.postmodel.postplot.partpointautoscale='on'; gui.flmodel{1}.postmodel.postplot.partpointscale='1'; gui.flmodel{1}.postmodel.postplot.partdroptype='once'; gui.flmodel{1}.postmodel.postplot.partdroptimes=''; gui.flmodel{1}.postmodel.postplot.partdropfreq=''; gui.flmodel{1}.postmodel.postplot.partbnd='stick'; gui.flmodel{1}.postmodel.postplot.partlinetype='line'; gui.flmodel{1}.postmodel.postplot.partradiusdata={'mises_smsld'}; gui.flmodel{1}.postmodel.postplot.partradiusdatacheck='off'; gui.flmodel{1}.postmodel.postplot.partradiusunit='Pa'; gui.flmodel{1}.postmodel.postplot.partradiusautoscale='on'; gui.flmodel{1}.postmodel.postplot.partradiusscale='1'; gui.flmodel{1}.postmodel.postplot.partmasslessrtol='0.001'; gui.flmodel{1}.postmodel.postplot.partmasslessatolmanual='off'; gui.flmodel{1}.postmodel.postplot.partmasslessatol={''}; gui.flmodel{1}.postmodel.postplot.partmasslessstepsizemanual='off'; gui.flmodel{1}.postmodel.postplot.partmasslesstendauto='on'; gui.flmodel{1}.postmodel.postplot.partmasslessmaxstepsauto='on'; gui.flmodel{1}.postmodel.postplot.partmasslessedgetol='0.001'; gui.flmodel{1}.postmodel.postplot.partmasslesstvar='partt'; gui.flmodel{1}.postmodel.postplot.partmasslessstatic='off'; gui.flmodel{1}.postmodel.postplot.partmasslessres='5'; gui.flmodel{1}.postmodel.postplot.partrtol='0.001'; gui.flmodel{1}.postmodel.postplot.partatolmanual='off'; gui.flmodel{1}.postmodel.postplot.partatol={'',''}; gui.flmodel{1}.postmodel.postplot.partstepsizemanual='off'; gui.flmodel{1}.postmodel.postplot.parttendauto='on'; gui.flmodel{1}.postmodel.postplot.partmaxstepsauto='on'; gui.flmodel{1}.postmodel.postplot.partedgetol='0.001'; gui.flmodel{1}.postmodel.postplot.partvelvar={'partu','partv','partw'}; gui.flmodel{1}.postmodel.postplot.parttvar='partt'; gui.flmodel{1}.postmodel.postplot.partstatic='off'; gui.flmodel{1}.postmodel.postplot.partres='5'; gui.flmodel{1}.postmodel.postplot.maxminplot='off'; gui.flmodel{1}.postmodel.postplot.maxminsubdata={'mises_smsld'}; gui.flmodel{1}.postmodel.postplot.maxminsubdatacheck='on'; gui.flmodel{1}.postmodel.postplot.maxminsubunit='Pa'; gui.flmodel{1}.postmodel.postplot.maxminbnddata={'disp_smsld'}; gui.flmodel{1}.postmodel.postplot.maxminbnddatacheck='off'; gui.flmodel{1}.postmodel.postplot.maxminbndunit='m'; gui.flmodel{1}.postmodel.postplot.maxminedgdata={'disp_smsld'}; gui.flmodel{1}.postmodel.postplot.maxminedgdatacheck='off'; gui.flmodel{1}.postmodel.postplot.maxminedgunit='m'; gui.flmodel{1}.postmodel.postplot.geom='on'; gui.flmodel{1}.postmodel.postplot.roughplot='off'; gui.flmodel{1}.postmodel.postplot.autorefine='on'; gui.flmodel{1}.postmodel.postplot.refine='5'; gui.flmodel{1}.postmodel.postplot.geomnum={'Geom1'}; gui.flmodel{1}.postmodel.postplot.phase='0'; gui.flmodel{1}.postmodel.postplot.frame='ale'; gui.flmodel{1}.postmodel.postplot.solnum='0'; gui.flmodel{1}.postmodel.postplot.selectvia='stored'; gui.flmodel{1}.postmodel.postplot.autotitle='on'; gui.flmodel{1}.postmodel.postplot.customtitle=''; gui.flmodel{1}.postmodel.postplot.smoothinternal='on'; gui.flmodel{1}.postmodel.postplot.useellogic='off'; gui.flmodel{1}.postmodel.postplot.ellogic=''; gui.flmodel{1}.postmodel.postplot.ellogictype='all'; gui.flmodel{1}.postmodel.postplot.deformplot='on'; gui.flmodel{1}.postmodel.postplot.deformsub='on'; gui.flmodel{1}.postmodel.postplot.deformbnd='on'; gui.flmodel{1}.postmodel.postplot.deformedg='on'; gui.flmodel{1}.postmodel.postplot.deformsubdata={'u','v','w'}; gui.flmodel{1}.postmodel.postplot.deformsubunit='m'; gui.flmodel{1}.postmodel.postplot.deformbnddata={'u','v','w'}; gui.flmodel{1}.postmodel.postplot.deformbndunit='m'; gui.flmodel{1}.postmodel.postplot.deformedgdata={'u','v','w'}; gui.flmodel{1}.postmodel.postplot.deformedgunit='m'; gui.flmodel{1}.postmodel.postplot.deformautoscale='on'; gui.flmodel{1}.postmodel.postplot.deformscale='2849.56665'; gui.flmodel{1}.postmodel.postplot.animate_solnum='0'; gui.flmodel{1}.postmodel.postplot.animate_selectvia='stored'; gui.flmodel{1}.postmodel.postplot.filetype='AVI'; gui.flmodel{1}.postmodel.postplot.width='640'; gui.flmodel{1}.postmodel.postplot.height='480'; gui.flmodel{1}.postmodel.postplot.fps='10'; gui.flmodel{1}.postmodel.postplot.statfunctype='full'; gui.flmodel{1}.postmodel.postplot.statnframes='11'; gui.flmodel{1}.postmodel.postplot.reverse='off'; gui.flmodel{1}.postmodel.postplot.movieinmatlab='off'; gui.flmodel{1}.postmodel.postplot.copyaxis='off'; gui.flmodel{1}.postmodel.intdata{1}.intdata={'disp_smsld'}; gui.flmodel{1}.postmodel.intdata{1}.intunit='m'; gui.flmodel{1}.postmodel.intdata{1}.phase='0'; gui.flmodel{1}.postmodel.intdata{1}.solnum='0'; gui.flmodel{1}.postmodel.intdata{1}.selectvia='stored'; gui.flmodel{1}.postmodel.intdata{2}.autoorder='on'; gui.flmodel{1}.postmodel.intdata{2}.multiplyexpr='off'; gui.flmodel{1}.postmodel.intdata{2}.intdata={'disp_smsld'}; gui.flmodel{1}.postmodel.intdata{2}.intunit='m^2'; gui.flmodel{1}.postmodel.intdata{2}.phase='0'; gui.flmodel{1}.postmodel.intdata{2}.frame='ale'; gui.flmodel{1}.postmodel.intdata{2}.solnum='0'; gui.flmodel{1}.postmodel.intdata{2}.selectvia='stored'; gui.flmodel{1}.postmodel.intdata{3}.autoorder='on'; gui.flmodel{1}.postmodel.intdata{3}.multiplyexpr='off'; gui.flmodel{1}.postmodel.intdata{3}.intdata={'disp_smsld'}; gui.flmodel{1}.postmodel.intdata{3}.intunit='m^3'; gui.flmodel{1}.postmodel.intdata{3}.phase='0'; gui.flmodel{1}.postmodel.intdata{3}.frame='ale'; gui.flmodel{1}.postmodel.intdata{3}.solnum='0'; gui.flmodel{1}.postmodel.intdata{3}.selectvia='stored'; gui.flmodel{1}.postmodel.intdata{4}.autoorder='on'; gui.flmodel{1}.postmodel.intdata{4}.multiplyexpr='off'; gui.flmodel{1}.postmodel.intdata{4}.intdata={'mises_smsld'}; gui.flmodel{1}.postmodel.intdata{4}.intunit='J'; gui.flmodel{1}.postmodel.intdata{4}.phase='0'; gui.flmodel{1}.postmodel.intdata{4}.frame='ale'; gui.flmodel{1}.postmodel.intdata{4}.solnum='0'; gui.flmodel{1}.postmodel.intdata{4}.selectvia='stored'; gui.flmodel{1}.postmodel.domainplot.colordata={'disp_smsld'}; gui.flmodel{1}.postmodel.domainplot.colorrangeauto='on'; gui.flmodel{1}.postmodel.domainplot.colorrangemin=''; gui.flmodel{1}.postmodel.domainplot.colorrangemax=''; gui.flmodel{1}.postmodel.domainplot.colorcont='on'; gui.flmodel{1}.postmodel.domainplot.colorunit='m'; gui.flmodel{1}.postmodel.domainplot.surfacesurfacemap='jet'; gui.flmodel{1}.postmodel.domainplot.surfacesurfacemapdepth='1024'; gui.flmodel{1}.postmodel.domainplot.surfacesurfacebar='on'; gui.flmodel{1}.postmodel.domainplot.surfacesurfaceusemap='on'; gui.flmodel{1}.postmodel.domainplot.surfacesurfacecolor='255,0,0'; gui.flmodel{1}.postmodel.domainplot.surfacesurfacecoloring='interp'; gui.flmodel{1}.postmodel.domainplot.surfacesurfacefill='fill'; gui.flmodel{1}.postmodel.domainplot.surfacexyaxisxaxistype='arc'; gui.flmodel{1}.postmodel.domainplot.surfacexyaxisuseexpr='off'; gui.flmodel{1}.postmodel.domainplot.xaxisdata={'disp_smsld'}; gui.flmodel{1}.postmodel.domainplot.xaxisunit='m'; gui.flmodel{1}.postmodel.domainplot.yaxisdata={'disp_smsld'}; gui.flmodel{1}.postmodel.domainplot.yaxisunit='m'; gui.flmodel{1}.postmodel.domainplot.extrusion='off'; gui.flmodel{1}.postmodel.domainplot.lineyaxisdata={'disp_smsld'}; gui.flmodel{1}.postmodel.domainplot.lineyaxiscont='on'; gui.flmodel{1}.postmodel.domainplot.lineyaxisunit='m'; gui.flmodel{1}.postmodel.domainplot.linexaxisxaxistype='arc'; gui.flmodel{1}.postmodel.domainplot.linexaxisuseexpr='off'; gui.flmodel{1}.postmodel.domainplot.linexaxisdata={'disp_smsld'}; gui.flmodel{1}.postmodel.domainplot.linexaxisunit='m'; gui.flmodel{1}.postmodel.domainplot.linelinestyle='solid'; gui.flmodel{1}.postmodel.domainplot.linelinecolor='cyclecolor'; gui.flmodel{1}.postmodel.domainplot.linelinemarker='none'; gui.flmodel{1}.postmodel.domainplot.linelegend='off'; gui.flmodel{1}.postmodel.domainplot.linelinelabels='off'; gui.flmodel{1}.postmodel.domainplot.linecolor='255,0,0'; gui.flmodel{1}.postmodel.domainplot.linesurfacemap='jet'; gui.flmodel{1}.postmodel.domainplot.linesurfacemapdepth='1024'; gui.flmodel{1}.postmodel.domainplot.linesurfacebar='on'; gui.flmodel{1}.postmodel.domainplot.linesurfaceusemap='on'; gui.flmodel{1}.postmodel.domainplot.linesurfacecolor='255,0,0'; gui.flmodel{1}.postmodel.domainplot.linesurfacecoloring='interp'; gui.flmodel{1}.postmodel.domainplot.linesurfacefill='fill'; gui.flmodel{1}.postmodel.domainplot.pointyaxisdata={'disp_smsld'}; gui.flmodel{1}.postmodel.domainplot.pointyaxisunit='m'; gui.flmodel{1}.postmodel.domainplot.pointxxaxistype=''; gui.flmodel{1}.postmodel.domainplot.pointxuseexpr='off'; gui.flmodel{1}.postmodel.domainplot.pointxdata={'disp_smsld'}; gui.flmodel{1}.postmodel.domainplot.pointxunit='m'; gui.flmodel{1}.postmodel.domainplot.pointlinestyle='solid'; gui.flmodel{1}.postmodel.domainplot.pointlinecolor='cyclecolor'; gui.flmodel{1}.postmodel.domainplot.pointlinemarker='none'; gui.flmodel{1}.postmodel.domainplot.pointlegend='off'; gui.flmodel{1}.postmodel.domainplot.pointlinelabels='off'; gui.flmodel{1}.postmodel.domainplot.pointcolor='255,0,0'; gui.flmodel{1}.postmodel.domainplot.crossdispcolor='255,0,0'; gui.flmodel{1}.postmodel.domainplot.phase='0'; gui.flmodel{1}.postmodel.domainplot.frame='ale'; gui.flmodel{1}.postmodel.domainplot.solnum='0'; gui.flmodel{1}.postmodel.domainplot.selectvia='stored'; gui.flmodel{1}.postmodel.domainplot.autotitle='on'; gui.flmodel{1}.postmodel.domainplot.customtitle=''; gui.flmodel{1}.postmodel.domainplot.autolabelx='on'; gui.flmodel{1}.postmodel.domainplot.customlabelx=''; gui.flmodel{1}.postmodel.domainplot.autolabely='on'; gui.flmodel{1}.postmodel.domainplot.customlabely=''; gui.flmodel{1}.postmodel.domainplot.axistype={'lin','lin'}; gui.flmodel{1}.postmodel.domainplot.smoothinternal='on'; gui.flmodel{1}.postmodel.domainplot.autorefine='on'; gui.flmodel{1}.postmodel.domainplot.refine='1'; gui.flmodel{1}.postmodel.domainplot.plottypeind='0'; gui.flmodel{1}.postmodel.crossplot.refine='1'; gui.flmodel{1}.postmodel.crossplot.autorefine='on'; gui.flmodel{1}.postmodel.crossplot.slicedata={'mises_smsld'}; gui.flmodel{1}.postmodel.crossplot.slicerangeauto='on'; gui.flmodel{1}.postmodel.crossplot.slicerangemin=''; gui.flmodel{1}.postmodel.crossplot.slicerangemax=''; gui.flmodel{1}.postmodel.crossplot.slicecont='on'; gui.flmodel{1}.postmodel.crossplot.sliceunit='Pa'; gui.flmodel{1}.postmodel.crossplot.slicecoord={'0','1','0','0','0','1','0','0','0'}; gui.flmodel{1}.postmodel.crossplot.slicegeombound='on'; gui.flmodel{1}.postmodel.crossplot.slicecoordsys='global'; gui.flmodel{1}.postmodel.crossplot.slicelevels='5'; gui.flmodel{1}.postmodel.crossplot.slicevectorlevels=''; gui.flmodel{1}.postmodel.crossplot.sliceisvector='off'; gui.flmodel{1}.postmodel.crossplot.sliceactive='off'; gui.flmodel{1}.postmodel.crossplot.slicesurfacemap='jet'; gui.flmodel{1}.postmodel.crossplot.slicesurfacemapdepth='1024'; gui.flmodel{1}.postmodel.crossplot.slicesurfacebar='on'; gui.flmodel{1}.postmodel.crossplot.slicesurfaceusemap='on'; gui.flmodel{1}.postmodel.crossplot.slicesurfacecolor='255,0,0'; gui.flmodel{1}.postmodel.crossplot.slicesurfacecoloring='interp'; gui.flmodel{1}.postmodel.crossplot.slicesurfacefill='fill'; gui.flmodel{1}.postmodel.crossplot.extrusion='off'; gui.flmodel{1}.postmodel.crossplot.lineyaxisdata={'mises_smsld'}; gui.flmodel{1}.postmodel.crossplot.lineyaxisunit='Pa'; gui.flmodel{1}.postmodel.crossplot.linexaxisxaxistype='arc'; gui.flmodel{1}.postmodel.crossplot.linexaxisuseexpr='off'; gui.flmodel{1}.postmodel.crossplot.linexaxisdata={'mises_smsld'}; gui.flmodel{1}.postmodel.crossplot.linexaxisunit='Pa'; gui.flmodel{1}.postmodel.crossplot.linelinestyle='solid'; gui.flmodel{1}.postmodel.crossplot.linelinecolor='cyclecolor'; gui.flmodel{1}.postmodel.crossplot.linelinemarker='none'; gui.flmodel{1}.postmodel.crossplot.linelegend='off'; gui.flmodel{1}.postmodel.crossplot.linelinelabels='off'; gui.flmodel{1}.postmodel.crossplot.linecolor='255,0,0'; gui.flmodel{1}.postmodel.crossplot.linesurfacemap='jet'; gui.flmodel{1}.postmodel.crossplot.linesurfacemapdepth='1024'; gui.flmodel{1}.postmodel.crossplot.linesurfacebar='on'; gui.flmodel{1}.postmodel.crossplot.linesurfaceusemap='on'; gui.flmodel{1}.postmodel.crossplot.linesurfacecolor='255,0,0'; gui.flmodel{1}.postmodel.crossplot.linesurfacecoloring='interp'; gui.flmodel{1}.postmodel.crossplot.linesurfacefill='fill'; gui.flmodel{1}.postmodel.crossplot.lineresolution='200'; gui.flmodel{1}.postmodel.crossplot.linecoord={'0','1','0','0','0','0'}; gui.flmodel{1}.postmodel.crossplot.pointyaxisdata={'mises_smsld'}; gui.flmodel{1}.postmodel.crossplot.pointyaxisunit='Pa'; gui.flmodel{1}.postmodel.crossplot.pointxxaxistype=''; gui.flmodel{1}.postmodel.crossplot.pointxuseexpr='off'; gui.flmodel{1}.postmodel.crossplot.pointxdata={'mises_smsld'}; gui.flmodel{1}.postmodel.crossplot.pointxunit='Pa'; gui.flmodel{1}.postmodel.crossplot.pointlinestyle='solid'; gui.flmodel{1}.postmodel.crossplot.pointlinecolor='cyclecolor'; gui.flmodel{1}.postmodel.crossplot.pointlinemarker='none'; gui.flmodel{1}.postmodel.crossplot.pointlegend='off'; gui.flmodel{1}.postmodel.crossplot.pointlinelabels='off'; gui.flmodel{1}.postmodel.crossplot.pointcolor='255,0,0'; gui.flmodel{1}.postmodel.crossplot.pointcoord={'0','0','0'}; gui.flmodel{1}.postmodel.crossplot.crossdispcolor='255,0,0'; gui.flmodel{1}.postmodel.crossplot.phase='0'; gui.flmodel{1}.postmodel.crossplot.frame='ale'; gui.flmodel{1}.postmodel.crossplot.solnum='0'; gui.flmodel{1}.postmodel.crossplot.selectvia='stored'; gui.flmodel{1}.postmodel.crossplot.autotitle='on'; gui.flmodel{1}.postmodel.crossplot.customtitle=''; gui.flmodel{1}.postmodel.crossplot.autolabelx='on'; gui.flmodel{1}.postmodel.crossplot.customlabelx=''; gui.flmodel{1}.postmodel.crossplot.autolabely='on'; gui.flmodel{1}.postmodel.crossplot.customlabely=''; gui.flmodel{1}.postmodel.crossplot.axistype={'lin','lin'}; gui.flmodel{1}.postmodel.crossplot.smoothinternal='on'; gui.flmodel{1}.postmodel.crossplot.plottypeind='1'; gui.flmodel{1}.postmodel.dataexport.pntdata={'disp_smsld'}; gui.flmodel{1}.postmodel.dataexport.pntunit='m'; gui.flmodel{1}.postmodel.dataexport.pntlocation='element'; gui.flmodel{1}.postmodel.dataexport.pntlagorder='2'; gui.flmodel{1}.postmodel.dataexport.edgdata={'disp_smsld'}; gui.flmodel{1}.postmodel.dataexport.edgcont='off'; gui.flmodel{1}.postmodel.dataexport.edgunit='m'; gui.flmodel{1}.postmodel.dataexport.edglocation='element'; gui.flmodel{1}.postmodel.dataexport.edglagorder='2'; gui.flmodel{1}.postmodel.dataexport.bnddata={'disp_smsld'}; gui.flmodel{1}.postmodel.dataexport.bndcont='off'; gui.flmodel{1}.postmodel.dataexport.bndunit='m'; gui.flmodel{1}.postmodel.dataexport.bndlocation='element'; gui.flmodel{1}.postmodel.dataexport.bndlagorder='2'; gui.flmodel{1}.postmodel.dataexport.subdata={'mises_smsld'}; gui.flmodel{1}.postmodel.dataexport.subcont='off'; gui.flmodel{1}.postmodel.dataexport.subunit='Pa'; gui.flmodel{1}.postmodel.dataexport.subxspacing='10'; gui.flmodel{1}.postmodel.dataexport.subxvectorspacing=''; gui.flmodel{1}.postmodel.dataexport.subxisvector='off'; gui.flmodel{1}.postmodel.dataexport.subyspacing='10'; gui.flmodel{1}.postmodel.dataexport.subyvectorspacing=''; gui.flmodel{1}.postmodel.dataexport.subyisvector='off'; gui.flmodel{1}.postmodel.dataexport.subzspacing='10'; gui.flmodel{1}.postmodel.dataexport.subzvectorspacing=''; gui.flmodel{1}.postmodel.dataexport.subzisvector='off'; gui.flmodel{1}.postmodel.dataexport.sublocation='element'; gui.flmodel{1}.postmodel.dataexport.sublagorder='2'; gui.flmodel{1}.postmodel.dataexport.phase='0'; gui.flmodel{1}.postmodel.dataexport.frame='ale'; gui.flmodel{1}.postmodel.dataexport.solnum='0'; gui.flmodel{1}.postmodel.dataexport.selectvia='stored'; gui.flmodel{1}.postmodel.dataexport.exportformat='ptd'; gui.flmodel{1}.postmodel.dataexport.exportedim='3'; gui.flmodel{1}.postmodel.dataexport.struct='off'; gui.flmodel{1}.postmodel.datadisplay.fullprecision='off'; gui.flmodel{1}.postmodel.datadisplay.phase='0'; gui.flmodel{1}.postmodel.datadisplay.frame='ale'; gui.flmodel{1}.postmodel.datadisplay.solnum='0'; gui.flmodel{1}.postmodel.datadisplay.selectvia='stored'; gui.flmodel{1}.postmodel.datadisplay.interpdata={'mises_smsld'}; gui.flmodel{1}.postmodel.datadisplay.interpunit='Pa'; gui.flmodel{1}.postmodel.datadisplay.coord={'0','0','0'}; gui.flmodel{1}.postmodel.globalplot.globalyaxisexprs={}; gui.flmodel{1}.postmodel.globalplot.globalyaxisexprsdisp={}; gui.flmodel{1}.postmodel.globalplot.linexaxisxaxistype=''; gui.flmodel{1}.postmodel.globalplot.linexaxisuseexpr='off'; gui.flmodel{1}.postmodel.globalplot.globallinestyle='solid'; gui.flmodel{1}.postmodel.globalplot.globallinecolor='cyclecolor'; gui.flmodel{1}.postmodel.globalplot.globallinemarker='none'; gui.flmodel{1}.postmodel.globalplot.globallegend='off'; gui.flmodel{1}.postmodel.globalplot.globallinelabels='off'; gui.flmodel{1}.postmodel.globalplot.globalcolor='255,0,0'; gui.flmodel{1}.postmodel.globalplot.autotitle='on'; gui.flmodel{1}.postmodel.globalplot.customtitle=''; gui.flmodel{1}.postmodel.globalplot.autolabelx='on'; gui.flmodel{1}.postmodel.globalplot.customlabelx=''; gui.flmodel{1}.postmodel.globalplot.autolabely='on'; gui.flmodel{1}.postmodel.globalplot.customlabely=''; gui.flmodel{1}.postmodel.globalplot.axistype={'lin','lin'}; gui.flmodel{1}.postmodel.globalplot.phase='0'; gui.flmodel{1}.postmodel.globalplot.solnum='0'; gui.flmodel{1}.postmodel.globalplot.selectvia='stored'; gui.flmodel{1}.postmodel.globalplot.plotimported='off'; gui.flmodel{1}.geommodel.suppressed{1}=[]; gui.flmodel{1}.geommodel.suppressed{2}=[]; gui.flmodel{1}.geommodel.suppressed{3}=[]; gui.flmodel{1}.geommodel.suppressed{4}=[]; gui.flmodel{1}.meshmodel.meshplot.subplot='off'; gui.flmodel{1}.meshmodel.meshplot.meshsubuseexpression='off'; gui.flmodel{1}.meshmodel.meshplot.meshsubcolor='192,192,192'; gui.flmodel{1}.meshmodel.meshplot.meshsubbordercheck='on'; gui.flmodel{1}.meshmodel.meshplot.meshwiresubplot='on'; gui.flmodel{1}.meshmodel.meshplot.meshwiresubcolor='0,0,0'; gui.flmodel{1}.meshmodel.meshplot.bndplot='on'; gui.flmodel{1}.meshmodel.meshplot.meshbnduseexpression='off'; gui.flmodel{1}.meshmodel.meshplot.meshbndcolor='128,128,128'; gui.flmodel{1}.meshmodel.meshplot.meshbndbordercheck='on'; gui.flmodel{1}.meshmodel.meshplot.meshwirebndplot='on'; gui.flmodel{1}.meshmodel.meshplot.meshwirebndcolor='0,0,0'; gui.flmodel{1}.meshmodel.meshplot.edgplot='off'; gui.flmodel{1}.meshmodel.meshplot.meshedgcolor='0,0,0'; gui.flmodel{1}.meshmodel.meshplot.useellogic='off'; gui.flmodel{1}.meshmodel.meshplot.ellogic=''; gui.flmodel{1}.meshmodel.meshplot.ellogictype='all'; gui.flmodel{1}.meshmodel.meshplot.meshkeepfraction='1'; gui.flmodel{1}.meshmodel.meshplot.meshkeeptype='random'; gui.flmodel{1}.meshmodel.meshplot.meshmap='jet'; gui.flmodel{1}.meshmodel.meshplot.meshmapdepth='1024'; gui.flmodel{1}.meshmodel.meshplot.meshbar='on'; uq ~ t��clear fem0 clear vrsn vrsn.name = 'COMSOL 3.3'; vrsn.ext = ''; vrsn.major = 0; vrsn.build = 405; vrsn.rcs = '$Name: $'; vrsn.date = '$Date: 2006/08/31 18:03:47 $'; fem0.version = vrsn; fem0.id = 43; fem0.geomdata = 'geom'; fem0.eqvars = 'on'; fem0.cplbndeq = 'on'; fem0.cplbndsh = 'off'; fem0.drawvalid = 'on'; fem0.geomvalid = 'on'; fem0.solvalid = 'on'; fem0.meshtime = 'TIME'; clear appl appl.mode.class = 'SmeSolid3'; appl.mode.type = 'cartesian'; appl.dim = {'u','v','w','p'}; appl.sdim = {'X','Y','Z'}; appl.name = 'smsld'; appl.module = 'SME'; appl.shape = {'shlag(2,''u'')','shlag(2,''v'')','shlag(2,''w'')'}; appl.gporder = 4; appl.cporder = 2; appl.sshape = 2; appl.border = 'off'; appl.assignsuffix = '_smsld'; appl.assign = {'Ex','EX_smsld','Ey','EY_smsld','Ez','EZ_smsld','Fx','FX_smsld', ... 'FxAmp','FXAmp_smsld','FxPh','FXPh_smsld','Fxg','FXg_smsld','Fy','FY_smsld', ... 'FyAmp','FYAmp_smsld','FyPh','FYPh_smsld','Fyg','FYg_smsld','Fz','FZ_smsld', ... 'FzAmp','FZAmp_smsld','FzPh','FZPh_smsld','Fzg','FZg_smsld','Gxy','GXY_smsld', ... 'Gxz','GXZ_smsld','Gyz','GYZ_smsld','Tax','TaX_smsld','Tay','TaY_smsld', ... 'Taz','TaZ_smsld','alphax','alphaX_smsld','alphay','alphaY_smsld', ... 'alphaz','alphaZ_smsld','ex','eX_smsld','exi','eXi_smsld','exy','eXY_smsld', ... 'exyi','eXYi_smsld','exz','eXZ_smsld','exzi','eXZi_smsld','ey','eY_smsld', ... 'eyi','eYi_smsld','eyz','eYZ_smsld','eyzi','eYZi_smsld','ez','eZ_smsld', ... 'ezi','eZi_smsld','nuxy','nuXY_smsld','nuxz','nuXZ_smsld','nuyz','nuYZ_smsld', ... 'sx','sX_smsld','sxi','sXi_smsld','sxy','sXY_smsld','sxyi','sXYi_smsld', ... 'sxz','sXZ_smsld','sxzi','sXZi_smsld','sy','sY_smsld','syi','sYi_smsld', ... 'syz','sYZ_smsld','syzi','sYZi_smsld','sz','sZ_smsld','szi','sZi_smsld'}; clear prop prop.elemdefault='Lag2'; prop.analysis='static'; prop.largedef='on'; prop.impl='weak'; prop.eigtype='freq'; prop.createframe='off'; prop.deformframe='ref'; prop.frame='ref'; clear weakconstr weakconstr.value = 'off'; weakconstr.dim = {'lm1','lm2','lm3'}; prop.weakconstr = weakconstr; appl.prop = prop; clear pnt pnt.Fx = {'0'}; pnt.FxAmp = {'1'}; pnt.FxPh = {'0'}; pnt.Fy = {'0'}; pnt.FyAmp = {'1'}; pnt.FyPh = {'0'}; pnt.Fz = {'0'}; pnt.FzAmp = {'1'}; pnt.FzPh = {'0'}; pnt.loadcoord = {'global'}; pnt.constrcoord = {'global'}; pnt.constrtype = {'standard'}; pnt.H = {{'0','0','0';'0','0','0';'0','0','0'}}; pnt.R = {{'0';'0';'0'}}; pnt.Rx = {'0'}; pnt.Hx = {'0'}; pnt.Ry = {'0'}; pnt.Hy = {'0'}; pnt.Rz = {'0'}; pnt.Hz = {'0'}; pnt.style = {{{'0'},{'0','0','0'}}}; pnt.ind = [1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1]; appl.pnt = pnt; clear edg edg.Fx = {'0'}; edg.FxAmp = {'1'}; edg.FxPh = {'0'}; edg.Fy = {'0'}; edg.FyAmp = {'1'}; edg.FyPh = {'0'}; edg.Fz = {'0'}; edg.FzAmp = {'1'}; edg.FzPh = {'0'}; edg.loadcoord = {'global'}; edg.constrcoord = {'global'}; edg.constrtype = {'standard'}; edg.H = {{'0','0','0';'0','0','0';'0','0','0'}}; edg.R = {{'0';'0';'0'}}; edg.Rx = {'0'}; edg.Hx = {'0'}; edg.Ry = {'0'}; edg.Hy = {'0'}; edg.Rz = {'0'}; edg.Hz = {'0'}; edg.style = {{{'0'},{'0','0','0'},{'solid'}}}; edg.ind = [1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1, ... 1,1,1,1,1,1]; appl.edg = edg; clear bnd bnd.name = {'Fluid load','',''}; bnd.Fx = {'-lm7*dvol_ale/dvol','0','0'}; bnd.FxAmp = {'1','1','1'}; bnd.FxPh = {'0','0','0'}; bnd.Fy = {'-lm8*dvol_ale/dvol','0','0'}; bnd.FyAmp = {'1','1','1'}; bnd.FyPh = {'0','0','0'}; bnd.Fz = {'-lm9*dvol_ale/dvol','0','0'}; bnd.FzAmp = {'1','1','1'}; bnd.FzPh = {'0','0','0'}; bnd.loadcoord = {'global','global','global'}; bnd.constrcoord = {'global','global','global'}; bnd.constrtype = {'standard','standard','standard'}; bnd.constrcond = {'free','free','fixed'}; bnd.velx = {'0','0','0'}; bnd.Hvelx = {'0','0','0'}; bnd.accx = {'0','0','0'}; bnd.Haccx = {'0','0','0'}; bnd.vely = {'0','0','0'}; bnd.Hvely = {'0','0','0'}; bnd.accy = {'0','0','0'}; bnd.Haccy = {'0','0','0'}; bnd.velz = {'0','0','0'}; bnd.Hvelz = {'0','0','0'}; bnd.accz = {'0','0','0'}; bnd.Haccz = {'0','0','0'}; bnd.H = {{'0','0','0';'0','0','0';'0','0','0'},{'0','0','0';'0','0','0'; ... '0','0','0'},{'0','0','0';'0','0','0';'0','0','0'}}; bnd.R = {{'0';'0';'0'},{'0';'0';'0'},{'0';'0';'0'}}; bnd.Rx = {'0','0','0'}; bnd.Hx = {'0','0','0'}; bnd.Ry = {'0','0','0'}; bnd.Hy = {'0','0','0'}; bnd.Rz = {'0','0','0'}; bnd.Hz = {'0','0','0'}; bnd.offset = {'0','0','0'}; bnd.pn = {'E_smsld/h*min(1e-4*5^auglagiter,0.1)','E_smsld/h*min(1e-4*5^auglagiter,0.1)', ... 'E_smsld/h*min(1e-4*5^auglagiter,0.1)'}; bnd.pt = {'E_smsld/h*min(1e-4*5^auglagiter,0.1)','E_smsld/h*min(1e-4*5^auglagiter,0.1)', ... 'E_smsld/h*min(1e-4*5^auglagiter,0.1)'}; bnd.frictiontype = {'nofriction','nofriction','nofriction'}; bnd.searchmethod = {'fast','fast','fast'}; bnd.searchdist = {'auto','auto','auto'}; bnd.mandist = {'1e-2','1e-2','1e-2'}; bnd.contacttol = {'auto','auto','auto'}; bnd.mantol = {'1e-6','1e-6','1e-6'}; bnd.dynfric = {'0','0','0'}; bnd.mustat = {'0','0','0'}; bnd.mudyn = {'0','0','0'}; bnd.dcfric = {'0','0','0'}; bnd.cohe = {'0','0','0'}; bnd.Ttmax = {'Inf','Inf','Inf'}; bnd.Tni = {'1e6','1e6','1e6'}; bnd.contact_oldi = {'0','0','0'}; bnd.Ttxi = {'1e6','1e6','1e6'}; bnd.xm_oldi = {'0','0','0'}; bnd.Ttyi = {'1e6','1e6','1e6'}; bnd.ym_oldi = {'0','0','0'}; bnd.Ttzi = {'1e6','1e6','1e6'}; bnd.zm_oldi = {'0','0','0'}; bnd.style = {{{'0'},{'217','217','217'}},{{'0'},{'0','0','255'}},{{'0'},{'0', ... '255','0'}}}; bnd.ind = [2,2,2,2,2,2,2,2,2,2,2,1,1,3,3,1,1,2]; appl.bnd = bnd; clear equ equ.shape = {[1;2;3],[1;2;3]}; equ.gporder = {{1;1;1;1},{1;1;1;1}}; equ.cporder = {{1;1;1;1},{1;1;1;1}}; equ.init = {{'0';'0';'0';'0'},{'0';'0';'0';'0'}}; equ.usage = {0,1}; equ.name = {'Fluid domain','Solid domain'}; equ.E = {'2.0e11','31100'}; equ.rho = {'7850','7850'}; equ.dampingtype = {'Rayleigh','Rayleigh'}; equ.alphadM = {'1','1'}; equ.betadK = {'0.001','0.001'}; equ.eta = {'0','0'}; equ.dinit = {{'0';'0';'0';'0'},{'0';'0';'0';'0'}}; equ.materialcoord = {'global','global'}; equ.materialmodel = {'iso','iso'}; equ.mixedform = {'0','0'}; equ.hypertype = {'neo_hookean','neo_hookean'}; equ.mu = {'8e5','8e5'}; equ.C10 = {'2e5','2e5'}; equ.C01 = {'2e5','2e5'}; equ.kappa = {'1e10','1e10'}; equ.hardeningmodel = {'iso','iso'}; equ.yieldtype = {'mises','mises'}; equ.isodata = {'tangent','tangent'}; equ.ETkin = {'2.0e10','2.0e10'}; equ.ETiso = {'2.0e10','2.0e10'}; equ.Sys = {'2.0e8','2.0e8'}; equ.Syfunc = {'mises_smsld','mises_smsld'}; equ.Syfunc_kin = {'misesKin_smsld','misesKin_smsld'}; equ.Shard = {'2.0e10/(1-2.0e10/2.0e11)*epe_smsld','2.0e10/(1-2.0e10/2.0e11)*epe_smsld'}; equ.ini_stress = {'0','0'}; equ.ini_strain = {'0','0'}; equ.alpha = {'1.2e-5','1.2e-5'}; equ.nu = {'0.33','.49'}; equ.Ex = {'2.0e11','2.0e11'}; equ.nuxy = {'0.33','0.33'}; equ.Gxy = {'7.52e10','7.52e10'}; equ.alphax = {'1.2e-5','1.2e-5'}; equ.sxi = {'0','0'}; equ.sxyi = {'0','0'}; equ.exi = {'0','0'}; equ.exyi = {'0','0'}; equ.Ey = {'2.0e11','2.0e11'}; equ.nuyz = {'0.33','0.33'}; equ.Gyz = {'7.52e10','7.52e10'}; equ.alphay = {'1.2e-5','1.2e-5'}; equ.syi = {'0','0'}; equ.syzi = {'0','0'}; equ.eyi = {'0','0'}; equ.eyzi = {'0','0'}; equ.Ez = {'2.0e11','2.0e11'}; equ.nuxz = {'0.33','0.33'}; equ.Gxz = {'7.52e10','7.52e10'}; equ.alphaz = {'1.2e-5','1.2e-5'}; equ.szi = {'0','0'}; equ.sxzi = {'0','0'}; equ.ezi = {'0','0'}; equ.exzi = {'0','0'}; equ.D = {{'2.0e11/((1+0.33)*(1-2*0.33))*(1-0.33)','2.0e11/((1+0.33)*(1-2*0.33))*0.33', ... '2.0e11/((1+0.33)*(1-2*0.33))*0.33','0','0','0';'2.0e11/((1+0.33)*(1-2*0.33))*0.33', ... '2.0e11/((1+0.33)*(1-2*0.33))*(1-0.33)','2.0e11/((1+0.33)*(1-2*0.33))*0.33', ... '0','0','0';'2.0e11/((1+0.33)*(1-2*0.33))*0.33','2.0e11/((1+0.33)*(1-2*0.33))*0.33', ... '2.0e11/((1+0.33)*(1-2*0.33))*(1-0.33)','0','0','0';'0','0','0','2.0e11/((1+0.33)*2)', ... '0','0';'0','0','0','0','2.0e11/((1+0.33)*2)','0';'0','0','0','0','0','2.0e11/((1+0.33)*2)'}, ... {'2.0e11/((1+0.33)*(1-2*0.33))*(1-0.33)','2.0e11/((1+0.33)*(1-2*0.33))*0.33', ... '2.0e11/((1+0.33)*(1-2*0.33))*0.33','0','0','0';'2.0e11/((1+0.33)*(1-2*0.33))*0.33', ... '2.0e11/((1+0.33)*(1-2*0.33))*(1-0.33)','2.0e11/((1+0.33)*(1-2*0.33))*0.33', ... '0','0','0';'2.0e11/((1+0.33)*(1-2*0.33))*0.33','2.0e11/((1+0.33)*(1-2*0.33))*0.33', ... '2.0e11/((1+0.33)*(1-2*0.33))*(1-0.33)','0','0','0';'0','0','0','2.0e11/((1+0.33)*2)', ... '0','0';'0','0','0','0','2.0e11/((1+0.33)*2)','0';'0','0','0','0','0','2.0e11/((1+0.33)*2)'}}; equ.alphavector = {{'1.2e-5';'1.2e-5';'1.2e-5';'0';'0';'0'},{'1.2e-5';'1.2e-5';'1.2e-5'; ... '0';'0';'0'}}; equ.Fx = {'0','0'}; equ.FxAmp = {'1','1'}; equ.FxPh = {'0','0'}; equ.Fy = {'0','0'}; equ.FyAmp = {'1','1'}; equ.FyPh = {'0','0'}; equ.Fz = {'0','0'}; equ.FzAmp = {'1','1'}; equ.FzPh = {'0','0'}; equ.loadcoord = {'global','global'}; equ.Tflag = {'0','0'}; equ.Temp = {'0','0'}; equ.Tempref = {'0','0'}; equ.constrcoord = {'global','global'}; equ.constrtype = {'standard','standard'}; equ.constrcond = {'free','free'}; equ.velx = {'0','0'}; equ.Hvelx = {'0','0'}; equ.accx = {'0','0'}; equ.Haccx = {'0','0'}; equ.vely = {'0','0'}; equ.Hvely = {'0','0'}; equ.accy = {'0','0'}; equ.Haccy = {'0','0'}; equ.velz = {'0','0'}; equ.Hvelz = {'0','0'}; equ.accz = {'0','0'}; equ.Haccz = {'0','0'}; equ.H = {{'0','0','0';'0','0','0';'0','0','0'},{'0','0','0';'0','0','0'; ... '0','0','0'}}; equ.R = {{'0';'0';'0'},{'0';'0';'0'}}; equ.Rx = {'0','0'}; equ.Hx = {'0','0'}; equ.Ry = {'0','0'}; equ.Hy = {'0','0'}; equ.Rz = {'0','0'}; equ.Hz = {'0','0'}; equ.dx = {'1','1'}; equ.Lx = {'dx_smsld','dx_smsld'}; equ.dy = {'1','1'}; equ.Ly = {'dy_smsld','dy_smsld'}; equ.dz = {'1','1'}; equ.Lz = {'dz_smsld','dz_smsld'}; equ.dr = {'1','1'}; equ.Lr = {'dr_smsld','dr_smsld'}; equ.R0 = {'1','1'}; equ.userPML = {{'0';'0';'0'},{'0';'0';'0'}}; equ.PMLtype = {'none','none'}; equ.coordOnx = {'0','0'}; equ.coordOny = {'0','0'}; equ.coordOnz = {'0','0'}; equ.x0 = {'0','0'}; equ.y0 = {'0','0'}; equ.z0 = {'0','0'}; equ.srcaxis = {{'0';'0';'0'},{'0';'0';'0'}}; equ.coord = {'global','global'}; equ.rOn = {'0','0'}; equ.style = {{{'0'},{'217','217','217'}},{{'0'},{'0','0','255'}}}; equ.ind = [1,1,2]; appl.equ = equ; appl.var = {'t_old_ini','-1'}; fem0.appl{1} = appl; clear appl appl.mode.class = 'MovingMesh'; appl.mode.type = 'cartesian'; appl.dim = {}; appl.sdim = {'X','Y','Z'}; appl.name = 'ale'; appl.shape = {'shlag(2,''lm4'')','shlag(2,''lm5'')','shlag(2,''lm6'')','shlag(2,''x'')','shlag(2,''y'')','shlag(2,''z'')'}; appl.gporder = {10,4}; appl.cporder = 2; appl.sshape = 2; appl.border = 'off'; appl.assignsuffix = '_ale'; clear prop prop.elemdefault='Lag2'; prop.smoothing='winslow'; prop.analysis='static'; prop.allowremesh='off'; prop.definesframe='ale'; prop.origrefframe='ale'; prop.frame='ref'; clear weakconstr weakconstr.value = 'non-ideal'; weakconstr.dim = {'lm4','lm5','lm6'}; prop.weakconstr = weakconstr; appl.prop = prop; clear bnd bnd.name = {'Fixed','','Structural displacement'}; bnd.constrcoord = {'global','global','global'}; bnd.type = {'def','def','def'}; bnd.deform = {{'0';'0';'0'},{'0';'0';'0'},{'u';'v';'w'}}; bnd.veldeform = {{'0';'0';'0'},{'0';'0';'0'},{'0';'0';'0'}}; bnd.defflag = {{'1';'1';'1'},{'0';'0';'0'},{'1';'1';'1'}}; bnd.veldefflag = {{'0';'0';'0'},{'0';'0';'0'},{'0';'0';'0'}}; bnd.weakconstr = {'1','1','1'}; bnd.wcshape = {[1;2;3],[1;2;3],[1;2;3]}; bnd.wcgporder = {{1;1;1},{1;1;1},{1;1;1}}; bnd.wcinit = {{'0';'0';'0'},{'0';'0';'0'},{'0';'0';'0'}}; bnd.style = {{{'0'},{'217','217','217'}},{{'0'},{'0','255','0'}},{{'0'},{'217', ... '217','217'}}}; bnd.ind = [1,1,1,1,1,1,1,1,1,2,1,3,3,2,2,3,3,1]; appl.bnd = bnd; clear equ equ.shape = {[4;5;6],[4;5;6]}; equ.gporder = {{2;2;2},{2;2;2}}; equ.cporder = {{1;1;1},{1;1;1}}; equ.init = {{'xinit_ale';'yinit_ale';'zinit_ale'},{'xinit_ale';'yinit_ale'; ... 'zinit_ale'}}; equ.usage = {1,1}; equ.name = {'Fluid domain','Solid domain'}; equ.type = {'free','phys'}; equ.physexpr = {{'0';'0';'0'},{'u';'v';'w'}}; equ.presexpr = {{'0';'0';'0'},{'0';'0';'0'}}; equ.dinit = {{'0';'0';'0'},{'0';'0';'0'}}; equ.style = {{{'0'},{'217','217','217'}},{{'0'},{'217','217','217'}}}; equ.ind = [1,1,2]; appl.equ = equ; fem0.appl{2} = appl; clear appl appl.mode.class = 'NavierStokes'; appl.mode.type = 'cartesian'; appl.dim = {'u2','v2','w2','p2','logk','logd','nxw','nyw','nzw'}; appl.sdim = {'x','y','z'}; appl.name = 'chns'; appl.module = 'CHEM'; appl.shape = {'shlag(2,''lm7'')','shlag(2,''lm8'')','shlag(2,''lm9'')','shlag(2,''u2'')','shlag(2,''v2'')','shlag(2,''w2'')','shlag(1,''p2'')'}; appl.gporder = {10,4,2}; appl.cporder = {2,1}; appl.sshape = 2; appl.border = 'off'; appl.assignsuffix = '_chns'; clear prop prop.elemdefault='Lagp2p1'; prop.analysis='static'; prop.stensor='full'; prop.cornersmoothing='Off'; prop.nisot='Off'; prop.turbmodel='None'; prop.nonnewtonian='Off'; prop.brinkmandef='Off'; prop.GLSdef='false'; prop.frame='ale'; clear weakconstr weakconstr.value = 'non-ideal'; weakconstr.dim = {'lm7','lm8','lm9','lm10','lm11'}; prop.weakconstr = weakconstr; appl.prop = prop; clear pnt pnt.p0 = {'0'}; pnt.pnton = {'0'}; pnt.style = {{{'0'},{'0','0','0'}}}; pnt.ind = [1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1]; appl.pnt = pnt; clear edg edg.style = {{{'0'},{'0','0','0'},{'solid'}}}; edg.ind = [1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1, ... 1,1,1,1,1,1]; appl.edg = edg; clear bnd bnd.name = {'','','','Structural velocity'}; bnd.type = {'uv','neutral','noslip','uv'}; bnd.u0 = {'.5','0','0','ut'}; bnd.v0 = {'0','0','0','vt'}; bnd.w0 = {'0','0','0','wt'}; bnd.p0 = {'0','0','0','0'}; bnd.k0 = {'0','0','0','0'}; bnd.d0 = {'0','0','0','0'}; bnd.TinType = {'LTIT','LTIT','LTIT','LTIT'}; bnd.LT = {'0.07','0.07','0.07','0.07'}; bnd.IT = {'0.054','0.054','0.054','0.054'}; bnd.uw = {'0','0','0','0'}; bnd.vw = {'0','0','0','0'}; bnd.ww = {'0','0','0','0'}; bnd.dwType = {'dw','dw','dw','dw'}; bnd.dw = {'h','h','h','h'}; bnd.yplus = {'100','100','100','100'}; bnd.relExpr = {{},{},{},{}}; bnd.U0 = {'0','0','0','0'}; bnd.V0 = {'0','0','0','0'}; bnd.p0_inl = {'0','0','0','0'}; bnd.L0 = {'1','1','1','1'}; bnd.flowtype = {'velocity','velocity','velocity','velocity'}; bnd.constr0 = {'0','0','0','0'}; bnd.weakconstr = {'1','1','1','1'}; bnd.wcshape = {[1;2;3],[1;2;3],[1;2;3],[1;2;3]}; bnd.wcgporder = {{1;1;1},{1;1;1},{1;1;1},{1;1;1}}; bnd.wcinit = {{'0';'0';'0';'0';'0'},{'0';'0';'0';'0';'0'},{'0';'0';'0';'0'; ... '0'},{'0';'0';'0';'0';'0'}}; bnd.style = {{{'0'},{'0','255','255'}},{{'0'},{'0','255','0'}},{{'0'},{'0', ... '0','255'}},{{'0'},{'217','217','217'}}}; bnd.ind = [1,2,2,2,2,2,3,3,3,2,3,3,3,3,3,3,3,2]; appl.bnd = bnd; clear equ equ.shape = {[4;5;6;7],[4;5;6;7]}; equ.gporder = {{2;2;2;3},{2;2;2;3}}; equ.cporder = {{1;1;1;2},{1;1;1;2}}; equ.init = {{'0';'0';'0';'0';'0';'0';'0';'0';'0'},{'0';'0';'0';'0';'0';'0'; ... '0';'0';'0'}}; equ.usage = {1,0}; equ.name = {'Fluid domain','Solid domain'}; equ.rho = {'1000','1'}; equ.eta = {'.001','1'}; equ.F_x = {'0','0'}; equ.F_y = {'0','0'}; equ.F_z = {'0','0'}; equ.relExpr = {{},{}}; equ.eta0 = {'1','1'}; equ.m = {'0','0'}; equ.n = {'0','0'}; equ.lambda = {'0','0'}; equ.eta_inf = {'0','0'}; equ.type_visc = {'power','power'}; equ.kappadv = {'0','0'}; equ.k = {'1','1'}; equ.brinkmaneqns = {'0','0'}; equ.idon = {{'0'},{'0'}}; equ.delid = {{'0.5'},{'0.5'}}; equ.sdon = {{'0'},{'0'}}; equ.sdtype = {{'pgc'},{'pgc'}}; equ.delsd = {{'0.25'},{'0.25'}}; equ.cdon = {{'0'},{'0'}}; equ.cdtype = {{'sc'},{'sc'}}; equ.delcd = {{'0.35'},{'0.35'}}; equ.pson = {{'0'},{'0'}}; equ.delps = {{'1'},{'1'}}; equ.idTon = {'1','1'}; equ.delidT = {'0.5','0.5'}; equ.sdTon = {'1','1'}; equ.sdTtype = {'ad','ad'}; equ.delsdT = {'0.5','0.5'}; equ.style = {{{'0'},{'217','217','217'}},{{'0'},{'217','217','217'}}}; equ.ind = [1,1,2]; appl.equ = equ; fem0.appl{3} = appl; fem0.geom = flbinary('fem44','geom','tutorial5.mph'); fem0.mesh = flbinary('fem44.0','mesh','tutorial5.mph'); fem0.sdim = {{'X','Y','Z'},{'x','y','z'}}; fem0.frame = {'ref','ale'}; fem0.shape = {'shlag(''order'',2,''basename'',''x'',''frame'',''ale'')','shlag(''order'',2,''basename'',''y'',''frame'',''ale'')','shlag(''order'',2,''basename'',''z'',''frame'',''ale'')','shlag(''order'',2,''basename'',''u2'',''frame'',''ale'')','shlag(''order'',2,''basename'',''v2'',''frame'',''ale'')','shlag(''order'',2,''basename'',''w2'',''frame'',''ale'')','shlag(''order'',1,''basename'',''p2'',''frame'',''ale'')','shlag(2,''u'')','shlag(2,''v'')','shlag(2,''w'')','shlag(2,''lm4'')','shlag(2,''lm5'')','shlag(2,''lm6'')','shlag(2,''lm7'')','shlag(2,''lm8'')','shlag(2,''lm9'')'}; fem0.gporder = {4,2,10}; fem0.cporder = {2,1}; clear sshape sshape.dvolname = 'dvol'; sshape.frame = 'ref'; sshape.sorder = 2; sshape.type = 'fixed'; fem0.sshape{1} = sshape; clear sshape sshape.dvolname = 'dvol_ale'; sshape.frame = 'ale'; sshape.sorder = 2; sshape.type = 'moving_abs'; sshape.refframe = 'ref'; fem0.sshape{2} = sshape; clear sshape sshape.dvolname = 'dvol_ale'; sshape.frame = 'ale'; sshape.sorder = 2; sshape.type = 'moving_rel'; sshape.sdimdofs = {'u','v','w'}; sshape.refframe = 'ref'; fem0.sshape{3} = sshape; fem0.simplify = 'on'; fem0.border = 1; fem0.form = 'weak'; clear units; units.basesystem = 'SI'; fem0.units = units; clear equ equ.shape = {[1;2;3;4;5;6;7],[8;9;10]}; equ.gporder = {{1;1;1;1;1;1;1;1;1;1;2},{1;1;1;1;1;1;1;1;1;1;2}}; equ.cporder = {{1;1;1;1;1;1;1;1;1;1;2},{1;1;1;1;1;1;1;1;1;1;2}}; equ.init = {{'0';'0';'0';'0';'xinit_ale';'yinit_ale';'zinit_ale';'0';'0'; ... '0';'0'},{'0';'0';'0';'0';'0';'0';'0';'0';'0';'0';'0'}}; equ.dinit = {{'0';'0';'0';'0';'0';'0';'0';'0';'0';'0';'0'},{'0';'0';'0';'0'; ... '0';'0';'0';'0';'0';'0';'0'}}; equ.weak = {{'0';'0';'0';'0';'-Xx*test(xx)-Xy*test(xy)-Xz*test(xz)';'-Yx*test(yx)-Yy*test(yy)-Yz*test(yz)'; ... '-Zx*test(zx)-Zy*test(zy)-Zz*test(zz)';'(-2*eta_chns*u2x+p2)*test(u2x)-eta_chns*(u2y+v2x)*test(u2y)-eta_chns*(u2z+w2x)*test(u2z)-rho_chns*(u2*u2x+v2*u2y+w2*u2z)*test(u2)'; ... '-eta_chns*(v2x+u2y)*test(v2x)+(-2*eta_chns*v2y+p2)*test(v2y)-eta_chns*(v2z+w2y)*test(v2z)-rho_chns*(u2*v2x+v2*v2y+w2*v2z)*test(v2)'; ... '-eta_chns*(w2x+u2z)*test(w2x)-eta_chns*(w2y+v2z)*test(w2y)+(-2*eta_chns*w2z+p2)*test(w2z)-rho_chns*(u2*w2x+v2*w2y+w2*w2z)*test(w2)'; ... '-divU_chns*test(p2)'},{'-eX_smsld_test*SX_smsld-eY_smsld_test*SY_smsld-eZ_smsld_test*SZ_smsld-2*eXY_smsld_test*SXY_smsld-2*eYZ_smsld_test*SYZ_smsld-2*eXZ_smsld_test*SXZ_smsld'; ... '0';'0';'0';'0';'0';'0';'0';'0';'0';'0'}}; equ.dweak = {{'0';'0';'0';'0';'0';'0';'0';'rho_chns*u2t*test(u2)';'rho_chns*v2t*test(v2)'; ... 'rho_chns*w2t*test(w2)';'0'},{'0';'0';'0';'0';'0';'0';'0';'0';'0';'0';'0'}}; equ.constr = {{'0';'0';'0';'0';'0';'0';'0';'0';'0';'0';'0'},{'0';'0';'0';'0'; ... '0';'0';'0';'0';'0';'0';'0'}}; equ.sshape = {[1;2],[1;3]}; equ.sshapedim = {{1;1;1;1;2;2;2;2;2;2;2},{1;1;1;1;3;3;3;3;3;3;3}}; equ.ind = [1,1,2]; equ.dim = {'u','v','w','p','x','y','z','u2','v2','w2','p2'}; equ.var = {'FXg_smsld',{'','0'}, ... 'FYg_smsld',{'','0'}, ... 'FZg_smsld',{'','0'}, ... 'disp_smsld',{'','sqrt(real(u)^2+real(v)^2+real(w)^2)'}, ... 'sX_smsld',{'','(F11_smsld*(SX_smsld*F11_smsld+SXY_smsld*F12_smsld+SXZ_smsld*F13_smsld)+F12_smsld*(SXY_smsld*F11_smsld+SY_smsld*F12_smsld+SYZ_smsld*F13_smsld)+F13_smsld*(SXZ_smsld*F11_smsld+SYZ_smsld*F12_smsld+SZ_smsld*F13_smsld))/J_smsld'}, ... 'sY_smsld',{'','(F21_smsld*(SX_smsld*F21_smsld+SXY_smsld*F22_smsld+SXZ_smsld*F23_smsld)+F22_smsld*(SXY_smsld*F21_smsld+SY_smsld*F22_smsld+SYZ_smsld*F23_smsld)+F23_smsld*(SXZ_smsld*F21_smsld+SYZ_smsld*F22_smsld+SZ_smsld*F23_smsld))/J_smsld'}, ... 'sZ_smsld',{'','(F31_smsld*(SX_smsld*F31_smsld+SXY_smsld*F32_smsld+SXZ_smsld*F33_smsld)+F32_smsld*(SXY_smsld*F31_smsld+SY_smsld*F32_smsld+SYZ_smsld*F33_smsld)+F33_smsld*(SXZ_smsld*F31_smsld+SYZ_smsld*F32_smsld+SZ_smsld*F33_smsld))/J_smsld'}, ... 'sXY_smsld',{'','(F11_smsld*(SX_smsld*F21_smsld+SXY_smsld*F22_smsld+SXZ_smsld*F23_smsld)+F12_smsld*(SXY_smsld*F21_smsld+SY_smsld*F22_smsld+SYZ_smsld*F23_smsld)+F13_smsld*(SXZ_smsld*F21_smsld+SYZ_smsld*F22_smsld+SZ_smsld*F23_smsld))/J_smsld'}, ... 'sYZ_smsld',{'','(F21_smsld*(SX_smsld*F31_smsld+SXY_smsld*F32_smsld+SXZ_smsld*F33_smsld)+F22_smsld*(SXY_smsld*F31_smsld+SY_smsld*F32_smsld+SYZ_smsld*F33_smsld)+F23_smsld*(SXZ_smsld*F31_smsld+SYZ_smsld*F32_smsld+SZ_smsld*F33_smsld))/J_smsld'}, ... 'sXZ_smsld',{'','(F11_smsld*(SX_smsld*F31_smsld+SXY_smsld*F32_smsld+SXZ_smsld*F33_smsld)+F12_smsld*(SXY_smsld*F31_smsld+SY_smsld*F32_smsld+SYZ_smsld*F33_smsld)+F13_smsld*(SXZ_smsld*F31_smsld+SYZ_smsld*F32_smsld+SZ_smsld*F33_smsld))/J_smsld'}, ... 'eX_smsld',{'','uX+0.5*(uX^2+vX^2+wX^2)'}, ... 'eY_smsld',{'','vY+0.5*(uY^2+vY^2+wY^2)'}, ... 'eZ_smsld',{'','wZ+0.5*(uZ^2+vZ^2+wZ^2)'}, ... 'eXY_smsld',{'','0.5*(uY+vX+uX*uY+vX*vY+wX*wY)'}, ... 'eYZ_smsld',{'','0.5*(vZ+wY+uY*uZ+vY*vZ+wY*wZ)'}, ... 'eXZ_smsld',{'','0.5*(uZ+wX+uX*uZ+vX*vZ+wX*wZ)'}, ... 'SX_smsld',{'','E_smsld*(1-nu_smsld)*eX_smsld/((1+nu_smsld)*(1-2*nu_smsld))+E_smsld*nu_smsld*eY_smsld/((1+nu_smsld)*(1-2*nu_smsld))+E_smsld*nu_smsld*eZ_smsld/((1+nu_smsld)*(1-2*nu_smsld))'}, ... 'SY_smsld',{'','E_smsld*nu_smsld*eX_smsld/((1+nu_smsld)*(1-2*nu_smsld))+E_smsld*(1-nu_smsld)*eY_smsld/((1+nu_smsld)*(1-2*nu_smsld))+E_smsld*nu_smsld*eZ_smsld/((1+nu_smsld)*(1-2*nu_smsld))'}, ... 'SZ_smsld',{'','E_smsld*nu_smsld*eX_smsld/((1+nu_smsld)*(1-2*nu_smsld))+E_smsld*nu_smsld*eY_smsld/((1+nu_smsld)*(1-2*nu_smsld))+E_smsld*(1-nu_smsld)*eZ_smsld/((1+nu_smsld)*(1-2*nu_smsld))'}, ... 'SXY_smsld',{'','E_smsld*eXY_smsld/(1+nu_smsld)'}, ... 'SYZ_smsld',{'','E_smsld*eYZ_smsld/(1+nu_smsld)'}, ... 'SXZ_smsld',{'','E_smsld*eXZ_smsld/(1+nu_smsld)'}, ... 'K_smsld',{'','E_smsld/(3*(1-2*nu_smsld))'}, ... 'G_smsld',{'','0.5*E_smsld/(1+nu_smsld)'}, ... 'mises_smsld',{'','sqrt(sX_smsld^2+sY_smsld^2+sZ_smsld^2-sX_smsld*sY_smsld-sY_smsld*sZ_smsld-sX_smsld*sZ_smsld+3*sXY_smsld^2+3*sYZ_smsld^2+3*sXZ_smsld^2)'}, ... 'Ws_smsld',{'','0.5*(sX_smsld*eX_smsld+sY_smsld*eY_smsld+sZ_smsld*eZ_smsld+2*sXY_smsld*eXY_smsld+2*sYZ_smsld*eYZ_smsld+2*sXZ_smsld*eXZ_smsld)'}, ... 'evol_smsld',{'','-1+Jel_smsld'}, ... 'F11_smsld',{'','1+uX'}, ... 'F12_smsld',{'','uY'}, ... 'F13_smsld',{'','uZ'}, ... 'F21_smsld',{'','vX'}, ... 'F22_smsld',{'','1+vY'}, ... 'F23_smsld',{'','vZ'}, ... 'F31_smsld',{'','wX'}, ... 'F32_smsld',{'','wY'}, ... 'F33_smsld',{'','1+wZ'}, ... 'detF_smsld',{'','F11_smsld*F22_smsld*F33_smsld+F12_smsld*F23_smsld*F31_smsld+F13_smsld*F21_smsld*F32_smsld-F11_smsld*F23_smsld*F32_smsld-F12_smsld*F21_smsld*F33_smsld-F13_smsld*F22_smsld*F31_smsld'}, ... 'J_smsld',{'','detF_smsld'}, ... 'Jel_smsld',{'','J_smsld'}, ... 'invF11_smsld',{'','(F22_smsld*F33_smsld-F23_smsld*F32_smsld)/detF_smsld'}, ... 'invF12_smsld',{'','(F13_smsld*F32_smsld-F12_smsld*F33_smsld)/detF_smsld'}, ... 'invF13_smsld',{'','(F12_smsld*F23_smsld-F13_smsld*F22_smsld)/detF_smsld'}, ... 'invF21_smsld',{'','(F31_smsld*F23_smsld-F21_smsld*F33_smsld)/detF_smsld'}, ... 'invF22_smsld',{'','(F11_smsld*F33_smsld-F31_smsld*F13_smsld)/detF_smsld'}, ... 'invF23_smsld',{'','(F21_smsld*F13_smsld-F11_smsld*F23_smsld)/detF_smsld'}, ... 'invF31_smsld',{'','(F21_smsld*F32_smsld-F31_smsld*F22_smsld)/detF_smsld'}, ... 'invF32_smsld',{'','(F31_smsld*F12_smsld-F11_smsld*F32_smsld)/detF_smsld'}, ... 'invF33_smsld',{'','(F11_smsld*F22_smsld-F21_smsld*F12_smsld)/detF_smsld'}, ... 'tresca_smsld',{'','max(max(abs(s1_smsld-s2_smsld),abs(s2_smsld-s3_smsld)),abs(s1_smsld-s3_smsld))'}, ... 'xinit_ale',{'X','X'}, ... 'yinit_ale',{'Y','Y'}, ... 'zinit_ale',{'Z','Z'}, ... 'dx_ale',{'x-X','x-X'}, ... 'dy_ale',{'y-Y','y-Y'}, ... 'dz_ale',{'z-Z','z-Z'}, ... 'U_chns',{'sqrt(u2^2+v2^2+w2^2)',''}, ... 'Vx_chns',{'w2y-v2z',''}, ... 'Vy_chns',{'u2z-w2x',''}, ... 'Vz_chns',{'v2x-u2y',''}, ... 'V_chns',{'sqrt(Vx_chns^2+Vy_chns^2+Vz_chns^2)',''}, ... 'divU_chns',{'u2x+v2y+w2z',''}, ... 'cellRe_chns',{'rho_chns*U_chns*h_ale/eta_chns',''}, ... 'res_u2_chns',{'rho_chns*(u2*u2x+v2*u2y+w2*u2z)+p2x-F_x_chns-eta_chns*(2*u2xx+u2yy+v2xy+u2zz+w2xz)', ... ''}, ... 'res_tst_u2_chns',{'nojac(rho_chns)*(nojac(u2)*u2x+nojac(v2)*u2y+nojac(w2)*u2z)+p2x-nojac(eta_chns)*(2*u2xx+u2yy+v2xy+u2zz+w2xz)', ... ''}, ... 'res_sc_u2_chns',{'rho_chns*(u2*u2x+v2*u2y+w2*u2z)+p2x-F_x_chns', ... ''}, ... 'res_v2_chns',{'rho_chns*(u2*v2x+v2*v2y+w2*v2z)+p2y-F_y_chns-eta_chns*(v2xx+u2yx+2*v2yy+v2zz+w2yz)', ... ''}, ... 'res_tst_v2_chns',{'nojac(rho_chns)*(nojac(u2)*v2x+nojac(v2)*v2y+nojac(w2)*v2z)+p2y-nojac(eta_chns)*(v2xx+u2yx+2*v2yy+v2zz+w2yz)', ... ''}, ... 'res_sc_v2_chns',{'rho_chns*(u2*v2x+v2*v2y+w2*v2z)+p2y-F_y_chns', ... ''}, ... 'res_w2_chns',{'rho_chns*(u2*w2x+v2*w2y+w2*w2z)+p2z-F_z_chns-eta_chns*(w2xx+u2zx+w2yy+v2zy+2*w2zz)', ... ''}, ... 'res_tst_w2_chns',{'nojac(rho_chns)*(nojac(u2)*w2x+nojac(v2)*w2y+nojac(w2)*w2z)+p2z-nojac(eta_chns)*(w2xx+u2zx+w2yy+v2zy+2*w2zz)', ... ''}, ... 'res_sc_w2_chns',{'rho_chns*(u2*w2x+v2*w2y+w2*w2z)+p2z-F_z_chns', ... ''}, ... 'beta_x_chns',{'rho_chns*u2',''}, ... 'beta_y_chns',{'rho_chns*v2',''}, ... 'beta_z_chns',{'rho_chns*w2',''}, ... 'Dm_chns',{'eta_chns',''}, ... 'da_chns',{'rho_chns',''},'E_smsld',{'','31100'}, ... 'rho_smsld',{'','7850'}, ... 'alphadM_smsld',{'','1'}, ... 'betadK_smsld',{'','0.001'}, ... 'eta_smsld',{'','0'}, ... 'mixedform_smsld',{'','0'}, ... 'mu_smsld',{'','800000'}, ... 'C10_smsld',{'','200000'}, ... 'C01_smsld',{'','200000'}, ... 'kappa_smsld',{'','10000000000'}, ... 'ETkin_smsld',{'','20000000000'}, ... 'ETiso_smsld',{'','20000000000'}, ... 'Sys_smsld',{'','200000000'}, ... 'Syfunc_smsld',{'','mises_smsld'}, ... 'Syfunc_kin_smsld',{'','misesKin_smsld'}, ... 'Shard_smsld',{'','22222222222.22222*epe_smsld'}, ... 'alpha_smsld',{'','1.2e-005'}, ... 'nu_smsld',{'','0.49'}, ... 'EX_smsld',{'','200000000000'}, ... 'nuXY_smsld',{'','0.33'}, ... 'GXY_smsld',{'','75200000000'}, ... 'alphaX_smsld',{'','1.2e-005'}, ... 'sXi_smsld',{'','0'}, ... 'sXYi_smsld',{'','0'}, ... 'eXi_smsld',{'','0'}, ... 'eXYi_smsld',{'','0'}, ... 'EY_smsld',{'','200000000000'}, ... 'nuYZ_smsld',{'','0.33'}, ... 'GYZ_smsld',{'','75200000000'}, ... 'alphaY_smsld',{'','1.2e-005'}, ... 'sYi_smsld',{'','0'}, ... 'sYZi_smsld',{'','0'}, ... 'eYi_smsld',{'','0'}, ... 'eYZi_smsld',{'','0'}, ... 'EZ_smsld',{'','200000000000'}, ... 'nuXZ_smsld',{'','0.33'}, ... 'GXZ_smsld',{'','75200000000'}, ... 'alphaZ_smsld',{'','1.2e-005'}, ... 'sZi_smsld',{'','0'}, ... 'sXZi_smsld',{'','0'}, ... 'eZi_smsld',{'','0'}, ... 'eXZi_smsld',{'','0'}, ... 'FX_smsld',{'','0'}, ... 'FXAmp_smsld',{'','1'}, ... 'FXPh_smsld',{'','0'}, ... 'FY_smsld',{'','0'}, ... 'FYAmp_smsld',{'','1'}, ... 'FYPh_smsld',{'','0'}, ... 'FZ_smsld',{'','0'}, ... 'FZAmp_smsld',{'','1'}, ... 'FZPh_smsld',{'','0'}, ... 'Temp_smsld',{'','0'}, ... 'Tempref_smsld',{'','0'}, ... 'velx_smsld',{'','0'}, ... 'accx_smsld',{'','0'}, ... 'vely_smsld',{'','0'}, ... 'accy_smsld',{'','0'}, ... 'velz_smsld',{'','0'}, ... 'accz_smsld',{'','0'}, ... 'Rx_smsld',{'','0'}, ... 'Ry_smsld',{'','0'}, ... 'Rz_smsld',{'','0'}, ... 'dx_smsld',{'','1'}, ... 'Lx_smsld',{'','dx_smsld'}, ... 'dy_smsld',{'','1'}, ... 'Ly_smsld',{'','dy_smsld'}, ... 'dz_smsld',{'','1'}, ... 'Lz_smsld',{'','dz_smsld'}, ... 'dr_smsld',{'','1'}, ... 'Lr_smsld',{'','dr_smsld'}, ... 'R0_smsld',{'','1'}, ... 'x0_smsld',{'','0'}, ... 'y0_smsld',{'','0'}, ... 'z0_smsld',{'','0'}, ... 'H11_smsld',{'','0'}, ... 'H12_smsld',{'','0'}, ... 'H13_smsld',{'','0'}, ... 'H21_smsld',{'','0'}, ... 'H22_smsld',{'','0'}, ... 'H23_smsld',{'','0'}, ... 'H31_smsld',{'','0'}, ... 'H32_smsld',{'','0'}, ... 'H33_smsld',{'','0'}, ... 'rho_chns',{'1000',''}, ... 'eta_chns',{'0.001',''}, ... 'F_x_chns',{'0',''}, ... 'F_y_chns',{'0',''}, ... 'F_z_chns',{'0',''}, ... 'eta0_chns',{'1',''}, ... 'm_chns',{'0',''}, ... 'n_chns',{'0',''}, ... 'lambda_chns',{'0',''}, ... 'eta_inf_chns',{'0',''}, ... 'type_visc_chns',{'power',''}, ... 'kappadv_chns',{'0',''}, ... 'k_chns',{'1',''}, ... 'brinkmaneqns_chns',{'0',''}}; equ.expr = {}; equ.bnd.weak = {{'0';'0';'0';'0';'0';'0';'0';'0';'0';'0';'0'}}; equ.bnd.gporder = {{1;1;1;1;1;1;1;1;1;1;2}}; equ.bnd.ind = [1,1,1]; equ.bnd.expr = {}; equ.lock = [0,0,0]; equ.mlock = {[0,0,0]}; fem0.equ = equ; clear bnd bnd.weak = {{'0';'0';'0';'0';'0';'0';'0';'-lm4_test*(-x+X)+x_test*lm4';'-lm5_test*(-y+Y)+y_test*lm5'; ... '-lm6_test*(-z+Z)+z_test*lm6';'0';'0';'0';'0';'-lm7_test*(-u2+u0_chns)+u2_test*lm7'; ... 'lm8_test*v2+v2_test*lm8';'lm9_test*w2+w2_test*lm9'},{'0';'0';'0';'0';'0'; ... '0';'0';'-lm4_test*(-x+X)+x_test*lm4';'-lm5_test*(-y+Y)+y_test*lm5';'-lm6_test*(-z+Z)+z_test*lm6'; ... '0';'0';'0';'0';'0';'0';'0'},{'0';'0';'0';'0';'0';'0';'0';'-lm4_test*(-x+X)+x_test*lm4'; ... '-lm5_test*(-y+Y)+y_test*lm5';'-lm6_test*(-z+Z)+z_test*lm6';'0';'0';'0'; ... '0';'lm7_test*u2+u2_test*lm7';'lm8_test*v2+v2_test*lm8';'lm9_test*w2+w2_test*lm9'}, ... {'0';'0';'0';'0';'0';'0';'0';'0';'0';'0';'0';'0';'0';'0';'0';'0';'0'},{'FXg_smsld*u_test'; ... 'FYg_smsld*v_test';'FZg_smsld*w_test';'0';'0';'0';'0';'-lm4_test*(u-x+X)+x_test*lm4'; ... '-lm5_test*(v-y+Y)+y_test*lm5';'-lm6_test*(w-z+Z)+z_test*lm6';'0';'0';'0'; ... '0';'lm7_test*u2+u2_test*lm7';'lm8_test*v2+v2_test*lm8';'lm9_test*w2+w2_test*lm9'}, ... {'0';'0';'0';'0';'0';'0';'0';'0';'0';'0';'0';'0';'0';'0';'0';'0';'0'},{'FXg_smsld*u_test'; ... 'FYg_smsld*v_test';'FZg_smsld*w_test';'0';'0';'0';'0';'-lm4_test*(u-x+X)+x_test*lm4'; ... '-lm5_test*(v-y+Y)+y_test*lm5';'-lm6_test*(w-z+Z)+z_test*lm6';'0';'0';'0'; ... '0';'lm7_test*u2+u2_test*lm7';'lm8_test*v2+v2_test*lm8';'lm9_test*w2+w2_test*lm9'}}; bnd.dweak = {{'0';'0';'0';'0';'0';'0';'0';'0';'0';'0';'0';'0';'0';'0';'0'; ... '0';'0'},{'0';'0';'0';'0';'0';'0';'0';'0';'0';'0';'0';'0';'0';'0';'0';'0'; ... '0'},{'0';'0';'0';'0';'0';'0';'0';'0';'0';'0';'0';'0';'0';'0';'0';'0';'0'}, ... {'0';'0';'0';'0';'0';'0';'0';'0';'0';'0';'0';'0';'0';'0';'0';'0';'0'},{'0'; ... '0';'0';'0';'0';'0';'0';'0';'0';'0';'0';'0';'0';'0';'0';'0';'0'},{'0';'0'; ... '0';'0';'0';'0';'0';'0';'0';'0';'0';'0';'0';'0';'0';'0';'0'},{'0';'0';'0'; ... '0';'0';'0';'0';'0';'0';'0';'0';'0';'0';'0';'0';'0';'0'}}; bnd.constr = {{'0';'0';'0';'0';'0';'0';'0';'0';'0';'0';'0';'0';'0';'0';'0'; ... '0';'0'},{'0';'0';'0';'0';'0';'0';'0';'0';'0';'0';'0';'0';'0';'0';'0';'0'; ... '0'},{'0';'0';'0';'0';'0';'0';'0';'0';'0';'0';'0';'0';'0';'0';'0';'0';'0'}, ... {'0';'0';'0';'0';'0';'0';'0';'0';'0';'0';'0';'0';'0';'0';'0';'0';'0'},{'0'; ... '0';'0';'0';'0';'0';'0';'0';'0';'0';'0';'0';'0';'0';'0';'0';'0'},{'-u'; ... '-v';'-w';'0';'0';'0';'0';'0';'0';'0';'0';'0';'0';'0';'0';'0';'0'},{'0'; ... '0';'0';'0';'0';'0';'0';'0';'0';'0';'0';'0';'0';'0';'0';'0';'0'}}; bnd.shape = {[1;2;3;11;12;13;4;5;6;7;14;15;16],[1;2;3;11;12;13;4;5;6;7],[1; ... 2;3;11;12;13;4;5;6;7;14;15;16],[1;2;3;4;5;6;7],[8;9;10;1;2;3;11;12;13;4; ... 5;6;7;14;15;16],[8;9;10],[8;9;10;1;2;3;11;12;13;4;5;6;7;14;15;16]}; bnd.sshape = {[1;2],[1;2],[1;2],[1;2],[1;2],[1;3],[1;2]}; bnd.sshapedim = {{1;1;1;1;2;2;2;2;2;2;2;2;2;2;2;2;2},{1;1;1;1;2;2;2;2;2;2;2;2; ... 2;2;2;2;2},{1;1;1;1;2;2;2;2;2;2;2;2;2;2;2;2;2},{1;1;1;1;2;2;2;2;2;2;2;2; ... 2;2;2;2;2},{1;1;1;1;2;2;2;2;2;2;2;2;2;2;2;2;2},{1;1;1;1;3;3;3;3;3;3;3;3; ... 3;3;3;3;3},{1;1;1;1;2;2;2;2;2;2;2;2;2;2;2;2;2}}; bnd.gporder = {{1;1;1;1;1;1;1;3;3;3;1;1;1;2;3;3;3},{1;1;1;1;1;1;1;3;3;3;1;1; ... 1;2;3;3;3},{1;1;1;1;1;1;1;3;3;3;1;1;1;2;3;3;3},{1;1;1;1;1;1;1;3;3;3;1;1; ... 1;2;3;3;3},{1;1;1;1;1;1;1;3;3;3;1;1;1;2;3;3;3},{1;1;1;1;1;1;1;3;3;3;1;1; ... 1;2;3;3;3},{1;1;1;1;1;1;1;3;3;3;1;1;1;2;3;3;3}}; bnd.cporder = {{1;1;1;1;1;1;1;1;1;1;1;1;1;2;1;1;1},{1;1;1;1;1;1;1;1;1;1;1;1; ... 1;2;1;1;1},{1;1;1;1;1;1;1;1;1;1;1;1;1;2;1;1;1},{1;1;1;1;1;1;1;1;1;1;1;1; ... 1;2;1;1;1},{1;1;1;1;1;1;1;1;1;1;1;1;1;2;1;1;1},{1;1;1;1;1;1;1;1;1;1;1;1; ... 1;2;1;1;1},{1;1;1;1;1;1;1;1;1;1;1;1;1;2;1;1;1}}; bnd.init = {{'';'';'';'';'';'';'';'0';'0';'0';'';'';'';'';'0';'0';'0'},{''; ... '';'';'';'';'';'';'0';'0';'0';'';'';'';'';'0';'0';'0'},{'';'';'';'';''; ... '';'';'0';'0';'0';'';'';'';'';'0';'0';'0'},{'';'';'';'';'';'';'';'0';'0'; ... '0';'';'';'';'';'0';'0';'0'},{'';'';'';'';'';'';'';'0';'0';'0';'';'';''; ... '';'0';'0';'0'},{'';'';'';'';'';'';'';'0';'0';'0';'';'';'';'';'0';'0';'0'}, ... {'';'';'';'';'';'';'';'0';'0';'0';'';'';'';'';'0';'0';'0'}}; bnd.dinit = {{'';'';'';'';'';'';'';'';'';'';'';'';'';'';'';'';''},{'';'';''; ... '';'';'';'';'';'';'';'';'';'';'';'';'';''},{'';'';'';'';'';'';'';'';''; ... '';'';'';'';'';'';'';''},{'';'';'';'';'';'';'';'';'';'';'';'';'';'';''; ... '';''},{'';'';'';'';'';'';'';'';'';'';'';'';'';'';'';'';''},{'';'';'';''; ... '';'';'';'';'';'';'';'';'';'';'';'';''},{'';'';'';'';'';'';'';'';'';''; ... '';'';'';'';'';'';''}}; bnd.ind = [1,2,2,2,2,2,3,3,3,4,3,5,5,6,6,7,7,2]; bnd.dim = {'u','v','w','p','x','y','z','lm4','lm5','lm6','u2','v2','w2', ... 'p2','lm7','lm8','lm9'}; bnd.var = {'FXg_smsld',{'','','','','FX_smsld','0','FX_smsld'}, ... 'FYg_smsld',{'','','','','FY_smsld','0','FY_smsld'}, ... 'FZg_smsld',{'','','','','FZ_smsld','0','FZ_smsld'}, ... 'disp_smsld',{'','','','','sqrt(real(u)^2+real(v)^2+real(w)^2)','sqrt(real(u)^2+real(v)^2+real(w)^2)', ... 'sqrt(real(u)^2+real(v)^2+real(w)^2)'}, ... 'TaX_smsld',{'','','','','(F11_smsld*SX_smsld+F12_smsld*SXY_smsld+F13_smsld*SXZ_smsld)*nX_smsld+(F11_smsld*SXY_smsld+F12_smsld*SY_smsld+F13_smsld*SYZ_smsld)*nY_smsld+(F11_smsld*SXZ_smsld+F12_smsld*SYZ_smsld+F13_smsld*SZ_smsld)*nZ_smsld', ... '(F11_smsld*SX_smsld+F12_smsld*SXY_smsld+F13_smsld*SXZ_smsld)*nX_smsld+(F11_smsld*SXY_smsld+F12_smsld*SY_smsld+F13_smsld*SYZ_smsld)*nY_smsld+(F11_smsld*SXZ_smsld+F12_smsld*SYZ_smsld+F13_smsld*SZ_smsld)*nZ_smsld', ... '(F11_smsld*SX_smsld+F12_smsld*SXY_smsld+F13_smsld*SXZ_smsld)*nX_smsld+(F11_smsld*SXY_smsld+F12_smsld*SY_smsld+F13_smsld*SYZ_smsld)*nY_smsld+(F11_smsld*SXZ_smsld+F12_smsld*SYZ_smsld+F13_smsld*SZ_smsld)*nZ_smsld'}, ... 'TaY_smsld',{'','','','','(F21_smsld*SX_smsld+F22_smsld*SXY_smsld+F23_smsld*SXZ_smsld)*nX_smsld+(F21_smsld*SXY_smsld+F22_smsld*SY_smsld+F23_smsld*SYZ_smsld)*nY_smsld+(F21_smsld*SXZ_smsld+F22_smsld*SYZ_smsld+F23_smsld*SZ_smsld)*nZ_smsld', ... '(F21_smsld*SX_smsld+F22_smsld*SXY_smsld+F23_smsld*SXZ_smsld)*nX_smsld+(F21_smsld*SXY_smsld+F22_smsld*SY_smsld+F23_smsld*SYZ_smsld)*nY_smsld+(F21_smsld*SXZ_smsld+F22_smsld*SYZ_smsld+F23_smsld*SZ_smsld)*nZ_smsld', ... '(F21_smsld*SX_smsld+F22_smsld*SXY_smsld+F23_smsld*SXZ_smsld)*nX_smsld+(F21_smsld*SXY_smsld+F22_smsld*SY_smsld+F23_smsld*SYZ_smsld)*nY_smsld+(F21_smsld*SXZ_smsld+F22_smsld*SYZ_smsld+F23_smsld*SZ_smsld)*nZ_smsld'}, ... 'TaZ_smsld',{'','','','','(F31_smsld*SX_smsld+F32_smsld*SXY_smsld+F33_smsld*SXZ_smsld)*nX_smsld+(F31_smsld*SXY_smsld+F32_smsld*SY_smsld+F33_smsld*SYZ_smsld)*nY_smsld+(F31_smsld*SXZ_smsld+F32_smsld*SYZ_smsld+F33_smsld*SZ_smsld)*nZ_smsld', ... '(F31_smsld*SX_smsld+F32_smsld*SXY_smsld+F33_smsld*SXZ_smsld)*nX_smsld+(F31_smsld*SXY_smsld+F32_smsld*SY_smsld+F33_smsld*SYZ_smsld)*nY_smsld+(F31_smsld*SXZ_smsld+F32_smsld*SYZ_smsld+F33_smsld*SZ_smsld)*nZ_smsld', ... '(F31_smsld*SX_smsld+F32_smsld*SXY_smsld+F33_smsld*SXZ_smsld)*nX_smsld+(F31_smsld*SXY_smsld+F32_smsld*SY_smsld+F33_smsld*SYZ_smsld)*nY_smsld+(F31_smsld*SXZ_smsld+F32_smsld*SYZ_smsld+F33_smsld*SZ_smsld)*nZ_smsld'}, ... 'K_x_chns',{'2*nx_chns*eta_chns*u2x+ny_chns*eta_chns*(u2y+v2x)+nz_chns*eta_chns*(u2z+w2x)', ... '2*nx_chns*eta_chns*u2x+ny_chns*eta_chns*(u2y+v2x)+nz_chns*eta_chns*(u2z+w2x)', ... '2*nx_chns*eta_chns*u2x+ny_chns*eta_chns*(u2y+v2x)+nz_chns*eta_chns*(u2z+w2x)', ... '2*nx_chns*eta_chns*u2x+ny_chns*eta_chns*(u2y+v2x)+nz_chns*eta_chns*(u2z+w2x)', ... '2*nx_chns*eta_chns*u2x+ny_chns*eta_chns*(u2y+v2x)+nz_chns*eta_chns*(u2z+w2x)', ... '','2*nx_chns*eta_chns*u2x+ny_chns*eta_chns*(u2y+v2x)+nz_chns*eta_chns*(u2z+w2x)'}, ... 'T_x_chns',{'-nx_chns*p2+2*nx_chns*eta_chns*u2x+ny_chns*eta_chns*(u2y+v2x)+nz_chns*eta_chns*(u2z+w2x)', ... '-nx_chns*p2+2*nx_chns*eta_chns*u2x+ny_chns*eta_chns*(u2y+v2x)+nz_chns*eta_chns*(u2z+w2x)', ... '-nx_chns*p2+2*nx_chns*eta_chns*u2x+ny_chns*eta_chns*(u2y+v2x)+nz_chns*eta_chns*(u2z+w2x)', ... '-nx_chns*p2+2*nx_chns*eta_chns*u2x+ny_chns*eta_chns*(u2y+v2x)+nz_chns*eta_chns*(u2z+w2x)', ... '-nx_chns*p2+2*nx_chns*eta_chns*u2x+ny_chns*eta_chns*(u2y+v2x)+nz_chns*eta_chns*(u2z+w2x)', ... '','-nx_chns*p2+2*nx_chns*eta_chns*u2x+ny_chns*eta_chns*(u2y+v2x)+nz_chns*eta_chns*(u2z+w2x)'}, ... 'K_y_chns',{'nx_chns*eta_chns*(v2x+u2y)+2*ny_chns*eta_chns*v2y+nz_chns*eta_chns*(v2z+w2y)', ... 'nx_chns*eta_chns*(v2x+u2y)+2*ny_chns*eta_chns*v2y+nz_chns*eta_chns*(v2z+w2y)', ... 'nx_chns*eta_chns*(v2x+u2y)+2*ny_chns*eta_chns*v2y+nz_chns*eta_chns*(v2z+w2y)', ... 'nx_chns*eta_chns*(v2x+u2y)+2*ny_chns*eta_chns*v2y+nz_chns*eta_chns*(v2z+w2y)', ... 'nx_chns*eta_chns*(v2x+u2y)+2*ny_chns*eta_chns*v2y+nz_chns*eta_chns*(v2z+w2y)', ... '','nx_chns*eta_chns*(v2x+u2y)+2*ny_chns*eta_chns*v2y+nz_chns*eta_chns*(v2z+w2y)'}, ... 'T_y_chns',{'-ny_chns*p2+nx_chns*eta_chns*(v2x+u2y)+2*ny_chns*eta_chns*v2y+nz_chns*eta_chns*(v2z+w2y)', ... '-ny_chns*p2+nx_chns*eta_chns*(v2x+u2y)+2*ny_chns*eta_chns*v2y+nz_chns*eta_chns*(v2z+w2y)', ... '-ny_chns*p2+nx_chns*eta_chns*(v2x+u2y)+2*ny_chns*eta_chns*v2y+nz_chns*eta_chns*(v2z+w2y)', ... '-ny_chns*p2+nx_chns*eta_chns*(v2x+u2y)+2*ny_chns*eta_chns*v2y+nz_chns*eta_chns*(v2z+w2y)', ... '-ny_chns*p2+nx_chns*eta_chns*(v2x+u2y)+2*ny_chns*eta_chns*v2y+nz_chns*eta_chns*(v2z+w2y)', ... '','-ny_chns*p2+nx_chns*eta_chns*(v2x+u2y)+2*ny_chns*eta_chns*v2y+nz_chns*eta_chns*(v2z+w2y)'}, ... 'K_z_chns',{'nx_chns*eta_chns*(w2x+u2z)+ny_chns*eta_chns*(w2y+v2z)+2*nz_chns*eta_chns*w2z', ... 'nx_chns*eta_chns*(w2x+u2z)+ny_chns*eta_chns*(w2y+v2z)+2*nz_chns*eta_chns*w2z', ... 'nx_chns*eta_chns*(w2x+u2z)+ny_chns*eta_chns*(w2y+v2z)+2*nz_chns*eta_chns*w2z', ... 'nx_chns*eta_chns*(w2x+u2z)+ny_chns*eta_chns*(w2y+v2z)+2*nz_chns*eta_chns*w2z', ... 'nx_chns*eta_chns*(w2x+u2z)+ny_chns*eta_chns*(w2y+v2z)+2*nz_chns*eta_chns*w2z', ... '','nx_chns*eta_chns*(w2x+u2z)+ny_chns*eta_chns*(w2y+v2z)+2*nz_chns*eta_chns*w2z'}, ... 'T_z_chns',{'-nz_chns*p2+nx_chns*eta_chns*(w2x+u2z)+ny_chns*eta_chns*(w2y+v2z)+2*nz_chns*eta_chns*w2z', ... '-nz_chns*p2+nx_chns*eta_chns*(w2x+u2z)+ny_chns*eta_chns*(w2y+v2z)+2*nz_chns*eta_chns*w2z', ... '-nz_chns*p2+nx_chns*eta_chns*(w2x+u2z)+ny_chns*eta_chns*(w2y+v2z)+2*nz_chns*eta_chns*w2z', ... '-nz_chns*p2+nx_chns*eta_chns*(w2x+u2z)+ny_chns*eta_chns*(w2y+v2z)+2*nz_chns*eta_chns*w2z', ... '-nz_chns*p2+nx_chns*eta_chns*(w2x+u2z)+ny_chns*eta_chns*(w2y+v2z)+2*nz_chns*eta_chns*w2z', ... '','-nz_chns*p2+nx_chns*eta_chns*(w2x+u2z)+ny_chns*eta_chns*(w2y+v2z)+2*nz_chns*eta_chns*w2z'},'FX_smsld',{'','','','','-lm7*dvol_ale/dvol','0','-lm7*dvol_ale/dvol'}, ... 'FXAmp_smsld',{'','','','','1','1','1'}, ... 'FXPh_smsld',{'','','','','0','0','0'}, ... 'FY_smsld',{'','','','','-lm8*dvol_ale/dvol','0','-lm8*dvol_ale/dvol'}, ... 'FYAmp_smsld',{'','','','','1','1','1'}, ... 'FYPh_smsld',{'','','','','0','0','0'}, ... 'FZ_smsld',{'','','','','-lm9*dvol_ale/dvol','0','-lm9*dvol_ale/dvol'}, ... 'FZAmp_smsld',{'','','','','1','1','1'}, ... 'FZPh_smsld',{'','','','','0','0','0'}, ... 'velx_smsld',{'','','','','0','0','0'}, ... 'accx_smsld',{'','','','','0','0','0'}, ... 'vely_smsld',{'','','','','0','0','0'}, ... 'accy_smsld',{'','','','','0','0','0'}, ... 'velz_smsld',{'','','','','0','0','0'}, ... 'accz_smsld',{'','','','','0','0','0'}, ... 'Rx_smsld',{'','','','','0','0','0'}, ... 'Ry_smsld',{'','','','','0','0','0'}, ... 'Rz_smsld',{'','','','','0','0','0'}, ... 'nX_smsld',{'','','','','unX','nX','dnX'}, ... 'nY_smsld',{'','','','','unY','nY','dnY'}, ... 'nZ_smsld',{'','','','','unZ','nZ','dnZ'}, ... 'H11_smsld',{'','','','','0','1','0'}, ... 'H12_smsld',{'','','','','0','0','0'}, ... 'H13_smsld',{'','','','','0','0','0'}, ... 'H21_smsld',{'','','','','0','0','0'}, ... 'H22_smsld',{'','','','','0','1','0'}, ... 'H23_smsld',{'','','','','0','0','0'}, ... 'H31_smsld',{'','','','','0','0','0'}, ... 'H32_smsld',{'','','','','0','0','0'}, ... 'H33_smsld',{'','','','','0','1','0'}, ... 'nX_ale',{'nX','nX','nX','nX','dnX','','unX'}, ... 'nY_ale',{'nY','nY','nY','nY','dnY','','unY'}, ... 'nZ_ale',{'nZ','nZ','nZ','nZ','dnZ','','unZ'}, ... 'u0_chns',{'0.5','0','0','0','0','','0'}, ... 'v0_chns',{'0','0','0','0','0','','0'}, ... 'w0_chns',{'0','0','0','0','0','','0'}, ... 'p0_chns',{'0','0','0','0','0','','0'}, ... 'k0_chns',{'0','0','0','0','0','','0'}, ... 'd0_chns',{'0','0','0','0','0','','0'}, ... 'LT_chns',{'0.07000000000000001','0.07000000000000001','0.07000000000000001', ... '0.07000000000000001','0.07000000000000001','','0.07000000000000001'}, ... 'IT_chns',{'0.054','0.054','0.054','0.054','0.054','','0.054'}, ... 'uw_chns',{'0','0','0','0','0','','0'}, ... 'vw_chns',{'0','0','0','0','0','','0'}, ... 'ww_chns',{'0','0','0','0','0','','0'}, ... 'dw_chns',{'h','h','h','h','h','','h'}, ... 'yplus_chns',{'100','100','100','100','100','','100'}, ... 'U0_chns',{'0','0','0','0','0','','0'}, ... 'V0_chns',{'0','0','0','0','0','','0'}, ... 'p0_inl_chns',{'0','0','0','0','0','','0'}, ... 'L0_chns',{'1','1','1','1','1','','1'}, ... 'nx_chns',{'nx','nx','nx','nx','dnx','','unx'}, ... 'ny_chns',{'ny','ny','ny','ny','dny','','uny'}, ... 'nz_chns',{'nz','nz','nz','nz','dnz','','unz'}}; bnd.expr = {}; bnd.lock = [0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0]; bnd.mlock = {[0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0]}; fem0.bnd = bnd; clear edg edg.weak = {{'0';'0';'0';'0';'0';'0';'0';'0';'0';'0';'0'},{'0';'0';'0';'0'; ... '0';'0';'0';'0';'0';'0';'0'}}; edg.dweak = {{'0';'0';'0';'0';'0';'0';'0';'0';'0';'0';'0'},{'0';'0';'0';'0'; ... '0';'0';'0';'0';'0';'0';'0'}}; edg.constr = {{'0';'0';'0';'0';'0';'0';'0';'0';'0';'0';'0'},{'0';'0';'0';'0'; ... '0';'0';'0';'0';'0';'0';'0'}}; edg.shape = {[1;2;3;4;5;6;7],[8;9;10;1;2;3;4;5;6;7]}; edg.sshape = {[1;2],[1;2]}; edg.sshapedim = {{1;1;1;1;2;2;2;2;2;2;2},{1;1;1;1;2;2;2;2;2;2;2}}; edg.gporder = {{1;1;1;1;1;1;1;1;1;1;2},{1;1;1;1;1;1;1;1;1;1;2}}; edg.cporder = {{1;1;1;1;1;1;1;1;1;1;2},{1;1;1;1;1;1;1;1;1;1;2}}; edg.init = {{'';'';'';'';'';'';'';'';'';'';''},{'';'';'';'';'';'';'';'';''; ... '';''}}; edg.dinit = {{'';'';'';'';'';'';'';'';'';'';''},{'';'';'';'';'';'';'';'';''; ... '';''}}; edg.ind = [1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,2,2,2,2,2,2,2,2,2,2, ... 2,2,1,1,1,1]; edg.dim = {'u','v','w','p','x','y','z','u2','v2','w2','p2'}; edg.var = {'FXg_smsld',{'','0'}, ... 'FYg_smsld',{'','0'}, ... 'FZg_smsld',{'','0'}, ... 'disp_smsld',{'','sqrt(real(u)^2+real(v)^2+real(w)^2)'},'FX_smsld',{'','0'}, ... 'FXAmp_smsld',{'','1'}, ... 'FXPh_smsld',{'','0'}, ... 'FY_smsld',{'','0'}, ... 'FYAmp_smsld',{'','1'}, ... 'FYPh_smsld',{'','0'}, ... 'FZ_smsld',{'','0'}, ... 'FZAmp_smsld',{'','1'}, ... 'FZPh_smsld',{'','0'}, ... 'Rx_smsld',{'','0'}, ... 'Ry_smsld',{'','0'}, ... 'Rz_smsld',{'','0'}, ... 'H11_smsld',{'','0'}, ... 'H12_smsld',{'','0'}, ... 'H13_smsld',{'','0'}, ... 'H21_smsld',{'','0'}, ... 'H22_smsld',{'','0'}, ... 'H23_smsld',{'','0'}, ... 'H31_smsld',{'','0'}, ... 'H32_smsld',{'','0'}, ... 'H33_smsld',{'','0'}}; edg.expr = {}; edg.lock = [0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0]; edg.mlock = {[0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0]}; fem0.edg = edg; clear pnt pnt.weak = {{'0';'0';'0';'0';'0';'0';'0';'0';'0';'0';'0';'0'},{'0';'0';'0'; ... '0';'0';'0';'0';'0';'0';'0';'0';'0'}}; pnt.dweak = {{'0';'0';'0';'0';'0';'0';'0';'0';'0';'0';'0';'0'},{'0';'0';'0'; ... '0';'0';'0';'0';'0';'0';'0';'0';'0'}}; pnt.constr = {{'0';'0';'0';'0';'0';'0';'0';'0';'0';'0';'0';'0'},{'0';'0';'0'; ... '0';'0';'0';'0';'0';'0';'0';'0';'0'}}; pnt.shape = {[1;2;3;4;5;6;7],[8;9;10;1;2;3;4;5;6;7]}; pnt.sshape = {[1;2],[1;2]}; pnt.sshapedim = {{1;1;1;1;2;2;2;2;2;2;2;2},{1;1;1;1;2;2;2;2;2;2;2;2}}; pnt.init = {{'';'';'';'';'';'';'';'';'';'';'';''},{'';'';'';'';'';'';'';''; ... '';'';'';''}}; pnt.dinit = {{'';'';'';'';'';'';'';'';'';'';'';''},{'';'';'';'';'';'';'';''; ... '';'';'';''}}; pnt.ind = [1,1,1,1,1,1,1,1,1,1,1,1,2,2,2,2,2,2,2,2,1,1,1,1]; pnt.dim = {'u','v','w','p','x','y','z','u2','v2','w2','p2','Pinl_chns'}; pnt.var = {'FXg_smsld',{'','0'}, ... 'FYg_smsld',{'','0'}, ... 'FZg_smsld',{'','0'}, ... 'disp_smsld',{'','sqrt(real(u)^2+real(v)^2+real(w)^2)'},'FX_smsld',{'','0'}, ... 'FXAmp_smsld',{'','1'}, ... 'FXPh_smsld',{'','0'}, ... 'FY_smsld',{'','0'}, ... 'FYAmp_smsld',{'','1'}, ... 'FYPh_smsld',{'','0'}, ... 'FZ_smsld',{'','0'}, ... 'FZAmp_smsld',{'','1'}, ... 'FZPh_smsld',{'','0'}, ... 'Rx_smsld',{'','0'}, ... 'Ry_smsld',{'','0'}, ... 'Rz_smsld',{'','0'}, ... 'H11_smsld',{'','0'}, ... 'H12_smsld',{'','0'}, ... 'H13_smsld',{'','0'}, ... 'H21_smsld',{'','0'}, ... 'H22_smsld',{'','0'}, ... 'H23_smsld',{'','0'}, ... 'H31_smsld',{'','0'}, ... 'H32_smsld',{'','0'}, ... 'H33_smsld',{'','0'}, ... 'p0_chns',{'0','0'}}; pnt.expr = {}; pnt.lock = [0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0]; pnt.mlock = {[0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0]}; fem0.pnt = pnt; fem0.var = {'t_old_ini_smsld','-1','ut','uTIME','utt','uTIMETIME','uXt','uXTIME','uYt','uYTIME','uZt','uZTIME','vt','vTIME','vtt','vTIMETIME','vXt','vXTIME','vYt','vYTIME','vZt','vZTIME','wt','wTIME','wtt','wTIMETIME','wXt','wXTIME','wYt','wYTIME','wZt','wZTIME','xt','xTIME','xtt','xTIMETIME','xxt','xxTIME','xyt','xyTIME','xzt','xzTIME','yt','yTIME','ytt','yTIMETIME','yxt','yxTIME','yyt','yyTIME','yzt','yzTIME','zt','zTIME','ztt','zTIMETIME','zxt','zxTIME','zyt','zyTIME','zzt','zzTIME','u2t','u2TIME-xTIME*u2x-yTIME*u2y-zTIME*u2z','u2tt','u2TIMETIME','u2xt','u2xTIME-xTIME*u2xx-yTIME*u2xy-zTIME*u2xz','u2yt','u2yTIME-xTIME*u2yx-yTIME*u2yy-zTIME*u2yz','u2zt','u2zTIME-xTIME*u2zx-yTIME*u2zy-zTIME*u2zz','v2t','v2TIME-xTIME*v2x-yTIME*v2y-zTIME*v2z','v2tt','v2TIMETIME','v2xt','v2xTIME-xTIME*v2xx-yTIME*v2xy-zTIME*v2xz','v2yt','v2yTIME-xTIME*v2yx-yTIME*v2yy-zTIME*v2yz','v2zt','v2zTIME-xTIME*v2zx-yTIME*v2zy-zTIME*v2zz','w2t','w2TIME-xTIME*w2x-yTIME*w2y-zTIME*w2z','w2tt','w2TIMETIME','w2xt','w2xTIME-xTIME*w2xx-yTIME*w2xy-zTIME*w2xz','w2yt','w2yTIME-xTIME*w2yx-yTIME*w2yy-zTIME*w2yz','w2zt','w2zTIME-xTIME*w2zx-yTIME*w2zy-zTIME*w2zz','p2t','p2TIME-xTIME*p2x-yTIME*p2y-zTIME*p2z','p2tt','p2TIMETIME','p2xt','p2xTIME-xTIME*p2xx-yTIME*p2xy-zTIME*p2xz','p2yt','p2yTIME-xTIME*p2yx-yTIME*p2yy-zTIME*p2yz','p2zt','p2zTIME-xTIME*p2zx-yTIME*p2zy-zTIME*p2zz'}; fem0.expr = {}; clear elemmph clear elem elem.elem = 'elpric'; elem.g = {'1'}; clear equ equ.tensor = {{'real(sX_smsld)','real(sXY_smsld)','real(sY_smsld)','real(sXZ_smsld)', ... 'real(sYZ_smsld)','real(sZ_smsld)'}}; equ.ind = {{'3'}}; elem.geomdim = {{{},{},{},equ}}; elem.basename = 's'; elem.postname = 'smsld'; elem.sdim = {'X','Y','Z'}; elemmph{1} = elem; clear elem elem.elem = 'elpric'; elem.g = {'1'}; clear equ equ.tensor = {{'real(eX_smsld)','real(eXY_smsld)','real(eY_smsld)','real(eXZ_smsld)', ... 'real(eYZ_smsld)','real(eZ_smsld)'}}; equ.ind = {{'3'}}; elem.geomdim = {{{},{},{},equ}}; elem.basename = 'e'; elem.postname = 'smsld'; elem.sdim = {'X','Y','Z'}; elemmph{2} = elem; clear elem elem.elem = 'elcplscalar'; elem.g = {'1'}; src = cell(1,1); clear equ equ.expr = {{'-qual_ale'}}; equ.ipoints = {{'4'}}; equ.ind = {{'1','2','3'}}; src{1} = {{},{},{},equ}; elem.src = src; geomdim = cell(1,1); geomdim{1} = {}; elem.geomdim = geomdim; elem.var = {'minqual_tmp1_ale'}; elem.global = {'1'}; elem.maxvars = {'1'}; elemmph{3} = elem; clear elem elem.elem = 'elconst'; elem.var = {'minqual1_ale','-minqual_tmp1_ale'}; elemmph{4} = elem; clear elem elem.elem = 'elcplextr'; elem.g = {'1'}; src = cell(1,1); src{1} = {{}}; elem.src = src; geomdim = cell(1,1); geomdim{1} = {{},{},{},{}}; elem.geomdim = geomdim; elem.var = {}; map = cell(1,0); elem.map = map; elemmph{5} = elem; clear elem elem.elem = 'elcplscalar'; elem.g = {'1'}; src = cell(1,1); src{1} = {{},{},{},{}}; elem.src = src; geomdim = cell(1,1); geomdim{1} = {{},{},{},{}}; elem.geomdim = geomdim; elem.var = {}; elem.global = {}; elemmph{6} = elem; fem0.elemmph = elemmph; clear draw draw.p.objs = {}; draw.p.name = {}; draw.c.objs = {}; draw.c.name = {}; draw.f.objs = {}; draw.f.name = {}; draw.s.objs = {}; draw.s.name = {}; fem0.draw = draw; fem0.const = {}; fem0.globalexpr = {}; clear fcns fem0.functions = {}; fem0.sol = flbinary('fem0','solution','tutorial5.mph'); fem0.xmcases = [0]; fem0.mcases = [0]; flbinary clear; fem0.rulingmode = 'smsld'; fem0.solform = 'weak'; sr com.femlab.geom.Geom3]�(�.� xr com.femlab.geom.Geom�6��{ L geomAssoct Lcom/femlab/geom/GeomAssoc;L versionq ~ xpw q ~ w H�ur [B���T� xp H� G e o m 3 �����|�=-C��6? �h㈵���'�>�?�'�>�?���������|�=�h㈵���'�>�?t^c���?���������|�=�h㈵��v^c���?�'�>�?���������|�=�h㈵��v^c���?t^c���?���������|�= ���������|�=!�rh��̻ -C��6*?���������|�= �'�>�?�'�>�?���������|�= �'�>�?v^c���?���������|�= v^c���?�'�>�?���������|�= v^c���?v^c���?���������|�= -C��6*? ���������|�=!�rh��̻-C��6*?-C��6*?���������|�=-C��6*? !�rh���;���������|�=-C��6*? -C��6*?���������|�=-C��6*?-C��6*?!�rh���;���������|�=-C��6*?-C��6*?-C��6*?���������|�=4 ���� � ���� � ���� � �?���� � �? ���� � �? �?���� � �?���� � �? ���� � �? �?���� � �? �?���� � �? ���� � �?���� � ���� � ���� � ���� � �?���� � �? ���� � �? �?���� � 333333�?433333�?���� � �? ���� � �?���� � �? �?���� � 333333�?gfffff�?���� � �? �?���� � �?���� � �? ���� � gfffff�?433333�?���� � �? �?���� � �?���� � �? ���� � gfffff�?gfffff�?���� � ���� � ���� � ���� � �?���� � �? ���� � �? �?���� � �? �?���� � �? ���� � �?���� � �? ���� � �?���� � �? �?���� � �? �?���� � �?���� � �? ���� � �? �?���� � �?���� � �? ���� � ���� � ���� � ���� � �?����! � �? �������� � �?���� � �? �������� � �?���� � �?����# � �? �������� � �? �������� � �?���� � �? �������� � �?���� � �? �������� � �?���� � �?����" � �?���� � �?���� � �?����$ � �?���� � �? �������� � �? �������� � �?���� � �?���� � �?���� � �?���� �4 �? � �? � �? ���� � �? ���� � �? � �? � �? ���� � �? ���� � �? � �? � �? � �? � �? ���� � �? ���� � �? ���� � " �? ���� � �? � �? � # �? ���� � $ �? ���� � ) �? � * �? � &