1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
|
% Copyright (C) 2011-2016 Yves Renard.
%
% This file is a part of GetFEM++
%
% GetFEM++ is free software; you can redistribute it and/or modify it
% under the terms of the GNU Lesser General Public License as published
% by the Free Software Foundation; either version 3 of the License, or
% (at your option) any later version along with the GCC Runtime Library
% Exception either version 3.1 or (at your option) any later version.
% This program is distributed in the hope that it will be useful, but
% WITHOUT ANY WARRANTY; without even the implied warranty of MERCHANTABILITY
% or FITNESS FOR A PARTICULAR PURPOSE. See the GNU Lesser General Public
% License and GCC Runtime Library Exception for more details.
% You should have received a copy of the GNU Lesser General Public License
% along with this program; if not, write to the Free Software Foundation,
% Inc., 51 Franklin St, Fifth Floor, Boston, MA 02110-1301, USA.
%
% Mesh generation with the experimental meshing procedure of Getfem which
% uses simple primitives to describe the mesh geometry.
%
% This program is used to check that matlab-getfem is working. This is also
% a good example of use of GetFEM++.
%
gf_workspace('clear all'); clear all;
N = 2; % dimension of the mesh
K = 2; % degree of the mesh (for curved boundaries)
if (N == 1)
mo = gf_mesher_object('ball', [0], 2);
fixed_vertices = [0];
h = 0.5;
elseif (N == 2)
if (0)
mo = gf_mesher_object('ball', [0 4], 2);
fixed_vertices = [0; 4];
h = 0.5;
else
mo1 = gf_mesher_object('ball', [0 1.1], 1);
mo2 = gf_mesher_object('rectangle', [-0.2 1.085], [0.2, 1.115]);
mo3 = gf_mesher_object('ball', [0.4 1.1], 0.2);
mo4 = gf_mesher_object('ball', [-0.4 1.1], 0.2);
mo5 = gf_mesher_object('union', mo2, mo3, mo4);
mo = gf_mesher_object('set minus', mo1, mo5);
fixed_vertices = [];
h = 0.05;
end
elseif (N == 3)
if (0)
mo1 = gf_mesher_object('ball', [0 0 1], 2);
mo2 = gf_mesher_object('ball', [0 0 -1], 2);
mo3 = gf_mesher_object('intersect', mo1, mo2);
mo4 = gf_mesher_object('ball', [0 0 0], 1.3);
mo5 = gf_mesher_object('union', mo4, mo3);
mo6 = gf_mesher_object('ball', [-1 0 0], 1);
mo = gf_mesher_object('set minus', mo5, mo6);
fixed_vertices = []; h = 0.3;
elseif (1)
alpha = pi/5;
L = 20;
R = L * tan(alpha) * 0.7;
mo1 = gf_mesher_object('cone', [0 0 0], [0 0 1], L, alpha);
mo2 = gf_mesher_object('cylinder', [0 0 L], [0, 0, 1], L, R);
mo = gf_mesher_object('union', mo1, mo2);
fixed_vertices = []; h = 2.1;
else
mo = gf_mesher_object('ball', [0 0 20], 20);
fixed_vertices = []; h = 3;
end
elseif (N == 4)
mo = gf_mesher_object('ball', [0 0 0 4], 2);
fixed_vertices = [0; 0; 0; 4];
h = 1;
else
error('It is not very reasonable to build a mesh in dimension greater than 4 !');
end;
tl = gf_util('trace level');
% gf_util('trace level', 2); % no trace for mesh generation
gf_util('trace level', 4); % all traces for mesh generation
m = gf_mesh('generate', mo, h, K, fixed_vertices);
gf_util('trace level', tl);
if (N <= 2)
gf_plot_mesh(m);
elseif (N == 3)
mf=gf_mesh_fem(m, 1);
gf_mesh_fem_set(mf, 'classical fem', K);
VAL = gf_mesh_fem_get(mf, 'eval', {'x+y+z'});
gf_plot(mf, VAL, 'mesh', 'on', 'cvlst', gf_mesh_get(mf,'outer faces'), 'refine', 4);
axis on; camlight; gf_colormap('chouette');
end;
% gf_mesh_get(m, 'save', 'name.mesh');
|
