boundary.geo

efg code with matlab

// density
h = 0.05;

Point(1) = {0,0,0,h};
Point(2) = {1,0,0,h};
Point(3) = {1,1,0,h};
Point(4) = {0,1,0,h};
Line(1) = {1,2};
Line(2) = {2,3};
Line(3) = {3,4};
Line(4) = {4,1};
Line Loop(5) = {1,2,3,4};
Plane Surface(1000) = {5};

// define physical entities for definition of displacement and force boundary conditions
Physical Line(20)={1}; // elements in the displacement boundary 
Physical Line(21)={4}; // elements in the traction boundary 
Physical Line(22)={3}; // elements in the traction boundary 

Physical Surface(30)={1000}; // elements in the interior domain