structuralmodel = createpde('structural','transient-solid');

gm = multicuboid(0.06,0.005,0.01);
structuralmodel.Geometry = gm;
pdegplot(structuralmodel,'FaceLabels','on','FaceAlpha',0.5)
view(50,20)

structuralProperties(structuralmodel,'YoungsModulus',210E9, ...
                                     'PoissonsRatio',0.3, ...
                                     'MassDensity',7800);

structuralBC(structuralmodel,'Face',5,'Constraint','fixed');

structuralBC(structuralmodel,'Face',3, ...
                             'YDisplacement',1E-4, ...
                             'Frequency',50);

generateMesh(structuralmodel,'Hmax',0.01);

structuralIC(structuralmodel,'Displacement',[0,0,0],'Velocity',[0,0,0]);

tlist = 0:0.002:0.2;
structuralresults = solve(structuralmodel,tlist);

stress = evaluateStress(structuralresults);

figure
pdeplot3D(structuralmodel,'ColorMapData',stress.sxx(:,end))
title('x-Direction Normal Stress in the Beam of the Last Time-Step')