function out = model % % pde2.m % % Model exported on Apr 1 2021, 15:10 by COMSOL 5.4.0.295. import com.comsol.model.* import com.comsol.model.util.* model = ModelUtil.create('Model'); model.modelPath('E:\academy\phononic_crystals\simulation\example'); model.label('pde2.mph'); model.param.set('lambda1', '1.23e11'); model.param.set('mu1', '8.2e10'); model.param.set('rho1', '7630'); model.param.set('freq1', '1'); model.param.set('A', '1'); model.param.set('lambda2', '4.213e9'); model.param.set('mu2', '1.331e9'); model.param.set('rho2', '1100'); model.component.create('comp1', true); model.component('comp1').geom.create('geom1', 2); model.component('comp1').mesh.create('mesh1'); model.component('comp1').geom('geom1').create('sq1', 'Square'); model.component('comp1').geom('geom1').feature('sq1').set('size', 'A'); model.component('comp1').geom('geom1').run; model.component('comp1').physics.create('c', 'CoefficientFormPDE', 'geom1'); model.component('comp1').physics('c').field('dimensionless').component({'u1'}); model.component('comp1').mesh('mesh1').create('ftri1', 'FreeTri'); model.component('comp1').view('view1').axis.set('xmin', -0.1039215624332428); model.component('comp1').view('view1').axis.set('xmax', 1.1039215326309204); model.component('comp1').view('view1').axis.set('ymin', -0.04999998211860657); model.component('comp1').view('view1').axis.set('ymax', 1.0499999523162842); model.component('comp1').physics('c').feature('cfeq1').set('a', 1); model.component('comp1').physics('c').feature('cfeq1').set('ea', 1); model.component('comp1').physics('c').feature('cfeq1').set('al', [1 1]); model.component('comp1').physics('c').feature('cfeq1').set('be', [1 1]); model.component('comp1').physics('c').feature('cfeq1').set('ga', [1 1]); model.component('comp1').mesh('mesh1').feature('size').set('hauto', 3); model.component('comp1').mesh('mesh1').run; model.study.create('std1'); model.study('std1').create('eigv', 'Eigenvalue'); model.sol.create('sol1'); model.sol('sol1').study('std1'); model.sol('sol1').attach('std1'); model.sol('sol1').create('st1', 'StudyStep'); model.sol('sol1').create('v1', 'Variables'); model.sol('sol1').create('e1', 'Eigenvalue'); model.sol('sol1').feature('e1').create('im1', 'InputMatrix'); model.sol('sol1').attach('std1'); model.sol('sol1').feature('e1').feature('im1').set('K', true); model.sol('sol1').runAll; out = model;