1 <!doctype HTML public "-//W3C//DTD HTML 4.0 Frameset//EN">
\r
6 <title>Quality Controls</title>
\r
7 <meta http-equiv="content-type" content="text/html; charset=windows-1252">
\r
8 <meta name="generator" content="RoboHelp by eHelp Corporation www.ehelp.com"><style type="text/css">
\r
10 p.whs1 { margin-top:0pt; margin-bottom:0pt; font-family:'Lucida Console' , monospace; }
\r
11 p.whs2 { margin-top:0pt; margin-bottom:0pt; }
\r
12 p.whs3 { margin-top:0pt; margin-bottom:0pt; font-family:'Times New Roman' , serif; }
\r
14 </style><script type="text/javascript" language="JavaScript">
\r
16 if ((navigator.appName == "Netscape") && (parseInt(navigator.appVersion) == 4))
\r
18 var strNSS = "<style type='text/css'>";
\r
19 strNSS += "p.whs1 {margin-top:1pt;margin-bottom:1pt; }";
\r
20 strNSS += "p.whs2 {margin-top:1pt;margin-bottom:1pt; }";
\r
21 strNSS += "p.whs3 {margin-top:1pt;margin-bottom:1pt; }";
\r
22 strNSS +="</style>";
\r
23 document.write(strNSS);
\r
27 <script type="text/javascript" language="JavaScript" title="WebHelpInlineScript">
\r
30 if (innerWidth != origWidth || innerHeight != origHeight)
\r
33 if ((parseInt(navigator.appVersion) == 4) && (navigator.appName == "Netscape")) {
\r
34 origWidth = innerWidth;
\r
35 origHeight = innerHeight;
\r
41 <style type="text/css">
\r
43 div.WebHelpPopupMenu { position:absolute; left:0px; top:0px; z-index:4; visibility:hidden; }
\r
44 p.WebHelpNavBar { text-align:right; }
\r
46 </style><script type="text/javascript" language="javascript1.2" src="whmsg.js"></script>
\r
47 <script type="text/javascript" language="javascript" src="whver.js"></script>
\r
48 <script type="text/javascript" language="javascript1.2" src="whproxy.js"></script>
\r
49 <script type="text/javascript" language="javascript1.2" src="whutils.js"></script>
\r
50 <script type="text/javascript" language="javascript1.2" src="whtopic.js"></script>
\r
51 <script type="text/javascript" language="javascript1.2">
\r
53 if (window.gbWhTopic)
\r
55 if (window.setRelStartPage)
\r
57 addTocInfo("MESH module\nTUI Scripts\nQuality Controls");
\r
58 addButton("show",BTN_IMG,"Show","","","","",0,0,"whd_show0.gif","whd_show2.gif","whd_show1.gif");
\r
59 addButton("hide",BTN_IMG,"Hide","","","","",0,0,"whd_hide0.gif","whd_hide2.gif","whd_hide1.gif");
\r
64 if (window.setRelStartPage)
\r
66 setRelStartPage("index.htm");
\r
76 document.location.reload();
\r
80 <body><script type="text/javascript" language="javascript1.2">
\r
82 if (window.writeIntopicBar)
\r
86 <h1>Quality Controls</h1>
\r
88 <h3><a name=bookmark>Free Borders</a></h3>
\r
90 <p class="whs1">import salome</p>
\r
92 <p class="whs1">import geompy</p>
\r
94 <p class="whs1"> </p>
\r
96 <p class="whs1">import smesh</p>
\r
98 <p class="whs2"> </p>
\r
100 <p class="whs2"># create open shell: a
\r
101 box without one plane</p>
\r
103 <p class="whs1"><span style="font-family: 'Lucida Console', monospace;">box
\r
104 = geompy.MakeBox(0., 0., 0., 20., 20., 15.)</span></p>
\r
106 <p class="whs1">FaceList = geompy.SubShapeAll(box,
\r
107 geompy.ShapeType["FACE"])</p>
\r
109 <p class="whs1">FaceList.remove(FaceList[5])</p>
\r
111 <p class="whs1">box = geompy.MakeShell(FaceList)</p>
\r
113 <p class="whs1">idbox = geompy.addToStudy(box,
\r
114 "box")</p>
\r
116 <p class="whs1"> </p>
\r
118 <p class="whs2"># create a mesh</p>
\r
120 <p class="whs1"><span style="font-family: 'Lucida Console', monospace;">mesh
\r
121 = smesh.Mesh(box, "Mesh_free_borders")</span></p>
\r
123 <p class="whs1">algo = mesh.Segment()</p>
\r
125 <p class="whs1">algo.NumberOfSegments(5)</p>
\r
127 <p class="whs1">algo = mesh.Triangle()</p>
\r
129 <p class="whs1">algo.MaxElementArea(20.)</p>
\r
131 <p class="whs1">mesh.Compute() </p>
\r
133 <p class="whs1"> </p>
\r
135 <p class="whs2"># criterion : free borders</p>
\r
137 <p class="whs1"><span style="font-family: 'Lucida Console', monospace;">aFilter
\r
138 = smesh.GetFilter(smesh.EDGE, smesh.FT_FreeBorders)</span><span style="margin-top: 0pt;
\r
139 margin-bottom: 0pt;
\r
140 font-family: 'Lucida Console', monospace;"> </span></p>
\r
142 <p class="whs1">anIds = mesh.GetIdsFromFilter(aFilter)</p>
\r
144 <p class="whs1"> </p>
\r
146 <p class="whs2"> </p>
\r
148 <p class="whs2"># print the result</p>
\r
150 <p class="whs1">print "Criterion:
\r
151 Free borders Nb = ", len(anIds)</p>
\r
153 <p class="whs1">j = 1</p>
\r
155 <p class="whs1">for i in range(len(anIds)):</p>
\r
157 <p class="whs1"> if
\r
158 j > 20: j = 1; print ""</p>
\r
160 <p class="whs1"> print
\r
163 <p class="whs1"> j
\r
166 <p class="whs1"> pass</p>
\r
168 <p class="whs1">print ""</p>
\r
170 <p class="whs2"> </p>
\r
172 <p class="whs2"># create a group</p>
\r
174 <p class="whs1">aGroup = mesh.CreateGroup(SMESH.EDGE,
\r
175 "Free borders")</p>
\r
177 <p class="whs1">aGroup.Add(anIds)</p>
\r
179 <p class="whs1"> </p>
\r
181 <p class="whs2"><span style="font-family: 'Lucida Console', monospace;">salome.sg.updateObjBrowser(1)</span>
\r
184 <h3><a name=bookmark1>Borders at Multiconnection</a></h3>
\r
186 <p class="whs1">import salome</p>
\r
188 <p class="whs1">import geompy</p>
\r
190 <p class="whs1"> </p>
\r
192 <p class="whs1">import smesh</p>
\r
194 <p class="whs1">import SMESH</p>
\r
196 <p class="whs1"> </p>
\r
198 <p class="whs3"># create open shell:
\r
199 a box without one plane</p>
\r
201 <p class="whs1">box = geompy.MakeBox(0.,
\r
202 0., 0., 20., 20., 15.)</p>
\r
204 <p class="whs1">FaceList = geompy.SubShapeAll(box,
\r
205 geompy.ShapeType["FACE"])</p>
\r
207 <p class="whs1">FaceList.remove(FaceList[5])</p>
\r
209 <p class="whs1">box = geompy.MakeShell(FaceList)</p>
\r
211 <p class="whs1">idbox = geompy.addToStudy(box,
\r
212 "box")</p>
\r
214 <p class="whs1"> </p>
\r
216 <p class="whs3"># create a mesh</p>
\r
218 <p class="whs1"><span style="font-family: 'Lucida Console', monospace;">mesh
\r
219 = smesh.Mesh(box, "Mesh_borders_at_multi-connections")</span></p>
\r
221 <p class="whs1">algo = mesh.Segment()</p>
\r
223 <p class="whs1">algo.NumberOfSegments(5)</p>
\r
225 <p class="whs1">algo = mesh.Triangle()</p>
\r
227 <p class="whs1">algo.MaxElementArea(20.)</p>
\r
229 <p class="whs1">mesh.Compute() </p>
\r
231 <p class="whs1"> </p>
\r
233 <p class="whs3"># Criterion : Borders
\r
234 at multi-connection</p>
\r
236 <p class="whs1">nb_conn = 2</p>
\r
238 <p class="whs1"> </p>
\r
240 <p class="whs1"><span style="font-family: 'Lucida Console', monospace;">aFilter
\r
241 = smesh.GetFilter(smesh.EDGE, smesh.FT_MultiConnection, smesh.FT_EqualTo,
\r
242 nb_conn)</span></p>
\r
244 <p class="whs1">anIds = mesh.GetIdsFromFilter(aFilter)</p>
\r
246 <p class="whs1"> </p>
\r
248 <p class="whs1"> </p>
\r
250 <p class="whs3"># print the result</p>
\r
252 <p class="whs1">print "Criterion:
\r
253 Borders at multi-connections Nb = ", len(anIds)</p>
\r
255 <p class="whs1">j = 1</p>
\r
257 <p class="whs1">for i in range(len(anIds)):</p>
\r
259 <p class="whs1"> if
\r
260 j > 20: j = 1; print ""</p>
\r
262 <p class="whs1"> print
\r
265 <p class="whs1"> j
\r
268 <p class="whs1"> pass</p>
\r
270 <p class="whs1">print ""</p>
\r
272 <p class="whs1"> </p>
\r
274 <p class="whs3"># create a group</p>
\r
276 <p class="whs1">aGroup = mesh.CreateGroup(SMESH.EDGE,
\r
277 "Borders at multi-connections")</p>
\r
279 <p class="whs1">aGroup.Add(anIds)</p>
\r
281 <p class="whs1">salome.sg.updateObjBrowser(1)
\r
284 <h3><a name=bookmark2>Length 1D</a></h3>
\r
286 <p class="whs1">import salome</p>
\r
288 <p class="whs1">import geompy</p>
\r
290 <p class="whs1"> </p>
\r
292 <p class="whs1">import smesh</p>
\r
294 <p class="whs2"> </p>
\r
296 <p class="whs2"># create open shell: a
\r
297 box without one plane</p>
\r
299 <p class="whs1">box = geompy.MakeBox(0.,
\r
300 0., 0., 20., 20., 15.)</p>
\r
302 <p class="whs1">FaceList = geompy.SubShapeAll(box,
\r
303 geompy.ShapeType["FACE"])</p>
\r
305 <p class="whs1">FaceList.remove(FaceList[5])</p>
\r
307 <p class="whs1">box = geompy.MakeShell(FaceList)</p>
\r
309 <p class="whs1">idbox = geompy.addToStudy(box,
\r
310 "box")</p>
\r
312 <p class="whs2"> </p>
\r
314 <p class="whs2"># create a mesh</p>
\r
316 <p class="whs1"><span style="font-family: 'Lucida Console', monospace;">mesh
\r
317 = smesh.Mesh(box, "Mesh_Length_1D")</span></p>
\r
319 <p class="whs1">algo = mesh.Segment()</p>
\r
321 <p class="whs1">algo.NumberOfSegments(5)</p>
\r
323 <p class="whs1">algo = mesh.Triangle()</p>
\r
325 <p class="whs1">algo.MaxElementArea(20.)</p>
\r
327 <p class="whs1">mesh.Compute() </p>
\r
329 <p class="whs2"> </p>
\r
331 <p class="whs2"># Criterion : Length >
\r
334 <p class="whs1">length_margin =
\r
337 <p class="whs1"> </p>
\r
339 <p class="whs1"><span style="font-family: 'Lucida Console', monospace;">aFilter
\r
340 = smesh.GetFilter(smesh.EDGE, smesh.FT_Length, smesh.FT_MoreThan, length_margin)</span></p>
\r
342 <p class="whs1">anIds = mesh.GetIdsFromFilter(aFilter)
\r
345 <p class="whs2"> </p>
\r
347 <p class="whs2"># print the result</p>
\r
349 <p class="whs1">print "Criterion:
\r
350 Edges length > ", length_margin, " Nb = ", len(anIds)</p>
\r
352 <p class="whs1">j = 1</p>
\r
354 <p class="whs1">for i in range(len(anIds)):</p>
\r
356 <p class="whs1"> if
\r
357 j > 20: j = 1; print ""</p>
\r
359 <p class="whs1"> print
\r
362 <p class="whs1"> j
\r
365 <p class="whs1"> pass</p>
\r
367 <p class="whs1">print ""</p>
\r
369 <p class="whs2"> </p>
\r
371 <p class="whs2"># create a group</p>
\r
373 <p class="whs1">aGroup = mesh.CreateGroup(SMESH.EDGE,
\r
374 "Edges with length > " + `length_margin`)</p>
\r
376 <p class="whs1">aGroup.Add(anIds)</p>
\r
378 <p class="whs1"> </p>
\r
380 <p class="whs2"><span style="font-family: 'Lucida Console', monospace;">salome.sg.updateObjBrowser(1)</span>
\r
383 <h3><a name=bookmark3>Free Edges</a></h3>
\r
385 <p class="whs1">import SMESH_mechanic</p>
\r
387 <p class="whs1"> </p>
\r
389 <p class="whs1">smesh =
\r
390 SMESH_mechanic.smesh</p>
\r
392 <p class="whs1">mesh =
\r
393 SMESH_mechanic.mesh</p>
\r
395 <p class="whs1">salome = SMESH_mechanic.salome</p>
\r
397 <p class="whs1"> </p>
\r
399 <p class="whs1">aFilterMgr = smesh.CreateFilterManager()</p>
\r
401 <p class="whs1"> </p>
\r
403 <p class="whs3"># Remove some elements
\r
404 to obtain free edges</p>
\r
406 <p class="whs3"># Criterion : AREA
\r
409 <p class="whs1">area_margin = 95.</p>
\r
411 <p class="whs1"> </p>
\r
413 <p class="whs1"><span style="font-family: 'Lucida Console', monospace;">aFilter
\r
414 = smesh.GetFilter(smesh.FACE, smesh.FT_Area, smesh.FT_MoreThan, area_margin)</span></p>
\r
416 <p class="whs1"> </p>
\r
418 <p class="whs1">anIds = mesh.GetIdsFromFilter(aFilter)</p>
\r
420 <p class="whs1"> </p>
\r
422 <p class="whs1">mesh.RemoveElements(anIds)</p>
\r
424 <p class="whs1"> </p>
\r
426 <p class="whs3"># Criterion : Free
\r
429 <p class="whs1">aBorders = mesh.GetFreeBorders()
\r
432 <p class="whs1"> </p>
\r
434 <p class="whs3"># create groups</p>
\r
436 <p class="whs1">aGroupF = mesh.CreateGroup(SMESH.FACE,
\r
437 "Faces with free edges")</p>
\r
439 <p class="whs1">aGroupN = mesh.CreateGroup(SMESH.NODE,
\r
440 "Nodes on free edges")</p>
\r
442 <p class="whs1"> </p>
\r
444 <p class="whs3"># fill groups with
\r
445 elements, corresponding to the criterion</p>
\r
447 <p class="whs1">print ""</p>
\r
449 <p class="whs1">print "Criterion:
\r
450 Free edges Nb = ", len(aBorders)</p>
\r
452 <p class="whs1">for i in range(len(aBorders)):</p>
\r
454 <p class="whs1"> aBorder
\r
457 <p class="whs1"> print
\r
458 "Face # ", aBorder.myElemId, " : Edge between nodes (",</p>
\r
460 <p class="whs1"> print
\r
461 aBorder.myPnt1, ", ", aBorder.myPnt2, ")"</p>
\r
463 <p class="whs1"> </p>
\r
465 <p class="whs1"> <span
\r
466 style="margin-top: 0pt;
\r
467 margin-bottom: 0pt;
\r
468 font-family: 'Lucida Console', monospace;">aGroupF.Add([aBorder.myElemId])</span></p>
\r
470 <p class="whs1"> aGroupN.Add([aBorder.myPnt1,
\r
471 aBorder.myPnt2])</p>
\r
473 <p class="whs1"> </p>
\r
475 <p class="whs1">salome.sg.updateObjBrowser(1)
\r
478 <p class="whs1"> </p>
\r
480 <h3><a name=bookmark4>Length 2D</a></h3>
\r
482 <p class="whs1">import salome</p>
\r
484 <p class="whs1">import geompy</p>
\r
486 <p class="whs1"> </p>
\r
488 <p class="whs1">import smesh</p>
\r
490 <p class="whs1"> </p>
\r
492 <p class="whs3"># create open shell:
\r
493 a box without one plane</p>
\r
495 <p class="whs1">box = geompy.MakeBox(0.,
\r
496 0., 0., 20., 20., 15.)</p>
\r
498 <p class="whs1">FaceList = geompy.SubShapeAll(box,
\r
499 geompy.ShapeType["FACE"])</p>
\r
501 <p class="whs1">FaceList.remove(FaceList[5])</p>
\r
503 <p class="whs1">box = geompy.MakeShell(FaceList)</p>
\r
505 <p class="whs1">idbox = geompy.addToStudy(box,
\r
506 "box")</p>
\r
508 <p class="whs1"> </p>
\r
510 <p class="whs3"># create a mesh</p>
\r
512 <p class="whs1"><span style="font-family: 'Lucida Console', monospace;">mesh
\r
513 = smesh.Mesh(box, "Mesh_Length_2D")</span></p>
\r
515 <p class="whs1">algo = mesh.Segment()</p>
\r
517 <p class="whs1">algo.NumberOfSegments(5)</p>
\r
519 <p class="whs1">algo = mesh.Triangle()</p>
\r
521 <p class="whs1">algo.MaxElementArea(20.)</p>
\r
523 <p class="whs1">mesh.Compute()</p>
\r
525 <p class="whs1"> </p>
\r
527 <p class="whs3"># Criterion : Length
\r
530 <p class="whs1">length_margin =
\r
533 <p class="whs1"> </p>
\r
535 <p class="whs1"><span style="font-family: 'Lucida Console', monospace;">aFilter
\r
536 = smesh.GetFilter(smesh.FACE, smesh.FT_Length2D, smesh.FT_MoreThan, length_margin)</span></p>
\r
538 <p class="whs1"> </p>
\r
540 <p class="whs1">anIds = mesh.GetIdsFromFilter(aFilter)</p>
\r
542 <p class="whs1"> </p>
\r
544 <p class="whs3"># print the result</p>
\r
546 <p class="whs1">print "Criterion:
\r
547 Edges length 2D > ", length_margin, " Nb = ", len(anIds)</p>
\r
549 <p class="whs1">j = 1</p>
\r
551 <p class="whs1">for i in range(len(anIds)):</p>
\r
553 <p class="whs1"> if
\r
554 j > 20: j = 1; print ""</p>
\r
556 <p class="whs1"> print
\r
559 <p class="whs1"> j
\r
562 <p class="whs1"> pass</p>
\r
564 <p class="whs1">print ""</p>
\r
566 <p class="whs1"> </p>
\r
568 <p class="whs3"># create a group</p>
\r
570 <p class="whs1"><span style="font-family: 'Lucida Console', monospace;">aGroup
\r
571 = mesh.CreateEmptyGroup(smesh.FACE, "Faces with length 2D > "
\r
572 + `length_margin`)</span></p>
\r
574 <p class="whs1">aGroup.Add(anIds)</p>
\r
576 <p class="whs1"> </p>
\r
578 <p class="whs1">salome.sg.updateObjBrowser(1)
\r
581 <p class="whs2"> </p>
\r
583 <h3><a name=bookmark5>Borders at Multiconnection 2D</a></h3>
\r
585 <p class="whs1">import salome</p>
\r
587 <p class="whs1">import geompy</p>
\r
589 <p class="whs1"> </p>
\r
591 <p class="whs1">import smesh</p>
\r
593 <p class="whs3"> </p>
\r
595 <p class="whs3"># create a compound
\r
596 of two glued boxes</p>
\r
598 <p class="whs1">box1 = geompy.MakeBox(0.,
\r
599 0., 0., 20., 20., 15.)</p>
\r
601 <p class="whs1">box2 = geompy.MakeTranslation(box1,
\r
604 <p class="whs1">comp = geompy.MakeCompound([box1,
\r
607 <p class="whs1">box = geompy.MakeGlueFaces(comp,
\r
610 <p class="whs1">idbox = geompy.addToStudy(box,
\r
611 "box")</p>
\r
613 <p class="whs1"> </p>
\r
615 <p class="whs3"># create a mesh</p>
\r
617 <p class="whs1"><span style="font-family: 'Lucida Console', monospace;">mesh
\r
618 = smesh.Mesh(box, "Box compound : 2D triangle mesh")</span></p>
\r
620 <p class="whs1">algo = mesh.Segment()</p>
\r
622 <p class="whs1">algo.NumberOfSegments(5)</p>
\r
624 <p class="whs1">algo = mesh.Triangle()</p>
\r
626 <p class="whs1">algo.MaxElementArea(20.)</p>
\r
628 <p class="whs1">mesh.Compute() </p>
\r
630 <p class="whs1"> </p>
\r
632 <p class="whs3"># Criterion : MULTI-CONNECTION
\r
635 <p class="whs1">nb_conn = 3</p>
\r
637 <p class="whs1"> </p>
\r
639 <p class="whs1"><span style="font-family: 'Lucida Console', monospace;">aFilter
\r
640 = smesh.GetFilter(smesh.FACE, smesh.FT_MultiConnection2D, smesh.FT_EqualTo,
\r
641 nb_conn)</span></p>
\r
643 <p class="whs1"> </p>
\r
645 <p class="whs1">anIds = mesh.GetIdsFromFilter(aFilter)
\r
648 <p class="whs1"> </p>
\r
650 <p class="whs3"># print the result</p>
\r
652 <p class="whs1">print "Criterion:
\r
653 Borders at multi-connection 2D = ", nb_conn, " Nb = ",
\r
656 <p class="whs1">j = 1</p>
\r
658 <p class="whs1">for i in range(len(anIds)):</p>
\r
660 <p class="whs1"> if
\r
661 j > 20: j = 1; print ""</p>
\r
663 <p class="whs1"> print
\r
666 <p class="whs1"> j
\r
669 <p class="whs1"> pass</p>
\r
671 <p class="whs1">print ""</p>
\r
673 <p class="whs1"> </p>
\r
675 <p class="whs3"># create a group</p>
\r
677 <p class="whs1"><span style="font-family: 'Lucida Console', monospace;">aGroup
\r
678 = mesh.CreateEmptyGroup(smesh.FACE, "Borders at multi-connection
\r
679 2D = " + `nb_conn`)</span></p>
\r
681 <p class="whs1">aGroup.Add(anIds)</p>
\r
683 <p class="whs1"> </p>
\r
685 <p class="whs1">salome.sg.updateObjBrowser(1)
\r
688 <p class="whs1"> </p>
\r
690 <h3><a name=bookmark6>Area</a></h3>
\r
692 <p class="whs1">import SMESH_mechanic</p>
\r
694 <p class="whs1"> </p>
\r
696 <p class="whs1">smesh =
\r
697 SMESH_mechanic.smesh</p>
\r
699 <p class="whs1">mesh =
\r
700 SMESH_mechanic.mesh</p>
\r
702 <p class="whs1">salome = SMESH_mechanic.salome</p>
\r
704 <p class="whs2"> </p>
\r
706 <p class="whs2"># Criterion : AREA >
\r
709 <p class="whs1">area_margin = 100.</p>
\r
711 <p class="whs1"> </p>
\r
713 <p class="whs1">aFilter = smesh.GetFilter(smesh.FACE,
\r
714 smesh.FT_Area, smesh.FT_MoreThan, area_margin)</p>
\r
716 <p class="whs1"> </p>
\r
718 <p class="whs1">anIds = mesh.GetIdsFromFilter(aFilter)
\r
721 <p class="whs2"> </p>
\r
723 <p class="whs2"># print the result</p>
\r
725 <p class="whs1">print "Criterion:
\r
726 Area > ", area_margin, " Nb = ", len(anIds)</p>
\r
728 <p class="whs1">j = 1</p>
\r
730 <p class="whs1">for i in range(len(anIds)):</p>
\r
732 <p class="whs1"> if
\r
733 j > 20: j = 1; print ""</p>
\r
735 <p class="whs1"> print
\r
738 <p class="whs1"> j
\r
741 <p class="whs1"> pass</p>
\r
743 <p class="whs1">print ""</p>
\r
745 <p class="whs1"> </p>
\r
747 <p class="whs2"># create a group</p>
\r
749 <p class="whs1">aGroup = mesh.CreateEmptyGroup(smesh.FACE,
\r
750 "Area > " + `area_margin`)</p>
\r
752 <p class="whs1">aGroup.Add(anIds)</p>
\r
754 <p class="whs1"> </p>
\r
756 <p class="whs2"><span style="font-family: 'Lucida Console', monospace;">salome.sg.updateObjBrowser(1)</span>
\r
759 <h3><a name=bookmark7>Taper</a></h3>
\r
761 <p class="whs1">import SMESH_mechanic</p>
\r
763 <p class="whs1"> </p>
\r
765 <p class="whs1">smesh =
\r
766 SMESH_mechanic.smesh</p>
\r
768 <p class="whs1">mesh =
\r
769 SMESH_mechanic.mesh</p>
\r
771 <p class="whs1">salome = SMESH_mechanic.salome</p>
\r
773 <p class="whs1"> </p>
\r
775 <p class="whs3"># Criterion : Taper
\r
778 <p class="whs1">taper_margin = 3e-20</p>
\r
780 <p class="whs1"> </p>
\r
782 <p class="whs1">aFilter = smesh.GetFilter(smesh.FACE,
\r
783 smesh.FT_Taper, smesh.FT_MoreThan, taper_margin)</p>
\r
785 <p class="whs1"> </p>
\r
787 <p class="whs1">anIds = mesh.GetIdsFromFilter(aFilter)
\r
790 <p class="whs1"> </p>
\r
792 <p class="whs3"># print the result</p>
\r
794 <p class="whs1">print "Criterion:
\r
795 Taper > ", taper_margin, " Nb = ", len(anIds)</p>
\r
797 <p class="whs1">j = 1</p>
\r
799 <p class="whs1">for i in range(len(anIds)):</p>
\r
801 <p class="whs1"> if
\r
802 j > 20: j = 1; print ""</p>
\r
804 <p class="whs1"> print
\r
807 <p class="whs1"> j
\r
810 <p class="whs1"> pass</p>
\r
812 <p class="whs1">print ""</p>
\r
814 <p class="whs1"> </p>
\r
816 <p class="whs3"># create a group</p>
\r
818 <p class="whs1"><span style="font-family: 'Lucida Console', monospace;">aGroup
\r
819 = mesh.CreateEmptyGroup(smesh.FACE, "Taper > " + `taper_margin`)</span></p>
\r
821 <p class="whs1">aGroup.Add(anIds)</p>
\r
823 <p class="whs1"> </p>
\r
825 <p class="whs1">salome.sg.updateObjBrowser(1)
\r
828 <h3><a name=bookmark8>Aspect Ratio</a></h3>
\r
830 <p class="whs1">import SMESH_mechanic</p>
\r
832 <p class="whs1"> </p>
\r
834 <p class="whs1">smesh =
\r
835 SMESH_mechanic.smesh</p>
\r
837 <p class="whs1">mesh =
\r
838 SMESH_mechanic.mesh</p>
\r
840 <p class="whs1">salome = SMESH_mechanic.salome</p>
\r
842 <p class="whs1"> </p>
\r
844 <p class="whs3"># Criterion : ASPECT
\r
847 <p class="whs1">ar_margin = 1.8</p>
\r
849 <p class="whs1"> </p>
\r
851 <p class="whs1"><span style="font-family: 'Lucida Console', monospace;">aFilter
\r
852 = smesh.GetFilter(smesh.FACE, smesh.FT_AspectRatio, smesh.FT_MoreThan,
\r
853 ar_margin)</span></p>
\r
855 <p class="whs1"> </p>
\r
857 <p class="whs1">anIds = mesh.GetIdsFromFilter(aFilter)
\r
860 <p class="whs1"> </p>
\r
862 <p class="whs3"># print the result</p>
\r
864 <p class="whs1">print "Criterion:
\r
865 Aspect Ratio > ", ar_margin, " Nb = ", len(anIds)</p>
\r
867 <p class="whs1">j = 1</p>
\r
869 <p class="whs1">for i in range(len(anIds)):</p>
\r
871 <p class="whs1"> if
\r
872 j > 20: j = 1; print ""</p>
\r
874 <p class="whs1"> print
\r
877 <p class="whs1"> j
\r
880 <p class="whs1"> pass</p>
\r
882 <p class="whs1">print ""</p>
\r
884 <p class="whs1"> </p>
\r
886 <p class="whs3"># create a group</p>
\r
888 <p class="whs1"><span style="font-family: 'Lucida Console', monospace;">aGroup
\r
889 = mesh.CreateEmptyGroup(smesh.FACE, "Aspect Ratio > " + `ar_margin`)</span></p>
\r
891 <p class="whs1">aGroup.Add(anIds)</p>
\r
893 <p class="whs1"> </p>
\r
895 <p class="whs1">salome.sg.updateObjBrowser(1)
\r
898 <p class="whs1"> </p>
\r
900 <h3><a name=bookmark9>Minimum Angle</a></h3>
\r
902 <p class="whs1">import SMESH_mechanic</p>
\r
904 <p class="whs1"> </p>
\r
906 <p class="whs1">smesh =
\r
907 SMESH_mechanic.smesh</p>
\r
909 <p class="whs1">mesh =
\r
910 SMESH_mechanic.mesh</p>
\r
912 <p class="whs1">salome = SMESH_mechanic.salome</p>
\r
914 <p class="whs1"> </p>
\r
916 <p class="whs3"># Criterion : MINIMUM
\r
919 <p class="whs1">min_angle = 35.</p>
\r
921 <p class="whs1"> </p>
\r
923 <p class="whs1">aFilter = smesh.GetFilter(smesh.FACE,
\r
924 smesh.FT_MinimumAngle, smesh.FT_LessThan, min_angle)</p>
\r
926 <p class="whs1"> </p>
\r
928 <p class="whs1">anIds = mesh.GetIdsFromFilter(aFilter)
\r
931 <p class="whs1"> </p>
\r
933 <p class="whs3"># print the result</p>
\r
935 <p class="whs1">print "Criterion:
\r
936 Minimum Angle < ", min_angle, " Nb = ", len(anIds)</p>
\r
938 <p class="whs1">j = 1</p>
\r
940 <p class="whs1">for i in range(len(anIds)):</p>
\r
942 <p class="whs1"> if
\r
943 j > 20: j = 1; print ""</p>
\r
945 <p class="whs1"> print
\r
948 <p class="whs1"> j
\r
951 <p class="whs1"> pass</p>
\r
953 <p class="whs1">print ""</p>
\r
955 <p class="whs1"> </p>
\r
957 <p class="whs3"># create a group</p>
\r
959 <p class="whs1"><span style="font-family: 'Lucida Console', monospace;">aGroup
\r
960 = mesh.CreateEmptyGroup(smesh.FACE, "Minimum Angle < " +
\r
961 `min_angle`)</span></p>
\r
963 <p class="whs1"> </p>
\r
965 <p class="whs1">aGroup.Add(anIds)</p>
\r
967 <p class="whs1"> </p>
\r
969 <p class="whs1">salome.sg.updateObjBrowser(1)
\r
972 <p class="whs1"> </p>
\r
974 <h3><a name=bookmark10>Warping</a></h3>
\r
976 <p class="whs1">import SMESH_mechanic</p>
\r
978 <p class="whs1"> </p>
\r
980 <p class="whs1">smesh =
\r
981 SMESH_mechanic.smesh</p>
\r
983 <p class="whs1">mesh =
\r
984 SMESH_mechanic.mesh</p>
\r
986 <p class="whs1">salome = SMESH_mechanic.salome</p>
\r
988 <p class="whs1"> </p>
\r
990 <p class="whs3"># Criterion : WARP
\r
991 ANGLE > 1e-15</p>
\r
993 <p class="whs1">wa_margin = 1e-15</p>
\r
995 <p class="whs1"> </p>
\r
997 <p class="whs1">aFilter = smesh.GetFilter(smesh.FACE,
\r
998 smesh.FT_Warping, smesh.FT_MoreThan, wa_margin)</p>
\r
1000 <p class="whs1"> </p>
\r
1002 <p class="whs1">anIds = mesh.GetIdsFromFilter(aFilter)
\r
1005 <p class="whs1"> </p>
\r
1007 <p class="whs3"># print the result</p>
\r
1009 <p class="whs1">print "Criterion:
\r
1010 Warp > ", wa_margin, " Nb = ", len(anIds)</p>
\r
1012 <p class="whs1">j = 1</p>
\r
1014 <p class="whs1">for i in range(len(anIds)):</p>
\r
1016 <p class="whs1"> if
\r
1017 j > 20: j = 1; print ""</p>
\r
1019 <p class="whs1"> print
\r
1022 <p class="whs1"> j
\r
1025 <p class="whs1"> pass</p>
\r
1027 <p class="whs1">print ""</p>
\r
1029 <p class="whs1"> </p>
\r
1031 <p class="whs3"># create a group</p>
\r
1033 <p class="whs1">aGroup = mesh.CreateEmptyGroup(smesh.FACE,
\r
1034 "Warp > " + `wa_margin`)</p>
\r
1036 <p class="whs1"> </p>
\r
1038 <p class="whs1">aGroup.Add(anIds)</p>
\r
1040 <p class="whs1"> </p>
\r
1042 <p class="whs1">salome.sg.updateObjBrowser(1)
\r
1045 <p class="whs1"> </p>
\r
1047 <h3><a name=bookmark12>Skew</a></h3>
\r
1049 <p class="whs1">import SMESH_mechanic</p>
\r
1051 <p class="whs1"> </p>
\r
1053 <p class="whs1">smesh =
\r
1054 SMESH_mechanic.smesh</p>
\r
1056 <p class="whs1">mesh =
\r
1057 SMESH_mechanic.mesh</p>
\r
1059 <p class="whs1">salome = SMESH_mechanic.salome</p>
\r
1061 <p class="whs1"> </p>
\r
1063 <p class="whs3"># Criterion : Skew
\r
1066 <p class="whs1">skew_margin = 38.</p>
\r
1068 <p class="whs1"> </p>
\r
1070 <p class="whs1">aFilter = smesh.GetFilter(smesh.FACE,
\r
1071 smesh.FT_Skew, smesh.FT_MoreThan, skew_margin)</p>
\r
1073 <p class="whs1"> </p>
\r
1075 <p class="whs1">anIds = mesh.GetIdsFromFilter(aFilter)</p>
\r
1077 <p class="whs1"> </p>
\r
1079 <p class="whs3"># print the result</p>
\r
1081 <p class="whs1">print "Criterion:
\r
1082 Skew > ", skew_margin, " Nb = ", len(anIds)</p>
\r
1084 <p class="whs1">j = 1</p>
\r
1086 <p class="whs1">for i in range(len(anIds)):</p>
\r
1088 <p class="whs1"> if
\r
1089 j > 20: j = 1; print ""</p>
\r
1091 <p class="whs1"> print
\r
1094 <p class="whs1"> j
\r
1097 <p class="whs1"> pass</p>
\r
1099 <p class="whs1">print ""</p>
\r
1101 <p class="whs1"> </p>
\r
1103 <p class="whs3"># create a group</p>
\r
1105 <p class="whs1"><span style="font-family: 'Lucida Console', monospace;">aGroup
\r
1106 = mesh.CreateEmptyGroup(smesh.FACE, "Skew > " + `skew_margin`)</span></p>
\r
1108 <p class="whs1">aGroup.Add(anIds)</p>
\r
1110 <p class="whs1"> </p>
\r
1112 <p class="whs1">salome.sg.updateObjBrowser(1)
\r
1115 <p class="whs1"> </p>
\r
1117 <h3><a name=bookmark11>Aspect Ratio 3D</a></h3>
\r
1119 <p class="whs1">import SMESH_mechanic_tetra
\r
1122 <p class="whs1"> </p>
\r
1124 <p class="whs1">smesh =
\r
1125 SMESH_mechanic_tetra.smesh</p>
\r
1127 <p class="whs1">mesh =
\r
1128 SMESH_mechanic_tetra.mesh</p>
\r
1130 <p class="whs1">salome = SMESH_mechanic_tetra.salome</p>
\r
1132 <p class="whs1"> </p>
\r
1134 <p class="whs3"># Criterion : ASPECT
\r
1135 RATIO 3D > 4.5</p>
\r
1137 <p class="whs1">ar_margin = 4.5</p>
\r
1139 <p class="whs1"> </p>
\r
1141 <p class="whs1">aFilter = smesh.GetFilter(smesh.VOLUME,
\r
1142 smesh.FT_AspectRatio3D, smesh.FT_MoreThan, ar_margin)</p>
\r
1144 <p class="whs1"> </p>
\r
1146 <p class="whs1">anIds = mesh.GetIdsFromFilter(aFilter)
\r
1149 <p class="whs1"> </p>
\r
1151 <p class="whs3"># print the result</p>
\r
1153 <p class="whs1">print "Criterion:
\r
1154 Aspect Ratio 3D > ", ar_margin, " Nb = ", len(anIds)</p>
\r
1156 <p class="whs1">j = 1</p>
\r
1158 <p class="whs1">for i in range(len(anIds)):</p>
\r
1160 <p class="whs1"> if
\r
1161 j > 20: j = 1; print ""</p>
\r
1163 <p class="whs1"> print
\r
1166 <p class="whs1"> j
\r
1169 <p class="whs1"> pass</p>
\r
1171 <p class="whs1">print ""</p>
\r
1173 <p class="whs1"> </p>
\r
1175 <p class="whs3"># create a group</p>
\r
1177 <p class="whs1">aGroup = mesh.CreateEmptyGroup(smesh.VOLUME,
\r
1178 "Aspect Ratio 3D > " + `ar_margin`)</p>
\r
1180 <p class="whs1"> </p>
\r
1182 <p class="whs1">aGroup.Add(anIds)</p>
\r
1184 <p class="whs1"> </p>
\r
1186 <p class="whs1">salome.sg.updateObjBrowser(1)
\r
1189 <h3><a name=Volume>Volume</a></h3>
\r
1191 <p class="whs1">import SMESH_mechanic_tetra</p>
\r
1193 <p class="whs1"> </p>
\r
1195 <p class="whs1">smesh =
\r
1196 SMESH_mechanic_tetra.smesh</p>
\r
1198 <p class="whs1">mesh =
\r
1199 SMESH_mechanic_tetra.mesh</p>
\r
1201 <p class="whs1">salome = SMESH_mechanic_tetra.salome</p>
\r
1203 <p class="whs2"> </p>
\r
1205 <p class="whs2"># Criterion : VOLUME <
\r
1208 <p class="whs1">volume_margin =
\r
1211 <p class="whs1"> </p>
\r
1213 <p class="whs1">aFilter = smesh.GetFilter(smesh.VOLUME,
\r
1214 smesh.FT_Volume3D, smesh.FT_LessThan, volume_margin)</p>
\r
1216 <p class="whs1"> </p>
\r
1218 <p class="whs1">anIds = mesh.GetIdsFromFilter(aFilter)
\r
1221 <p class="whs2"> </p>
\r
1223 <p class="whs2"># print the result</p>
\r
1225 <p class="whs1">print ""</p>
\r
1227 <p class="whs1">print "Criterion:
\r
1228 Volume < ", volume_margin, " Nb = ", len(anIds)</p>
\r
1230 <p class="whs1">j = 1</p>
\r
1232 <p class="whs1">for i in range(len(anIds)):</p>
\r
1234 <p class="whs1"> if
\r
1235 j > 20: j = 1; print ""</p>
\r
1237 <p class="whs1"> print
\r
1240 <p class="whs1"> j
\r
1243 <p class="whs1"> pass</p>
\r
1245 <p class="whs1">print ""</p>
\r
1247 <p class="whs2"> </p>
\r
1249 <p class="whs2"># create a group</p>
\r
1251 <p class="whs1">aGroup = mesh.CreateEmptyGroup(smesh.VOLUME,
\r
1252 "Volume < " + `volume_margin`)</p>
\r
1254 <p class="whs1"> </p>
\r
1256 <p class="whs1">aGroup.Add(anIds)</p>
\r
1258 <p class="whs2"> </p>
\r
1260 <p class="whs2"><span style="font-family: 'Lucida Console', monospace;">salome.sg.updateObjBrowser(1)
\r
1263 <script type="text/javascript" language="javascript1.2">
\r
1265 if (window.writeIntopicBar)
\r
1266 writeIntopicBar(0);
\r