% COMSOL Multiphysics Model M-file


%Initialize
flclear fem

% COMSOL version
clear vrsn
vrsn.name = 'COMSOL 3.4';
vrsn.ext = '';
vrsn.major = 0;
vrsn.build = 250;
vrsn.rcs = '$Name:  $';
vrsn.date = '$Date: 2008/01/24 15:34:30 $';
fem.version = vrsn;





% Construct Geometry

%outer support
g8=cylinder3('35e-3','140e-3','pos',{'0','0','0'},'axis',{'0','0','1'},'rot','0');
g9=cylinder3('25e-3','57.5e-3','pos',{'0','0','10e-3'},'axis',{'0','0','1'},'rot','0');
g10=mirror(g9,[0,0,70e-3],[0,0,1]);
g11=cylinder3('15e-3','100e-3','pos',{'0','0','20e-3'},'axis',{'0','0','1'},'rot','0');
g12=geomcomp({g8,g9,g10,g11},'ns',{'CYL1','CYL2','CO1','CYL3'},'sf','CYL1-CYL2-CO1-CYL3','face','none','edge','all');
g15=cylinder3('5e-3','1','pos',{'0','0','-10e-3'},'axis',{'0','0','1'},'rot','0');

g16=geomcomp({g12,g15},'ns',{'CO2','CYL2'},'sf','CO2-CYL2','face','none','edge','all');

%cavity and endcaps
g13=cylinder3('15e-3','100e-3','pos',{'0','0','20e-3'},'axis',{'0','0','1'},'rot','0');
g17=cylinder3('5e-3','1','pos',{'0','0','0'},'axis',{'0','0','1'},'rot','0');

g18=geomcomp({g13,g17},'ns',{'CYL1','CYL2'},'sf','CYL1-CYL2','face','none','edge','all');
g1=cylinder3('10e-3','5e-3','pos',{'0','0','15e-3'},'axis',{'0','0','1'},'rot','0');
g6=cylinder3('10e-3','5e-3','pos',{'0','0','120e-3'},'axis',{'0','0','1'},'rot','0');

%legs
g19=cylinder3('5e-3','17e-3','pos',{'25e-3','0','-17e-3'},'axis',{'0','0','1'},'rot','0');
g21=block3('10e-3','10e-3','10e-3','base','center','pos',{'25e-3','0','-22e-3'},'axis',{'0','0','1'},'rot','0');
[g24,g25]=geomcopy({g19,g21});
g24=rotate(g24,2.0943951023931953,[0,0,1],[0,0,0]);
g25=rotate(g25,2.0943951023931953,[0,0,1],[0,0,0]);
[g26,g27]=geomcopy({g19,g21});
g26=rotate(g26,-2.0943951023931953,[0,0,1],[0,0,0]);
g27=rotate(g27,-2.0943951023931953,[0,0,1],[0,0,0]);


% Geometry objects
clear s
s.objs={g21,g16,g27,g18,g19,g25,g26,g24,g1,g6};
s.name={'BLK1','CO1','BLK3','CO2','CYL1','BLK2','CYL3','CYL2','CYL4', ...
  'CYL5'};
s.tags={'g21','g16','g27','g18','g19','g25','g26','g24','g1','g6'};

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





% Make materials
% Library materials
clear lib
lib.mat{1}.name='ULE 7972';
lib.mat{1}.varname='mat1';
lib.mat{1}.variables.nu='0.17';
lib.mat{1}.variables.E='67.6[GPa]';
lib.mat{1}.variables.CTE='0';
lib.mat{1}.variables.TD='0.0079[cm^2/s]';
lib.mat{1}.variables.Gxy='29.0[GPa]';
lib.mat{1}.variables.alpha='0';
lib.mat{1}.variables.C='797[J/(kg*K)]';
lib.mat{1}.variables.kappa='34.1[GPa]';
lib.mat{1}.variables.epsilon='0.05';
lib.mat{1}.variables.Syt='49.8[MPa]';
lib.mat{1}.variables.rho='2.21[g/cm^3]';
lib.mat{1}.variables.k='1.31[W/(m*K)]';
lib.mat{2}.name='Zerodur';
lib.mat{2}.varname='mat2';
lib.mat{2}.variables.nu='0.243';
lib.mat{2}.variables.E='90.3[GPa]';
lib.mat{2}.variables.CTE='0';
lib.mat{2}.variables.TD='0.72*10^-6[m^2/s]';
lib.mat{2}.variables.C='0.8[J/(g*K)]';
lib.mat{2}.variables.rho='2.53[g/cm^3]';
lib.mat{2}.variables.k='1.46[W/(m*K)]';

fem.lib = lib;







% Initialize physics

clear appl
appl.mode.class = 'FlSolid3';
appl.gporder = 4;
appl.cporder = 2;
appl.assignsuffix = '_sld';
clear equ
equ.nu = 'mat1_nu';
equ.E = 'mat1_E';
equ.rho = 'mat1_rho';
equ.Fz = {'-21680*1',0};
equ.Hy = {0,1};
equ.Hx = {0,1};
equ.Hz = {0,1};
equ.ind = [1,2,2,1,1,1,1,1,2,1];
appl.equ = equ;
fem.appl{1} = appl;
fem.sdim = {'x','y','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);