% COMSOL Multiphysics Model M-file
% Generated by COMSOL 3.5 (COMSOL 3.5.0.494, $Date: 2008/09/19 16:09:48 $)

flclear fem

% COMSOL version
clear vrsn
vrsn.name = 'COMSOL 3.5';
vrsn.ext = '';
vrsn.major = 0;
vrsn.build = 494;
vrsn.rcs = '$Name:  $';
vrsn.date = '$Date: 2008/09/19 16:09:48 $';
fem.version = vrsn;

% Geometry
g1=rect2(1.4,0.8,'base','corner','pos',[0.4,0]);
g2=rect2('0.08','0.05','base','corner','pos',{'0','-0.02'},'rot','0');
g3=rect2(0.01,0.02,'base','corner','pos',[0.015,-0.0050]);
g4=rect2('0.05','0.003','base','corner','pos',{'0','-0.0030'},'rot','0');
g5=rect2(0.015,0.015,'base','corner','pos',[0.01,0.0050]);
g6=rect2('0.0015','0.015','base','corner','pos',{'0',num2str(lift)},'rot','0');
g7=rect2(0.015,0.015,'base','corner','pos',[0.01,0.0050]);
g8=rect2('0.0015','0.005','base','corner','pos',{'0.0015',num2str(lift)},'rot','0');
g9=rect2(0.02,0.015,'base','corner','pos',[0.0050,0.01]);
g10=rect2('0.0020','0.0050','base','corner','pos',{'0.0030',num2str(lift)},'rot','0');
g11=rect2(0.01,0.0050,'base','corner','pos',[0.01,0.015]);
g12=rect2('0.0020','0.0070','base','corner','pos',{'0.005',num2str(lift)},'rot','0');

% Analyzed geometry
clear s
s.objs={g2,g4,g6,g8,g10,g12};
s.name={'R1','R2','R3','R4','R5','R6'};
s.tags={'g2','g4','g6','g8','g10','g12'};

fem.draw=struct('s',s);
fem.geom=geomcsg(fem);

% Constants
fem.const = {'J','1E9'};

% Constants
fem.const = {'J','1E9'};


% Initialize mesh
fem.mesh=meshinit(fem, ...
                  'hauto',5, ...
                  'hmaxsub',[1,0.001,2,0.001,3,0.001,4,0.001,5,0.001,6,0.001]);

% (Default values are not included)

% Application mode 1
clear appl
appl.mode.class = 'AzimuthalCurrents';
appl.mode.type = 'axi';
appl.module = 'ACDC';
appl.sshape = 2;
appl.border = 'on';
appl.assignsuffix = '_emqa';
clear prop
prop.analysis='transient';
appl.prop = prop;
clear bnd
bnd.type = {'tH0','cont','ax'};
bnd.ind = [3,1,3,2,3,2,3,2,3,2,1,2,2,2,2,2,2,2,2,2,2,2,2,2,1];
appl.bnd = bnd;
clear equ
equ.sigma = {0,3.77e7,0,6e7,6e7};
equ.mur = {1,1,2000,1,1};
equ.Jephi = {0,0,0,0,'J*flc2hs(t-1.5e-3,0.5e-3)'};
equ.ind = [1,2,3,4,3,5];
appl.equ = equ;
fem.appl{1} = appl;
fem.sdim = {'r','z'};
fem.frame = {'ref'};
fem.border = 1;
clear units;
units.basesystem = 'SI';
fem.units = units;


% ODE Settings
clear ode
clear units;
units.basesystem = 'SI';
ode.units = units;
fem.ode=ode;

% Multiphysics
fem=multiphysics(fem);

% Extend mesh
fem.xmesh=meshextend(fem);

% Solve problem
fem.sol=femtime(fem, ...
                'solcomp',{'Aphi'}, ...
                'outcomp',{'Aphi','Aphit'}, ...
                'blocksize','auto', ...
                'tlist',[0:4e-5:4e-3], ...
                'tout','tlist', ...
                'tsteps','strict', ...
                'atol',{'0.010'});

% Save current fem structure for restart purposes
fem0=fem;