]> SALOME platform Git repositories - tools/medcoupling.git/blob - doc/tutorial/medcoupling_2Dpolygon.rst
Salome HOME
Merge branch 'master' of https://codev-tuleap.cea.fr/plugins/git/salome/medcoupling
[tools/medcoupling.git] / doc / tutorial / medcoupling_2Dpolygon.rst
1
2 Meshing a 2D geometry with polygons
3 -----------------------------------
4
5 Objective
6 ~~~~~~~~~
7
8 In this example we propose to build a mesh on a 2D geometry. We first create triangles in the center of the geometry and then hexagons around them. The result will be saved in a med file, and then visualized with the SMESH module of Salome.
9
10 .. image:: images/Mesh_polygons2D.jpg
11
12
13 Beginning of implementation
14 ~~~~~~~~~~~~~~~~~~~~~~~~~~~
15
16 To implement this exercice we use the python language script and import the MEDCoupling and MEDLoader parts of the MED module::
17
18         from MEDCoupling import *
19         from MEDLoader import *
20
21         from math import *
22
23 Then we must instantiate a meshing object::
24
25         mesh=MEDCouplingUMesh.New()
26         mesh.setMeshDimension(2)
27         mesh.allocateCells(numberOfCells)
28         mesh.setName("MaFleur")
29
30
31 You must define 2 variables for space dimension and total number of nodes::
32
33         numberOfNodes = 25
34         numberOfCells = 12
35         spaceDimension = 2
36
37
38 You must define the coordinates of the nodes of the central hexagon::
39
40         X = [1.,0.5,-0.5,-1.,-0.5,0.5]
41         Y = [0.,sqrt(3.)/2.,sqrt(3.)/2.,0.,-sqrt(3.)/2.,-sqrt(3.)/2.]
42
43 Definition of nodes coordinates
44 ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
45
46 Define the coordinates of the origin and the nodes of the 2D central polygon::
47
48         coordinates = []
49         # origin
50         coordinates.append(0.)
51         coordinates.append(0.)
52
53         # central polygon
54         for i in range(6):
55           coordinates.append(X[i])
56           coordinates.append(Y[i])
57
58 Define the coordinates of the extra nodes of the hexagons around the central hexagon::
59
60         # second couron
61         xc = 1.5
62         yc = - sqrt(3.)/2.
63
64         for i in range(6):
65           ...
66
67         myCoords=DataArrayDouble.New()
68         myCoords.setValues(coordinates,numberOfNodes,2)
69         mesh.setCoords(myCoords)
70
71
72 Definition of triangles connectivity
73 ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
74
75 For each triangle of the mesh, you have to give its connectivity: the list of the nodes which belong to the triangles::
76
77         connectivity = []
78         for i in range(6):
79                 connectivity.append(0)
80                 connectivity.append(...)
81                 connectivity.append(...)
82         for i in range(6):
83                 mesh.insertNextCell(NORM_TRI3,3,connectivity[3*i:3*(i+1)])
84                 pass
85
86
87 Definition of hexagons connectivity
88 ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
89
90 For each hexagon of the mesh, you have to give its connectivity: the list of the nodes which belong to the hexagons::
91
92         connectivity = []
93         for i in range(6):
94                 start = i%6+1
95                 connectivity.append(...)
96         for i in range(6):
97                 mesh.insertNextCell(NORM_POLYGON,6,connectivity[6*i:6*(i+1)])
98                 pass
99         mesh.checkConsistencyLight()
100
101 Saving the mesh in a med file
102 ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
103
104 You have to create a med file with the MED driver::
105
106         medFileName = "MEDCoupling_Fleur.med"
107         MEDLoader.WriteUMesh(medFileName,mesh,True)
108
109 Visualize the mesh with SMESH module of Salome
110 ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
111
112 Launch Salome platform, then select SMESH module and import your MED file. First You can see the number of elements of your mesh. For that, select your mesh in the object browser, set select Mesh menu and "Advanced Mesh Info" option. Verify the number of triangles and the number of hexagons. To visualize your mesh: click right bottom on your mesh and select "Show" option.
113
114 Solution
115 ~~~~~~~~
116
117 :ref:`python_testMEDCoupling2D_solution`