No subject

[Nicolas MOES]

Voila le geo qui pose probleme.

ta1 = 0.5;
ta2 = 2;
nbpt1 = 6;
nbpt2 = 12;
a = 1.0;
L = 5.0;
angle = Pi/4.0;
p = a*Cos(angle);
normal = -5;

Point(1) = {0.0,0.0,0.0,ta1};
Point(2) = {a,0.0,0.0,ta1};
Point(3) = {p,p,0.,ta1};
Point(4) = {0.0,a,0.0,ta1};
Point(5) = {L,0.0,0.0,ta2};
Point(6) = {L,L,0.0,ta2};
Point(7) = {0.0,L,0.0,ta2};

Line(1) = {1,2};
Line(2) = {2,5};
Line(3) = {5,6};
Line(4) = {6,7}; 
Line(5) = {7,4}; 
Line(6) = {4,1}; 
Line(7) = {1,3};
Line(8) = {3,6};
Circle(9) = {2,1,3};
Circle(10) = {3,1,4}; 

Line Loop(11) = {-9,-1,7};
Plane Surface(12) = {11};
Line Loop(13) = {8,-3,-2,9};
Plane Surface(14) = {13};
Line Loop(15) = {4,5,-10,8};
Plane Surface(16) = {15};
Line Loop(17) = {6,7,10};
Plane Surface(18) = {17};

Transfinite Line{1,7} = nbpt1 Using Power 1.0;
Transfinite Line{9} = nbpt1 Using Power 1.0;
Transfinite Surface{12} = {1,2,3};

Transfinite Line{2,8} = nbpt1 Using Power 1.0;
Transfinite Line{3,9} = nbpt2 Using Power 1.0;
Transfinite Surface{14} = {2,5,6,3};

Transfinite Line{5,8} = nbpt1 Using Power 1.0;
Transfinite Line{4,10} = nbpt2 Using Power 1.0;
Transfinite Surface{16} = {4,3,6,7};

Transfinite Line{6,7} = nbpt1 Using Power 1.0;
Transfinite Line{10} = nbpt2 Using Power 1.0;
Transfinite Surface{18} = {1,3,4};

Recombine Surface{12,18,16,14};

/* Extrude surface 12 */

Extrude(12,{0,0,normal});

Transfinite Line{23,27,24,7,1,21,20} = nbpt1 Using Power 1.0;
Transfinite Line{9,19} = nbpt2 Using Power 1.0;
Transfinite Surface{26} = {2,8,10,3};
Transfinite Surface{12} = {1,2,3};
Transfinite Surface{31} = {1,3,9,10};
Transfinite Surface{29} = {2,1,9,8};
Transfinite Surface{32} = {9,8,10};
Recombine Surface{29,26,31,32,12};

Surface Loop(1000) = {31,-26,32,-29,-12};
Complex Volume(2000) = {1000};

/* subtlety: first and fourth must be the degenerate ones for */
/* the prismatic hexahedral element (both faces in ccw orientation) */

Transfinite Volume{2000} = {9,8,10,1,2,3};

Extrude(18,{0,0,normal});

Transfinite Line{27,2001,21,24,2005,6,7} = nbpt1 Using Power 1.0;
Transfinite Line{2003,10} = nbpt2 Using Power 1.0;
Transfinite Surface{2012} = {9,10,11};
Transfinite Surface{18} = {1,3,4};
Transfinite Surface{2007} = {9,1,4,11};
Transfinite Surface{31} = {9,1,3,10};
Transfinite Surface{2011} = {11,10,3,4};
Recombine Surface{2012,18,2011,2007,31};

Surface Loop(2013) = {31,-18,2007,2012,2011};
Complex Volume(3000) = {2013};

/* subtlety: first and fourth must be the degenerate ones for */
/* the prismatic hexahedral element (both faces in ccw orientation) */

Transfinite Volume{3000} = {9,10,11,1,3,4};

Extrude(14,{0,0,normal});

Transfinite Line{2,3009,3003,23,24,8,3006,3001} = nbpt1 Using Power 1.0;
Transfinite Line{9,3,19,3002} = nbpt2 Using Power 1.0;
Transfinite Surface{26} = {2,8,10,3};
Transfinite Surface{3013} = {2,5,14,8};
Transfinite Surface{14} = {3,2,5,6};
Transfinite Surface{3016} = {10,8,14,13};
Transfinite Surface{3008} = {3,6,13,10};
Transfinite Surface{3011} = {5,14,13,6};

Recombine Surface{3013,3016,26,3008,14,3011};

Surface Loop(3017) = {26,-14,3008,-3011,3016,-3013};
Complex Volume(3018) = {3017};
Transfinite Volume{3018} = {2,5,14,8,3,6,13,10};

Extrude(16,{0,0,normal});

Transfinite Line{5,3020,8,3024,3006,3001,24,2005} = nbpt1 Using Power 1.0;
Transfinite Line{2003,10,3019,4} = nbpt2 Using Power 1.0;

Transfinite Surface{3028} = {4,11,16,7};
Transfinite Surface{3026} = {7,16,13,6};
Transfinite Surface{3008} = {10,3,6,13};
Transfinite Surface{2011} = {4,3,10,11};
Transfinite Surface{3033} = {10,13,16,11};
Transfinite Surface{16} = {4,3,6,7};

Recombine Surface{3008,3026,3028,2011,3033,16};

Surface Loop(3034) = {2011,16,-3026,-3028,-3033,-3008};
Complex Volume(3035) = {3034};
Transfinite Volume{3035} = {4,3,10,11,7,6,13,16};

/* Physical Definitions */

/* 1/8th Symmetry */
/* 12,14,16,18 Top surface z = 0 (Traction +1) */
/* 3016,3033 Part of bottom surface z = -5 (u_z = 0) */
/* 3013,29 Symmetry surface u_y = 0 */
/* 3028,2007 Symmetry surface u_x = 0 */

Physical Surface(4000) = {12,14,16,18,3016,3033,3013,29,3028,2007};
Physical Surface(4001) = {32,2012}; /* crack surface on z = -5 */
Physical Volume(5000) = {2000,3000,3018,3035};
Recombine Surface{2012,18,2007,2011,31};