1 <!doctype HTML public "-//W3C//DTD HTML 4.0 Frameset//EN">
\r
6 <title>Transforming Meshes</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 h3.whs1 { margin-top:0pt; margin-bottom:0pt; }
\r
11 h4.whs2 { margin-top:0pt; margin-bottom:0pt; }
\r
12 p.whs3 { font-family:'Lucida Console' , monospace; margin-top:0pt; margin-bottom:0pt; }
\r
13 p.whs4 { margin-top:0pt; margin-bottom:0pt; }
\r
14 p.whs5 { margin-top:0pt; margin-bottom:0pt; font-family:'Lucida Console' , monospace; }
\r
16 </style><script type="text/javascript" language="JavaScript">
\r
18 if ((navigator.appName == "Netscape") && (parseInt(navigator.appVersion) == 4))
\r
20 var strNSS = "<style type='text/css'>";
\r
21 strNSS += "h3.whs1 {margin-top:1pt;margin-bottom:1pt; }";
\r
22 strNSS += "h4.whs2 {margin-top:1pt;margin-bottom:1pt; }";
\r
23 strNSS += "p.whs3 {margin-top:1pt;margin-bottom:1pt; }";
\r
24 strNSS += "p.whs4 {margin-top:1pt;margin-bottom:1pt; }";
\r
25 strNSS += "p.whs5 {margin-top:1pt;margin-bottom:1pt; }";
\r
26 strNSS +="</style>";
\r
27 document.write(strNSS);
\r
31 <script type="text/javascript" language="JavaScript" title="WebHelpInlineScript">
\r
34 if (innerWidth != origWidth || innerHeight != origHeight)
\r
37 if ((parseInt(navigator.appVersion) == 4) && (navigator.appName == "Netscape")) {
\r
38 origWidth = innerWidth;
\r
39 origHeight = innerHeight;
\r
45 <style type="text/css">
\r
47 div.WebHelpPopupMenu { position:absolute; left:0px; top:0px; z-index:4; visibility:hidden; }
\r
49 </style><script type="text/javascript" language="javascript1.2" src="whmsg.js"></script>
\r
50 <script type="text/javascript" language="javascript" src="whver.js"></script>
\r
51 <script type="text/javascript" language="javascript1.2" src="whproxy.js"></script>
\r
52 <script type="text/javascript" language="javascript1.2" src="whutils.js"></script>
\r
53 <script type="text/javascript" language="javascript1.2" src="whtopic.js"></script>
\r
55 <body><script type="text/javascript" language="javascript1.2">
\r
57 if (window.gbWhTopic)
\r
59 if (window.addTocInfo)
\r
61 addTocInfo("MESH module\nTUI Scripts\nTransforming Meshes");
\r
62 addButton("show",BTN_TEXT,"Show","","","","",0,0,"","","");
\r
65 if (window.writeBtnStyle)
\r
68 if (window.writeIntopicBar)
\r
71 if (window.setRelStartPage)
\r
73 setRelStartPage("smesh.htm");
\r
82 document.location.reload();
\r
85 <h1>Transforming Meshes</h1>
\r
87 <h3 class="whs1">Transforming Meshes</h3>
\r
89 <h4 class="whs2"> </h4>
\r
91 <h4 class="whs2"><a name=bookmark>Translation</a></h4>
\r
93 <p class="whs3"> </p>
\r
95 <p class="whs3">import SMESH</p>
\r
97 <p class="whs3">import SMESH_mechanic</p>
\r
99 <p class="whs3"> </p>
\r
101 <p class="whs3">mesh =
\r
102 SMESH_mechanic.mesh</p>
\r
104 <p class="whs3">salome = SMESH_mechanic.salome</p>
\r
106 <p class="whs4"> </p>
\r
108 <p class="whs4"># translate a mesh</p>
\r
110 <p class="whs3">anEditor = mesh.GetMeshEditor()</p>
\r
112 <p class="whs3">point = SMESH.PointStruct(-150, -150,
\r
115 <p class="whs3">vector = SMESH.DirStruct(point) </p>
\r
117 <p class="whs3">anEditor.TranslateObject(mesh, vector,
\r
120 <p class="whs4"> </p>
\r
122 <p class="whs4"><span style="font-family: 'Lucida Console', monospace;">salome.sg.updateObjBrowser(1)</span>
\r
125 <p class="whs4"> </p>
\r
127 <h4 class="whs2"><a name=bookmark1>Rotation</a></h4>
\r
129 <p class="whs3"> </p>
\r
131 <p class="whs3">import SMESH</p>
\r
133 <p class="whs3">import SMESH_mechanic</p>
\r
135 <p class="whs3"> </p>
\r
137 <p class="whs3">import math</p>
\r
139 <p class="whs3">mesh =
\r
140 SMESH_mechanic.mesh</p>
\r
142 <p class="whs3">salome = SMESH_mechanic.salome</p>
\r
144 <p class="whs4"> </p>
\r
146 <p class="whs4"># rotate a mesh </p>
\r
148 <p class="whs3">anEditor = mesh.GetMeshEditor()</p>
\r
150 <p class="whs3">axisXYZ = SMESH.AxisStruct(0, 0, 0, 5,
\r
153 <p class="whs3">angle180 = 1.5*math.pi</p>
\r
155 <p class="whs3">anEditor.RotateObject(mesh, axisXYZ, angle180,
\r
158 <p class="whs3"> </p>
\r
160 <p class="whs5">salome.sg.updateObjBrowser(1)</p>
\r
162 <p class="whs4"> </p>
\r
164 <p class="whs4"> </p>
\r
166 <h4 class="whs2"><a name=bookmark2>Symmetry</a></h4>
\r
168 <p class="whs3"> </p>
\r
170 <p class="whs3">import SMESH</p>
\r
172 <p class="whs3">import SMESH_mechanic</p>
\r
174 <p class="whs3"> </p>
\r
176 <p class="whs3">import math</p>
\r
178 <p class="whs3">mesh =
\r
179 SMESH_mechanic.mesh</p>
\r
181 <p class="whs3">salome = SMESH_mechanic.salome</p>
\r
183 <p class="whs4"> </p>
\r
185 <p class="whs4"># create a symmetrical
\r
186 copy of the mesh mirrored through a point</p>
\r
188 <p class="whs3">anEditor = mesh.GetMeshEditor()</p>
\r
190 <p class="whs3">anEditor.MirrorObject(mesh, SMESH.AxisStruct(0,
\r
191 0, 0, 0, 0, 0), SMESH.SMESH_MeshEditor.POINT, 1) </p>
\r
193 <p class="whs3"> </p>
\r
195 <p class="whs4"><span style="font-family: 'Lucida Console', monospace;">salome.sg.updateObjBrowser(1)</span>
\r
198 <p class="whs4"> </p>
\r
200 <p class="whs4"> </p>
\r
202 <h4 class="whs2"><a name=bookmark3>Merging
\r
205 <p class="whs3"> </p>
\r
207 <p class="whs3">import SMESH</p>
\r
209 <p class="whs3">import SMESH_mechanic</p>
\r
211 <p class="whs3"> </p>
\r
213 <p class="whs3">smesh =
\r
214 SMESH_mechanic.smesh</p>
\r
216 <p class="whs3">mesh =
\r
217 SMESH_mechanic.mesh</p>
\r
219 <p class="whs3">salome = SMESH_mechanic.salome</p>
\r
221 <p class="whs4"> </p>
\r
223 <p class="whs4"># merge nodes</p>
\r
225 <p class="whs3">aMeshEditor = mesh.GetMeshEditor()</p>
\r
227 <p class="whs3">Tolerance = 25.0</p>
\r
229 <p class="whs3"> </p>
\r
231 <p class="whs3">GroupsOfNodes = aMeshEditor.FindCoincidentNodes(Tolerance)</p>
\r
233 <p class="whs3">aMeshEditor.MergeNodes(GroupsOfNodes)</p>
\r
235 <p class="whs3"> </p>
\r
237 <p class="whs4"><span style="font-family: 'Lucida Console', monospace;">salome.sg.updateObjBrowser(1)</span>
\r
240 <p class="whs4"> </p>
\r
242 <h4><a name=bookmark8>Merging Elements</a></h4>
\r
244 <p class="whs5">import geompy</p>
\r
246 <p class="whs5">import smesh</p>
\r
248 <p class="whs5">import salome</p>
\r
250 <p class="whs5">import SMESH</p>
\r
252 <p class="whs4"> </p>
\r
254 <p class="whs4"># create a face to be meshed</p>
\r
256 <p class="whs5">px = geompy.MakeVertex(100.,
\r
260 <p class="whs5">py = geompy.MakeVertex(0.
\r
264 <p class="whs5">pz = geompy.MakeVertex(0.
\r
265 , 0. ,
\r
268 <p class="whs5"> </p>
\r
270 <p class="whs5">vxy = geompy.MakeVector(px,
\r
273 <p class="whs5">arc = geompy.MakeArc(py,
\r
276 <p class="whs5"> </p>
\r
278 <p class="whs5">wire = geompy.MakeWire([vxy,
\r
281 <p class="whs5">isPlanarFace = 1</p>
\r
283 <p class="whs5"> </p>
\r
285 <p class="whs5">face1 = geompy.MakeFace(wire,
\r
288 <p class="whs5">id_face1 = geompy.addToStudy(face1,
\r
289 "Face1")</p>
\r
291 <p class="whs4"> </p>
\r
293 <p class="whs4"># create a 2D mesh on the
\r
296 <p class="whs5">trias = smesh.Mesh(face1,
\r
297 "Face : 2D mesh")</p>
\r
299 <p class="whs5"> </p>
\r
301 <p class="whs5">algo = trias.Segment()</p>
\r
303 <p class="whs5">algo.NumberOfSegments(6)</p>
\r
305 <p class="whs5"> </p>
\r
307 <p class="whs5">algo = trias.Triangle()</p>
\r
309 <p class="whs5">algo.LengthFromEdges()</p>
\r
311 <p class="whs5"> </p>
\r
313 <p class="whs5">trias.Compute()</p>
\r
315 <p class="whs5"> </p>
\r
317 <p class="whs5">tri_mesh = trias.GetMesh()</p>
\r
319 <p class="whs4"> </p>
\r
321 <p class="whs4"># create a path mesh</p>
\r
323 <p class="whs5">px1 = geompy.MakeVertex(100.,
\r
327 <p class="whs5">py1 = geompy.MakeVertex(-100.
\r
331 <p class="whs5">pz1 = geompy.MakeVertex(0.
\r
332 , 0. ,
\r
335 <p class="whs5"> </p>
\r
337 <p class="whs5">circle = geompy.MakeCircleThreePnt(py1,
\r
340 <p class="whs5">id_circle = geompy.addToStudy(circle,
\r
341 "Path")</p>
\r
343 <p class="whs5"> </p>
\r
345 <p class="whs5">circlemesh = smesh.Mesh(circle,
\r
346 "Path mesh")</p>
\r
348 <p class="whs5"> </p>
\r
350 <p class="whs5">algo = circlemesh.Segment()</p>
\r
352 <p class="whs5">algo.NumberOfSegments(10)</p>
\r
354 <p class="whs5"> </p>
\r
356 <p class="whs5">circlemesh.Compute()</p>
\r
358 <p class="whs4"> </p>
\r
360 <p class="whs4"># extrusion of the mesh</p>
\r
362 <p class="whs5">aMeshEditor = tri_mesh.GetMeshEditor()</p>
\r
364 <p class="whs5">aMeshEditor.ExtrusionAlongPathObject(tri_mesh,
\r
365 circlemesh.GetMesh(), circle, 1, 0, [], 0, SMESH.PointStruct(0, 0, 0))</p>
\r
367 <p class="whs4"> </p>
\r
369 <p class="whs4"># merge nodes</p>
\r
371 <p class="whs5">print "Number
\r
372 of nodes before MergeNodes:", tri_mesh.NbNodes()</p>
\r
374 <p class="whs5"> </p>
\r
376 <p class="whs5">tolerance = 0.001</p>
\r
378 <p class="whs5">array_of_nodes_groups
\r
379 = aMeshEditor.FindCoincidentNodes(tolerance)</p>
\r
381 <p class="whs5">aMeshEditor.MergeNodes(array_of_nodes_groups)</p>
\r
383 <p class="whs5"> </p>
\r
385 <p class="whs5">print "Number
\r
386 of nodes after MergeNodes:", tri_mesh.NbNodes()</p>
\r
388 <p class="whs5">print ""</p>
\r
390 <p class="whs5">print "Number
\r
391 of elements before MergeEqualElements:"</p>
\r
393 <p class="whs5">print "Edges
\r
394 :
\r
395 ", tri_mesh.NbEdges()</p>
\r
397 <p class="whs5">print "Triangles
\r
399 tri_mesh.NbTriangles()</p>
\r
401 <p class="whs5">print "Quadrangles:
\r
402 ", tri_mesh.NbQuadrangles()</p>
\r
404 <p class="whs5">print "Volumes
\r
405 :
\r
406 ", tri_mesh.NbVolumes()</p>
\r
408 <p class="whs4"> </p>
\r
410 <p class="whs4"># merge elements</p>
\r
412 <p class="whs5">aMeshEditor.MergeEqualElements()</p>
\r
414 <p class="whs5"> </p>
\r
416 <p class="whs5">print "Number
\r
417 of elements after MergeEqualElements:"</p>
\r
419 <p class="whs5">print "Edges
\r
420 :
\r
421 ", tri_mesh.NbEdges()</p>
\r
423 <p class="whs5">print "Triangles
\r
425 tri_mesh.NbTriangles()</p>
\r
427 <p class="whs5">print "Quadrangles:
\r
428 ", tri_mesh.NbQuadrangles()</p>
\r
430 <p class="whs5">print "Volumes
\r
431 :
\r
432 ", tri_mesh.NbVolumes()</p>
\r
434 <p class="whs5"> </p>
\r
436 <p class="whs4"><span style="font-family: 'Lucida Console', monospace;">salome.sg.updateObjBrowser(1)</span>
\r
439 <h3>Sewing Meshes</h3>
\r
441 <h4><a name=bookmark4>Sew Meshes Border to Border</a></h4>
\r
443 <p class="whs5">import salome</p>
\r
445 <p class="whs5">import geompy</p>
\r
447 <p class="whs5">import smesh</p>
\r
449 <p class="whs4"> </p>
\r
451 <p class="whs4"># create two faces of a
\r
454 <p class="whs5">box1 = geompy.MakeBox(0.,
\r
455 0., -10., 30., 20., 25.)</p>
\r
457 <p class="whs5">subShapeList1 =
\r
458 geompy.SubShapeAll(box1, geompy.ShapeType["FACE"])</p>
\r
460 <p class="whs5"> </p>
\r
462 <p class="whs5">box2 = geompy.MakeBox(0.,
\r
463 5., 0., 20., 20., 15.)</p>
\r
465 <p class="whs5">subShapeList2 =
\r
466 geompy.SubShapeAll(box2, geompy.ShapeType["FACE"])</p>
\r
468 <p class="whs5">EdgesList = geompy.SubShapeAll(subShapeList2[
\r
469 1 ], geompy.ShapeType["EDGE"])</p>
\r
471 <p class="whs5"> </p>
\r
473 <p class="whs5">aComp = geompy.MakeCompound(
\r
474 [subShapeList1[ 2 ], subShapeList2[ 1 ]] )</p>
\r
476 <p class="whs5">idComp = geompy.addToStudy(
\r
477 aComp, "Two faces" )</p>
\r
479 <p class="whs5">aCompobj =
\r
480 salome.IDToObject( idComp )</p>
\r
482 <p class="whs5"> </p>
\r
484 <p class="whs4"># create a mesh on two
\r
487 <p class="whs4"> </p>
\r
489 <p class="whs5">mesh = smesh.Mesh(aCompobj,
\r
490 "Two faces : quadrangle mesh")</p>
\r
492 <p class="whs5">algo = mesh.Segment()</p>
\r
494 <p class="whs5"> </p>
\r
496 <p class="whs4"># define "NumberOfSegments"
\r
497 hypothesis to cut an edge in a fixed number of segments</p>
\r
499 <p class="whs5">algo.NumberOfSegments(9)</p>
\r
501 <p class="whs4"> </p>
\r
503 <p class="whs4"># create a quadrangle 2D
\r
504 algorithm for faces</p>
\r
506 <p class="whs5">mesh.Quadrangle()</p>
\r
508 <p class="whs4"> </p>
\r
510 <p class="whs4"># create a local hypothesis</p>
\r
512 <p class="whs5">algo = mesh.Segment(EdgesList[2])</p>
\r
514 <p class="whs4"> </p>
\r
516 <p class="whs4"># define "Arithmetic1D"
\r
517 hypothesis to cut an edge in several segments with increasing arithmetic
\r
520 <p class="whs5">algo.Arithmetic1D(1,
\r
523 <p class="whs4"> </p>
\r
525 <p class="whs4"># define "Propagation"
\r
526 hypothesis that propagates all other hypothesis on all edges on the opposite
\r
527 side in case of quadrangular faces</p>
\r
529 <p class="whs5">algo.Propagation()</p>
\r
531 <p class="whs5">mesh.Compute()</p>
\r
533 <p class="whs4"> </p>
\r
535 <p class="whs4"># sew free borders</p>
\r
537 <p class="whs5">anEditor = mesh.GetMesh().GetMeshEditor()</p>
\r
539 <p class="whs5">anEditor.SewBorderToSide(5,
\r
540 45, 6, 113, 109, 0, 0)</p>
\r
542 <p class="whs5"> </p>
\r
544 <p class="whs4"><span style="font-family: 'Lucida Console', monospace;">salome.sg.updateObjBrowser(1)</span>
\r
547 <p class="whs4"> </p>
\r
549 <h4><a name=bookmark5>Sew Conform Free Borders</a></h4>
\r
551 <p class="whs5">import salome</p>
\r
553 <p class="whs5">import geompy</p>
\r
555 <p class="whs5">import smesh</p>
\r
557 <p class="whs4"> </p>
\r
559 <p class="whs4"># create two faces of the
\r
562 <p class="whs5">box1 = geompy.MakeBox(0.,
\r
563 0., -10., 20., 20., 15.)</p>
\r
565 <p class="whs5">subShapeList1 =
\r
566 geompy.SubShapeAll(box1, geompy.ShapeType["FACE"])</p>
\r
568 <p class="whs5"> </p>
\r
570 <p class="whs5">box2 = geompy.MakeBox(0.,
\r
571 5., 0., 20., 20., 15.)</p>
\r
573 <p class="whs5">subShapeList2 =
\r
574 geompy.SubShapeAll(box2, geompy.ShapeType["FACE"])</p>
\r
576 <p class="whs5">EdgesList = geompy.SubShapeAll(subShapeList2[
\r
577 1 ], geompy.ShapeType["EDGE"])</p>
\r
579 <p class="whs5"> </p>
\r
581 <p class="whs5">aComp = geompy.MakeCompound(
\r
582 [subShapeList1[ 2 ], subShapeList2[ 1 ]] )</p>
\r
584 <p class="whs5">idComp = geompy.addToStudy(
\r
585 aComp, "Two faces" )</p>
\r
587 <p class="whs5">aCompobj =
\r
588 salome.IDToObject( idComp )</p>
\r
590 <p class="whs4"> </p>
\r
592 <p class="whs4"># create a mesh on two
\r
595 <p class="whs5">mesh = smesh.Mesh(aCompobj,
\r
596 "Two faces : quadrangle mesh")</p>
\r
598 <p class="whs5">algo = mesh.Segment()</p>
\r
600 <p class="whs4"> </p>
\r
602 <p class="whs4"># define "NumberOfSegments"
\r
603 hypothesis to cut an edge in a fixed number of segments</p>
\r
605 <p class="whs5">algo.NumberOfSegments(9)</p>
\r
607 <p class="whs4"> </p>
\r
609 <p class="whs4"># create a quadrangle 2D
\r
610 algorithm for faces</p>
\r
612 <p class="whs5">mesh.Quadrangle()</p>
\r
614 <p class="whs4"> </p>
\r
616 <p class="whs4"># create a local hypothesis</p>
\r
618 <p class="whs5">algo = mesh.Segment(EdgesList[2])</p>
\r
620 <p class="whs4"> </p>
\r
622 <p class="whs4"># define "Arithmetic1D"
\r
623 hypothesis to cut an edge in several segments with increasing arithmetic
\r
626 <p class="whs5">algo.Arithmetic1D(1,
\r
629 <p class="whs4"> </p>
\r
631 <p class="whs4"># define "Propagation"
\r
632 hypothesis that propagate all other hypothesis on all edges on<span style="margin-top: 0pt; margin-bottom: 0pt;">
\r
633 the opposite side in case of quadrangular faces</span></p>
\r
635 <p class="whs5">algo.Propagation()</p>
\r
637 <p class="whs5">mesh.Compute()</p>
\r
639 <p class="whs4"> </p>
\r
641 <p class="whs4"># sew free borders</p>
\r
643 <p class="whs5">anEditor = mesh.GetMesh().GetMeshEditor()</p>
\r
645 <p class="whs5">anEditor.SewConformFreeBorders(5,
\r
648 <p class="whs4"><span style="font-family: 'Lucida Console', monospace;">salome.sg.updateObjBrowser(1)</span>
\r
651 <p class="whs4"> </p>
\r
653 <h4><a name=bookmark6>Sew Free Borders</a></h4>
\r
655 <p class="whs5">import salome</p>
\r
657 <p class="whs5">import geompy</p>
\r
659 <p class="whs5">import smesh</p>
\r
661 <p class="whs4"> </p>
\r
663 <p class="whs4"># create two faces of the
\r
666 <p class="whs5">box1 = geompy.MakeBox(0.,
\r
667 0., 0., 20., 20., 15.)</p>
\r
669 <p class="whs5">subShapeList1 =
\r
670 geompy.SubShapeAll(box1, geompy.ShapeType["FACE"])</p>
\r
672 <p class="whs4"> </p>
\r
674 <p class="whs5">box2 = geompy.MakeBox(0.,
\r
675 5., 0., 20., 20., 15.)</p>
\r
677 <p class="whs5">subShapeList2 =
\r
678 geompy.SubShapeAll(box2, geompy.ShapeType["FACE"])</p>
\r
680 <p class="whs5">EdgesList = geompy.SubShapeAll(subShapeList2[
\r
681 1 ], geompy.ShapeType["EDGE"])</p>
\r
683 <p class="whs5"> </p>
\r
685 <p class="whs5">aComp = geompy.MakeCompound(
\r
686 [subShapeList1[ 2 ], subShapeList2[ 1 ]] )</p>
\r
688 <p class="whs5">idComp = geompy.addToStudy(
\r
689 aComp, "Two faces" )</p>
\r
691 <p class="whs5">aCompobj =
\r
692 salome.IDToObject( idComp )</p>
\r
694 <p class="whs5"> </p>
\r
696 <p class="whs4"># create a mesh on two
\r
699 <p class="whs5">mesh = smesh.Mesh(aCompobj,
\r
700 "Two faces : quadrangle mesh")</p>
\r
702 <p class="whs5">algo = mesh.Segment()</p>
\r
704 <p class="whs4"> </p>
\r
706 <p class="whs4"># define "NumberOfSegments"
\r
707 hypothesis to cut an edge in a fixed number of segments</p>
\r
709 <p class="whs5">algo.NumberOfSegments(4)</p>
\r
711 <p class="whs4"> </p>
\r
713 <p class="whs4"># creates a quadrangle
\r
714 2D algorithm for faces</p>
\r
716 <p class="whs5">mesh.Quadrangle()</p>
\r
718 <p class="whs4"> </p>
\r
720 <p class="whs4"># create a local hypothesis</p>
\r
722 <p class="whs5">algo = mesh.Segment(EdgesList[2])</p>
\r
724 <p class="whs4"> </p>
\r
726 <p class="whs4"># define "Arithmetic1D"
\r
727 hypothesis to cut an edge in several segments with increasing
\r
728 arithmetic length</p>
\r
730 <p class="whs5">algo.Arithmetic1D(1,
\r
733 <p class="whs4"> </p>
\r
735 <p class="whs4"># define "Propagation"
\r
736 hypothesis that propagate all other hypothesis on all edges on<span style="margin-top: 0pt; margin-bottom: 0pt;">
\r
737 the opposite side in case of quadrangular faces</span></p>
\r
739 <p class="whs5">algo.Propagation()</p>
\r
741 <p class="whs5">mesh.Compute()</p>
\r
743 <p class="whs4"> </p>
\r
745 <p class="whs4"># sew free borders</p>
\r
747 <p class="whs5">anEditor = mesh.GetMesh().GetMeshEditor()</p>
\r
749 <p class="whs5">anEditor.SewFreeBorders(6,
\r
750 21, 5, 1, 12, 3, 0, 0)</p>
\r
752 <p class="whs5"> </p>
\r
754 <p class="whs4"><span style="font-family: 'Lucida Console', monospace;">salome.sg.updateObjBrowser(1)</span>
\r
757 <p class="whs4"> </p>
\r
759 <h4><a name=bookmark7>Sew Side Elements</a></h4>
\r
761 <p class="whs5">import salome</p>
\r
763 <p class="whs5">import geompy</p>
\r
765 <p class="whs5">import smesh</p>
\r
767 <p class="whs5"> </p>
\r
769 <p class="whs4"># create two faces of the
\r
772 <p class="whs5">box1 = geompy.MakeBox(0.,
\r
776 <p class="whs5">box2 = geompy.MakeBox(0.,
\r
777 15., 0., 20., 25., 10.)</p>
\r
779 <p class="whs5">EdgesList = geompy.SubShapeAll(box2,
\r
780 geompy.ShapeType["EDGE"])</p>
\r
782 <p class="whs4"> </p>
\r
784 <p class="whs5">aComp = geompy.MakeCompound(
\r
787 <p class="whs5">idComp = geompy.addToStudy(
\r
788 aComp, "Two faces" )</p>
\r
790 <p class="whs5">aCompobj =
\r
791 salome.IDToObject( idComp )</p>
\r
793 <p class="whs5"> </p>
\r
795 <p class="whs4"># create a mesh on two
\r
798 <p class="whs5">mesh = smesh.Mesh(aCompobj,
\r
799 "Two faces : quadrangle mesh")</p>
\r
801 <p class="whs5">algo = mesh.Segment()</p>
\r
803 <p class="whs4"> </p>
\r
805 <p class="whs4"># define "NumberOfSegments"
\r
806 hypothesis to cut an edge in a fixed number of segments</p>
\r
808 <p class="whs5">algo.NumberOfSegments(2)</p>
\r
810 <p class="whs4"> </p>
\r
812 <p class="whs4"># create a quadrangle 2D
\r
813 algorithm for faces</p>
\r
815 <p class="whs5">mesh.Quadrangle()</p>
\r
817 <p class="whs4"> </p>
\r
819 <p class="whs4"># create a local hypothesis</p>
\r
821 <p class="whs5">algo = mesh.Segment(EdgesList[8])</p>
\r
823 <p class="whs4"> </p>
\r
825 <p class="whs4"># define "Arithmetic1D"
\r
826 hypothesis to cut an edge in several segments with increasing arithmetic
\r
829 <p class="whs5">algo.NumberOfSegments(4)</p>
\r
831 <p class="whs4"> </p>
\r
833 <p class="whs4"># define "Propagation"
\r
834 hypothesis that propagates all other hypothesis on all edges on <span
\r
835 style="margin-top: 0pt; margin-bottom: 0pt;">the opposite side in case
\r
836 of quadrangular faces</span></p>
\r
838 <p class="whs5">algo.Propagation()</p>
\r
840 <p class="whs5">mesh.Compute()</p>
\r
842 <p class="whs4"> </p>
\r
844 <p class="whs4"># sew free borders</p>
\r
846 <p class="whs5">anEditor = mesh.GetMesh().GetMeshEditor()</p>
\r
848 <p class="whs5">anEditor.SewSideElements([69,
\r
849 70, 71, 72], [91, 92, 89, 90], 8, 38, 23, 58)</p>
\r
851 <p class="whs4"><span style="font-family: 'Lucida Console', monospace;">salome.sg.updateObjBrowser(1)</span>
\r
854 <script type="text/javascript" language="javascript1.2">
\r
856 if (window.writeIntopicBar)
\r
857 writeIntopicBar(0);
\r