Skip to content
Snippets Groups Projects
Commit 7a1bcc72 authored by Christophe Geuzaine's avatar Christophe Geuzaine
Browse files

obsolete

parent 609b6bfb
No related branches found
No related tags found
No related merge requests found
// Ce fichier est la description geometrique utilisee par GMSH pour
// Ce fichier est la description geometrique utilisee par GMSH pour
// generer le fichier "four.msh".
//
// Copyright (C) 2003 Ali ABAKAR <abakar@edf.fr>
......@@ -17,10 +17,15 @@ ProfondeurMarche = 1. ; LongueurMarche = 2. ;
Include "Param_test.geo";
// Densité de maillage (taille des mailles autour d'un point géomètrique)
lfour = LongueurFour/NbElLongF ;
// Description géomètrique du four
// Création des points
pf1 = newp; Point(pf1) = { LongueurFour, LargeurFour/2., -ProfondeurMarche, lfour};
pf2 = newp; Point(pf2) = { LongueurFour-LongueurMarche, LargeurFour/2., -ProfondeurMarche, lfour};
pf3 = newp; Point(pf3) = { LongueurFour-LongueurMarche, LargeurFour/2., 0., lfour};
......@@ -70,13 +75,13 @@ lf18 = newreg; Line(lf18) = {pf5, pf12};
// Définition des contour(s) du dessus du four
DessusFour[0] = newreg; Line Loop(DessusFour[0]) = {lf9,lf10,lf11,lf12};
DessusFour[0] = newreg; Line Loop(DessusFour[0]) = {lf9,lf10,lf11,lf12};
// Définitions des contours des bords du four
// Bord en y = LargeurFour/2.
BordFour[0] = newreg; Line Loop(BordFour[0]) = {lf1,lf14,lf5,lf15,-lf9, -lf13};
BordFour[0] = newreg; Line Loop(BordFour[0]) = {lf1,lf14,lf5,lf15,-lf9, -lf13};
// Bord en x = 0.
......@@ -88,7 +93,7 @@ BordFour[2] = newreg; Line Loop(BordFour[2]) = {lf7,-lf17,lf3,lf18,-lf11,-lf16};
// Bord en x = LongueurFour
BordFour[3] = newreg; Line Loop(BordFour[3]) = {lf4, lf13,-lf12,-lf18};
BordFour[3] = newreg; Line Loop(BordFour[3]) = {lf4, lf13,-lf12,-lf18};
// Définitions des contour(s) du fond du Four
......@@ -118,7 +123,7 @@ Include "Electrode_8.geo" ;
// Rayon et Hauteur Electrode
r = 0.076/2.; l = 0.35;
r = 0.076/2.; l = 0.35;
// Vecteur axe electrode
......@@ -139,7 +144,7 @@ Call Electrode;
t = t + 1 ;
R = 0.5 ;
L = 0.8 ;
yc = 3.5 ; zc = 0 ;
yc = 3.5 ; zc = 0 ;
Call Electrode;
......@@ -167,7 +172,7 @@ Call Electrode;
t = t + 1 ;
R = 0.5 ;
L = 0.8 ;
yc = 3.5 ; zc = 0 ;
yc = 3.5 ; zc = 0 ;
Call Electrode;
......@@ -195,7 +200,7 @@ Call Electrode;
t = t + 1 ;
R = 0.5 ;
L = 0.8 ;
yc = 3.5 ; zc = 0 ;
yc = 3.5 ; zc = 0 ;
Call Electrode;
......@@ -214,7 +219,7 @@ xc = xc + 2.5 ;
// ELECTRODE I
t = t+1;
R = 1. ;
R = 1. ;
yc = 0 ; zc = 0 ;
Call Electrode;
......@@ -224,7 +229,7 @@ Call Electrode;
t = t + 1 ;
R = 0.5 ;
L = 0.8 ;
yc = 3.5 ; zc = 0 ;
yc = 3.5 ; zc = 0 ;
Call Electrode;
......@@ -253,7 +258,7 @@ Call Electrode;
// ELECTRODE II
t = t + 1 ;
yc = 0.9 ; zc = 0 ;
yc = 0.9 ; zc = 0 ;
Call Electrode;
......@@ -314,7 +319,7 @@ Call Electrode;
Plane Surface(FondFour[0]) = { FondFour[0]};
Plane Surface(FondFour[1]) = { FondFour[1],ContBaseBoiteElect[]};
Plane Surface(FondFour[2]) = { FondFour[2]};
Plane Surface(FondFour[2]) = { FondFour[2]};
// Surface dessus du Four - trace(s) Electrode(s)
......@@ -331,7 +336,7 @@ Plane Surface(BordFour[3]) = { BordFour[3]};
// -------------------------------------------------
FrontExtFour = newreg;
Surface Loop(FrontExtFour) = {FondFour[],DessusFour[],BordFour[]};
Surface Loop(FrontExtFour) = {FondFour[],DessusFour[],BordFour[]};
Verre = newreg ;
Volume(Verre) = {FrontExtFour,FrontElect[]};
......
......@@ -14,8 +14,3 @@ NbBoxSup = 3.;
NbBoxInf = 3.;
NbBoxGen = 2.;
// Densité de maillage (taille des mailles autour d'un point géomètrique)
lfour = LongueurFour/NbElLongF ;
// Extrusion and rotation. Christophe Geuzaine, Feb 2002.
lc = 0.005;
lc2 = 4*lc;
Point(1)={0.09,0,0,lc};
Point(2)={0.11,0,0,lc};
Point(3)={0,0.09,0,lc};
Point(4)={0,0.11,0,lc};
Point(5)={0.2,0.09,0,lc};
Point(6)={0.2,0.11,0,lc};
Point(7)={0.09,0.2,0,lc};
Point(8)={0.11,0.2,0,lc};
Point(9)={0.09,0.09,0,lc2};
Point(10)={0.09,0.11,0,lc2};
Point(11)={0.11,0.11,0,lc2};
Point(12)={0.11,0.09,0,lc2};
Line(1) = {1,2};
Line(2) = {2,12};
Line(3) = {12,5};
Line(4) = {5,6};
Line(5) = {6,11};
Line(6) = {11,8};
Line(7) = {8,7};
Line(8) = {7,10};
Line(9) = {10,4};
Line(10) = {4,3};
Line(11) = {3,9};
Line(12) = {9,1};
Line Loop(13) = {2,3,4,5,6,7,8,9,10,11,12,1};
Plane Surface(14) = {13};
turns = 1/3;
zz = 0.01;
cc = 0.1;
Extrude Surface {news-1, {0,0,2*zz}, {0,0,1} , {cc,cc,0} , 0}
{ Layers {2}; Recombine; };
For j In {1:3}
Extrude Surface {news-1, {0,0,zz}, {0,0,1} , {cc,cc,0} , Pi/(40*(4-j))}
{ Layers {1}; Recombine; };
EndFor
For j In {1:turns*8}
Extrude Surface {news-1, {0,0,10*zz}, {0,0,1} , {cc,cc,0} , Pi/4}
{ Layers {10}; Recombine; };
EndFor
For j In {1:3}
Extrude Surface {news-1, {0,0,zz}, {0,0,1} , {cc,cc,0} , Pi/(40*j)}
{ Layers {1}; Recombine; };
EndFor
Extrude Surface {news-1, {0,0,2*zz}, {0,0,1} , {cc,cc,0} , 0}
{ Layers {2}; Recombine; };
p = newp;
Point(p+1000)={0.1, -0.02, 0,lc};
Point(p+1001)={0.1, -0.01, 0,lc};
Point(p+1002)={0.1, -0.02, 0.01,lc};
Point(p+1003)={0.1, -0.03, 0,lc};
Point(p+1004)={0.1, -0.02, -0.01,lc};
Circle(1635) = {p+1001,p+1000,p+1002};
Circle(1636) = {p+1002,p+1000,p+1003};
Circle(1637) = {p+1003,p+1000,p+1004};
Circle(1638) = {p+1004,p+1000,p+1001};
Line Loop(1639) = {1635,1636,1637,1638};
Plane Surface(1640) = {1639};
turns = 2;
zz = 0.01/4;
For j In {1:8*turns}
Extrude Surface {news-1, {0,0,10*zz}, {0,0,1} , {cc,cc,0} , Pi/4}
{ Layers {10}; Recombine; };
EndFor
Point(p+4000)={0.1, -0.02 + 0.24, 0, lc};
Point(p+4001)={0.1, -0.01 + 0.24, 0, lc};
Point(p+4002)={0.1, -0.02 + 0.24, 0.01, lc};
Point(p+4003)={0.1, -0.03 + 0.24, 0, lc};
Point(p+4004)={0.1, -0.02 + 0.24, -0.01, lc};
Circle(4635) = {p+4001,p+4000,p+4002};
Circle(4636) = {p+4002,p+4000,p+4003};
Circle(4637) = {p+4003,p+4000,p+4004};
Circle(4638) = {p+4004,p+4000,p+4001};
Line Loop(4639) = {4635,4636,4637,4638};
Plane Surface(4640) = {4639};
For j In {1:8*turns}
Extrude Surface {news-1, {0,0,10*zz}, {0,0,1} , {cc,cc,0} , Pi/4}
{ Layers {10}; Recombine; };
EndFor
// options
Mesh.PartitionByExtrusion=1; // forces partitioning by columns
Mesh.NbPartitions=20;
Mesh.ColorCarousel = 3;
Mesh.SurfaceFaces= 1;
Mesh.VolumeEdges = 0;
Mesh.VolumeFaces = 1;
Mesh.Color.Zero = {255,0,0};
Mesh.Color.One = {0,0,255};
Mesh.Color.Two = {0,255,0};
Mesh.Color.Three = {255,255,0};
Mesh.Color.Four = {255,0,255};
Mesh.Color.Six = {255,128,0};
Mesh.Color.Seven = {64,128,32};
lc = 0.1;
Point(1)={0.0,0.0,0.0,lc};
Point(2)={1.0,0.0,0.0,lc};
Line(1)={1,2};
Extrude Line { 1, {0,0.1,0.0} }
{ Layers { {10,10}, {10,20}, {0.2,1} };
Recombine;
};
Extrude Point { 1, {0,-0.2,-0.3} }
{ Layers { {10,10}, {11,21}, {0.2,1} }; };
Extrude Surface{ 5, {0,0.1,0.6}}
{ Layers { {10,10}, {110,210}, {0.2,1} }; Recombine; };
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment