function out = model
%
% single_solid_mechanics_physics.m
%
% Model exported on Oct 28 2019, 20:05 by COMSOL 5.4.0.246.

import com.comsol.model.*
import com.comsol.model.util.*

model = ModelUtil.create('Model');

model.modelPath('C:\Users\sayinde\Desktop\exampletest');

model.label('single_solid_mechanics_physics.mph');

model.component.create('comp1', true);

model.component('comp1').geom.create('geom1', 2);

model.component('comp1').mesh.create('mesh1');

model.component('comp1').geom('geom1').create('r1', 'Rectangle');
model.component('comp1').geom('geom1').feature('r1').set('pos', [-0.45 0]);
model.component('comp1').geom('geom1').feature('r1').set('size', [1 0.5]);
model.component('comp1').geom('geom1').run;

model.component('comp1').material.create('mat1', 'Common');
model.component('comp1').material('mat1').propertyGroup.create('Enu', 'Young''s modulus and Poisson''s ratio');

model.component('comp1').physics.create('solid', 'SolidMechanics', 'geom1');
model.component('comp1').physics('solid').create('fix1', 'Fixed', 1);
model.component('comp1').physics('solid').feature('fix1').selection.set([1]);
model.component('comp1').physics('solid').create('pl1', 'PointLoad', 0);
model.component('comp1').physics('solid').feature('pl1').selection.set([4]);

model.component('comp1').mesh('mesh1').create('map1', 'Map');
model.component('comp1').mesh('mesh1').feature('map1').create('dis1', 'Distribution');
model.component('comp1').mesh('mesh1').feature('map1').create('dis2', 'Distribution');
model.component('comp1').mesh('mesh1').feature('map1').feature('dis1').selection.set([2 3]);
model.component('comp1').mesh('mesh1').feature('map1').feature('dis2').selection.set([1 4]);

model.component('comp1').view('view1').axis.set('xmin', -0.47002750635147095);
model.component('comp1').view('view1').axis.set('xmax', 0.6031083464622498);
model.component('comp1').view('view1').axis.set('ymin', -0.2756078243255615);
model.component('comp1').view('view1').axis.set('ymax', 0.709986686706543);

model.component('comp1').material('mat1').propertyGroup('Enu').set('youngsmodulus', '3');
model.component('comp1').material('mat1').propertyGroup('Enu').set('poissonsratio', '0.3');

model.component('comp1').physics('solid').prop('ShapeProperty').set('order_displacement', 1);
model.component('comp1').physics('solid').feature('pl1').set('Fp', [0; -10; 0]);

model.component('comp1').mesh('mesh1').feature('map1').feature('dis1').set('numelem', 2);
model.component('comp1').mesh('mesh1').feature('map1').feature('dis2').set('numelem', 1);
model.component('comp1').mesh('mesh1').run;

model.study.create('std1');
model.study('std1').create('stat', 'Stationary');

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('s1', 'Stationary');
model.sol('sol1').feature('s1').create('fc1', 'FullyCoupled');
model.sol('sol1').feature('s1').feature.remove('fcDef');

model.result.create('pg1', 'PlotGroup2D');
model.result('pg1').create('surf1', 'Surface');
model.result('pg1').feature('surf1').create('def', 'Deform');
model.result.export.create('data1', 'Data');

model.sol('sol1').attach('std1');
model.sol('sol1').feature('s1').feature('aDef').set('cachepattern', true);
model.sol('sol1').runAll;

model.result('pg1').label('Stress (solid)');
model.result('pg1').feature('surf1').set('expr', 'u');
model.result('pg1').feature('surf1').set('descr', 'Displacement field, X component');
model.result('pg1').feature('surf1').set('const', {'solid.refpntx' '0' 'Reference point for moment computation, x coordinate'; 'solid.refpnty' '0' 'Reference point for moment computation, y coordinate'; 'solid.refpntz' '0' 'Reference point for moment computation, z coordinate'});
model.result('pg1').feature('surf1').set('colortable', 'RainbowLight');
model.result('pg1').feature('surf1').set('resolution', 'normal');
model.result('pg1').feature('surf1').feature('def').set('scale', 9.286117660543049E-4);
model.result('pg1').feature('surf1').feature('def').set('scaleactive', false);
model.result.export('data1').set('expr', {'u'});
model.result.export('data1').set('unit', {'m'});
model.result.export('data1').set('descr', {'Displacement field, X component'});
model.result.export('data1').set('const', {'solid.refpntx' '0' 'Reference point for moment computation, x coordinate'; 'solid.refpnty' '0' 'Reference point for moment computation, y coordinate'; 'solid.refpntz' '0' 'Reference point for moment computation, z coordinate'});


 U = mphgetu(model, 'soltag', 'sol1');
 
save result


out = model;