%
% BlockSinWave.m
%
% Model exported on Jan 26 2012, 10:52 by COMSOL 4.2.0.187.

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

model = ModelUtil.create('Model');

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

model.name('Block1.mph');

model.param.set('W0', '66[Hz]');

model.modelNode.create('mod1');

model.geom.create('geom1', 3);
model.geom('geom1').feature.create('blk1', 'Block');
model.geom('geom1').feature('blk1').set('size', {'3[mm]' '2[mm]' '2[mm]'});
model.geom('geom1').run;

model.material.create('mat1');

model.physics.create('solid', 'SolidMechanics', 'geom1');
model.physics('solid').feature.create('disp1', 'Displacement2', 2);
model.physics('solid').feature('disp1').selection.set([1]);

model.mesh.create('mesh1', 'geom1');
model.mesh('mesh1').feature.create('ftet1', 'FreeTet');

model.material('mat1').propertyGroup('def').set('youngsmodulus', '2e11');
model.material('mat1').propertyGroup('def').set('poissonsratio', '.29');
model.material('mat1').propertyGroup('def').set('density', '7850');

model.physics('solid').feature('disp1').set('U0', {'0.00001*sin(W0*2*pi*t)'; '0'; '0'});
model.physics('solid').feature('disp1').set('Direction', {'1'; '0'; '0'});

model.mesh('mesh1').feature('size').set('hauto', 6);
model.mesh('mesh1').run;

model.study.create('std1');
model.study('std1').feature.create('time', 'Transient');

model.sol.create('sol1');
model.sol('sol1').study('std1');
model.sol('sol1').attach('std1');
model.sol('sol1').feature.create('st1', 'StudyStep');
model.sol('sol1').feature.create('v1', 'Variables');
model.sol('sol1').feature.create('t1', 'Time');
model.sol('sol1').feature('t1').feature.create('fc1', 'FullyCoupled');
model.sol('sol1').feature('t1').feature.remove('fcDef');

model.result.dataset.create('cpt1', 'CutPoint3D');

model.study('std1').feature('time').set('tlist', 'range(0,0.0001,4/W0)');

model.sol('sol1').feature('st1').name('Compile Equations: Time Dependent');
model.sol('sol1').feature('st1').set('studystep', 'time');
model.sol('sol1').feature('v1').set('control', 'time');
model.sol('sol1').feature('v1').feature('mod1_u').set('scalemethod', 'manual');
model.sol('sol1').feature('v1').feature('mod1_u').set('scaleval', '1e-2*0.00412310562561766');
model.sol('sol1').feature('t1').set('tlist', 'range(0,0.0001,4/W0)');
model.sol('sol1').feature('t1').set('timemethod', 'genalpha');
model.sol('sol1').feature('t1').set('tstepsgenalpha', 'strict');
model.sol('sol1').runAll;

model.result.dataset('cpt1').set('pointx', '0');
model.result.dataset('cpt1').set('pointy', '1[mm]');
model.result.dataset('cpt1').set('pointz', '1[mm]');

out = model;

 A=mpheval(model,{'u','t'},'edim',2,'selection',4)
 uu=A.d1(:,13);
 tt=A.d2(:,13);
 hold on
 plot(tt,uu)