%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% ALL NEW %%%%%%%%%%%%%%%%%%%%%% fclose all; close all; clear all; disp('# starting'); i=0; lstart=6E-6; lstop =20.0E-6; lstep =.5E-6; num = (lstop-lstart)/lstep+1 out=[]; s1='C:\Users\Matthias Imboden\Documents\Research_Mohanty_Lab\Dissipation Review Article\ClampingLosses\'; s2='output_3DModel_6-20_um_vfine_smallBase.dat'; sn=[s1,s2]; for l1 = [lstart: lstep: lstop] % COMSOL Multiphysics Model M-file % Generated by COMSOL 3.5a (COMSOL 3.5.0.603, $Date: 2008/12/03 17:02:19 $) flclear fem % COMSOL version clear vrsn vrsn.name = 'COMSOL 3.5'; vrsn.ext = 'a'; vrsn.major = 0; vrsn.build = 603; vrsn.rcs = '$Name: $'; vrsn.date = '$Date: 2008/12/03 17:02:19 $'; fem.version = vrsn; % Geometry g1=block3('1.5e-6','1.5e-6','1e-6','base','center','pos',{'0','0','0'},'axis',{'0','0','1'},'rot','0'); g1=move(g1,[-0.75e-6,0,-0.5e-6]); g2=block3(l1,'.5e-6','.3e-6','base','center','pos',{'0','0','0'},'axis',{'0','0','1'},'rot','0'); g2=move(g2,[l1/2,0,-.15e-6]); % Analyzed geometry clear s s.objs={g1,g2}; s.name={'BLK1','BLK2'}; s.tags={'g1','g2'}; fem.draw=struct('s',s); fem.geom=geomcsg(fem); % Initialize mesh fem.mesh=meshinit(fem, ... 'hauto',5); % Refine mesh fem.mesh=meshrefine(fem, ... 'mcase',0, ... 'rmethod','longest'); % Refine mesh fem.mesh=meshrefine(fem, ... 'mcase',0, ... 'rmethod','longest'); % Refine mesh fem.mesh=meshrefine(fem, ... 'mcase',0, ... 'rmethod','longest'); %%%%%%%%%%%%%%%%%%%%%%%%% fine mesh % Refine mesh fem.mesh=meshrefine(fem, ... 'mcase',0, ... 'rmethod','longest'); % Refine mesh %fem.mesh=meshrefine(fem, ... % 'mcase',0, ... % 'rmethod','longest'); % (Default values are not included) % Application mode 1 clear appl appl.mode.class = 'SmeSolid3'; appl.module = 'MEMS'; appl.gporder = 4; appl.cporder = 2; appl.sshape = 2; appl.assignsuffix = '_smsld'; clear prop prop.analysis='eigen'; appl.prop = prop; clear bnd bnd.constrcond = {'free','fixed'}; bnd.ind = [2,2,2,1,2,1,1,1,1,1,1,1]; appl.bnd = bnd; clear equ equ.nu = .2; equ.rho = 3500; equ.E = 900e9; 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; % 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','v','u'}, ... 'outcomp',{'w','v','u'}, ... 'rowscale','off', ... 'blocksize','auto', ... 'neigs',2, ... 'linsolver','spooles'); % Save current fem structure for restart purposes fem0=fem; % Plot solution %postplot(fem, ... % 'slicedata',{'mises_smsld','cont','internal','unit','Pa'}, ... % 'slicexspacing',5, ... % 'sliceyspacing',0, ... % 'slicezspacing',0, ... % 'slicemap','Rainbow', ... % 'solnum',1, ... % 'title','eigfreq_smsld(1)=9.748376e6 Slice: von Mises stress [Pa]', ... % 'grid','on', ... % 'campos',[-3.734829576007982E-5,-5.114937476534861E-5,3.5973143604420465E-5], ... % 'camtarget',[1.899999915622175E-6,0,-1.249999968422344E-6], ... % 'camup',[0,0,1], ... % 'camva',7.219154077388316); % Plot solution %postplot(fem, ... % 'tetdata',{'mises_smsld','cont','internal','unit','Pa'}, ... % 'tetmap','Rainbow', ... % 'tetkeep',1, ... % 'tetkeeptype','random', ... % 'deformsub',{'u','v','w'}, ... % 'solnum',1, ... % 'title','eigfreq_smsld(1)=9.748376e6 Subdomain: von Mises stress [Pa] Deformation: Displacement', ... % 'grid','on', ... % 'campos',[-3.734829576007982E-5,-5.114937476534861E-5,3.5973143604420465E-5], ... % 'camtarget',[1.899999915622175E-6,0,-1.249999968422344E-6], ... % 'camup',[0,0,1], ... % 'camva',7.2191540773884295); % Plot solution %postplot(fem, ... % 'tetdata',{'mises_smsld','cont','internal','unit','Pa'}, ... % 'tetmap','Rainbow', ... % 'tetkeep',1, ... % 'tetkeeptype','random', ... % 'deformsub',{'u','v','w'}, ... % 'solnum','end', ... % 'title','eigfreq_smsld(2)=1.61479e7 Subdomain: von Mises stress [Pa] Deformation: Displacement', ... % 'grid','on', ... % 'campos',[-3.734829576007982E-5,-5.114937476534861E-5,3.5973143604420465E-5], ... % 'camtarget',[1.899999915622175E-6,0,-1.249999968422344E-6], ... % 'camup',[0,0,1], ... % 'camva',7.2191540773884295); % Plot solution %postplot(fem, ... % 'tetdata',{'mises_smsld','cont','internal','unit','Pa'}, ... % 'tetmap','Rainbow', ... % 'tetkeep',1, ... % 'tetkeeptype','random', ... % 'deformsub',{'u','v','w'}, ... % 'solnum','end', ... % 'title','eigfreq_smsld(2)=1.61479e7 Subdomain: von Mises stress [Pa] Deformation: Displacement', ... % 'grid','on', ... % 'campos',[-3.734829576007982E-5,-5.114937476534861E-5,3.5973143604420465E-5], ... % 'camtarget',[1.899999915622175E-6,0,-1.249999968422344E-6], ... % 'camup',[0,0,1], ... % 'camva',7.2191540773884295); % Integrate I1=postint(fem,'Ws_smsld', ... 'unit','J', ... 'recover','off', ... 'dl',2, ... 'solnum',1); % Integrate I2=postint(fem,'Ws_smsld', ... 'unit','J', ... 'recover','off', ... 'dl',1, ... 'solnum',1); % Integrate I3=postint(fem,'Ws_smsld', ... 'unit','J', ... 'recover','off', ... 'dl',2, ... 'solnum','end'); % Integrate I4=postint(fem,'Ws_smsld', ... 'unit','J', ... 'recover','off', ... 'dl',1, ... 'solnum','end'); l1 %f0=postmin(fem,'eigfreq_smsld') f0=-imag(fem.sol.lambda(1))/2/pi f1=-imag(fem.sol.lambda(2))/2/pi out(i+1,1)=l1; % beam length ? out(i+1,2)=f0; % out of plane resoannce frequnecy out(i+1,3)=I1; % beam strain energy out(i+1,4)=I2; % base strain energy out(i+1,5)=I2/(I1+I2); % loss factor out(i+1,6)=f1; % in plane resoannce frequnecy out(i+1,7)=I3; % beam strain energy out(i+1,8)=I4; % base strain energy out(i+1,9)=I4/(I3+I4); % loss factor out(i+1,10)=1/l1; % 1/L out(i+1,11)=1/l1^2; % 1/L^2 (to see if frequency scales) out(i+1,12)=1/l1^3; % 1/L^3 out(i+1,13)=1/l1^5; % 1/L^5 (to see if dissipation scales) i=i+1; disp('% completed:'); i/num*100 end; save(sn,'out','-ASCII','-double'); %% ending disp('# ending');