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

new examples

parent 8bf3b327
No related branches found
No related tags found
No related merge requests found
Show whitespace changes
Inline Side-by-side
Showing
with 359 additions and 135 deletions
demos/ex01.geo → demos/demo-01.geo
+ 10
− 43
View file @ 0d6d2c3d
......@@ -14,7 +14,12 @@ lc = 0.007;
/* README: This defines the point '1', located at x=0, y=0, z=0, with
an associated characteristic length 'lc'. The characteristic length
sets the target mesh size at the point. */
sets the target mesh size at the point.
The same convention is applied for all command: "If a number
defines a new entity, it is enclosed between parentheses. If a
number refers to a previously defined entity, it is enclosed
between braces." */
Point(1) = {0, 0, 0, lc};
......@@ -72,8 +77,9 @@ Color{
toggle the visibility of the particular view. A right mouse click
will give you access to the view's options:
- Reload: reload the file from which the view was loaded
- Remove: removes the view
- Duplicate: makes a formal copy of the view
- Duplicate: makes a copy of the view (without duplicating the data)
- Lighting: activates/deactivates lighting for the view
- Scale: gives access to the scale menu (range definition, iso-value
choice, ...)
......@@ -93,47 +99,8 @@ Color{
/* README: The 'msh' file format is the native output file format for
Gmsh. The file is divided in several sections (enclosed in $KEY and
$ENDKEY pairs). Two fields are important: $NOD/$ENDNOD defines the
nodes and $ELM/$ENDELM defines the elements.
The syntax is as follows:
$NOD
number-of-nodes
node-number x-coord y-coord z-coord
...
$ENDNOD
$ELM
number-of-elements
elm-number elm-type elm-region unused number-of-nodes node-numbers
...
$ENDELM
All the syntactic variables stand for integers except x-coord,
y-coord and z-coord which stand for floating point values. The
elm-type value defines the geometrical type for the element:
elm-type:
1 Line (2 nodes, 1 edge).
2 Triangle (3 nodes, 3 edges).
3 Quadrangle (4 nodes, 4 edges).
4 Tetrahedron (4 nodes, 6 edges, 4 facets).
5 Hexahedron (8 nodes, 12 edges, 6 facets).
6 Prism (6 nodes, 9 edges, 5 facets).
7 Pyramid (5 nodes, 8 edges, 5 facets).
15 Point (1 node).
The elm-region value is the number of the physical entity to which
the element belongs.
The 'unv' is the universal file format standard.
*/
/* README: The format of post-processing views is detailled in 'view01.pos' */
Gmsh. See the FORMAT file for an explanation of all the formats
used by Gmsh. */
......
demos/demo-02.geo 0 → 100644
+ 38
− 0
View file @ 0d6d2c3d
/*
Gmsh demo file (C) 2000 C. Geuzaine, J.-F. Remacle
Trivial 3D example
All important comments are marked with "README"
*/
// Let's include the first example
Include "demo-01.geo" ;
// We can extrude the surface to create a box of height h
h = 0.1 ;
Extrude Surface { 6, {0,0,h} };
// and specify a different characteristic length to some of the
// automatically created nodes
Characteristic Length{5,14} = lc * 5 ;
// And define a volume
/* README: In the same way one defines line loops to build surfaces,
one must define surface loops to build volumes. The following
volume is very simple, without holes (only one surface loop in the
volume) */
Surface Loop(10022) = {10021,10008,6,10012,10016,10020};
Volume(10023) = {10022};
// Save volumic elements of volume 10023 with the associate number '1'
Physical Volume (1) = {41} ;
demos/demo-03.geo 0 → 100644
+ 32
− 0
View file @ 0d6d2c3d
/*
Gmsh demo file (C) 2000 C. Geuzaine, J.-F. Remacle
Trivial 3D example, with extrusion
All important comments are marked with "README"
*/
// Let's include the first example, again
Include "demo-01.geo" ;
// In the same way as in demo-02.geo, we can extrude the surface to
// create a box of height h. But here we also extrude the 2D mesh into
// two layers (of 2 and 4 elements), with volume numbers 9000 and 9001
// and repective heights of 0.33*h and 0.67*h
h = 0.1 ;
Extrude Surface { 6, {0,0,h} } { Layers { {2,4}, {9000,9001}, {0.33,1} } ; };
// The extrusion can be combined with a rotation, and the extruded 3D
// mesh can be recombined into prisms (wedges). All rotations are
// specified by an axis direction ({0,1,0}), an axis point ({0,0,0})
// and a rotation angle (Pi/2)
Extrude Surface { 10021, {0,1,0} , {-0.1,0,0.1} , -Pi/2 } {
Recombine ; Layers { {7}, {9002}, {1} } ;
};
demos/ex06.geo → demos/demo-04.geo
+ 14
− 14
View file @ 0d6d2c3d
......@@ -103,27 +103,27 @@ Line Loop(122) = {32,-22,-27,21};Ruled Surface(123) = {122};
Line Loop(124) = {33,-23,-28,22};Ruled Surface(125) = {124};
Line Loop(126) = {34,-24,-29,23};Ruled Surface(127) = {126};
Surface Loop(128) = {93,-73,-55,95,-91}; Complex Volume(129) = {128}; /* int */
Surface Loop(130) = {107,-75,-97,95,57,121}; Complex Volume(131) = {130}; /* int b */
Surface Loop(132) = {105,-65,-97,-83,-93}; Complex Volume(133) = {132}; /* int h */
Surface Loop(134) = {99,-111,77,123,59,107}; Complex Volume(135) = {134}; /* shell b */
Surface Loop(136) = {99,-109,67,105,85}; Complex Volume(137) = {136}; /* shell h */
Surface Loop(138) = {113,79,-101,-111,-125,-61};Complex Volume(139) = {138}; /* ext b */
Surface Loop(140) = {115,-69,-101,-87,-109}; Complex Volume(141) = {140}; /* ext h */
Surface Loop(142) = {103,-117,-81,113,127,63}; Complex Volume(143) = {142}; /* inf b */
Surface Loop(144) = {89,-119,71,103,115}; Complex Volume(145) = {144}; /* inf h */
Surface Loop(128) = {93,-73,-55,95,-91}; Volume(129) = {128}; /* int */
Surface Loop(130) = {107,-75,-97,95,57,121}; Volume(131) = {130}; /* int b */
Surface Loop(132) = {105,-65,-97,-83,-93}; Volume(133) = {132}; /* int h */
Surface Loop(134) = {99,-111,77,123,59,107}; Volume(135) = {134}; /* shell b */
Surface Loop(136) = {99,-109,67,105,85}; Volume(137) = {136}; /* shell h */
Surface Loop(138) = {113,79,-101,-111,-125,-61}; Volume(139) = {138}; /* ext b */
Surface Loop(140) = {115,-69,-101,-87,-109}; Volume(141) = {140}; /* ext h */
Surface Loop(142) = {103,-117,-81,113,127,63}; Volume(143) = {142}; /* inf b */
Surface Loop(144) = {89,-119,71,103,115}; Volume(145) = {144}; /* inf h */
/* README: Tranfinite line commands explicitly specify the number of
points and their repartition. A minus sign in the argument list of
the transfinite command will produce the reversed mesh. */
Transfinite Line{35,21,22,23,24,38,17,18,19,20} = nbpt_phi Using Power 1.0;
Transfinite Line{35,21,22,23,24,38,17,18,19,20} = nbpt_phi ;
Transfinite Line{31,26,48,44,42} = nbpt_int Using Progression 0.95;
Transfinite Line{41,37,36,9,11,1,3,13,5,15,7} = nbpt_arc1 Using Power 1.0;
Transfinite Line{30,25,40,39,10,2,12,4,14,6,16,8} = nbpt_arc2 Using Power 1.0;
Transfinite Line{32,27,49,45,43} = nbpt_shell Using Power 1.0;
Transfinite Line{41,37,36,9,11,1,3,13,5,15,7} = nbpt_arc1 ;
Transfinite Line{30,25,40,39,10,2,12,4,14,6,16,8} = nbpt_arc2 ;
Transfinite Line{32,27,49,45,43} = nbpt_shell ;
Transfinite Line{33,28,46,50,52} = nbpt_far Using Progression 1.05 ;
Transfinite Line{34,29,51,47,53} = nbpt_inf Using Progression 1.2;
Transfinite Line{34,29,51,47,53} = nbpt_inf Using Progression 0.01;
/* README: *All* 2D and 3D transfinite entities are defined in respect
to points. The ordering of the points defines the ordering of the
......
demos/ex02.geo → demos/demo-05.geo
+ 63
− 0
View file @ 0d6d2c3d
......@@ -6,47 +6,35 @@
All important comments are marked with "README"
*/
// -clscale 0.2 ==> 169932 Nodes, 1005779 Elements in 1095.95 s
/* README: You can specify a global scaling factor for all
characteristic lengths in the geometry with the command line option
'-clscale'
With -clscale 1, this example produces a mesh of 2126 Nodes and
10221 Elements (in 3.61706 seconds on my DEC alpha 666MHz)
With -clscale 0.2 (i.e. all characteristic lengths a divided by 5),
this example produces a mesh of 169932 Nodes, and 1005779 (in
1095.95 s) */
lcar1 = .1;
lcar2 = .0005;
Point(1) = {0.5,0.5,0.5,lcar2};
Point(2) = {0.5,0.5,0,lcar1};
Point(3) = {0,0.5,0.5,lcar1};
Point(4) = {0,0,0.5,lcar1};
Point(5) = {0.5,0,0.5,lcar1};
Point(6) = {0.5,0,0,lcar1};
Point(7) = {0,0.5,0,lcar1};
Point(8) = {0,1,0,lcar1};
Point(9) = {1,1,0,lcar1};
Point(10) = {0,0,1,lcar1};
Point(11) = {0,1,1,lcar1};
Point(12) = {1,1,1,lcar1};
Point(13) = {1,0,1,lcar1};
Point(14) = {1,0,0,lcar1};
Line(1) = {8,9};
Line(2) = {9,12};
Line(3) = {12,11};
Line(4) = {11,8};
Line(5) = {9,14};
Line(6) = {14,13};
Line(7) = {13,12};
Line(8) = {11,10};
Line(9) = {10,13};
Line(10) = {10,4};
Line(11) = {4,5};
Line(12) = {5,6};
Line(13) = {6,2};
Line(14) = {2,1};
Line(15) = {1,3};
Line(16) = {3,7};
Line(17) = {7,2};
Line(18) = {3,4};
Line(19) = {5,1};
Line(20) = {7,8};
Line(21) = {6,14};
Point(1) = {0.5,0.5,0.5,lcar2}; Point(2) = {0.5,0.5,0,lcar1};
Point(3) = {0,0.5,0.5,lcar1}; Point(4) = {0,0,0.5,lcar1};
Point(5) = {0.5,0,0.5,lcar1}; Point(6) = {0.5,0,0,lcar1};
Point(7) = {0,0.5,0,lcar1}; Point(8) = {0,1,0,lcar1};
Point(9) = {1,1,0,lcar1}; Point(10) = {0,0,1,lcar1};
Point(11) = {0,1,1,lcar1}; Point(12) = {1,1,1,lcar1};
Point(13) = {1,0,1,lcar1}; Point(14) = {1,0,0,lcar1};
Line(1) = {8,9}; Line(2) = {9,12}; Line(3) = {12,11};
Line(4) = {11,8}; Line(5) = {9,14}; Line(6) = {14,13};
Line(7) = {13,12}; Line(8) = {11,10}; Line(9) = {10,13};
Line(10) = {10,4}; Line(11) = {4,5}; Line(12) = {5,6};
Line(13) = {6,2}; Line(14) = {2,1}; Line(15) = {1,3};
Line(16) = {3,7}; Line(17) = {7,2}; Line(18) = {3,4};
Line(19) = {5,1}; Line(20) = {7,8}; Line(21) = {6,14};
Line Loop(22) = {11,19,15,18}; Plane Surface(23) = {22};
Line Loop(24) = {16,17,14,15}; Plane Surface(25) = {24};
......@@ -58,12 +46,8 @@ Line Loop(34) = {7,3,8,9}; Plane Surface(35) = {34};
Line Loop(36) = {10,-18,16,20,-4,8}; Plane Surface(37) = {36};
Line Loop(38) = {-14,-13,-12,19}; Plane Surface(39) = {38};
/* README: In the same way one defines line loops to build surfaces,
one must define surface loops to build volumes. The following
volume is very simple, without holes (only one surface loop) */
Surface Loop(40) = {35,31,29,37,33,23,39,25,27};
Complex Volume(41) = {40};
Volume(41) = {40};
// Save volumic elements in region 1000
......
demos/ex03.geo → demos/demo-06.geo
+ 4
− 4
View file @ 0d6d2c3d
/*
Gmsh demo file (C) 2000 C. Geuzaine, J.-F. Remacle
Simple 3D example
Simple 3D example, with customization commands
All important comments are marked with "README"
*/
......@@ -130,13 +130,13 @@ Ruled Surface(76) = {75};
// Definition of volumes
Surface Loop(77) = {60,-58,54,-44,-56};
Complex Volume(78) = {77};
Volume(78) = {77};
Surface Loop(79) = {64,-68,-74,-62,-48,50,66,76};
Complex Volume(80) = {79};
Volume(80) = {79};
Surface Loop(81) = {64,-68,-74,-62,70,-46,-72,58,-60,42,-40,-52,76};
Complex Volume(82) = {81};
Volume(82) = {81};
// Some definitions for saving the mesh
......
demos/ex04.geo → demos/demo-07.geo
+ 1
− 1
View file @ 0d6d2c3d
/*
Gmsh demo file (C) 2000 C. Geuzaine, J.-F. Remacle
2D example with hole
Trivial 2D example, with hole
All important comments are marked with "README"
*/
......
demos/ex05.geo → demos/demo-08.geo
+ 24
− 28
View file @ 0d6d2c3d
/*
Gmsh demo file (C) 2000 C. Geuzaine, J.-F. Remacle
More complex 2D example
More complex 2D example with holes
All important comments are marked with "README"
*/
......@@ -112,22 +112,14 @@ R2 = r2r ;
*ALL* variables in Gmsh are global. Including a file is similar to
paste its content where the include command is located. */
i = 1 ; th = th0r + (i - 1) * dthr ;
Include "ex05.i1" ;
i = 2 ; th = th0r + (i - 1) * dthr ;
Include "ex05.i1" ;
i = 3 ; th = th0r + (i - 1) * dthr ;
Include "ex05.i1" ;
i = 4 ; th = th0r + (i - 1) * dthr ;
Include "ex05.i1" ;
i = 5 ; th = th0r + (i - 1) * dthr ;
Include "ex05.i1" ;
i = 6 ; th = th0r + (i - 1) * dthr ;
Include "ex05.i1" ;
i = 7 ; th = th0r + (i - 1) * dthr ;
Include "ex05.i1" ;
i = 8 ; th = th0r + (i - 1) * dthr ;
Include "ex05.i1" ;
i = 1 ; th = th0r + (i - 1) * dthr ; Include "demo-08.i1" ;
i = 2 ; th = th0r + (i - 1) * dthr ; Include "demo-08.i1" ;
i = 3 ; th = th0r + (i - 1) * dthr ; Include "demo-08.i1" ;
i = 4 ; th = th0r + (i - 1) * dthr ; Include "demo-08.i1" ;
i = 5 ; th = th0r + (i - 1) * dthr ; Include "demo-08.i1" ;
i = 6 ; th = th0r + (i - 1) * dthr ; Include "demo-08.i1" ;
i = 7 ; th = th0r + (i - 1) * dthr ; Include "demo-08.i1" ;
i = 8 ; th = th0r + (i - 1) * dthr ; Include "demo-08.i1" ;
// 9 stator poles
......@@ -142,31 +134,31 @@ E2 = e2s ;
H1 = h1s ;
i = 1 ; th = th0s + (i - 1) * dths ;
Include "ex05.i2" ;
Include "demo-08.i2" ;
PP2 = p1 ; PP3 = p9 ;
i = 2 ; th = th0s + (i - 1) * dths ;
Include "ex05.i2" ;
Include "demo-08.i2" ;
PP4 = p1 ; PP5 = p9 ;
i = 3 ; th = th0s + (i - 1) * dths ;
Include "ex05.i2" ;
Include "demo-08.i2" ;
PP6 = p1 ; PP7 = p9 ;
i = 4 ; th = th0s + (i - 1) * dths ;
Include "ex05.i2" ;
Include "demo-08.i2" ;
PP8 = p1 ; PP9 = p9 ;
i = 5 ; th = th0s + (i - 1) * dths ;
Include "ex05.i2" ;
Include "demo-08.i2" ;
PP10 = p1 ; PP11 = p9 ;
i = 6 ; th = th0s + (i - 1) * dths ;
Include "ex05.i2" ;
Include "demo-08.i2" ;
PP12 = p1 ; PP13 = p9 ;
i = 7 ; th = th0s + (i - 1) * dths ;
Include "ex05.i2" ;
Include "demo-08.i2" ;
PP14 = p1 ; PP15 = p9 ;
i = 8 ; th = th0s + (i - 1) * dths ;
Include "ex05.i2" ;
Include "demo-08.i2" ;
PP16 = p1 ; PP17 = p9 ;
i = 9 ; th = th0s + (i - 1) * dths ;
Include "ex05.i2" ;
Include "demo-08.i2" ;
PP18 = p1 ; PP19 = p9 ;
lin1 = newreg ; Line(lin1) = {PP1 , PP2 } ;
......@@ -184,10 +176,14 @@ lin1 = newreg ; Line(lin1) = {PP19, PPB } ;
Line Loop(145) = {8,-2,6,7};
Plane Surface(146) = {145,68,61,54,47,40,33,26,19};
Line Loop(147) = {-7,9,133,-74,134,-81,135,-88,136,-95,137,-102,138,-109,139,-116,140,-123,141,-130,142,10};
Line Loop(147) = {-7,9,133,-74,134,-81,135,-88,136,-95,137,-102,
138,-109,139,-116,140,-123,141,-130,142,10};
Plane Surface(148) = {147};
Line Loop(149) = {70,71,72,73,134,77,78,79,80,135,84,85,86,87,136,91,92,93,94,137,98,99,100,101,138,105,106,107,108,139,112,113,114,115,140,119,120,121,122,141,126,127,128,129,142,-13,-12,-11,133};
Line Loop(149) = {70,71,72,73,134,77,78,79,80,135,84,85,86,87,136,
91,92,93,94,137,98,99,100,101,138,105,106,107,108,
139,112,113,114,115,140,119,120,121,122,141,126,127,
128,129,142,-13,-12,-11,133};
Plane Surface(150) = {149};
/* README: One should define physical regions to specify what to
......
demos/ex05.i1 → demos/demo-08.i1
+ 0
− 0
View file @ 0d6d2c3d
File moved
demos/ex05.i2 → demos/demo-08.i2
+ 0
− 0
View file @ 0d6d2c3d
File moved
demos/ex07.geo → demos/demo-09.geo
+ 6
− 3
View file @ 0d6d2c3d
/*
Gmsh demo file (C) 2000 C. Geuzaine, J.-F. Remacle
3D example with (too) large characteristic lengths
3D example with large characteristic lengths
All important comments are marked with "README"
*/
......@@ -429,7 +429,7 @@ Complex Volume(608) = {607};
Surface Loop(6001) = {213,-239,-182,180,186,237,231,-229,227,-233,-188,-201,-199,-219,132,138,136,134,112,156,-114,162,-160,158,-140,101,-146,-144,-142,154,148,-105,-152,-150};
Complex Volume(6002) = {6001};
Delete { Line(169); }
Delete { Line{169} ; }
Line(6003) = {83,87};
Line(6004) = {87,85};
......@@ -480,6 +480,9 @@ Physical Volume (R1) = 602;
Physical Volume (R2) = 604;
Physical Volume (SOU) = 6032;
Physical Surface (CLD0) = {231,229,223,237,227, 199,201,219,209,207,203,205, 180,182,188,176,178,174,172 , 6026,6024,6028,6030};
Physical Surface (CLD0) = {231,229,223,237,227,
199,201,219,209,207,203,205,
180,182,188,176,178,174,172,
6026,6024,6028,6030};
Physical Surface (CLTEM) = {211,6022};
demos/ex09.geo → demos/demo-10.geo
+ 125
− 0
View file @ 0d6d2c3d
......@@ -35,62 +35,58 @@ Circle(6) = {3,2,5};
r = 2*3.14159/5;
Point(10) = { (l1 + r2) * Cos(r/2) , (l1 + r2) * Sin(r/2), 0.0, lc};
Rotate({0.0,0.0,1.0},{0.0,0.0,0.0},r) {
// Remember, all rotations are specified by the axis direction
// ({0,0,1}), an axis point ({0,0,0}) and a rotation angle (r)
Rotate {{0.0,0.0,1.0},{0.0,0.0,0.0},r} {
Duplicata {
Line(1);
Line(2);
Line(3);
Line(4);
Line(5);
Line(6);
Point(10);
Line{1};
Line{2};
Line{3};
Line{4};
Line{5};
Line{6};
Point{10};
}
}
Rotate({0.0,0.0,1.0},{0.0,0.0,0.0},2*r) {
Rotate{{0.0,0.0,1.0},{0.0,0.0,0.0},2*r} {
Duplicata {
Line(1);
Line(2);
Line(3);
Line(4);
Line(5);
Line(6);
Point(10);
Line{1};
Line{2};
Line{3};
Line{4};
Line{5};
Line{6};
Point{10};
}
}
Rotate({0.0,0.0,1.0},{0.0,0.0,0.0},3*r) {
Rotate{{0.0,0.0,1.0},{0.0,0.0,0.0},3*r} {
Duplicata {
Line(1);
Line(2);
Line(3);
Line(4);
Line(5);
Line(6);
Point(10);
Line{1};
Line{2};
Line{3};
Line{4};
Line{5};
Line{6};
Point{10};
}
}
Rotate({0.0,0.0,1.0},{0.0,0.0,0.0},4*r) {
Rotate{{0.0,0.0,1.0},{0.0,0.0,0.0},4*r} {
Duplicata {
Line(1);
Line(2);
Line(3);
Line(4);
Line(5);
Line(6);
Point(10);
Line{1};
Line{2};
Line{3};
Line{4};
Line{5};
Line{6};
Point{10};
}
}
Coherence;
Circle(31) = {16,10,9};
Circle(32) = {31,25,14};
Circle(33) = {46,40,29};
Circle(34) = {61,55,44};
Circle(35) = {8,70,59};
Point(newp) = {rint,0,0,lc};
Point(newp) = {rext,0,0,lc};
Point(newp) = {-rint,0,0,lc};
......@@ -99,32 +95,31 @@ Point(newp) = {0,rint,0,lc};
Point(newp) = {0,rext,0,lc};
Point(newp) = {0,-rint,0,lc};
Point(newp) = {0,-rext,0,lc};
Circle(36) = {75,1,73};
Circle(37) = {73,1,77};
Circle(38) = {77,1,71};
Circle(39) = {71,1,75};
Circle(40) = {76,1,74};
Circle(41) = {74,1,78};
Circle(42) = {78,1,72};
Circle(43) = {72,1,76};
Line Loop(44) = {-8,-7,11,12,10,9,31,-2,-1,5,6,4,3,35,-26,-25,29,30,28,27,34,-20,-19,23,24,22,21,33,-14,-13,17,18,16,15,32};
Line Loop(45) = {40,41,42,43};
Plane Surface(46) = {44,45};
Line Loop(47) = {36,37,38,39};
Plane Surface(48) = {45,47};
Extrude(46, {0,0,.2});
Coherence;
Extrude(48, {0,0,.2});
Coherence;
Extrude(238, {0,0,.2});
Circle(31) = {8,118,97};
Circle(32) = {20,10,9};
Circle(33) = {47,37,16};
Circle(34) = {74,64,43};
Circle(35) = {101,91,70};
Circle(36) = {119,1,123};
Circle(37) = {123,1,121};
Circle(38) = {121,1,125};
Circle(39) = {125,1,119};
Circle(40) = {124,1,122};
Circle(41) = {122,1,126};
Circle(42) = {126,1,120};
Circle(43) = {120,1,124};
Line Loop(44) = {36,37,38,39};
Plane Surface(45) = {44};
Line Loop(46) = {43,40,41,42};
Plane Surface(47) = {46,44};
Line Loop(48) = {-26,-25,29,30,28,27,35,-20,-19,23,24,22,21,34,
-14,-13,17,18,16,15,33,-8,-7,11,12,10,9,32,-2,
-1,5,6,4,3,31};
Plane Surface(49) = {48,46};
Extrude Surface {45, {0,0,0.2}};
Coherence;
Surface Loop(274) = {235,-48,46,88,-189,-207,91,-94,-97,-100,-103,-106,109,112,-115,-118,-121,-124,-127,130,133,-136,-139,-142,-145,-148,151,154,-157,-160,-163,-166,-169,172,175,-178,-181,-184,-187,247,250,253,255,273,264,229,232,267,270,272,237};
Complex Volume(275) = {274};
Physical Volume(1000) = 275;
Surface Loop(72) = {45,58,62,66,70,71};
Complex Volume(73) = {72};
demos/ex12.geo → demos/demo-12.geo
+ 1
− 1
View file @ 0d6d2c3d
......@@ -6,7 +6,7 @@
All important comments are marked with "README"
*/
Include "ex12.a" ;
Include "demo-12.i1" ;
/* Post-Processing meshes */
......
demos/ex12.a → demos/demo-12.i1
+ 8
− 8
View file @ 0d6d2c3d
......@@ -29,36 +29,36 @@ z8 = 150 ; dy8 = 0. ;
p0 = 2. ; p1 = 30. ;
Include "ex12.b" ;
Include "demo-12.i2" ;
/* ------------ */
/* Fil 1a, b, c */
/* ------------ */
x = x1a ; y = y1a ; i_p = 100 ; i_l = 10 ;
Include "ex12.c" ;
Include "demo-12.i3" ;
x = x1b ; y = y1b ; i_p = 200 ; i_l = 20 ;
Include "ex12.c" ;
Include "demo-12.i3" ;
x = x1c ; y = y1c ; i_p = 300 ; i_l = 30 ;
Include "ex12.c" ;
Include "demo-12.i3" ;
/* ------------ */
/* Fil 2a, b, c */
/* ------------ */
x = x2a ; y = y2a ; i_p = 400 ; i_l = 40 ;
Include "ex12.c" ;
Include "demo-12.i3" ;
x = x2b ; y = y2b ; i_p = 500 ; i_l = 50 ;
Include "ex12.c" ;
Include "demo-12.i3" ;
x = x2c ; y = y2c ; i_p = 600 ; i_l = 60 ;
Include "ex12.c" ;
Include "demo-12.i3" ;
/* ------ */
/* Fil g */
/* ------ */
x = xg ; y = yg ; i_p = 700 ; i_l = 70 ;
Include "ex12.c" ;
Include "demo-12.i3" ;
/* --------------------------------------------- */
......
demos/ex12.b → demos/demo-12.i2
+ 11
− 11
View file @ 0d6d2c3d
......@@ -41,33 +41,33 @@ pt0 = 1. ; pt1 = 1. ; pt2 = 1. ; pt3 = 1. ; pt4 = 1. ; pt5 = 1. ;
pt6 = 1. ; pt7 = 1. ; pt8 = 1. ; pt9 = 1. ;
ci = c0 ; hi = h0 ; pti = pt0 ; i_p = 0 ;
Include "ex12.d" ;
Include "demo-12.i4" ;
ci = c1 ; hi = h1 ; pti = pt1 ; i_p = 10 ;
Include "ex12.d" ;
Include "ex12.e" ;
Include "demo-12.i4" ;
Include "demo-12.i5" ;
ci = c2 ; hi = h2 ; pti = pt2 ; i_p = 20 ;
Include "ex12.d" ;
Include "ex12.e" ;
Include "demo-12.i4" ;
Include "demo-12.i5" ;
ci = c3 ; hi = h3 ; pti = pt3 ; i_p = 30 ;
Include "ex12.d" ;
Include "demo-12.i4" ;
ci = c4 ; hi = h4 ; pti = pt4 ; i_p = 40 ;
Include "ex12.d" ;
Include "demo-12.i4" ;
ci = c5 ; hi = h5 ; pti = pt5 ; i_p = 50 ;
Include "ex12.d" ;
Include "demo-12.i4" ;
ci = c6 ; hi = h6 ; pti = pt6 ; i_p = 60 ;
Include "ex12.d" ;
Include "demo-12.i4" ;
ci = c7 ; hi = h7 ; pti = pt7 ; i_p = 70 ;
Include "ex12.d" ;
Include "demo-12.i4" ;
ci = c8 ; hi = h8 ; pti = pt8 ; i_p = 80 ;
Include "ex12.d" ;
Include "demo-12.i4" ;
ci = c9 ; hi = h9 ; pti = pt9 ; i_p = 90 ;
Point(91) = { ci , hi , 0. , pti } ;
......
demos/ex12.c → demos/demo-12.i3
+ 0
− 0
View file @ 0d6d2c3d
File moved
demos/ex12.d → demos/demo-12.i4
+ 0
− 0
View file @ 0d6d2c3d
File moved
demos/ex12.e → demos/demo-12.i5
+ 0
− 0
View file @ 0d6d2c3d
File moved
demos/ex13.geo → demos/demo-13.geo
+ 22
− 22
View file @ 0d6d2c3d
......@@ -78,28 +78,28 @@ l20 = 1404.e-3/2-ll/2 ; r20 = 12.e-3/2 ; d20 = 1241.7e-3 ; e20 = em+dc+d20*Sin(t
l21 = 1648.e-3/2-ll/2 ; r21 = 12.e-3/2 ; d21 = 1467.7e-3 ; e21 = em+dc+d21*Sin(t) ; s21 = f3*r21 ;
l22 = 1934.e-3/2-ll/2 ; r22 = 12.e-3/2 ; d22 = 1734.3e-3 ; e22 = em+dc+d22*Sin(t) ; s22 = f3*r22 ;
dx = d1 ; rx = r1 ; sx = s1 ; lx = l1 ; e = e1 ; x = 100; Include "ex13.bar" ;
dx = d2 ; rx = r2 ; sx = s2 ; lx = l2 ; e =-e2 ; x = 200; Include "ex13.bar" ;
dx = d3 ; rx = r3 ; sx = s3 ; lx = l3 ; e = e3 ; x = 300; Include "ex13.bar" ;
dx = d4 ; rx = r4 ; sx = s4 ; lx = l4 ; e =-e4 ; x = 400; Include "ex13.bar" ;
dx = d5 ; rx = r5 ; sx = s5 ; lx = l5 ; e = e5 ; x = 500; Include "ex13.bar" ;
dx = d6 ; rx = r6 ; sx = s6 ; lx = l6 ; e =-e6 ; x = 600; Include "ex13.bar" ;
dx = d7 ; rx = r7 ; sx = s7 ; lx = l7 ; e = e7 ; x = 700; Include "ex13.bar" ;
dx = d8 ; rx = r8 ; sx = s8 ; lx = l8 ; e =-e8 ; x = 800; Include "ex13.bar" ;
dx = d9 ; rx = r9 ; sx = s9 ; lx = l9 ; e = e9 ; x = 900; Include "ex13.bar" ;
dx = d10; rx = r10; sx = s10; lx = l10; e =-e10; x =1000; Include "ex13.bar" ;
dx = d11; rx = r11; sx = s11; lx = l11; e = e11; x =1100; Include "ex13.bar" ;
dx = d12; rx = r12; sx = s12; lx = l12; e =-e12; x =1200; Include "ex13.bar" ;
dx = d13; rx = r13; sx = s13; lx = l13; e = e13; x =1300; Include "ex13.bar" ;
dx = d14; rx = r14; sx = s14; lx = l14; e =-e14; x =1400; Include "ex13.bar" ;
dx = d15; rx = r15; sx = s15; lx = l15; e = e15; x =1500; Include "ex13.bar" ;
dx = d16; rx = r16; sx = s16; lx = l16; e =-e16; x =1600; Include "ex13.bar" ;
dx = d17; rx = r17; sx = s17; lx = l17; e = e17; x =1700; Include "ex13.bar" ;
dx = d18; rx = r18; sx = s18; lx = l18; e =-e18; x =1800; Include "ex13.bar" ;
dx = d19; rx = r19; sx = s19; lx = l19; e = e19; x =1900; Include "ex13.bar" ;
dx = d20; rx = r20; sx = s20; lx = l20; e =-e20; x =2000; Include "ex13.bar" ;
dx = d21; rx = r21; sx = s21; lx = l21; e = e21; x =2100; Include "ex13.bar" ;
dx = d22; rx = r22; sx = s22; lx = l22; e =-e22; x =2200; Include "ex13.bar" ;
dx = d1 ; rx = r1 ; sx = s1 ; lx = l1 ; e = e1 ; x = 100; Include "demo-13.i1" ;
dx = d2 ; rx = r2 ; sx = s2 ; lx = l2 ; e =-e2 ; x = 200; Include "demo-13.i1" ;
dx = d3 ; rx = r3 ; sx = s3 ; lx = l3 ; e = e3 ; x = 300; Include "demo-13.i1" ;
dx = d4 ; rx = r4 ; sx = s4 ; lx = l4 ; e =-e4 ; x = 400; Include "demo-13.i1" ;
dx = d5 ; rx = r5 ; sx = s5 ; lx = l5 ; e = e5 ; x = 500; Include "demo-13.i1" ;
dx = d6 ; rx = r6 ; sx = s6 ; lx = l6 ; e =-e6 ; x = 600; Include "demo-13.i1" ;
dx = d7 ; rx = r7 ; sx = s7 ; lx = l7 ; e = e7 ; x = 700; Include "demo-13.i1" ;
dx = d8 ; rx = r8 ; sx = s8 ; lx = l8 ; e =-e8 ; x = 800; Include "demo-13.i1" ;
dx = d9 ; rx = r9 ; sx = s9 ; lx = l9 ; e = e9 ; x = 900; Include "demo-13.i1" ;
dx = d10; rx = r10; sx = s10; lx = l10; e =-e10; x =1000; Include "demo-13.i1" ;
dx = d11; rx = r11; sx = s11; lx = l11; e = e11; x =1100; Include "demo-13.i1" ;
dx = d12; rx = r12; sx = s12; lx = l12; e =-e12; x =1200; Include "demo-13.i1" ;
dx = d13; rx = r13; sx = s13; lx = l13; e = e13; x =1300; Include "demo-13.i1" ;
dx = d14; rx = r14; sx = s14; lx = l14; e =-e14; x =1400; Include "demo-13.i1" ;
dx = d15; rx = r15; sx = s15; lx = l15; e = e15; x =1500; Include "demo-13.i1" ;
dx = d16; rx = r16; sx = s16; lx = l16; e =-e16; x =1600; Include "demo-13.i1" ;
dx = d17; rx = r17; sx = s17; lx = l17; e = e17; x =1700; Include "demo-13.i1" ;
dx = d18; rx = r18; sx = s18; lx = l18; e =-e18; x =1800; Include "demo-13.i1" ;
dx = d19; rx = r19; sx = s19; lx = l19; e = e19; x =1900; Include "demo-13.i1" ;
dx = d20; rx = r20; sx = s20; lx = l20; e =-e20; x =2000; Include "demo-13.i1" ;
dx = d21; rx = r21; sx = s21; lx = l21; e = e21; x =2100; Include "demo-13.i1" ;
dx = d22; rx = r22; sx = s22; lx = l22; e =-e22; x =2200; Include "demo-13.i1" ;
/*
......
demos/ex13.bar → demos/demo-13.i1
+ 0
− 0
View file @ 0d6d2c3d
File moved
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment