sr!com.femlab.server.ModelFileHeaderD%LtagstLcom/femlab/util/FlStringList;Ltypesq~LvrsntLcom/femlab/util/FlVersion;xpwsrcom.femlab.util.FlVersion%/B = IbuildImajorLdatetLjava/lang/String;Lextq~Lnameq~Lrcsq~L reactionExtq~L reactionNameq~L scriptExtq~L scriptNameq~xpwtCOMSOL Script 1.2tt#COMSOL Reaction Engineering Lab 1.4q~t COMSOL 3.4q~wt $Name: $t$Date: 2007/10/10 16:07:51 $xur[Ljava.lang.String;V{Gxpt modelinfotxfemtguitg8tg20tg1tfem1t mfileinfouq~ q~t femstructt guistructtdrawq~q~tgeomq~xsrcom.femlab.api.client.ModelInfo^%Ldescrq~LdocURLq~[imaget[Bxpwptpxuq~ tclear xfem clear vrsn vrsn.name = 'COMSOL 3.4'; vrsn.ext = ''; vrsn.major = 0; vrsn.build = 248; vrsn.rcs = '$Name: $'; vrsn.date = '$Date: 2007/10/10 16:07:51 $'; xfem.version = vrsn; xfem.id = 1; xfem.geomdata = 'geom'; xfem.eqvars = 'on'; xfem.cplbndeq = 'on'; xfem.cplbndsh = 'off'; xfem.drawvalid = 'on'; xfem.geomvalid = 'on'; xfem.solvalid = 'off'; xfem.linshape = 'on'; xfem.linshapetol = 0.1; xfem.meshtime = 't'; clear appl appl.mode.class = 'Darcy'; appl.mode.type = 'axi'; appl.dim = {'p'}; appl.sdim = {'r','phi','z'}; appl.name = 'esdl'; appl.module = 'ES'; appl.shape = {'shlag(2,''p'')'}; appl.gporder = 4; appl.cporder = 2; appl.sshape = 2; appl.border = 'on'; appl.assignsuffix = '_esdl'; clear prop prop.elemdefault='Lag2'; prop.variable='H'; prop.analysis='time'; prop.frame='ref'; clear weakconstr weakconstr.value = 'off'; weakconstr.dim = {'lm1'}; prop.weakconstr = weakconstr; prop.constrtype='ideal'; appl.prop = prop; clear pnt pnt.p0 = {{}}; pnt.tp0 = {{}}; pnt.tpf = {{}}; pnt.N0 = {{}}; pnt.tN0 = {{}}; pnt.tNf = {{}}; pnt.style = {{{'0'},{'0','0','0'}}}; pnt.ind = [1,1,1,1,1,1,1,1]; appl.pnt = pnt; clear bnd bnd.name = {'','','','',''}; bnd.H0 = {'0','0','0','0','0'}; bnd.Hp0 = {'0','0','0','0','0'}; bnd.p0 = {'0','0','0','0','0'}; bnd.N0 = {'0','0','1e-10','0','0'}; bnd.Hb = {'0','0','0','0','0'}; bnd.Hpb = {'0','0','0','0','0'}; bnd.pb = {'0','0','0','0','0'}; bnd.Rb = {'0','0','0','0','0'}; bnd.Db = {'0','0','0','0','0'}; bnd.type = {'Na','Nz','N0','cont','H0'}; bnd.style = {{{'0'},{'0','255','0'},{'solid'}},{{'0'},{'0','0','0'},{'solid'}}, ... {{'0'},{'255','0','255'},{'solid'}},{{'0'},{'0','0','255'},{'solid'}},{{'0'}, ... {'0','0','255'},{'solid'}}}; bnd.ind = [1,2,1,3,1,4,5,2,2]; appl.bnd = bnd; clear equ equ.shape = {[1]}; equ.gporder = {{1}}; equ.cporder = {{1}}; equ.init = {{'0'}}; equ.usage = {1}; equ.Stype = {'userdef'}; equ.S = {'S_alpha'}; equ.thetas = {'0.25'}; equ.chip = {'1'}; equ.chif = {'1'}; equ.Ks = {'K_s'}; equ.kaps = {'1'}; equ.Ktype = {'conductivity'}; equ.Ar = {{'1','0';'0','1'}}; equ.Aron = {'0'}; equ.rhof = {'rho_f'}; equ.eta = {'0.001'}; equ.Qs = {'-alpha_b*(er_t_smaxi+ez_t_smaxi)'}; equ.CSs = {'1'}; equ.CKs = {'1'}; equ.CQs = {'1'}; equ.style = {{{'0'},{'193','193','193'}}}; equ.ind = [1,1]; appl.equ = equ; appl.var = {'tscale','1e-5', ... 'g','9.81', ... 'D','0'}; xfem.appl{1} = appl; clear appl appl.mode.class = 'SmeAxialSolid'; appl.mode.type = 'axi'; appl.dim = {'uor2','w2','p2'}; appl.sdim = {'r','phi','z'}; appl.name = 'smaxi'; appl.module = 'SME'; appl.shape = {'shlag(2,''uor2'')','shlag(2,''w2'')','shlag(1,''p2'')'}; appl.gporder = 4; appl.cporder = 2; appl.sshape = 2; appl.border = 'off'; appl.assignsuffix = '_smaxi'; appl.assign = {'Ephi','Ephi_smaxi','Er','Er_smaxi','Ez','Ez_smaxi','Fr','Fr_smaxi', ... 'FrAmp','FrAmp_smaxi','FrPh','FrPh_smaxi','Frg','Frg_smaxi','Fz','Fz_smaxi', ... 'FzAmp','FzAmp_smaxi','FzPh','FzPh_smaxi','Fzg','Fzg_smaxi','Grz','Grz_smaxi', ... 'Tar','Tar_smaxi','Taz','Taz_smaxi','alphaphi','alphaphi_smaxi','alphar', ... 'alphar_smaxi','alphaz','alphaz_smaxi','ephi','ephi_smaxi','ephii', ... 'ephii_smaxi','er','er_smaxi','eri','eri_smaxi','erz','erz_smaxi','ez', ... 'ez_smaxi','ezi','ezi_smaxi','nuphiz','nuphiz_smaxi','nurphi','nurphi_smaxi', ... 'nurz','nurz_smaxi','sphi','sphi_smaxi','sphii','sphii_smaxi','sr','sr_smaxi', ... 'sri','sri_smaxi','srz','srz_smaxi','sz','sz_smaxi','szi','szi_smaxi', ... 'uaxir','uaxir_smaxi','uaxiz','uaxiz_smaxi'}; clear prop prop.elemdefault='Lag2'; prop.analysis='time'; prop.largedef='off'; prop.eigtype='freq'; prop.createframe='off'; prop.deformframe='ref'; prop.frame='ref'; clear weakconstr weakconstr.value = 'off'; weakconstr.dim = {'lm2','lm3'}; prop.weakconstr = weakconstr; prop.constrtype='ideal'; appl.prop = prop; clear pnt pnt.Fr = {'0'}; pnt.FrPh = {'0'}; pnt.Fz = {'0'}; pnt.FzPh = {'0'}; pnt.loadcoord = {'global'}; pnt.constrcoord = {'global'}; pnt.constrtype = {'standard'}; pnt.H = {{'0','0';'0','0'}}; pnt.R = {{'0';'0'}}; pnt.Rr = {'0'}; pnt.Hr = {'0'}; pnt.Rz = {'0'}; pnt.Hz = {'0'}; pnt.style = {{{'0'},{'0','0','0'}}}; pnt.ind = [1,1,1,1,1,1,1,1]; appl.pnt = pnt; clear bnd bnd.name = {'','',''}; bnd.Fr = {'0','0','0'}; bnd.FrPh = {'0','0','0'}; bnd.Fz = {'0','0','0'}; bnd.FzPh = {'0','0','0'}; bnd.P = {'0','0','0'}; bnd.PPh = {'0','0','0'}; bnd.loadcond = {'distr_force','distr_force','distr_force'}; bnd.loadcoord = {'global','global','global'}; bnd.constrcoord = {'global','global','global'}; bnd.constrtype = {'standard','standard','standard'}; bnd.constrcond = {'free','displacement','displacement'}; bnd.velr = {'0','0','0'}; bnd.Hvelr = {'0','0','0'}; bnd.accr = {'0','0','0'}; bnd.Haccr = {'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'}}; bnd.R = {{'0';'0'},{'0';'0'},{'0';'0'}}; bnd.Rr = {'0','0','0'}; bnd.Hr = {'0','1','1'}; bnd.Rz = {'0','0','0'}; bnd.Hz = {'0','0','1'}; bnd.offset = {'0','0','0'}; bnd.pn = {'E_smaxi/h*min(1e-3*5^auglagiter,1)','E_smaxi/h*min(1e-3*5^auglagiter,1)', ... 'E_smaxi/h*min(1e-3*5^auglagiter,1)'}; bnd.pt = {'E_smaxi/h*min(1e-3*5^auglagiter,1)','E_smaxi/h*min(1e-3*5^auglagiter,1)', ... 'E_smaxi/h*min(1e-3*5^auglagiter,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 = {'1e-3','1e-3','1e-3'}; bnd.contact_oldi = {'0','0','0'}; bnd.Ttri = {'1e-3','1e-3','1e-3'}; bnd.rm_oldi = {'0','0','0'}; bnd.Ttzi = {'1e-3','1e-3','1e-3'}; bnd.zm_oldi = {'0','0','0'}; bnd.style = {{{'0'},{'0','0','0'},{'solid'}},{{'0'},{'0','0','255'},{'solid'}}, ... {{'0'},{'255','0','255'},{'solid'}}}; bnd.ind = [2,3,2,1,2,1,1,3,3]; appl.bnd = bnd; clear equ equ.shape = {[1;2;3]}; equ.gporder = {{1;1;1}}; equ.cporder = {{1;1;1}}; equ.init = {{'0';'0';'0'}}; equ.usage = {1}; equ.E = {'E'}; equ.rho = {'rho_s'}; equ.dampingtype = {'nodamping'}; equ.alphadM = {'1'}; equ.betadK = {'0.001'}; equ.eta = {'0'}; equ.dinit = {{'0';'0';'0'}}; equ.materialcoord = {'global'}; equ.materialmodel = {'iso'}; equ.mixedform = {'1'}; equ.hypertype = {'neo_hookean'}; equ.mu = {'8e5'}; equ.C10 = {'2e5'}; equ.C01 = {'2e5'}; equ.kappa = {'1e10'}; equ.hardeningmodel = {'iso'}; equ.yieldtype = {'mises'}; equ.isodata = {'tangent'}; equ.ETkin = {'2.0e10'}; equ.ETiso = {'2.0e10'}; equ.Sys = {'2.0e8'}; equ.Syfunc = {'mises_smaxi'}; equ.Syfunc_kin = {'misesKin_smaxi'}; equ.Shard = {'2.0e10[Pa]/(1-2.0e10[Pa]/2.0e11[Pa])*epe_smaxi'}; equ.ini_stress = {'1'}; equ.ini_strain = {'0'}; equ.srzi = {'0'}; equ.erzi = {'0'}; equ.sri = {'0'}; equ.eri = {'0'}; equ.sphii = {'0'}; equ.ephii = {'0'}; equ.szi = {'0'}; equ.ezi = {'0'}; equ.alpha = {'0'}; equ.nu = {'nu'}; equ.Grz = {'7.52e10'}; equ.nurphi = {'0.33'}; equ.alphar = {'1.2e-5'}; equ.Er = {'2.0e11'}; equ.nuphiz = {'0.33'}; equ.alphaphi = {'1.2e-5'}; equ.Ephi = {'2.0e11'}; equ.nurz = {'0.33'}; equ.alphaz = {'1.2e-5'}; equ.Ez = {'2.0e11'}; 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';'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';'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','2.0e11/((1+0.33)*2)'}}; equ.alphavector = {{'1.2e-5';'1.2e-5';'1.2e-5';'0'}}; equ.Fr = {'F_f2sr'}; equ.FrPh = {'0'}; equ.Fz = {'F_f2sz'}; equ.FzPh = {'0'}; equ.loadcoord = {'global'}; equ.Tflag = {'0'}; equ.Temp = {'0'}; equ.Tempref = {'0'}; equ.constrcoord = {'global'}; equ.constrtype = {'standard'}; equ.constrcond = {'displacement'}; equ.velr = {'0'}; equ.Hvelr = {'0'}; equ.accr = {'0'}; equ.Haccr = {'0'}; equ.velz = {'0'}; equ.Hvelz = {'0'}; equ.accz = {'0'}; equ.Haccz = {'0'}; equ.H = {{'0','0';'0','0'}}; equ.R = {{'0';'0'}}; equ.Rr = {'0'}; equ.Hr = {'0'}; equ.Rz = {'0'}; equ.Hz = {'0'}; equ.dr = {'dr_guess_smaxi'}; equ.Lr = {'dr_smaxi'}; equ.dz = {'dz_guess_smaxi'}; equ.Lz = {'dz_smaxi'}; equ.R0 = {'R0_guess_smaxi'}; equ.X0 = {'0'}; equ.Y0 = {'0'}; equ.Z0 = {'Z0_guess_smaxi'}; equ.userPML = {{'r';'z'}}; equ.PMLtype = {'none'}; equ.coordOnr = {'0'}; equ.coordOnz = {'0'}; equ.r0 = {'0'}; equ.z0 = {'0'}; equ.coord = {'0'}; equ.rOn = {'0'}; equ.n = {'1'}; equ.style = {{{'0'},{'0','0','255'}}}; equ.ind = [1,1]; appl.equ = equ; appl.var = {'t_old_ini','-1'}; xfem.appl{2} = appl; xfem.geom = flbinary('fem1','geom','surficial_caprock_model.mph'); xfem.sdim = {'r','z'}; xfem.frame = {'ref'}; xfem.shape = {'shlag(2,''p'')','shlag(2,''uor2'')','shlag(2,''w2'')','shlag(1,''p2'')'}; xfem.gporder = 4; xfem.cporder = 2; xfem.sshape = 2; xfem.simplify = 'off'; xfem.border = 1; xfem.outform = 'general'; xfem.form = 'general'; clear units; units.basesystem = 'SI'; xfem.units = units; clear equ equ.shape = {[1;2;3;4],[1;2;3;4]}; equ.gporder = {{1;1;1;1},{1;1;1;1}}; equ.cporder = {{1;1;1;1},{1;1;1;1}}; equ.init = {{'(((0)-D_esdl)*rhof_esdl*g_esdl)';'(0)';'(0)';'(0)'},{'(((0)-D_esdl)*rhof_esdl*g_esdl)'; ... '(0)';'(0)';'(0)'}}; equ.dinit = {{'0';'(0)';'(0)';'(0)'},{'0';'(0)';'(0)';'(0)'}}; equ.weak = {{'0';'-(er_smaxi_test*sr_smaxi+ephi_smaxi_test*sphi_smaxi+ez_smaxi_test*sz_smaxi+2*erz_smaxi_test*srz_smaxi)*r+Frg_smaxi*uor2_test*r^2'; ... 'Fzg_smaxi*w2_test*r';'r*p2_test*(p2/1/K_smaxi+(+(1)*(er_smaxi-0*alpha_smaxi*(Temp_smaxi-Tempref_smaxi)/1-0*eri_smaxi/1)+(1)*(ephi_smaxi-0*alpha_smaxi*(Temp_smaxi-Tempref_smaxi)/1-0*ephii_smaxi/1)+(1)*(ez_smaxi-0*alpha_smaxi*(Temp_smaxi-Tempref_smaxi)/1-0*ezi_smaxi/1)+(0)*(erz_smaxi-0/1-0*erzi_smaxi/1))-1*-(sri_smaxi+sphii_smaxi+szi_smaxi)/3/1/K_smaxi)'}, ... {'0';'-(er_smaxi_test*sr_smaxi+ephi_smaxi_test*sphi_smaxi+ez_smaxi_test*sz_smaxi+2*erz_smaxi_test*srz_smaxi)*r+Frg_smaxi*uor2_test*r^2'; ... 'Fzg_smaxi*w2_test*r';'r*p2_test*(p2/1/K_smaxi+(+(1)*(er_smaxi-0*alpha_smaxi*(Temp_smaxi-Tempref_smaxi)/1-0*eri_smaxi/1)+(1)*(ephi_smaxi-0*alpha_smaxi*(Temp_smaxi-Tempref_smaxi)/1-0*ephii_smaxi/1)+(1)*(ez_smaxi-0*alpha_smaxi*(Temp_smaxi-Tempref_smaxi)/1-0*ezi_smaxi/1)+(0)*(erz_smaxi-0/1-0*erzi_smaxi/1))-1*-(sri_smaxi+sphii_smaxi+szi_smaxi)/3/1/K_smaxi)'}}; equ.dweak = {{'0';'rho_smaxi*(uor2_test*uor2tt*r^3+w2_test*w2tt*r)';'0';'0'}, ... {'0';'rho_smaxi*(uor2_test*uor2tt*r^3+w2_test*w2tt*r)';'0';'0'}}; equ.constr = {{'0';'0-((((0)*(1)+(0)*(0)))*r*uor2+(((0)*(0)+(0)*(1)))*w2)'; ... '0-((((0)*(1)+(0)*(0)))*r*uor2+(((0)*(0)+(0)*(1)))*w2)';'0'},{'0';'0-((((0)*(1)+(0)*(0)))*r*uor2+(((0)*(0)+(0)*(1)))*w2)'; ... '0-((((0)*(1)+(0)*(0)))*r*uor2+(((0)*(0)+(0)*(1)))*w2)';'0'}}; equ.constrf = {{'0';'0';'0';'0'},{'0';'0';'0';'0'}}; equ.c = {{{'-diff(-(r*K_esdl/rhof_esdl/g_esdl)*pr+(-r*K_esdl*diff(D_esdl,r)),pr)', ... '-diff(-(r*K_esdl/rhof_esdl/g_esdl)*pr+(-r*K_esdl*diff(D_esdl,r)),pz)'; ... '-diff(-(r*K_esdl/rhof_esdl/g_esdl)*pz+(-r*K_esdl*diff(D_esdl,z)),pr)', ... '-diff(-(r*K_esdl/rhof_esdl/g_esdl)*pz+(-r*K_esdl*diff(D_esdl,z)),pz)'}, ... {'-diff(-(r*K_esdl/rhof_esdl/g_esdl)*pr+(-r*K_esdl*diff(D_esdl,r)),uor2r)', ... '-diff(-(r*K_esdl/rhof_esdl/g_esdl)*pr+(-r*K_esdl*diff(D_esdl,r)),uor2z)'; ... '-diff(-(r*K_esdl/rhof_esdl/g_esdl)*pz+(-r*K_esdl*diff(D_esdl,z)),uor2r)', ... '-diff(-(r*K_esdl/rhof_esdl/g_esdl)*pz+(-r*K_esdl*diff(D_esdl,z)),uor2z)'}, ... {'-diff(-(r*K_esdl/rhof_esdl/g_esdl)*pr+(-r*K_esdl*diff(D_esdl,r)),w2r)', ... '-diff(-(r*K_esdl/rhof_esdl/g_esdl)*pr+(-r*K_esdl*diff(D_esdl,r)),w2z)'; ... '-diff(-(r*K_esdl/rhof_esdl/g_esdl)*pz+(-r*K_esdl*diff(D_esdl,z)),w2r)', ... '-diff(-(r*K_esdl/rhof_esdl/g_esdl)*pz+(-r*K_esdl*diff(D_esdl,z)),w2z)'}, ... {'-diff(-(r*K_esdl/rhof_esdl/g_esdl)*pr+(-r*K_esdl*diff(D_esdl,r)),p2r)', ... '-diff(-(r*K_esdl/rhof_esdl/g_esdl)*pr+(-r*K_esdl*diff(D_esdl,r)),p2z)'; ... '-diff(-(r*K_esdl/rhof_esdl/g_esdl)*pz+(-r*K_esdl*diff(D_esdl,z)),p2r)', ... '-diff(-(r*K_esdl/rhof_esdl/g_esdl)*pz+(-r*K_esdl*diff(D_esdl,z)),p2z)'}; ... {'-0'},{'-0'},{'-0'},{'-0'};{'-0'},{'-0'},{'-0'},{'-0'};{'-0'},{'-0'},{'-0'}, ... {'-0'}},{{'-diff(-(r*K_esdl/rhof_esdl/g_esdl)*pr+(-r*K_esdl*diff(D_esdl,r)),pr)', ... '-diff(-(r*K_esdl/rhof_esdl/g_esdl)*pr+(-r*K_esdl*diff(D_esdl,r)),pz)'; ... '-diff(-(r*K_esdl/rhof_esdl/g_esdl)*pz+(-r*K_esdl*diff(D_esdl,z)),pr)', ... '-diff(-(r*K_esdl/rhof_esdl/g_esdl)*pz+(-r*K_esdl*diff(D_esdl,z)),pz)'}, ... {'-diff(-(r*K_esdl/rhof_esdl/g_esdl)*pr+(-r*K_esdl*diff(D_esdl,r)),uor2r)', ... '-diff(-(r*K_esdl/rhof_esdl/g_esdl)*pr+(-r*K_esdl*diff(D_esdl,r)),uor2z)'; ... '-diff(-(r*K_esdl/rhof_esdl/g_esdl)*pz+(-r*K_esdl*diff(D_esdl,z)),uor2r)', ... '-diff(-(r*K_esdl/rhof_esdl/g_esdl)*pz+(-r*K_esdl*diff(D_esdl,z)),uor2z)'}, ... {'-diff(-(r*K_esdl/rhof_esdl/g_esdl)*pr+(-r*K_esdl*diff(D_esdl,r)),w2r)', ... '-diff(-(r*K_esdl/rhof_esdl/g_esdl)*pr+(-r*K_esdl*diff(D_esdl,r)),w2z)'; ... '-diff(-(r*K_esdl/rhof_esdl/g_esdl)*pz+(-r*K_esdl*diff(D_esdl,z)),w2r)', ... '-diff(-(r*K_esdl/rhof_esdl/g_esdl)*pz+(-r*K_esdl*diff(D_esdl,z)),w2z)'}, ... {'-diff(-(r*K_esdl/rhof_esdl/g_esdl)*pr+(-r*K_esdl*diff(D_esdl,r)),p2r)', ... '-diff(-(r*K_esdl/rhof_esdl/g_esdl)*pr+(-r*K_esdl*diff(D_esdl,r)),p2z)'; ... '-diff(-(r*K_esdl/rhof_esdl/g_esdl)*pz+(-r*K_esdl*diff(D_esdl,z)),p2r)', ... '-diff(-(r*K_esdl/rhof_esdl/g_esdl)*pz+(-r*K_esdl*diff(D_esdl,z)),p2z)'}; ... {'-0'},{'-0'},{'-0'},{'-0'};{'-0'},{'-0'},{'-0'},{'-0'};{'-0'},{'-0'},{'-0'}, ... {'-0'}}}; equ.a = {{'-diff(+(r*Qs_esdl),p)','-diff(+(r*Qs_esdl),uor2)','-diff(+(r*Qs_esdl),w2)', ... '-diff(+(r*Qs_esdl),p2)';'-0','-0','-0','-0';'-0','-0','-0','-0';'-0','-0', ... '-0','-0'},{'-diff(+(r*Qs_esdl),p)','-diff(+(r*Qs_esdl),uor2)','-diff(+(r*Qs_esdl),w2)', ... '-diff(+(r*Qs_esdl),p2)';'-0','-0','-0','-0';'-0','-0','-0','-0';'-0','-0', ... '-0','-0'}}; equ.f = {{'+(r*Qs_esdl)';'0';'0';'0'},{'+(r*Qs_esdl)';'0';'0';'0'}}; equ.ea = {{'0','0','0','0';'0','0','0','0';'0','0','0','0';'0','0','0', ... '0'},{'0','0','0','0';'0','0','0','0';'0','0','0','0';'0','0','0','0'}}; equ.da = {{'r*S_esdl','0','0','0';'0','0','0','0';'0','0','0','0';'0','0', ... '0','0'},{'r*S_esdl','0','0','0';'0','0','0','0';'0','0','0','0';'0','0', ... '0','0'}}; equ.al = {{{'-diff(-(r*K_esdl/rhof_esdl/g_esdl)*pr+(-r*K_esdl*diff(D_esdl,r)),p)'; ... '-diff(-(r*K_esdl/rhof_esdl/g_esdl)*pz+(-r*K_esdl*diff(D_esdl,z)),p)'}, ... {'-diff(-(r*K_esdl/rhof_esdl/g_esdl)*pr+(-r*K_esdl*diff(D_esdl,r)),uor2)'; ... '-diff(-(r*K_esdl/rhof_esdl/g_esdl)*pz+(-r*K_esdl*diff(D_esdl,z)),uor2)'}, ... {'-diff(-(r*K_esdl/rhof_esdl/g_esdl)*pr+(-r*K_esdl*diff(D_esdl,r)),w2)'; ... '-diff(-(r*K_esdl/rhof_esdl/g_esdl)*pz+(-r*K_esdl*diff(D_esdl,z)),w2)'}, ... {'-diff(-(r*K_esdl/rhof_esdl/g_esdl)*pr+(-r*K_esdl*diff(D_esdl,r)),p2)'; ... '-diff(-(r*K_esdl/rhof_esdl/g_esdl)*pz+(-r*K_esdl*diff(D_esdl,z)),p2)'}; ... {'-0';'-0'},{'-0';'-0'},{'-0';'-0'},{'-0';'-0'};{'-0';'-0'},{'-0';'-0'}, ... {'-0';'-0'},{'-0';'-0'};{'-0';'-0'},{'-0';'-0'},{'-0';'-0'},{'-0';'-0'}}, ... {{'-diff(-(r*K_esdl/rhof_esdl/g_esdl)*pr+(-r*K_esdl*diff(D_esdl,r)),p)'; ... '-diff(-(r*K_esdl/rhof_esdl/g_esdl)*pz+(-r*K_esdl*diff(D_esdl,z)),p)'}, ... {'-diff(-(r*K_esdl/rhof_esdl/g_esdl)*pr+(-r*K_esdl*diff(D_esdl,r)),uor2)'; ... '-diff(-(r*K_esdl/rhof_esdl/g_esdl)*pz+(-r*K_esdl*diff(D_esdl,z)),uor2)'}, ... {'-diff(-(r*K_esdl/rhof_esdl/g_esdl)*pr+(-r*K_esdl*diff(D_esdl,r)),w2)'; ... '-diff(-(r*K_esdl/rhof_esdl/g_esdl)*pz+(-r*K_esdl*diff(D_esdl,z)),w2)'}, ... {'-diff(-(r*K_esdl/rhof_esdl/g_esdl)*pr+(-r*K_esdl*diff(D_esdl,r)),p2)'; ... '-diff(-(r*K_esdl/rhof_esdl/g_esdl)*pz+(-r*K_esdl*diff(D_esdl,z)),p2)'}; ... {'-0';'-0'},{'-0';'-0'},{'-0';'-0'},{'-0';'-0'};{'-0';'-0'},{'-0';'-0'}, ... {'-0';'-0'},{'-0';'-0'};{'-0';'-0'},{'-0';'-0'},{'-0';'-0'},{'-0';'-0'}}}; equ.be = {{{'-diff(+(r*Qs_esdl),pr)';'-diff(+(r*Qs_esdl),pz)'},{'-diff(+(r*Qs_esdl),uor2r)'; ... '-diff(+(r*Qs_esdl),uor2z)'},{'-diff(+(r*Qs_esdl),w2r)';'-diff(+(r*Qs_esdl),w2z)'}, ... {'-diff(+(r*Qs_esdl),p2r)';'-diff(+(r*Qs_esdl),p2z)'};{'-0';'-0'},{'-0'; ... '-0'},{'-0';'-0'},{'-0';'-0'};{'-0';'-0'},{'-0';'-0'},{'-0';'-0'},{'-0'; ... '-0'};{'-0';'-0'},{'-0';'-0'},{'-0';'-0'},{'-0';'-0'}},{{'-diff(+(r*Qs_esdl),pr)'; ... '-diff(+(r*Qs_esdl),pz)'},{'-diff(+(r*Qs_esdl),uor2r)';'-diff(+(r*Qs_esdl),uor2z)'}, ... {'-diff(+(r*Qs_esdl),w2r)';'-diff(+(r*Qs_esdl),w2z)'},{'-diff(+(r*Qs_esdl),p2r)'; ... '-diff(+(r*Qs_esdl),p2z)'};{'-0';'-0'},{'-0';'-0'},{'-0';'-0'},{'-0';'-0'}; ... {'-0';'-0'},{'-0';'-0'},{'-0';'-0'},{'-0';'-0'};{'-0';'-0'},{'-0';'-0'}, ... {'-0';'-0'},{'-0';'-0'}}}; equ.ga = {{{'-(r*K_esdl/rhof_esdl/g_esdl)*pr+(-r*K_esdl*diff(D_esdl,r))'; ... '-(r*K_esdl/rhof_esdl/g_esdl)*pz+(-r*K_esdl*diff(D_esdl,z))'};{'0';'0'}; ... {'0';'0'};{'0';'0'}},{{'-(r*K_esdl/rhof_esdl/g_esdl)*pr+(-r*K_esdl*diff(D_esdl,r))'; ... '-(r*K_esdl/rhof_esdl/g_esdl)*pz+(-r*K_esdl*diff(D_esdl,z))'};{'0';'0'}; ... {'0';'0'};{'0';'0'}}}; equ.sshape = {[1],[1]}; equ.sshapedim = {{1;1;1;1},{1;1;1;1}}; equ.ind = [1,2]; equ.dim = {'p','uor2','w2','p2'}; equ.var = {'S_esdl',{'(S_alpha)/rhof_esdl/g_esdl*CSs_esdl','(S_alpha)/rhof_esdl/g_esdl*CSs_esdl'}, ... 'Qs_esdl',{'(-alpha_b*(er_t_smaxi+ez_t_smaxi))*CQs_esdl','(-alpha_b*(er_t_smaxi+ez_t_smaxi))*CQs_esdl'}, ... 'K_esdl',{'Ks_esdl*CKs_esdl','Ks_esdl*CKs_esdl'}, ... 'Krr_esdl',{'K_esdl','K_esdl'}, ... 'Krz_esdl',{'0','0'}, ... 'Kzr_esdl',{'0','0'}, ... 'Kzz_esdl',{'K_esdl','K_esdl'}, ... 'kap_esdl',{'kaps_esdl*CKs_esdl','kaps_esdl*CKs_esdl'}, ... 'kaprr_esdl',{'kap_esdl','kap_esdl'}, ... 'kaprz_esdl',{'0','0'}, ... 'kapzr_esdl',{'0','0'}, ... 'kapzz_esdl',{'kap_esdl','kap_esdl'}, ... 'H_esdl',{'p/rhof_esdl/g_esdl+D_esdl','p/rhof_esdl/g_esdl+D_esdl'}, ... 'gradP_esdl',{'sqrt(+pr^2+pz^2)','sqrt(+pr^2+pz^2)'}, ... 'Hr_esdl',{'diff(H_esdl,r)','diff(H_esdl,r)'}, ... 'Hz_esdl',{'diff(H_esdl,z)','diff(H_esdl,z)'}, ... 'gradH_esdl',{'sqrt(+Hr_esdl^2+Hz_esdl^2)','sqrt(+Hr_esdl^2+Hz_esdl^2)'}, ... 'u_esdl',{'1/rhof_esdl/g_esdl*(-Krr_esdl*(pr+diff(rhof_esdl*g_esdl*D_esdl,r))-Krz_esdl*(pz+diff(rhof_esdl*g_esdl*D_esdl,z)))', ... '1/rhof_esdl/g_esdl*(-Krr_esdl*(pr+diff(rhof_esdl*g_esdl*D_esdl,r))-Krz_esdl*(pz+diff(rhof_esdl*g_esdl*D_esdl,z)))'}, ... 'v_esdl',{'1/rhof_esdl/g_esdl*(-Kzr_esdl*(pr+diff(rhof_esdl*g_esdl*D_esdl,r))-Kzz_esdl*(pz+diff(rhof_esdl*g_esdl*D_esdl,z)))', ... '1/rhof_esdl/g_esdl*(-Kzr_esdl*(pr+diff(rhof_esdl*g_esdl*D_esdl,r))-Kzz_esdl*(pz+diff(rhof_esdl*g_esdl*D_esdl,z)))'}, ... 'U_esdl',{'sqrt(+(u_esdl)^2+(v_esdl)^2)','sqrt(+(u_esdl)^2+(v_esdl)^2)'}, ... 'Frg_smaxi',{'Fr_smaxi','Fr_smaxi'}, ... 'Fzg_smaxi',{'Fz_smaxi','Fz_smaxi'}, ... 'disp_smaxi',{'sqrt(real(uaxi_smaxi)^2+real(w2)^2)','sqrt(real(uaxi_smaxi)^2+real(w2)^2)'}, ... 'uaxi_smaxi',{'uor2*r','uor2*r'}, ... 'uaxir_smaxi',{'uor2r*r+uor2','uor2r*r+uor2'}, ... 'uaxiz_smaxi',{'uor2z*r','uor2z*r'}, ... 'uaxi_t_smaxi',{'uor2t*r','uor2t*r'}, ... 'sr_smaxi',{'(4/3*G_smaxi)*(1*er_smaxi-0*eri_smaxi-0*alpha_smaxi*(Temp_smaxi-Tempref_smaxi))+1*sri_smaxi+-(1*sri_smaxi+1*sphii_smaxi+1*szi_smaxi)/3+(-2/3*G_smaxi)*(1*ephi_smaxi-0*ephii_smaxi-0*alpha_smaxi*(Temp_smaxi-Tempref_smaxi))+(-2/3*G_smaxi)*(1*ez_smaxi-0*ezi_smaxi-0*alpha_smaxi*(Temp_smaxi-Tempref_smaxi))+0*(1*erz_smaxi-0*erzi_smaxi-0*0*(Temp_smaxi-Tempref_smaxi))-1*p2', ... '(4/3*G_smaxi)*(1*er_smaxi-0*eri_smaxi-0*alpha_smaxi*(Temp_smaxi-Tempref_smaxi))+1*sri_smaxi+-(1*sri_smaxi+1*sphii_smaxi+1*szi_smaxi)/3+(-2/3*G_smaxi)*(1*ephi_smaxi-0*ephii_smaxi-0*alpha_smaxi*(Temp_smaxi-Tempref_smaxi))+(-2/3*G_smaxi)*(1*ez_smaxi-0*ezi_smaxi-0*alpha_smaxi*(Temp_smaxi-Tempref_smaxi))+0*(1*erz_smaxi-0*erzi_smaxi-0*0*(Temp_smaxi-Tempref_smaxi))-1*p2'}, ... 'sz_smaxi',{'(-2/3*G_smaxi)*(1*er_smaxi-0*eri_smaxi-0*alpha_smaxi*(Temp_smaxi-Tempref_smaxi))+1*szi_smaxi+-(1*sri_smaxi+1*sphii_smaxi+1*szi_smaxi)/3+(-2/3*G_smaxi)*(1*ephi_smaxi-0*ephii_smaxi-0*alpha_smaxi*(Temp_smaxi-Tempref_smaxi))+(4/3*G_smaxi)*(1*ez_smaxi-0*ezi_smaxi-0*alpha_smaxi*(Temp_smaxi-Tempref_smaxi))+0*(1*erz_smaxi-0*erzi_smaxi-0*0*(Temp_smaxi-Tempref_smaxi))-1*p2', ... '(-2/3*G_smaxi)*(1*er_smaxi-0*eri_smaxi-0*alpha_smaxi*(Temp_smaxi-Tempref_smaxi))+1*szi_smaxi+-(1*sri_smaxi+1*sphii_smaxi+1*szi_smaxi)/3+(-2/3*G_smaxi)*(1*ephi_smaxi-0*ephii_smaxi-0*alpha_smaxi*(Temp_smaxi-Tempref_smaxi))+(4/3*G_smaxi)*(1*ez_smaxi-0*ezi_smaxi-0*alpha_smaxi*(Temp_smaxi-Tempref_smaxi))+0*(1*erz_smaxi-0*erzi_smaxi-0*0*(Temp_smaxi-Tempref_smaxi))-1*p2'}, ... 'sphi_smaxi',{'(-2/3*G_smaxi)*(1*er_smaxi-0*eri_smaxi-0*alpha_smaxi*(Temp_smaxi-Tempref_smaxi))+1*sphii_smaxi+-(1*sri_smaxi+1*sphii_smaxi+1*szi_smaxi)/3+(4/3*G_smaxi)*(1*ephi_smaxi-0*ephii_smaxi-0*alpha_smaxi*(Temp_smaxi-Tempref_smaxi))+(-2/3*G_smaxi)*(1*ez_smaxi-0*ezi_smaxi-0*alpha_smaxi*(Temp_smaxi-Tempref_smaxi))+0*(1*erz_smaxi-0*erzi_smaxi-0*0*(Temp_smaxi-Tempref_smaxi))-1*p2', ... '(-2/3*G_smaxi)*(1*er_smaxi-0*eri_smaxi-0*alpha_smaxi*(Temp_smaxi-Tempref_smaxi))+1*sphii_smaxi+-(1*sri_smaxi+1*sphii_smaxi+1*szi_smaxi)/3+(4/3*G_smaxi)*(1*ephi_smaxi-0*ephii_smaxi-0*alpha_smaxi*(Temp_smaxi-Tempref_smaxi))+(-2/3*G_smaxi)*(1*ez_smaxi-0*ezi_smaxi-0*alpha_smaxi*(Temp_smaxi-Tempref_smaxi))+0*(1*erz_smaxi-0*erzi_smaxi-0*0*(Temp_smaxi-Tempref_smaxi))-1*p2'}, ... 'srz_smaxi',{'0*(1*er_smaxi-0*eri_smaxi-0*alpha_smaxi*(Temp_smaxi-Tempref_smaxi))+1*srzi_smaxi+0*(1*ephi_smaxi-0*ephii_smaxi-0*alpha_smaxi*(Temp_smaxi-Tempref_smaxi))+0*(1*ez_smaxi-0*ezi_smaxi-0*alpha_smaxi*(Temp_smaxi-Tempref_smaxi))+(2*G_smaxi)*(1*erz_smaxi-0*erzi_smaxi-0*0*(Temp_smaxi-Tempref_smaxi))-0*p2', ... '0*(1*er_smaxi-0*eri_smaxi-0*alpha_smaxi*(Temp_smaxi-Tempref_smaxi))+1*srzi_smaxi+0*(1*ephi_smaxi-0*ephii_smaxi-0*alpha_smaxi*(Temp_smaxi-Tempref_smaxi))+0*(1*ez_smaxi-0*ezi_smaxi-0*alpha_smaxi*(Temp_smaxi-Tempref_smaxi))+(2*G_smaxi)*(1*erz_smaxi-0*erzi_smaxi-0*0*(Temp_smaxi-Tempref_smaxi))-0*p2'}, ... 'er_smaxi',{'((uor2r*r+uor2))','((uor2r*r+uor2))'}, ... 'ez_smaxi',{'w2z','w2z'}, ... 'ephi_smaxi',{'uor2','uor2'}, ... 'erz_smaxi',{'((uor2z*r)+w2r)/2','((uor2z*r)+w2r)/2'}, ... 'sr_t_smaxi',{'(4/3*G_smaxi)*(er_t_smaxi-0-0*(0-0))+0+-(0+0+0)/3+(-2/3*G_smaxi)*(ephi_t_smaxi-0-0*(0-0))+(-2/3*G_smaxi)*(ez_t_smaxi-0-0*(0-0))+0*(erz_t_smaxi-0-0*(0-0))-1*p2t', ... '(4/3*G_smaxi)*(er_t_smaxi-0-0*(0-0))+0+-(0+0+0)/3+(-2/3*G_smaxi)*(ephi_t_smaxi-0-0*(0-0))+(-2/3*G_smaxi)*(ez_t_smaxi-0-0*(0-0))+0*(erz_t_smaxi-0-0*(0-0))-1*p2t'}, ... 'sz_t_smaxi',{'(-2/3*G_smaxi)*(er_t_smaxi-0-0*(0-0))+0+-(0+0+0)/3+(-2/3*G_smaxi)*(ephi_t_smaxi-0-0*(0-0))+(4/3*G_smaxi)*(ez_t_smaxi-0-0*(0-0))+0*(erz_t_smaxi-0-0*(0-0))-1*p2t', ... '(-2/3*G_smaxi)*(er_t_smaxi-0-0*(0-0))+0+-(0+0+0)/3+(-2/3*G_smaxi)*(ephi_t_smaxi-0-0*(0-0))+(4/3*G_smaxi)*(ez_t_smaxi-0-0*(0-0))+0*(erz_t_smaxi-0-0*(0-0))-1*p2t'}, ... 'sphi_t_smaxi',{'(-2/3*G_smaxi)*(er_t_smaxi-0-0*(0-0))+0+-(0+0+0)/3+(4/3*G_smaxi)*(ephi_t_smaxi-0-0*(0-0))+(-2/3*G_smaxi)*(ez_t_smaxi-0-0*(0-0))+0*(erz_t_smaxi-0-0*(0-0))-1*p2t', ... '(-2/3*G_smaxi)*(er_t_smaxi-0-0*(0-0))+0+-(0+0+0)/3+(4/3*G_smaxi)*(ephi_t_smaxi-0-0*(0-0))+(-2/3*G_smaxi)*(ez_t_smaxi-0-0*(0-0))+0*(erz_t_smaxi-0-0*(0-0))-1*p2t'}, ... 'srz_t_smaxi',{'0*(er_t_smaxi-0-0*(0-0))+0+0*(ephi_t_smaxi-0-0*(0-0))+0*(ez_t_smaxi-0-0*(0-0))+(2*G_smaxi)*(erz_t_smaxi-0-0*(0-0))-0*p2t', ... '0*(er_t_smaxi-0-0*(0-0))+0+0*(ephi_t_smaxi-0-0*(0-0))+0*(ez_t_smaxi-0-0*(0-0))+(2*G_smaxi)*(erz_t_smaxi-0-0*(0-0))-0*p2t'}, ... 'er_t_smaxi',{'(uor2rt*r+uor2t)','(uor2rt*r+uor2t)'}, ... 'ez_t_smaxi',{'w2zt','w2zt'}, ... 'ephi_t_smaxi',{'uor2t','uor2t'}, ... 'erz_t_smaxi',{'(uor2zt*r+w2rt)/2','(uor2zt*r+w2rt)/2'}, ... 'K_smaxi',{'E_smaxi/(3*(1-2*nu_smaxi))','E_smaxi/(3*(1-2*nu_smaxi))'}, ... 'G_smaxi',{'E_smaxi/(2*(1+nu_smaxi))','E_smaxi/(2*(1+nu_smaxi))'}, ... 'mises_smaxi',{'sqrt(sr_smaxi^2+sphi_smaxi^2+sz_smaxi^2-sr_smaxi*sphi_smaxi-sphi_smaxi*sz_smaxi-sr_smaxi*sz_smaxi+3*srz_smaxi^2)', ... 'sqrt(sr_smaxi^2+sphi_smaxi^2+sz_smaxi^2-sr_smaxi*sphi_smaxi-sphi_smaxi*sz_smaxi-sr_smaxi*sz_smaxi+3*srz_smaxi^2)'}, ... 'Ws_smaxi',{'0.5*(er_smaxi*sr_smaxi+ephi_smaxi*sphi_smaxi+ez_smaxi*sz_smaxi+2*erz_smaxi*srz_smaxi)', ... '0.5*(er_smaxi*sr_smaxi+ephi_smaxi*sphi_smaxi+ez_smaxi*sz_smaxi+2*erz_smaxi*srz_smaxi)'}, ... 'evol_smaxi',{'er_smaxi+ephi_smaxi+ez_smaxi','er_smaxi+ephi_smaxi+ez_smaxi'}, ... 'tresca_smaxi',{'max(max(abs(s1_smaxi-s2_smaxi),abs(s2_smaxi-s3_smaxi)),abs(s1_smaxi-s3_smaxi))', ... 'max(max(abs(s1_smaxi-s2_smaxi),abs(s2_smaxi-s3_smaxi)),abs(s1_smaxi-s3_smaxi))'},'thetas_esdl',{'0.25','0.25'}, ... 'chip_esdl',{'1','1'}, ... 'chif_esdl',{'1','1'}, ... 'Ks_esdl',{'K_s','K_s'}, ... 'kaps_esdl',{'1','1'}, ... 'rhof_esdl',{'rho_f','rho_f'}, ... 'eta_esdl',{'0.001','0.001'}, ... 'CSs_esdl',{'1','1'}, ... 'CKs_esdl',{'1','1'}, ... 'CQs_esdl',{'1','1'}, ... 'E_smaxi',{'E','E'}, ... 'rho_smaxi',{'rho_s','rho_s'}, ... 'alphadM_smaxi',{'1','1'}, ... 'betadK_smaxi',{'0.001','0.001'}, ... 'eta_smaxi',{'0','0'}, ... 'mixedform_smaxi',{'1','1'}, ... 'mu_smaxi',{'8e5','8e5'}, ... 'C10_smaxi',{'2e5','2e5'}, ... 'C01_smaxi',{'2e5','2e5'}, ... 'kappa_smaxi',{'1e10','1e10'}, ... 'ETkin_smaxi',{'2.0e10','2.0e10'}, ... 'ETiso_smaxi',{'2.0e10','2.0e10'}, ... 'Sys_smaxi',{'2.0e8','2.0e8'}, ... 'Syfunc_smaxi',{'mises_smaxi','mises_smaxi'}, ... 'Syfunc_kin_smaxi',{'misesKin_smaxi','misesKin_smaxi'}, ... 'Shard_smaxi',{'2.0e10[Pa]/(1-2.0e10[Pa]/2.0e11[Pa])*epe_smaxi','2.0e10[Pa]/(1-2.0e10[Pa]/2.0e11[Pa])*epe_smaxi'}, ... 'srzi_smaxi',{'0','0'}, ... 'erzi_smaxi',{'0','0'}, ... 'sri_smaxi',{'0','0'}, ... 'eri_smaxi',{'0','0'}, ... 'sphii_smaxi',{'0','0'}, ... 'ephii_smaxi',{'0','0'}, ... 'szi_smaxi',{'0','0'}, ... 'ezi_smaxi',{'0','0'}, ... 'alpha_smaxi',{'0','0'}, ... 'nu_smaxi',{'nu','nu'}, ... 'Grz_smaxi',{'7.52e10','7.52e10'}, ... 'nurphi_smaxi',{'0.33','0.33'}, ... 'alphar_smaxi',{'1.2e-5','1.2e-5'}, ... 'Er_smaxi',{'2.0e11','2.0e11'}, ... 'nuphiz_smaxi',{'0.33','0.33'}, ... 'alphaphi_smaxi',{'1.2e-5','1.2e-5'}, ... 'Ephi_smaxi',{'2.0e11','2.0e11'}, ... 'nurz_smaxi',{'0.33','0.33'}, ... 'alphaz_smaxi',{'1.2e-5','1.2e-5'}, ... 'Ez_smaxi',{'2.0e11','2.0e11'}, ... 'Fr_smaxi',{'F_f2sr','F_f2sr'}, ... 'FrPh_smaxi',{'0','0'}, ... 'Fz_smaxi',{'F_f2sz','F_f2sz'}, ... 'FzPh_smaxi',{'0','0'}, ... 'Temp_smaxi',{'0','0'}, ... 'Tempref_smaxi',{'0','0'}, ... 'velr_smaxi',{'0','0'}, ... 'accr_smaxi',{'0','0'}, ... 'velz_smaxi',{'0','0'}, ... 'accz_smaxi',{'0','0'}, ... 'Rr_smaxi',{'0','0'}, ... 'Rz_smaxi',{'0','0'}, ... 'dr_smaxi',{'dr_guess_smaxi','dr_guess_smaxi'}, ... 'Lr_smaxi',{'dr_smaxi','dr_smaxi'}, ... 'dz_smaxi',{'dz_guess_smaxi','dz_guess_smaxi'}, ... 'Lz_smaxi',{'dz_smaxi','dz_smaxi'}, ... 'R0_smaxi',{'R0_guess_smaxi','R0_guess_smaxi'}, ... 'X0_smaxi',{'0','0'}, ... 'Y0_smaxi',{'0','0'}, ... 'Z0_smaxi',{'Z0_guess_smaxi','Z0_guess_smaxi'}, ... 'r0_smaxi',{'0','0'}, ... 'z0_smaxi',{'0','0'}, ... 'n_smaxi',{'1','1'}, ... 'H11_smaxi',{'((0)*(1)+(0)*(0))','((0)*(1)+(0)*(0))'}, ... 'H12_smaxi',{'((0)*(0)+(0)*(1))','((0)*(0)+(0)*(1))'}, ... 'H21_smaxi',{'((0)*(1)+(0)*(0))','((0)*(1)+(0)*(0))'}, ... 'H22_smaxi',{'((0)*(0)+(0)*(1))','((0)*(0)+(0)*(1))'}}; equ.expr = {'S_alpha',{'0','0'}, ... 'K_s',{'1e-5','1e-7'}, ... 'nu',{'0.25','0.25'}, ... 'E',{'50[MPa]','50[MPa]'}, ... 'C_s2f',{'alpha_b*(er_t_smaxi+ez_t_smaxi)','alpha_b*(er_t_smaxi+ez_t_smaxi)'}}; equ.bnd.weak = {{'0';'0';'0';'0'}}; equ.bnd.gporder = {{1;1;1;1}}; equ.bnd.ind = [1,1]; equ.bnd.var = {}; equ.bnd.expr = {}; equ.lock = [0,0]; equ.mlock = {[0,0]}; xfem.equ = equ; clear bnd bnd.weak = {{'0';'0*uor2_test*r^2*1';'0*w2_test*r*1';'0'},{'0';'0*uor2_test*r^2*1'; ... '0*w2_test*r*1';'0'},{'0';'0';'0';'0'},{'0';'0';'0';'0'},{'0';'0*uor2_test*r^2*1'; ... '0*w2_test*r*1';'0'}}; bnd.dweak = {{'0';'0';'0';'0'},{'0';'0';'0';'0'},{'0';'0';'0';'0'},{'0';'0'; ... '0';'0'},{'0';'0';'0';'0'}}; bnd.constr = {{'0';'0-(((((1))*(1)+(0)*(0)))*r*uor2+((((1))*(0)+(0)*(1)))*w2)'; ... '0-((((0)*(1)+(0)*(0)))*r*uor2+(((0)*(0)+(0)*(1)))*w2)';'0'},{'0';'0-(((((1))*(1)+(0)*(0)))*r*uor2+((((1))*(0)+(0)*(1)))*w2)'; ... '0-((((0)*(1)+((1))*(0)))*r*uor2+(((0)*(0)+((1))*(1)))*w2)';'0'},{'0';'0'; ... '0';'0'},{'0';'0';'0';'0'},{'0';'0-((((0)*(1)+(0)*(0)))*r*uor2+(((0)*(0)+(0)*(1)))*w2)'; ... '0-((((0)*(1)+(0)*(0)))*r*uor2+(((0)*(0)+(0)*(1)))*w2)';'0'}}; bnd.constrf = {{'0';'test(-r*uor2)';'0';'0'},{'0';'test(-r*uor2)';'test(-w2)'; ... '0'},{'0';'0';'0';'0'},{'0';'0';'0';'0'},{'0';'0';'0';'0'}}; bnd.q = {{'-0','-0','-0','-0';'-0','-0','-0','-0';'-0','-0','-0','-0'; ... '-0','-0','-0','-0'},{'-0','-0','-0','-0';'-0','-0','-0','-0';'-0','-0', ... '-0','-0';'-0','-0','-0','-0'},{'-diff(+(r*N0_esdl),p)','-diff(+(r*N0_esdl),uor2)', ... '-diff(+(r*N0_esdl),w2)','-diff(+(r*N0_esdl),p2)';'-0','-0','-0','-0';'-0', ... '-0','-0','-0';'-0','-0','-0','-0'},{'-0','-0','-0','-0';'-0','-0','-0', ... '-0';'-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.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','-0';'-0','-0','-0','-0'},{'-0','-0','-0','-0';'-0','-0','-0','-0'; ... '-0','-0','-0','-0';'-0','-0','-0','-0'},{'-0','-0','-0','-0';'-0','-0', ... '-0','-0';'-0','-0','-0','-0';'-0','-0','-0','-0'},{'-diff(-(1)*p+(rhof_esdl*g_esdl*(0-D_esdl)),p)', ... '-diff(-(1)*p+(rhof_esdl*g_esdl*(0-D_esdl)),uor2)','-diff(-(1)*p+(rhof_esdl*g_esdl*(0-D_esdl)),w2)', ... '-diff(-(1)*p+(rhof_esdl*g_esdl*(0-D_esdl)),p2)';'-0','-0','-0','-0';'-0', ... '-0','-0','-0';'-0','-0','-0','-0'}}; bnd.g = {{'0';'0';'0';'0'},{'0';'0';'0';'0'},{'+(r*N0_esdl)';'0';'0';'0'}, ... {'0';'0';'0';'0'},{'0';'0';'0';'0'}}; bnd.r = {{'0';'0';'0';'0'},{'0';'0';'0';'0'},{'0';'0';'0';'0'},{'0';'0'; ... '0';'0'},{'-(1)*p+(rhof_esdl*g_esdl*(0-D_esdl))';'0';'0';'0'}}; bnd.shape = {[1;2;3;4],[1;2;3;4],[1;2;3;4],[1;2;3;4],[1;2;3;4]}; bnd.sshape = {[1],[1],[1],[1],[1]}; bnd.sshapedim = {{1;1;1;1},{1;1;1;1},{1;1;1;1},{1;1;1;1},{1;1;1;1}}; bnd.gporder = {{1;1;1;1},{1;1;1;1},{1;1;1;1},{1;1;1;1},{1;1;1;1}}; bnd.cporder = {{1;1;1;1},{1;1;1;1},{1;1;1;1},{1;1;1;1},{1;1;1;1}}; bnd.init = {{'';'';'';''},{'';'';'';''},{'';'';'';''},{'';'';'';''},{'';''; ... '';''}}; bnd.dinit = {{'';'';'';''},{'';'';'';''},{'';'';'';''},{'';'';'';''},{'';''; ... '';''}}; bnd.ind = [1,2,1,3,1,4,5,2,2]; bnd.dim = {'p','uor2','w2','p2'}; bnd.var = {'nU_esdl',{'+(u_esdl*nr_esdl)+(v_esdl*nz_esdl)','+(u_esdl*nr_esdl)+(v_esdl*nz_esdl)', ... '+(u_esdl*nr_esdl)+(v_esdl*nz_esdl)','+(u_esdl*nr_esdl)+(v_esdl*nz_esdl)', ... '+(u_esdl*nr_esdl)+(v_esdl*nz_esdl)'}, ... 'flux_esdl',{'+(u_esdl*nr_esdl)+(v_esdl*nz_esdl)','+(u_esdl*nr_esdl)+(v_esdl*nz_esdl)', ... '+(u_esdl*nr_esdl)+(v_esdl*nz_esdl)','+(u_esdl*nr_esdl)+(v_esdl*nz_esdl)', ... '+(u_esdl*nr_esdl)+(v_esdl*nz_esdl)'}, ... 'Frg_smaxi',{'0','0','0','0','0'}, ... 'Fzg_smaxi',{'0','0','0','0','0'}, ... 'disp_smaxi',{'sqrt(real(uaxi_smaxi)^2+real(w2)^2)','sqrt(real(uaxi_smaxi)^2+real(w2)^2)', ... 'sqrt(real(uaxi_smaxi)^2+real(w2)^2)','sqrt(real(uaxi_smaxi)^2+real(w2)^2)', ... 'sqrt(real(uaxi_smaxi)^2+real(w2)^2)'}, ... 'uaxi_smaxi',{'uor2*r','uor2*r','uor2*r','uor2*r','uor2*r'}, ... 'uaxir_smaxi',{'uor2r*r+uor2','uor2r*r+uor2','uor2r*r+uor2','uor2r*r+uor2', ... 'uor2r*r+uor2'}, ... 'uaxiz_smaxi',{'uor2z*r','uor2z*r','uor2z*r','uor2z*r','uor2z*r'}, ... 'uaxi_t_smaxi',{'uor2t*r','uor2t*r','uor2t*r','uor2t*r','uor2t*r'}, ... 'Tar_smaxi',{'((sr_smaxi)*(nr_smaxi)+(srz_smaxi)*(nz_smaxi))','((sr_smaxi)*(nr_smaxi)+(srz_smaxi)*(nz_smaxi))', ... '((sr_smaxi)*(nr_smaxi)+(srz_smaxi)*(nz_smaxi))','((sr_smaxi)*(nr_smaxi)+(srz_smaxi)*(nz_smaxi))', ... '((sr_smaxi)*(nr_smaxi)+(srz_smaxi)*(nz_smaxi))'}, ... 'Taz_smaxi',{'((srz_smaxi)*(nr_smaxi)+(sz_smaxi)*(nz_smaxi))','((srz_smaxi)*(nr_smaxi)+(sz_smaxi)*(nz_smaxi))', ... '((srz_smaxi)*(nr_smaxi)+(sz_smaxi)*(nz_smaxi))','((srz_smaxi)*(nr_smaxi)+(sz_smaxi)*(nz_smaxi))', ... '((srz_smaxi)*(nr_smaxi)+(sz_smaxi)*(nz_smaxi))'},'H0_esdl',{'0','0','0','0','0'}, ... 'Hp0_esdl',{'0','0','0','0','0'}, ... 'p0_esdl',{'0','0','0','0','0'}, ... 'N0_esdl',{'0','0','1e-10','0','0'}, ... 'Hb_esdl',{'0','0','0','0','0'}, ... 'Hpb_esdl',{'0','0','0','0','0'}, ... 'pb_esdl',{'0','0','0','0','0'}, ... 'Rb_esdl',{'0','0','0','0','0'}, ... 'Db_esdl',{'0','0','0','0','0'}, ... 'nr_esdl',{'nr','nr','nr','nr','nr'}, ... 'nz_esdl',{'nz','nz','nz','nz','nz'}, ... 'Fr_smaxi',{'0','0','0','0','0'}, ... 'FrPh_smaxi',{'0','0','0','0','0'}, ... 'Fz_smaxi',{'0','0','0','0','0'}, ... 'FzPh_smaxi',{'0','0','0','0','0'}, ... 'P_smaxi',{'0','0','0','0','0'}, ... 'PPh_smaxi',{'0','0','0','0','0'}, ... 'velr_smaxi',{'0','0','0','0','0'}, ... 'accr_smaxi',{'0','0','0','0','0'}, ... 'velz_smaxi',{'0','0','0','0','0'}, ... 'accz_smaxi',{'0','0','0','0','0'}, ... 'Rr_smaxi',{'0','0','0','0','0'}, ... 'Rz_smaxi',{'0','0','0','0','0'}, ... 'nr_smaxi',{'nr','nr','nr','nr','nr'}, ... 'nz_smaxi',{'nz','nz','nz','nz','nz'}, ... 'H11_smaxi',{'(((1))*(1)+(0)*(0))','(((1))*(1)+(0)*(0))','((0)*(1)+(0)*(0))', ... '((0)*(1)+(0)*(0))','((0)*(1)+(0)*(0))'}, ... 'H12_smaxi',{'(((1))*(0)+(0)*(1))','(((1))*(0)+(0)*(1))','((0)*(0)+(0)*(1))', ... '((0)*(0)+(0)*(1))','((0)*(0)+(0)*(1))'}, ... 'H21_smaxi',{'((0)*(1)+(0)*(0))','((0)*(1)+((1))*(0))','((0)*(1)+(0)*(0))', ... '((0)*(1)+(0)*(0))','((0)*(1)+(0)*(0))'}, ... 'H22_smaxi',{'((0)*(0)+(0)*(1))','((0)*(0)+((1))*(1))','((0)*(0)+(0)*(1))', ... '((0)*(0)+(0)*(1))','((0)*(0)+(0)*(1))'}}; bnd.expr = {}; bnd.lock = [0,0,0,0,0,0,0,0,0]; bnd.mlock = {[0,0,0,0,0,0,0,0,0]}; xfem.bnd = bnd; clear pnt pnt.weak = {{'0';'(0*uor2_test*r)/(2*pi)*1';'(0*w2_test)/(2*pi)*1';'0'},{'0'; ... '(0*uor2_test*r)/(2*pi)*1';'(0*w2_test)/(2*pi)*1';'0'}}; pnt.dweak = {{'0';'0';'0';'0'},{'0';'0';'0';'0'}}; pnt.constr = {{'0';'0-((((0)*(1)+(0)*(0)))*r*uor2+(((0)*(0)+(0)*(1)))*w2)'; ... '0-((((0)*(1)+(0)*(0)))*r*uor2+(((0)*(0)+(0)*(1)))*w2)';'0'},{'0';'0-((((0)*(1)+(0)*(0)))*r*uor2+(((0)*(0)+(0)*(1)))*w2)'; ... '0-((((0)*(1)+(0)*(0)))*r*uor2+(((0)*(0)+(0)*(1)))*w2)';'0'}}; pnt.constrf = {{'0';'0';'0';'0'},{'0';'0';'0';'0'}}; pnt.shape = {[1;2;3;4],[1;2;3;4]}; pnt.sshape = {[1],[1]}; pnt.sshapedim = {{1;1;1;1},{1;1;1;1}}; pnt.init = {{'';'';'';''},{'';'';'';''}}; pnt.dinit = {{'';'';'';''},{'';'';'';''}}; pnt.ind = [1,1,2,1,1,1,2,1]; pnt.dim = {'p','uor2','w2','p2'}; pnt.var = {'rhof_esdl',{'(+(rho_f))/1','(+(rho_f)+(rho_f))/2'}, ... 'Frg_smaxi',{'0','0'}, ... 'Fzg_smaxi',{'0','0'}, ... 'disp_smaxi',{'sqrt(real(uaxi_smaxi)^2+real(w2)^2)','sqrt(real(uaxi_smaxi)^2+real(w2)^2)'}, ... 'uaxi_smaxi',{'uor2*r','uor2*r'}, ... 'uaxir_smaxi',{'uor2r*r+uor2','uor2r*r+uor2'}, ... 'uaxiz_smaxi',{'uor2z*r','uor2z*r'}, ... 'uaxi_t_smaxi',{'uor2t*r','uor2t*r'},'Fr_smaxi',{'0','0'}, ... 'FrPh_smaxi',{'0','0'}, ... 'Fz_smaxi',{'0','0'}, ... 'FzPh_smaxi',{'0','0'}, ... 'Rr_smaxi',{'0','0'}, ... 'Rz_smaxi',{'0','0'}, ... 'H11_smaxi',{'((0)*(1)+(0)*(0))','((0)*(1)+(0)*(0))'}, ... 'H12_smaxi',{'((0)*(0)+(0)*(1))','((0)*(0)+(0)*(1))'}, ... 'H21_smaxi',{'((0)*(1)+(0)*(0))','((0)*(1)+(0)*(0))'}, ... 'H22_smaxi',{'((0)*(0)+(0)*(1))','((0)*(0)+(0)*(1))'}}; pnt.expr = {}; pnt.lock = [0,0,0,0,0,0,0,0]; pnt.mlock = {[0,0,0,0,0,0,0,0]}; xfem.pnt = pnt; xfem.var = {'tscale_esdl','1e-5','g_esdl','9.81','D_esdl','0','t_old_ini_smaxi','-1'}; xfem.expr = {'F_f2sr','-alpha_b*rho_f*g_esdl*Hr_esdl', ... 'F_f2sz','-alpha_b*rho_f*g_esdl*Hz_esdl'}; clear descr descr.expr= {'F_f2sr','Fluid-to-structure coupling expression, r direction','F_f2sz','Fluid-to-structure coupling expression, z direction'}; xfem.descr = descr; clear elemmph clear elem elem.elem = 'elpric'; elem.g = {'1'}; clear equ equ.tensor = {{'real(sr_smaxi)','0','real(sphi_smaxi)','real(srz_smaxi)','0', ... 'real(sz_smaxi)'}}; equ.ind = {{'1','2'}}; elem.geomdim = {{{},{},equ}}; elem.basename = 's'; elem.postname = 'smaxi'; elem.sdim = {'r','phi','z'}; elemmph{1} = elem; clear elem elem.elem = 'elpric'; elem.g = {'1'}; clear equ equ.tensor = {{'real(er_smaxi)','0','real(ephi_smaxi)','real(erz_smaxi)','0', ... 'real(ez_smaxi)'}}; equ.ind = {{'1','2'}}; elem.geomdim = {{{},{},equ}}; elem.basename = 'e'; elem.postname = 'smaxi'; elem.sdim = {'r','phi','z'}; elemmph{2} = elem; xfem.elemmph = elemmph; clear draw draw.p.objs = {}; draw.p.name = {}; draw.c.objs = {flbinary('g20','draw','surficial_caprock_model.mph'),flbinary('g8','draw','surficial_caprock_model.mph')}; draw.c.name = {'B2','B1'}; draw.s.objs = {flbinary('g1','draw','surficial_caprock_model.mph')}; draw.s.name = {'R1'}; xfem.draw = draw; xfem.const = {'rho_f','1000[kg/m^3]','rho_s','2780[kg/m^3]','alpha_b','1'}; xfem.globalexpr = {}; clear fcns xfem.functions = {}; descr = xfem.descr; descr.const= {'rho_f','Fluid density','alpha_b','Biot-Willis coefficient','rho_s','Solids density'}; xfem.descr = descr; xfem.xmcases = [0]; xfem.mcases = [0]; flbinary clear; xfem.rulingmode = 'esdl'; xfem.solform = 'weak'; clear ode clear units; units.basesystem = 'SI'; ode.units = units; xfem.ode=ode;uq~ tgui.solvemodel.toutcomp='off'; gui.solvemodel.currsolver='time'; gui.solvemodel.solveroption=''; gui.solvemodel.postsolver='stationary'; gui.solvemodel.nonlin='auto'; gui.solvemodel.ntol='1.0E-6'; gui.solvemodel.maxiter='25'; gui.solvemodel.maxsegiter='100'; gui.solvemodel.manualdamp='off'; gui.solvemodel.damping='on'; 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='0.0010'; gui.solvemodel.augmaxiter='25'; gui.solvemodel.augsolver='lumped'; gui.solvemodel.atol='0.0010'; gui.solvemodel.rtol='0.01'; gui.solvemodel.tlist='0:164000:164000*2000'; 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.minorder='1'; gui.solvemodel.maxstep='1.0'; gui.solvemodel.timeusestopcond='off'; gui.solvemodel.timestopcond=''; 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='2'; gui.solvemodel.tpmult='1.7'; gui.solvemodel.tpworst='0.5'; gui.solvemodel.tpfract='0.5'; gui.solvemodel.autosolver='on'; gui.solvemodel.varcomp=''; gui.solvemodel.oldvarcomp=''; gui.solvemodel.manualoptlim='off'; gui.solvemodel.manualhessupd='off'; gui.solvemodel.qpsolver='cholesky'; gui.solvemodel.gradient='adjoint'; gui.solvemodel.hessupd='10'; gui.solvemodel.totitlim='10000'; gui.solvemodel.itlim='500'; gui.solvemodel.majitlim='1000'; gui.solvemodel.opttol='1.0e-6'; gui.solvemodel.feastol='1.0e-6'; gui.solvemodel.majfeastol='1.0e-6'; gui.solvemodel.solcomp='p2,p,uor2,w2'; gui.solvemodel.outcomp='p2,p,uor2,w2'; 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.pardreorder='nd'; gui.solvemodel.linsolvernode.pardrreorder='on'; gui.solvemodel.linsolvernode.pivotperturb='1.0E-8'; gui.solvemodel.linsolvernode.errorchk='on'; gui.solvemodel.linsolvernode.errorchkd='off'; gui.solvemodel.linsolvernode.termination='tol'; 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.prefuntype='left'; gui.solvemodel.linsolvernode.iluiter='1'; 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.vankablocked='off'; gui.solvemodel.linsolvernode.sorblocked='off'; gui.solvemodel.linsolvernode.sorvecdof=''; 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.solversegmodel.seggrps{1}.segcomp='p'; gui.solvemodel.solversegmodel.seggrps{1}.ntol='1e-3'; gui.solvemodel.solversegmodel.seggrps{1}.linsolvernode.currlinsolver='umfpack'; gui.solvemodel.solversegmodel.seggrps{1}.linsolvernode.type='linsolver'; gui.solvemodel.solversegmodel.seggrps{1}.linsolvernode.droptol='0.0'; gui.solvemodel.solversegmodel.seggrps{1}.linsolvernode.thresh='0.1'; gui.solvemodel.solversegmodel.seggrps{1}.linsolvernode.umfalloc='0.7'; gui.solvemodel.solversegmodel.seggrps{1}.linsolvernode.preorder='mmd'; gui.solvemodel.solversegmodel.seggrps{1}.linsolvernode.preroworder='on'; gui.solvemodel.solversegmodel.seggrps{1}.linsolvernode.pardreorder='nd'; gui.solvemodel.solversegmodel.seggrps{1}.linsolvernode.pardrreorder='on'; gui.solvemodel.solversegmodel.seggrps{1}.linsolvernode.pivotperturb='1.0E-8'; gui.solvemodel.solversegmodel.seggrps{1}.linsolvernode.errorchk='on'; gui.solvemodel.solversegmodel.seggrps{1}.linsolvernode.errorchkd='off'; gui.solvemodel.solversegmodel.seggrps{1}.linsolvernode.termination='tol'; gui.solvemodel.solversegmodel.seggrps{1}.linsolvernode.iter='2'; gui.solvemodel.solversegmodel.seggrps{1}.linsolvernode.itol='1.0E-6'; gui.solvemodel.solversegmodel.seggrps{1}.linsolvernode.rhob='400.0'; gui.solvemodel.solversegmodel.seggrps{1}.linsolvernode.maxlinit='10000'; gui.solvemodel.solversegmodel.seggrps{1}.linsolvernode.prefuntype='left'; gui.solvemodel.solversegmodel.seggrps{1}.linsolvernode.iluiter='1'; gui.solvemodel.solversegmodel.seggrps{1}.linsolvernode.itrestart='50'; gui.solvemodel.solversegmodel.seggrps{1}.linsolvernode.seconditer='1'; gui.solvemodel.solversegmodel.seggrps{1}.linsolvernode.relax='1.0'; gui.solvemodel.solversegmodel.seggrps{1}.linsolvernode.amgauto='3'; gui.solvemodel.solversegmodel.seggrps{1}.linsolvernode.mglevels='6'; gui.solvemodel.solversegmodel.seggrps{1}.linsolvernode.mgcycle='v'; gui.solvemodel.solversegmodel.seggrps{1}.linsolvernode.maxcoarsedof='5000'; gui.solvemodel.solversegmodel.seggrps{1}.linsolvernode.oocmemory='128.0'; gui.solvemodel.solversegmodel.seggrps{1}.linsolvernode.oocfilename=''; gui.solvemodel.solversegmodel.seggrps{1}.linsolvernode.modified='off'; gui.solvemodel.solversegmodel.seggrps{1}.linsolvernode.fillratio='2.0'; gui.solvemodel.solversegmodel.seggrps{1}.linsolvernode.respectpattern='on'; gui.solvemodel.solversegmodel.seggrps{1}.linsolvernode.droptype='droptol'; gui.solvemodel.solversegmodel.seggrps{1}.linsolvernode.vankavars=''; gui.solvemodel.solversegmodel.seggrps{1}.linsolvernode.vankasolv='gmres'; gui.solvemodel.solversegmodel.seggrps{1}.linsolvernode.vankatol='0.02'; gui.solvemodel.solversegmodel.seggrps{1}.linsolvernode.vankarestart='100'; gui.solvemodel.solversegmodel.seggrps{1}.linsolvernode.vankarelax='0.8'; gui.solvemodel.solversegmodel.seggrps{1}.linsolvernode.vankablocked='off'; gui.solvemodel.solversegmodel.seggrps{1}.linsolvernode.sorblocked='off'; gui.solvemodel.solversegmodel.seggrps{1}.linsolvernode.sorvecdof=''; gui.solvemodel.solversegmodel.seggrps{1}.linsolvernode.mgauto='shape'; gui.solvemodel.solversegmodel.seggrps{1}.linsolvernode.rmethod='regular'; gui.solvemodel.solversegmodel.seggrps{1}.linsolvernode.coarseassem='on'; gui.solvemodel.solversegmodel.seggrps{1}.linsolvernode.meshscale='2'; gui.solvemodel.solversegmodel.seggrps{1}.linsolvernode.mgautolevels='2'; gui.solvemodel.solversegmodel.seggrps{1}.linsolvernode.mgkeep='off'; gui.solvemodel.solversegmodel.seggrps{1}.linsolvernode.mggeom='Geom1'; gui.solvemodel.solversegmodel.seggrps{1}.linsolvernode.mcase0='on'; gui.solvemodel.solversegmodel.seggrps{1}.linsolvernode.mgassem0='on'; gui.solvemodel.solversegmodel.seggrps{2}.segcomp='uor w'; gui.solvemodel.solversegmodel.seggrps{2}.ntol='1e-3'; gui.solvemodel.solversegmodel.seggrps{2}.linsolvernode.currlinsolver='umfpack'; gui.solvemodel.solversegmodel.seggrps{2}.linsolvernode.type='linsolver'; gui.solvemodel.solversegmodel.seggrps{2}.linsolvernode.droptol='0.0'; gui.solvemodel.solversegmodel.seggrps{2}.linsolvernode.thresh='0.1'; gui.solvemodel.solversegmodel.seggrps{2}.linsolvernode.umfalloc='0.7'; gui.solvemodel.solversegmodel.seggrps{2}.linsolvernode.preorder='mmd'; gui.solvemodel.solversegmodel.seggrps{2}.linsolvernode.preroworder='on'; gui.solvemodel.solversegmodel.seggrps{2}.linsolvernode.pardreorder='nd'; gui.solvemodel.solversegmodel.seggrps{2}.linsolvernode.pardrreorder='on'; gui.solvemodel.solversegmodel.seggrps{2}.linsolvernode.pivotperturb='1.0E-8'; gui.solvemodel.solversegmodel.seggrps{2}.linsolvernode.errorchk='on'; gui.solvemodel.solversegmodel.seggrps{2}.linsolvernode.errorchkd='off'; gui.solvemodel.solversegmodel.seggrps{2}.linsolvernode.termination='tol'; gui.solvemodel.solversegmodel.seggrps{2}.linsolvernode.iter='2'; gui.solvemodel.solversegmodel.seggrps{2}.linsolvernode.itol='1.0E-6'; gui.solvemodel.solversegmodel.seggrps{2}.linsolvernode.rhob='400.0'; gui.solvemodel.solversegmodel.seggrps{2}.linsolvernode.maxlinit='10000'; gui.solvemodel.solversegmodel.seggrps{2}.linsolvernode.prefuntype='left'; gui.solvemodel.solversegmodel.seggrps{2}.linsolvernode.iluiter='1'; gui.solvemodel.solversegmodel.seggrps{2}.linsolvernode.itrestart='50'; gui.solvemodel.solversegmodel.seggrps{2}.linsolvernode.seconditer='1'; gui.solvemodel.solversegmodel.seggrps{2}.linsolvernode.relax='1.0'; gui.solvemodel.solversegmodel.seggrps{2}.linsolvernode.amgauto='3'; gui.solvemodel.solversegmodel.seggrps{2}.linsolvernode.mglevels='6'; gui.solvemodel.solversegmodel.seggrps{2}.linsolvernode.mgcycle='v'; gui.solvemodel.solversegmodel.seggrps{2}.linsolvernode.maxcoarsedof='5000'; gui.solvemodel.solversegmodel.seggrps{2}.linsolvernode.oocmemory='128.0'; gui.solvemodel.solversegmodel.seggrps{2}.linsolvernode.oocfilename=''; gui.solvemodel.solversegmodel.seggrps{2}.linsolvernode.modified='off'; gui.solvemodel.solversegmodel.seggrps{2}.linsolvernode.fillratio='2.0'; gui.solvemodel.solversegmodel.seggrps{2}.linsolvernode.respectpattern='on'; gui.solvemodel.solversegmodel.seggrps{2}.linsolvernode.droptype='droptol'; gui.solvemodel.solversegmodel.seggrps{2}.linsolvernode.vankavars=''; gui.solvemodel.solversegmodel.seggrps{2}.linsolvernode.vankasolv='gmres'; gui.solvemodel.solversegmodel.seggrps{2}.linsolvernode.vankatol='0.02'; gui.solvemodel.solversegmodel.seggrps{2}.linsolvernode.vankarestart='100'; gui.solvemodel.solversegmodel.seggrps{2}.linsolvernode.vankarelax='0.8'; gui.solvemodel.solversegmodel.seggrps{2}.linsolvernode.vankablocked='off'; gui.solvemodel.solversegmodel.seggrps{2}.linsolvernode.sorblocked='off'; gui.solvemodel.solversegmodel.seggrps{2}.linsolvernode.sorvecdof=''; gui.solvemodel.solversegmodel.seggrps{2}.linsolvernode.mgauto='shape'; gui.solvemodel.solversegmodel.seggrps{2}.linsolvernode.rmethod='regular'; gui.solvemodel.solversegmodel.seggrps{2}.linsolvernode.coarseassem='on'; gui.solvemodel.solversegmodel.seggrps{2}.linsolvernode.meshscale='2'; gui.solvemodel.solversegmodel.seggrps{2}.linsolvernode.mgautolevels='2'; gui.solvemodel.solversegmodel.seggrps{2}.linsolvernode.mgkeep='off'; gui.solvemodel.solversegmodel.seggrps{2}.linsolvernode.mggeom='Geom1'; gui.solvemodel.solversegmodel.seggrps{2}.linsolvernode.mcase0='on'; gui.solvemodel.solversegmodel.seggrps{2}.linsolvernode.mgassem0='on'; gui.solvemodel.solversegmodel.seggrps{3}.segcomp='uor2 w2 p2'; gui.solvemodel.solversegmodel.seggrps{3}.ntol='1e-3'; gui.solvemodel.solversegmodel.seggrps{3}.linsolvernode.currlinsolver='umfpack'; gui.solvemodel.solversegmodel.seggrps{3}.linsolvernode.type='linsolver'; gui.solvemodel.solversegmodel.seggrps{3}.linsolvernode.droptol='0.0'; gui.solvemodel.solversegmodel.seggrps{3}.linsolvernode.thresh='0.1'; gui.solvemodel.solversegmodel.seggrps{3}.linsolvernode.umfalloc='0.7'; gui.solvemodel.solversegmodel.seggrps{3}.linsolvernode.preorder='mmd'; gui.solvemodel.solversegmodel.seggrps{3}.linsolvernode.preroworder='on'; gui.solvemodel.solversegmodel.seggrps{3}.linsolvernode.pardreorder='nd'; gui.solvemodel.solversegmodel.seggrps{3}.linsolvernode.pardrreorder='on'; gui.solvemodel.solversegmodel.seggrps{3}.linsolvernode.pivotperturb='1.0E-8'; gui.solvemodel.solversegmodel.seggrps{3}.linsolvernode.errorchk='on'; gui.solvemodel.solversegmodel.seggrps{3}.linsolvernode.errorchkd='off'; gui.solvemodel.solversegmodel.seggrps{3}.linsolvernode.termination='tol'; gui.solvemodel.solversegmodel.seggrps{3}.linsolvernode.iter='2'; gui.solvemodel.solversegmodel.seggrps{3}.linsolvernode.itol='1.0E-6'; gui.solvemodel.solversegmodel.seggrps{3}.linsolvernode.rhob='400.0'; gui.solvemodel.solversegmodel.seggrps{3}.linsolvernode.maxlinit='10000'; gui.solvemodel.solversegmodel.seggrps{3}.linsolvernode.prefuntype='left'; gui.solvemodel.solversegmodel.seggrps{3}.linsolvernode.iluiter='1'; gui.solvemodel.solversegmodel.seggrps{3}.linsolvernode.itrestart='50'; gui.solvemodel.solversegmodel.seggrps{3}.linsolvernode.seconditer='1'; gui.solvemodel.solversegmodel.seggrps{3}.linsolvernode.relax='1.0'; gui.solvemodel.solversegmodel.seggrps{3}.linsolvernode.amgauto='3'; gui.solvemodel.solversegmodel.seggrps{3}.linsolvernode.mglevels='6'; gui.solvemodel.solversegmodel.seggrps{3}.linsolvernode.mgcycle='v'; gui.solvemodel.solversegmodel.seggrps{3}.linsolvernode.maxcoarsedof='5000'; gui.solvemodel.solversegmodel.seggrps{3}.linsolvernode.oocmemory='128.0'; gui.solvemodel.solversegmodel.seggrps{3}.linsolvernode.oocfilename=''; gui.solvemodel.solversegmodel.seggrps{3}.linsolvernode.modified='off'; gui.solvemodel.solversegmodel.seggrps{3}.linsolvernode.fillratio='2.0'; gui.solvemodel.solversegmodel.seggrps{3}.linsolvernode.respectpattern='on'; gui.solvemodel.solversegmodel.seggrps{3}.linsolvernode.droptype='droptol'; gui.solvemodel.solversegmodel.seggrps{3}.linsolvernode.vankavars=''; gui.solvemodel.solversegmodel.seggrps{3}.linsolvernode.vankasolv='gmres'; gui.solvemodel.solversegmodel.seggrps{3}.linsolvernode.vankatol='0.02'; gui.solvemodel.solversegmodel.seggrps{3}.linsolvernode.vankarestart='100'; gui.solvemodel.solversegmodel.seggrps{3}.linsolvernode.vankarelax='0.8'; gui.solvemodel.solversegmodel.seggrps{3}.linsolvernode.vankablocked='off'; gui.solvemodel.solversegmodel.seggrps{3}.linsolvernode.sorblocked='off'; gui.solvemodel.solversegmodel.seggrps{3}.linsolvernode.sorvecdof=''; gui.solvemodel.solversegmodel.seggrps{3}.linsolvernode.mgauto='shape'; gui.solvemodel.solversegmodel.seggrps{3}.linsolvernode.rmethod='regular'; gui.solvemodel.solversegmodel.seggrps{3}.linsolvernode.coarseassem='on'; gui.solvemodel.solversegmodel.seggrps{3}.linsolvernode.meshscale='2'; gui.solvemodel.solversegmodel.seggrps{3}.linsolvernode.mgautolevels='2'; gui.solvemodel.solversegmodel.seggrps{3}.linsolvernode.mgkeep='off'; gui.solvemodel.solversegmodel.seggrps{3}.linsolvernode.mggeom='Geom1'; gui.solvemodel.solversegmodel.seggrps{3}.linsolvernode.mcase0='on'; gui.solvemodel.solversegmodel.seggrps{3}.linsolvernode.mgassem0='on'; gui.solvemodel.solversegmodel.segorder=[1,2,3]; gui.solvemodel.solversegmodel.segdamp=[0.5,0.5,0.5]; gui.solvemodel.solversegmodel.subiter=[1,1,1]; gui.solvemodel.solversegmodel.manualsteps='off'; gui.solvemodel.solversegmodel.llimitdof=''; gui.solvemodel.solversegmodel.llimitval=''; 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='1035819'; gui.registry.saved_version='COMSOL 3.4.0.248'; 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='1454226810866'; gui.registry.info_creationdate='1362358641673'; gui.registry.info_modelresult=''; gui.registry.spice_netlist=''; gui.registry.spice_forceac='off'; gui.registry.spice_includeparams='off'; 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='boundary'; gui.flmodel{1}.currappl='0'; gui.flmodel{1}.axis.xmin='-50.0'; gui.flmodel{1}.axis.xmax='1050.0'; gui.flmodel{1}.axis.ymin='-1061.0'; gui.flmodel{1}.axis.ymax='61.0'; gui.flmodel{1}.axis.zmin='-1.0'; gui.flmodel{1}.axis.zmax='1.0'; gui.flmodel{1}.axis.xspacing='10.0'; gui.flmodel{1}.axis.yspacing='10.0'; gui.flmodel{1}.axis.zspacing='0.2'; gui.flmodel{1}.axis.extrax=''; gui.flmodel{1}.axis.extray=''; gui.flmodel{1}.axis.extraz=''; gui.flmodel{1}.camera.xmin='-0.55'; gui.flmodel{1}.camera.xmax='0.55'; gui.flmodel{1}.camera.ymin='-0.5609999999999999'; gui.flmodel{1}.camera.ymax='0.5609999999999999'; gui.flmodel{1}.camera.camposx='0.5'; gui.flmodel{1}.camera.camposy='-0.5'; gui.flmodel{1}.camera.camposz='5.609999999999999'; gui.flmodel{1}.camera.camtargetx='0.5'; gui.flmodel{1}.camera.camtargety='-0.5'; gui.flmodel{1}.camera.camtargetz='0.0'; gui.flmodel{1}.camera.camupx='0.0'; gui.flmodel{1}.camera.camupy='1.0'; gui.flmodel{1}.camera.camupz='0.0'; 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='0.5'; gui.flmodel{1}.lightmodel.headlight.ypos='-0.5'; gui.flmodel{1}.lightmodel.headlight.zpos='5.61'; 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='off'; gui.flmodel{1}.registry.axis_csys='on'; gui.flmodel{1}.registry.grid_visible='on'; gui.flmodel{1}.registry.grid_auto='off'; 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='0.0010'; gui.flmodel{1}.registry.render_mesh='off'; 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='5'; gui.flmodel{1}.meshparam.usehauto='on'; gui.flmodel{1}.meshparam.hmax=''; gui.flmodel{1}.meshparam.hmaxfact='1'; gui.flmodel{1}.meshparam.hcurve='0.3'; gui.flmodel{1}.meshparam.hgrad='1.3'; gui.flmodel{1}.meshparam.hcutoff='0.001'; gui.flmodel{1}.meshparam.hnarrow='1'; gui.flmodel{1}.meshparam.hpnt='10'; gui.flmodel{1}.meshparam.xscale='1.0'; gui.flmodel{1}.meshparam.yscale='1.0'; gui.flmodel{1}.meshparam.jiggle='on'; gui.flmodel{1}.meshparam.mcase='0'; gui.flmodel{1}.meshparam.rmethod='regular'; gui.flmodel{1}.meshparam.hmaxvtx={'','','','','','','',''}; gui.flmodel{1}.meshparam.hgradvtx={'','','','','','','',''}; gui.flmodel{1}.meshparam.hmaxedg={'','100','','0.5','','','2','100','30'}; gui.flmodel{1}.meshparam.hcutoffedg={'','','','','','','','',''}; gui.flmodel{1}.meshparam.hcurveedg={'','','','','','','','',''}; gui.flmodel{1}.meshparam.hgradedg={'','','','','','','','',''}; gui.flmodel{1}.meshparam.hgradsub={'',''}; gui.flmodel{1}.meshparam.methodsub={'tri','tri'}; gui.flmodel{1}.meshparam.hmaxsub={'100','5'}; gui.flmodel{1}.postmodel.postplot.triplot='on'; gui.flmodel{1}.postmodel.postplot.tridata={'w2'}; gui.flmodel{1}.postmodel.postplot.trirangeauto='on'; gui.flmodel{1}.postmodel.postplot.trirangemin='0.0'; gui.flmodel{1}.postmodel.postplot.trirangemax='3.539723330024922E-7'; gui.flmodel{1}.postmodel.postplot.tricont='on'; gui.flmodel{1}.postmodel.postplot.triunit='m'; gui.flmodel{1}.postmodel.postplot.triheightdata={'H_esdl'}; gui.flmodel{1}.postmodel.postplot.triheightdatacheck='off'; gui.flmodel{1}.postmodel.postplot.triheightunit='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.contplot='off'; gui.flmodel{1}.postmodel.postplot.contdata={'H_esdl'}; gui.flmodel{1}.postmodel.postplot.contcont='on'; gui.flmodel{1}.postmodel.postplot.contunit='m'; gui.flmodel{1}.postmodel.postplot.contheightdata={'H_esdl'}; gui.flmodel{1}.postmodel.postplot.contheightdatacheck='off'; gui.flmodel{1}.postmodel.postplot.contheightunit='m'; gui.flmodel{1}.postmodel.postplot.contcolordata={'H_esdl'}; gui.flmodel{1}.postmodel.postplot.contcolorrangeauto='on'; gui.flmodel{1}.postmodel.postplot.contcolorrangemin=''; gui.flmodel{1}.postmodel.postplot.contcolorrangemax=''; gui.flmodel{1}.postmodel.postplot.contcolordatacheck='off'; gui.flmodel{1}.postmodel.postplot.contcolorunit='m'; gui.flmodel{1}.postmodel.postplot.contmap='cool'; gui.flmodel{1}.postmodel.postplot.contmapdepth='1024'; gui.flmodel{1}.postmodel.postplot.contbar='on'; gui.flmodel{1}.postmodel.postplot.contusemap='on'; gui.flmodel{1}.postmodel.postplot.contcolor='255,0,0'; gui.flmodel{1}.postmodel.postplot.contlevels='20'; gui.flmodel{1}.postmodel.postplot.contvectorlevels=''; gui.flmodel{1}.postmodel.postplot.contisvector='off'; gui.flmodel{1}.postmodel.postplot.contlabel='off'; gui.flmodel{1}.postmodel.postplot.contfill='off'; gui.flmodel{1}.postmodel.postplot.linplot='off'; gui.flmodel{1}.postmodel.postplot.lindata={'H_esdl'}; 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.linheightdata={'H_esdl'}; gui.flmodel{1}.postmodel.postplot.linheightdatacheck='off'; gui.flmodel{1}.postmodel.postplot.linheightunit='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={'Hr_esdl','Hz_esdl'}; gui.flmodel{1}.postmodel.postplot.arrowunit='1'; gui.flmodel{1}.postmodel.postplot.arrowbnddata={'uaxi_smaxi','w2'}; gui.flmodel{1}.postmodel.postplot.arrowbndunit='m'; gui.flmodel{1}.postmodel.postplot.arrowheightdata={'H_esdl'}; gui.flmodel{1}.postmodel.postplot.arrowheightdatacheck='off'; gui.flmodel{1}.postmodel.postplot.arrowheightunit='m'; gui.flmodel{1}.postmodel.postplot.arrowxspacing='15'; gui.flmodel{1}.postmodel.postplot.arrowxvectorspacing=''; gui.flmodel{1}.postmodel.postplot.arrowxisvector='off'; gui.flmodel{1}.postmodel.postplot.arrowyspacing='15'; gui.flmodel{1}.postmodel.postplot.arrowyvectorspacing=''; gui.flmodel{1}.postmodel.postplot.arrowyisvector='off'; gui.flmodel{1}.postmodel.postplot.arrowtype='arrow'; 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.princdata={'s1_smaxi','s1r_smaxi','s1phi_smaxi','s1z_smaxi','s2_smaxi','s2r_smaxi','s2phi_smaxi','s2z_smaxi','s3_smaxi','s3r_smaxi','s3phi_smaxi','s3z_smaxi'}; gui.flmodel{1}.postmodel.postplot.princunit='Pa'; gui.flmodel{1}.postmodel.postplot.princheightdata={'H_esdl'}; gui.flmodel{1}.postmodel.postplot.princheightdatacheck='off'; gui.flmodel{1}.postmodel.postplot.princheightunit='m'; gui.flmodel{1}.postmodel.postplot.princxspacing='8'; gui.flmodel{1}.postmodel.postplot.princxvectorspacing=''; gui.flmodel{1}.postmodel.postplot.princxisvector='off'; gui.flmodel{1}.postmodel.postplot.princyspacing='8'; gui.flmodel{1}.postmodel.postplot.princyvectorspacing=''; gui.flmodel{1}.postmodel.postplot.princyisvector='off'; gui.flmodel{1}.postmodel.postplot.princtype='arrow'; 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={'Hr_esdl','Hz_esdl'}; gui.flmodel{1}.postmodel.postplot.flowunit='1'; gui.flmodel{1}.postmodel.postplot.flowuseexpression='off'; gui.flmodel{1}.postmodel.postplot.flowcolor='255,0,0'; gui.flmodel{1}.postmodel.postplot.flowcolordata={'H_esdl'}; gui.flmodel{1}.postmodel.postplot.flowcolorunit='m'; 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.flowheightdata={'H_esdl'}; gui.flmodel{1}.postmodel.postplot.flowheightdatacheck='off'; gui.flmodel{1}.postmodel.postplot.flowheightunit='m'; 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.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.05'; gui.flmodel{1}.postmodel.postplot.flowlinesvel='20'; gui.flmodel{1}.postmodel.postplot.flowseedmanual='off'; gui.flmodel{1}.postmodel.postplot.flowseed1=''; gui.flmodel{1}.postmodel.postplot.flowseed2=''; 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=''; gui.flmodel{1}.postmodel.postplot.partmass='1'; gui.flmodel{1}.postmodel.postplot.partforce={'','',''}; gui.flmodel{1}.postmodel.postplot.part_massless_flowdata={'Hr_esdl','Hz_esdl'}; gui.flmodel{1}.postmodel.postplot.part_massless_flowunit='1'; 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.partstartptssel='partstart_explicit'; gui.flmodel{1}.postmodel.postplot.partstartdl=''; gui.flmodel{1}.postmodel.postplot.partstartedim1levels='10'; gui.flmodel{1}.postmodel.postplot.partstartedim1vectorlevels=''; gui.flmodel{1}.postmodel.postplot.partstartedim1isvector='off'; gui.flmodel{1}.postmodel.postplot.explicitcoord={'0','0'}; gui.flmodel{1}.postmodel.postplot.partuseexpression='off'; gui.flmodel{1}.postmodel.postplot.partcolor='255,0,0'; gui.flmodel{1}.postmodel.postplot.partcolordata={'H_esdl'}; gui.flmodel{1}.postmodel.postplot.partcolorunit='m'; 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.partmasslessrtol='0.001'; gui.flmodel{1}.postmodel.postplot.partmasslessatolmanual='off'; gui.flmodel{1}.postmodel.postplot.partmasslessatol={''}; gui.flmodel{1}.postmodel.postplot.partmasslessatolexpanded={''}; 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.partatolexpanded={'',''}; 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={'H_esdl'}; gui.flmodel{1}.postmodel.postplot.maxminsubdatacheck='on'; gui.flmodel{1}.postmodel.postplot.maxminsubunit='m'; gui.flmodel{1}.postmodel.postplot.maxminbnddata={'H_esdl'}; gui.flmodel{1}.postmodel.postplot.maxminbnddatacheck='off'; gui.flmodel{1}.postmodel.postplot.maxminbndunit='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='3'; gui.flmodel{1}.postmodel.postplot.geomnum={'Geom1'}; gui.flmodel{1}.postmodel.postplot.phase='0'; 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.deformsubdata={'uaxi_smaxi','w2'}; gui.flmodel{1}.postmodel.postplot.deformsubunit='m'; gui.flmodel{1}.postmodel.postplot.deformbnddata={'uaxi_smaxi','w2'}; gui.flmodel{1}.postmodel.postplot.deformbndunit='m'; gui.flmodel{1}.postmodel.postplot.deformautoscale='on'; gui.flmodel{1}.postmodel.postplot.deformscale='2.980428e8'; 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={'H_esdl'}; 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={'H_esdl'}; gui.flmodel{1}.postmodel.intdata{2}.intunit='m^2'; gui.flmodel{1}.postmodel.intdata{2}.phase='0'; 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={'H_esdl'}; gui.flmodel{1}.postmodel.intdata{3}.intunit='m^3'; gui.flmodel{1}.postmodel.intdata{3}.phase='0'; gui.flmodel{1}.postmodel.intdata{3}.solnum='0'; gui.flmodel{1}.postmodel.intdata{3}.selectvia='stored'; gui.flmodel{1}.postmodel.domainplot.colordata={'H_esdl'}; 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.extrusion='off'; gui.flmodel{1}.postmodel.domainplot.lineyaxisdata={'H_esdl'}; 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={'H_esdl'}; 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={'H_esdl'}; 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={'H_esdl'}; 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.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.extrusion='off'; gui.flmodel{1}.postmodel.crossplot.lineyaxisdata={'w2'}; gui.flmodel{1}.postmodel.crossplot.lineyaxisunit='m'; gui.flmodel{1}.postmodel.crossplot.linexaxisxaxistype='arc'; gui.flmodel{1}.postmodel.crossplot.linexaxisuseexpr='off'; gui.flmodel{1}.postmodel.crossplot.linexaxisdata={'H_esdl'}; gui.flmodel{1}.postmodel.crossplot.linexaxisunit='m'; 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='50'; gui.flmodel{1}.postmodel.crossplot.linecoord={'0','300','10','10'}; gui.flmodel{1}.postmodel.crossplot.linelevels='5'; gui.flmodel{1}.postmodel.crossplot.linevectorlevels=''; gui.flmodel{1}.postmodel.crossplot.lineisvector='off'; gui.flmodel{1}.postmodel.crossplot.lineactive='off'; gui.flmodel{1}.postmodel.crossplot.pointyaxisdata={'w2'}; gui.flmodel{1}.postmodel.crossplot.pointyaxisunit='m'; gui.flmodel{1}.postmodel.crossplot.pointxxaxistype=''; gui.flmodel{1}.postmodel.crossplot.pointxuseexpr='off'; gui.flmodel{1}.postmodel.crossplot.pointxdata={'H_esdl'}; gui.flmodel{1}.postmodel.crossplot.pointxunit='m'; 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','10'}; gui.flmodel{1}.postmodel.crossplot.crossdispcolor='255,0,0'; gui.flmodel{1}.postmodel.crossplot.phase='0'; 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.plottypeind='0'; gui.flmodel{1}.postmodel.dataexport.pntdata={'H_esdl'}; 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.bnddata={'H_esdl'}; 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={'H_esdl'}; gui.flmodel{1}.postmodel.dataexport.subcont='off'; gui.flmodel{1}.postmodel.dataexport.subunit='m'; 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.sublocation='element'; gui.flmodel{1}.postmodel.dataexport.sublagorder='2'; gui.flmodel{1}.postmodel.dataexport.phase='0'; 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='2'; 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.solnum='0'; gui.flmodel{1}.postmodel.datadisplay.selectvia='stored'; gui.flmodel{1}.postmodel.datadisplay.interpdata={'H_esdl'}; gui.flmodel{1}.postmodel.datadisplay.interpunit='m'; gui.flmodel{1}.postmodel.datadisplay.coord={'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}.geommodel.suppressed{1}=[]; gui.flmodel{1}.geommodel.suppressed{2}=[]; gui.flmodel{1}.geommodel.suppressed{3}=[]; gui.flmodel{1}.workplane.type='0'; gui.flmodel{1}.workplane.wrkpln='0,1,0,0,0,1,0,0,0'; gui.flmodel{1}.workplane.localsys='0,1,0,0,0,0,1,0,0,0,0,1'; gui.flmodel{1}.workplane.model2d='Geom1'; gui.flmodel{1}.workplane.quicktype='10'; gui.flmodel{1}.workplane.parameter='0'; gui.flmodel{1}.workplane.zdir='up'; gui.flmodel{1}.meshmodel.meshplot.subplot='on'; gui.flmodel{1}.meshmodel.meshplot.meshsubuseexpression='off'; gui.flmodel{1}.meshmodel.meshplot.meshsubcolor='128,128,128'; gui.flmodel{1}.meshmodel.meshplot.meshsubbordercheck='off'; gui.flmodel{1}.meshmodel.meshplot.meshwiresubplot='on'; gui.flmodel{1}.meshmodel.meshplot.meshwiresubcolor='192,192,192'; gui.flmodel{1}.meshmodel.meshplot.bndplot='on'; gui.flmodel{1}.meshmodel.meshplot.meshbndcolor='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'; srcom.femlab.geom.Geom2Vc/Oxrcom.femlab.geom.Geom͹6{L geomAssoctLcom/femlab/geom/GeomAssoc;Lversionq~xpwq~wur[BTxpGeom2|=-C6?@@? BezierCurve?@@? AssocAttrib VectorInt$ xwxsq~$wq~wuq~(Geom2|=-C6?$$@$? BezierCurve$?$@$? AssocAttrib VectorInt$ xwxsrcom.femlab.geom.Rectȉ,E6DlxDlyLlxExprt!Lcom/femlab/geom/Geom$Expression;LlyExprq~-xrcom.femlab.geom.Prim2Z~DrotLbaseq~Lconstrq~[post[D[posExprt"[Lcom/femlab/geom/Geom$Expression;LrotExprq~-xq~$wq~w|uq~(|Geom2|=-C6?@@@@$@$@???? BezierCurve?@@? BezierCurve@@?@@$@? BezierCurve@@$@?$@? BezierCurve$@?? AssocAttrib VectorInt, xwxwtcornerur[D>cZxpur"[Lcom.femlab.geom.Geom$Expression;\2YUxpsrcom.femlab.geom.Geom$Expression OpI_hDnumScaleLexprStrq~Lthis$0tLcom/femlab/geom/Geom;xpwpw?xsq~8wpw?xwsq~8wpw?FR9xxw @@sq~8wpw?xw@sq~8wpw?xxsq~$wq~wuq~(Geom2|=-C6?$$@$@$@@@@@@$@ ?????????  BezierCurve?$? BezierCurve?@@? BezierCurve$?߼? BezierCurve$?$@$? BezierCurve߼?$@? BezierCurve?@@? BezierCurve$@?@@$@? BezierCurve@@?@@#=? BezierCurve@@#=?@@$@? AssocAttrib VectorInt@  xwxsrcom.femlab.api.client.MFileInfo3$$LfemNameq~[historyt[Ljava/lang/String;[mfileTagsAndTypest[[Ljava/lang/String;[ resetHistoryq~B[ storedNamesq~BLversionq~xpwsq~wq~q~q~ q~t COMSOL 3.4twt $Name: $t$Date: 2007/10/10 16:07:51 $xuq~ t~% COMSOL Multiphysics Model M-file % Generated by COMSOL 3.4 (COMSOL 3.4.0.248, $Date: 2007/10/10 16:07:51 $) flclear fem % COMSOL version clear vrsn vrsn.name = 'COMSOL 3.4'; vrsn.ext = ''; vrsn.major = 0; vrsn.build = 248; vrsn.rcs = '$Name: $'; vrsn.date = '$Date: 2007/10/10 16:07:51 $'; fem.version = vrsn; % Geometry g1=rect2('400','400','base','corner','pos',{'0','-400'},'rot','0'); parr={point2(0,-10)}; g2=geomcoerce('point',parr); % Analyzed geometry clear p s p.objs={g2}; p.name={'PT1'}; p.tags={'g2'}; s.objs={g1}; s.name={'R1'}; s.tags={'g1'}; fem.draw=struct('p',p,'s',s); fem.geom=geomcsg(fem); parr={point2(1,-10)}; g3=geomcoerce('point',parr); carr={curve2([0,1],[-10,-10],[1,1])}; g4=geomcoerce('curve',carr); g5=curve2([1,1],[0,-10]); % Constants fem.const = {'rho_f','1000[kg/m^3]', ... 'rho_s','2750[kg/m^3]', ... 'alpha_b','1'}; % Constants fem.const = {'rho_f','1000[kg/m^3]', ... 'rho_s','2750[kg/m^3]', ... 'alpha_b','1'}; % Geometry % Analyzed geometry clear p c s p.objs={g2,g3}; p.name={'PT1','PT2'}; p.tags={'g2','g3'}; c.objs={g4,g5}; c.name={'B1','B2'}; c.tags={'g4','g5'}; s.objs={g1}; s.name={'R1'}; s.tags={'g1'}; fem.draw=struct('p',p,'c',c,'s',s); fem.geom=geomcsg(fem); % Constants fem.const = {'rho_f','1000[kg/m^3]', ... 'rho_s','2750[kg/m^3]', ... 'alpha_b','1'}; % Constants fem.const = {'rho_f','1000[kg/m^3]', ... 'rho_s','2750[kg/m^3]', ... 'alpha_b','1'}; % Constants fem.const = {'rho_f','1000[kg/m^3]', ... 'rho_s','1920[kg/m^3]', ... 'alpha_b','1'}; % Constants fem.const = {'rho_f','1000[kg/m^3]', ... 'rho_s','1920[kg/m^3]', ... 'alpha_b','1'}; % Geometry % Analyzed geometry clear p c s p.objs={g2,g3}; p.name={'PT1','PT2'}; p.tags={'g2','g3'}; c.objs={g4,g5}; c.name={'B1','B2'}; c.tags={'g4','g5'}; s.objs={g1}; s.name={'R1'}; s.tags={'g1'}; fem.draw=struct('p',p,'c',c,'s',s); fem.geom=geomcsg(fem); % COMSOL Multiphysics Model M-file % Generated by COMSOL 3.4 (COMSOL 3.4.0.248, $Date: 2007/10/10 16:07:51 $) % Geometry g6=rect2(280,400,'base','corner','pos',[0,-400]); g7=rect2(220,400,'base','corner','pos',[0,-400]); % Analyzed geometry clear p c s p.objs={g2,g3}; p.name={'PT1','PT2'}; p.tags={'g2','g3'}; c.objs={g5,g4}; c.name={'B2','B1'}; c.tags={'g5','g4'}; s.objs={g7}; s.name={'R1'}; s.tags={'g7'}; fem.draw=struct('p',p,'c',c,'s',s); fem.geom=geomcsg(fem); g8=rect2(210,220,'base','corner','pos',[0,-220]); % Analyzed geometry clear p c s p.objs={g2,g3}; p.name={'PT1','PT2'}; p.tags={'g2','g3'}; c.objs={g5,g4}; c.name={'B2','B1'}; c.tags={'g5','g4'}; s.objs={g8}; s.name={'R1'}; s.tags={'g8'}; fem.draw=struct('p',p,'c',c,'s',s); fem.geom=geomcsg(fem); g9=rect2(140,220,'base','corner','pos',[0,-220]); % Analyzed geometry clear p c s p.objs={g2,g3}; p.name={'PT1','PT2'}; p.tags={'g2','g3'}; c.objs={g5,g4}; c.name={'B2','B1'}; c.tags={'g5','g4'}; s.objs={g9}; s.name={'R1'}; s.tags={'g9'}; fem.draw=struct('p',p,'c',c,'s',s); fem.geom=geomcsg(fem); % COMSOL Multiphysics Model M-file % Generated by COMSOL 3.4 (COMSOL 3.4.0.248, $Date: 2007/10/10 16:07:51 $) % Geometry g1=rect2(110,130,'base','corner','pos',[0,-130]); % Analyzed geometry clear p c s p.objs={g2,g3}; p.name={'PT1','PT2'}; p.tags={'g2','g3'}; c.objs={g5,g4}; c.name={'B2','B1'}; c.tags={'g5','g4'}; s.objs={g1}; s.name={'R1'}; s.tags={'g1'}; fem.draw=struct('p',p,'c',c,'s',s); fem.geom=geomcsg(fem); % COMSOL Multiphysics Model M-file % Generated by COMSOL 3.4 (COMSOL 3.4.0.248, $Date: 2007/10/10 16:07:51 $) % Geometry g6=rect2(100,110,'base','corner','pos',[0,-110]); % Analyzed geometry clear p c s p.objs={g2,g3}; p.name={'PT1','PT2'}; p.tags={'g2','g3'}; c.objs={g5,g4}; c.name={'B2','B1'}; c.tags={'g5','g4'}; s.objs={g6}; s.name={'R1'}; s.tags={'g6'}; fem.draw=struct('p',p,'c',c,'s',s); fem.geom=geomcsg(fem); g7=rect2(95,110,'base','corner','pos',[0,-110]); % Analyzed geometry clear p c s p.objs={g2,g3}; p.name={'PT1','PT2'}; p.tags={'g2','g3'}; c.objs={g5,g4}; c.name={'B2','B1'}; c.tags={'g5','g4'}; s.objs={g7}; s.name={'R1'}; s.tags={'g7'}; fem.draw=struct('p',p,'c',c,'s',s); fem.geom=geomcsg(fem); g8=rect2(95,105,'base','corner','pos',[0,-105]); % Analyzed geometry clear p c s p.objs={g2,g3}; p.name={'PT1','PT2'}; p.tags={'g2','g3'}; c.objs={g5,g4}; c.name={'B2','B1'}; c.tags={'g5','g4'}; s.objs={g8}; s.name={'R1'}; s.tags={'g8'}; fem.draw=struct('p',p,'c',c,'s',s); fem.geom=geomcsg(fem); g9=rect2(100,105,'base','corner','pos',[0,-105]); % Analyzed geometry clear p c s p.objs={g2,g3}; p.name={'PT1','PT2'}; p.tags={'g2','g3'}; c.objs={g5,g4}; c.name={'B2','B1'}; c.tags={'g5','g4'}; s.objs={g9}; s.name={'R1'}; s.tags={'g9'}; fem.draw=struct('p',p,'c',c,'s',s); fem.geom=geomcsg(fem); g10=rect2(97.5,107.5,'base','corner','pos',[0,-107.5]); % Analyzed geometry clear p c s p.objs={g2,g3}; p.name={'PT1','PT2'}; p.tags={'g2','g3'}; c.objs={g5,g4}; c.name={'B2','B1'}; c.tags={'g5','g4'}; s.objs={g10}; s.name={'R1'}; s.tags={'g10'}; fem.draw=struct('p',p,'c',c,'s',s); fem.geom=geomcsg(fem); g11=rect2(105,107.5,'base','corner','pos',[0,-107.5]); g12=curve2([0,0.1],[-10,-10]); g13=curve2([0.1,1],[-10,-10]); g14=curve2([0,2],[-12,-12]); g15=curve2([0,2],[-12,0]); g16=curve2([2,2],[-12,0]); % Analyzed geometry clear p c s p.objs={g2,g3}; p.name={'PT1','PT2'}; p.tags={'g2','g3'}; c.objs={g5,g12,g13,g14,g16}; c.name={'B2','B1','B3','B4','B5'}; c.tags={'g5','g12','g13','g14','g16'}; s.objs={g11}; s.name={'R1'}; s.tags={'g11'}; fem.draw=struct('p',p,'c',c,'s',s); fem.geom=geomcsg(fem); g17=rect2(105,105,'base','corner','pos',[0,-105]); % Analyzed geometry clear p c s p.objs={g2,g3}; p.name={'PT1','PT2'}; p.tags={'g2','g3'}; c.objs={g5,g12,g13,g14,g16}; c.name={'B2','B1','B3','B4','B5'}; c.tags={'g5','g12','g13','g14','g16'}; s.objs={g17}; s.name={'R1'}; s.tags={'g17'}; fem.draw=struct('p',p,'c',c,'s',s); fem.geom=geomcsg(fem); g18=rect2(105,110,'base','corner','pos',[0,-110]); % Analyzed geometry clear p c s p.objs={g2,g3}; p.name={'PT1','PT2'}; p.tags={'g2','g3'}; c.objs={g5,g12,g13,g14,g16}; c.name={'B2','B1','B3','B4','B5'}; c.tags={'g5','g12','g13','g14','g16'}; s.objs={g18}; s.name={'R1'}; s.tags={'g18'}; fem.draw=struct('p',p,'c',c,'s',s); fem.geom=geomcsg(fem); g19=rect2(107.5,110,'base','corner','pos',[0,-110]); % Analyzed geometry clear p c s p.objs={g2,g3}; p.name={'PT1','PT2'}; p.tags={'g2','g3'}; c.objs={g5,g12,g13,g14,g16}; c.name={'B2','B1','B3','B4','B5'}; c.tags={'g5','g12','g13','g14','g16'}; s.objs={g19}; s.name={'R1'}; s.tags={'g19'}; fem.draw=struct('p',p,'c',c,'s',s); fem.geom=geomcsg(fem); g20=rect2(107.5,102.5,'base','corner','pos',[0,-102.5]); % Analyzed geometry clear p c s p.objs={g2,g3}; p.name={'PT1','PT2'}; p.tags={'g2','g3'}; c.objs={g5,g12,g13,g14,g16}; c.name={'B2','B1','B3','B4','B5'}; c.tags={'g5','g12','g13','g14','g16'}; s.objs={g20}; s.name={'R1'}; s.tags={'g20'}; fem.draw=struct('p',p,'c',c,'s',s); fem.geom=geomcsg(fem); g21=rect2(102.5,102.5,'base','corner','pos',[0,-102.5]); % Analyzed geometry clear p c s p.objs={g2,g3}; p.name={'PT1','PT2'}; p.tags={'g2','g3'}; c.objs={g5,g12,g13,g14,g16}; c.name={'B2','B1','B3','B4','B5'}; c.tags={'g5','g12','g13','g14','g16'}; s.objs={g21}; s.name={'R1'}; s.tags={'g21'}; fem.draw=struct('p',p,'c',c,'s',s); fem.geom=geomcsg(fem); g22=rect2(100,100,'base','corner','pos',[0,-100]); % Analyzed geometry clear p c s p.objs={g2,g3}; p.name={'PT1','PT2'}; p.tags={'g2','g3'}; c.objs={g5,g12,g13,g14,g16}; c.name={'B2','B1','B3','B4','B5'}; c.tags={'g5','g12','g13','g14','g16'}; s.objs={g22}; s.name={'R1'}; s.tags={'g22'}; fem.draw=struct('p',p,'c',c,'s',s); fem.geom=geomcsg(fem); g23=rect2(95,100,'base','corner','pos',[0,-100]); % Analyzed geometry clear p c s p.objs={g2,g3}; p.name={'PT1','PT2'}; p.tags={'g2','g3'}; c.objs={g5,g12,g13,g14,g16}; c.name={'B2','B1','B3','B4','B5'}; c.tags={'g5','g12','g13','g14','g16'}; s.objs={g23}; s.name={'R1'}; s.tags={'g23'}; fem.draw=struct('p',p,'c',c,'s',s); fem.geom=geomcsg(fem); g24=rect2(100,100,'base','corner','pos',[0,-100]); % Analyzed geometry clear p c s p.objs={g2,g3}; p.name={'PT1','PT2'}; p.tags={'g2','g3'}; c.objs={g5,g12,g13,g14,g16}; c.name={'B2','B1','B3','B4','B5'}; c.tags={'g5','g12','g13','g14','g16'}; s.objs={g24}; s.name={'R1'}; s.tags={'g24'}; fem.draw=struct('p',p,'c',c,'s',s); fem.geom=geomcsg(fem); g25=rect2(102.5,105,'base','corner','pos',[0,-105]); g26=rect2(105,105,'base','corner','pos',[0,-105]); % Analyzed geometry clear p c s p.objs={g2,g3}; p.name={'PT1','PT2'}; p.tags={'g2','g3'}; c.objs={g5,g12,g13,g14,g16}; c.name={'B2','B1','B3','B4','B5'}; c.tags={'g5','g12','g13','g14','g16'}; s.objs={g26}; s.name={'R1'}; s.tags={'g26'}; fem.draw=struct('p',p,'c',c,'s',s); fem.geom=geomcsg(fem); g27=rect2(100,100,'base','corner','pos',[0,-100]); % Analyzed geometry clear p c s p.objs={g2,g3}; p.name={'PT1','PT2'}; p.tags={'g2','g3'}; c.objs={g5,g12,g13,g14,g16}; c.name={'B2','B1','B3','B4','B5'}; c.tags={'g5','g12','g13','g14','g16'}; s.objs={g27}; s.name={'R1'}; s.tags={'g27'}; fem.draw=struct('p',p,'c',c,'s',s); fem.geom=geomcsg(fem); % COMSOL Multiphysics Model M-file % Generated by COMSOL 3.4 (COMSOL 3.4.0.248, $Date: 2007/10/10 16:07:51 $) % Constants fem.const = {'rho_f','1000[kg/m^3]', ... 'rho_s','1920[kg/m^3]', ... 'alpha_b','1'}; % COMSOL Multiphysics Model M-file % Generated by COMSOL 3.4 (COMSOL 3.4.0.248, $Date: 2007/10/10 16:07:51 $) % Constants fem.const = {'rho_f','1000[kg/m^3]', ... 'rho_s','2780[kg/m^3]', ... 'alpha_b','1'}; % Constants fem.const = {'rho_f','1000[kg/m^3]', ... 'rho_s','2780[kg/m^3]', ... 'alpha_b','1'}; % COMSOL Multiphysics Model M-file % Generated by COMSOL 3.4 (COMSOL 3.4.0.248, $Date: 2007/10/10 16:07:51 $) % Geometry g1=rect2(100,97.5,'base','corner','pos',[0,-97.5]); % Analyzed geometry clear p c s p.objs={g2,g3}; p.name={'PT1','PT2'}; p.tags={'g2','g3'}; c.objs={g16,g12,g13,g5,g14}; c.name={'B5','B1','B3','B2','B4'}; c.tags={'g16','g12','g13','g5','g14'}; s.objs={g1}; s.name={'R1'}; s.tags={'g1'}; fem.draw=struct('p',p,'c',c,'s',s); fem.geom=geomcsg(fem); g4=rect2(100,110,'base','corner','pos',[0,-110]); % Analyzed geometry clear p c s p.objs={g2,g3}; p.name={'PT1','PT2'}; p.tags={'g2','g3'}; c.objs={g16,g12,g13,g5,g14}; c.name={'B5','B1','B3','B2','B4'}; c.tags={'g16','g12','g13','g5','g14'}; s.objs={g4}; s.name={'R1'}; s.tags={'g4'}; fem.draw=struct('p',p,'c',c,'s',s); fem.geom=geomcsg(fem); g6=rect2(102.5,115,'base','corner','pos',[0,-115]); % Analyzed geometry clear p c s p.objs={g2,g3}; p.name={'PT1','PT2'}; p.tags={'g2','g3'}; c.objs={g16,g12,g13,g5,g14}; c.name={'B5','B1','B3','B2','B4'}; c.tags={'g16','g12','g13','g5','g14'}; s.objs={g6}; s.name={'R1'}; s.tags={'g6'}; fem.draw=struct('p',p,'c',c,'s',s); fem.geom=geomcsg(fem); g7=rect2(105,115,'base','corner','pos',[0,-115]); g8=rect2(110,115,'base','corner','pos',[0,-115]); % Analyzed geometry clear p c s p.objs={g2,g3}; p.name={'PT1','PT2'}; p.tags={'g2','g3'}; c.objs={g16,g12,g13,g5,g14}; c.name={'B5','B1','B3','B2','B4'}; c.tags={'g16','g12','g13','g5','g14'}; s.objs={g8}; s.name={'R1'}; s.tags={'g8'}; fem.draw=struct('p',p,'c',c,'s',s); fem.geom=geomcsg(fem); g9=rect2(105,115,'base','corner','pos',[0,-115]); % Analyzed geometry clear p c s p.objs={g2,g3}; p.name={'PT1','PT2'}; p.tags={'g2','g3'}; c.objs={g16,g12,g13,g5,g14}; c.name={'B5','B1','B3','B2','B4'}; c.tags={'g16','g12','g13','g5','g14'}; s.objs={g9}; s.name={'R1'}; s.tags={'g9'}; fem.draw=struct('p',p,'c',c,'s',s); fem.geom=geomcsg(fem); g10=rect2(105,120,'base','corner','pos',[0,-120]); % Analyzed geometry clear p c s p.objs={g2,g3}; p.name={'PT1','PT2'}; p.tags={'g2','g3'}; c.objs={g16,g12,g13,g5,g14}; c.name={'B5','B1','B3','B2','B4'}; c.tags={'g16','g12','g13','g5','g14'}; s.objs={g10}; s.name={'R1'}; s.tags={'g10'}; fem.draw=struct('p',p,'c',c,'s',s); fem.geom=geomcsg(fem); g11=rect2(102.5,120,'base','corner','pos',[0,-120]); % Analyzed geometry clear p c s p.objs={g2,g3}; p.name={'PT1','PT2'}; p.tags={'g2','g3'}; c.objs={g16,g12,g13,g5,g14}; c.name={'B5','B1','B3','B2','B4'}; c.tags={'g16','g12','g13','g5','g14'}; s.objs={g11}; s.name={'R1'}; s.tags={'g11'}; fem.draw=struct('p',p,'c',c,'s',s); fem.geom=geomcsg(fem); g15=rect2(102.5,125,'base','corner','pos',[0,-125]); % Analyzed geometry clear p c s p.objs={g2,g3}; p.name={'PT1','PT2'}; p.tags={'g2','g3'}; c.objs={g16,g12,g13,g5,g14}; c.name={'B5','B1','B3','B2','B4'}; c.tags={'g16','g12','g13','g5','g14'}; s.objs={g15}; s.name={'R1'}; s.tags={'g15'}; fem.draw=struct('p',p,'c',c,'s',s); fem.geom=geomcsg(fem); g17=rect2(102.5,150,'base','corner','pos',[0,-150]); % Analyzed geometry clear p c s p.objs={g2,g3}; p.name={'PT1','PT2'}; p.tags={'g2','g3'}; c.objs={g16,g12,g13,g5,g14}; c.name={'B5','B1','B3','B2','B4'}; c.tags={'g16','g12','g13','g5','g14'}; s.objs={g17}; s.name={'R1'}; s.tags={'g17'}; fem.draw=struct('p',p,'c',c,'s',s); fem.geom=geomcsg(fem); % COMSOL Multiphysics Model M-file % Generated by COMSOL 3.4 (COMSOL 3.4.0.248, $Date: 2007/10/10 16:07:51 $) % Geometry g1=curve2([0,1],[-9.5,-9.5]); g4=curve2([1,1],[-9.5,-10.5]); g6=curve2([0,1],[-10.5,-10.5]); g7=curve2([0,5],[-7.5,-7.5]); g8=curve2([5,5],[-7.5,-12.5]); g9=curve2([0,5],[-12.5,-12.5]); g11=curve2([0,10],[-5,-5]); g15=curve2([10,10],[-5,-15]); g18=curve2([0,10],[-15,-15]); g19=curve2([0,15],[-2.5,-2.5]); g20=curve2([15,15],[-2.5,-17.5]); g21=curve2([0,15],[-17.5,-17.5]); % Analyzed geometry clear p c s p.objs={g2}; p.name={'PT1'}; p.tags={'g2'}; c.objs={g12,g1,g4,g6,g7,g8,g9,g11,g15,g18,g19,g20,g21}; c.name={'B1','B2','B3','B4','B5','B6','B7','B8','B9','B10','B11', ... 'B12','B13'}; c.tags={'g12','g1','g4','g6','g7','g8','g9','g11','g15','g18','g19', ... 'g20','g21'}; s.objs={g17}; s.name={'R1'}; s.tags={'g17'}; fem.draw=struct('p',p,'c',c,'s',s); fem.geom=geomcsg(fem); g22=rect2(95,105,'base','corner','pos',[0,-105]); % Analyzed geometry clear p c s p.objs={g2}; p.name={'PT1'}; p.tags={'g2'}; c.objs={g12,g1,g4,g6,g7,g8,g9,g11,g15,g18,g19,g20,g21}; c.name={'B1','B2','B3','B4','B5','B6','B7','B8','B9','B10','B11', ... 'B12','B13'}; c.tags={'g12','g1','g4','g6','g7','g8','g9','g11','g15','g18','g19', ... 'g20','g21'}; s.objs={g22}; s.name={'R1'}; s.tags={'g22'}; fem.draw=struct('p',p,'c',c,'s',s); fem.geom=geomcsg(fem); g23=rect2(90,107.5,'base','corner','pos',[0,-107.5]); % Analyzed geometry clear p c s p.objs={g2}; p.name={'PT1'}; p.tags={'g2'}; c.objs={g12,g1,g4,g6,g7,g8,g9,g11,g15,g18,g19,g20,g21}; c.name={'B1','B2','B3','B4','B5','B6','B7','B8','B9','B10','B11', ... 'B12','B13'}; c.tags={'g12','g1','g4','g6','g7','g8','g9','g11','g15','g18','g19', ... 'g20','g21'}; s.objs={g23}; s.name={'R1'}; s.tags={'g23'}; fem.draw=struct('p',p,'c',c,'s',s); fem.geom=geomcsg(fem); g24=rect2(87.5,107.5,'base','corner','pos',[0,-107.5]); % Analyzed geometry clear p c s p.objs={g2}; p.name={'PT1'}; p.tags={'g2'}; c.objs={g12,g1,g4,g6,g7,g8,g9,g11,g15,g18,g19,g20,g21}; c.name={'B1','B2','B3','B4','B5','B6','B7','B8','B9','B10','B11', ... 'B12','B13'}; c.tags={'g12','g1','g4','g6','g7','g8','g9','g11','g15','g18','g19', ... 'g20','g21'}; s.objs={g24}; s.name={'R1'}; s.tags={'g24'}; fem.draw=struct('p',p,'c',c,'s',s); fem.geom=geomcsg(fem); g25=rect2(92.5,107.5,'base','corner','pos',[0,-107.5]); % Analyzed geometry clear p c s p.objs={g2}; p.name={'PT1'}; p.tags={'g2'}; c.objs={g12,g1,g4,g6,g7,g8,g9,g11,g15,g18,g19,g20,g21}; c.name={'B1','B2','B3','B4','B5','B6','B7','B8','B9','B10','B11', ... 'B12','B13'}; c.tags={'g12','g1','g4','g6','g7','g8','g9','g11','g15','g18','g19', ... 'g20','g21'}; s.objs={g25}; s.name={'R1'}; s.tags={'g25'}; fem.draw=struct('p',p,'c',c,'s',s); fem.geom=geomcsg(fem); g26=rect2(87.5,107.5,'base','corner','pos',[0,-107.5]); g27=rect2(85,107.5,'base','corner','pos',[0,-107.5]); % Analyzed geometry clear p c s p.objs={g2}; p.name={'PT1'}; p.tags={'g2'}; c.objs={g12,g1,g4,g6,g7,g8,g9,g11,g15,g18,g19,g20,g21}; c.name={'B1','B2','B3','B4','B5','B6','B7','B8','B9','B10','B11', ... 'B12','B13'}; c.tags={'g12','g1','g4','g6','g7','g8','g9','g11','g15','g18','g19', ... 'g20','g21'}; s.objs={g27}; s.name={'R1'}; s.tags={'g27'}; fem.draw=struct('p',p,'c',c,'s',s); fem.geom=geomcsg(fem); g28=rect2(85,105,'base','corner','pos',[0,-105]); % Analyzed geometry clear p c s p.objs={g2}; p.name={'PT1'}; p.tags={'g2'}; c.objs={g12,g1,g4,g6,g7,g8,g9,g11,g15,g18,g19,g20,g21}; c.name={'B1','B2','B3','B4','B5','B6','B7','B8','B9','B10','B11', ... 'B12','B13'}; c.tags={'g12','g1','g4','g6','g7','g8','g9','g11','g15','g18','g19', ... 'g20','g21'}; s.objs={g28}; s.name={'R1'}; s.tags={'g28'}; fem.draw=struct('p',p,'c',c,'s',s); fem.geom=geomcsg(fem); g29=rect2(85,110,'base','corner','pos',[0,-110]); % Analyzed geometry clear p c s p.objs={g2}; p.name={'PT1'}; p.tags={'g2'}; c.objs={g12,g1,g4,g6,g7,g8,g9,g11,g15,g18,g19,g20,g21}; c.name={'B1','B2','B3','B4','B5','B6','B7','B8','B9','B10','B11', ... 'B12','B13'}; c.tags={'g12','g1','g4','g6','g7','g8','g9','g11','g15','g18','g19', ... 'g20','g21'}; s.objs={g29}; s.name={'R1'}; s.tags={'g29'}; fem.draw=struct('p',p,'c',c,'s',s); fem.geom=geomcsg(fem); g30=rect2(85,112.5,'base','corner','pos',[0,-112.5]); % Analyzed geometry clear p c s p.objs={g2}; p.name={'PT1'}; p.tags={'g2'}; c.objs={g12,g1,g4,g6,g7,g8,g9,g11,g15,g18,g19,g20,g21}; c.name={'B1','B2','B3','B4','B5','B6','B7','B8','B9','B10','B11', ... 'B12','B13'}; c.tags={'g12','g1','g4','g6','g7','g8','g9','g11','g15','g18','g19', ... 'g20','g21'}; s.objs={g30}; s.name={'R1'}; s.tags={'g30'}; fem.draw=struct('p',p,'c',c,'s',s); fem.geom=geomcsg(fem); g31=rect2(85,117.5,'base','corner','pos',[0,-117.5]); % Analyzed geometry clear p c s p.objs={g2}; p.name={'PT1'}; p.tags={'g2'}; c.objs={g12,g1,g4,g6,g7,g8,g9,g11,g15,g18,g19,g20,g21}; c.name={'B1','B2','B3','B4','B5','B6','B7','B8','B9','B10','B11', ... 'B12','B13'}; c.tags={'g12','g1','g4','g6','g7','g8','g9','g11','g15','g18','g19', ... 'g20','g21'}; s.objs={g31}; s.name={'R1'}; s.tags={'g31'}; fem.draw=struct('p',p,'c',c,'s',s); fem.geom=geomcsg(fem); g32=rect2(90,117.5,'base','corner','pos',[0,-117.5]); % Analyzed geometry clear p c s p.objs={g2}; p.name={'PT1'}; p.tags={'g2'}; c.objs={g12,g1,g4,g6,g7,g8,g9,g11,g15,g18,g19,g20,g21}; c.name={'B1','B2','B3','B4','B5','B6','B7','B8','B9','B10','B11', ... 'B12','B13'}; c.tags={'g12','g1','g4','g6','g7','g8','g9','g11','g15','g18','g19', ... 'g20','g21'}; s.objs={g32}; s.name={'R1'}; s.tags={'g32'}; fem.draw=struct('p',p,'c',c,'s',s); fem.geom=geomcsg(fem); g33=rect2(90,95,'base','corner','pos',[0,-95]); % Analyzed geometry clear p c s p.objs={g2}; p.name={'PT1'}; p.tags={'g2'}; c.objs={g12,g1,g4,g6,g7,g8,g9,g11,g15,g18,g19,g20,g21}; c.name={'B1','B2','B3','B4','B5','B6','B7','B8','B9','B10','B11', ... 'B12','B13'}; c.tags={'g12','g1','g4','g6','g7','g8','g9','g11','g15','g18','g19', ... 'g20','g21'}; s.objs={g33}; s.name={'R1'}; s.tags={'g33'}; fem.draw=struct('p',p,'c',c,'s',s); fem.geom=geomcsg(fem); g34=rect2(90,117.5,'base','corner','pos',[0,-117.5]); g35=rect2(85,117.5,'base','corner','pos',[0,-117.5]); % Analyzed geometry clear p c s p.objs={g2}; p.name={'PT1'}; p.tags={'g2'}; c.objs={g12,g1,g4,g6,g7,g8,g9,g11,g15,g18,g19,g20,g21}; c.name={'B1','B2','B3','B4','B5','B6','B7','B8','B9','B10','B11', ... 'B12','B13'}; c.tags={'g12','g1','g4','g6','g7','g8','g9','g11','g15','g18','g19', ... 'g20','g21'}; s.objs={g35}; s.name={'R1'}; s.tags={'g35'}; fem.draw=struct('p',p,'c',c,'s',s); fem.geom=geomcsg(fem); % COMSOL Multiphysics Model M-file % Generated by COMSOL 3.4 (COMSOL 3.4.0.248, $Date: 2007/10/10 16:07:51 $) % Geometry g3=curve2([0,1],[-10,-10]); % Analyzed geometry clear c s c.objs={g9,g21,g7,g8,g11,g20,g18,g19,g15,g3}; c.name={'B7','B13','B5','B6','B8','B12','B10','B11','B9','B1'}; c.tags={'g9','g21','g7','g8','g11','g20','g18','g19','g15','g3'}; s.objs={g35}; s.name={'R1'}; s.tags={'g35'}; fem.draw=struct('c',c,'s',s); fem.geom=geomcsg(fem); % COMSOL Multiphysics Model M-file % Generated by COMSOL 3.4 (COMSOL 3.4.0.248, $Date: 2007/10/10 16:07:51 $) % Geometry g1=rect2(100,117.5,'base','corner','pos',[0,-117.5]); % Analyzed geometry clear c s c.objs={g21,g9,g7,g20,g11,g8,g18,g19,g3,g15}; c.name={'B13','B7','B5','B12','B8','B6','B10','B11','B1','B9'}; c.tags={'g21','g9','g7','g20','g11','g8','g18','g19','g3','g15'}; s.objs={g1}; s.name={'R1'}; s.tags={'g1'}; fem.draw=struct('c',c,'s',s); fem.geom=geomcsg(fem); % COMSOL Multiphysics Model M-file % Generated by COMSOL 3.4 (COMSOL 3.4.0.248, $Date: 2007/10/10 16:07:51 $) % Geometry g2=rect2(100,127.5,'base','corner','pos',[0,-117.5]); g4=curve2([0,100],[0,0]); % Analyzed geometry clear c s c.objs={g4}; c.name={'B1'}; c.tags={'g4'}; s.objs={g2}; s.name={'R1'}; s.tags={'g2'}; fem.draw=struct('c',c,'s',s); fem.geom=geomcsg(fem); g5=curve2([0,10],[-10,-10]); [g6]=geomcopy({g5}); g10=curve2([0,1],[-10,-10]); % Analyzed geometry clear c s c.objs={g4,g10}; c.name={'B1','B2'}; c.tags={'g4','g10'}; s.objs={g2}; s.name={'R1'}; s.tags={'g2'}; fem.draw=struct('c',c,'s',s); fem.geom=geomcsg(fem); % COMSOL Multiphysics Model M-file % Generated by COMSOL 3.4 (COMSOL 3.4.0.248, $Date: 2007/10/10 16:07:51 $) % Constants fem.const = {'rho_f','1000[kg/m^3]', ... 'rho_s','2780[kg/m^3]', ... 'alpha_b','1'}; % Constants fem.const = {'rho_f','1000[kg/m^3]', ... 'rho_s','2780[kg/m^3]', ... 'alpha_b','1'}; % Geometry g1=rect2(105,142.5,'base','corner','pos',[0,-132.5]); g3=rect2(200,210,'base','corner','pos',[0,-200]); % Analyzed geometry clear c s c.objs={g10,g4}; c.name={'B2','B1'}; c.tags={'g10','g4'}; s.objs={g3}; s.name={'R1'}; s.tags={'g3'}; fem.draw=struct('c',c,'s',s); fem.geom=geomcsg(fem); [g5]=geomcopy({g4}); g6=curve2([0,200],[0,0]); % Analyzed geometry clear c s c.objs={g10,g6}; c.name={'B2','B1'}; c.tags={'g10','g6'}; s.objs={g3}; s.name={'R1'}; s.tags={'g3'}; fem.draw=struct('c',c,'s',s); fem.geom=geomcsg(fem); % COMSOL Multiphysics Model M-file % Generated by COMSOL 3.4 (COMSOL 3.4.0.248, $Date: 2007/10/10 16:07:51 $) % Constants fem.const = {'rho_f','1000[kg/m^3]', ... 'rho_s','2780[kg/m^3]', ... 'alpha_b','1'}; % Constants fem.const = {'rho_f','1000[kg/m^3]', ... 'rho_s','2780[kg/m^3]', ... 'alpha_b','1'}; % COMSOL Multiphysics Model M-file % Generated by COMSOL 3.4 (COMSOL 3.4.0.248, $Date: 2007/10/10 16:07:51 $) % Geometry g1=rect2(300,210,'base','corner','pos',[0,-200]); g6=move(g6,[2.5,0]); [g2]=geomcopy({g6}); g4=curve2([0,300],[0,0]); % Analyzed geometry clear c s c.objs={g10,g4}; c.name={'B2','B1'}; c.tags={'g10','g4'}; s.objs={g1}; s.name={'R1'}; s.tags={'g1'}; fem.draw=struct('c',c,'s',s); fem.geom=geomcsg(fem); g5=rect2(600,410,'base','corner','pos',[0,-400]); [g7]=geomcopy({g4}); g8=curve2([0,600],[0,0]); % Analyzed geometry clear c s c.objs={g10,g8}; c.name={'B2','B1'}; c.tags={'g10','g8'}; s.objs={g5}; s.name={'R1'}; s.tags={'g5'}; fem.draw=struct('c',c,'s',s); fem.geom=geomcsg(fem); g9=rect2(1002.5,520,'base','corner','pos',[0,-510]); [g11]=geomcopy({g8}); g12=curve2([0,1000],[0,0]); g13=rect2(1000,520,'base','corner','pos',[0,-510]); % Analyzed geometry clear c s c.objs={g10,g12}; c.name={'B2','B1'}; c.tags={'g10','g12'}; s.objs={g13}; s.name={'R1'}; s.tags={'g13'}; fem.draw=struct('c',c,'s',s); fem.geom=geomcsg(fem); g14=rect2(1500,520,'base','corner','pos',[0,-510]); [g15]=geomcopy({g12}); g16=curve2([0,1500],[0,0]); % Analyzed geometry clear c s c.objs={g10,g16}; c.name={'B2','B1'}; c.tags={'g10','g16'}; s.objs={g14}; s.name={'R1'}; s.tags={'g14'}; fem.draw=struct('c',c,'s',s); fem.geom=geomcsg(fem); g17=rect2(1500,1017.5,'base','corner','pos',[0,-1007.5]); % Analyzed geometry clear c s c.objs={g10,g16}; c.name={'B2','B1'}; c.tags={'g10','g16'}; s.objs={g17}; s.name={'R1'}; s.tags={'g17'}; fem.draw=struct('c',c,'s',s); fem.geom=geomcsg(fem); % COMSOL Multiphysics Model M-file % Generated by COMSOL 3.4 (COMSOL 3.4.0.248, $Date: 2007/10/10 16:07:51 $) % Geometry g1=rect2(1500,510,'base','corner','pos',[0,-500]); % Analyzed geometry clear c s c.objs={g10,g16}; c.name={'B2','B1'}; c.tags={'g10','g16'}; s.objs={g1}; s.name={'R1'}; s.tags={'g1'}; fem.draw=struct('c',c,'s',s); fem.geom=geomcsg(fem); % COMSOL Multiphysics Model M-file % Generated by COMSOL 3.4 (COMSOL 3.4.0.248, $Date: 2007/10/10 16:07:51 $) % Geometry g3=rect2(1500,262.5,'base','corner','pos',[0,-252.5]); % Analyzed geometry clear c s c.objs={g10,g16}; c.name={'B2','B1'}; c.tags={'g10','g16'}; s.objs={g3}; s.name={'R1'}; s.tags={'g3'}; fem.draw=struct('c',c,'s',s); fem.geom=geomcsg(fem); g4=rect2(800,262.5,'base','corner','pos',[0,-252.5]); g5=curve2([0,800],[0,0]); % Analyzed geometry clear c s c.objs={g10,g5}; c.name={'B2','B1'}; c.tags={'g10','g5'}; s.objs={g4}; s.name={'R1'}; s.tags={'g4'}; fem.draw=struct('c',c,'s',s); fem.geom=geomcsg(fem); g6=rect2(800,510,'base','corner','pos',[0,-500]); % Analyzed geometry clear c s c.objs={g10,g5}; c.name={'B2','B1'}; c.tags={'g10','g5'}; s.objs={g6}; s.name={'R1'}; s.tags={'g6'}; fem.draw=struct('c',c,'s',s); fem.geom=geomcsg(fem); g7=rect2(800,230,'base','corner','pos',[0,-220]); g8=rect2(810,170,'base','corner','pos',[0,-160]); g9=rect2(800,170,'base','corner','pos',[0,-160]); % Analyzed geometry clear c s c.objs={g10,g5}; c.name={'B2','B1'}; c.tags={'g10','g5'}; s.objs={g9}; s.name={'R1'}; s.tags={'g9'}; fem.draw=struct('c',c,'s',s); fem.geom=geomcsg(fem); % COMSOL Multiphysics Model M-file % Generated by COMSOL 3.4 (COMSOL 3.4.0.248, $Date: 2007/10/10 16:07:51 $) % Geometry g1=rect2(1200,170,'base','corner','pos',[0,-160]); g2=curve2([0,1200],[0,0]); % Analyzed geometry clear c s c.objs={g10,g2}; c.name={'B2','B1'}; c.tags={'g10','g2'}; s.objs={g1}; s.name={'R1'}; s.tags={'g1'}; fem.draw=struct('c',c,'s',s); fem.geom=geomcsg(fem); g3=rect2(1200,472.5,'base','corner','pos',[0,-462.5]); % Analyzed geometry clear c s c.objs={g10,g2}; c.name={'B2','B1'}; c.tags={'g10','g2'}; s.objs={g3}; s.name={'R1'}; s.tags={'g3'}; fem.draw=struct('c',c,'s',s); fem.geom=geomcsg(fem); g4=rect2(1500,472.5,'base','corner','pos',[0,-462.5]); g6=curve2([0,1500],[0,0]); % Analyzed geometry clear c s c.objs={g10,g6}; c.name={'B2','B1'}; c.tags={'g10','g6'}; s.objs={g4}; s.name={'R1'}; s.tags={'g4'}; fem.draw=struct('c',c,'s',s); fem.geom=geomcsg(fem); g7=rect2(1000,472.5,'base','corner','pos',[0,-462.5]); g8=curve2([0,1000],[0,0]); % Analyzed geometry clear c s c.objs={g10,g8}; c.name={'B2','B1'}; c.tags={'g10','g8'}; s.objs={g7}; s.name={'R1'}; s.tags={'g7'}; fem.draw=struct('c',c,'s',s); fem.geom=geomcsg(fem); g7=move(g7,[0,-3.75]); g11=rect2(800,407.5,'base','corner','pos',[0,-401.25]); [g12]=geomcopy({g8}); g13=curve2([0,800],[0,0]); % Analyzed geometry clear c s c.objs={g10,g13}; c.name={'B2','B1'}; c.tags={'g10','g13'}; s.objs={g11}; s.name={'R1'}; s.tags={'g11'}; fem.draw=struct('c',c,'s',s); fem.geom=geomcsg(fem); g11=move(g11,[0,3.75]); % Analyzed geometry clear c s c.objs={g10,g13}; c.name={'B2','B1'}; c.tags={'g10','g13'}; s.objs={g11}; s.name={'R1'}; s.tags={'g11'}; fem.draw=struct('c',c,'s',s); fem.geom=geomcsg(fem); g14=rect2(1000,407.5,'base','corner','pos',[0,-397.5]); [g15]=geomcopy({g13}); g16=curve2([0,1000],[0,0]); % Analyzed geometry clear c s c.objs={g10,g16}; c.name={'B2','B1'}; c.tags={'g10','g16'}; s.objs={g14}; s.name={'R1'}; s.tags={'g14'}; fem.draw=struct('c',c,'s',s); fem.geom=geomcsg(fem); g17=rect2(800,407.5,'base','corner','pos',[0,-397.5]); g18=curve2([0,800],[0,0]); % Analyzed geometry clear c s c.objs={g10,g18}; c.name={'B2','B1'}; c.tags={'g10','g18'}; s.objs={g17}; s.name={'R1'}; s.tags={'g17'}; fem.draw=struct('c',c,'s',s); fem.geom=geomcsg(fem); [g19]=geomcopy({g10}); g20=curve2([0,10],[-10,-10]); % Analyzed geometry clear c s c.objs={g18,g20}; c.name={'B1','B2'}; c.tags={'g18','g20'}; s.objs={g17}; s.name={'R1'}; s.tags={'g17'}; fem.draw=struct('c',c,'s',s); fem.geom=geomcsg(fem); g21=rect2(1000,512.5,'base','corner','pos',[0,-502.5]); g22=curve2([0,1000],[0,0]); % Analyzed geometry clear c s c.objs={g20,g22}; c.name={'B2','B1'}; c.tags={'g20','g22'}; s.objs={g21}; s.name={'R1'}; s.tags={'g21'}; fem.draw=struct('c',c,'s',s); fem.geom=geomcsg(fem); g21=move(g21,[0,2.5]); g23=rect2(1200,512.5,'base','corner','pos',[0,-500]); [g24]=geomcopy({g22}); g25=curve2([0,1200],[0,0]); % Analyzed geometry clear c s c.objs={g20,g25}; c.name={'B2','B1'}; c.tags={'g20','g25'}; s.objs={g23}; s.name={'R1'}; s.tags={'g23'}; fem.draw=struct('c',c,'s',s); fem.geom=geomcsg(fem); g26=rect2(1200,510,'base','corner','pos',[0,-500]); g27=curve2([0,1200],[0,0]); % Analyzed geometry clear c s c.objs={g20,g27}; c.name={'B2','B1'}; c.tags={'g20','g27'}; s.objs={g26}; s.name={'R1'}; s.tags={'g26'}; fem.draw=struct('c',c,'s',s); fem.geom=geomcsg(fem); % COMSOL Multiphysics Model M-file % Generated by COMSOL 3.4 (COMSOL 3.4.0.248, $Date: 2007/10/10 16:07:51 $) % Geometry g1=rect2(1000,507.5,'base','corner','pos',[0,-497.5]); g2=rect2(1000,1210,'base','corner','pos',[0,-1200]); g5=rect2(990,1110,'base','corner','pos',[0,-1100]); g6=rect2(1000,1110,'base','corner','pos',[0,-1100]); [g7]=geomcopy({g27}); g8=curve2([0,1000],[0,0]); % Analyzed geometry clear c s c.objs={g20,g8}; c.name={'B2','B1'}; c.tags={'g20','g8'}; s.objs={g6}; s.name={'R1'}; s.tags={'g6'}; fem.draw=struct('c',c,'s',s); fem.geom=geomcsg(fem); % COMSOL Multiphysics Model M-file % Generated by COMSOL 3.4 (COMSOL 3.4.0.248, $Date: 2007/10/10 16:07:51 $) % Geometry g1=rect2(1000,1020,'base','corner','pos',[0,-1010]); % Analyzed geometry clear c s c.objs={g8,g20}; c.name={'B1','B2'}; c.tags={'g8','g20'}; s.objs={g1}; s.name={'R1'}; s.tags={'g1'}; fem.draw=struct('c',c,'s',s); fem.geom=geomcsg(fem); % COMSOL Multiphysics Model M-file % Generated by COMSOL 3.4 (COMSOL 3.4.0.248, $Date: 2007/10/10 16:07:51 $) % (Default values are not included) % Application mode 1 clear appl appl.mode.class = 'Darcy'; appl.mode.type = 'axi'; appl.module = 'ES'; appl.border = 'on'; appl.assignsuffix = '_esdl'; clear prop prop.variable='H'; prop.analysis='time'; appl.prop = prop; clear bnd bnd.N0 = {0,0,1e-10,0,0}; bnd.type = {'Na','Nz','N0','cont','H0'}; bnd.ind = [1,2,1,3,1,4,5,2,2]; appl.bnd = bnd; clear equ equ.Ks = 'K_s'; equ.rhof = 'rho_f'; equ.Qs = '-alpha_b*(er_t_smaxi+ez_t_smaxi)'; equ.S = 'S_alpha'; equ.ind = [1,1]; appl.equ = equ; appl.var = {'g','9.81'}; fem.appl{1} = appl; % Application mode 2 clear appl appl.mode.class = 'SmeAxialSolid'; appl.mode.type = 'axi'; appl.dim = {'uor2','w2','p2'}; appl.module = 'SME'; appl.shape = {'shlag(2,''uor2'')','shlag(2,''w2'')','shlag(1,''p2'')'}; appl.gporder = 4; appl.cporder = 2; appl.assignsuffix = '_smaxi'; clear prop prop.analysis='time'; clear weakconstr weakconstr.value = 'off'; weakconstr.dim = {'lm2','lm3'}; prop.weakconstr = weakconstr; appl.prop = prop; clear bnd bnd.Hr = {0,1,1}; bnd.Hz = {0,0,1}; bnd.constrcond = {'free','displacement','displacement'}; bnd.ind = [2,3,2,1,2,1,1,3,3]; appl.bnd = bnd; clear equ equ.nu = 'nu'; equ.rho = 'rho_s'; equ.ini_stress = 1; equ.mixedform = 1; equ.dampingtype = 'nodamping'; equ.E = 'E'; equ.alpha = 0; equ.Fz = 'F_f2sz'; equ.Fr = 'F_f2sr'; equ.constrcond = 'displacement'; equ.ind = [1,1]; appl.equ = equ; fem.appl{2} = appl; fem.sdim = {'r','z'}; fem.frame = {'ref'}; % Simplify expressions fem.simplify = 'off'; fem.border = 1; fem.outform = 'general'; clear units; units.basesystem = 'SI'; fem.units = units; % Subdomain settings clear equ equ.ind = [1,2]; equ.dim = {'p','uor2','w2','p2'}; % Subdomain expressions equ.expr = {'S_alpha',0, ... 'K_s',{1e-5,1e-7}, ... 'nu',0.25, ... 'E','50[MPa]', ... 'C_s2f','alpha_b*(er_t_smaxi+ez_t_smaxi)'}; fem.equ = equ; % Scalar expressions fem.expr = {'F_f2sr','-alpha_b*rho_f*g_esdl*Hr_esdl', ... 'F_f2sz','-alpha_b*rho_f*g_esdl*Hz_esdl'}; % Descriptions clear descr descr.expr= {'F_f2sr','Fluid-to-structure coupling expression, r direction','F_f2sz','Fluid-to-structure coupling expression, z direction'}; fem.descr = descr; % Descriptions descr = fem.descr; descr.const= {'rho_f','Fluid density','alpha_b','Biot-Willis coefficient','rho_s','Solids density'}; fem.descr = descr; % ODE Settings clear ode clear units; units.basesystem = 'SI'; ode.units = units; fem.ode=ode; % Multiphysics fem=multiphysics(fem); uq~ tt% COMSOL Multiphysics Model M-file % Generated by COMSOL 3.4 (COMSOL 3.4.0.248, $Date: 2007/10/10 16:07:51 $) flclear fem % COMSOL version clear vrsn vrsn.name = 'COMSOL 3.4'; vrsn.ext = ''; vrsn.major = 0; vrsn.build = 248; vrsn.rcs = '$Name: $'; vrsn.date = '$Date: 2007/10/10 16:07:51 $'; fem.version = vrsn; % Geometry g1=rect2('400','400','base','corner','pos',{'0','-400'},'rot','0'); parr={point2(0,-10)}; g2=geomcoerce('point',parr); % Analyzed geometry clear p s p.objs={g2}; p.name={'PT1'}; p.tags={'g2'}; s.objs={g1}; s.name={'R1'}; s.tags={'g1'}; fem.draw=struct('p',p,'s',s); fem.geom=geomcsg(fem); parr={point2(1,-10)}; g3=geomcoerce('point',parr); carr={curve2([0,1],[-10,-10],[1,1])}; g4=geomcoerce('curve',carr); g5=curve2([1,1],[0,-10]); % Constants fem.const = {'rho_f','1000[kg/m^3]', ... 'rho_s','2750[kg/m^3]', ... 'alpha_b','1'}; % Constants fem.const = {'rho_f','1000[kg/m^3]', ... 'rho_s','2750[kg/m^3]', ... 'alpha_b','1'}; % Geometry % Analyzed geometry clear p c s p.objs={g2,g3}; p.name={'PT1','PT2'}; p.tags={'g2','g3'}; c.objs={g4,g5}; c.name={'B1','B2'}; c.tags={'g4','g5'}; s.objs={g1}; s.name={'R1'}; s.tags={'g1'}; fem.draw=struct('p',p,'c',c,'s',s); fem.geom=geomcsg(fem); % Constants fem.const = {'rho_f','1000[kg/m^3]', ... 'rho_s','2750[kg/m^3]', ... 'alpha_b','1'}; % Constants fem.const = {'rho_f','1000[kg/m^3]', ... 'rho_s','2750[kg/m^3]', ... 'alpha_b','1'}; % Constants fem.const = {'rho_f','1000[kg/m^3]', ... 'rho_s','1920[kg/m^3]', ... 'alpha_b','1'}; % Constants fem.const = {'rho_f','1000[kg/m^3]', ... 'rho_s','1920[kg/m^3]', ... 'alpha_b','1'}; % Geometry % Analyzed geometry clear p c s p.objs={g2,g3}; p.name={'PT1','PT2'}; p.tags={'g2','g3'}; c.objs={g4,g5}; c.name={'B1','B2'}; c.tags={'g4','g5'}; s.objs={g1}; s.name={'R1'}; s.tags={'g1'}; fem.draw=struct('p',p,'c',c,'s',s); fem.geom=geomcsg(fem); % COMSOL Multiphysics Model M-file % Generated by COMSOL 3.4 (COMSOL 3.4.0.248, $Date: 2007/10/10 16:07:51 $) % Geometry g6=rect2(280,400,'base','corner','pos',[0,-400]); g7=rect2(220,400,'base','corner','pos',[0,-400]); % Analyzed geometry clear p c s p.objs={g2,g3}; p.name={'PT1','PT2'}; p.tags={'g2','g3'}; c.objs={g5,g4}; c.name={'B2','B1'}; c.tags={'g5','g4'}; s.objs={g7}; s.name={'R1'}; s.tags={'g7'}; fem.draw=struct('p',p,'c',c,'s',s); fem.geom=geomcsg(fem); g8=rect2(210,220,'base','corner','pos',[0,-220]); % Analyzed geometry clear p c s p.objs={g2,g3}; p.name={'PT1','PT2'}; p.tags={'g2','g3'}; c.objs={g5,g4}; c.name={'B2','B1'}; c.tags={'g5','g4'}; s.objs={g8}; s.name={'R1'}; s.tags={'g8'}; fem.draw=struct('p',p,'c',c,'s',s); fem.geom=geomcsg(fem); g9=rect2(140,220,'base','corner','pos',[0,-220]); % Analyzed geometry clear p c s p.objs={g2,g3}; p.name={'PT1','PT2'}; p.tags={'g2','g3'}; c.objs={g5,g4}; c.name={'B2','B1'}; c.tags={'g5','g4'}; s.objs={g9}; s.name={'R1'}; s.tags={'g9'}; fem.draw=struct('p',p,'c',c,'s',s); fem.geom=geomcsg(fem); % COMSOL Multiphysics Model M-file % Generated by COMSOL 3.4 (COMSOL 3.4.0.248, $Date: 2007/10/10 16:07:51 $) % Geometry g1=rect2(110,130,'base','corner','pos',[0,-130]); % Analyzed geometry clear p c s p.objs={g2,g3}; p.name={'PT1','PT2'}; p.tags={'g2','g3'}; c.objs={g5,g4}; c.name={'B2','B1'}; c.tags={'g5','g4'}; s.objs={g1}; s.name={'R1'}; s.tags={'g1'}; fem.draw=struct('p',p,'c',c,'s',s); fem.geom=geomcsg(fem); % COMSOL Multiphysics Model M-file % Generated by COMSOL 3.4 (COMSOL 3.4.0.248, $Date: 2007/10/10 16:07:51 $) % Geometry g6=rect2(100,110,'base','corner','pos',[0,-110]); % Analyzed geometry clear p c s p.objs={g2,g3}; p.name={'PT1','PT2'}; p.tags={'g2','g3'}; c.objs={g5,g4}; c.name={'B2','B1'}; c.tags={'g5','g4'}; s.objs={g6}; s.name={'R1'}; s.tags={'g6'}; fem.draw=struct('p',p,'c',c,'s',s); fem.geom=geomcsg(fem); g7=rect2(95,110,'base','corner','pos',[0,-110]); % Analyzed geometry clear p c s p.objs={g2,g3}; p.name={'PT1','PT2'}; p.tags={'g2','g3'}; c.objs={g5,g4}; c.name={'B2','B1'}; c.tags={'g5','g4'}; s.objs={g7}; s.name={'R1'}; s.tags={'g7'}; fem.draw=struct('p',p,'c',c,'s',s); fem.geom=geomcsg(fem); g8=rect2(95,105,'base','corner','pos',[0,-105]); % Analyzed geometry clear p c s p.objs={g2,g3}; p.name={'PT1','PT2'}; p.tags={'g2','g3'}; c.objs={g5,g4}; c.name={'B2','B1'}; c.tags={'g5','g4'}; s.objs={g8}; s.name={'R1'}; s.tags={'g8'}; fem.draw=struct('p',p,'c',c,'s',s); fem.geom=geomcsg(fem); g9=rect2(100,105,'base','corner','pos',[0,-105]); % Analyzed geometry clear p c s p.objs={g2,g3}; p.name={'PT1','PT2'}; p.tags={'g2','g3'}; c.objs={g5,g4}; c.name={'B2','B1'}; c.tags={'g5','g4'}; s.objs={g9}; s.name={'R1'}; s.tags={'g9'}; fem.draw=struct('p',p,'c',c,'s',s); fem.geom=geomcsg(fem); g10=rect2(97.5,107.5,'base','corner','pos',[0,-107.5]); % Analyzed geometry clear p c s p.objs={g2,g3}; p.name={'PT1','PT2'}; p.tags={'g2','g3'}; c.objs={g5,g4}; c.name={'B2','B1'}; c.tags={'g5','g4'}; s.objs={g10}; s.name={'R1'}; s.tags={'g10'}; fem.draw=struct('p',p,'c',c,'s',s); fem.geom=geomcsg(fem); g11=rect2(105,107.5,'base','corner','pos',[0,-107.5]); g12=curve2([0,0.1],[-10,-10]); g13=curve2([0.1,1],[-10,-10]); g14=curve2([0,2],[-12,-12]); g15=curve2([0,2],[-12,0]); g16=curve2([2,2],[-12,0]); % Analyzed geometry clear p c s p.objs={g2,g3}; p.name={'PT1','PT2'}; p.tags={'g2','g3'}; c.objs={g5,g12,g13,g14,g16}; c.name={'B2','B1','B3','B4','B5'}; c.tags={'g5','g12','g13','g14','g16'}; s.objs={g11}; s.name={'R1'}; s.tags={'g11'}; fem.draw=struct('p',p,'c',c,'s',s); fem.geom=geomcsg(fem); g17=rect2(105,105,'base','corner','pos',[0,-105]); % Analyzed geometry clear p c s p.objs={g2,g3}; p.name={'PT1','PT2'}; p.tags={'g2','g3'}; c.objs={g5,g12,g13,g14,g16}; c.name={'B2','B1','B3','B4','B5'}; c.tags={'g5','g12','g13','g14','g16'}; s.objs={g17}; s.name={'R1'}; s.tags={'g17'}; fem.draw=struct('p',p,'c',c,'s',s); fem.geom=geomcsg(fem); g18=rect2(105,110,'base','corner','pos',[0,-110]); % Analyzed geometry clear p c s p.objs={g2,g3}; p.name={'PT1','PT2'}; p.tags={'g2','g3'}; c.objs={g5,g12,g13,g14,g16}; c.name={'B2','B1','B3','B4','B5'}; c.tags={'g5','g12','g13','g14','g16'}; s.objs={g18}; s.name={'R1'}; s.tags={'g18'}; fem.draw=struct('p',p,'c',c,'s',s); fem.geom=geomcsg(fem); g19=rect2(107.5,110,'base','corner','pos',[0,-110]); % Analyzed geometry clear p c s p.objs={g2,g3}; p.name={'PT1','PT2'}; p.tags={'g2','g3'}; c.objs={g5,g12,g13,g14,g16}; c.name={'B2','B1','B3','B4','B5'}; c.tags={'g5','g12','g13','g14','g16'}; s.objs={g19}; s.name={'R1'}; s.tags={'g19'}; fem.draw=struct('p',p,'c',c,'s',s); fem.geom=geomcsg(fem); g20=rect2(107.5,102.5,'base','corner','pos',[0,-102.5]); % Analyzed geometry clear p c s p.objs={g2,g3}; p.name={'PT1','PT2'}; p.tags={'g2','g3'}; c.objs={g5,g12,g13,g14,g16}; c.name={'B2','B1','B3','B4','B5'}; c.tags={'g5','g12','g13','g14','g16'}; s.objs={g20}; s.name={'R1'}; s.tags={'g20'}; fem.draw=struct('p',p,'c',c,'s',s); fem.geom=geomcsg(fem); g21=rect2(102.5,102.5,'base','corner','pos',[0,-102.5]); % Analyzed geometry clear p c s p.objs={g2,g3}; p.name={'PT1','PT2'}; p.tags={'g2','g3'}; c.objs={g5,g12,g13,g14,g16}; c.name={'B2','B1','B3','B4','B5'}; c.tags={'g5','g12','g13','g14','g16'}; s.objs={g21}; s.name={'R1'}; s.tags={'g21'}; fem.draw=struct('p',p,'c',c,'s',s); fem.geom=geomcsg(fem); g22=rect2(100,100,'base','corner','pos',[0,-100]); % Analyzed geometry clear p c s p.objs={g2,g3}; p.name={'PT1','PT2'}; p.tags={'g2','g3'}; c.objs={g5,g12,g13,g14,g16}; c.name={'B2','B1','B3','B4','B5'}; c.tags={'g5','g12','g13','g14','g16'}; s.objs={g22}; s.name={'R1'}; s.tags={'g22'}; fem.draw=struct('p',p,'c',c,'s',s); fem.geom=geomcsg(fem); g23=rect2(95,100,'base','corner','pos',[0,-100]); % Analyzed geometry clear p c s p.objs={g2,g3}; p.name={'PT1','PT2'}; p.tags={'g2','g3'}; c.objs={g5,g12,g13,g14,g16}; c.name={'B2','B1','B3','B4','B5'}; c.tags={'g5','g12','g13','g14','g16'}; s.objs={g23}; s.name={'R1'}; s.tags={'g23'}; fem.draw=struct('p',p,'c',c,'s',s); fem.geom=geomcsg(fem); g24=rect2(100,100,'base','corner','pos',[0,-100]); % Analyzed geometry clear p c s p.objs={g2,g3}; p.name={'PT1','PT2'}; p.tags={'g2','g3'}; c.objs={g5,g12,g13,g14,g16}; c.name={'B2','B1','B3','B4','B5'}; c.tags={'g5','g12','g13','g14','g16'}; s.objs={g24}; s.name={'R1'}; s.tags={'g24'}; fem.draw=struct('p',p,'c',c,'s',s); fem.geom=geomcsg(fem); g25=rect2(102.5,105,'base','corner','pos',[0,-105]); g26=rect2(105,105,'base','corner','pos',[0,-105]); % Analyzed geometry clear p c s p.objs={g2,g3}; p.name={'PT1','PT2'}; p.tags={'g2','g3'}; c.objs={g5,g12,g13,g14,g16}; c.name={'B2','B1','B3','B4','B5'}; c.tags={'g5','g12','g13','g14','g16'}; s.objs={g26}; s.name={'R1'}; s.tags={'g26'}; fem.draw=struct('p',p,'c',c,'s',s); fem.geom=geomcsg(fem); g27=rect2(100,100,'base','corner','pos',[0,-100]); % Analyzed geometry clear p c s p.objs={g2,g3}; p.name={'PT1','PT2'}; p.tags={'g2','g3'}; c.objs={g5,g12,g13,g14,g16}; c.name={'B2','B1','B3','B4','B5'}; c.tags={'g5','g12','g13','g14','g16'}; s.objs={g27}; s.name={'R1'}; s.tags={'g27'}; fem.draw=struct('p',p,'c',c,'s',s); fem.geom=geomcsg(fem); % COMSOL Multiphysics Model M-file % Generated by COMSOL 3.4 (COMSOL 3.4.0.248, $Date: 2007/10/10 16:07:51 $) % Constants fem.const = {'rho_f','1000[kg/m^3]', ... 'rho_s','1920[kg/m^3]', ... 'alpha_b','1'}; % COMSOL Multiphysics Model M-file % Generated by COMSOL 3.4 (COMSOL 3.4.0.248, $Date: 2007/10/10 16:07:51 $) % Constants fem.const = {'rho_f','1000[kg/m^3]', ... 'rho_s','2780[kg/m^3]', ... 'alpha_b','1'}; % Constants fem.const = {'rho_f','1000[kg/m^3]', ... 'rho_s','2780[kg/m^3]', ... 'alpha_b','1'}; % COMSOL Multiphysics Model M-file % Generated by COMSOL 3.4 (COMSOL 3.4.0.248, $Date: 2007/10/10 16:07:51 $) % Geometry g1=rect2(100,97.5,'base','corner','pos',[0,-97.5]); % Analyzed geometry clear p c s p.objs={g2,g3}; p.name={'PT1','PT2'}; p.tags={'g2','g3'}; c.objs={g16,g12,g13,g5,g14}; c.name={'B5','B1','B3','B2','B4'}; c.tags={'g16','g12','g13','g5','g14'}; s.objs={g1}; s.name={'R1'}; s.tags={'g1'}; fem.draw=struct('p',p,'c',c,'s',s); fem.geom=geomcsg(fem); g4=rect2(100,110,'base','corner','pos',[0,-110]); % Analyzed geometry clear p c s p.objs={g2,g3}; p.name={'PT1','PT2'}; p.tags={'g2','g3'}; c.objs={g16,g12,g13,g5,g14}; c.name={'B5','B1','B3','B2','B4'}; c.tags={'g16','g12','g13','g5','g14'}; s.objs={g4}; s.name={'R1'}; s.tags={'g4'}; fem.draw=struct('p',p,'c',c,'s',s); fem.geom=geomcsg(fem); g6=rect2(102.5,115,'base','corner','pos',[0,-115]); % Analyzed geometry clear p c s p.objs={g2,g3}; p.name={'PT1','PT2'}; p.tags={'g2','g3'}; c.objs={g16,g12,g13,g5,g14}; c.name={'B5','B1','B3','B2','B4'}; c.tags={'g16','g12','g13','g5','g14'}; s.objs={g6}; s.name={'R1'}; s.tags={'g6'}; fem.draw=struct('p',p,'c',c,'s',s); fem.geom=geomcsg(fem); g7=rect2(105,115,'base','corner','pos',[0,-115]); g8=rect2(110,115,'base','corner','pos',[0,-115]); % Analyzed geometry clear p c s p.objs={g2,g3}; p.name={'PT1','PT2'}; p.tags={'g2','g3'}; c.objs={g16,g12,g13,g5,g14}; c.name={'B5','B1','B3','B2','B4'}; c.tags={'g16','g12','g13','g5','g14'}; s.objs={g8}; s.name={'R1'}; s.tags={'g8'}; fem.draw=struct('p',p,'c',c,'s',s); fem.geom=geomcsg(fem); g9=rect2(105,115,'base','corner','pos',[0,-115]); % Analyzed geometry clear p c s p.objs={g2,g3}; p.name={'PT1','PT2'}; p.tags={'g2','g3'}; c.objs={g16,g12,g13,g5,g14}; c.name={'B5','B1','B3','B2','B4'}; c.tags={'g16','g12','g13','g5','g14'}; s.objs={g9}; s.name={'R1'}; s.tags={'g9'}; fem.draw=struct('p',p,'c',c,'s',s); fem.geom=geomcsg(fem); g10=rect2(105,120,'base','corner','pos',[0,-120]); % Analyzed geometry clear p c s p.objs={g2,g3}; p.name={'PT1','PT2'}; p.tags={'g2','g3'}; c.objs={g16,g12,g13,g5,g14}; c.name={'B5','B1','B3','B2','B4'}; c.tags={'g16','g12','g13','g5','g14'}; s.objs={g10}; s.name={'R1'}; s.tags={'g10'}; fem.draw=struct('p',p,'c',c,'s',s); fem.geom=geomcsg(fem); g11=rect2(102.5,120,'base','corner','pos',[0,-120]); % Analyzed geometry clear p c s p.objs={g2,g3}; p.name={'PT1','PT2'}; p.tags={'g2','g3'}; c.objs={g16,g12,g13,g5,g14}; c.name={'B5','B1','B3','B2','B4'}; c.tags={'g16','g12','g13','g5','g14'}; s.objs={g11}; s.name={'R1'}; s.tags={'g11'}; fem.draw=struct('p',p,'c',c,'s',s); fem.geom=geomcsg(fem); g15=rect2(102.5,125,'base','corner','pos',[0,-125]); % Analyzed geometry clear p c s p.objs={g2,g3}; p.name={'PT1','PT2'}; p.tags={'g2','g3'}; c.objs={g16,g12,g13,g5,g14}; c.name={'B5','B1','B3','B2','B4'}; c.tags={'g16','g12','g13','g5','g14'}; s.objs={g15}; s.name={'R1'}; s.tags={'g15'}; fem.draw=struct('p',p,'c',c,'s',s); fem.geom=geomcsg(fem); g17=rect2(102.5,150,'base','corner','pos',[0,-150]); % Analyzed geometry clear p c s p.objs={g2,g3}; p.name={'PT1','PT2'}; p.tags={'g2','g3'}; c.objs={g16,g12,g13,g5,g14}; c.name={'B5','B1','B3','B2','B4'}; c.tags={'g16','g12','g13','g5','g14'}; s.objs={g17}; s.name={'R1'}; s.tags={'g17'}; fem.draw=struct('p',p,'c',c,'s',s); fem.geom=geomcsg(fem); % COMSOL Multiphysics Model M-file % Generated by COMSOL 3.4 (COMSOL 3.4.0.248, $Date: 2007/10/10 16:07:51 $) % Geometry g1=curve2([0,1],[-9.5,-9.5]); g4=curve2([1,1],[-9.5,-10.5]); g6=curve2([0,1],[-10.5,-10.5]); g7=curve2([0,5],[-7.5,-7.5]); g8=curve2([5,5],[-7.5,-12.5]); g9=curve2([0,5],[-12.5,-12.5]); g11=curve2([0,10],[-5,-5]); g15=curve2([10,10],[-5,-15]); g18=curve2([0,10],[-15,-15]); g19=curve2([0,15],[-2.5,-2.5]); g20=curve2([15,15],[-2.5,-17.5]); g21=curve2([0,15],[-17.5,-17.5]); % Analyzed geometry clear p c s p.objs={g2}; p.name={'PT1'}; p.tags={'g2'}; c.objs={g12,g1,g4,g6,g7,g8,g9,g11,g15,g18,g19,g20,g21}; c.name={'B1','B2','B3','B4','B5','B6','B7','B8','B9','B10','B11', ... 'B12','B13'}; c.tags={'g12','g1','g4','g6','g7','g8','g9','g11','g15','g18','g19', ... 'g20','g21'}; s.objs={g17}; s.name={'R1'}; s.tags={'g17'}; fem.draw=struct('p',p,'c',c,'s',s); fem.geom=geomcsg(fem); g22=rect2(95,105,'base','corner','pos',[0,-105]); % Analyzed geometry clear p c s p.objs={g2}; p.name={'PT1'}; p.tags={'g2'}; c.objs={g12,g1,g4,g6,g7,g8,g9,g11,g15,g18,g19,g20,g21}; c.name={'B1','B2','B3','B4','B5','B6','B7','B8','B9','B10','B11', ... 'B12','B13'}; c.tags={'g12','g1','g4','g6','g7','g8','g9','g11','g15','g18','g19', ... 'g20','g21'}; s.objs={g22}; s.name={'R1'}; s.tags={'g22'}; fem.draw=struct('p',p,'c',c,'s',s); fem.geom=geomcsg(fem); g23=rect2(90,107.5,'base','corner','pos',[0,-107.5]); % Analyzed geometry clear p c s p.objs={g2}; p.name={'PT1'}; p.tags={'g2'}; c.objs={g12,g1,g4,g6,g7,g8,g9,g11,g15,g18,g19,g20,g21}; c.name={'B1','B2','B3','B4','B5','B6','B7','B8','B9','B10','B11', ... 'B12','B13'}; c.tags={'g12','g1','g4','g6','g7','g8','g9','g11','g15','g18','g19', ... 'g20','g21'}; s.objs={g23}; s.name={'R1'}; s.tags={'g23'}; fem.draw=struct('p',p,'c',c,'s',s); fem.geom=geomcsg(fem); g24=rect2(87.5,107.5,'base','corner','pos',[0,-107.5]); % Analyzed geometry clear p c s p.objs={g2}; p.name={'PT1'}; p.tags={'g2'}; c.objs={g12,g1,g4,g6,g7,g8,g9,g11,g15,g18,g19,g20,g21}; c.name={'B1','B2','B3','B4','B5','B6','B7','B8','B9','B10','B11', ... 'B12','B13'}; c.tags={'g12','g1','g4','g6','g7','g8','g9','g11','g15','g18','g19', ... 'g20','g21'}; s.objs={g24}; s.name={'R1'}; s.tags={'g24'}; fem.draw=struct('p',p,'c',c,'s',s); fem.geom=geomcsg(fem); g25=rect2(92.5,107.5,'base','corner','pos',[0,-107.5]); % Analyzed geometry clear p c s p.objs={g2}; p.name={'PT1'}; p.tags={'g2'}; c.objs={g12,g1,g4,g6,g7,g8,g9,g11,g15,g18,g19,g20,g21}; c.name={'B1','B2','B3','B4','B5','B6','B7','B8','B9','B10','B11', ... 'B12','B13'}; c.tags={'g12','g1','g4','g6','g7','g8','g9','g11','g15','g18','g19', ... 'g20','g21'}; s.objs={g25}; s.name={'R1'}; s.tags={'g25'}; fem.draw=struct('p',p,'c',c,'s',s); fem.geom=geomcsg(fem); g26=rect2(87.5,107.5,'base','corner','pos',[0,-107.5]); g27=rect2(85,107.5,'base','corner','pos',[0,-107.5]); % Analyzed geometry clear p c s p.objs={g2}; p.name={'PT1'}; p.tags={'g2'}; c.objs={g12,g1,g4,g6,g7,g8,g9,g11,g15,g18,g19,g20,g21}; c.name={'B1','B2','B3','B4','B5','B6','B7','B8','B9','B10','B11', ... 'B12','B13'}; c.tags={'g12','g1','g4','g6','g7','g8','g9','g11','g15','g18','g19', ... 'g20','g21'}; s.objs={g27}; s.name={'R1'}; s.tags={'g27'}; fem.draw=struct('p',p,'c',c,'s',s); fem.geom=geomcsg(fem); g28=rect2(85,105,'base','corner','pos',[0,-105]); % Analyzed geometry clear p c s p.objs={g2}; p.name={'PT1'}; p.tags={'g2'}; c.objs={g12,g1,g4,g6,g7,g8,g9,g11,g15,g18,g19,g20,g21}; c.name={'B1','B2','B3','B4','B5','B6','B7','B8','B9','B10','B11', ... 'B12','B13'}; c.tags={'g12','g1','g4','g6','g7','g8','g9','g11','g15','g18','g19', ... 'g20','g21'}; s.objs={g28}; s.name={'R1'}; s.tags={'g28'}; fem.draw=struct('p',p,'c',c,'s',s); fem.geom=geomcsg(fem); g29=rect2(85,110,'base','corner','pos',[0,-110]); % Analyzed geometry clear p c s p.objs={g2}; p.name={'PT1'}; p.tags={'g2'}; c.objs={g12,g1,g4,g6,g7,g8,g9,g11,g15,g18,g19,g20,g21}; c.name={'B1','B2','B3','B4','B5','B6','B7','B8','B9','B10','B11', ... 'B12','B13'}; c.tags={'g12','g1','g4','g6','g7','g8','g9','g11','g15','g18','g19', ... 'g20','g21'}; s.objs={g29}; s.name={'R1'}; s.tags={'g29'}; fem.draw=struct('p',p,'c',c,'s',s); fem.geom=geomcsg(fem); g30=rect2(85,112.5,'base','corner','pos',[0,-112.5]); % Analyzed geometry clear p c s p.objs={g2}; p.name={'PT1'}; p.tags={'g2'}; c.objs={g12,g1,g4,g6,g7,g8,g9,g11,g15,g18,g19,g20,g21}; c.name={'B1','B2','B3','B4','B5','B6','B7','B8','B9','B10','B11', ... 'B12','B13'}; c.tags={'g12','g1','g4','g6','g7','g8','g9','g11','g15','g18','g19', ... 'g20','g21'}; s.objs={g30}; s.name={'R1'}; s.tags={'g30'}; fem.draw=struct('p',p,'c',c,'s',s); fem.geom=geomcsg(fem); g31=rect2(85,117.5,'base','corner','pos',[0,-117.5]); % Analyzed geometry clear p c s p.objs={g2}; p.name={'PT1'}; p.tags={'g2'}; c.objs={g12,g1,g4,g6,g7,g8,g9,g11,g15,g18,g19,g20,g21}; c.name={'B1','B2','B3','B4','B5','B6','B7','B8','B9','B10','B11', ... 'B12','B13'}; c.tags={'g12','g1','g4','g6','g7','g8','g9','g11','g15','g18','g19', ... 'g20','g21'}; s.objs={g31}; s.name={'R1'}; s.tags={'g31'}; fem.draw=struct('p',p,'c',c,'s',s); fem.geom=geomcsg(fem); g32=rect2(90,117.5,'base','corner','pos',[0,-117.5]); % Analyzed geometry clear p c s p.objs={g2}; p.name={'PT1'}; p.tags={'g2'}; c.objs={g12,g1,g4,g6,g7,g8,g9,g11,g15,g18,g19,g20,g21}; c.name={'B1','B2','B3','B4','B5','B6','B7','B8','B9','B10','B11', ... 'B12','B13'}; c.tags={'g12','g1','g4','g6','g7','g8','g9','g11','g15','g18','g19', ... 'g20','g21'}; s.objs={g32}; s.name={'R1'}; s.tags={'g32'}; fem.draw=struct('p',p,'c',c,'s',s); fem.geom=geomcsg(fem); g33=rect2(90,95,'base','corner','pos',[0,-95]); % Analyzed geometry clear p c s p.objs={g2}; p.name={'PT1'}; p.tags={'g2'}; c.objs={g12,g1,g4,g6,g7,g8,g9,g11,g15,g18,g19,g20,g21}; c.name={'B1','B2','B3','B4','B5','B6','B7','B8','B9','B10','B11', ... 'B12','B13'}; c.tags={'g12','g1','g4','g6','g7','g8','g9','g11','g15','g18','g19', ... 'g20','g21'}; s.objs={g33}; s.name={'R1'}; s.tags={'g33'}; fem.draw=struct('p',p,'c',c,'s',s); fem.geom=geomcsg(fem); g34=rect2(90,117.5,'base','corner','pos',[0,-117.5]); g35=rect2(85,117.5,'base','corner','pos',[0,-117.5]); % Analyzed geometry clear p c s p.objs={g2}; p.name={'PT1'}; p.tags={'g2'}; c.objs={g12,g1,g4,g6,g7,g8,g9,g11,g15,g18,g19,g20,g21}; c.name={'B1','B2','B3','B4','B5','B6','B7','B8','B9','B10','B11', ... 'B12','B13'}; c.tags={'g12','g1','g4','g6','g7','g8','g9','g11','g15','g18','g19', ... 'g20','g21'}; s.objs={g35}; s.name={'R1'}; s.tags={'g35'}; fem.draw=struct('p',p,'c',c,'s',s); fem.geom=geomcsg(fem); % COMSOL Multiphysics Model M-file % Generated by COMSOL 3.4 (COMSOL 3.4.0.248, $Date: 2007/10/10 16:07:51 $) % Geometry g3=curve2([0,1],[-10,-10]); % Analyzed geometry clear c s c.objs={g9,g21,g7,g8,g11,g20,g18,g19,g15,g3}; c.name={'B7','B13','B5','B6','B8','B12','B10','B11','B9','B1'}; c.tags={'g9','g21','g7','g8','g11','g20','g18','g19','g15','g3'}; s.objs={g35}; s.name={'R1'}; s.tags={'g35'}; fem.draw=struct('c',c,'s',s); fem.geom=geomcsg(fem); % COMSOL Multiphysics Model M-file % Generated by COMSOL 3.4 (COMSOL 3.4.0.248, $Date: 2007/10/10 16:07:51 $) % Geometry g1=rect2(100,117.5,'base','corner','pos',[0,-117.5]); % Analyzed geometry clear c s c.objs={g21,g9,g7,g20,g11,g8,g18,g19,g3,g15}; c.name={'B13','B7','B5','B12','B8','B6','B10','B11','B1','B9'}; c.tags={'g21','g9','g7','g20','g11','g8','g18','g19','g3','g15'}; s.objs={g1}; s.name={'R1'}; s.tags={'g1'}; fem.draw=struct('c',c,'s',s); fem.geom=geomcsg(fem); % COMSOL Multiphysics Model M-file % Generated by COMSOL 3.4 (COMSOL 3.4.0.248, $Date: 2007/10/10 16:07:51 $) % Geometry g2=rect2(100,127.5,'base','corner','pos',[0,-117.5]); g4=curve2([0,100],[0,0]); % Analyzed geometry clear c s c.objs={g4}; c.name={'B1'}; c.tags={'g4'}; s.objs={g2}; s.name={'R1'}; s.tags={'g2'}; fem.draw=struct('c',c,'s',s); fem.geom=geomcsg(fem); g5=curve2([0,10],[-10,-10]); [g6]=geomcopy({g5}); g10=curve2([0,1],[-10,-10]); % Analyzed geometry clear c s c.objs={g4,g10}; c.name={'B1','B2'}; c.tags={'g4','g10'}; s.objs={g2}; s.name={'R1'}; s.tags={'g2'}; fem.draw=struct('c',c,'s',s); fem.geom=geomcsg(fem); % COMSOL Multiphysics Model M-file % Generated by COMSOL 3.4 (COMSOL 3.4.0.248, $Date: 2007/10/10 16:07:51 $) % Constants fem.const = {'rho_f','1000[kg/m^3]', ... 'rho_s','2780[kg/m^3]', ... 'alpha_b','1'}; % Constants fem.const = {'rho_f','1000[kg/m^3]', ... 'rho_s','2780[kg/m^3]', ... 'alpha_b','1'}; % Geometry g1=rect2(105,142.5,'base','corner','pos',[0,-132.5]); g3=rect2(200,210,'base','corner','pos',[0,-200]); % Analyzed geometry clear c s c.objs={g10,g4}; c.name={'B2','B1'}; c.tags={'g10','g4'}; s.objs={g3}; s.name={'R1'}; s.tags={'g3'}; fem.draw=struct('c',c,'s',s); fem.geom=geomcsg(fem); [g5]=geomcopy({g4}); g6=curve2([0,200],[0,0]); % Analyzed geometry clear c s c.objs={g10,g6}; c.name={'B2','B1'}; c.tags={'g10','g6'}; s.objs={g3}; s.name={'R1'}; s.tags={'g3'}; fem.draw=struct('c',c,'s',s); fem.geom=geomcsg(fem); % COMSOL Multiphysics Model M-file % Generated by COMSOL 3.4 (COMSOL 3.4.0.248, $Date: 2007/10/10 16:07:51 $) % Constants fem.const = {'rho_f','1000[kg/m^3]', ... 'rho_s','2780[kg/m^3]', ... 'alpha_b','1'}; % Constants fem.const = {'rho_f','1000[kg/m^3]', ... 'rho_s','2780[kg/m^3]', ... 'alpha_b','1'}; % COMSOL Multiphysics Model M-file % Generated by COMSOL 3.4 (COMSOL 3.4.0.248, $Date: 2007/10/10 16:07:51 $) % Geometry g1=rect2(300,210,'base','corner','pos',[0,-200]); g6=move(g6,[2.5,0]); [g2]=geomcopy({g6}); g4=curve2([0,300],[0,0]); % Analyzed geometry clear c s c.objs={g10,g4}; c.name={'B2','B1'}; c.tags={'g10','g4'}; s.objs={g1}; s.name={'R1'}; s.tags={'g1'}; fem.draw=struct('c',c,'s',s); fem.geom=geomcsg(fem); g5=rect2(600,410,'base','corner','pos',[0,-400]); [g7]=geomcopy({g4}); g8=curve2([0,600],[0,0]); % Analyzed geometry clear c s c.objs={g10,g8}; c.name={'B2','B1'}; c.tags={'g10','g8'}; s.objs={g5}; s.name={'R1'}; s.tags={'g5'}; fem.draw=struct('c',c,'s',s); fem.geom=geomcsg(fem); g9=rect2(1002.5,520,'base','corner','pos',[0,-510]); [g11]=geomcopy({g8}); g12=curve2([0,1000],[0,0]); g13=rect2(1000,520,'base','corner','pos',[0,-510]); % Analyzed geometry clear c s c.objs={g10,g12}; c.name={'B2','B1'}; c.tags={'g10','g12'}; s.objs={g13}; s.name={'R1'}; s.tags={'g13'}; fem.draw=struct('c',c,'s',s); fem.geom=geomcsg(fem); g14=rect2(1500,520,'base','corner','pos',[0,-510]); [g15]=geomcopy({g12}); g16=curve2([0,1500],[0,0]); % Analyzed geometry clear c s c.objs={g10,g16}; c.name={'B2','B1'}; c.tags={'g10','g16'}; s.objs={g14}; s.name={'R1'}; s.tags={'g14'}; fem.draw=struct('c',c,'s',s); fem.geom=geomcsg(fem); g17=rect2(1500,1017.5,'base','corner','pos',[0,-1007.5]); % Analyzed geometry clear c s c.objs={g10,g16}; c.name={'B2','B1'}; c.tags={'g10','g16'}; s.objs={g17}; s.name={'R1'}; s.tags={'g17'}; fem.draw=struct('c',c,'s',s); fem.geom=geomcsg(fem); % COMSOL Multiphysics Model M-file % Generated by COMSOL 3.4 (COMSOL 3.4.0.248, $Date: 2007/10/10 16:07:51 $) % Geometry g1=rect2(1500,510,'base','corner','pos',[0,-500]); % Analyzed geometry clear c s c.objs={g10,g16}; c.name={'B2','B1'}; c.tags={'g10','g16'}; s.objs={g1}; s.name={'R1'}; s.tags={'g1'}; fem.draw=struct('c',c,'s',s); fem.geom=geomcsg(fem); % COMSOL Multiphysics Model M-file % Generated by COMSOL 3.4 (COMSOL 3.4.0.248, $Date: 2007/10/10 16:07:51 $) % Geometry g3=rect2(1500,262.5,'base','corner','pos',[0,-252.5]); % Analyzed geometry clear c s c.objs={g10,g16}; c.name={'B2','B1'}; c.tags={'g10','g16'}; s.objs={g3}; s.name={'R1'}; s.tags={'g3'}; fem.draw=struct('c',c,'s',s); fem.geom=geomcsg(fem); g4=rect2(800,262.5,'base','corner','pos',[0,-252.5]); g5=curve2([0,800],[0,0]); % Analyzed geometry clear c s c.objs={g10,g5}; c.name={'B2','B1'}; c.tags={'g10','g5'}; s.objs={g4}; s.name={'R1'}; s.tags={'g4'}; fem.draw=struct('c',c,'s',s); fem.geom=geomcsg(fem); g6=rect2(800,510,'base','corner','pos',[0,-500]); % Analyzed geometry clear c s c.objs={g10,g5}; c.name={'B2','B1'}; c.tags={'g10','g5'}; s.objs={g6}; s.name={'R1'}; s.tags={'g6'}; fem.draw=struct('c',c,'s',s); fem.geom=geomcsg(fem); g7=rect2(800,230,'base','corner','pos',[0,-220]); g8=rect2(810,170,'base','corner','pos',[0,-160]); g9=rect2(800,170,'base','corner','pos',[0,-160]); % Analyzed geometry clear c s c.objs={g10,g5}; c.name={'B2','B1'}; c.tags={'g10','g5'}; s.objs={g9}; s.name={'R1'}; s.tags={'g9'}; fem.draw=struct('c',c,'s',s); fem.geom=geomcsg(fem); % COMSOL Multiphysics Model M-file % Generated by COMSOL 3.4 (COMSOL 3.4.0.248, $Date: 2007/10/10 16:07:51 $) % Geometry g1=rect2(1200,170,'base','corner','pos',[0,-160]); g2=curve2([0,1200],[0,0]); % Analyzed geometry clear c s c.objs={g10,g2}; c.name={'B2','B1'}; c.tags={'g10','g2'}; s.objs={g1}; s.name={'R1'}; s.tags={'g1'}; fem.draw=struct('c',c,'s',s); fem.geom=geomcsg(fem); g3=rect2(1200,472.5,'base','corner','pos',[0,-462.5]); % Analyzed geometry clear c s c.objs={g10,g2}; c.name={'B2','B1'}; c.tags={'g10','g2'}; s.objs={g3}; s.name={'R1'}; s.tags={'g3'}; fem.draw=struct('c',c,'s',s); fem.geom=geomcsg(fem); g4=rect2(1500,472.5,'base','corner','pos',[0,-462.5]); g6=curve2([0,1500],[0,0]); % Analyzed geometry clear c s c.objs={g10,g6}; c.name={'B2','B1'}; c.tags={'g10','g6'}; s.objs={g4}; s.name={'R1'}; s.tags={'g4'}; fem.draw=struct('c',c,'s',s); fem.geom=geomcsg(fem); g7=rect2(1000,472.5,'base','corner','pos',[0,-462.5]); g8=curve2([0,1000],[0,0]); % Analyzed geometry clear c s c.objs={g10,g8}; c.name={'B2','B1'}; c.tags={'g10','g8'}; s.objs={g7}; s.name={'R1'}; s.tags={'g7'}; fem.draw=struct('c',c,'s',s); fem.geom=geomcsg(fem); g7=move(g7,[0,-3.75]); g11=rect2(800,407.5,'base','corner','pos',[0,-401.25]); [g12]=geomcopy({g8}); g13=curve2([0,800],[0,0]); % Analyzed geometry clear c s c.objs={g10,g13}; c.name={'B2','B1'}; c.tags={'g10','g13'}; s.objs={g11}; s.name={'R1'}; s.tags={'g11'}; fem.draw=struct('c',c,'s',s); fem.geom=geomcsg(fem); g11=move(g11,[0,3.75]); % Analyzed geometry clear c s c.objs={g10,g13}; c.name={'B2','B1'}; c.tags={'g10','g13'}; s.objs={g11}; s.name={'R1'}; s.tags={'g11'}; fem.draw=struct('c',c,'s',s); fem.geom=geomcsg(fem); g14=rect2(1000,407.5,'base','corner','pos',[0,-397.5]); [g15]=geomcopy({g13}); g16=curve2([0,1000],[0,0]); % Analyzed geometry clear c s c.objs={g10,g16}; c.name={'B2','B1'}; c.tags={'g10','g16'}; s.objs={g14}; s.name={'R1'}; s.tags={'g14'}; fem.draw=struct('c',c,'s',s); fem.geom=geomcsg(fem); g17=rect2(800,407.5,'base','corner','pos',[0,-397.5]); g18=curve2([0,800],[0,0]); % Analyzed geometry clear c s c.objs={g10,g18}; c.name={'B2','B1'}; c.tags={'g10','g18'}; s.objs={g17}; s.name={'R1'}; s.tags={'g17'}; fem.draw=struct('c',c,'s',s); fem.geom=geomcsg(fem); [g19]=geomcopy({g10}); g20=curve2([0,10],[-10,-10]); % Analyzed geometry clear c s c.objs={g18,g20}; c.name={'B1','B2'}; c.tags={'g18','g20'}; s.objs={g17}; s.name={'R1'}; s.tags={'g17'}; fem.draw=struct('c',c,'s',s); fem.geom=geomcsg(fem); g21=rect2(1000,512.5,'base','corner','pos',[0,-502.5]); g22=curve2([0,1000],[0,0]); % Analyzed geometry clear c s c.objs={g20,g22}; c.name={'B2','B1'}; c.tags={'g20','g22'}; s.objs={g21}; s.name={'R1'}; s.tags={'g21'}; fem.draw=struct('c',c,'s',s); fem.geom=geomcsg(fem); g21=move(g21,[0,2.5]); g23=rect2(1200,512.5,'base','corner','pos',[0,-500]); [g24]=geomcopy({g22}); g25=curve2([0,1200],[0,0]); % Analyzed geometry clear c s c.objs={g20,g25}; c.name={'B2','B1'}; c.tags={'g20','g25'}; s.objs={g23}; s.name={'R1'}; s.tags={'g23'}; fem.draw=struct('c',c,'s',s); fem.geom=geomcsg(fem); g26=rect2(1200,510,'base','corner','pos',[0,-500]); g27=curve2([0,1200],[0,0]); % Analyzed geometry clear c s c.objs={g20,g27}; c.name={'B2','B1'}; c.tags={'g20','g27'}; s.objs={g26}; s.name={'R1'}; s.tags={'g26'}; fem.draw=struct('c',c,'s',s); fem.geom=geomcsg(fem); % COMSOL Multiphysics Model M-file % Generated by COMSOL 3.4 (COMSOL 3.4.0.248, $Date: 2007/10/10 16:07:51 $) % Geometry g1=rect2(1000,507.5,'base','corner','pos',[0,-497.5]); g2=rect2(1000,1210,'base','corner','pos',[0,-1200]); g5=rect2(990,1110,'base','corner','pos',[0,-1100]); g6=rect2(1000,1110,'base','corner','pos',[0,-1100]); [g7]=geomcopy({g27}); g8=curve2([0,1000],[0,0]); % Analyzed geometry clear c s c.objs={g20,g8}; c.name={'B2','B1'}; c.tags={'g20','g8'}; s.objs={g6}; s.name={'R1'}; s.tags={'g6'}; fem.draw=struct('c',c,'s',s); fem.geom=geomcsg(fem); % COMSOL Multiphysics Model M-file % Generated by COMSOL 3.4 (COMSOL 3.4.0.248, $Date: 2007/10/10 16:07:51 $) % Geometry g1=rect2(1000,1020,'base','corner','pos',[0,-1010]); % Analyzed geometry clear c s c.objs={g8,g20}; c.name={'B1','B2'}; c.tags={'g8','g20'}; s.objs={g1}; s.name={'R1'}; s.tags={'g1'}; fem.draw=struct('c',c,'s',s); fem.geom=geomcsg(fem); tclear mfile clear vrsn vrsn.name = 'COMSOL 3.4'; vrsn.ext = ''; vrsn.major = 0; vrsn.build = 248; vrsn.rcs = '$Name: $'; vrsn.date = '$Date: 2007/10/10 16:07:51 $'; mfile.version=vrsn; mfile.fem='fem'; mfile.stored={'fem0','fem1'}; mfile.tags={}; mfile.types={}; x