% 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
g9=flbinary('g9','draw',flbinaryfile);
draw.s.objs = {g9};
draw.s.name = {'CO2'};
draw.s.tags = {'g9'};
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);
g5=flbinary('g5','draw',flbinaryfile);
g3=flbinary('g3','draw',flbinaryfile);
draw.s.objs = {g1,g25,g5,g3};
draw.s.name = {'CO1','CO2','R1','CO3'};
draw.s.tags = {'g1','g25','g5','g3'};
fem.draw = draw;

% Geometry 2

% Geometry objects
clear s
s.objs={g1,g25,g3,g5};
s.name={'CO1','CO2','CO3','R1'};
s.tags={'g1','g25','g3','g5'};

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.assignsuffix = '_smpz3d';
clear prop
prop.analysis='time';
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];
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];
appl.edg = edg;
clear bnd
bnd.electrictype = {'V','V','nD0','nD0','V0'};
bnd.constrcoord = {'global','local','global','global','local'};
bnd.wcshape = [1;2;3;4];
bnd.V0 = {'40*cos(0.1*t)','40*cos(0.1*t)',0,0,'40*sin(60*t)'};
bnd.Hz = {0,1,1,0,1};
bnd.constrcond = {'free','free','displacement','free','displacement'};
bnd.wcgporder = 2;
bnd.ind = [1,2,3,4,1,5,5,5,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_'}};
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_'}};
equ.gporder = 3;
equ.rho = 'mat1_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_'}};
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_'}};
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_'}};
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_'}};
equ.ind = [1];
appl.equ = equ;
appl.var = {'freq','80e6', ...
  'Qtot','Qes_tot_smpz3d'};
fem.appl{1} = 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,1,1,1,1];
bnd.dim = {'u','v','w','V','lm1','lm2','lm3','lm4'};
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'},{{},'v'},{{},'w'},{{},'V'}};
bnd.map = {{'1','1'},{'1','1'},{'1','1'},{'1','1'}};
bnd.ind = {{'1','3','4','5','6','7','8','9','10'},{'2'}};
src{1} = {{},{},bnd,{}};
elem.src = src;
geomdim = cell(1,1);
clear bnd
bnd.map = {{{},'2'},{{},'2'},{{},'2'},{{},'2'}};
bnd.ind = {{'1','2','3','4','6','7','8','9','10'},{'5'}};
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 = {'3','4','15','16'};
submap.dg = '1';
submap.dv = {'1','2','13','14'};
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-(u)','pconstr2-(v)','pconstr3-(w)','pconstr4-(V)'}};
bnd.cpoints = {{'2','2','2','2'}};
bnd.ind = {{'5'}};
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);