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

*** empty log message ***

parent df449c5a
Branches
Tags
No related merge requests found
Mesh.CharacteristicLengthFactor=2;
//begin "my example"
lowPr = 0.1;
Point(1) = {0,0,0,lowPr};
Point(2) = {1,0,0,lowPr};
Point(3) = {1,1,0,lowPr};
Point(4) = {0,1,0,lowPr};
Point(5) = {0.5,0.5,0,0};
// points to see where the DistMax and DistMin are
Point(6) = {0.05+0.5,0.5,0,lowPr};
Point(7) = {0.2+0.5,0.5,0,lowPr};
Point(8) = {0.3+0.5,0.5,0,lowPr};
Point(9) = {0.5+0.5,0.5,0,lowPr};
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(6) = {5};
Field[1] = Attractor;
Field[1].NodesList = {5};
Field[2] = Threshold;
Field[2].IField = 1;
Field[2].LcMin = 0.001;
Field[2].LcMax = 0.005;
Field[2].DistMin = 0.05;
Field[2].DistMax = 0.2;
Field[2].StopAtDistMax = 1; // don't impose element size outside DistMax
//Field[2].Sigmoid = 1; // use sigmoid instead of linear interpol
Field[3] = Attractor;
Field[3].NodesList = {5};
Field[4] = Threshold;
Field[4].IField = 3;
Field[4].LcMin = 0.005;
Field[4].LcMax = 0.1;
Field[4].DistMin = 0.3;
Field[4].DistMax = 0.5;
Field[5] = Min;
Field[5].FieldsList = {2,4};
Background Field = 5;
// Options
Mesh.RecombineAlgo = 2 ;
Mesh.MshFileVersion = 1;
// Variables
ch_length = 0.4 ;
// Noeuds
p1 = newp; Point(p1) = {0, 0, 0, ch_length} ;
p2 = newp; Point(p2) = {2, 0, 0, ch_length} ;
p3 = newp; Point(p3) = {2, 5, 0, ch_length} ;
p4 = newp; Point(p4) = {0, 5, 0, ch_length} ;
p5 = newp; Point(p5) = {0.5, 2, 0, ch_length} ;
p6 = newp; Point(p6) = {1, 1.6, 0, ch_length} ;
p7 = newp; Point(p7) = {1.2, 2.6, 0, ch_length} ;
// Lignes
l1 = newl; Line(l1) = {p1, p2} ;
l2 = newl; Line(l2) = {p2, p3} ;
l3 = newl; Line(l3) = {p3, p4} ;
l4 = newl; Line(l4) = {p4, p1} ;
l5 = newl; Line(l5) = {p6, p7} ;
// Surface Dalle
ll1 = newll; Line Loop (ll1) = {l1, l2, l3, l4} ;
s1 = news; Plane Surface (s1) = {ll1} ;
// Affinage des lignes
/*
Transfinite Line {l1} = 1 + 2 / ch_length ;
Transfinite Line {l2} = 1 + 5 / ch_length ;
Transfinite Line {l3} = 1 + 2 / ch_length ;
Transfinite Line {l4} = 1 + 5 / ch_length ;
Transfinite Line {l5} = 1 + 1.3 / ch_length ;
*/
Point{p5} In Surface{s1} ;
Line{l5} In Surface{s1} ;
// Transformation des triangles en quadrangles
Recombine Surface {s1} ;
$MeshFormat
2 0 8
$EndMeshFormat
$Nodes
4
1 0 0 0
2 1 0 0
3 1 1 0
4 0 1 0
$EndNodes
$Elements
1
1 3 3 0 1 0 1 2 3 4
$EndElements
$NodeData
1
"quad"
1
0.0
3
0
1
4
1 1
2 2
3 3
4 4
$EndNodeData
$MeshFormat
2 0 8
$EndMeshFormat
$Nodes
27
1 0 0 0
2 1 0 0
3 1 1 0
4 0 1 0
5 1 1 1
6 1 0 1
7 0 0 1
8 0 1 1
9 0.4999999999986718 1 0
10 1 0.5000000000013305 0
11 0.5000000000013305 0 0
12 0 0.4999999999986718 0
13 1 0.5000000000013305 1
14 0.5000000000013305 0 1
15 0 0.4999999999986718 1
16 0.4999999999986718 1 1
17 1 1 0.4999999999986718
18 1 0 0.4999999999986718
19 0 0 0.4999999999986718
20 0 1 0.4999999999986718
21 0.5 0.5 0
22 1 0.5 0.5
23 0.5 0 0.5
24 0 0.5 0.5
25 0.5 1 0.5
26 0.5 0.5 1
27 0.5 0.5 0.5
$EndNodes
$Elements
38
1 15 3 0 1 0 1
2 15 3 0 2 0 2
3 15 3 0 3 0 3
4 15 3 0 4 0 4
5 15 3 0 5 0 5
6 15 3 0 6 0 6
7 15 3 0 10 0 7
8 15 3 0 14 0 8
9 8 3 0 1 0 4 3 9
10 8 3 0 2 0 3 2 10
11 8 3 0 3 0 2 1 11
12 8 3 0 4 0 1 4 12
13 8 3 0 8 0 5 6 13
14 8 3 0 9 0 6 7 14
15 8 3 0 10 0 7 8 15
16 8 3 0 11 0 8 5 16
17 8 3 0 13 0 3 5 17
18 8 3 0 14 0 2 6 18
19 8 3 0 18 0 1 7 19
20 8 3 0 22 0 4 8 20
21 9 3 0 6 0 3 2 4 10 21 9
22 9 3 0 6 0 2 1 4 11 12 21
23 9 3 0 15 0 2 6 5 18 13 22
24 9 3 0 15 0 3 2 5 10 22 17
25 9 3 0 19 0 6 2 7 18 23 14
26 9 3 0 19 0 2 1 7 11 19 23
27 9 3 0 23 0 1 4 7 12 24 19
28 9 3 0 23 0 4 8 7 20 15 24
29 9 3 0 27 0 4 3 5 9 17 25
30 9 3 0 27 0 8 4 5 20 25 16
31 9 3 0 28 0 7 5 6 26 13 14
32 9 3 0 28 0 7 8 5 15 16 26
33 11 3 0 1 0 7 4 8 2 24 20 15 23 27 21
34 11 3 0 1 0 5 4 3 2 25 9 17 22 10 21
35 11 3 0 1 0 6 5 2 7 13 22 18 14 23 26
36 11 3 0 1 0 7 4 2 1 24 21 23 19 11 12
37 11 3 0 1 0 5 4 2 8 25 21 22 16 27 20
38 11 3 0 1 0 5 8 2 7 16 27 22 26 23 15
$EndElements
$NodeData
1
"2nd order view"
1
0.0
3
0
1
27
1 2
2 2
3 2
4 2
5 2
6 2
7 2
8 2
9 1
10 1
11 1
12 1
13 1
14 1
15 1
16 1
17 1
18 1
19 1
20 1
21 1
22 1
23 1
24 1
25 1
26 1
27 1
$EndNodeData
$MeshFormat
2.0 0 8
$EndMeshFormat
$Nodes
4
1 0 0 0
2 0 1 0
3 1 0 0
4 1 1 0
$EndNodes
$Elements
1
1 3 0 1 3 4 2
$EndElements
$NodeData
1
"A scalar view"
1
10.0
3
0
1
4
1 0.0
2 0.0
3 1.0
4 1.0
$EndNodeData
$MeshFormat
2.0 0 8
$EndMeshFormat
$Nodes
6
1 0 0 0
2 0 1.2 0
3 1.2 0 0
4 1.2 1.2 0
5 0.6 0 0
6 0.6 1.2 0
$EndNodes
$Elements
2
1 3 0 1 2 6 5
2 3 0 3 4 6 5
$EndElements
$NodeData
1
"A scalar view"
1
11.0
3
1
1
6
1 0.0
2 0.0
3 2.0
4 2.0
5 2.0
6 2.0
$EndNodeData
$NodeData
1
"A scalar view"
1
12.0
3
2
1
6
1 0.5
2 0.2
3 2.2
4 1.2
5 0.7
6 2.0
$EndNodeData
Source diff could not be displayed: it is too large. Options to address this: view the blob.
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment