// Gmsh project created on Fri Dec 07 15:55:29 2018

// #1, anisotropic material
// Material constants are arranged as:
// temperature density
// c11 c12 c13 c14 c15 c16
// c22 c23 c24 c25 c26
// c33 c34 c35 c36
// c44 c45 c46
// c55 c56
// c66
Physical Point("1 2 1
0 0
8.29 5.49 5.6 0 0 0
8.29 5.6 0 0 0
8.49 0 0 0
1.5 0 0
1.5 0
1.5
0 0 0 0 0 0 0
") = {};
// Material name: 1 -- MAT1

Point(1) = {-0.5, -0.5, 0, 1.0};
Point(2) = {0.5, -0.5, 0, 1.0};
Point(3) = {0.5, 0.5, 0, 1.0};
Point(4) = {-0.5, 0.5, 0, 1.0};
Point(5) = {-0.5, -0.5, 1, 1.0};
Point(6) = {0.5, -0.5, 1, 1.0};
Point(7) = {0.5, 0.5, 1, 1.0};
Point(8) = {-0.5, 0.5, 1, 1.0};
Line(2) = {1, 2};
Line(3) = {2, 3};
Line(4) = {3, 4};
Line(5) = {4, 1};
Line(6) = {5, 6};
Line(7) = {6, 7};
Line(8) = {7, 8};
Line(9) = {8, 5};
Line(10) = {1, 5};
Line(11) = {8, 4};
Line(12) = {2, 6};
Line(13) = {3, 7};
Point(9) = {0, -0.125, 0, 1.0};
Point(10) = {0.125, +0.125, 0, 1.0};
Point(11) = {-0.125, +0.125, 0, 1.0};
Point(12) = {0, -0.125, 1, 1.0};
Point(13) = {0.125, 0.125, 1, 1.0};
Point(14) = {-0.125, 0.125, 1, 1.0};
Line(14) = {9, 10};
Line(15) = {10, 11};
Line(16) = {11, 9};
Line(17) = {12, 13};
Line(18) = {13, 14};
Line(19) = {14, 12};
Line(20) = {9, 12};
Line(21) = {10, 13};
Line(22) = {11, 14};Line Loop(23) = {10, 6, -12, -2};
Plane Surface(24) = {23};
Line Loop(25) = {5, 10, -9, 11};
Plane Surface(26) = {25};
Line Loop(27) = {13, 8, 11, -4};
Plane Surface(28) = {27};
Line Loop(29) = {12, 7, -13, -3};
Plane Surface(30) = {29};
Line Loop(31) = {7, 8, 9, 6};
Line Loop(32) = {17, 18, 19};
Plane Surface(33) = {31, 32};
Line Loop(34) = {18, -22, -15, 21};
Plane Surface(35) = {34};
Line Loop(36) = {19, -20, -16, 22};
Plane Surface(37) = {36};
Line Loop(38) = {21, -17, -20, 14};
Plane Surface(39) = {38};
Line Loop(40) = {4, 5, 2, 3};
Line Loop(41) = {16, 14, 15};
Plane Surface(42) = {40, 41};
Surface Loop(43) = {33, 30, 24, 26, 42, 28, 35, 37, 39};
Volume(44) = {43};
Physical Volume(1) = {44};

Periodic Surface{30} = {26} Translate {1,0,0};
Periodic Surface{24} = {28} Translate {0,-1,0};
Periodic Surface{33} = {42} Translate {0,0,1};
