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(0.8,0.4,'base','corner','pos',[-0.6,-0.2]); g2=rect2('1','1','base','corner','pos',{'0','0'},'rot','0'); % Analyzed geometry clear s s.objs={g2}; s.name={'R1'}; s.tags={'g2'}; fem.draw=struct('s',s); fem.geom=geomcsg(fem); % Initialize mesh fem.mesh=meshinit(fem, ... 'hauto',5); % (Default values are not included) % Application mode 1 clear appl appl.mode.class = 'SmeMindlin'; appl.module = 'SME'; appl.gporder = 4; appl.cporder = 1; appl.assignsuffix = '_smdrm'; clear prop prop.analysis='eigen'; appl.prop = prop; clear equ equ.nu = 0.3; equ.thickness = 0.1; equ.rho = 7800; equ.ind = [1]; appl.equ = equ; fem.appl{1} = appl; 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=femeig(fem, ... 'solcomp',{'w','thx','thy','thn_smdrm'}, ... 'outcomp',{'thx','w','thy','thn_smdrm'}, ... 'rowscale','off', ... 'blocksize','auto'); % Save current fem structure for restart purposes fem0=fem; % Plot solution postplot(fem, ... 'tridata',{'w','cont','internal','unit','m'}, ... 'trimap','jet(1024)', ... 'solnum',1, ... 'title','eigfreq_smdrm(1)=0 (1) Surface: z-displacement [m]', ... 'axis',[-0.23994845360824751,1.2399484536082475,-0.05,1.05]); % Plot solution postplot(fem, ... 'tridata',{'w','cont','internal','unit','m'}, ... 'trimap','jet(1024)', ... 'solnum',4, ... 'title','eigfreq_smdrm(4)=312.404097 Surface: z-displacement [m]', ... 'axis',[-0.2020969855832241,1.2020969855832242,-0.05,1.05]); % Plot solution postplot(fem, ... 'tridata',{'w','cont','internal','unit','m'}, ... 'triz',{'w','unit','m'}, ... 'trimap','jet(1024)', ... 'solnum',4, ... 'title','eigfreq_smdrm(4)=312.404097 Surface: z-displacement [m] Height: z-displacement [m]', ... 'grid','on'); % Plot solution postplot(fem, ... 'tridata',{'w','cont','internal','unit','m'}, ... 'triz',{'w','unit','m'}, ... 'trimap','jet(1024)', ... 'solnum',5, ... 'title','eigfreq_smdrm(5)=463.097957 Surface: z-displacement [m] Height: z-displacement [m]', ... 'grid','on'); % Using multiple geometries flclear xfem % Keep global fields fields={'version','const','functions','xmesh','sol','globalexpr','event'}; for ii=1:length(fields) if isfield(fem,fields{ii}) xfem.(fields{ii})=fem.(fields{ii}); fem=rmfield(fem,fields{ii}); end end xfem.fem{1}=fem; % Geometry 1 carr={curve2([0.4,0.4],[1,0],[1,1])}; g3=geomcoerce('curve',carr); gg=geomedit(g3); gg{1}=beziercurve2([0.5,0.5],[1,0],[1,1]); g4=geomedit(g3,gg); fem=xfem.fem{1}; % 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); % Geometry 2 g5=curve3([0.5,0.5],[0,1],[0,0]); xfem.fem{1}=fem; flclear fem % Analyzed geometry clear c c.objs={g5}; c.name={'B1'}; c.tags={'g5'}; fem.draw=struct('c',c); fem.geom=geomcsg(fem); % Initialize mesh for geometry 2 fem.mesh=meshinit(fem, ... 'hauto',5); xfem.fem{2}=fem; fem=xfem.fem{1}; % Initialize mesh for geometry 1 fem.mesh=meshinit(fem, ... 'hauto',5); % Initialize mesh for geometry 1 fem.mesh=meshinit(fem, ... 'hauto',5, ... 'hnumedg',{4,10}); xfem.fem{1}=fem; % (Default values are not included) fem=xfem.fem{1}; % Application mode 1 clear appl appl.mode.class = 'SmeMindlin'; appl.module = 'SME'; appl.gporder = 4; appl.cporder = 1; appl.assignsuffix = '_smdrm'; clear prop prop.analysis='eigen'; appl.prop = prop; clear equ equ.nu = 0.3; equ.thickness = 0.1; equ.rho = 7800; equ.ind = [1,1]; appl.equ = equ; fem.appl{1} = appl; fem.frame = {'ref'}; fem.border = 1; clear units; units.basesystem = 'SI'; fem.units = units; % Coupling variable elements clear elemcpl % Extrusion coupling variables clear elem elem.elem = 'elcplextr'; elem.g = {'1','2'}; src = cell(1,2); clear bnd bnd.expr = {{{},'w'},{{},'thx'},{{},'thy'}}; bnd.map = {{'1','1'},{'1','1'},{'1','1'}}; bnd.ind = {{'1','2','3','5','6','7'},{'4'}}; src{1} = {{},bnd,{}}; src{2} = {}; elem.src = src; geomdim = cell(1,2); geomdim{1} = {}; clear edg edg.map = {{'2'},{'2'},{'2'}}; edg.ind = {{'1'}}; geomdim{2} = {{},edg,{},{}}; elem.geomdim = geomdim; elem.var = {'w1','thx1','thy1'}; map = cell(1,2); clear submap submap.type = 'unit'; map{1} = submap; clear submap submap.type = 'linear'; submap.sg = '2'; submap.sv = {'2','1'}; submap.dg = '1'; submap.dv = {'4','3'}; map{2} = submap; elem.map = map; elemcpl{1} = elem; fem.elemcpl = elemcpl; xfem.fem{1} = fem; fem=xfem.fem{2}; % Application mode 1 clear appl appl.mode.class = 'Sme3DEulerBeam'; appl.module = 'SME'; appl.gporder = 6; appl.cporder = 1; appl.assignsuffix = '_smeul3d'; clear prop prop.analysis='eigen'; appl.prop = prop; clear edg edg.nu = 0.3; edg.Izz = '0.1^4/12'; edg.rho = 7800; edg.constrcoord = 'local'; edg.Hz = 1; edg.Rthy = 'thy1'; edg.Rthx = 'thx1'; edg.J = '2*0.1^4/12'; edg.Rz = 'w1'; edg.Iyy = '0.1^4/12'; edg.constrcond = 'displacement'; edg.Hthx = 1; edg.Hthy = 1; edg.ind = [1]; appl.edg = edg; fem.appl{1} = appl; fem.frame = {'ref'}; fem.border = 1; clear units; units.basesystem = 'SI'; fem.units = units; xfem.fem{2} = fem; % ODE Settings clear ode clear units; units.basesystem = 'SI'; ode.units = units; xfem.ode=ode; % Multiphysics xfem=multiphysics(xfem); % Extend mesh xfem.xmesh=meshextend(xfem); % Evaluate initial value init = asseminit(xfem,'blocksize','auto'); % Update model xfem.sol = init; % Save current fem structure for restart purposes fem0=xfem; % Plot solution postplot(xfem, ... 'tridata',{'w','cont','internal','unit','m'}, ... 'triz',{'w','unit','m'}, ... 'trimap','jet(1024)', ... 'title','Time=0 Surface: z-displacement [m] Height: z-displacement [m]', ... 'grid','on'); % (Default values are not included) fem=xfem.fem{1}; % Application mode 1 clear appl appl.mode.class = 'SmeMindlin'; appl.module = 'SME'; appl.gporder = 4; appl.cporder = 1; appl.assignsuffix = '_smdrm'; clear prop prop.analysis='eigen'; appl.prop = prop; clear equ equ.nu = 0.3; equ.thickness = 0.1; equ.rho = 7800; equ.ind = [1,1]; appl.equ = equ; fem.appl{1} = appl; fem.frame = {'ref'}; fem.border = 1; clear units; units.basesystem = 'SI'; fem.units = units; % Coupling variable elements clear elemcpl % Extrusion coupling variables clear elem elem.elem = 'elcplextr'; elem.g = {'1','2'}; src = cell(1,2); clear bnd bnd.expr = {{{},'w'},{{},'thx'},{{},'thy'}}; bnd.map = {{'1','1'},{'1','1'},{'1','1'}}; bnd.ind = {{'1','2','3','5','6','7'},{'4'}}; src{1} = {{},bnd,{}}; src{2} = {}; elem.src = src; geomdim = cell(1,2); geomdim{1} = {}; clear edg edg.map = {{'2'},{'2'},{'2'}}; edg.ind = {{'1'}}; geomdim{2} = {{},edg,{},{}}; elem.geomdim = geomdim; elem.var = {'w1','thx1','thy1'}; map = cell(1,2); clear submap submap.type = 'unit'; map{1} = submap; clear submap submap.type = 'linear'; submap.sg = '2'; submap.sv = {'2','1'}; submap.dg = '1'; submap.dv = {'4','3'}; map{2} = submap; elem.map = map; elemcpl{1} = elem; fem.elemcpl = elemcpl; xfem.fem{1} = fem; fem=xfem.fem{2}; % Application mode 1 clear appl appl.mode.class = 'Sme3DEulerBeam'; appl.module = 'SME'; appl.gporder = 6; appl.cporder = 1; appl.assignsuffix = '_smeul3d'; clear prop prop.analysis='eigen'; appl.prop = prop; clear edg edg.nu = 0.3; edg.Izz = '0.1^4/12'; edg.rho = 7800; edg.constrcoord = 'local'; edg.Hz = 1; edg.Rthy = 'thy1'; edg.Rthx = 'thx1'; edg.J = '2*0.1^4/12'; edg.Rz = 'w1'; edg.Iyy = '0.1^4/12'; edg.constrcond = 'displacement'; edg.Hthx = 1; edg.Hthy = 1; edg.ind = [1]; appl.edg = edg; fem.appl{1} = appl; fem.frame = {'ref'}; fem.border = 1; clear units; units.basesystem = 'SI'; fem.units = units; xfem.fem{2} = fem; % ODE Settings clear ode clear units; units.basesystem = 'SI'; ode.units = units; xfem.ode=ode; % Multiphysics xfem=multiphysics(xfem); % Extend mesh xfem.xmesh=meshextend(xfem); % Solve problem xfem.sol=femeig(xfem, ... 'init',fem0.sol, ... 'solcomp',{'w','thx','thy','v','u','thz','thn_smdrm'}, ... 'outcomp',{'w','thx','v','thy','u','thz','thn_smdrm'}, ... 'rowscale','off', ... 'blocksize','auto', ... 'neigs',10); % Save current fem structure for restart purposes fem0=xfem; % Plot solution postplot(xfem, ... 'tridata',{'w','cont','internal','unit','m'}, ... 'triz',{'w','unit','m'}, ... 'trimap','jet(1024)', ... 'solnum',1, ... 'title','eigfreq_smdrm(1)=0 (1) Surface: z-displacement [m] Height: z-displacement [m]', ... 'grid','on'); % Plot solution postplot(xfem, ... 'tridata',{'w','cont','internal','unit','m'}, ... 'triz',{'w','unit','m'}, ... 'trimap','jet(1024)', ... 'solnum',1, ... 'title','eigfreq_smdrm(1)=0 (1) Surface: z-displacement [m] Height: z-displacement [m]', ... 'grid','on'); % Plot solution postplot(xfem, ... 'tridata',{'w','cont','internal','unit','m'}, ... 'triz',{'w','unit','m'}, ... 'trimap','jet(1024)', ... 'solnum',6, ... 'title','eigfreq_smdrm(6)=361.354068 Surface: z-displacement [m] Height: z-displacement [m]', ... 'grid','on'); % COMSOL Multiphysics Model M-file % Generated by COMSOL 3.5 (COMSOL 3.5.0.494, $Date: 2008/09/19 16:09:48 $) % (Default values are not included) fem=xfem.fem{1}; % Application mode 1 clear appl appl.mode.class = 'SmeMindlin'; appl.module = 'SME'; appl.gporder = 4; appl.cporder = 1; appl.assignsuffix = '_smdrm'; clear prop prop.analysis='eigen'; appl.prop = prop; clear equ equ.nu = 0.3; equ.thickness = 0.1; equ.rho = 7800; equ.ind = [1,1]; appl.equ = equ; fem.appl{1} = appl; fem.frame = {'ref'}; fem.border = 1; clear units; units.basesystem = 'SI'; fem.units = units; % Coupling variable elements clear elemcpl % Extrusion coupling variables clear elem elem.elem = 'elcplextr'; elem.g = {'1','2'}; src = cell(1,2); clear bnd bnd.expr = {{{},'w'},{{},'thx'},{{},'thy'}}; bnd.map = {{'1','1'},{'1','1'},{'1','1'}}; bnd.ind = {{'1','2','3','5','6','7'},{'4'}}; src{1} = {{},bnd,{}}; src{2} = {}; elem.src = src; geomdim = cell(1,2); geomdim{1} = {}; clear edg edg.map = {{'2'},{'2'},{'2'}}; edg.ind = {{'1'}}; geomdim{2} = {{},edg,{},{}}; elem.geomdim = geomdim; elem.var = {'w1','thx1','thy1'}; map = cell(1,2); clear submap submap.type = 'unit'; map{1} = submap; clear submap submap.type = 'linear'; submap.sg = '2'; submap.sv = {'2','1'}; submap.dg = '1'; submap.dv = {'4','3'}; map{2} = submap; elem.map = map; elemcpl{1} = elem; fem.elemcpl = elemcpl; xfem.fem{1} = fem; fem=xfem.fem{2}; % Application mode 1 clear appl appl.mode.class = 'Sme3DEulerBeam'; appl.module = 'SME'; appl.gporder = 6; appl.cporder = 1; appl.assignsuffix = '_smeul3d'; clear prop prop.analysis='eigen'; appl.prop = prop; clear edg edg.nu = 0.3; edg.Izz = '0.1^4/12'; edg.rho = 7800; edg.constrcoord = 'local'; edg.Hz = 1; edg.Rthy = 'thy1'; edg.Rthx = 'thx1'; edg.J = '2*0.1^4/12'; edg.Rz = 'w1'; edg.Iyy = '0.1^4/12'; edg.constrcond = 'displacement'; edg.Hthx = 1; edg.Hthy = 1; edg.ind = [1]; appl.edg = edg; fem.appl{1} = appl; fem.frame = {'ref'}; fem.border = 1; clear units; units.basesystem = 'SI'; fem.units = units; xfem.fem{2} = fem; % ODE Settings clear ode clear units; units.basesystem = 'SI'; ode.units = units; xfem.ode=ode; % Multiphysics xfem=multiphysics(xfem); % Extend mesh xfem.xmesh=meshextend(xfem); % Solve problem xfem.sol=femeig(xfem, ... 'solcomp',{'thx','w','v','thy','thz','u','thn_smdrm'}, ... 'outcomp',{'thx','w','thy','v','thz','u','thn_smdrm'}, ... 'rowscale','off', ... 'blocksize','auto', ... 'neigs',10); % Save current fem structure for restart purposes fem0=xfem; % Plot solution postplot(xfem, ... 'tridata',{'w','cont','internal','unit','m'}, ... 'triz',{'w','unit','m'}, ... 'trimap','jet(1024)', ... 'solnum',1, ... 'title','eigfreq_smdrm(1)=0 (1) Surface: z-displacement [m] Height: z-displacement [m]', ... 'grid','on'); % Plot solution postplot(xfem, ... 'tridata',{'w','cont','internal','unit','m'}, ... 'triz',{'w','unit','m'}, ... 'trimap','jet(1024)', ... 'solnum',6, ... 'title','eigfreq_smdrm(6)=347.586047 Surface: z-displacement [m] Height: z-displacement [m]', ... 'grid','on'); % Plot solution postplot(xfem, ... 'tridata',{'w','cont','internal','unit','m'}, ... 'triz',{'w','unit','m'}, ... 'trimap','jet(1024)', ... 'solnum',7, ... 'title','eigfreq_smdrm(7)=414.00202 Surface: z-displacement [m] Height: z-displacement [m]', ... 'grid','on'); % Plot solution postplot(xfem, ... 'tridata',{'w','cont','internal','unit','m'}, ... 'triz',{'w','unit','m'}, ... 'trimap','jet(1024)', ... 'solnum',8, ... 'title','eigfreq_smdrm(8)=452.292124 Surface: z-displacement [m] Height: z-displacement [m]', ... 'grid','on'); % Geometry 2 [g1]=geomcopy({g5}); g3=curve3([0.5,0.5],[0,1],[0,0]); fem=xfem.fem{2}; % Analyzed geometry clear c c.objs={g3}; c.name={'B1'}; c.tags={'g3'}; fem.draw=struct('c',c); fem.geom=geomcsg(fem); % Initialize mesh for geometry 2 fem.mesh=meshinit(fem, ... 'hauto',5); xfem.fem{2}=fem; % (Default values are not included) fem=xfem.fem{1}; % Application mode 1 clear appl appl.mode.class = 'SmeMindlin'; appl.module = 'SME'; appl.gporder = 4; appl.cporder = 1; appl.assignsuffix = '_smdrm'; clear prop prop.analysis='eigen'; appl.prop = prop; clear equ equ.nu = 0.3; equ.thickness = 0.1; equ.rho = 7800; equ.ind = [1,1]; appl.equ = equ; fem.appl{1} = appl; fem.frame = {'ref'}; fem.border = 1; clear units; units.basesystem = 'SI'; fem.units = units; % Coupling variable elements clear elemcpl % Extrusion coupling variables clear elem elem.elem = 'elcplextr'; elem.g = {'1','2'}; src = cell(1,2); clear bnd bnd.expr = {{{},'w'},{{},'thx'},{{},'thy'}}; bnd.map = {{'1','1'},{'1','1'},{'1','1'}}; bnd.ind = {{'1','2','3','5','6','7'},{'4'}}; src{1} = {{},bnd,{}}; src{2} = {}; elem.src = src; geomdim = cell(1,2); geomdim{1} = {}; clear edg edg.map = {{'2'},{'2'},{'2'}}; edg.ind = {{'1'}}; geomdim{2} = {{},edg,{},{}}; elem.geomdim = geomdim; elem.var = {'w1','thx1','thy1'}; map = cell(1,2); clear submap submap.type = 'unit'; map{1} = submap; clear submap submap.type = 'linear'; submap.sg = '2'; submap.sv = {'2','1'}; submap.dg = '1'; submap.dv = {'4','3'}; map{2} = submap; elem.map = map; elemcpl{1} = elem; fem.elemcpl = elemcpl; xfem.fem{1} = fem; fem=xfem.fem{2}; % Application mode 1 clear appl appl.mode.class = 'Sme3DEulerBeam'; appl.module = 'SME'; appl.gporder = 6; appl.cporder = 1; appl.assignsuffix = '_smeul3d'; clear prop prop.analysis='eigen'; appl.prop = prop; clear edg edg.nu = 0.3; edg.Izz = '0.1^4/12'; edg.rho = 7800; edg.Hz = 1; edg.Rthy = 'thy1'; edg.Rthx = 'thx1'; edg.J = '2*0.1^4/12'; edg.Rz = 'w1'; edg.Iyy = '0.1^4/12'; edg.constrcond = 'displacement'; edg.Hthx = 1; edg.Hthy = 1; edg.ind = [1]; appl.edg = edg; fem.appl{1} = appl; fem.frame = {'ref'}; fem.border = 1; clear units; units.basesystem = 'SI'; fem.units = units; xfem.fem{2} = fem; % ODE Settings clear ode clear units; units.basesystem = 'SI'; ode.units = units; xfem.ode=ode; % Multiphysics xfem=multiphysics(xfem); % Extend mesh xfem.xmesh=meshextend(xfem); % Solve problem xfem.sol=femeig(xfem, ... 'solcomp',{'thx','w','v','thy','thz','u','thn_smdrm'}, ... 'outcomp',{'thx','w','thy','v','thz','u','thn_smdrm'}, ... 'rowscale','off', ... 'blocksize','auto', ... 'neigs',10); % Save current fem structure for restart purposes fem0=xfem; % Plot solution postplot(xfem, ... 'tridata',{'w','cont','internal','unit','m'}, ... 'triz',{'w','unit','m'}, ... 'trimap','jet(1024)', ... 'solnum',1, ... 'title','eigfreq_smdrm(1)=0 (1) Surface: z-displacement [m] Height: z-displacement [m]', ... 'grid','on'); % Plot solution postplot(xfem, ... 'tridata',{'w','cont','internal','unit','m'}, ... 'triz',{'w','unit','m'}, ... 'trimap','jet(1024)', ... 'solnum',7, ... 'title','eigfreq_smdrm(7)=268.184646 Surface: z-displacement [m] Height: z-displacement [m]', ... 'grid','on'); % Plot solution postplot(xfem, ... 'tridata',{'w','cont','internal','unit','m'}, ... 'triz',{'w','unit','m'}, ... 'trimap','jet(1024)', ... 'solnum',8, ... 'title','eigfreq_smdrm(8)=448.250006 Surface: z-displacement [m] Height: z-displacement [m]', ... 'grid','on'); % Plot solution postplot(xfem, ... 'tridata',{'w','cont','internal','unit','m'}, ... 'triz',{'w','unit','m'}, ... 'trimap','jet(1024)', ... 'solnum',9, ... 'title','eigfreq_smdrm(9)=513.676378 Surface: z-displacement [m] Height: z-displacement [m]', ... 'grid','on'); % Plot solution postplot(xfem, ... 'tridata',{'w','cont','internal','unit','m'}, ... 'triz',{'w','unit','m'}, ... 'trimap','jet(1024)', ... 'solnum','end', ... 'title','eigfreq_smdrm(10)=563.555664 Surface: z-displacement [m] Height: z-displacement [m]', ... 'grid','on'); % (Default values are not included) fem=xfem.fem{1}; % Application mode 1 clear appl appl.mode.class = 'SmeMindlin'; appl.module = 'SME'; appl.gporder = 4; appl.cporder = 1; appl.assignsuffix = '_smdrm'; clear prop prop.analysis='eigen'; appl.prop = prop; clear equ equ.nu = 0.3; equ.thickness = 0.1; equ.rho = 7800; equ.ind = [1,1]; appl.equ = equ; fem.appl{1} = appl; fem.frame = {'ref'}; fem.border = 1; clear units; units.basesystem = 'SI'; fem.units = units; % Coupling variable elements clear elemcpl % Extrusion coupling variables clear elem elem.elem = 'elcplextr'; elem.g = {'1','2'}; src = cell(1,2); clear bnd bnd.expr = {{{},'w'},{{},'thx'},{{},'thy'}}; bnd.map = {{'1','1'},{'1','1'},{'1','1'}}; bnd.ind = {{'1','2','3','5','6','7'},{'4'}}; src{1} = {{},bnd,{}}; src{2} = {}; elem.src = src; geomdim = cell(1,2); geomdim{1} = {}; clear edg edg.map = {{'2'},{'2'},{'2'}}; edg.ind = {{'1'}}; geomdim{2} = {{},edg,{},{}}; elem.geomdim = geomdim; elem.var = {'w1','thx1','thy1'}; map = cell(1,2); clear submap submap.type = 'unit'; map{1} = submap; clear submap submap.type = 'linear'; submap.sg = '2'; submap.sv = {'2','1'}; submap.dg = '1'; submap.dv = {'4','3'}; map{2} = submap; elem.map = map; elemcpl{1} = elem; fem.elemcpl = elemcpl; xfem.fem{1} = fem; fem=xfem.fem{2}; % Application mode 1 clear appl appl.mode.class = 'Sme3DEulerBeam'; appl.module = 'SME'; appl.gporder = 6; appl.cporder = 1; appl.assignsuffix = '_smeul3d'; clear prop prop.analysis='eigen'; appl.prop = prop; clear edg edg.nu = 0.3; edg.Izz = '0.1^4/12'; edg.rho = 7800; edg.Hz = 1; edg.Rthy = 'thy1'; edg.Rthx = 'thx1'; edg.J = '2*0.1^4/12'; edg.Rz = 'w1'; edg.Iyy = '0.1^4/12'; edg.constrcond = 'displacement'; edg.Hthx = 1; edg.Hthy = 1; edg.ind = [1]; appl.edg = edg; fem.appl{1} = appl; fem.frame = {'ref'}; fem.border = 1; clear units; units.basesystem = 'SI'; fem.units = units; xfem.fem{2} = fem; % ODE Settings clear ode clear units; units.basesystem = 'SI'; ode.units = units; xfem.ode=ode; % Multiphysics xfem=multiphysics(xfem); % Extend mesh xfem.xmesh=meshextend(xfem); % Solve problem xfem.sol=femeig(xfem, ... 'solcomp',{'w','thx','thy','v','u','thz','thn_smdrm'}, ... 'outcomp',{'w','thx','v','thy','u','thz','thn_smdrm'}, ... 'rowscale','off', ... 'blocksize','auto', ... 'neigs',20); % Save current fem structure for restart purposes fem0=xfem; % Plot solution postplot(xfem, ... 'tridata',{'w','cont','internal','unit','m'}, ... 'triz',{'w','unit','m'}, ... 'trimap','jet(1024)', ... 'solnum',1, ... 'title','eigfreq_smdrm(1)=0 (1) Surface: z-displacement [m] Height: z-displacement [m]', ... 'grid','on'); % Plot solution postplot(xfem, ... 'tridata',{'w','cont','internal','unit','m'}, ... 'triz',{'w','unit','m'}, ... 'trimap','jet(1024)', ... 'solnum',8, ... 'title','eigfreq_smdrm(8)=321.9655 Surface: z-displacement [m] Height: z-displacement [m]', ... 'grid','on'); % Plot solution postplot(xfem, ... 'tridata',{'w','cont','internal','unit','m'}, ... 'triz',{'w','unit','m'}, ... 'trimap','jet(1024)', ... 'solnum',9, ... 'title','eigfreq_smdrm(9)=446.812651 Surface: z-displacement [m] Height: z-displacement [m]', ... 'grid','on'); % Plot solution postplot(xfem, ... 'tridata',{'w','cont','internal','unit','m'}, ... 'triz',{'w','unit','m'}, ... 'trimap','jet(1024)', ... 'solnum',10, ... 'title','eigfreq_smdrm(10)=518.951796 Surface: z-displacement [m] Height: z-displacement [m]', ... 'grid','on'); % Plot solution postplot(xfem, ... 'tridata',{'w','cont','internal','unit','m'}, ... 'triz',{'w','unit','m'}, ... 'trimap','jet(1024)', ... 'solnum',11, ... 'title','eigfreq_smdrm(11)=577.111255 Surface: z-displacement [m] Height: z-displacement [m]', ... 'grid','on'); % Plot solution postplot(xfem, ... 'tridata',{'w','cont','internal','unit','m'}, ... 'triz',{'w','unit','m'}, ... 'trimap','jet(1024)', ... 'solnum',12, ... 'title','eigfreq_smdrm(12)=736.581964 Surface: z-displacement [m] Height: z-displacement [m]', ... 'grid','on'); % (Default values are not included) fem=xfem.fem{1}; % Application mode 1 clear appl appl.mode.class = 'SmeMindlin'; appl.module = 'SME'; appl.gporder = 4; appl.cporder = 1; appl.assignsuffix = '_smdrm'; clear prop prop.analysis='eigen'; appl.prop = prop; clear equ equ.nu = 0.3; equ.thickness = 0.1; equ.rho = 7800; equ.ind = [1,1]; appl.equ = equ; fem.appl{1} = appl; fem.frame = {'ref'}; fem.border = 1; clear units; units.basesystem = 'SI'; fem.units = units; % Coupling variable elements clear elemcpl % Extrusion coupling variables clear elem elem.elem = 'elcplextr'; elem.g = {'1','2'}; src = cell(1,2); clear bnd bnd.expr = {{{},'w'},{{},'thx'},{{},'thy'}}; bnd.map = {{'1','1'},{'1','1'},{'1','1'}}; bnd.ind = {{'1','2','3','5','6','7'},{'4'}}; src{1} = {{},bnd,{}}; src{2} = {}; elem.src = src; geomdim = cell(1,2); geomdim{1} = {}; clear edg edg.map = {{'2'},{'2'},{'2'}}; edg.ind = {{'1'}}; geomdim{2} = {{},edg,{},{}}; elem.geomdim = geomdim; elem.var = {'w1','thx1','thy1'}; map = cell(1,2); clear submap submap.type = 'unit'; map{1} = submap; clear submap submap.type = 'linear'; submap.sg = '2'; submap.sv = {'2','1'}; submap.dg = '1'; submap.dv = {'4','3'}; map{2} = submap; elem.map = map; elemcpl{1} = elem; fem.elemcpl = elemcpl; xfem.fem{1} = fem; fem=xfem.fem{2}; % Application mode 1 clear appl appl.mode.class = 'Sme3DEulerBeam'; appl.module = 'SME'; appl.gporder = 6; appl.cporder = 1; appl.assignsuffix = '_smeul3d'; clear prop prop.analysis='eigen'; appl.prop = prop; clear edg edg.nu = 0.3; edg.Izz = '0.1^4/12'; edg.rho = 7800; edg.Hz = 1; edg.Rthy = 'thy1'; edg.Rthx = 'thx1'; edg.J = '2*0.1^4/12'; edg.Rz = 'w1'; edg.Iyy = '0.1^4/12'; edg.constrcond = 'displacement'; edg.Hthx = 1; edg.Hthy = 1; edg.ind = [1]; appl.edg = edg; fem.appl{1} = appl; fem.frame = {'ref'}; fem.border = 1; clear units; units.basesystem = 'SI'; fem.units = units; xfem.fem{2} = fem; % ODE Settings clear ode clear units; units.basesystem = 'SI'; ode.units = units; xfem.ode=ode; % Multiphysics xfem=multiphysics(xfem); % Extend mesh xfem.xmesh=meshextend(xfem); % Solve problem xfem.sol=femeig(xfem, ... 'solcomp',{'w','thx','thy','v','u','thz','thn_smdrm'}, ... 'outcomp',{'w','thx','v','thy','u','thz','thn_smdrm'}, ... 'rowscale','off', ... 'blocksize','auto', ... 'neigs',20); % Save current fem structure for restart purposes fem0=xfem; % Plot solution postplot(xfem, ... 'tridata',{'w','cont','internal','unit','m'}, ... 'triz',{'w','unit','m'}, ... 'trimap','jet(1024)', ... 'solnum',1, ... 'title','eigfreq_smdrm(1)=0 (1) Surface: z-displacement [m] Height: z-displacement [m]', ... 'grid','on'); % Plot solution postplot(xfem, ... 'tridata',{'w','cont','internal','unit','m'}, ... 'triz',{'w','unit','m'}, ... 'trimap','jet(1024)', ... 'solnum',8, ... 'title','eigfreq_smdrm(8)=319.844163 Surface: z-displacement [m] Height: z-displacement [m]', ... 'grid','on'); % Plot solution postplot(xfem, ... 'tridata',{'w','cont','internal','unit','m'}, ... 'triz',{'w','unit','m'}, ... 'trimap','jet(1024)', ... 'solnum',9, ... 'title','eigfreq_smdrm(9)=363.043932 (1) Surface: z-displacement [m] Height: z-displacement [m]', ... 'grid','on'); % (Default values are not included) fem=xfem.fem{1}; % Application mode 1 clear appl appl.mode.class = 'SmeMindlin'; appl.module = 'SME'; appl.gporder = 4; appl.cporder = 1; appl.assignsuffix = '_smdrm'; clear prop prop.analysis='eigen'; appl.prop = prop; clear equ equ.nu = 0.3; equ.thickness = 0.1; equ.rho = 7800; equ.ind = [1,1]; appl.equ = equ; fem.appl{1} = appl; fem.frame = {'ref'}; fem.border = 1; clear units; units.basesystem = 'SI'; fem.units = units; % Coupling variable elements clear elemcpl % Extrusion coupling variables clear elem elem.elem = 'elcplextr'; elem.g = {'1','2'}; src = cell(1,2); clear bnd bnd.expr = {{{},'w'},{{},'thx'},{{},'thy'}}; bnd.map = {{'1','1'},{'1','1'},{'1','1'}}; bnd.ind = {{'1','2','3','5','6','7'},{'4'}}; src{1} = {{},bnd,{}}; src{2} = {}; elem.src = src; geomdim = cell(1,2); geomdim{1} = {}; clear edg edg.map = {{'2'},{'2'},{'2'}}; edg.ind = {{'1'}}; geomdim{2} = {{},edg,{},{}}; elem.geomdim = geomdim; elem.var = {'w1','thx1','thy1'}; map = cell(1,2); clear submap submap.type = 'unit'; map{1} = submap; clear submap submap.type = 'linear'; submap.sg = '2'; submap.sv = {'2','1'}; submap.dg = '1'; submap.dv = {'4','3'}; map{2} = submap; elem.map = map; elemcpl{1} = elem; fem.elemcpl = elemcpl; xfem.fem{1} = fem; fem=xfem.fem{2}; % Application mode 1 clear appl appl.mode.class = 'Sme3DEulerBeam'; appl.module = 'SME'; appl.gporder = 6; appl.cporder = 1; appl.assignsuffix = '_smeul3d'; clear prop prop.analysis='eigen'; appl.prop = prop; clear edg edg.nu = 0.3; edg.Izz = '0.1^4/12'; edg.rho = 7800; edg.Hz = 1; edg.Rthy = 'thy1'; edg.Rthx = 'thx1'; edg.J = '2*0.1^4/12'; edg.Rz = 'w1'; edg.Iyy = '0.1^4/12'; edg.constrcond = 'displacement'; edg.Hthx = 1; edg.Hthy = 1; edg.ind = [1]; appl.edg = edg; fem.appl{1} = appl; fem.frame = {'ref'}; fem.border = 1; clear units; units.basesystem = 'SI'; fem.units = units; xfem.fem{2} = fem; % ODE Settings clear ode clear units; units.basesystem = 'SI'; ode.units = units; xfem.ode=ode; % Multiphysics xfem=multiphysics(xfem); % Extend mesh xfem.xmesh=meshextend(xfem); % Solve problem xfem.sol=femeig(xfem, ... 'solcomp',{'w','thx','thy','v','u','thz','thn_smdrm'}, ... 'outcomp',{'w','thx','v','thy','u','thz','thn_smdrm'}, ... 'rowscale','off', ... 'blocksize','auto', ... 'neigs',20); % Save current fem structure for restart purposes fem0=xfem; % Plot solution postplot(xfem, ... 'tridata',{'w','cont','internal','unit','m'}, ... 'triz',{'w','unit','m'}, ... 'trimap','jet(1024)', ... 'solnum',1, ... 'title','eigfreq_smdrm(1)=0 (1) Surface: z-displacement [m] Height: z-displacement [m]', ... 'grid','on'); % Plot solution postplot(xfem, ... 'tridata',{'w','cont','internal','unit','m'}, ... 'triz',{'w','unit','m'}, ... 'trimap','jet(1024)', ... 'solnum',7, ... 'title','eigfreq_smdrm(7)=313.698433 Surface: z-displacement [m] Height: z-displacement [m]', ... 'grid','on'); % Plot solution postplot(xfem, ... 'tridata',{'w','cont','internal','unit','m'}, ... 'triz',{'w','unit','m'}, ... 'trimap','jet(1024)', ... 'solnum',8, ... 'title','eigfreq_smdrm(8)=451.247551 Surface: z-displacement [m] Height: z-displacement [m]', ... 'grid','on'); % Plot solution postplot(xfem, ... 'tridata',{'w','cont','internal','unit','m'}, ... 'triz',{'w','unit','m'}, ... 'trimap','jet(1024)', ... 'solnum',9, ... 'title','eigfreq_smdrm(9)=521.759498 Surface: z-displacement [m] Height: z-displacement [m]', ... 'grid','on'); % Plot solution postplot(xfem, ... 'tridata',{'w','cont','internal','unit','m'}, ... 'triz',{'w','unit','m'}, ... 'trimap','jet(1024)', ... 'solnum',7, ... 'title','eigfreq_smdrm(7)=313.698433 Surface: z-displacement [m] Height: z-displacement [m]', ... 'grid','on'); % Plot solution postplot(xfem, ... 'tridata',{'w','cont','internal','unit','m'}, ... 'triz',{'w','unit','m'}, ... 'trimap','jet(1024)', ... 'solnum',8, ... 'title','eigfreq_smdrm(8)=451.247551 Surface: z-displacement [m] Height: z-displacement [m]', ... 'grid','on'); % Plot solution postplot(xfem, ... 'tridata',{'w','cont','internal','unit','m'}, ... 'triz',{'w','unit','m'}, ... 'trimap','jet(1024)', ... 'solnum',9, ... 'title','eigfreq_smdrm(9)=521.759498 Surface: z-displacement [m] Height: z-displacement [m]', ... 'grid','on'); % Plot solution postplot(xfem, ... 'tridata',{'w','cont','internal','unit','m'}, ... 'triz',{'w','unit','m'}, ... 'trimap','jet(1024)', ... 'solnum',11, ... 'title','eigfreq_smdrm(11)=750.527356 Surface: z-displacement [m] Height: z-displacement [m]', ... 'grid','on'); % COMSOL Multiphysics Model M-file % Generated by COMSOL 3.5 (COMSOL 3.5.0.494, $Date: 2008/09/19 16:09:48 $) % (Default values are not included) fem=xfem.fem{1}; % Application mode 1 clear appl appl.mode.class = 'SmeMindlin'; appl.module = 'SME'; appl.gporder = 4; appl.cporder = 1; appl.border = 'on'; appl.assignsuffix = '_smdrm'; clear prop prop.analysis='eigen'; appl.prop = prop; clear bnd bnd.Rthy = {0,'thy2'}; bnd.Rthx = {0,'thx2'}; bnd.constrcoord = {'local','global'}; bnd.Hz = {0,1}; bnd.Rz = {0,'w2'}; bnd.Hthx = {0,1}; bnd.Hthy = {0,1}; bnd.ind = [1,1,1,2,1,1,1]; appl.bnd = bnd; clear equ equ.nu = 0.3; equ.thickness = 0.1; equ.rho = 7800; equ.ind = [1,1]; appl.equ = equ; fem.appl{1} = appl; fem.frame = {'ref'}; fem.border = 1; clear units; units.basesystem = 'SI'; fem.units = units; % Coupling variable elements clear elemcpl % Extrusion coupling variables clear elem elem.elem = 'elcplextr'; elem.g = {'1','2'}; src = cell(1,2); clear bnd bnd.expr = {{{},'w'},{{},'thx'},{{},'thy'}}; bnd.map = {{'1','1'},{'1','1'},{'1','1'}}; bnd.ind = {{'1','2','3','5','6','7'},{'4'}}; src{1} = {{},bnd,{}}; src{2} = {}; elem.src = src; geomdim = cell(1,2); geomdim{1} = {}; clear edg edg.map = {{'2'},{'2'},{'2'}}; edg.ind = {{'1'}}; geomdim{2} = {{},edg,{},{}}; elem.geomdim = geomdim; elem.var = {'w1','thx1','thy1'}; map = cell(1,2); clear submap submap.type = 'unit'; map{1} = submap; clear submap submap.type = 'linear'; submap.sg = '2'; submap.sv = {'2','1'}; submap.dg = '1'; submap.dv = {'4','3'}; map{2} = submap; elem.map = map; elemcpl{1} = elem; fem.elemcpl = elemcpl; xfem.fem{1} = fem; fem=xfem.fem{2}; % Application mode 1 clear appl appl.mode.class = 'Sme3DEulerBeam'; appl.module = 'SME'; appl.gporder = 6; appl.cporder = 1; appl.assignsuffix = '_smeul3d'; clear prop prop.analysis='eigen'; appl.prop = prop; clear edg edg.nu = 0.3; edg.Izz = '0.1^4/12'; edg.rho = 7800; edg.Hz = 1; edg.Rthy = 'thy1'; edg.Rthx = 'thx1'; edg.J = '2*0.1^4/12'; edg.Rz = 'w1'; edg.Iyy = '0.1^4/12'; edg.constrcond = 'displacement'; edg.Hthx = 1; edg.Hthy = 1; edg.ind = [1]; appl.edg = edg; fem.appl{1} = appl; fem.frame = {'ref'}; fem.border = 1; clear units; units.basesystem = 'SI'; fem.units = units; % Coupling variable elements clear elemcpl % Extrusion coupling variables clear elem elem.elem = 'elcplextr'; elem.g = {'1','2'}; src = cell(1,2); src{1} = {}; clear edg edg.expr = {{'w'},{'thx'},{'thy'}}; edg.map = {{'1'},{'1'},{'1'}}; edg.ind = {{'1'}}; src{2} = {{},edg,{},{}}; elem.src = src; geomdim = cell(1,2); clear bnd bnd.map = {{{},'2'},{{},'2'},{{},'2'}}; bnd.ind = {{'1','2','3','5','6','7'},{'4'}}; geomdim{1} = {{},bnd,{}}; geomdim{2} = {}; elem.geomdim = geomdim; elem.var = {'w2','thx2','thy2'}; map = cell(1,2); clear submap submap.type = 'unit'; map{1} = submap; clear submap submap.type = 'linear'; submap.sg = '1'; submap.sv = {'4','3'}; submap.dg = '2'; submap.dv = {'2','1'}; map{2} = submap; elem.map = map; elemcpl{1} = elem; fem.elemcpl = elemcpl; xfem.fem{2} = fem; % ODE Settings clear ode clear units; units.basesystem = 'SI'; ode.units = units; xfem.ode=ode; % Multiphysics xfem=multiphysics(xfem); % Extend mesh xfem.xmesh=meshextend(xfem); % Solve problem xfem.sol=femeig(xfem, ... 'solcomp',{'w','thx','thy','v','u','thz','thn_smdrm'}, ... 'outcomp',{'w','thx','v','thy','u','thz','thn_smdrm'}, ... 'rowscale','off', ... 'blocksize','auto', ... 'neigs',20); % Save current fem structure for restart purposes fem0=xfem; % Plot solution postplot(xfem, ... 'tridata',{'w','cont','internal','unit','m'}, ... 'triz',{'w','unit','m'}, ... 'trimap','jet(1024)', ... 'solnum',1, ... 'title','eigfreq_smdrm(1)=0 (1) Surface: z-displacement [m] Height: z-displacement [m]', ... 'grid','on'); % Plot solution postplot(xfem, ... 'tridata',{'w','cont','internal','unit','m'}, ... 'triz',{'w','unit','m'}, ... 'trimap','jet(1024)', ... 'solnum',7, ... 'title','eigfreq_smdrm(7)=221.196567 (1) Surface: z-displacement [m] Height: z-displacement [m]', ... 'grid','on'); % Plot solution postplot(xfem, ... 'tridata',{'w','cont','internal','unit','m'}, ... 'triz',{'w','unit','m'}, ... 'trimap','jet(1024)', ... 'solnum',9, ... 'title','eigfreq_smdrm(9)=307.829196 Surface: z-displacement [m] Height: z-displacement [m]', ... 'grid','on'); % Plot solution postplot(xfem, ... 'tridata',{'w','cont','internal','unit','m'}, ... 'triz',{'w','unit','m'}, ... 'trimap','jet(1024)', ... 'solnum',10, ... 'title','eigfreq_smdrm(10)=425.878522 Surface: z-displacement [m] Height: z-displacement [m]', ... 'grid','on'); % Plot solution postplot(xfem, ... 'tridata',{'w','cont','internal','unit','m'}, ... 'triz',{'w','unit','m'}, ... 'trimap','jet(1024)', ... 'solnum',9, ... 'title','eigfreq_smdrm(9)=307.829196 Surface: z-displacement [m] Height: z-displacement [m]', ... 'grid','on'); % Plot solution postplot(xfem, ... 'tridata',{'w','cont','internal','unit','m'}, ... 'triz',{'w','unit','m'}, ... 'trimap','jet(1024)', ... 'solnum',10, ... 'title','eigfreq_smdrm(10)=425.878522 Surface: z-displacement [m] Height: z-displacement [m]', ... 'grid','on'); % Plot solution postplot(xfem, ... 'tridata',{'w','cont','internal','unit','m'}, ... 'triz',{'w','unit','m'}, ... 'trimap','jet(1024)', ... 'solnum',11, ... 'title','eigfreq_smdrm(11)=613.73189 Surface: z-displacement [m] Height: z-displacement [m]', ... 'grid','on'); % (Default values are not included) fem=xfem.fem{1}; % Application mode 1 clear appl appl.mode.class = 'SmeMindlin'; appl.module = 'SME'; appl.gporder = 4; appl.cporder = 1; appl.border = 'on'; appl.assignsuffix = '_smdrm'; clear prop prop.analysis='eigen'; appl.prop = prop; clear bnd bnd.Rthy = {0,'thy2'}; bnd.Rthx = {0,'thx2'}; bnd.constrcoord = {'local','global'}; bnd.Hz = {0,1}; bnd.Rz = {0,'w2'}; bnd.Hthx = {0,1}; bnd.Hthy = {0,1}; bnd.ind = [1,1,1,2,1,1,1]; appl.bnd = bnd; clear equ equ.nu = 0.3; equ.thickness = 0.1; equ.rho = 7800; equ.ind = [1,1]; appl.equ = equ; fem.appl{1} = appl; fem.frame = {'ref'}; fem.border = 1; clear units; units.basesystem = 'SI'; fem.units = units; % Coupling variable elements clear elemcpl % Extrusion coupling variables clear elem elem.elem = 'elcplextr'; elem.g = {'1','2'}; src = cell(1,2); clear bnd bnd.expr = {{{},'w'},{{},'thx'},{{},'thy'}}; bnd.map = {{'1','1'},{'1','1'},{'1','1'}}; bnd.ind = {{'1','2','3','5','6','7'},{'4'}}; src{1} = {{},bnd,{}}; src{2} = {}; elem.src = src; geomdim = cell(1,2); geomdim{1} = {}; clear edg edg.map = {{'2'},{'2'},{'2'}}; edg.ind = {{'1'}}; geomdim{2} = {{},edg,{},{}}; elem.geomdim = geomdim; elem.var = {'w1','thx1','thy1'}; map = cell(1,2); clear submap submap.type = 'unit'; map{1} = submap; clear submap submap.type = 'linear'; submap.sg = '2'; submap.sv = {'2','1'}; submap.dg = '1'; submap.dv = {'4','3'}; map{2} = submap; elem.map = map; elemcpl{1} = elem; fem.elemcpl = elemcpl; xfem.fem{1} = fem; fem=xfem.fem{2}; % Application mode 1 clear appl appl.mode.class = 'Sme3DEulerBeam'; appl.module = 'SME'; appl.gporder = 6; appl.cporder = 1; appl.assignsuffix = '_smeul3d'; clear prop prop.analysis='eigen'; appl.prop = prop; clear edg edg.nu = 0.3; edg.Izz = '0.1^4/12'; edg.rho = 7800; edg.Hz = 1; edg.Rthy = 'thy1'; edg.Rthx = 'thx1'; edg.J = '2*0.1^4/12'; edg.Rz = 'w1'; edg.Iyy = '0.1^4/12'; edg.constrcond = 'displacement'; edg.Hthx = 1; edg.Hthy = 1; edg.ind = [1]; appl.edg = edg; fem.appl{1} = appl; fem.frame = {'ref'}; fem.border = 1; clear units; units.basesystem = 'SI'; fem.units = units; % Coupling variable elements clear elemcpl % Extrusion coupling variables clear elem elem.elem = 'elcplextr'; elem.g = {'1','2'}; src = cell(1,2); src{1} = {}; clear edg edg.expr = {{'w'},{'thx'},{'thy'}}; edg.map = {{'1'},{'1'},{'1'}}; edg.ind = {{'1'}}; src{2} = {{},edg,{},{}}; elem.src = src; geomdim = cell(1,2); clear bnd bnd.map = {{{},'2'},{{},'2'},{{},'2'}}; bnd.ind = {{'1','2','3','5','6','7'},{'4'}}; geomdim{1} = {{},bnd,{}}; geomdim{2} = {}; elem.geomdim = geomdim; elem.var = {'w2','thx2','thy2'}; map = cell(1,2); clear submap submap.type = 'unit'; map{1} = submap; clear submap submap.type = 'linear'; submap.sg = '1'; submap.sv = {'4','3'}; submap.dg = '2'; submap.dv = {'2','1'}; map{2} = submap; elem.map = map; elemcpl{1} = elem; fem.elemcpl = elemcpl; xfem.fem{2} = fem; % ODE Settings clear ode clear units; units.basesystem = 'SI'; ode.units = units; xfem.ode=ode; % Multiphysics xfem=multiphysics(xfem); % Extend mesh xfem.xmesh=meshextend(xfem); % Solve problem xfem.sol=femeig(xfem, ... 'solcomp',{'w','thx','thy','thn_smdrm'}, ... 'outcomp',{'w','thx','v','thy','u','thz','thn_smdrm'}, ... 'rowscale','off', ... 'blocksize','auto', ... 'neigs',20); % Save current fem structure for restart purposes fem0=xfem; % Plot solution postplot(xfem, ... 'tridata',{'w','cont','internal','unit','m'}, ... 'triz',{'w','unit','m'}, ... 'trimap','jet(1024)', ... 'solnum',1, ... 'title','eigfreq_smdrm(1)=0 (1) Surface: z-displacement [m] Height: z-displacement [m]', ... 'grid','on'); % Plot solution postplot(xfem, ... 'tridata',{'w','cont','internal','unit','m'}, ... 'triz',{'w','unit','m'}, ... 'trimap','jet(1024)', ... 'solnum',5, ... 'title','eigfreq_smdrm(5)=260.851107 Surface: z-displacement [m] Height: z-displacement [m]', ... 'grid','on'); % Plot solution postplot(xfem, ... 'tridata',{'w','cont','internal','unit','m'}, ... 'triz',{'w','unit','m'}, ... 'trimap','jet(1024)', ... 'solnum',6, ... 'title','eigfreq_smdrm(6)=318.529411 Surface: z-displacement [m] Height: z-displacement [m]', ... 'grid','on'); % Plot solution postplot(xfem, ... 'tridata',{'w','cont','internal','unit','m'}, ... 'triz',{'w','unit','m'}, ... 'trimap','jet(1024)', ... 'solnum',7, ... 'title','eigfreq_smdrm(7)=447.379825 Surface: z-displacement [m] Height: z-displacement [m]', ... 'grid','on'); % Plot solution postplot(xfem, ... 'tridata',{'w','cont','internal','unit','m'}, ... 'triz',{'w','unit','m'}, ... 'trimap','jet(1024)', ... 'solnum',8, ... 'title','eigfreq_smdrm(8)=677.85519 Surface: z-displacement [m] Height: z-displacement [m]', ... 'grid','on'); % COMSOL Multiphysics Model M-file % Generated by COMSOL 3.5 (COMSOL 3.5.0.494, $Date: 2008/09/19 16:09:48 $) % Plot solution postplot(xfem, ... 'tridata',{'w','cont','internal','unit','m'}, ... 'triz',{'w','unit','m'}, ... 'trimap','jet(1024)', ... 'solnum',8, ... 'title','eigfreq_smdrm(8)=677.85519 Surface: z-displacement [m] Height: z-displacement [m]', ... 'grid','on'); % Plot solution postplot(xfem, ... 'tridata',{'w','cont','internal','unit','m'}, ... 'triz',{'w','unit','m'}, ... 'trimap','jet(1024)', ... 'solnum',5, ... 'title','eigfreq_smdrm(5)=260.851107 Surface: z-displacement [m] Height: z-displacement [m]', ... 'grid','on'); % (Default values are not included) fem=xfem.fem{1}; % Application mode 1 clear appl appl.mode.class = 'SmeMindlin'; appl.module = 'SME'; appl.gporder = 4; appl.cporder = 1; appl.assignsuffix = '_smdrm'; clear prop prop.analysis='eigen'; appl.prop = prop; clear bnd bnd.Rthy = {0,'thy2'}; bnd.Rthx = {0,'thx2'}; bnd.constrcoord = {'local','global'}; bnd.Hz = {0,1}; bnd.Rz = {0,'w2'}; bnd.Hthx = {0,1}; bnd.Hthy = {0,1}; bnd.ind = [1,1,1,2,1,1,1]; appl.bnd = bnd; clear equ equ.nu = 0.3; equ.thickness = 0.1; equ.rho = 7800; equ.ind = [1,1]; appl.equ = equ; fem.appl{1} = appl; fem.frame = {'ref'}; fem.border = 1; clear units; units.basesystem = 'SI'; fem.units = units; % Coupling variable elements clear elemcpl % Extrusion coupling variables clear elem elem.elem = 'elcplextr'; elem.g = {'1','2'}; src = cell(1,2); clear bnd bnd.expr = {{{},'w'},{{},'thx'},{{},'thy'}}; bnd.map = {{'1','1'},{'1','1'},{'1','1'}}; bnd.ind = {{'1','2','3','5','6','7'},{'4'}}; src{1} = {{},bnd,{}}; src{2} = {}; elem.src = src; geomdim = cell(1,2); geomdim{1} = {}; clear edg edg.map = {{'2'},{'2'},{'2'}}; edg.ind = {{'1'}}; geomdim{2} = {{},edg,{},{}}; elem.geomdim = geomdim; elem.var = {'w1','thx1','thy1'}; map = cell(1,2); clear submap submap.type = 'unit'; map{1} = submap; clear submap submap.type = 'linear'; submap.sg = '2'; submap.sv = {'2','1'}; submap.dg = '1'; submap.dv = {'4','3'}; map{2} = submap; elem.map = map; elemcpl{1} = elem; fem.elemcpl = elemcpl; xfem.fem{1} = fem; fem=xfem.fem{2}; % Application mode 1 clear appl appl.mode.class = 'Sme3DEulerBeam'; appl.module = 'SME'; appl.gporder = 6; appl.cporder = 1; appl.assignsuffix = '_smeul3d'; clear prop prop.analysis='eigen'; appl.prop = prop; clear edg edg.nu = 0.3; edg.Izz = '0.1^4/12'; edg.rho = 7800; edg.Hz = 1; edg.Rthy = 'thy1'; edg.Rthx = 'thx1'; edg.J = '2*0.1^4/12'; edg.Rz = 'w1'; edg.Iyy = '0.1^4/12'; edg.constrcond = 'norot'; edg.Hthx = 1; edg.Hthy = 1; edg.ind = [1]; appl.edg = edg; fem.appl{1} = appl; fem.frame = {'ref'}; fem.border = 1; clear units; units.basesystem = 'SI'; fem.units = units; % Coupling variable elements clear elemcpl % Extrusion coupling variables clear elem elem.elem = 'elcplextr'; elem.g = {'1','2'}; src = cell(1,2); src{1} = {}; clear edg edg.expr = {{'w'},{'thx'},{'thy'}}; edg.map = {{'1'},{'1'},{'1'}}; edg.ind = {{'1'}}; src{2} = {{},edg,{},{}}; elem.src = src; geomdim = cell(1,2); clear bnd bnd.map = {{{},'2'},{{},'2'},{{},'2'}}; bnd.ind = {{'1','2','3','5','6','7'},{'4'}}; geomdim{1} = {{},bnd,{}}; geomdim{2} = {}; elem.geomdim = geomdim; elem.var = {'w2','thx2','thy2'}; map = cell(1,2); clear submap submap.type = 'unit'; map{1} = submap; clear submap submap.type = 'linear'; submap.sg = '1'; submap.sv = {'4','3'}; submap.dg = '2'; submap.dv = {'2','1'}; map{2} = submap; elem.map = map; elemcpl{1} = elem; fem.elemcpl = elemcpl; xfem.fem{2} = fem; % ODE Settings clear ode clear units; units.basesystem = 'SI'; ode.units = units; xfem.ode=ode; % Multiphysics xfem=multiphysics(xfem); % Extend mesh xfem.xmesh=meshextend(xfem); % Solve problem xfem.sol=femeig(xfem, ... 'solcomp',{'thx','w','thy','thn_smdrm'}, ... 'outcomp',{'thx','w','thy','v','thz','u','thn_smdrm'}, ... 'rowscale','off', ... 'blocksize','auto', ... 'neigs',20); % Save current fem structure for restart purposes fem0=xfem; % Plot solution postplot(xfem, ... 'tridata',{'w','cont','internal','unit','m'}, ... 'triz',{'w','unit','m'}, ... 'trimap','jet(1024)', ... 'solnum',5, ... 'title','eigfreq_smdrm(5)=308.606761 Surface: z-displacement [m] Height: z-displacement [m]', ... 'grid','on'); % Plot solution postplot(xfem, ... 'tridata',{'w','cont','internal','unit','m'}, ... 'triz',{'w','unit','m'}, ... 'trimap','jet(1024)', ... 'solnum',5, ... 'title','eigfreq_smdrm(5)=308.606761 Surface: z-displacement [m] Height: z-displacement [m]', ... 'grid','on'); % Plot solution postplot(xfem, ... 'tridata',{'w','cont','internal','unit','m'}, ... 'triz',{'w','unit','m'}, ... 'trimap','jet(1024)', ... 'solnum',6, ... 'title','eigfreq_smdrm(6)=482.054525 Surface: z-displacement [m] Height: z-displacement [m]', ... 'grid','on'); % Plot solution postplot(xfem, ... 'tridata',{'w','cont','internal','unit','m'}, ... 'triz',{'w','unit','m'}, ... 'trimap','jet(1024)', ... 'solnum',7, ... 'title','eigfreq_smdrm(7)=560.584814 Surface: z-displacement [m] Height: z-displacement [m]', ... 'grid','on'); % (Default values are not included) fem=xfem.fem{1}; % Application mode 1 clear appl appl.mode.class = 'SmeMindlin'; appl.module = 'SME'; appl.gporder = 4; appl.cporder = 1; appl.assignsuffix = '_smdrm'; clear prop prop.analysis='eigen'; appl.prop = prop; clear bnd bnd.Rthy = {0,'thy2'}; bnd.Rthx = {0,'thx2'}; bnd.constrcoord = {'local','global'}; bnd.Hz = {0,1}; bnd.Rz = {0,'w2'}; bnd.Hthx = {0,1}; bnd.Hthy = {0,1}; bnd.ind = [1,1,1,2,1,1,1]; appl.bnd = bnd; clear equ equ.nu = 0.3; equ.thickness = 0.1; equ.rho = 7800; equ.ind = [1,1]; appl.equ = equ; fem.appl{1} = appl; fem.frame = {'ref'}; fem.border = 1; clear units; units.basesystem = 'SI'; fem.units = units; % Coupling variable elements clear elemcpl % Extrusion coupling variables clear elem elem.elem = 'elcplextr'; elem.g = {'1','2'}; src = cell(1,2); clear bnd bnd.expr = {{{},'w'},{{},'thx'},{{},'thy'}}; bnd.map = {{'1','1'},{'1','1'},{'1','1'}}; bnd.ind = {{'1','2','3','5','6','7'},{'4'}}; src{1} = {{},bnd,{}}; src{2} = {}; elem.src = src; geomdim = cell(1,2); geomdim{1} = {}; clear edg edg.map = {{'2'},{'2'},{'2'}}; edg.ind = {{'1'}}; geomdim{2} = {{},edg,{},{}}; elem.geomdim = geomdim; elem.var = {'w1','thx1','thy1'}; map = cell(1,2); clear submap submap.type = 'unit'; map{1} = submap; clear submap submap.type = 'linear'; submap.sg = '2'; submap.sv = {'2','1'}; submap.dg = '1'; submap.dv = {'4','3'}; map{2} = submap; elem.map = map; elemcpl{1} = elem; fem.elemcpl = elemcpl; xfem.fem{1} = fem; fem=xfem.fem{2}; % Application mode 1 clear appl appl.mode.class = 'Sme3DEulerBeam'; appl.module = 'SME'; appl.gporder = 6; appl.cporder = 1; appl.assignsuffix = '_smeul3d'; clear prop prop.analysis='eigen'; appl.prop = prop; clear edg edg.nu = 0.3; edg.Izz = '0.1^4/12'; edg.rho = 7800; edg.Hz = 1; edg.Rthy = 'thy1'; edg.Rthx = 'thx1'; edg.J = '2*0.1^4/12'; edg.Rz = 'w1'; edg.Iyy = '0.1^4/12'; edg.constrcond = 'norot'; edg.Hthx = 1; edg.Hthy = 1; edg.ind = [1]; appl.edg = edg; fem.appl{1} = appl; fem.frame = {'ref'}; fem.border = 1; clear units; units.basesystem = 'SI'; fem.units = units; % Coupling variable elements clear elemcpl % Extrusion coupling variables clear elem elem.elem = 'elcplextr'; elem.g = {'1','2'}; src = cell(1,2); src{1} = {}; clear edg edg.expr = {{'w'},{'thx'},{'thy'}}; edg.map = {{'1'},{'1'},{'1'}}; edg.ind = {{'1'}}; src{2} = {{},edg,{},{}}; elem.src = src; geomdim = cell(1,2); clear bnd bnd.map = {{{},'2'},{{},'2'},{{},'2'}}; bnd.ind = {{'1','2','3','5','6','7'},{'4'}}; geomdim{1} = {{},bnd,{}}; geomdim{2} = {}; elem.geomdim = geomdim; elem.var = {'w2','thx2','thy2'}; map = cell(1,2); clear submap submap.type = 'unit'; map{1} = submap; clear submap submap.type = 'linear'; submap.sg = '1'; submap.sv = {'4','3'}; submap.dg = '2'; submap.dv = {'2','1'}; map{2} = submap; elem.map = map; elemcpl{1} = elem; fem.elemcpl = elemcpl; xfem.fem{2} = fem; % ODE Settings clear ode clear units; units.basesystem = 'SI'; ode.units = units; xfem.ode=ode; % Multiphysics xfem=multiphysics(xfem); % Extend mesh xfem.xmesh=meshextend(xfem); % Solve problem xfem.sol=femeig(xfem, ... 'solcomp',{'w','thx','thy','v','thz','u','thn_smdrm'}, ... 'outcomp',{'w','thx','v','thy','u','thz','thn_smdrm'}, ... 'rowscale','off', ... 'blocksize','auto', ... 'neigs',20); % Save current fem structure for restart purposes fem0=xfem; % Plot solution postplot(xfem, ... 'tridata',{'w','cont','internal','unit','m'}, ... 'triz',{'w','unit','m'}, ... 'trimap','jet(1024)', ... 'solnum',1, ... 'title','eigfreq_smdrm(1)=0 (1) Surface: z-displacement [m] Height: z-displacement [m]', ... 'grid','on'); % Plot solution postplot(xfem, ... 'tridata',{'w','cont','internal','unit','m'}, ... 'triz',{'w','unit','m'}, ... 'trimap','jet(1024)', ... 'solnum',7, ... 'title','eigfreq_smdrm(7)=311.249003 Surface: z-displacement [m] Height: z-displacement [m]', ... 'grid','on'); % Plot solution postplot(xfem, ... 'tridata',{'w','cont','internal','unit','m'}, ... 'triz',{'w','unit','m'}, ... 'trimap','jet(1024)', ... 'solnum',8, ... 'title','eigfreq_smdrm(8)=485.349 Surface: z-displacement [m] Height: z-displacement [m]', ... 'grid','on'); % Plot solution postplot(xfem, ... 'tridata',{'w','cont','internal','unit','m'}, ... 'triz',{'w','unit','m'}, ... 'trimap','jet(1024)', ... 'solnum',9, ... 'title','eigfreq_smdrm(9)=585.097347 Surface: z-displacement [m] Height: z-displacement [m]', ... 'grid','on'); % (Default values are not included) fem=xfem.fem{1}; % Application mode 1 clear appl appl.mode.class = 'SmeMindlin'; appl.module = 'SME'; appl.gporder = 4; appl.cporder = 1; appl.assignsuffix = '_smdrm'; clear prop prop.analysis='eigen'; appl.prop = prop; clear bnd bnd.Rthy = {0,'thy2'}; bnd.Rthx = {0,'thx2'}; bnd.constrcoord = {'local','global'}; bnd.Hz = {0,1}; bnd.Rz = {0,'w2'}; bnd.Hthx = {0,1}; bnd.Hthy = {0,1}; bnd.ind = [1,1,1,2,1,1,1]; appl.bnd = bnd; clear equ equ.nu = 0.3; equ.thickness = 0.1; equ.rho = 7800; equ.ind = [1,1]; appl.equ = equ; fem.appl{1} = appl; fem.frame = {'ref'}; fem.border = 1; clear units; units.basesystem = 'SI'; fem.units = units; % Coupling variable elements clear elemcpl % Extrusion coupling variables clear elem elem.elem = 'elcplextr'; elem.g = {'1','2'}; src = cell(1,2); clear bnd bnd.expr = {{{},'w'},{{},'thx'},{{},'thy'}}; bnd.map = {{'1','1'},{'1','1'},{'1','1'}}; bnd.ind = {{'1','2','3','5','6','7'},{'4'}}; src{1} = {{},bnd,{}}; src{2} = {}; elem.src = src; geomdim = cell(1,2); geomdim{1} = {}; clear edg edg.map = {{'2'},{'2'},{'2'}}; edg.ind = {{'1'}}; geomdim{2} = {{},edg,{},{}}; elem.geomdim = geomdim; elem.var = {'w1','thx1','thy1'}; map = cell(1,2); clear submap submap.type = 'unit'; map{1} = submap; clear submap submap.type = 'linear'; submap.sg = '2'; submap.sv = {'2','1'}; submap.dg = '1'; submap.dv = {'4','3'}; map{2} = submap; elem.map = map; elemcpl{1} = elem; fem.elemcpl = elemcpl; xfem.fem{1} = fem; fem=xfem.fem{2}; % Application mode 1 clear appl appl.mode.class = 'Sme3DEulerBeam'; appl.module = 'SME'; appl.gporder = 6; appl.cporder = 1; appl.assignsuffix = '_smeul3d'; clear prop prop.analysis='eigen'; appl.prop = prop; clear edg edg.nu = 0.3; edg.Izz = '0.1^4/12'; edg.rho = 7800; edg.Hz = 1; edg.Rthy = 'thy1'; edg.Rthx = 'thx1'; edg.J = '2*0.1^4/12'; edg.Rz = 'w1'; edg.Iyy = '0.1^4/12'; edg.constrcond = 'displacement'; edg.Hthx = 1; edg.Hthy = 1; edg.ind = [1]; appl.edg = edg; fem.appl{1} = appl; fem.frame = {'ref'}; fem.border = 1; clear units; units.basesystem = 'SI'; fem.units = units; % Coupling variable elements clear elemcpl % Extrusion coupling variables clear elem elem.elem = 'elcplextr'; elem.g = {'1','2'}; src = cell(1,2); src{1} = {}; clear edg edg.expr = {{'w'},{'thx'},{'thy'}}; edg.map = {{'1'},{'1'},{'1'}}; edg.ind = {{'1'}}; src{2} = {{},edg,{},{}}; elem.src = src; geomdim = cell(1,2); clear bnd bnd.map = {{{},'2'},{{},'2'},{{},'2'}}; bnd.ind = {{'1','2','3','5','6','7'},{'4'}}; geomdim{1} = {{},bnd,{}}; geomdim{2} = {}; elem.geomdim = geomdim; elem.var = {'w2','thx2','thy2'}; map = cell(1,2); clear submap submap.type = 'unit'; map{1} = submap; clear submap submap.type = 'linear'; submap.sg = '1'; submap.sv = {'4','3'}; submap.dg = '2'; submap.dv = {'2','1'}; map{2} = submap; elem.map = map; elemcpl{1} = elem; fem.elemcpl = elemcpl; xfem.fem{2} = fem; % ODE Settings clear ode clear units; units.basesystem = 'SI'; ode.units = units; xfem.ode=ode; % Multiphysics xfem=multiphysics(xfem); % Extend mesh xfem.xmesh=meshextend(xfem); % Solve problem xfem.sol=femeig(xfem, ... 'solcomp',{'w','thx','v','thy','u','thz','thn_smdrm'}, ... 'outcomp',{'thx','w','thy','v','u','thz','thn_smdrm'}, ... 'rowscale','off', ... 'blocksize','auto', ... 'neigs',20); % Save current fem structure for restart purposes fem0=xfem; % (Default values are not included) fem=xfem.fem{1}; % Application mode 1 clear appl appl.mode.class = 'SmeMindlin'; appl.module = 'SME'; appl.gporder = 4; appl.cporder = 1; appl.assignsuffix = '_smdrm'; clear prop prop.analysis='eigen'; appl.prop = prop; clear bnd bnd.Rthy = {0,'thy2'}; bnd.Rthx = {0,'thx2'}; bnd.constrcoord = {'local','global'}; bnd.Hz = {0,1}; bnd.Rz = {0,'w2'}; bnd.Hthx = {0,1}; bnd.Hthy = {0,1}; bnd.ind = [1,1,1,2,1,1,1]; appl.bnd = bnd; clear equ equ.nu = 0.3; equ.thickness = 0.1; equ.rho = 7800; equ.ind = [1,1]; appl.equ = equ; fem.appl{1} = appl; fem.frame = {'ref'}; fem.border = 1; clear units; units.basesystem = 'SI'; fem.units = units; % Coupling variable elements clear elemcpl % Extrusion coupling variables clear elem elem.elem = 'elcplextr'; elem.g = {'1','2'}; src = cell(1,2); clear bnd bnd.expr = {{{},'w'},{{},'thx'},{{},'thy'}}; bnd.map = {{'1','1'},{'1','1'},{'1','1'}}; bnd.ind = {{'1','2','3','5','6','7'},{'4'}}; src{1} = {{},bnd,{}}; src{2} = {}; elem.src = src; geomdim = cell(1,2); geomdim{1} = {}; clear edg edg.map = {{'2'},{'2'},{'2'}}; edg.ind = {{'1'}}; geomdim{2} = {{},edg,{},{}}; elem.geomdim = geomdim; elem.var = {'w1','thx1','thy1'}; map = cell(1,2); clear submap submap.type = 'unit'; map{1} = submap; clear submap submap.type = 'linear'; submap.sg = '2'; submap.sv = {'2','1'}; submap.dg = '1'; submap.dv = {'4','3'}; map{2} = submap; elem.map = map; elemcpl{1} = elem; fem.elemcpl = elemcpl; xfem.fem{1} = fem; fem=xfem.fem{2}; % Application mode 1 clear appl appl.mode.class = 'Sme3DEulerBeam'; appl.module = 'SME'; appl.gporder = 6; appl.cporder = 1; appl.assignsuffix = '_smeul3d'; clear prop prop.analysis='eigen'; appl.prop = prop; clear edg edg.nu = 0.3; edg.Izz = '0.1^4/12'; edg.rho = 7800; edg.Hz = 1; edg.Rthy = 'thy1'; edg.Rthx = 'thx1'; edg.J = '2*0.1^4/12'; edg.Rz = 'w1'; edg.Iyy = '0.1^4/12'; edg.constrcond = 'displacement'; edg.Hthx = 1; edg.Hthy = 1; edg.ind = [1]; appl.edg = edg; fem.appl{1} = appl; fem.frame = {'ref'}; fem.border = 1; clear units; units.basesystem = 'SI'; fem.units = units; % Coupling variable elements clear elemcpl % Extrusion coupling variables clear elem elem.elem = 'elcplextr'; elem.g = {'1','2'}; src = cell(1,2); src{1} = {}; clear edg edg.expr = {{'w'},{'thx'},{'thy'}}; edg.map = {{'1'},{'1'},{'1'}}; edg.ind = {{'1'}}; src{2} = {{},edg,{},{}}; elem.src = src; geomdim = cell(1,2); clear bnd bnd.map = {{{},'2'},{{},'2'},{{},'2'}}; bnd.ind = {{'1','2','3','5','6','7'},{'4'}}; geomdim{1} = {{},bnd,{}}; geomdim{2} = {}; elem.geomdim = geomdim; elem.var = {'w2','thx2','thy2'}; map = cell(1,2); clear submap submap.type = 'unit'; map{1} = submap; clear submap submap.type = 'linear'; submap.sg = '1'; submap.sv = {'4','3'}; submap.dg = '2'; submap.dv = {'2','1'}; map{2} = submap; elem.map = map; elemcpl{1} = elem; fem.elemcpl = elemcpl; xfem.fem{2} = fem; % ODE Settings clear ode clear units; units.basesystem = 'SI'; ode.units = units; xfem.ode=ode; % Multiphysics xfem=multiphysics(xfem); % Extend mesh xfem.xmesh=meshextend(xfem); % Solve problem xfem.sol=femeig(xfem, ... 'solcomp',{'thx','w','thy','v','thz','u','thn_smdrm'}, ... 'outcomp',{'w','thx','v','thy','thz','u','thn_smdrm'}, ... 'rowscale','off', ... 'blocksize','auto', ... 'neigs',20); % Save current fem structure for restart purposes fem0=xfem; % Plot solution postplot(xfem, ... 'tridata',{'w','cont','internal','unit','m'}, ... 'triz',{'w','unit','m'}, ... 'trimap','jet(1024)', ... 'solnum',1, ... 'title','eigfreq_smdrm(1)=0 (1) Surface: z-displacement [m] Height: z-displacement [m]', ... 'grid','on'); % Plot solution postplot(xfem, ... 'tridata',{'w','cont','internal','unit','m'}, ... 'triz',{'w','unit','m'}, ... 'trimap','jet(1024)', ... 'solnum',9, ... 'title','eigfreq_smdrm(9)=237.447096 Surface: z-displacement [m] Height: z-displacement [m]', ... 'grid','on'); % Plot solution postplot(xfem, ... 'tridata',{'w','cont','internal','unit','m'}, ... 'triz',{'w','unit','m'}, ... 'trimap','jet(1024)', ... 'solnum',10, ... 'title','eigfreq_smdrm(10)=468.583656 (1) Surface: z-displacement [m] Height: z-displacement [m]', ... 'grid','on'); % Plot solution postplot(xfem, ... 'tridata',{'w','cont','internal','unit','m'}, ... 'triz',{'w','unit','m'}, ... 'trimap','jet(1024)', ... 'solnum',14, ... 'title','eigfreq_smdrm(14)=515.376445 Surface: z-displacement [m] Height: z-displacement [m]', ... 'grid','on'); % (Default values are not included) fem=xfem.fem{1}; % Application mode 1 clear appl appl.mode.class = 'SmeMindlin'; appl.module = 'SME'; appl.gporder = 4; appl.cporder = 1; appl.assignsuffix = '_smdrm'; clear prop prop.analysis='eigen'; appl.prop = prop; clear bnd bnd.Rthy = {0,'thy2'}; bnd.Rthx = {0,'thx2'}; bnd.constrcoord = {'local','global'}; bnd.Hz = {0,1}; bnd.Rz = {0,'w2'}; bnd.Hthx = {0,1}; bnd.Hthy = {0,1}; bnd.ind = [1,1,1,2,1,1,1]; appl.bnd = bnd; clear equ equ.nu = 0.3; equ.thickness = 0.1; equ.rho = 7800; equ.ind = [1,1]; appl.equ = equ; fem.appl{1} = appl; fem.frame = {'ref'}; fem.border = 1; clear units; units.basesystem = 'SI'; fem.units = units; % Coupling variable elements clear elemcpl % Extrusion coupling variables clear elem elem.elem = 'elcplextr'; elem.g = {'1','2'}; src = cell(1,2); clear bnd bnd.expr = {{{},'w'},{{},'thx'},{{},'thy'}}; bnd.map = {{'1','1'},{'1','1'},{'1','1'}}; bnd.ind = {{'1','2','3','5','6','7'},{'4'}}; src{1} = {{},bnd,{}}; src{2} = {}; elem.src = src; geomdim = cell(1,2); geomdim{1} = {}; clear edg edg.map = {{'2'},{'2'},{'2'}}; edg.ind = {{'1'}}; geomdim{2} = {{},edg,{},{}}; elem.geomdim = geomdim; elem.var = {'w1','thx1','thy1'}; map = cell(1,2); clear submap submap.type = 'unit'; map{1} = submap; clear submap submap.type = 'linear'; submap.sg = '2'; submap.sv = {'2','1'}; submap.dg = '1'; submap.dv = {'4','3'}; map{2} = submap; elem.map = map; elemcpl{1} = elem; fem.elemcpl = elemcpl; xfem.fem{1} = fem; fem=xfem.fem{2}; % Application mode 1 clear appl appl.mode.class = 'Sme3DEulerBeam'; appl.module = 'SME'; appl.gporder = 6; appl.cporder = 1; appl.assignsuffix = '_smeul3d'; clear prop prop.analysis='eigen'; appl.prop = prop; clear edg edg.nu = 0.3; edg.Izz = '0.1^4/12'; edg.rho = 7800; edg.Hz = 1; edg.Rthy = 'thy1'; edg.Rthx = 'thx1'; edg.J = '2*0.1^4/12'; edg.Rz = 'w1'; edg.Iyy = '0.1^4/12'; edg.constrcond = 'displacement'; edg.Hthx = 1; edg.Hthy = 1; edg.ind = [1]; appl.edg = edg; fem.appl{1} = appl; fem.frame = {'ref'}; fem.border = 1; clear units; units.basesystem = 'SI'; fem.units = units; % Coupling variable elements clear elemcpl % Extrusion coupling variables clear elem elem.elem = 'elcplextr'; elem.g = {'1','2'}; src = cell(1,2); src{1} = {}; src{2} = {{},{},{},{}}; elem.src = src; geomdim = cell(1,2); geomdim{1} = {{},{},{}}; geomdim{2} = {}; elem.geomdim = geomdim; elem.var = {}; map = cell(1,0); elem.map = map; elemcpl{1} = elem; fem.elemcpl = elemcpl; xfem.fem{2} = fem; % ODE Settings clear ode clear units; units.basesystem = 'SI'; ode.units = units; xfem.ode=ode; % Multiphysics xfem=multiphysics(xfem); % Extend mesh xfem.xmesh=meshextend(xfem); % Solve problem xfem.sol=femeig(xfem, ... 'solcomp',{'w','thx','thy','v','u','thz','thn_smdrm'}, ... 'outcomp',{'w','thx','v','thy','u','thz','thn_smdrm'}, ... 'rowscale','off', ... 'blocksize','auto', ... 'neigs',20); % Save current fem structure for restart purposes fem0=xfem; % (Default values are not included) fem=xfem.fem{1}; % Application mode 1 clear appl appl.mode.class = 'SmeMindlin'; appl.module = 'SME'; appl.gporder = 4; appl.cporder = 1; appl.assignsuffix = '_smdrm'; clear prop prop.analysis='eigen'; appl.prop = prop; clear bnd bnd.Rthy = {0,'thy2'}; bnd.Rthx = {0,'thx2'}; bnd.constrcoord = {'local','global'}; bnd.Hz = {0,1}; bnd.Rz = {0,'w2'}; bnd.Hthx = {0,1}; bnd.Hthy = {0,1}; bnd.ind = [1,1,1,2,1,1,1]; appl.bnd = bnd; clear equ equ.nu = 0.3; equ.thickness = 0.1; equ.rho = 7800; equ.ind = [1,1]; appl.equ = equ; fem.appl{1} = appl; fem.frame = {'ref'}; fem.border = 1; clear units; units.basesystem = 'SI'; fem.units = units; % Coupling variable elements clear elemcpl % Extrusion coupling variables clear elem elem.elem = 'elcplextr'; elem.g = {'1','2'}; src = cell(1,2); clear bnd bnd.expr = {{{},'w'},{{},'thx'},{{},'thy'}}; bnd.map = {{'1','1'},{'1','1'},{'1','1'}}; bnd.ind = {{'1','2','3','5','6','7'},{'4'}}; src{1} = {{},bnd,{}}; src{2} = {}; elem.src = src; geomdim = cell(1,2); geomdim{1} = {}; clear edg edg.map = {{'2'},{'2'},{'2'}}; edg.ind = {{'1'}}; geomdim{2} = {{},edg,{},{}}; elem.geomdim = geomdim; elem.var = {'w1','thx1','thy1'}; map = cell(1,2); clear submap submap.type = 'unit'; map{1} = submap; clear submap submap.type = 'linear'; submap.sg = '2'; submap.sv = {'2','1'}; submap.dg = '1'; submap.dv = {'4','3'}; map{2} = submap; elem.map = map; elemcpl{1} = elem; fem.elemcpl = elemcpl; xfem.fem{1} = fem; fem=xfem.fem{2}; % Application mode 1 clear appl appl.mode.class = 'Sme3DEulerBeam'; appl.module = 'SME'; appl.gporder = 6; appl.cporder = 1; appl.assignsuffix = '_smeul3d'; clear prop prop.analysis='eigen'; appl.prop = prop; clear edg edg.nu = 0.3; edg.Izz = '0.1^4/12'; edg.rho = 7800; edg.Hz = 1; edg.Rthy = 'thy1'; edg.Rthx = 'thx1'; edg.J = '2*0.1^4/12'; edg.Rz = 'w1'; edg.Iyy = '0.1^4/12'; edg.constrcond = 'displacement'; edg.Hthx = 1; edg.Hthy = 1; edg.ind = [1]; appl.edg = edg; fem.appl{1} = appl; fem.frame = {'ref'}; fem.border = 1; clear units; units.basesystem = 'SI'; fem.units = units; % Coupling variable elements clear elemcpl % Extrusion coupling variables clear elem elem.elem = 'elcplextr'; elem.g = {'1','2'}; src = cell(1,2); src{1} = {}; src{2} = {{},{},{},{}}; elem.src = src; geomdim = cell(1,2); geomdim{1} = {{},{},{}}; geomdim{2} = {}; elem.geomdim = geomdim; elem.var = {}; map = cell(1,0); elem.map = map; elemcpl{1} = elem; fem.elemcpl = elemcpl; xfem.fem{2} = fem; % ODE Settings clear ode clear units; units.basesystem = 'SI'; ode.units = units; xfem.ode=ode; % Multiphysics xfem=multiphysics(xfem); % Extend mesh xfem.xmesh=meshextend(xfem); % Solve problem xfem.sol=femeig(xfem, ... 'solcomp',{'thx','w','v','thy','thz','u','thn_smdrm'}, ... 'outcomp',{'thx','w','thy','v','thz','u','thn_smdrm'}, ... 'rowscale','off', ... 'blocksize','auto', ... 'neigs',20); % Save current fem structure for restart purposes fem0=xfem; % Plot solution postplot(xfem, ... 'tridata',{'w','cont','internal','unit','m'}, ... 'triz',{'w','unit','m'}, ... 'trimap','jet(1024)', ... 'solnum',1, ... 'title','eigfreq_smdrm(1)=0 (1) Surface: z-displacement [m] Height: z-displacement [m]', ... 'grid','on'); % Plot solution postplot(xfem, ... 'tridata',{'w','cont','internal','unit','m'}, ... 'triz',{'w','unit','m'}, ... 'trimap','jet(1024)', ... 'solnum',8, ... 'title','eigfreq_smdrm(8)=312.885349 Surface: z-displacement [m] Height: z-displacement [m]', ... 'grid','on'); % Plot solution postplot(xfem, ... 'tridata',{'w','cont','internal','unit','m'}, ... 'triz',{'w','unit','m'}, ... 'trimap','jet(1024)', ... 'solnum',9, ... 'title','eigfreq_smdrm(9)=417.090264 Surface: z-displacement [m] Height: z-displacement [m]', ... 'grid','on'); % Plot solution postplot(xfem, ... 'tridata',{'w','cont','internal','unit','m'}, ... 'triz',{'w','unit','m'}, ... 'trimap','jet(1024)', ... 'solnum',10, ... 'title','eigfreq_smdrm(10)=495.534267 Surface: z-displacement [m] Height: z-displacement [m]', ... 'grid','on'); % Plot solution postplot(xfem, ... 'tridata',{'w','cont','internal','unit','m'}, ... 'triz',{'w','unit','m'}, ... 'trimap','jet(1024)', ... 'solnum',8, ... 'title','eigfreq_smdrm(8)=312.885349 Surface: z-displacement [m] Height: z-displacement [m]', ... 'grid','on'); % Plot solution postplot(xfem, ... 'tridata',{'w','cont','internal','unit','m'}, ... 'triz',{'w','unit','m'}, ... 'trimap','jet(1024)', ... 'solnum',9, ... 'title','eigfreq_smdrm(9)=417.090264 Surface: z-displacement [m] Height: z-displacement [m]', ... 'grid','on');