]> SALOME platform Git repositories - modules/smesh.git/blob - doc/salome/gui/SMESH/input/tui_quality_controls.doc
Salome HOME
remove old files
[modules/smesh.git] / doc / salome / gui / SMESH / input / tui_quality_controls.doc
1 /*!
2
3 \page tui_quality_controls_page Quality Controls
4
5 <br>
6 \anchor tui_free_borders
7 <h2>Free Borders</h2>
8
9 \code
10 import salome
11 import geompy
12
13 import smesh
14
15 # create open shell: a box without one plane
16 box = geompy.MakeBox(0., 0., 0., 20., 20., 15.)
17 FaceList = geompy.SubShapeAll(box, geompy.ShapeType["FACE"])
18 FaceList.remove(FaceList[5])
19 box = geompy.MakeShell(FaceList)
20 idbox = geompy.addToStudy(box, "box")
21
22 # create a mesh
23 mesh = smesh.Mesh(box, "Mesh_free_borders")
24 algo = mesh.Segment()
25 algo.NumberOfSegments(5)
26 algo = mesh.Triangle()
27 algo.MaxElementArea(20.)
28 mesh.Compute() 
29
30 # criterion : free borders
31 aFilter = smesh.GetFilter(smesh.EDGE, smesh.FT_FreeBorders) 
32 anIds = mesh.GetIdsFromFilter(aFilter)
33
34 # print the result
35 print "Criterion: Free borders Nb = ", len(anIds)
36 j = 1
37 for i in range(len(anIds)):
38   if j > 20: j = 1; print ""
39   print anIds[i],
40   j = j + 1
41   pass
42 print ""
43
44 # create a group
45 aGroup = mesh.CreateGroup(SMESH.EDGE, "Free borders")
46 aGroup.Add(anIds)
47
48 salome.sg.updateObjBrowser(1)
49 \endcode
50
51 <br>
52 \anchor tui_borders_at_multiconnection
53 <h2>Borders at Multiconnection</h2>
54
55 \code
56 import salome
57 import geompy
58
59 import smesh
60 import SMESH
61
62 # create open shell: a box without one plane
63 box = geompy.MakeBox(0., 0., 0., 20., 20., 15.)
64 FaceList = geompy.SubShapeAll(box, geompy.ShapeType["FACE"])
65 FaceList.remove(FaceList[5])
66 box = geompy.MakeShell(FaceList)
67 idbox = geompy.addToStudy(box, "box")
68
69 # create a mesh
70 mesh = smesh.Mesh(box, "Mesh_borders_at_multi-connections")
71 algo = mesh.Segment()
72 algo.NumberOfSegments(5)
73 algo = mesh.Triangle()
74 algo.MaxElementArea(20.)
75 mesh.Compute() 
76
77 # Criterion : Borders at multi-connection
78 nb_conn = 2
79
80 aFilter = smesh.GetFilter(smesh.EDGE, smesh.FT_MultiConnection, smesh.FT_EqualTo, nb_conn)
81 anIds = mesh.GetIdsFromFilter(aFilter)
82
83 # print the result
84 print "Criterion: Borders at multi-connections Nb = ", len(anIds)
85 j = 1
86 for i in range(len(anIds)):
87   if j > 20: j = 1; print ""
88   print anIds[i],
89   j = j + 1
90   pass
91 print ""
92
93 # create a group
94 aGroup = mesh.CreateGroup(SMESH.EDGE, "Borders at multi-connections")
95 aGroup.Add(anIds)
96
97 salome.sg.updateObjBrowser(1)
98 \endcode
99
100 <br>
101 \anchor tui_length_1d
102 <h2>Length 1D</h2>
103
104 \code
105 import salome
106 import geompy
107
108 import smesh
109
110 # create open shell: a box without one plane
111 box = geompy.MakeBox(0., 0., 0., 20., 20., 15.)
112 FaceList = geompy.SubShapeAll(box, geompy.ShapeType["FACE"])
113 FaceList.remove(FaceList[5])
114 box = geompy.MakeShell(FaceList)
115 idbox = geompy.addToStudy(box, "box")
116
117 # create a mesh
118 mesh = smesh.Mesh(box, "Mesh_Length_1D")
119 algo = mesh.Segment()
120 algo.NumberOfSegments(5)
121 algo = mesh.Triangle()
122 algo.MaxElementArea(20.)
123 mesh.Compute() 
124
125 # Criterion : Length > 3.
126 length_margin = 3.
127
128 aFilter = smesh.GetFilter(smesh.EDGE, smesh.FT_Length, smesh.FT_MoreThan, length_margin)
129 anIds = mesh.GetIdsFromFilter(aFilter) 
130
131 # print the result
132 print "Criterion: Edges length > ", length_margin, " Nb = ", len(anIds)
133 j = 1
134 for i in range(len(anIds)):
135   if j > 20: j = 1; print ""
136   print anIds[i],
137   j = j + 1
138   pass
139 print ""
140
141 # create a group
142 aGroup = mesh.CreateGroup(SMESH.EDGE, "Edges with length > " + `length_margin`)
143 aGroup.Add(anIds)
144
145 salome.sg.updateObjBrowser(1)
146 \endcode
147
148 <br>
149 \anchor tui_free_edges
150 <h2>Free Edges</h2>
151
152 \code
153 import SMESH_mechanic
154
155 smesh  = SMESH_mechanic.smesh
156 mesh   = SMESH_mechanic.mesh
157 salome = SMESH_mechanic.salome
158
159 aFilterMgr = smesh.CreateFilterManager()
160
161 # Remove some elements to obtain free edges
162 # Criterion : AREA > 95.
163 area_margin = 95.
164
165 aFilter = smesh.GetFilter(smesh.FACE, smesh.FT_Area, smesh.FT_MoreThan, area_margin)
166
167 anIds = mesh.GetIdsFromFilter(aFilter)
168
169 mesh.RemoveElements(anIds)
170
171 # Criterion : Free Edges
172 aBorders = mesh.GetFreeBorders() 
173
174 # create groups
175 aGroupF = mesh.CreateGroup(SMESH.FACE, "Faces with free edges")
176 aGroupN = mesh.CreateGroup(SMESH.NODE, "Nodes on free edges")
177
178 # fill groups with elements, corresponding to the criterion
179 print ""
180 print "Criterion: Free edges Nb = ", len(aBorders)
181 for i in range(len(aBorders)):
182   aBorder = aBorders[i]
183   print "Face # ", aBorder.myElemId, " : Edge between nodes (",
184   print aBorder.myPnt1, ", ", aBorder.myPnt2, ")"
185
186   aGroupF.Add([aBorder.myElemId])
187   aGroupN.Add([aBorder.myPnt1, aBorder.myPnt2])
188
189 salome.sg.updateObjBrowser(1)
190 \endcode
191
192 <br>
193 \anchor tui_length_2d
194 <h2>Length 2D</h2>
195
196 \code
197 import salome
198 import geompy
199
200 import smesh
201
202 # create open shell: a box without one plane
203 box = geompy.MakeBox(0., 0., 0., 20., 20., 15.)
204 FaceList = geompy.SubShapeAll(box, geompy.ShapeType["FACE"])
205 FaceList.remove(FaceList[5])
206 box = geompy.MakeShell(FaceList)
207 idbox = geompy.addToStudy(box, "box")
208
209 # create a mesh
210 mesh = smesh.Mesh(box, "Mesh_Length_2D")
211 algo = mesh.Segment()
212 algo.NumberOfSegments(5)
213 algo = mesh.Triangle()
214 algo.MaxElementArea(20.)
215 mesh.Compute()
216
217 # Criterion : Length 2D > 5.7
218 length_margin = 5.7
219
220 aFilter = smesh.GetFilter(smesh.FACE, smesh.FT_Length2D, smesh.FT_MoreThan, length_margin)
221
222 anIds = mesh.GetIdsFromFilter(aFilter)
223
224 # print the result
225 print "Criterion: Edges length 2D > ", length_margin, " Nb = ", len(anIds)
226 j = 1
227 for i in range(len(anIds)):
228   if j > 20: j = 1; print ""
229   print anIds[i],
230   j = j + 1
231   pass
232 print ""
233
234 # create a group
235 aGroup = mesh.CreateEmptyGroup(smesh.FACE, "Faces with length 2D > " + `length_margin`)
236 aGroup.Add(anIds)
237
238 salome.sg.updateObjBrowser(1)
239 \endcode
240
241 <br>
242 \anchor tui_borders_at_multiconnection_2d
243 <h2>Borders at Multiconnection 2D</h2>
244
245 \code
246 import salome
247 import geompy
248
249 import smesh
250
251 # create a compound of two glued boxes
252 box1 = geompy.MakeBox(0., 0., 0., 20., 20., 15.)
253 box2 = geompy.MakeTranslation(box1, 0., 20., 0)
254 comp = geompy.MakeCompound([box1, box2])
255 box = geompy.MakeGlueFaces(comp, 0.000001)
256 idbox = geompy.addToStudy(box, "box")
257
258 # create a mesh
259 mesh = smesh.Mesh(box, "Box compound : 2D triangle mesh")
260 algo = mesh.Segment()
261 algo.NumberOfSegments(5)
262 algo = mesh.Triangle()
263 algo.MaxElementArea(20.)
264 mesh.Compute() 
265
266 # Criterion : MULTI-CONNECTION 2D = 3
267 nb_conn = 3
268
269 aFilter = smesh.GetFilter(smesh.FACE, smesh.FT_MultiConnection2D, smesh.FT_EqualTo, nb_conn)
270
271 anIds = mesh.GetIdsFromFilter(aFilter) 
272
273 # print the result
274 print "Criterion: Borders at multi-connection 2D = ", nb_conn, " Nb = ", len(anIds)
275 j = 1
276 for i in range(len(anIds)):
277   if j > 20: j = 1; print ""
278   print anIds[i],
279   j = j + 1
280   pass
281 print ""
282
283 # create a group
284 aGroup = mesh.CreateEmptyGroup(smesh.FACE, "Borders at multi-connection 2D = " + `nb_conn`)
285 aGroup.Add(anIds)
286
287 salome.sg.updateObjBrowser(1)
288 \endcode
289
290 <br>
291 \anchor tui_area
292 <h2>Area</h2>
293
294 \code
295 import SMESH_mechanic
296
297 smesh  = SMESH_mechanic.smesh
298 mesh   = SMESH_mechanic.mesh
299 salome = SMESH_mechanic.salome
300
301 # Criterion : AREA > 100.
302 area_margin = 100.
303
304 aFilter = smesh.GetFilter(smesh.FACE, smesh.FT_Area, smesh.FT_MoreThan, area_margin)
305
306 anIds = mesh.GetIdsFromFilter(aFilter) 
307
308 # print the result
309 print "Criterion: Area > ", area_margin, " Nb = ", len(anIds)
310 j = 1
311 for i in range(len(anIds)):
312   if j > 20: j = 1; print ""
313   print anIds[i],
314   j = j + 1
315   pass
316 print ""
317
318 # create a group
319 aGroup = mesh.CreateEmptyGroup(smesh.FACE, "Area > " + `area_margin`)
320 aGroup.Add(anIds)
321
322 salome.sg.updateObjBrowser(1)  
323 \endcode
324
325 <br>
326 \anchor tui_taper
327 <h2>Taper</h2>
328
329 \code
330 import SMESH_mechanic
331
332 smesh  = SMESH_mechanic.smesh
333 mesh   = SMESH_mechanic.mesh
334 salome = SMESH_mechanic.salome
335
336 # Criterion : Taper > 3e-20
337 taper_margin = 3e-20
338
339 aFilter = smesh.GetFilter(smesh.FACE, smesh.FT_Taper, smesh.FT_MoreThan, taper_margin)
340
341 anIds = mesh.GetIdsFromFilter(aFilter) 
342
343 # print the result
344 print "Criterion: Taper > ", taper_margin, " Nb = ", len(anIds)
345 j = 1
346 for i in range(len(anIds)):
347   if j > 20: j = 1; print ""
348   print anIds[i],
349   j = j + 1
350   pass
351 print ""
352
353 # create a group
354 aGroup = mesh.CreateEmptyGroup(smesh.FACE, "Taper > " + `taper_margin`)
355 aGroup.Add(anIds)
356
357 salome.sg.updateObjBrowser(1)
358 \endcode
359
360 <br>
361 \anchor tui_aspect_ratio
362 <h2>Aspect Ratio</h2>
363
364 \code
365 import SMESH_mechanic
366
367 smesh  = SMESH_mechanic.smesh
368 mesh   = SMESH_mechanic.mesh
369 salome = SMESH_mechanic.salome
370
371 # Criterion : ASPECT RATIO > 1.8
372 ar_margin = 1.8
373
374 aFilter = smesh.GetFilter(smesh.FACE, smesh.FT_AspectRatio, smesh.FT_MoreThan, ar_margin)
375
376 anIds = mesh.GetIdsFromFilter(aFilter) 
377
378 # print the result
379 print "Criterion: Aspect Ratio > ", ar_margin, " Nb = ", len(anIds)
380 j = 1
381 for i in range(len(anIds)):
382   if j > 20: j = 1; print ""
383   print anIds[i],
384   j = j + 1
385   pass
386 print ""
387
388 # create a group
389 aGroup = mesh.CreateEmptyGroup(smesh.FACE, "Aspect Ratio > " + `ar_margin`)
390 aGroup.Add(anIds)
391
392 salome.sg.updateObjBrowser(1)
393 \endcode
394
395 <br>
396 \anchor tui_minimum_angle
397 <h2>Minimum Angle</h2>
398
399 \code
400 import SMESH_mechanic
401
402 smesh  = SMESH_mechanic.smesh
403 mesh   = SMESH_mechanic.mesh
404 salome = SMESH_mechanic.salome
405
406 # Criterion : MINIMUM ANGLE < 35.
407 min_angle = 35.
408
409 aFilter = smesh.GetFilter(smesh.FACE, smesh.FT_MinimumAngle, smesh.FT_LessThan, min_angle)
410
411 anIds = mesh.GetIdsFromFilter(aFilter) 
412
413 # print the result
414 print "Criterion: Minimum Angle < ", min_angle, " Nb = ", len(anIds)
415 j = 1
416 for i in range(len(anIds)):
417   if j > 20: j = 1; print ""
418   print anIds[i],
419   j = j + 1
420   pass
421 print ""
422
423 # create a group
424 aGroup = mesh.CreateEmptyGroup(smesh.FACE, "Minimum Angle < " + `min_angle`)
425
426 aGroup.Add(anIds)
427
428 salome.sg.updateObjBrowser(1)
429 \endcode
430
431 <br>
432 \anchor tui_warping
433 <h2>Warping</h2>
434
435 \code
436 import SMESH_mechanic
437
438 smesh  = SMESH_mechanic.smesh
439 mesh   = SMESH_mechanic.mesh
440 salome = SMESH_mechanic.salome
441
442 # Criterion : WARP ANGLE > 1e-15
443 wa_margin = 1e-15
444
445 aFilter = smesh.GetFilter(smesh.FACE, smesh.FT_Warping, smesh.FT_MoreThan, wa_margin)
446
447 anIds = mesh.GetIdsFromFilter(aFilter) 
448
449 # print the result
450 print "Criterion: Warp > ", wa_margin, " Nb = ", len(anIds)
451 j = 1
452 for i in range(len(anIds)):
453   if j > 20: j = 1; print ""
454   print anIds[i],
455   j = j + 1
456   pass
457 print ""
458
459 # create a group
460 aGroup = mesh.CreateEmptyGroup(smesh.FACE, "Warp > " + `wa_margin`)
461
462 aGroup.Add(anIds)
463
464 salome.sg.updateObjBrowser(1) 
465 \endcode
466
467 <br>
468 \anchor tui_skew
469 <h2>Skew</h2>
470
471 \code
472 import SMESH_mechanic
473
474 smesh  = SMESH_mechanic.smesh
475 mesh   = SMESH_mechanic.mesh
476 salome = SMESH_mechanic.salome
477
478 # Criterion : Skew > 38.
479 skew_margin = 38.
480
481 aFilter = smesh.GetFilter(smesh.FACE, smesh.FT_Skew, smesh.FT_MoreThan, skew_margin)
482
483 anIds = mesh.GetIdsFromFilter(aFilter)
484
485 # print the result
486 print "Criterion: Skew > ", skew_margin, " Nb = ", len(anIds)
487 j = 1
488 for i in range(len(anIds)):
489   if j > 20: j = 1; print ""
490   print anIds[i],
491   j = j + 1
492   pass
493 print ""
494
495 # create a group
496 aGroup = mesh.CreateEmptyGroup(smesh.FACE, "Skew > " + `skew_margin`)
497 aGroup.Add(anIds)
498
499 salome.sg.updateObjBrowser(1)
500 \endcode
501
502 <br>
503 \anchor tui_aspect_ratio_3d
504 <h2>Aspect Ratio 3D</h2>
505
506 \code
507 import SMESH_mechanic_tetra 
508
509 smesh  = SMESH_mechanic_tetra.smesh
510 mesh   = SMESH_mechanic_tetra.mesh
511 salome = SMESH_mechanic_tetra.salome
512
513 # Criterion : ASPECT RATIO 3D > 4.5
514 ar_margin = 4.5
515
516 aFilter = smesh.GetFilter(smesh.VOLUME, smesh.FT_AspectRatio3D, smesh.FT_MoreThan, ar_margin)
517
518 anIds = mesh.GetIdsFromFilter(aFilter) 
519
520 # print the result
521 print "Criterion: Aspect Ratio 3D > ", ar_margin, " Nb = ", len(anIds)
522 j = 1
523 for i in range(len(anIds)):
524   if j > 20: j = 1; print ""
525   print anIds[i],
526   j = j + 1
527   pass
528 print ""
529
530 # create a group
531 aGroup = mesh.CreateEmptyGroup(smesh.VOLUME, "Aspect Ratio 3D > " + `ar_margin`)
532
533 aGroup.Add(anIds)
534
535 salome.sg.updateObjBrowser(1)
536 \endcode
537
538 <br>
539 \anchor tui_volume
540 <h2>Volume</h2>
541
542 \code
543 import SMESH_mechanic_tetra
544
545 smesh  = SMESH_mechanic_tetra.smesh
546 mesh   = SMESH_mechanic_tetra.mesh
547 salome = SMESH_mechanic_tetra.salome
548
549 # Criterion : VOLUME < 7.
550 volume_margin = 7.
551
552 aFilter = smesh.GetFilter(smesh.VOLUME, smesh.FT_Volume3D, smesh.FT_LessThan, volume_margin)
553
554 anIds = mesh.GetIdsFromFilter(aFilter) 
555
556 # print the result
557 print ""
558 print "Criterion: Volume < ", volume_margin, " Nb = ", len(anIds)
559 j = 1
560 for i in range(len(anIds)):
561   if j > 20: j = 1; print ""
562   print anIds[i],
563   j = j + 1
564   pass
565 print ""
566
567 # create a group
568 aGroup = mesh.CreateEmptyGroup(smesh.VOLUME, "Volume < " + `volume_margin`)
569
570 aGroup.Add(anIds)
571
572 salome.sg.updateObjBrowser(1) 
573 \endcode
574
575 */