clear; close all; clc; model = createpde('structural', 'static-solid'); importGeometry(model, 'piston.stl'); pdegplot(model); figure(); generateMesh(model, 'Hmax', 10, 'Hmin', 0.05); pdemesh(model); figure(); structuralProperties(model, 'YoungsModulus', 210E6, 'PoissonsRatio', 0.28); structuralBC(model, 'Face', [19 20], 'Constraint', 'fixed'); structuralBoundaryLoad(model, 'Face', 17, 'Pressure', 1); R = solve(model); pdeplot3D(model, 'ColorMapData', R.VonMisesStress);