% COMSOL Multiphysics Model M-file % Generated by COMSOL 3.5a (COMSOL 3.5.0.608, $Date: 2009/05/11 07:38:49 $) % Some geometry objects are stored in a separate file. % The name of this file is given by the variable 'flbinaryfile'. flclear xfem % COMSOL version clear vrsn vrsn.name = 'COMSOL 3.5'; vrsn.ext = 'a'; vrsn.major = 0; vrsn.build = 608; vrsn.rcs = '$Name: v35ap $'; vrsn.date = '$Date: 2009/05/11 07:38:49 $'; xfem.version = vrsn; flbinaryfile=''; % Constants xfem.const = {'rho','1000[ kg/(m^3)]'}; flclear fem % Geometry 1 clear draw g2=flbinary('g2','draw',flbinaryfile); g6=flbinary('g6','draw',flbinaryfile); g26=flbinary('g26','draw',flbinaryfile); g4=flbinary('g4','draw',flbinaryfile); draw.s.objs = {g2,g6,g26,g4}; draw.s.name = {'EXT1','EXT4','EXT2','EXT3'}; draw.s.tags = {'g2','g6','g26','g4'}; fem.draw = draw; fem.geom = geomcsg(fem); xfem.fem{1}=fem; flclear fem % Geometry 2 clear draw g1=flbinary('g1','draw',flbinaryfile); g25=flbinary('g25','draw',flbinaryfile); g3=flbinary('g3','draw',flbinaryfile); draw.s.objs = {g1,g25,g3}; draw.s.name = {'CO1','CO2','CO3'}; draw.s.tags = {'g1','g25','g3'}; fem.draw = draw; % Geometry 2 % Geometry objects clear s s.objs={g1,g25,g3}; s.name={'CO1','CO2','CO3'}; s.tags={'g1','g25','g3'}; fem.draw=struct('s',s); xfem.fem{2}=fem; % (Default values are not included) fem=xfem.fem{1}; % Application mode 1 clear appl appl.mode.class = 'PiezoSolid3'; appl.module = 'MEMS'; appl.shape = {'shlag(2,''lm1'')','shlag(2,''lm2'')','shlag(2,''lm3'')','shlag(2,''lm4'')','shlag(2,''u'')','shlag(2,''v'')','shlag(2,''w'')','shlag(2,''V'')'}; appl.gporder = {30,10,4}; appl.cporder = 2; appl.border = 'on'; appl.assignsuffix = '_smpz3d'; clear prop prop.analysis='freq'; clear weakconstr weakconstr.value = 'on'; weakconstr.dim = {'lm1','lm2','lm3','lm4'}; prop.weakconstr = weakconstr; appl.prop = prop; clear pnt pnt.wcshape = [1;2;3;4]; pnt.ind = [1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1, ... 1,1]; appl.pnt = pnt; clear edg edg.wcshape = [1;2;3;4]; edg.ind = [1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1, ... 1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1]; appl.edg = edg; clear bnd bnd.electrictype = {'nD0','nD0','V','V','cont'}; bnd.wcshape = [1;2;3;4]; bnd.V0 = {0,0,-1000,1000,0}; bnd.constrcond = {'sym','free','sym','sym','free'}; bnd.wcgporder = 2; bnd.ind = [1,2,2,2,1,3,2,2,4,2,5,5,2,5,5,2,1,4,2,2,3,2,2,1,1,5,1]; appl.bnd = bnd; clear equ equ.epsilonrT = {{'mat1_epsilonrT_1_1_','mat1_epsilonrT_1_2_','mat1_epsilonrT_1_3_'; ... 'mat1_epsilonrT_2_1_','mat1_epsilonrT_2_2_','mat1_epsilonrT_2_3_';'mat1_epsilonrT_3_1_', ... 'mat1_epsilonrT_3_2_','mat1_epsilonrT_3_3_'},{3130,0,0;0,3130, ... 0;0,0,3400},{3130,0,0;0,3130,0;0,0,3400}}; equ.epsilonrS = {{'mat1_epsilonrS_1_1_','mat1_epsilonrS_1_2_','mat1_epsilonrS_1_3_'; ... 'mat1_epsilonrS_2_1_','mat1_epsilonrS_2_2_','mat1_epsilonrS_2_3_';'mat1_epsilonrS_3_1_', ... 'mat1_epsilonrS_3_2_','mat1_epsilonrS_3_3_'},{1704.40,0,0;0,1704.40, ... 0;0,0,1433.61},{1704.40,0,0;0,1704.40,0;0,0,1433.61}}; equ.gporder = 3; equ.epsilonr = {1,'mat2_epsilonr','mat2_epsilonr'}; equ.rho = {'mat1_rho','mat2_rho','mat2_rho'}; equ.cE = {{'mat1_cE_1_1_','mat1_cE_1_2_','mat1_cE_1_3_','mat1_cE_1_4_', ... 'mat1_cE_1_5_','mat1_cE_1_6_';'mat1_cE_2_1_','mat1_cE_2_2_','mat1_cE_2_3_', ... 'mat1_cE_2_4_','mat1_cE_2_5_','mat1_cE_2_6_';'mat1_cE_3_1_','mat1_cE_3_2_', ... 'mat1_cE_3_3_','mat1_cE_3_4_','mat1_cE_3_5_','mat1_cE_3_6_';'mat1_cE_4_1_', ... 'mat1_cE_4_2_','mat1_cE_4_3_','mat1_cE_4_4_','mat1_cE_4_5_','mat1_cE_4_6_'; ... 'mat1_cE_5_1_','mat1_cE_5_2_','mat1_cE_5_3_','mat1_cE_5_4_','mat1_cE_5_5_', ... 'mat1_cE_5_6_';'mat1_cE_6_1_','mat1_cE_6_2_','mat1_cE_6_3_','mat1_cE_6_4_', ... 'mat1_cE_6_5_','mat1_cE_6_6_'},{1.27205e11,8.02122e10,8.46702e10, ... 0,0,0;8.02122e10,1.27205e11,8.46702e10,0,0,0;8.46702e10, ... 8.46702e10,1.17436e11,0,0,0;0,0,0,2.29886e10,0,0; ... 0,0,0,0,2.29886e10,0;0,0,0,0,0,2.34742e10},{1.27205e11, ... 8.02122e10,8.46702e10,0,0,0;8.02122e10,1.27205e11,8.46702e10, ... 0,0,0;8.46702e10,8.46702e10,1.17436e11,0,0,0;0,0, ... 0,2.29886e10,0,0;0,0,0,0,2.29886e10,0;0,0,0,0, ... 0,2.34742e10}}; equ.wcshape = [1;2;3;4]; equ.shape = [5;6;7;8]; equ.wcgporder = 3; equ.sE = {{'mat1_sE_1_1_','mat1_sE_1_2_','mat1_sE_1_3_','mat1_sE_1_4_', ... 'mat1_sE_1_5_','mat1_sE_1_6_';'mat1_sE_2_1_','mat1_sE_2_2_','mat1_sE_2_3_', ... 'mat1_sE_2_4_','mat1_sE_2_5_','mat1_sE_2_6_';'mat1_sE_3_1_','mat1_sE_3_2_', ... 'mat1_sE_3_3_','mat1_sE_3_4_','mat1_sE_3_5_','mat1_sE_3_6_';'mat1_sE_4_1_', ... 'mat1_sE_4_2_','mat1_sE_4_3_','mat1_sE_4_4_','mat1_sE_4_5_','mat1_sE_4_6_'; ... 'mat1_sE_5_1_','mat1_sE_5_2_','mat1_sE_5_3_','mat1_sE_5_4_','mat1_sE_5_5_', ... 'mat1_sE_5_6_';'mat1_sE_6_1_','mat1_sE_6_2_','mat1_sE_6_3_','mat1_sE_6_4_', ... 'mat1_sE_6_5_','mat1_sE_6_6_'},{16.5e-12,-4.78e-12,-8.45e-12,0, ... 0,0;-4.78e-12,16.5e-12,-8.45e-12,0,0,0;-8.45e-12,-8.45e-12, ... 20.7e-12,0,0,0;0,0,0,43.5e-12,0,0;0,0,0,0,43.5e-12, ... 0;0,0,0,0,0,42.6e-12},{16.5e-12,-4.78e-12,-8.45e-12, ... 0,0,0;-4.78e-12,16.5e-12,-8.45e-12,0,0,0;-8.45e-12, ... -8.45e-12,20.7e-12,0,0,0;0,0,0,43.5e-12,0,0;0,0, ... 0,0,43.5e-12,0;0,0,0,0,0,42.6e-12}}; equ.sigma = {5.99e7,'mat2_sigma','mat2_sigma'}; equ.d = {{'mat1_d_1_1_','mat1_d_1_2_','mat1_d_1_3_','mat1_d_1_4_','mat1_d_1_5_', ... 'mat1_d_1_6_';'mat1_d_2_1_','mat1_d_2_2_','mat1_d_2_3_','mat1_d_2_4_','mat1_d_2_5_', ... 'mat1_d_2_6_';'mat1_d_3_1_','mat1_d_3_2_','mat1_d_3_3_','mat1_d_3_4_','mat1_d_3_5_', ... 'mat1_d_3_6_'},{0,0,0,0,741e-12,0;0,0,0,741e-12,0, ... 0;-274e-12,-274e-12,593e-12,0,0,0},{0,0,0,0,741e-12, ... 0;0,0,0,741e-12,0,0;-274e-12,-274e-12,593e-12,0,0, ... 0}}; equ.e = {{'mat1_e_1_1_','mat1_e_1_2_','mat1_e_1_3_','mat1_e_1_4_','mat1_e_1_5_', ... 'mat1_e_1_6_';'mat1_e_2_1_','mat1_e_2_2_','mat1_e_2_3_','mat1_e_2_4_','mat1_e_2_5_', ... 'mat1_e_2_6_';'mat1_e_3_1_','mat1_e_3_2_','mat1_e_3_3_','mat1_e_3_4_','mat1_e_3_5_', ... 'mat1_e_3_6_'},{0,0,0,0,17.0345,0;0,0,0,17.0345,0, ... 0;-6.62281,-6.62281,23.2403,0,0,0},{0,0,0,0,17.0345, ... 0;0,0,0,17.0345,0,0;-6.62281,-6.62281,23.2403,0,0, ... 0}}; equ.usage = {1,1,0}; equ.ind = [1,2,3,2]; appl.equ = equ; appl.var = {'freq','12e6'}; fem.appl{1} = appl; % Application mode 2 clear appl appl.mode.class = 'GeneralLaminarFlow'; appl.dim = {'u4','v4','w4','p4','phi2','psi2','nxw2','nyw2','nzw2'}; appl.name = 'mmglf2'; appl.module = 'MEMS'; appl.shape = {'shlag(2,''lm12'')','shlag(2,''lm13'')','shlag(2,''lm14'')','shlag(1,''lm15'')','shlag(2,''u4'')','shlag(2,''v4'')','shlag(2,''w4'')','shlag(1,''p4'')'}; appl.gporder = {10,4,2}; appl.cporder = {2,1}; appl.assignsuffix = '_mmglf2'; clear prop prop.analysis='static'; clear weakconstr weakconstr.value = 'on'; weakconstr.dim = {'lm12','lm13','lm14','lm15','lm16','lm17'}; prop.weakconstr = weakconstr; appl.prop = prop; clear bnd bnd.opentype = {'ntotstress','ntotstress','ntotstress','novisc','ntotstress'}; bnd.type = {'int','inlet','walltype','walltype','outlet'}; bnd.ww = {0,0,0,'w_t_smpz3d',0}; bnd.wcshape = [1;2;3;4]; bnd.vw = {0,0,0,'v_t_smpz3d',0}; bnd.U0in = {1,1e-10,1,1,1}; bnd.wcgporder = 1; bnd.intype = {'uv','p','uv','uv','uv'}; bnd.walltype = {'noslip','noslip','noslip','lwall','noslip'}; bnd.uw = {0,0,0,'u_t_smpz3d',0}; bnd.ind = [1,1,1,1,1,1,1,1,1,1,2,3,4,3,3,1,1,1,1,1,1,1,1,1,1,5,1]; appl.bnd = bnd; clear equ equ.eta = {1e-3,1e-8}; equ.gporder = {{2;2;2;3}}; equ.rho = {1e3,1}; equ.cporder = {{1;1;1;2}}; equ.shape = [5;6;7;8]; equ.usage = {0,1}; equ.ind = [1,1,2,1]; appl.equ = equ; fem.appl{2} = appl; fem.frame = {'ref'}; fem.border = 1; clear units; units.basesystem = 'SI'; fem.units = units; % Boundary settings clear bnd bnd.ind = [1,1,1,1,1,1,2,1,1,1,1,1,2,1,1,1,1,1,2,1,1,1,1,1,1,1,1]; bnd.dim = {'u','v','w','V','lm1','lm2','lm3','lm4','u4','v4','w4','p4', ... 'lmx_mmglf2','lmy_mmglf2','lmz_mmglf2','lm12','lm13','lm14','lm15'}; % Boundary expressions bnd.expr = {'Z_acceleration',{'','w_tt_smpz3d'}, ... 'X_acceleration',{'','u_tt_smpz3d'}, ... 'Y_acceleration',{'','v_tt_smpz3d'}}; fem.bnd = bnd; % Coupling variable elements clear elemcpl % Extrusion coupling variables clear elem elem.elem = 'elcplextr'; elem.g = {'1'}; src = cell(1,1); clear bnd bnd.expr = {{{},'u','0*u'},{{},'v','0*v'},{{},'w','w'},{{},'V','0*V'}}; bnd.map = {{'1','1','1'},{'1','1','1'},{'1','1','1'},{'1','1','1'}}; bnd.ind = {{'1','3','4','5','7','8','9','10','11','12','13','14','15', ... '16','17','18','19','20','21','22','23','24','25','26','27'},{'2'},{'6'}}; src{1} = {{},{},bnd,{}}; elem.src = src; geomdim = cell(1,1); clear bnd bnd.map = {{{},'2'},{{},'2'},{{},'2'},{{},'2'}}; bnd.ind = {{'1','2','3','4','5','6','7','8','10','11','12','13','14', ... '15','16','17','18','19','20','21','22','24','25','26','27'},{'9','23'}}; geomdim{1} = {{},{},bnd,{}}; elem.geomdim = geomdim; elem.var = {'pconstr1','pconstr2','pconstr3','pconstr4'}; map = cell(1,2); clear submap submap.type = 'unit'; map{1} = submap; clear submap submap.type = 'linear'; submap.sg = '1'; submap.sv = {'5','6','21','22','15','16','31','32'}; submap.dg = '1'; submap.dv = {'3','4','19','20','1','2','17','18'}; map{2} = submap; elem.map = map; elemcpl{1} = elem; % Point constraint variables (used for periodic conditions) clear elem elem.elem = 'elpconstr'; elem.g = {'1'}; clear bnd bnd.constr = {{'pconstr1-(0)','pconstr2-(0)','pconstr3-(1e-019*w)','pconstr4-(0)'}, ... {'pconstr1-(u)','pconstr2-(v)','pconstr3-(w)','pconstr4-(V)'}}; bnd.cpoints = {{'2','2','2','2'},{'2','2','2','2'}}; bnd.ind = {{'9'},{'23'}}; elem.geomdim = {{{},{},bnd,{}}}; elemcpl{2} = elem; fem.elemcpl = elemcpl; xfem.fem{1} = fem; fem=xfem.fem{2}; fem.sdim = {'x','y'}; fem.border = 1; clear units; units.basesystem = 'SI'; fem.units = units; xfem.fem{2} = fem; % Descriptions clear descr descr.const= {'rho','Density'}; xfem.descr = descr; % Library materials clear lib lib.mat{1}.name='Lithium Niobate'; lib.mat{1}.varname='mat1'; lib.mat{1}.variables.epsilonrT={{'84','0','0'},{'84','0'},{'30'}}; lib.mat{1}.variables.d={{'0[C/N]','0[C/N]','0[C/N]','0[C/N]','6.8e-011[C/N]','-4.2e-011[C/N]'},{'-2.1e-011[C/N]','2.1e-011[C/N]','0[C/N]','6.8e-011[C/N]','0[C/N]','0[C/N]'},{'-1e-012[C/N]','-1e-012[C/N]','6e-012[C/N]','0[C/N]','0[C/N]','0[C/N]'}}; lib.mat{1}.variables.e={{'0[C/m^2]','0[C/m^2]','0[C/m^2]','0[C/m^2]','3.69594[C/m^2]','-2.53384[C/m^2]'},{'-2.53764[C/m^2]','2.53764[C/m^2]','0[C/m^2]','3.69548[C/m^2]','0[C/m^2]','0[C/m^2]'},{'0.193644[C/m^2]','0.193644[C/m^2]','1.30863[C/m^2]','0[C/m^2]','0[C/m^2]','0[C/m^2]'}}; lib.mat{1}.variables.epsilonrS={{'43.6','0','0'},{'43.6','0'},{'29.16'}}; lib.mat{1}.variables.sE={{'5.78e-012[1/Pa]','-1.01e-012[1/Pa]','-1.47e-012[1/Pa]','-1.02e-012[1/Pa]','0[1/Pa]','0[1/Pa]'},{'5.78e-012[1/Pa]','-1.47e-012[1/Pa]','1.02e-012[1/Pa]','0[1/Pa]','0[1/Pa]'},{'5.02e-012[1/Pa]','0[1/Pa]','0[1/Pa]','0[1/Pa]'},{'1.7e-011[1/Pa]','0[1/Pa]','0[1/Pa]'},{'1.7e-011[1/Pa]','-2.04e-012[1/Pa]'},{'1.36e-011[1/Pa]'}}; lib.mat{1}.variables.rho='4700[kg/m^3]'; lib.mat{1}.variables.cE={{'2.02897e+011[Pa]','5.29177e+010[Pa]','7.49098e+010[Pa]','8.99874e+009[Pa]','0[Pa]','0[Pa]'},{'2.02897e+011[Pa]','7.49098e+010[Pa]','-8.99874e+009[Pa]','0[Pa]','0[Pa]'},{'2.43075e+011[Pa]','0[Pa]','0[Pa]','0[Pa]'},{'5.99034e+010[Pa]','0[Pa]','0[Pa]'},{'5.99018e+010[Pa]','8.98526e+009[Pa]'},{'7.48772e+010[Pa]'}}; lib.mat{2}.name='Copper'; lib.mat{2}.varname='mat2'; lib.mat{2}.variables.alphares='3.9e-3[1/K]'; lib.mat{2}.variables.mur='1'; lib.mat{2}.variables.sigma='5.998e7[S/m]'; lib.mat{2}.variables.epsilonr='1'; lib.mat{2}.variables.C='385[J/(kg*K)]'; lib.mat{2}.variables.epsilon='0.5'; lib.mat{2}.variables.res0='1.72e-8[ohm*m]'; lib.mat{2}.variables.rho='8700[kg/m^3]'; lib.mat{2}.variables.k='400[W/(m*K)]'; lib.mat{2}.variables.T0='273.15[K]'; lib.matgroups{1}.name='Electromagnetic (AC/DC)'; lib.matgroups{1}.variables={'alphares','T0','res0','Br','normfH','fH','normfB','fB'}; lib.matgroups{1}.descr={'Temperature coefficient','Reference temperature','Resistivity at reference temperature','Remanent flux density','Nonlinear magnetic field, norm','Nonlinear magnetic field','Nonlinear magnetic flux density, norm','Nonlinear magnetic flux density'}; xfem.lib = lib; % ODE Settings clear ode clear units; units.basesystem = 'SI'; ode.units = units; xfem.ode=ode; % Multiphysics xfem=multiphysics(xfem);