-<!doctype HTML public "-//W3C//DTD HTML 4.0 Frameset//EN">\r
-\r
-<html>\r
-\r
-<head>\r
-<title>Modifying Meshes</title>\r
-<meta http-equiv="content-type" content="text/html; charset=windows-1252">\r
-<meta name="generator" content="RoboHelp by eHelp Corporation www.ehelp.com"><style>\r
-<!--\r
-P { margin-top:0pt; margin-bottom:0pt; }\r
-LI.kadov-P { }\r
--->\r
-</style><style type="text/css">\r
-<!--\r
-p.whs1 { margin-top:0pt; margin-bottom:0pt; font-weight:bold; color:#ff0000; }\r
-p.whs2 { margin-top:0pt; margin-bottom:0pt; }\r
-p.whs3 { margin-top:0pt; margin-bottom:0pt; font-family:'Lucida Console' , monospace; }\r
-p.whs4 { margin-top:0pt; margin-bottom:0pt; font-family:'Times New Roman' , serif; }\r
-h4.whs5 { margin-top:0pt; margin-bottom:0pt; }\r
-p.whs6 { margin-top:0.5pt; margin-bottom:0pt; font-family:'Lucida Console' , monospace; }\r
-p.whs7 { margin-top:0.5pt; margin-bottom:0pt; }\r
-p.whs8 { margin-top:0.5pt; margin-bottom:0pt; font-family:'Times New Roman' , serif; }\r
-p.whs9 { font-family:'Lucida Console' , monospace; }\r
-p.whs10 { font-family:'Times New Roman' , serif; }\r
-p.whs11 { margin-left:36px; font-family:'Lucida Console' , monospace; }\r
-p.whs12 { margin-left:36px; }\r
-p.whs13 { margin-left:192px; font-family:'Lucida Console' , monospace; }\r
-p.whs14 { font-family:'Lucida Console' , monospace; margin-left:0px; }\r
--->\r
-</style><script type="text/javascript" language="JavaScript">\r
-<!--\r
-if ((navigator.appName == "Netscape") && (parseInt(navigator.appVersion) == 4))\r
-{\r
- var strNSS = "<style type='text/css'>";\r
- strNSS += "p.whs1 {margin-top:1pt;margin-bottom:1pt; }";\r
- strNSS += "p.whs2 {margin-top:1pt;margin-bottom:1pt; }";\r
- strNSS += "p.whs3 {margin-top:1pt;margin-bottom:1pt; }";\r
- strNSS += "p.whs4 {margin-top:1pt;margin-bottom:1pt; }";\r
- strNSS += "h4.whs5 {margin-top:1pt;margin-bottom:1pt; }";\r
- strNSS += "p.whs6 {margin-top:1pt;margin-bottom:1pt; }";\r
- strNSS += "p.whs7 {margin-top:1pt;margin-bottom:1pt; }";\r
- strNSS += "p.whs8 {margin-top:1pt;margin-bottom:1pt; }";\r
- strNSS += "p.whs14 {margin-left:1pt; }";\r
- strNSS +="</style>";\r
- document.write(strNSS);\r
-}\r
-//-->\r
-</script>\r
-<script type="text/javascript" language="JavaScript" title="WebHelpInlineScript">\r
-<!--\r
-function reDo() {\r
- if (innerWidth != origWidth || innerHeight != origHeight)\r
- location.reload();\r
-}\r
-if ((parseInt(navigator.appVersion) == 4) && (navigator.appName == "Netscape")) {\r
- origWidth = innerWidth;\r
- origHeight = innerHeight;\r
- onresize = reDo;\r
-}\r
-onerror = null; \r
-//-->\r
-</script>\r
-<style type="text/css">\r
-<!--\r
-div.WebHelpPopupMenu { position:absolute; left:0px; top:0px; z-index:4; visibility:hidden; }\r
-p.WebHelpNavBar { text-align:right; }\r
--->\r
-</style><script type="text/javascript" language="javascript1.2" src="whmsg.js"></script>\r
-<script type="text/javascript" language="javascript" src="whver.js"></script>\r
-<script type="text/javascript" language="javascript1.2" src="whproxy.js"></script>\r
-<script type="text/javascript" language="javascript1.2" src="whutils.js"></script>\r
-<script type="text/javascript" language="javascript1.2" src="whtopic.js"></script>\r
-<script type="text/javascript" language="javascript1.2">\r
-<!--\r
-if (window.gbWhTopic)\r
-{\r
- if (window.setRelStartPage)\r
- {\r
- addTocInfo("MESH module\nTUI Scripts\nModifying Meshes");\r
-addButton("show",BTN_IMG,"Show","","","","",0,0,"whd_show0.gif","whd_show2.gif","whd_show1.gif");\r
-addButton("hide",BTN_IMG,"Hide","","","","",0,0,"whd_hide0.gif","whd_hide2.gif","whd_hide1.gif");\r
-\r
- }\r
-\r
-\r
- if (window.setRelStartPage)\r
- {\r
- setRelStartPage("index.htm");\r
-\r
- autoSync(1);\r
- sendSyncInfo();\r
- sendAveInfoOut();\r
- }\r
-\r
-}\r
-else\r
- if (window.gbIE4)\r
- document.location.reload();\r
-//-->\r
-</script>\r
-</head>\r
-<body><script type="text/javascript" language="javascript1.2">\r
-<!--\r
-if (window.writeIntopicBar)\r
- writeIntopicBar(4);\r
-//-->\r
-</script>\r
-<h1>Modifying Meshes</h1>\r
-\r
-<h3><a name=bookmark>Adding Nodes and Elements</a></h3>\r
-\r
-<p class="whs1"># Attention! The scripts for Adding nodes \r
- and Elements have \r
- been written using the old approach basing on direct usage of SMESH idl \r
- interface.</p>\r
-\r
-<p class="whs1"># For the moment smesh package doesn't provide \r
- methods to add nodes and elements. </p>\r
-\r
-<p class="whs2"><span style="font-weight: bold; color: #ff0000;"><font color=red ><B># \r
- In the next SALOME version the scripts will be updated to use only the \r
- commands from smesh package.</B></font></span></p>\r
-\r
-<h4>Add Node</h4>\r
-\r
-<p class="whs3"><span style="font-family: 'Lucida Console', monospace;">import \r
- SMESH</span></p>\r
-\r
-<p class="whs3">import SMESH_mechanic</p>\r
-\r
-<p class="whs3"> </p>\r
-\r
-<p class="whs3">mesh = SMESH_mechanic.mesh</p>\r
-\r
-<p class="whs3">aMeshEditor = mesh.GetMeshEditor()</p>\r
-\r
-<p class="whs3"> </p>\r
-\r
-<p class="whs4"># add node</p>\r
-\r
-<p class="whs3">new_id = aMeshEditor.AddNode(50, \r
- 10, 0)</p>\r
-\r
-<p class="whs3">print ""</p>\r
-\r
-<p class="whs3">if new_id == 0: \r
- print "KO node addition."</p>\r
-\r
-<p class="whs3">else: print \r
- "New Node has been added with ID ", new_id </p>\r
-\r
-<p class="whs3"> </p>\r
-\r
-<h4>Add Edge</h4>\r
-\r
-<p class="whs3"><span style="font-family: 'Lucida Console', monospace;">import \r
- SMESH</span></p>\r
-\r
-<p class="whs3">import SMESH_mechanic</p>\r
-\r
-<p class="whs3"> </p>\r
-\r
-<p class="whs3">mesh = SMESH_mechanic.mesh</p>\r
-\r
-<p class="whs3">aMeshEditor = mesh.GetMeshEditor()</p>\r
-\r
-<p class="whs3">print ""</p>\r
-\r
-<p class="whs3"> </p>\r
-\r
-<p class="whs4"># add node</p>\r
-\r
-<p class="whs3">n1 = aMeshEditor.AddNode(50, \r
- 10, 0)</p>\r
-\r
-<p class="whs3">if n1 == 0: print \r
- "KO node addition."</p>\r
-\r
-<p class="whs3"> </p>\r
-\r
-<p class="whs4"># add edge</p>\r
-\r
-<p class="whs3">e1 = aMeshEditor.AddEdge([n1, \r
- 38])</p>\r
-\r
-<p class="whs3">if e1 == 0: print \r
- "KO edge addition."</p>\r
-\r
-<p class="whs3">else: print \r
- "New Edge has been added with ID ", e1 </p>\r
-\r
-<p class="whs3"> </p>\r
-\r
-<h4>Add Triangle</h4>\r
-\r
-<p class="whs3"><span style="font-family: 'Lucida Console', monospace;">import \r
- SMESH</span></p>\r
-\r
-<p class="whs3">import SMESH_mechanic</p>\r
-\r
-<p class="whs3"> </p>\r
-\r
-<p class="whs3">mesh = SMESH_mechanic.mesh</p>\r
-\r
-<p class="whs3">aMeshEditor = mesh.GetMeshEditor()</p>\r
-\r
-<p class="whs3">print ""</p>\r
-\r
-<p class="whs3"> </p>\r
-\r
-<p class="whs4"># add node</p>\r
-\r
-<p class="whs3">n1 = aMeshEditor.AddNode(50, \r
- 10, 0)</p>\r
-\r
-<p class="whs3">if n1 == 0: print \r
- "KO node addition."</p>\r
-\r
-<p class="whs3"> </p>\r
-\r
-<p class="whs4"># add triangle</p>\r
-\r
-<p class="whs3">t1 = aMeshEditor.AddFace([n1, \r
- 38, 39])</p>\r
-\r
-<p class="whs3">if t1 == 0: print \r
- "KO triangle addition."</p>\r
-\r
-<p class="whs3">else: print \r
- "New Triangle has been added with ID ", t1 </p>\r
-\r
-<h4>Add Quadrangle</h4>\r
-\r
-<p class="whs3">import SMESH</p>\r
-\r
-<p class="whs3">import SMESH_mechanic</p>\r
-\r
-<p class="whs3"> </p>\r
-\r
-<p class="whs3">mesh = SMESH_mechanic.mesh</p>\r
-\r
-<p class="whs3">aMeshEditor = mesh.GetMeshEditor()</p>\r
-\r
-<p class="whs3">print ""</p>\r
-\r
-<p class="whs2"> </p>\r
-\r
-<p class="whs2"># add node</p>\r
-\r
-<p class="whs3">n1 = aMeshEditor.AddNode(50, \r
- 10, 0)</p>\r
-\r
-<p class="whs3">if n1 == 0: print \r
- "KO node addition."</p>\r
-\r
-<p class="whs3"> </p>\r
-\r
-<p class="whs3">n2 = aMeshEditor.AddNode(40, \r
- 20, 0)</p>\r
-\r
-<p class="whs3">if n2 == 0: print \r
- "KO node addition."</p>\r
-\r
-<p class="whs3"> </p>\r
-\r
-<p class="whs2"># add quadrangle</p>\r
-\r
-<p class="whs3">q1 = aMeshEditor.AddFace([n2, \r
- n1, 38, 39])</p>\r
-\r
-<p class="whs3">if q1 == 0: print \r
- "KO quadrangle addition."</p>\r
-\r
-<p class="whs2"><span style="font-family: 'Lucida Console', monospace;">else: \r
- print \r
- "New Quadrangle has been added with ID ", q1</span> </p>\r
-\r
-<h4>Add Tetrahedron</h4>\r
-\r
-<p class="whs3">import SMESH </p>\r
-\r
-<p class="whs3">import SMESH_mechanic</p>\r
-\r
-<p class="whs3"> </p>\r
-\r
-<p class="whs3">mesh = SMESH_mechanic.mesh</p>\r
-\r
-<p class="whs3">aMeshEditor = mesh.GetMeshEditor()</p>\r
-\r
-<p class="whs3">print ""</p>\r
-\r
-<p class="whs3"> </p>\r
-\r
-<p class="whs4"># add node</p>\r
-\r
-<p class="whs3">n1 = aMeshEditor.AddNode(50, \r
- 10, 0)</p>\r
-\r
-<p class="whs3">if n1 == 0: print \r
- "KO node addition."</p>\r
-\r
-<p class="whs3"> </p>\r
-\r
-<p class="whs4"># add tetrahedron</p>\r
-\r
-<p class="whs3">t1 = aMeshEditor.AddVolume([n1, \r
- 38, 39, 246])</p>\r
-\r
-<p class="whs3">if t1 == 0: print \r
- "KO tetrahedron addition."</p>\r
-\r
-<p class="whs2"><span style="font-family: 'Lucida Console', monospace;">else: \r
- print \r
- "New Tetrahedron has been added with ID ", t1</span> </p>\r
-\r
-<h4>Add Hexahedron</h4>\r
-\r
-<p class="whs3"><span style="font-family: 'Lucida Console', monospace;">import \r
- SMESH</span></p>\r
-\r
-<p class="whs3">import SMESH_mechanic</p>\r
-\r
-<p class="whs3"> </p>\r
-\r
-<p class="whs3">mesh = SMESH_mechanic.mesh</p>\r
-\r
-<p class="whs3">aMeshEditor = mesh.GetMeshEditor()</p>\r
-\r
-<p class="whs3">print ""</p>\r
-\r
-<p class="whs3"> </p>\r
-\r
-<p class="whs4"># add nodes</p>\r
-\r
-<p class="whs3">nId1 = aMeshEditor.AddNode(50, \r
- 10, 0)</p>\r
-\r
-<p class="whs3">nId2 = aMeshEditor.AddNode(47, \r
- 12, 0)</p>\r
-\r
-<p class="whs3">nId3 = aMeshEditor.AddNode(50, \r
- 10, 10)</p>\r
-\r
-<p class="whs3">nId4 = aMeshEditor.AddNode(47, \r
- 12, 10)</p>\r
-\r
-<p class="whs3"> </p>\r
-\r
-<p class="whs3">if nId1 == 0 or \r
- nId2 == 0 or nId3 == 0 or nId4 == 0: print "KO node addition."</p>\r
-\r
-<p class="whs3"> </p>\r
-\r
-<p class="whs4"># add hexahedron</p>\r
-\r
-<p class="whs3">vId = aMeshEditor.AddVolume([nId2, \r
- nId1, 38, 39, nId4, nId3, 245, 246])</p>\r
-\r
-<p class="whs3">if vId == 0: print \r
- "KO Hexahedron addition."</p>\r
-\r
-<p class="whs3">else: print \r
- "New Hexahedron has been added with ID ", vId </p>\r
-\r
-<p class="whs3"> </p>\r
-\r
-<h4>Add Polygon</h4>\r
-\r
-<p class="whs3"><span style="font-family: 'Lucida Console', monospace;">import \r
- math</span></p>\r
-\r
-<p class="whs3">import salome</p>\r
-\r
-<p class="whs3">import smesh</p>\r
-\r
-<p class="whs3"> </p>\r
-\r
-<p class="whs4"># create an empty mesh \r
- structure</p>\r
-\r
-<p class="whs3">gen = smesh.smesh</p>\r
-\r
-<p class="whs3">mesh = gen.CreateEmptyMesh()</p>\r
-\r
-<p class="whs3">MeshEditor = mesh.GetMeshEditor()</p>\r
-\r
-<p class="whs3"> </p>\r
-\r
-<p class="whs4"># a method to build \r
- a polygonal mesh element with <nb_vert> angles:</p>\r
-\r
-<p class="whs3">def MakePolygon \r
- (a_mesh, x0, y0, z0, radius, nb_vert):</p>\r
-\r
-<p class="whs3"> al \r
- = 2.0 * math.pi / nb_vert</p>\r
-\r
-<p class="whs3"> node_ids \r
- = []</p>\r
-\r
-<p class="whs3"> </p>\r
-\r
-<p class="whs4"> # \r
- Create nodes for a polygon</p>\r
-\r
-<p class="whs3"> for \r
- ii in range(nb_vert):</p>\r
-\r
-<p class="whs3"> nid \r
- = MeshEditor.AddNode(x0 + radius * math.cos(ii*al),</p>\r
-\r
-<p class="whs3"> y0 \r
- + radius * math.sin(ii*al),</p>\r
-\r
-<p class="whs3"> z0)</p>\r
-\r
-<p class="whs3"> node_ids.append(nid)</p>\r
-\r
-<p class="whs3"> pass</p>\r
-\r
-<p class="whs3"> </p>\r
-\r
-<p class="whs4"> # \r
- Create a polygon</p>\r
-\r
-<p class="whs3"> return \r
- MeshEditor.AddPolygonalFace(node_ids)</p>\r
-\r
-<p class="whs3"> </p>\r
-\r
-<p class="whs4"># Create three polygons</p>\r
-\r
-<p class="whs3">f1 = MakePolygon(mesh, \r
- 0, 0, 0, \r
- 30, 13)</p>\r
-\r
-<p class="whs3">f2 = MakePolygon(mesh, \r
- 0, 0, 10, 21, 9)</p>\r
-\r
-<p class="whs3">f3 = MakePolygon(mesh, \r
- 0, 0, 20, 13, 6)</p>\r
-\r
-<p class="whs3"> </p>\r
-\r
-<p class="whs3">salome.sg.updateObjBrowser(1) \r
- </p>\r
-\r
-<h4>Add Polyhedron</h4>\r
-\r
-<p class="whs3"><span style="font-family: 'Lucida Console', monospace;">import \r
- salome</span></p>\r
-\r
-<p class="whs3">import math</p>\r
-\r
-<p class="whs3">import smesh</p>\r
-\r
-<p class="whs3"> </p>\r
-\r
-<p class="whs4"># create an empty mesh \r
- structure</p>\r
-\r
-<p class="whs3">gen = smesh.smesh</p>\r
-\r
-<p class="whs3">mesh = gen.CreateEmptyMesh()</p>\r
-\r
-<p class="whs3">MeshEditor = mesh.GetMeshEditor()</p>\r
-\r
-<p class="whs3"> </p>\r
-\r
-<p class="whs4"># Create nodes for \r
- 12-hedron with pentagonal faces</p>\r
-\r
-<p class="whs3">al = 2 * math.pi \r
- / 5.0</p>\r
-\r
-<p class="whs3">cosal = math.cos(al)</p>\r
-\r
-<p class="whs3">aa = 13</p>\r
-\r
-<p class="whs3">rr = aa / (2.0 * \r
- math.sin(al/2.0))</p>\r
-\r
-<p class="whs3">dr = 2.0 * rr * \r
- cosal</p>\r
-\r
-<p class="whs3">r1 = rr + dr</p>\r
-\r
-<p class="whs3">dh = rr * math.sqrt(2.0 \r
- * (1.0 - cosal * (1.0 + 2.0 * cosal)))</p>\r
-\r
-<p class="whs3">hh = 2.0 * dh - \r
- dr * (rr*(cosal - 1) + (rr + dr)*(math.cos(al/2) - 1)) / dh</p>\r
-\r
-<p class="whs3"> </p>\r
-\r
-<p class="whs3">dd = [] <span style="font-family: 'Times New Roman', serif;"># \r
- top</span></p>\r
-\r
-<p class="whs3">cc = [] <span style="font-family: 'Times New Roman', serif;"># \r
- below top</span></p>\r
-\r
-<p class="whs3">bb = [] <span style="font-family: 'Times New Roman', serif;"># \r
- above bottom</span></p>\r
-\r
-<p class="whs3">aa = [] <span style="font-family: 'Times New Roman', serif;"># \r
- bottom</span></p>\r
-\r
-<p class="whs3"> </p>\r
-\r
-<p class="whs3">for i in range(5):</p>\r
-\r
-<p class="whs3"> cos_bot \r
- = math.cos(i*al)</p>\r
-\r
-<p class="whs3"> sin_bot \r
- = math.sin(i*al)</p>\r
-\r
-<p class="whs3"> </p>\r
-\r
-<p class="whs3"> <span \r
- style="margin-top: 0pt;\r
- margin-bottom: 0pt;\r
- font-family: 'Lucida Console', monospace;">cos_top \r
- = math.cos(i*al + al/2.0)</span></p>\r
-\r
-<p class="whs3"> sin_top \r
- = math.sin(i*al + al/2.0)</p>\r
-\r
-<p class="whs3"> </p>\r
-\r
-<p class="whs3"> nd \r
- = MeshEditor.AddNode(rr * cos_top, rr * sin_top, hh ) \r
- <span style="font-family: 'Times New Roman', serif;"># top</span></p>\r
-\r
-<p class="whs3"> nc \r
- = MeshEditor.AddNode(r1 * cos_top, r1 * sin_top, hh - dh) <span style="font-family: 'Times New Roman', serif;"># \r
- below top</span></p>\r
-\r
-<p class="whs3"> nb \r
- = MeshEditor.AddNode(r1 * cos_bot, r1 * sin_bot, dh) \r
- <span style="font-family: 'Times New Roman', serif;"># above bottom</span></p>\r
-\r
-<p class="whs3"> na \r
- = MeshEditor.AddNode(rr * cos_bot, rr * sin_bot, 0) \r
- <span style="font-family: 'Times New Roman', serif;"># bottom</span></p>\r
-\r
-<p class="whs3"> dd.append(nd) \r
- <span style="font-family: 'Times New Roman', serif;"># top</span></p>\r
-\r
-<p class="whs3"> cc.append(nc) \r
- <span style="font-family: 'Times New Roman', serif;"># below top</span></p>\r
-\r
-<p class="whs3"> bb.append(nb) \r
- <span style="font-family: 'Times New Roman', serif;"># above bottom</span></p>\r
-\r
-<p class="whs3"> aa.append(na) \r
- <span style="font-family: 'Times New Roman', serif;"># bottom</span></p>\r
-\r
-<p class="whs3"> </p>\r
-\r
-<p class="whs3"> pass</p>\r
-\r
-<p class="whs3"> </p>\r
-\r
-<p class="whs4"># Create a polyhedral \r
- volume (12-hedron with pentagonal faces)</p>\r
-\r
-<p class="whs3">MeshEditor.AddPolyhedralVolume([dd[0], \r
- dd[1], dd[2], dd[3], dd[4], <span style="font-family: 'Times New Roman', serif;"># \r
- top</span></p>\r
-\r
-<p class="whs3"> dd[0], \r
- cc[0], bb[1], cc[1], dd[1], # \r
- -</p>\r
-\r
-<p class="whs3"> dd[1], \r
- cc[1], bb[2], cc[2], dd[2], # \r
- -</p>\r
-\r
-<p class="whs3"> dd[2], \r
- cc[2], bb[3], cc[3], dd[3], <span style="font-family: 'Times New Roman', serif;"># \r
- - below top</span></p>\r
-\r
-<p class="whs3"> dd[3], \r
- cc[3], bb[4], cc[4], dd[4], # \r
- -</p>\r
-\r
-<p class="whs3"> dd[4], \r
- cc[4], bb[0], cc[0], dd[0], # \r
- -</p>\r
-\r
-<p class="whs3"> aa[4], \r
- bb[4], cc[4], bb[0], aa[0], # \r
- .</p>\r
-\r
-<p class="whs3"> aa[3], \r
- bb[3], cc[3], bb[4], aa[4], # \r
- .</p>\r
-\r
-<p class="whs3"> aa[2], \r
- bb[2], cc[2], bb[3], aa[3], <span style="font-family: 'Times New Roman', serif;"># \r
- . above bottom</span></p>\r
-\r
-<p class="whs3"> aa[1], \r
- bb[1], cc[1], bb[2], aa[2], # \r
- .</p>\r
-\r
-<p class="whs3"> aa[0], \r
- bb[0], cc[0], bb[1], aa[1], # \r
- .</p>\r
-\r
-<p class="whs3"> aa[0], \r
- aa[1], aa[2], aa[3], aa[4]], <span style="font-family: 'Times New Roman', serif;"># \r
- bottom</span></p>\r
-\r
-<p class="whs3"> [5,5,5,5,5,5,5,5,5,5,5,5])</p>\r
-\r
-<p class="whs3">salome.sg.updateObjBrowser(1) \r
- </p>\r
-\r
-<p class="whs3"> </p>\r
-\r
-<h3><a name=bookmark1>Removing Nodes and Elements</a></h3>\r
-\r
-<p class="whs1"><span style="font-weight: bold; color: #ff0000;"><font color=red ><B># \r
- Attention! This script has been written using the old approach basing \r
- on direct usage of SMESH idl interface.</B></font></span></p>\r
-\r
-<p class="whs1"># For the moment smesh package doesn't provide \r
- methods to remove nodes and elements. </p>\r
-\r
-<p class="whs1"># In the next SALOME version the scripts will \r
- be updated to use only the commands from smesh package.</p>\r
-\r
-<h4 class="whs5"> </h4>\r
-\r
-<h4 class="whs5">Removing Nodes</h4>\r
-\r
-<p class="whs2"> </p>\r
-\r
-<p class="whs3"><span style="font-family: 'Lucida Console', monospace;">import \r
- SMESH_mechanic</span></p>\r
-\r
-<p class="whs3"> </p>\r
-\r
-<p class="whs3">mesh = SMESH_mechanic.mesh</p>\r
-\r
-<p class="whs3">aMeshEditor = mesh.GetMeshEditor()</p>\r
-\r
-<p class="whs3"> </p>\r
-\r
-<p class="whs4"># remove nodes #246 \r
- and #255</p>\r
-\r
-<p class="whs3">res = aMeshEditor.RemoveNodes([246, \r
- 255])</p>\r
-\r
-<p class="whs3">if res == 1: print \r
- "Nodes removing is OK!"</p>\r
-\r
-<p class="whs3">else: print \r
- "KO nodes removing." </p>\r
-\r
-<p class="whs2"> </p>\r
-\r
-<h4>Removing Elements</h4>\r
-\r
-<p class="whs3"><span style="font-family: 'Lucida Console', monospace;">import \r
- SMESH_mechanic</span></p>\r
-\r
-<p class="whs3">mesh = SMESH_mechanic.mesh</p>\r
-\r
-<p class="whs3"> </p>\r
-\r
-<p class="whs3">anEditor = mesh.GetMeshEditor()</p>\r
-\r
-<p class="whs3"> </p>\r
-\r
-<p class="whs4"># remove three elements: \r
- #850, #859 and #814</p>\r
-\r
-<p class="whs3">res = anEditor.RemoveElements([850, \r
- 859, 814])</p>\r
-\r
-<p class="whs3">if res == 1: print \r
- "Elements removing is OK!"</p>\r
-\r
-<p class="whs3">else: print \r
- "KO Elements removing." </p>\r
-\r
-<p class="whs3"> </p>\r
-\r
-<h3><a name=bookmark2>Renumbering Nodes and Elements</a></h3>\r
-\r
-<p class="whs1"><span style="font-weight: bold; color: #ff0000;"><font color=red ><B># \r
- Attention! This script has been written using the old approach basing \r
- on direct usage of SMESH idl interface.</B></font></span></p>\r
-\r
-<p class="whs1"># For the moment smesh package doesn't provide \r
- methods to renumber nodes and elements. </p>\r
-\r
-<p class="whs1"># In the next SALOME version the scripts will \r
- be updated to use only the commands from smesh package.</p>\r
-\r
-<p class="whs6"> </p>\r
-\r
-<p class="whs6">import SMESH_mechanic</p>\r
-\r
-<p class="whs6"> </p>\r
-\r
-<p class="whs6">mesh = SMESH_mechanic.mesh</p>\r
-\r
-<p class="whs6">anEditor = mesh.GetMeshEditor()</p>\r
-\r
-<p class="whs7"> </p>\r
-\r
-<p class="whs7"><span style="font-family: 'Lucida Console', monospace;">anEditor.RenumberNodes()</span> \r
- </p>\r
-\r
-<h3><a name=bookmark3>Moving Nodes</a></h3>\r
-\r
-<p class="whs1"><span style="font-weight: bold; color: #ff0000;"><font color=red ><B># \r
- Attention! This script has been written using the old approach basing \r
- on direct usage of SMESH idl interface.</B></font></span></p>\r
-\r
-<p class="whs1"># For the moment smesh package doesn't provide \r
- methods to move nodes. </p>\r
-\r
-<p class="whs1"># In the next SALOME version the scripts will \r
- be updated to use only the commands from smesh package.</p>\r
-\r
-<p class="whs3"> </p>\r
-\r
-<p class="whs3">import SMESH</p>\r
-\r
-<p class="whs3">import SMESH_mechanic</p>\r
-\r
-<p class="whs3"> </p>\r
-\r
-<p class="whs3">mesh = SMESH_mechanic.mesh</p>\r
-\r
-<p class="whs3">aMeshEditor = mesh.GetMeshEditor()</p>\r
-\r
-<p class="whs2"> </p>\r
-\r
-<p class="whs2"># move node #38</p>\r
-\r
-<p class="whs3">aMeshEditor.MoveNode(38, \r
- 20., 10., 0.) </p>\r
-\r
-<p class="whs2"> </p>\r
-\r
-<h3><a name=bookmark4>Diagonal Inversion</a></h3>\r
-\r
-<p class="whs1"><span style="font-weight: bold; color: #ff0000;"><font color=red ><B># \r
- Attention! This script has been written using the old approach basing \r
- on direct usage of SMESH idl interface.</B></font></span></p>\r
-\r
-<p class="whs1"># For the moment smesh package doesn't provide \r
- methods to produce a diagonal inversion. </p>\r
-\r
-<p class="whs1"># In the next SALOME version the scripts will \r
- be updated to use only the commands from smesh package.</p>\r
-\r
-<p class="whs7"> </p>\r
-\r
-<p class="whs6">import salome</p>\r
-\r
-<p class="whs6">import smesh</p>\r
-\r
-<p class="whs6"> </p>\r
-\r
-<p class="whs8"># create an empty mesh \r
- structure</p>\r
-\r
-<p class="whs6">gen = smesh.smesh</p>\r
-\r
-<p class="whs6">mesh = gen.CreateEmptyMesh()</p>\r
-\r
-<p class="whs6">aMeshEditor = mesh.GetMeshEditor()</p>\r
-\r
-<p class="whs6"> </p>\r
-\r
-<p class="whs8"># create the following \r
- mesh:</p>\r
-\r
-<p class="whs6"># .----.----.----.</p>\r
-\r
-<p class="whs6"># | /| \r
- /| \r
- /|</p>\r
-\r
-<p class="whs6"># | / \r
- | / | / \r
- |</p>\r
-\r
-<p class="whs6"># | / | \r
- / | / |</p>\r
-\r
-<p class="whs6"># |/ |/ \r
- |/ \r
- |</p>\r
-\r
-<p class="whs6"># .----.----.----.</p>\r
-\r
-<p class="whs6"> </p>\r
-\r
-<p class="whs6">bb = [0, 0, 0, 0]</p>\r
-\r
-<p class="whs6">tt = [0, 0, 0, 0]</p>\r
-\r
-<p class="whs6">ff = [0, 0, 0, 0, \r
- 0, 0]</p>\r
-\r
-<p class="whs6"> </p>\r
-\r
-<p class="whs6">bb[0] = aMeshEditor.AddNode( \r
- 0., 0., 0.)</p>\r
-\r
-<p class="whs6">bb[1] = aMeshEditor.AddNode(10., \r
- 0., 0.)</p>\r
-\r
-<p class="whs6">bb[2] = aMeshEditor.AddNode(20., \r
- 0., 0.)</p>\r
-\r
-<p class="whs6">bb[3] = aMeshEditor.AddNode(30., \r
- 0., 0.)</p>\r
-\r
-<p class="whs6"> </p>\r
-\r
-<p class="whs6">tt[0] = aMeshEditor.AddNode( \r
- 0., 15., 0.)</p>\r
-\r
-<p class="whs6">tt[1] = aMeshEditor.AddNode(10., \r
- 15., 0.)</p>\r
-\r
-<p class="whs6">tt[2] = aMeshEditor.AddNode(20., \r
- 15., 0.)</p>\r
-\r
-<p class="whs6">tt[3] = aMeshEditor.AddNode(30., \r
- 15., 0.)</p>\r
-\r
-<p class="whs6"> </p>\r
-\r
-<p class="whs6">ff[0] = aMeshEditor.AddFace([bb[0], \r
- bb[1], tt[1]])</p>\r
-\r
-<p class="whs6">ff[1] = aMeshEditor.AddFace([bb[0], \r
- tt[1], tt[0]])</p>\r
-\r
-<p class="whs6">ff[2] = aMeshEditor.AddFace([bb[1], \r
- bb[2], tt[2]])</p>\r
-\r
-<p class="whs6">ff[3] = aMeshEditor.AddFace([bb[1], \r
- tt[2], tt[1]])</p>\r
-\r
-<p class="whs6">ff[4] = aMeshEditor.AddFace([bb[2], \r
- bb[3], tt[3]])</p>\r
-\r
-<p class="whs6">ff[5] = aMeshEditor.AddFace([bb[2], \r
- tt[3], tt[2]])</p>\r
-\r
-<p class="whs6"> </p>\r
-\r
-<p class="whs8"># inverse the diagonal \r
- bb[1] - tt[2]</p>\r
-\r
-<p class="whs6">print "\nDiagonal \r
- inversion ... ",</p>\r
-\r
-<p class="whs6">res = aMeshEditor.InverseDiag(bb[1], \r
- tt[2])</p>\r
-\r
-<p class="whs6">if not res: print \r
- "failed!"</p>\r
-\r
-<p class="whs6">else: print \r
- "done."</p>\r
-\r
-<p class="whs7"> </p>\r
-\r
-<p class="whs7"><span style="font-family: 'Lucida Console', monospace;">salome.sg.updateObjBrowser(1)</span> \r
- </p>\r
-\r
-<h3><a name=bookmark5>Uniting two Triangles</a></h3>\r
-\r
-<p class="whs1"><span style="font-weight: bold; color: #ff0000;"><font color=red ><B># \r
- Attention! This script has been written using the old approach basing \r
- on direct usage of SMESH idl interface.</B></font></span></p>\r
-\r
-<p class="whs1"># For the moment smesh package doesn't provide \r
- methods to unite two triangles. </p>\r
-\r
-<p class="whs1"># In the next SALOME version the scripts will \r
- be updated to use only the commands from smesh package.</p>\r
-\r
-<p> </p>\r
-\r
-<p class="whs9">import salome</p>\r
-\r
-<p class="whs9">import smesh</p>\r
-\r
-<p> </p>\r
-\r
-<p># create an empty mesh structure</p>\r
-\r
-<p class="whs9">gen = smesh.smesh</p>\r
-\r
-<p class="whs9">mesh = gen.CreateEmptyMesh()</p>\r
-\r
-<p class="whs9">aMeshEditor = mesh.GetMeshEditor()</p>\r
-\r
-<p> </p>\r
-\r
-<p># create the following mesh:</p>\r
-\r
-<p class="whs9"># .----.----.----.</p>\r
-\r
-<p class="whs9"># | /| \r
- /| \r
- /|</p>\r
-\r
-<p class="whs9"># | / \r
- | / | / \r
- |</p>\r
-\r
-<p class="whs9"># | / | \r
- / | / |</p>\r
-\r
-<p class="whs9"># |/ |/ \r
- |/ \r
- |</p>\r
-\r
-<p class="whs9"># .----.----.----.</p>\r
-\r
-<p> </p>\r
-\r
-<p class="whs9">bb = [0, 0, 0, 0]</p>\r
-\r
-<p class="whs9">tt = [0, 0, 0, 0]</p>\r
-\r
-<p class="whs9">ff = [0, 0, 0, 0, \r
- 0, 0]</p>\r
-\r
-<p class="whs9"> </p>\r
-\r
-<p class="whs9">bb[0] = aMeshEditor.AddNode( \r
- 0., 0., 0.)</p>\r
-\r
-<p class="whs9">bb[1] = aMeshEditor.AddNode(10., \r
- 0., 0.)</p>\r
-\r
-<p class="whs9">bb[2] = aMeshEditor.AddNode(20., \r
- 0., 0.)</p>\r
-\r
-<p class="whs9">bb[3] = aMeshEditor.AddNode(30., \r
- 0., 0.)</p>\r
-\r
-<p class="whs9"> </p>\r
-\r
-<p class="whs9">tt[0] = aMeshEditor.AddNode( \r
- 0., 15., 0.)</p>\r
-\r
-<p class="whs9">tt[1] = aMeshEditor.AddNode(10., \r
- 15., 0.)</p>\r
-\r
-<p class="whs9">tt[2] = aMeshEditor.AddNode(20., \r
- 15., 0.)</p>\r
-\r
-<p class="whs9">tt[3] = aMeshEditor.AddNode(30., \r
- 15., 0.)</p>\r
-\r
-<p class="whs9"> </p>\r
-\r
-<p class="whs9">ff[0] = aMeshEditor.AddFace([bb[0], \r
- bb[1], tt[1]])</p>\r
-\r
-<p class="whs9">ff[1] = aMeshEditor.AddFace([bb[0], \r
- tt[1], tt[0]])</p>\r
-\r
-<p class="whs9">ff[2] = aMeshEditor.AddFace([bb[1], \r
- bb[2], tt[2]])</p>\r
-\r
-<p class="whs9">ff[3] = aMeshEditor.AddFace([bb[1], \r
- tt[2], tt[1]])</p>\r
-\r
-<p class="whs9">ff[4] = aMeshEditor.AddFace([bb[2], \r
- bb[3], tt[3]])</p>\r
-\r
-<p class="whs9">ff[5] = aMeshEditor.AddFace([bb[2], \r
- tt[3], tt[2]])</p>\r
-\r
-<p> </p>\r
-\r
-<p># delete the diagonal bb[1] - tt[2]</p>\r
-\r
-<p class="whs9">print "\nUnite \r
- two triangles ... ",</p>\r
-\r
-<p class="whs9">res = aMeshEditor.DeleteDiag(bb[1], \r
- tt[2])</p>\r
-\r
-<p class="whs9">if not res: print \r
- "failed!"</p>\r
-\r
-<p class="whs9">else: print \r
- "done."</p>\r
-\r
-<p class="whs9"> </p>\r
-\r
-<p><span style="font-family: 'Lucida Console', monospace;">salome.sg.updateObjBrowser(1)</span> \r
- </p>\r
-\r
-<h3><a name=bookmark6>Uniting a Set of Triangles</a></h3>\r
-\r
-<p class="whs1"><span style="font-weight: bold; color: #ff0000;"><font color=red ><B># \r
- Attention! This script has been written using the old approach basing \r
- on direct usage of SMESH idl interface.</B></font></span></p>\r
-\r
-<p class="whs1"># For the moment smesh package doesn't provide \r
- methods to unite a set of triangles. </p>\r
-\r
-<p class="whs1"># In the next SALOME version the scripts will \r
- be updated to use only the commands from smesh package.</p>\r
-\r
-<p class="whs10"> </p>\r
-\r
-<p class="whs9">import salome</p>\r
-\r
-<p class="whs9">import smesh</p>\r
-\r
-<p class="whs9"> </p>\r
-\r
-<p class="whs10"># create an empty mesh \r
- structure</p>\r
-\r
-<p class="whs9">gen = smesh.smesh</p>\r
-\r
-<p class="whs9">mesh = gen.CreateEmptyMesh()</p>\r
-\r
-<p class="whs9">aMeshEditor = mesh.GetMeshEditor()</p>\r
-\r
-<p class="whs9"> </p>\r
-\r
-<p class="whs10"># create the following \r
- mesh:</p>\r
-\r
-<p class="whs9"># .----.----.----.</p>\r
-\r
-<p class="whs9"># | /| \r
- /| \r
- /|</p>\r
-\r
-<p class="whs9"># | / \r
- | / | / \r
- |</p>\r
-\r
-<p class="whs9"># | / | \r
- / | / |</p>\r
-\r
-<p class="whs9"># |/ |/ \r
- |/ \r
- |</p>\r
-\r
-<p class="whs9"># .----.----.----.</p>\r
-\r
-<p class="whs9"> </p>\r
-\r
-<p class="whs9">bb = [0, 0, 0, 0]</p>\r
-\r
-<p class="whs9">tt = [0, 0, 0, 0]</p>\r
-\r
-<p class="whs9">ff = [0, 0, 0, 0, \r
- 0, 0]</p>\r
-\r
-<p class="whs9"> </p>\r
-\r
-<p class="whs9">bb[0] = aMeshEditor.AddNode( \r
- 0., 0., 0.)</p>\r
-\r
-<p class="whs9">bb[1] = aMeshEditor.AddNode(10., \r
- 0., 0.)</p>\r
-\r
-<p class="whs9">bb[2] = aMeshEditor.AddNode(20., \r
- 0., 0.)</p>\r
-\r
-<p class="whs9">bb[3] = aMeshEditor.AddNode(30., \r
- 0., 0.)</p>\r
-\r
-<p class="whs9"> </p>\r
-\r
-<p class="whs9">tt[0] = aMeshEditor.AddNode( \r
- 0., 15., 0.)</p>\r
-\r
-<p class="whs9">tt[1] = aMeshEditor.AddNode(10., \r
- 15., 0.)</p>\r
-\r
-<p class="whs9">tt[2] = aMeshEditor.AddNode(20., \r
- 15., 0.)</p>\r
-\r
-<p class="whs9">tt[3] = aMeshEditor.AddNode(30., \r
- 15., 0.)</p>\r
-\r
-<p class="whs9"> </p>\r
-\r
-<p class="whs9">ff[0] = aMeshEditor.AddFace([bb[0], \r
- bb[1], tt[1]])</p>\r
-\r
-<p class="whs9">ff[1] = aMeshEditor.AddFace([bb[0], \r
- tt[1], tt[0]])</p>\r
-\r
-<p class="whs9">ff[2] = aMeshEditor.AddFace([bb[1], \r
- bb[2], tt[2]])</p>\r
-\r
-<p class="whs9">ff[3] = aMeshEditor.AddFace([bb[1], \r
- tt[2], tt[1]])</p>\r
-\r
-<p class="whs9">ff[4] = aMeshEditor.AddFace([bb[2], \r
- bb[3], tt[3]])</p>\r
-\r
-<p class="whs9">ff[5] = aMeshEditor.AddFace([bb[2], \r
- tt[3], tt[2]])</p>\r
-\r
-<p class="whs9"> </p>\r
-\r
-<p class="whs10"># unite a set of triangles</p>\r
-\r
-<p class="whs9">aFilterMgr = smesh.smesh.CreateFilterManager()</p>\r
-\r
-<p class="whs9">aFunctor = aFilterMgr.CreateMinimumAngle()</p>\r
-\r
-<p class="whs9">print "\nUnite \r
- a set of triangles ... ",</p>\r
-\r
-<p class="whs9">res = aMeshEditor.TriToQuad([ff[2], \r
- ff[3], ff[4], ff[5]], aFunctor, 60.)</p>\r
-\r
-<p class="whs9">if not res: print \r
- "failed!"</p>\r
-\r
-<p class="whs9">else: print \r
- "done."</p>\r
-\r
-<p class="whs9"> </p>\r
-\r
-<p class="whs9">salome.sg.updateObjBrowser(1) \r
- </p>\r
-\r
-<h3><a name=bookmark12>Orientation</a></h3>\r
-\r
-<p class="whs1"><span style="font-weight: bold; color: #ff0000;"><font color=red ><B># \r
- Attention! This script has been written using the old approach basing \r
- on direct usage of SMESH idl interface.</B></font></span></p>\r
-\r
-<p class="whs1"># For the moment smesh package doesn't provide \r
- methods to change orientation. </p>\r
-\r
-<p class="whs1"># In the next SALOME version the scripts will \r
- be updated to use only the commands from smesh package.</p>\r
-\r
-<p class="whs9"> </p>\r
-\r
-<p class="whs9">import salome</p>\r
-\r
-<p class="whs9">import smesh</p>\r
-\r
-<p class="whs9"> </p>\r
-\r
-<p class="whs10"># SMESH module</p>\r
-\r
-<p class="whs9">gen = smesh.smesh</p>\r
-\r
-<p class="whs9">mesh = gen.CreateEmptyMesh()</p>\r
-\r
-<p class="whs9">MeshEditor = mesh.GetMeshEditor()</p>\r
-\r
-<p class="whs9"> </p>\r
-\r
-<p class="whs10"># build five quadrangles:</p>\r
-\r
-<p class="whs9">dx = 10</p>\r
-\r
-<p class="whs9">dy = 20</p>\r
-\r
-<p class="whs9"> </p>\r
-\r
-<p class="whs9">n1 = \r
- MeshEditor.AddNode(0.0 * dx, 0, 0)</p>\r
-\r
-<p class="whs9">n2 = \r
- MeshEditor.AddNode(1.0 * dx, 0, 0)</p>\r
-\r
-<p class="whs9">n3 = \r
- MeshEditor.AddNode(2.0 * dx, 0, 0)</p>\r
-\r
-<p class="whs9">n4 = \r
- MeshEditor.AddNode(3.0 * dx, 0, 0)</p>\r
-\r
-<p class="whs9">n5 = \r
- MeshEditor.AddNode(4.0 * dx, 0, 0)</p>\r
-\r
-<p class="whs9">n6 = \r
- MeshEditor.AddNode(5.0 * dx, 0, 0)</p>\r
-\r
-<p class="whs9">n7 = \r
- MeshEditor.AddNode(0.0 * dx, dy, 0)</p>\r
-\r
-<p class="whs9">n8 = \r
- MeshEditor.AddNode(1.0 * dx, dy, 0)</p>\r
-\r
-<p class="whs9">n9 = \r
- MeshEditor.AddNode(2.0 * dx, dy, 0)</p>\r
-\r
-<p class="whs9">n10 = MeshEditor.AddNode(3.0 \r
- * dx, dy, 0)</p>\r
-\r
-<p class="whs9">n11 = MeshEditor.AddNode(4.0 \r
- * dx, dy, 0)</p>\r
-\r
-<p class="whs9">n12 = MeshEditor.AddNode(5.0 \r
- * dx, dy, 0)</p>\r
-\r
-<p class="whs9"> </p>\r
-\r
-<p class="whs9">f1 = MeshEditor.AddFace([n1, \r
- n2, n8 , n7 ])</p>\r
-\r
-<p class="whs9">f2 = MeshEditor.AddFace([n2, \r
- n3, n9 , n8 ])</p>\r
-\r
-<p class="whs9">f3 = MeshEditor.AddFace([n3, \r
- n4, n10, n9 ])</p>\r
-\r
-<p class="whs9">f4 = MeshEditor.AddFace([n4, \r
- n5, n11, n10])</p>\r
-\r
-<p class="whs9">f5 = MeshEditor.AddFace([n5, \r
- n6, n12, n11])</p>\r
-\r
-<p class="whs9"> </p>\r
-\r
-<p class="whs10"># Change the orientation \r
- of the second and the fourth faces.</p>\r
-\r
-<p class="whs9">MeshEditor.Reorient([2, \r
- 4])</p>\r
-\r
-<p class="whs9"> </p>\r
-\r
-<p class="whs9">salome.sg.updateObjBrowser(1) \r
- </p>\r
-\r
-<p class="whs9"> </p>\r
-\r
-<h3><a name=bookmark7>Cutting Quadrangles</a></h3>\r
-\r
-<p class="whs1"><span style="font-weight: bold; color: #ff0000;"><font color=red ><B># \r
- Attention! This script has been written using the old approach basing \r
- on direct usage of SMESH idl interface.</B></font></span></p>\r
-\r
-<p class="whs1"># For the moment smesh package doesn't provide \r
- methods to cut quadrangles. </p>\r
-\r
-<p class="whs1"># In the next SALOME version the scripts will \r
- be updated to use only the commands from smesh package.</p>\r
-\r
-<p class="whs9"> </p>\r
-\r
-<p class="whs9">import SMESH</p>\r
-\r
-<p class="whs9">import SMESH_mechanic</p>\r
-\r
-<p class="whs9"> </p>\r
-\r
-<p class="whs9">smesh = SMESH_mechanic.smesh</p>\r
-\r
-<p class="whs9">mesh = \r
- SMESH_mechanic.mesh</p>\r
-\r
-<p class="whs9"> </p>\r
-\r
-<p class="whs10"># cut two quadrangles: \r
- 405 and 406</p>\r
-\r
-<p class="whs9">aFilterMgr = smesh.CreateFilterManager()</p>\r
-\r
-<p class="whs9">aFunctor = aFilterMgr.CreateMinimumAngle()</p>\r
-\r
-<p class="whs9"> </p>\r
-\r
-<p class="whs9">aMeshEditor = mesh.GetMeshEditor()</p>\r
-\r
-<p class="whs9">aMeshEditor.QuadToTri([405, \r
- 406], aFunctor) </p>\r
-\r
-<p class="whs10"> </p>\r
-\r
-<h3><a name=bookmark8>Smoothing</a></h3>\r
-\r
-<p class="whs1"><span style="font-weight: bold; color: #ff0000;"><font color=red ><B># \r
- Attention! This script has been written using the old approach basing \r
- on direct usage of SMESH idl interface.</B></font></span></p>\r
-\r
-<p class="whs1"># For the moment smesh package doesn't provide \r
- methods to perform smoothing. </p>\r
-\r
-<p class="whs1"># In the next SALOME version the scripts will \r
- be updated to use only the commands from smesh package.</p>\r
-\r
-<p class="whs1"> </p>\r
-\r
-<p class="whs9">import salome</p>\r
-\r
-<p class="whs9">import geompy</p>\r
-\r
-<p class="whs9"> </p>\r
-\r
-<p class="whs9">import SMESH</p>\r
-\r
-<p class="whs9">import SMESH_mechanic</p>\r
-\r
-<p class="whs9"> </p>\r
-\r
-<p class="whs9">mesh = SMESH_mechanic.mesh</p>\r
-\r
-<p class="whs9">aMeshEditor = mesh.GetMeshEditor()</p>\r
-\r
-<p class="whs9"> </p>\r
-\r
-<p class="whs10"># select the top face</p>\r
-\r
-<p class="whs9">faces = geompy.SubShapeAllSorted(SMESH_mechanic.shape_mesh, \r
- geompy.ShapeType["FACE"])</p>\r
-\r
-<p class="whs9">face = faces[3]</p>\r
-\r
-<p class="whs9">geompy.addToStudyInFather(SMESH_mechanic.shape_mesh, \r
- face, "face planar with hole")</p>\r
-\r
-<p class="whs9"> </p>\r
-\r
-<p class="whs10"># create a group of faces \r
- to be smoothed</p>\r
-\r
-<p class="whs9">GroupSmooth = mesh.CreateGroupFromGEOM(SMESH.FACE, \r
- "Group of faces (smooth)", face)</p>\r
-\r
-<p class="whs9"> </p>\r
-\r
-<p class="whs10"># perform smoothing</p>\r
-\r
-<p class="whs10"># boolean SmoothObject(Object, \r
- IDsOfFixedNodes, MaxNbOfIterations, MaxAspectRatio, Method)</p>\r
-\r
-<p class="whs9">res = aMeshEditor.SmoothObject(GroupSmooth, \r
- [], 20, 2., SMESH.SMESH_MeshEditor.CENTROIDAL_SMOOTH)</p>\r
-\r
-<p class="whs9"> </p>\r
-\r
-<p class="whs9">print "\nSmoothing \r
- ... ",</p>\r
-\r
-<p class="whs9">if not res: print \r
- "failed!"</p>\r
-\r
-<p class="whs9">else: print \r
- "done."</p>\r
-\r
-<p class="whs9"> </p>\r
-\r
-<p class="whs9">salome.sg.updateObjBrowser(1) \r
- </p>\r
-\r
-<p class="whs9"> </p>\r
-\r
-<h3><a name=bookmark9>Extrusion</a></h3>\r
-\r
-<p class="whs1"><span style="font-weight: bold; color: #ff0000;"><font color=red ><B># \r
- Attention! This script has been written using the old approach basing \r
- on direct usage of SMESH idl interface.</B></font></span></p>\r
-\r
-<p class="whs1"># For the moment smesh package doesn't provide \r
- methods to perform extrusion. </p>\r
-\r
-<p class="whs1"># In the next SALOME version the scripts will \r
- be updated to use only the commands from smesh package.</p>\r
-\r
-<p class="whs9"> </p>\r
-\r
-<p class="whs9">import salome</p>\r
-\r
-<p class="whs9">import geompy</p>\r
-\r
-<p class="whs9"> </p>\r
-\r
-<p class="whs9">import SMESH</p>\r
-\r
-<p class="whs9">import SMESH_mechanic</p>\r
-\r
-<p class="whs9"> </p>\r
-\r
-<p class="whs9">mesh = SMESH_mechanic.mesh</p>\r
-\r
-<p class="whs9">aMeshEditor = mesh.GetMeshEditor()</p>\r
-\r
-<p class="whs9"> </p>\r
-\r
-<p class="whs10"># select the top face</p>\r
-\r
-<p class="whs9">faces = geompy.SubShapeAllSorted(SMESH_mechanic.shape_mesh, \r
- geompy.ShapeType["FACE"])</p>\r
-\r
-<p class="whs9">face = faces[7]</p>\r
-\r
-<p class="whs9">geompy.addToStudyInFather(SMESH_mechanic.shape_mesh, \r
- face, "face circular top")</p>\r
-\r
-<p class="whs9"> </p>\r
-\r
-<p class="whs10"># create a vector for \r
- extrusion</p>\r
-\r
-<p class="whs9">point = SMESH.PointStruct(0., \r
- 0., 5.)</p>\r
-\r
-<p class="whs9">vector = SMESH.DirStruct(point)</p>\r
-\r
-<p class="whs9"> </p>\r
-\r
-<p class="whs10"># create a group to be \r
- extruded</p>\r
-\r
-<p class="whs9">GroupTri = mesh.CreateGroupFromGEOM(SMESH.FACE, \r
- "Group of faces (extrusion)", face)</p>\r
-\r
-<p class="whs9"> </p>\r
-\r
-<p class="whs10"># perform extrusion of \r
- the group</p>\r
-\r
-<p class="whs9">aMeshEditor.ExtrusionSweepObject(GroupTri, \r
- vector, 5)</p>\r
-\r
-<p class="whs10"> </p>\r
-\r
-<p class="whs10"><span style="font-family: 'Lucida Console', monospace;">salome.sg.updateObjBrowser(1)</span> \r
- </p>\r
-\r
-<h3><a name=bookmark10>Extrusion along a Path</a></h3>\r
-\r
-<p class="whs1"><span style="font-weight: bold; color: #ff0000;"><font color=red ><B># \r
- Attention! This script has been written using the old approach basing \r
- on direct usage of SMESH idl interface.</B></font></span></p>\r
-\r
-<p class="whs1"># For the moment smesh package doesn't provide \r
- methods to perform extrusion along a path. </p>\r
-\r
-<p class="whs1"># In the next SALOME version the scripts will \r
- be updated to use only the commands from smesh package.</p>\r
-\r
-<p class="whs9"> </p>\r
-\r
-<p class="whs9">import math</p>\r
-\r
-<p class="whs9">import salome</p>\r
-\r
-<p class="whs9"> </p>\r
-\r
-<p class="whs10"># Geometry</p>\r
-\r
-<p class="whs9">import geompy</p>\r
-\r
-<p class="whs10"> </p>\r
-\r
-<p class="whs10"># 1. Create points</p>\r
-\r
-<p class="whs9">points = [[0, 0], \r
- [50, 30], [50, 110], [0, 150], [-80, 150], [-130, 70], [-130, -20]]</p>\r
-\r
-<p class="whs9"> </p>\r
-\r
-<p class="whs9">iv = 1</p>\r
-\r
-<p class="whs9">vertices = []</p>\r
-\r
-<p class="whs9">for point in points:</p>\r
-\r
-<p class="whs11">vert \r
- = geompy.MakeVertex(point[0], point[1], 0)</p>\r
-\r
-<p class="whs11">geompy.addToStudy(vert, \r
- "Vertex_" + `iv`)</p>\r
-\r
-<p class="whs11">vertices.append(vert)</p>\r
-\r
-<p class="whs11">iv \r
- += 1</p>\r
-\r
-<p class="whs11">pass</p>\r
-\r
-<p class="whs10"> </p>\r
-\r
-<p class="whs10"># 2. Create edges and \r
- wires</p>\r
-\r
-<p class="whs9">Edge_straight = geompy.MakeEdge(vertices[0], \r
- vertices[4])</p>\r
-\r
-<p class="whs9">Edge_bezierrr = geompy.MakeBezier(vertices)</p>\r
-\r
-<p class="whs9">Wire_polyline = geompy.MakePolyline(vertices)</p>\r
-\r
-<p class="whs9">Edge_Circle = \r
- geompy.MakeCircleThreePnt(vertices[0], vertices[1], vertices[2])</p>\r
-\r
-<p class="whs9"> </p>\r
-\r
-<p class="whs9">geompy.addToStudy(Edge_straight, \r
- "Edge_straight")</p>\r
-\r
-<p class="whs9">geompy.addToStudy(Edge_bezierrr, \r
- "Edge_bezierrr")</p>\r
-\r
-<p class="whs9">geompy.addToStudy(Wire_polyline, \r
- "Wire_polyline")</p>\r
-\r
-<p class="whs9">geompy.addToStudy(Edge_Circle \r
- , "Edge_Circle")</p>\r
-\r
-<p class="whs10"> </p>\r
-\r
-<p class="whs10"># 3. Explode wire on \r
- edges, as they will be used for mesh extrusion</p>\r
-\r
-<p class="whs9">Wire_polyline_edges \r
- = geompy.SubShapeAll(Wire_polyline, geompy.ShapeType["EDGE"])</p>\r
-\r
-<p class="whs9">for ii in range(len(Wire_polyline_edges)):</p>\r
-\r
-<p class="whs11">geompy.addToStudyInFather(Wire_polyline, \r
- Wire_polyline_edges[ii], "Edge_" + `ii + 1`)</p>\r
-\r
-<p class="whs11">pass</p>\r
-\r
-<p class="whs12"> </p>\r
-\r
-<p class="whs10"># Mesh</p>\r
-\r
-<p class="whs9">import smesh</p>\r
-\r
-<p class="whs9">import SMESH</p>\r
-\r
-<p class="whs9"> </p>\r
-\r
-<p class="whs9">gen = smesh.smesh</p>\r
-\r
-<p class="whs9"> </p>\r
-\r
-<p class="whs9">smeshgui = salome.ImportComponentGUI("SMESH")</p>\r
-\r
-<p class="whs9">smeshgui.Init(salome.myStudyId)</p>\r
-\r
-<p class="whs10"> </p>\r
-\r
-<p class="whs10"># 1D algorithm and three \r
- 1D hypotheses</p>\r
-\r
-<p class="whs9">Wire_discretisation \r
- = gen.CreateHypothesis('Regular_1D', 'libStdMeshersEngine.so')</p>\r
-\r
-<p class="whs9"> </p>\r
-\r
-<p class="whs9">Nb_Segments_3 = gen.CreateHypothesis('NumberOfSegments', \r
- 'libStdMeshersEngine.so')</p>\r
-\r
-<p class="whs9">Nb_Segments_7 = gen.CreateHypothesis('NumberOfSegments', \r
- 'libStdMeshersEngine.so')</p>\r
-\r
-<p class="whs9">Nb_Segments_8 = gen.CreateHypothesis('NumberOfSegments', \r
- 'libStdMeshersEngine.so')</p>\r
-\r
-<p class="whs9"> </p>\r
-\r
-<p class="whs9">Nb_Segments_3.SetNumberOfSegments(3)</p>\r
-\r
-<p class="whs9">Nb_Segments_7.SetNumberOfSegments(7)</p>\r
-\r
-<p class="whs9">Nb_Segments_8.SetNumberOfSegments(8)</p>\r
-\r
-<p class="whs10"> </p>\r
-\r
-<p class="whs10"># Mesh given shape with \r
- given 1d hypothesis</p>\r
-\r
-<p class="whs9">def Mesh1D(shape1d, \r
- hyp1d, name):</p>\r
-\r
-<p class="whs11">mesh1d_tool \r
- = smesh.Mesh(shape1d)</p>\r
-\r
-<p class="whs11">mesh1d \r
- = mesh1d_tool.GetMesh()</p>\r
-\r
-<p class="whs11">status \r
- = mesh1d.AddHypothesis(shape1d, hyp1d)</p>\r
-\r
-<p class="whs11">status \r
- = mesh1d.AddHypothesis(shape1d, Wire_discretisation)</p>\r
-\r
-<p class="whs11">isDone \r
- = mesh1d_tool.Compute()</p>\r
-\r
-<p class="whs11">if \r
- not isDone: print 'Mesh ', name, ': computation failed'</p>\r
-\r
-<p class="whs11">return \r
- mesh1d</p>\r
-\r
-<p class="whs9"> </p>\r
-\r
-<p class="whs10"># Create a mesh with \r
- six nodes, seven edges and two quadrangle faces</p>\r
-\r
-<p class="whs9">def MakeQuadMesh2(mesh_name):</p>\r
-\r
-<p class="whs11">quad_1 \r
- = gen.CreateEmptyMesh()</p>\r
-\r
-<p class="whs11">smeshgui.SetName(salome.ObjectToID(quad_1), \r
- mesh_name)</p>\r
-\r
-<p class="whs11">editor_1 \r
- = quad_1.GetMeshEditor()</p>\r
-\r
-<p class="whs12"> </p>\r
-\r
-<p class="whs12"># six nodes</p>\r
-\r
-<p class="whs11">n1 \r
- = editor_1.AddNode(0, 20, 10)</p>\r
-\r
-<p class="whs11">n2 \r
- = editor_1.AddNode(0, 40, 10)</p>\r
-\r
-<p class="whs11">n3 \r
- = editor_1.AddNode(0, 40, 30)</p>\r
-\r
-<p class="whs11">n4 \r
- = editor_1.AddNode(0, 20, 30)</p>\r
-\r
-<p class="whs11">n5 \r
- = editor_1.AddNode(0, 0, \r
- 30)</p>\r
-\r
-<p class="whs11">n6 \r
- = editor_1.AddNode(0, 0, \r
- 10)</p>\r
-\r
-<p class="whs12"> </p>\r
-\r
-<p class="whs12"># seven edges</p>\r
-\r
-<p class="whs11">editor_1.AddEdge([n1, \r
- n2]) # 1</p>\r
-\r
-<p class="whs11">editor_1.AddEdge([n2, \r
- n3]) # 2</p>\r
-\r
-<p class="whs11">editor_1.AddEdge([n3, \r
- n4]) # 3</p>\r
-\r
-<p class="whs11">editor_1.AddEdge([n4, \r
- n1]) # 4</p>\r
-\r
-<p class="whs11">editor_1.AddEdge([n4, \r
- n5]) # 5</p>\r
-\r
-<p class="whs11">editor_1.AddEdge([n5, \r
- n6]) # 6</p>\r
-\r
-<p class="whs11">editor_1.AddEdge([n6, \r
- n1]) # 7</p>\r
-\r
-<p class="whs12"> </p>\r
-\r
-<p class="whs12"># two quadrangle faces</p>\r
-\r
-<p class="whs11">editor_1.AddFace([n1, \r
- n2, n3, n4]) # 8</p>\r
-\r
-<p class="whs11">editor_1.AddFace([n1, \r
- n4, n5, n6]) # 9</p>\r
-\r
-<p class="whs11">return \r
- [quad_1, editor_1, [1,2,3,4,5,6,7], [8,9]]</p>\r
-\r
-<p class="whs10"> </p>\r
-\r
-<p class="whs10"># Path meshes</p>\r
-\r
-<p class="whs9">Edge_straight_mesh \r
- = Mesh1D(Edge_straight, Nb_Segments_7, "Edge_straight")</p>\r
-\r
-<p class="whs9">Edge_bezierrr_mesh \r
- = Mesh1D(Edge_bezierrr, Nb_Segments_7, "Edge_bezierrr")</p>\r
-\r
-<p class="whs9">Wire_polyline_mesh \r
- = Mesh1D(Wire_polyline, Nb_Segments_3, "Wire_polyline")</p>\r
-\r
-<p class="whs9">Edge_Circle_mesh = \r
- Mesh1D(Edge_Circle , \r
- Nb_Segments_8, "Edge_Circle")</p>\r
-\r
-<p class="whs10"> </p>\r
-\r
-<p class="whs10"># Initial meshes (to \r
- be extruded)</p>\r
-\r
-<p class="whs9">[quad_1, editor_1, \r
- ee_1, ff_1] = MakeQuadMesh2("quad_1")</p>\r
-\r
-<p class="whs9">[quad_2, editor_2, \r
- ee_2, ff_2] = MakeQuadMesh2("quad_2")</p>\r
-\r
-<p class="whs9">[quad_3, editor_3, \r
- ee_3, ff_3] = MakeQuadMesh2("quad_3")</p>\r
-\r
-<p class="whs9">[quad_4, editor_4, \r
- ee_4, ff_4] = MakeQuadMesh2("quad_4")</p>\r
-\r
-<p class="whs9">[quad_5, editor_5, \r
- ee_5, ff_5] = MakeQuadMesh2("quad_5")</p>\r
-\r
-<p class="whs9">[quad_6, editor_6, \r
- ee_6, ff_6] = MakeQuadMesh2("quad_6")</p>\r
-\r
-<p class="whs9">[quad_7, editor_7, \r
- ee_7, ff_7] = MakeQuadMesh2("quad_7")</p>\r
-\r
-<p class="whs10"> </p>\r
-\r
-<p class="whs10"># ExtrusionAlongPath</p>\r
-\r
-<p class="whs10"># IDsOfElements, PathMesh, \r
- PathShape, NodeStart,</p>\r
-\r
-<p class="whs10"># HasAngles, Angles, \r
- HasRefPoint, RefPoint</p>\r
-\r
-<p class="whs10"> </p>\r
-\r
-<p class="whs9">refPoint = SMESH.PointStruct(0, \r
- 0, 0)</p>\r
-\r
-<p class="whs9"> </p>\r
-\r
-<p class="whs9">a10 = 10.0*math.pi/180.0</p>\r
-\r
-<p class="whs9">a45 = 45.0*math.pi/180.0</p>\r
-\r
-<p class="whs10"> </p>\r
-\r
-<p class="whs10"># 1. Extrusion of two \r
- mesh edges along a straight path</p>\r
-\r
-<p class="whs9">error = editor_1.ExtrusionAlongPath([1,2], \r
- Edge_straight_mesh, Edge_straight, 1,</p>\r
-\r
-<p class="whs13"> 0, \r
- [], 0, refPoint)</p>\r
-\r
-<p class="whs10"> </p>\r
-\r
-<p class="whs10"># 2. Extrusion of one \r
- mesh edge along a curved path</p>\r
-\r
-<p class="whs9">error = editor_2.ExtrusionAlongPath([2], \r
- Edge_bezierrr_mesh, Edge_bezierrr, 1,</p>\r
-\r
-<p class="whs13"> 0, \r
- [], 0, refPoint)</p>\r
-\r
-<p class="whs10"> </p>\r
-\r
-<p class="whs10"># 3. Extrusion of one \r
- mesh edge along a curved path with usage of angles </p>\r
-\r
-<p class="whs9">error = editor_3.ExtrusionAlongPath([2], \r
- Edge_bezierrr_mesh, Edge_bezierrr, 1,</p>\r
-\r
-<p class="whs13"> 1, \r
- [a45, a45, a45, 0, -a45, -a45, -a45], 0, refPoint)</p>\r
-\r
-<p class="whs10"> </p>\r
-\r
-<p class="whs10"># 4. Extrusion of one \r
- mesh edge along the path, which is a part of a meshed wire</p>\r
-\r
-<p class="whs9">error = editor_4.ExtrusionAlongPath([4], \r
- Wire_polyline_mesh, Wire_polyline_edges[0], 1,</p>\r
-\r
-<p class="whs13"> 1, \r
- [a10, a10, a10], 0, refPoint)</p>\r
-\r
-<p class="whs10"> </p>\r
-\r
-<p class="whs10"># 5. Extrusion of two \r
- mesh faces along the path, which is a part of a meshed wire</p>\r
-\r
-<p class="whs9">error = editor_5.ExtrusionAlongPath(ff_5 \r
- , Wire_polyline_mesh, Wire_polyline_edges[2], 4,</p>\r
-\r
-<p class="whs13"> 0, \r
- [], 0, refPoint)</p>\r
-\r
-<p class="whs10"> </p>\r
-\r
-<p class="whs10"># 6. Extrusion of two \r
- mesh faces along a closed path</p>\r
-\r
-<p class="whs9">error = editor_6.ExtrusionAlongPath(ff_6 \r
- , Edge_Circle_mesh, Edge_Circle, 1,</p>\r
-\r
-<p class="whs13"> 0, \r
- [], 0, refPoint)</p>\r
-\r
-<p class="whs10"> </p>\r
-\r
-<p class="whs10"># 7. Extrusion of two \r
- mesh faces along a closed path with usage of angles </p>\r
-\r
-<p class="whs9">error = editor_7.ExtrusionAlongPath(ff_7, \r
- Edge_Circle_mesh, Edge_Circle, 1,</p>\r
-\r
-<p class="whs13"> 1, \r
- [a45, -a45, a45, -a45, a45, -a45, a45, -a45], 0, refPoint</p>\r
-\r
-<p class="whs13"> </p>\r
-\r
-<p class="whs14">salome.sg.updateObjBrowser(1)</p>\r
-\r
-<h3><a name=bookmark11>Revolution</a></h3>\r
-\r
-<p class="whs1"><span style="font-weight: bold; color: #ff0000;"><font color=red ><B># \r
- Attention! This script has been written using the old approach basing \r
- on direct usage of SMESH idl interface.</B></font></span></p>\r
-\r
-<p class="whs1"># For the moment smesh package doesn't provide \r
- methods to perform revolution. </p>\r
-\r
-<p class="whs1"># In the next SALOME version the scripts will \r
- be updated to use only the commands from smesh package.</p>\r
-\r
-<p class="whs1"> </p>\r
-\r
-<p class="whs9">import math</p>\r
-\r
-<p class="whs9">import SMESH</p>\r
-\r
-<p class="whs9">import SMESH_mechanic</p>\r
-\r
-<p class="whs9"> </p>\r
-\r
-<p class="whs9">mesh = SMESH_mechanic.mesh</p>\r
-\r
-<p class="whs9">aMeshEditor = mesh.GetMeshEditor()</p>\r
-\r
-<p class="whs9"> </p>\r
-\r
-<p class="whs10"># create a group of faces \r
- to be revolved</p>\r
-\r
-<p class="whs9">FacesRotate = [492, \r
- 493, 502, 503]</p>\r
-\r
-<p class="whs9">GroupRotate = mesh.CreateGroup(SMESH.FACE,"Group \r
- of faces (rotate)")</p>\r
-\r
-<p class="whs9">GroupRotate.Add(FacesRotate)</p>\r
-\r
-<p class="whs10"> </p>\r
-\r
-<p class="whs10"># define revolution angle \r
- and axis</p>\r
-\r
-<p class="whs9">angle45 = 45 * math.pi \r
- / 180</p>\r
-\r
-<p class="whs9">axisXYZ = SMESH.AxisStruct(-38.3128, \r
- -73.3658, -23.321, -13.3402, -13.3265, 6.66632)</p>\r
-\r
-<p class="whs9"> </p>\r
-\r
-<p class="whs10"># perform revolution \r
- of an object</p>\r
-\r
-<p class="whs9">aMeshEditor.RotationSweepObject(GroupRotate, \r
- axisXYZ, angle45, 4, 1e-5) </p>\r
-\r
-<p class="whs10"> </p>\r
-\r
-<h3><a name=bookmark13>Pattern Mapping</a></h3>\r
-\r
-<p class="whs1"><span style="font-weight: bold; color: #ff0000;"><font color=red ><B># \r
- Attention! This script has been written using the old approach basing \r
- on direct usage of SMESH idl interface.</B></font></span></p>\r
-\r
-<p class="whs1"># For the moment smesh package doesn't provide \r
- methods to perform pattern mapping. </p>\r
-\r
-<p class="whs1"># In the next SALOME version the scripts will \r
- be updated to use only the commands from smesh package.</p>\r
-\r
-<p class="whs9"> </p>\r
-\r
-<p class="whs9">import geompy</p>\r
-\r
-<p class="whs9">import smesh</p>\r
-\r
-<p class="whs9"> </p>\r
-\r
-<p class="whs10"># define the geometry</p>\r
-\r
-<p class="whs9">Box_1 = geompy.MakeBoxDXDYDZ(200., \r
- 200., 200.)</p>\r
-\r
-<p class="whs9">geompy.addToStudy(Box_1, \r
- "Box_1")</p>\r
-\r
-<p class="whs9"> </p>\r
-\r
-<p class="whs9">faces = geompy.SubShapeAll(Box_1, \r
- geompy.ShapeType["FACE"])</p>\r
-\r
-<p class="whs9">Face_1 = faces[0]</p>\r
-\r
-<p class="whs9">Face_2 = faces[1]</p>\r
-\r
-<p class="whs9">geompy.addToStudyInFather(Box_1, \r
- Face_1, "Face_1")</p>\r
-\r
-<p class="whs9">geompy.addToStudyInFather(Box_1, \r
- Face_2, "Face_2")</p>\r
-\r
-<p class="whs9"> </p>\r
-\r
-<p class="whs10"># build a quadrangle \r
- mesh 3x3 on Face_1</p>\r
-\r
-<p class="whs9">Mesh_1 = smesh.Mesh(Face_1)</p>\r
-\r
-<p class="whs9">algo1D = Mesh_1.Segment()</p>\r
-\r
-<p class="whs9">algo1D.NumberOfSegments(3)</p>\r
-\r
-<p class="whs9">Mesh_1.Quadrangle()</p>\r
-\r
-<p class="whs9"> </p>\r
-\r
-<p class="whs9">isDone = Mesh_1.Compute()</p>\r
-\r
-<p class="whs9">if not isDone: print \r
- 'Mesh Mesh_1 : computation failed'</p>\r
-\r
-<p class="whs9"> </p>\r
-\r
-<p class="whs10"># build a triangle mesh \r
- on Face_2</p>\r
-\r
-<p class="whs9">Mesh_2 = smesh.Mesh(Face_2)</p>\r
-\r
-<p class="whs9">algo1D = Mesh_2.Segment()</p>\r
-\r
-<p class="whs9">algo1D.NumberOfSegments(1)</p>\r
-\r
-<p class="whs9">algo2D = Mesh_2.Triangle()</p>\r
-\r
-<p class="whs9">algo2D.MaxElementArea(240)</p>\r
-\r
-<p class="whs9">isDone = Mesh_2.Compute()</p>\r
-\r
-<p class="whs9">if not isDone: print \r
- 'Mesh Mesh_2 : computation failed'</p>\r
-\r
-<p class="whs9"> </p>\r
-\r
-<p class="whs10"># create a pattern</p>\r
-\r
-<p class="whs9">pattern = smesh.smesh.GetPattern()</p>\r
-\r
-<p class="whs9">isDone = pattern.LoadFromFace(Mesh_2.GetMesh(), \r
- Face_2, 0)</p>\r
-\r
-<p class="whs9">if (isDone != 1): \r
- print 'LoadFromFace :', pattern.GetErrorCode()</p>\r
-\r
-<p class="whs9"> </p>\r
-\r
-<p class="whs10"># apply the pattern to \r
- a face of the first mesh</p>\r
-\r
-<p class="whs9">pattern.ApplyToMeshFaces(Mesh_1.GetMesh(), \r
- [17], 0, 0)</p>\r
-\r
-<p class="whs9">isDone = pattern.MakeMesh(Mesh_1.GetMesh(), \r
- 0, 0)</p>\r
-\r
-<p class="whs9">if (isDone != 1): \r
- print 'MakeMesh :', pattern.GetErrorCode() </p>\r
-\r
-<p class="whs9"> </p>\r
-\r
-<p class="whs10"> </p>\r
-\r
-<script type="text/javascript" language="javascript1.2">\r
-<!--\r
-if (window.writeIntopicBar)\r
- writeIntopicBar(0);\r
-//-->\r
-</script>\r
-</body>\r
-</html>\r
+<!doctype HTML public "-//W3C//DTD HTML 4.0 Frameset//EN">
+
+<html>
+
+<head>
+<title>Modifying Meshes</title>
+<meta http-equiv="content-type" content="text/html; charset=windows-1252">
+<meta name="generator" content="RoboHelp by eHelp Corporation www.ehelp.com"><style>
+<!--
+P { margin-top:0pt; margin-bottom:0pt; }
+LI.kadov-P { }
+-->
+</style><style type="text/css">
+<!--
+p.whs1 { margin-top:0pt; margin-bottom:0pt; font-family:'Lucida Console' , monospace; }
+p.whs2 { margin-top:0pt; margin-bottom:0pt; font-family:'Times New Roman' , serif; }
+p.whs3 { margin-top:0pt; margin-bottom:0pt; }
+h4.whs4 { margin-top:0pt; margin-bottom:0pt; }
+p.whs5 { margin-top:0.5pt; margin-bottom:0pt; font-family:'Lucida Console' , monospace; }
+p.whs6 { margin-top:0.5pt; margin-bottom:0pt; }
+p.whs7 { margin-top:0.5pt; margin-bottom:0pt; font-family:'Times New Roman' , serif; }
+p.whs8 { font-family:'Lucida Console' , monospace; }
+p.whs9 { font-family:'Times New Roman' , serif; }
+p.whs10 { margin-left:36px; font-family:'Lucida Console' , monospace; }
+p.whs11 { margin-left:36px; }
+p.whs12 { margin-left:0px; font-family:'Times New Roman' , serif; }
+p.whs13 { font-family:'Lucida Console' , monospace; margin-left:0px; }
+-->
+</style><script type="text/javascript" language="JavaScript">
+<!--
+if ((navigator.appName == "Netscape") && (parseInt(navigator.appVersion) == 4))
+{
+ var strNSS = "<style type='text/css'>";
+ strNSS += "p.whs1 {margin-top:1pt;margin-bottom:1pt; }";
+ strNSS += "p.whs2 {margin-top:1pt;margin-bottom:1pt; }";
+ strNSS += "p.whs3 {margin-top:1pt;margin-bottom:1pt; }";
+ strNSS += "h4.whs4 {margin-top:1pt;margin-bottom:1pt; }";
+ strNSS += "p.whs5 {margin-top:1pt;margin-bottom:1pt; }";
+ strNSS += "p.whs6 {margin-top:1pt;margin-bottom:1pt; }";
+ strNSS += "p.whs7 {margin-top:1pt;margin-bottom:1pt; }";
+ strNSS += "p.whs12 {margin-left:1pt; }";
+ strNSS += "p.whs13 {margin-left:1pt; }";
+ strNSS +="</style>";
+ document.write(strNSS);
+}
+//-->
+</script>
+<script type="text/javascript" language="JavaScript" title="WebHelpInlineScript">
+<!--
+function reDo() {
+ if (innerWidth != origWidth || innerHeight != origHeight)
+ location.reload();
+}
+if ((parseInt(navigator.appVersion) == 4) && (navigator.appName == "Netscape")) {
+ origWidth = innerWidth;
+ origHeight = innerHeight;
+ onresize = reDo;
+}
+onerror = null;
+//-->
+</script>
+<style type="text/css">
+<!--
+div.WebHelpPopupMenu { position:absolute; left:0px; top:0px; z-index:4; visibility:hidden; }
+p.WebHelpNavBar { text-align:right; }
+-->
+</style><script type="text/javascript" language="javascript1.2" src="whmsg.js"></script>
+<script type="text/javascript" language="javascript" src="whver.js"></script>
+<script type="text/javascript" language="javascript1.2" src="whproxy.js"></script>
+<script type="text/javascript" language="javascript1.2" src="whutils.js"></script>
+<script type="text/javascript" language="javascript1.2" src="whtopic.js"></script>
+<script type="text/javascript" language="javascript1.2">
+<!--
+if (window.gbWhTopic)
+{
+ if (window.setRelStartPage)
+ {
+ addTocInfo("MESH module\nTUI Scripts\nModifying Meshes");
+addButton("show",BTN_IMG,"Show","","","","",0,0,"whd_show0.gif","whd_show2.gif","whd_show1.gif");
+addButton("hide",BTN_IMG,"Hide","","","","",0,0,"whd_hide0.gif","whd_hide2.gif","whd_hide1.gif");
+
+ }
+
+
+ if (window.setRelStartPage)
+ {
+ setRelStartPage("index.htm");
+
+ autoSync(1);
+ sendSyncInfo();
+ sendAveInfoOut();
+ }
+
+}
+else
+ if (window.gbIE4)
+ document.location.reload();
+//-->
+</script>
+</head>
+<body><script type="text/javascript" language="javascript1.2">
+<!--
+if (window.writeIntopicBar)
+ writeIntopicBar(4);
+//-->
+</script>
+<h1>Modifying Meshes</h1>
+
+<h3><a name=bookmark>Adding Nodes and Elements</a></h3>
+
+<h4>Add Node</h4>
+
+<p class="whs1">import SMESH_mechanic</p>
+
+<p class="whs1"> </p>
+
+<p class="whs1">mesh = SMESH_mechanic.mesh</p>
+
+<p class="whs1"> </p>
+
+<p class="whs2"># add node</p>
+
+<p class="whs1"><span style="font-family: 'Lucida Console', monospace;">new_id
+ = mesh.AddNode(50, 10, 0)</span></p>
+
+<p class="whs1">print ""</p>
+
+<p class="whs1">if new_id == 0:
+ print "KO node addition."</p>
+
+<p class="whs1">else: print
+ "New Node has been added with ID ", new_id </p>
+
+<p class="whs1"> </p>
+
+<h4>Add Edge</h4>
+
+<p class="whs1">import SMESH_mechanic</p>
+
+<p class="whs1"> </p>
+
+<p class="whs1">mesh = SMESH_mechanic.mesh</p>
+
+<p class="whs1">print ""</p>
+
+<p class="whs1"> </p>
+
+<p class="whs2"># add node</p>
+
+<p class="whs1"><span style="font-family: 'Lucida Console', monospace;">n1
+ = mesh.AddNode(50, 10, 0)</span></p>
+
+<p class="whs1">if n1 == 0: print
+ "KO node addition." </p>
+
+<p class="whs1"> </p>
+
+<p class="whs2"># add edge</p>
+
+<p class="whs1"><span style="font-family: 'Lucida Console', monospace;">e1
+ = mesh.AddEdge([n1, 38])</span></p>
+
+<p class="whs1">if e1 == 0: print
+ "KO edge addition."</p>
+
+<p class="whs1">else: print
+ "New Edge has been added with ID ", e1 </p>
+
+<h4>Add Triangle</h4>
+
+<p class="whs1">import SMESH_mechanic</p>
+
+<p class="whs1"> </p>
+
+<p class="whs1">mesh = SMESH_mechanic.mesh</p>
+
+<p class="whs1">print ""</p>
+
+<p class="whs1"> </p>
+
+<p class="whs2"># add node</p>
+
+<p class="whs1">n1 = mesh.AddNode(50,
+ 10, 0)</p>
+
+<p class="whs1">if n1 == 0: print
+ "KO node addition."</p>
+
+<p class="whs1"> </p>
+
+<p class="whs2"># add triangle</p>
+
+<p class="whs1">t1 = mesh.AddFace([n1,
+ 38, 39])</p>
+
+<p class="whs1">if t1 == 0: print
+ "KO triangle addition."</p>
+
+<p class="whs1">else: print
+ "New Triangle has been added with ID ", t1 </p>
+
+<h4>Add Quadrangle</h4>
+
+<p class="whs1">import SMESH_mechanic</p>
+
+<p class="whs1"> </p>
+
+<p class="whs1">mesh = SMESH_mechanic.mesh</p>
+
+<p class="whs1">print ""</p>
+
+<p class="whs3"> </p>
+
+<p class="whs3"># add node</p>
+
+<p class="whs1">n1 = mesh.AddNode(50,
+ 10, 0)</p>
+
+<p class="whs1">if n1 == 0: print
+ "KO node addition."</p>
+
+<p class="whs1"> </p>
+
+<p class="whs1">n2 = mesh.AddNode(40,
+ 20, 0)</p>
+
+<p class="whs1">if n2 == 0: print
+ "KO node addition."</p>
+
+<p class="whs1"> </p>
+
+<p class="whs3"># add quadrangle</p>
+
+<p class="whs1">q1 = mesh.AddFace([n2,
+ n1, 38, 39])</p>
+
+<p class="whs1">if q1 == 0: print
+ "KO quadrangle addition."</p>
+
+<p class="whs3"><span style="font-family: 'Lucida Console', monospace;">else:
+ print
+ "New Quadrangle has been added with ID ", q1</span> </p>
+
+<h4>Add Tetrahedron</h4>
+
+<p class="whs1">import SMESH_mechanic</p>
+
+<p class="whs1"> </p>
+
+<p class="whs1">mesh = SMESH_mechanic.mesh</p>
+
+<p class="whs1">print ""</p>
+
+<p class="whs1"> </p>
+
+<p class="whs2"># add node</p>
+
+<p class="whs1">n1 = mesh.AddNode(50,
+ 10, 0)</p>
+
+<p class="whs1">if n1 == 0: print
+ "KO node addition."</p>
+
+<p class="whs1"> </p>
+
+<p class="whs2"># add tetrahedron</p>
+
+<p class="whs1">t1 = mesh.AddVolume([n1,
+ 38, 39, 246])</p>
+
+<p class="whs1">if t1 == 0: print
+ "KO tetrahedron addition."</p>
+
+<p class="whs3"><span style="font-family: 'Lucida Console', monospace;">else:
+ print
+ "New Tetrahedron has been added with ID ", t1</span> </p>
+
+<h4>Add Hexahedron</h4>
+
+<p class="whs1">import SMESH_mechanic</p>
+
+<p class="whs1"> </p>
+
+<p class="whs1">mesh = SMESH_mechanic.mesh</p>
+
+<p class="whs1">print ""</p>
+
+<p class="whs1"> </p>
+
+<p class="whs2"># add nodes</p>
+
+<p class="whs1">nId1 = mesh.AddNode(50,
+ 10, 0)</p>
+
+<p class="whs1">nId2 = mesh.AddNode(47,
+ 12, 0)</p>
+
+<p class="whs1">nId3 = mesh.AddNode(50,
+ 10, 10)</p>
+
+<p class="whs1">nId4 = mesh.AddNode(47,
+ 12, 10)</p>
+
+<p class="whs1"> </p>
+
+<p class="whs1">if nId1 == 0 or
+ nId2 == 0 or nId3 == 0 or nId4 == 0: print "KO node addition."</p>
+
+<p class="whs1"> </p>
+
+<p class="whs2"># add hexahedron</p>
+
+<p class="whs1">vId = mesh.AddVolume([nId2,
+ nId1, 38, 39, nId4, nId3, 245, 246])</p>
+
+<p class="whs1">if vId == 0: print
+ "KO Hexahedron addition."</p>
+
+<p class="whs1">else: print
+ "New Hexahedron has been added with ID ", vId </p>
+
+<p class="whs1"> </p>
+
+<h4>Add Polygon</h4>
+
+<p class="whs1"><span style="font-family: 'Lucida Console', monospace;">import
+ math</span></p>
+
+<p class="whs1">import salome</p>
+
+<p class="whs1"> </p>
+
+<p class="whs1">import smesh</p>
+
+<p class="whs1"> </p>
+
+<p class="whs2"># create an empty mesh
+ structure</p>
+
+<p class="whs1">mesh = smesh.Mesh()
+ </p>
+
+<p class="whs1"> </p>
+
+<p class="whs2"># a method to build
+ a polygonal mesh element with <nb_vert> angles:</p>
+
+<p class="whs1">def MakePolygon
+ (a_mesh, x0, y0, z0, radius, nb_vert):</p>
+
+<p class="whs1"> al
+ = 2.0 * math.pi / nb_vert</p>
+
+<p class="whs1"> node_ids
+ = []</p>
+
+<p class="whs1"> </p>
+
+<p class="whs2"> #
+ Create nodes for a polygon</p>
+
+<p class="whs1"> for
+ ii in range(nb_vert):</p>
+
+<p class="whs1"> nid
+ = mesh.AddNode(x0 + radius * math.cos(ii*al),</p>
+
+<p class="whs1"> y0
+ + radius * math.sin(ii*al),</p>
+
+<p class="whs1"> z0)</p>
+
+<p class="whs1"> node_ids.append(nid)</p>
+
+<p class="whs1"> pass</p>
+
+<p class="whs1"> </p>
+
+<p class="whs2"> #
+ Create a polygon</p>
+
+<p class="whs1"> return
+ mesh.AddPolygonalFace(node_ids)</p>
+
+<p class="whs1"> </p>
+
+<p class="whs2"># Create three polygons</p>
+
+<p class="whs1">f1 = MakePolygon(mesh,
+ 0, 0, 0,
+ 30, 13)</p>
+
+<p class="whs1">f2 = MakePolygon(mesh,
+ 0, 0, 10, 21, 9)</p>
+
+<p class="whs1">f3 = MakePolygon(mesh,
+ 0, 0, 20, 13, 6)</p>
+
+<p class="whs1"> </p>
+
+<p class="whs1">salome.sg.updateObjBrowser(1)
+ </p>
+
+<h4>Add Polyhedron</h4>
+
+<p class="whs1"><span style="font-family: 'Lucida Console', monospace;">import
+ salome</span></p>
+
+<p class="whs1">import math</p>
+
+<p class="whs1"> </p>
+
+<p class="whs2"># create an empty mesh
+ structure</p>
+
+<p class="whs1">mesh = smesh.Mesh()
+ </p>
+
+<p class="whs1"> </p>
+
+<p class="whs2"># Create nodes for
+ 12-hedron with pentagonal faces</p>
+
+<p class="whs1">al = 2 * math.pi
+ / 5.0</p>
+
+<p class="whs1">cosal = math.cos(al)</p>
+
+<p class="whs1">aa = 13</p>
+
+<p class="whs1">rr = aa / (2.0 *
+ math.sin(al/2.0))</p>
+
+<p class="whs1">dr = 2.0 * rr *
+ cosal</p>
+
+<p class="whs1">r1 = rr + dr</p>
+
+<p class="whs1">dh = rr * math.sqrt(2.0
+ * (1.0 - cosal * (1.0 + 2.0 * cosal)))</p>
+
+<p class="whs1">hh = 2.0 * dh -
+ dr * (rr*(cosal - 1) + (rr + dr)*(math.cos(al/2) - 1)) / dh</p>
+
+<p class="whs1"> </p>
+
+<p class="whs1">dd = [] <span style="font-family: 'Times New Roman', serif;">#
+ top</span></p>
+
+<p class="whs1">cc = [] <span style="font-family: 'Times New Roman', serif;">#
+ below top</span></p>
+
+<p class="whs1">bb = [] <span style="font-family: 'Times New Roman', serif;">#
+ above bottom</span></p>
+
+<p class="whs1">aa = [] <span style="font-family: 'Times New Roman', serif;">#
+ bottom</span></p>
+
+<p class="whs1"> </p>
+
+<p class="whs1">for i in range(5):</p>
+
+<p class="whs1"> cos_bot
+ = math.cos(i*al)</p>
+
+<p class="whs1"> sin_bot
+ = math.sin(i*al)</p>
+
+<p class="whs1"> </p>
+
+<p class="whs1"> <span
+ style="margin-top: 0pt;
+ margin-bottom: 0pt;
+ font-family: 'Lucida Console', monospace;">cos_top
+ = math.cos(i*al + al/2.0)</span></p>
+
+<p class="whs1"> sin_top
+ = math.sin(i*al + al/2.0)</p>
+
+<p class="whs1"> </p>
+
+<p class="whs1"> nd
+ = mesh.AddNode(rr * cos_top, rr * sin_top, hh )
+ <span style="font-family: 'Times New Roman', serif;"># top</span></p>
+
+<p class="whs1"> nc
+ = mesh.AddNode(r1 * cos_top, r1 * sin_top, hh - dh) <span style="font-family: 'Times New Roman', serif;">#
+ below top</span></p>
+
+<p class="whs1"> nb
+ = mesh.AddNode(r1 * cos_bot, r1 * sin_bot, dh)
+ <span style="font-family: 'Times New Roman', serif;"># above bottom</span></p>
+
+<p class="whs1"> na
+ = mesh.AddNode(rr * cos_bot, rr * sin_bot, 0)
+ <span style="font-family: 'Times New Roman', serif;"># bottom</span></p>
+
+<p class="whs1"> dd.append(nd)
+ <span style="font-family: 'Times New Roman', serif;"># top</span></p>
+
+<p class="whs1"> cc.append(nc)
+ <span style="font-family: 'Times New Roman', serif;"># below top</span></p>
+
+<p class="whs1"> bb.append(nb)
+ <span style="font-family: 'Times New Roman', serif;"># above bottom</span></p>
+
+<p class="whs1"> aa.append(na)
+ <span style="font-family: 'Times New Roman', serif;"># bottom</span></p>
+
+<p class="whs1"> </p>
+
+<p class="whs1"> pass</p>
+
+<p class="whs1"> </p>
+
+<p class="whs2"># Create a polyhedral
+ volume (12-hedron with pentagonal faces)</p>
+
+<p class="whs1">MeshEditor.AddPolyhedralVolume([dd[0],
+ dd[1], dd[2], dd[3], dd[4], <span style="font-family: 'Times New Roman', serif;">#
+ top</span></p>
+
+<p class="whs1"> dd[0],
+ cc[0], bb[1], cc[1], dd[1], #
+ -</p>
+
+<p class="whs1"> dd[1],
+ cc[1], bb[2], cc[2], dd[2], #
+ -</p>
+
+<p class="whs1"> dd[2],
+ cc[2], bb[3], cc[3], dd[3], <span style="font-family: 'Times New Roman', serif;">#
+ - below top</span></p>
+
+<p class="whs1"> dd[3],
+ cc[3], bb[4], cc[4], dd[4], #
+ -</p>
+
+<p class="whs1"> dd[4],
+ cc[4], bb[0], cc[0], dd[0], #
+ -</p>
+
+<p class="whs1"> aa[4],
+ bb[4], cc[4], bb[0], aa[0], #
+ .</p>
+
+<p class="whs1"> aa[3],
+ bb[3], cc[3], bb[4], aa[4], #
+ .</p>
+
+<p class="whs1"> aa[2],
+ bb[2], cc[2], bb[3], aa[3], <span style="font-family: 'Times New Roman', serif;">#
+ . above bottom</span></p>
+
+<p class="whs1"> aa[1],
+ bb[1], cc[1], bb[2], aa[2], #
+ .</p>
+
+<p class="whs1"> aa[0],
+ bb[0], cc[0], bb[1], aa[1], #
+ .</p>
+
+<p class="whs1"> aa[0],
+ aa[1], aa[2], aa[3], aa[4]], <span style="font-family: 'Times New Roman', serif;">#
+ bottom</span></p>
+
+<p class="whs1"> [5,5,5,5,5,5,5,5,5,5,5,5])</p>
+
+<p class="whs1">salome.sg.updateObjBrowser(1)
+ </p>
+
+<p class="whs1"> </p>
+
+<h3><a name=bookmark1>Removing Nodes and Elements</a></h3>
+
+<h4 class="whs4">Removing Nodes</h4>
+
+<p class="whs3"> </p>
+
+<p class="whs1"><span style="font-family: 'Lucida Console', monospace;">import
+ SMESH_mechanic</span></p>
+
+<p class="whs1"> </p>
+
+<p class="whs1">mesh = SMESH_mechanic.mesh</p>
+
+<p class="whs1"> </p>
+
+<p class="whs2"># remove nodes #246
+ and #255</p>
+
+<p class="whs1">res = mesh.RemoveNodes([246,
+ 255])</p>
+
+<p class="whs1">if res == 1: print
+ "Nodes removing is OK!"</p>
+
+<p class="whs1">else: print
+ "KO nodes removing." </p>
+
+<p class="whs3"> </p>
+
+<h4>Removing Elements</h4>
+
+<p class="whs1">import SMESH_mechanic</p>
+
+<p class="whs1"> </p>
+
+<p class="whs1">mesh = SMESH_mechanic.mesh</p>
+
+<p class="whs1"> </p>
+
+<p class="whs2"># remove three elements:
+ #850, #859 and #814</p>
+
+<p class="whs1">res = mesh.RemoveElements([850,
+ 859, 814])</p>
+
+<p class="whs1">if res == 1: print
+ "Elements removing is OK!"</p>
+
+<p class="whs1">else: print
+ "KO Elements removing." </p>
+
+<p class="whs1"> </p>
+
+<h3><a name=bookmark2>Renumbering Nodes and Elements</a></h3>
+
+<p class="whs5">import SMESH_mechanic</p>
+
+<p class="whs5"> </p>
+
+<p class="whs5">mesh = SMESH_mechanic.mesh</p>
+
+<p class="whs6"> </p>
+
+<p class="whs6"><span style="font-family: 'Lucida Console', monospace;">mesh.RenumberNodes()</span>
+ </p>
+
+<h3><a name=bookmark3>Moving Nodes</a></h3>
+
+<p class="whs1">import SMESH_mechanic</p>
+
+<p class="whs1"> </p>
+
+<p class="whs1">mesh = SMESH_mechanic.mesh</p>
+
+<p class="whs3"> </p>
+
+<p class="whs3"># move node #38</p>
+
+<p class="whs1">mesh.MoveNode(38,
+ 20., 10., 0.) </p>
+
+<p class="whs3"> </p>
+
+<h3><a name=bookmark4>Diagonal Inversion</a></h3>
+
+<p class="whs5">import salome</p>
+
+<p class="whs5">import smesh</p>
+
+<p class="whs5"> </p>
+
+<p class="whs7"># create an empty mesh
+ structure</p>
+
+<p class="whs5">mesh = smesh.Mesh()
+ </p>
+
+<p class="whs5"> </p>
+
+<p class="whs7"># create the following
+ mesh:</p>
+
+<p class="whs5"># .----.----.----.</p>
+
+<p class="whs5"># | /|
+ /|
+ /|</p>
+
+<p class="whs5"># | /
+ | / | /
+ |</p>
+
+<p class="whs5"># | / |
+ / | / |</p>
+
+<p class="whs5"># |/ |/
+ |/
+ |</p>
+
+<p class="whs5"># .----.----.----.</p>
+
+<p class="whs5"> </p>
+
+<p class="whs5">bb = [0, 0, 0, 0]</p>
+
+<p class="whs5">tt = [0, 0, 0, 0]</p>
+
+<p class="whs5">ff = [0, 0, 0, 0,
+ 0, 0]</p>
+
+<p class="whs5"> </p>
+
+<p class="whs5">bb[0] = mesh.AddNode(
+ 0., 0., 0.)</p>
+
+<p class="whs5">bb[1] = mesh.AddNode(10.,
+ 0., 0.)</p>
+
+<p class="whs5">bb[2] = mesh.AddNode(20.,
+ 0., 0.)</p>
+
+<p class="whs5">bb[3] = mesh.AddNode(30.,
+ 0., 0.)</p>
+
+<p class="whs5"> </p>
+
+<p class="whs5">tt[0] = mesh.AddNode(
+ 0., 15., 0.)</p>
+
+<p class="whs5">tt[1] = mesh.AddNode(10.,
+ 15., 0.)</p>
+
+<p class="whs5">tt[2] = mesh.AddNode(20.,
+ 15., 0.)</p>
+
+<p class="whs5">tt[3] = mesh.AddNode(30.,
+ 15., 0.)</p>
+
+<p class="whs5"> </p>
+
+<p class="whs5">ff[0] = mesh.AddFace([bb[0],
+ bb[1], tt[1]])</p>
+
+<p class="whs5">ff[1] = mesh.AddFace([bb[0],
+ tt[1], tt[0]])</p>
+
+<p class="whs5">ff[2] = mesh.AddFace([bb[1],
+ bb[2], tt[2]])</p>
+
+<p class="whs5">ff[3] = mesh.AddFace([bb[1],
+ tt[2], tt[1]])</p>
+
+<p class="whs5">ff[4] = mesh.AddFace([bb[2],
+ bb[3], tt[3]])</p>
+
+<p class="whs5">ff[5] = mesh.AddFace([bb[2],
+ tt[3], tt[2]])</p>
+
+<p class="whs5"> </p>
+
+<p class="whs7"># inverse the diagonal
+ bb[1] - tt[2]</p>
+
+<p class="whs5">print "\nDiagonal
+ inversion ... ",</p>
+
+<p class="whs5">res = mesh.InverseDiag(bb[1],
+ tt[2])</p>
+
+<p class="whs5">if not res: print
+ "failed!"</p>
+
+<p class="whs5">else: print
+ "done."</p>
+
+<p class="whs6"> </p>
+
+<p class="whs6"><span style="font-family: 'Lucida Console', monospace;">salome.sg.updateObjBrowser(1)</span>
+ </p>
+
+<h3><a name=bookmark5>Uniting two Triangles</a></h3>
+
+<p class="whs8">import salome</p>
+
+<p class="whs8">import smesh</p>
+
+<p> </p>
+
+<p># create an empty mesh structure</p>
+
+<p class="whs8">mesh = smesh.Mesh()
+ </p>
+
+<p> </p>
+
+<p># create the following mesh:</p>
+
+<p class="whs8"># .----.----.----.</p>
+
+<p class="whs8"># | /|
+ /|
+ /|</p>
+
+<p class="whs8"># | /
+ | / | /
+ |</p>
+
+<p class="whs8"># | / |
+ / | / |</p>
+
+<p class="whs8"># |/ |/
+ |/
+ |</p>
+
+<p class="whs8"># .----.----.----.</p>
+
+<p> </p>
+
+<p class="whs8">bb = [0, 0, 0, 0]</p>
+
+<p class="whs8">tt = [0, 0, 0, 0]</p>
+
+<p class="whs8">ff = [0, 0, 0, 0,
+ 0, 0]</p>
+
+<p class="whs8"> </p>
+
+<p class="whs8"><span style="font-family: 'Lucida Console', monospace;">bb[0]
+ = mesh.AddNode( 0., 0., 0.)</span></p>
+
+<p class="whs8">bb[1] = mesh.AddNode(10.,
+ 0., 0.)</p>
+
+<p class="whs8">bb[2] = mesh.AddNode(20.,
+ 0., 0.)</p>
+
+<p class="whs8">bb[3] = mesh.AddNode(30.,
+ 0., 0.)</p>
+
+<p class="whs8"> </p>
+
+<p class="whs8">tt[0] = mesh.AddNode(
+ 0., 15., 0.)</p>
+
+<p class="whs8">tt[1] = mesh.AddNode(10.,
+ 15., 0.)</p>
+
+<p class="whs8">tt[2] = mesh.AddNode(20.,
+ 15., 0.)</p>
+
+<p class="whs8">tt[3] = mesh.AddNode(30.,
+ 15., 0.)</p>
+
+<p class="whs8"> </p>
+
+<p class="whs8">ff[0] = mesh.AddFace([bb[0],
+ bb[1], tt[1]])</p>
+
+<p class="whs8">ff[1] = mesh.AddFace([bb[0],
+ tt[1], tt[0]])</p>
+
+<p class="whs8">ff[2] = mesh.AddFace([bb[1],
+ bb[2], tt[2]])</p>
+
+<p class="whs8">ff[3] = mesh.AddFace([bb[1],
+ tt[2], tt[1]])</p>
+
+<p class="whs8">ff[4] = mesh.AddFace([bb[2],
+ bb[3], tt[3]])</p>
+
+<p class="whs8">ff[5] = mesh.AddFace([bb[2],
+ tt[3], tt[2]]) </p>
+
+<p> </p>
+
+<p># delete the diagonal bb[1] - tt[2]</p>
+
+<p class="whs8">print "\nUnite
+ two triangles ... ",</p>
+
+<p class="whs8">res = mesh.DeleteDiag(bb[1],
+ tt[2])</p>
+
+<p class="whs8">if not res: print
+ "failed!"</p>
+
+<p class="whs8">else: print
+ "done."</p>
+
+<p class="whs8"> </p>
+
+<p><span style="font-family: 'Lucida Console', monospace;">salome.sg.updateObjBrowser(1)</span>
+ </p>
+
+<h3><a name=bookmark6>Uniting a Set of Triangles</a></h3>
+
+<p class="whs8">import salome</p>
+
+<p class="whs8">import smesh</p>
+
+<p class="whs8"> </p>
+
+<p class="whs9"># create an empty mesh
+ structure</p>
+
+<p class="whs8">mesh = smesh.Mesh()
+ </p>
+
+<p class="whs8"> </p>
+
+<p class="whs9"># create the following
+ mesh:</p>
+
+<p class="whs8"># .----.----.----.</p>
+
+<p class="whs8"># | /|
+ /|
+ /|</p>
+
+<p class="whs8"># | /
+ | / | /
+ |</p>
+
+<p class="whs8"># | / |
+ / | / |</p>
+
+<p class="whs8"># |/ |/
+ |/
+ |</p>
+
+<p class="whs8"># .----.----.----.</p>
+
+<p class="whs8"> </p>
+
+<p class="whs8">bb = [0, 0, 0, 0]</p>
+
+<p class="whs8">tt = [0, 0, 0, 0]</p>
+
+<p class="whs8">ff = [0, 0, 0, 0,
+ 0, 0]</p>
+
+<p class="whs8"> </p>
+
+<p class="whs8"><span style="font-family: 'Lucida Console', monospace;">bb[0]
+ = mesh.AddNode( 0., 0., 0.)</span></p>
+
+<p class="whs8">bb[1] = mesh.AddNode(10.,
+ 0., 0.)</p>
+
+<p class="whs8">bb[2] = mesh.AddNode(20.,
+ 0., 0.)</p>
+
+<p class="whs8">bb[3] = mesh.AddNode(30.,
+ 0., 0.)</p>
+
+<p class="whs8"> </p>
+
+<p class="whs8">tt[0] = mesh.AddNode(
+ 0., 15., 0.)</p>
+
+<p class="whs8">tt[1] = mesh.AddNode(10.,
+ 15., 0.)</p>
+
+<p class="whs8">tt[2] = mesh.AddNode(20.,
+ 15., 0.)</p>
+
+<p class="whs8">tt[3] = mesh.AddNode(30.,
+ 15., 0.)</p>
+
+<p class="whs8"> </p>
+
+<p class="whs8">ff[0] = mesh.AddFace([bb[0],
+ bb[1], tt[1]])</p>
+
+<p class="whs8">ff[1] = mesh.AddFace([bb[0],
+ tt[1], tt[0]])</p>
+
+<p class="whs8">ff[2] = mesh.AddFace([bb[1],
+ bb[2], tt[2]])</p>
+
+<p class="whs8">ff[3] = mesh.AddFace([bb[1],
+ tt[2], tt[1]])</p>
+
+<p class="whs8">ff[4] = mesh.AddFace([bb[2],
+ bb[3], tt[3]])</p>
+
+<p class="whs8">ff[5] = mesh.AddFace([bb[2],
+ tt[3], tt[2]])</p>
+
+<p class="whs8"> </p>
+
+<p class="whs9"># unite a set of triangles</p>
+
+<p class="whs8"><span style="font-family: 'Lucida Console', monospace;">print
+ "\nUnite a set of triangles ... ",</span></p>
+
+<p class="whs8">res = mesh.TriToQuad([ff[2],
+ ff[3], ff[4], ff[5]], smesh.FT_MinimumAngle, 60.)</p>
+
+<p class="whs8">if not res: print
+ "failed!"</p>
+
+<p class="whs8">else: print
+ "done."</p>
+
+<p class="whs8"> </p>
+
+<p class="whs8">salome.sg.updateObjBrowser(1)
+ </p>
+
+<h3><a name=bookmark12>Orientation</a></h3>
+
+<p class="whs8">import salome</p>
+
+<p class="whs8">import smesh</p>
+
+<p class="whs8"> </p>
+
+<p># create an empty mesh structure</p>
+
+<p class="whs8">mesh = smesh.Mesh()
+ </p>
+
+<p class="whs8"> </p>
+
+<p class="whs9"># build five quadrangles:</p>
+
+<p class="whs8">dx = 10</p>
+
+<p class="whs8">dy = 20</p>
+
+<p class="whs8"> </p>
+
+<p class="whs8"><span style="font-family: 'Lucida Console', monospace;">n1
+ = mesh.AddNode(0.0
+ * dx, 0, 0)</span></p>
+
+<p class="whs8">n2 =
+ mesh.AddNode(1.0 * dx, 0, 0)</p>
+
+<p class="whs8">n3 =
+ mesh.AddNode(2.0 * dx, 0, 0)</p>
+
+<p class="whs8">n4 =
+ mesh.AddNode(3.0 * dx, 0, 0)</p>
+
+<p class="whs8">n5 =
+ mesh.AddNode(4.0 * dx, 0, 0)</p>
+
+<p class="whs8">n6 =
+ mesh.AddNode(5.0 * dx, 0, 0)</p>
+
+<p class="whs8">n7 =
+ mesh.AddNode(0.0 * dx, dy, 0)</p>
+
+<p class="whs8">n8 =
+ mesh.AddNode(1.0 * dx, dy, 0)</p>
+
+<p class="whs8">n9 =
+ mesh.AddNode(2.0 * dx, dy, 0)</p>
+
+<p class="whs8">n10 = mesh.AddNode(3.0
+ * dx, dy, 0)</p>
+
+<p class="whs8">n11 = mesh.AddNode(4.0
+ * dx, dy, 0)</p>
+
+<p class="whs8">n12 = mesh.AddNode(5.0
+ * dx, dy, 0)</p>
+
+<p class="whs8"> </p>
+
+<p class="whs8">f1 = mesh.AddFace([n1,
+ n2, n8 , n7 ])</p>
+
+<p class="whs8">f2 = mesh.AddFace([n2,
+ n3, n9 , n8 ])</p>
+
+<p class="whs8">f3 = mesh.AddFace([n3,
+ n4, n10, n9 ])</p>
+
+<p class="whs8">f4 = mesh.AddFace([n4,
+ n5, n11, n10])</p>
+
+<p class="whs8">f5 = mesh.AddFace([n5,
+ n6, n12, n11]) </p>
+
+<p class="whs8"> </p>
+
+<p class="whs9"># Change the orientation
+ of the second and the fourth faces.</p>
+
+<p class="whs8">mesh.Reorient([2,
+ 4])</p>
+
+<p class="whs8"> </p>
+
+<p class="whs8">salome.sg.updateObjBrowser(1)
+ </p>
+
+<p class="whs8"> </p>
+
+<h3><a name=bookmark7>Cutting Quadrangles</a></h3>
+
+<p class="whs8">import SMESH_mechanic</p>
+
+<p class="whs8"> </p>
+
+<p class="whs8">smesh = SMESH_mechanic.smesh</p>
+
+<p class="whs8">mesh =
+ SMESH_mechanic.mesh</p>
+
+<p class="whs8"> </p>
+
+<p class="whs9"># cut two quadrangles:
+ 405 and 406</p>
+
+<p class="whs8">mesh.QuadToTri([405,
+ 406], smesh.FT_MinimumAngle) </p>
+
+<p class="whs9"> </p>
+
+<h3><a name=bookmark8>Smoothing</a></h3>
+
+<p class="whs8">import salome</p>
+
+<p class="whs8">import geompy</p>
+
+<p class="whs8"> </p>
+
+<p class="whs8">import SMESH_mechanic</p>
+
+<p class="whs8"> </p>
+
+<p class="whs8"><span style="font-family: 'Lucida Console', monospace;">smesh
+ = SMESH_mechanic.smesh</span></p>
+
+<p class="whs8">mesh = SMESH_mechanic.mesh</p>
+
+<p class="whs8"> </p>
+
+<p class="whs9"># select the top face</p>
+
+<p class="whs8"><span style="font-family: 'Lucida Console', monospace;">faces
+ = geompy.SubShapeAllSorted(SMESH_mechanic.shape_mesh, geompy.ShapeType["FACE"])</span></p>
+
+<p class="whs8">face = faces[3]</p>
+
+<p class="whs8">geompy.addToStudyInFather(SMESH_mechanic.shape_mesh,
+ face, "face planar with hole")</p>
+
+<p class="whs8"> </p>
+
+<p class="whs8"># create a group of
+ faces to be smoothed</p>
+
+<p class="whs8">GroupSmooth = mesh.GroupOnGeom(face,
+ "Group of faces (smooth)", smesh.FACE)</p>
+
+<p class="whs8"> </p>
+
+<p class="whs8"># perform smoothing</p>
+
+<p class="whs8"> </p>
+
+<p class="whs8"># boolean SmoothObject(Object,
+ IDsOfFixedNodes, MaxNbOfIterations, MaxAspectRatio, Method)</p>
+
+<p class="whs8">res = mesh.SmoothObject(GroupSmooth,
+ [], 20, 2., smesh.CENTROIDAL_SMOOTH)</p>
+
+<p class="whs8">print "\nSmoothing
+ ... ",</p>
+
+<p class="whs8">if not res: print
+ "failed!"</p>
+
+<p class="whs8">else: print
+ "done."</p>
+
+<p class="whs8"> </p>
+
+<p class="whs8">salome.sg.updateObjBrowser(1)
+ </p>
+
+<p class="whs8"> </p>
+
+<h3><a name=bookmark9>Extrusion</a></h3>
+
+<p class="whs8">import salome</p>
+
+<p class="whs8">import geompy</p>
+
+<p class="whs8"> </p>
+
+<p class="whs8">import SMESH_mechanic</p>
+
+<p class="whs8"> </p>
+
+<p class="whs8"><span style="font-family: 'Lucida Console', monospace;">smesh
+ = SMESH_mechanic.smesh</span></p>
+
+<p class="whs8">mesh = SMESH_mechanic.mesh
+ </p>
+
+<p class="whs8"> </p>
+
+<p class="whs9"># select the top face</p>
+
+<p class="whs8">faces = geompy.SubShapeAllSorted(SMESH_mechanic.shape_mesh,
+ geompy.ShapeType["FACE"])</p>
+
+<p class="whs8">face = faces[7]</p>
+
+<p class="whs8">geompy.addToStudyInFather(SMESH_mechanic.shape_mesh,
+ face, "face circular top")</p>
+
+<p class="whs8"> </p>
+
+<p class="whs8"># create a vector
+ for extrusion</p>
+
+<p class="whs8">point = smesh.PointStruct(0.,
+ 0., 5.)</p>
+
+<p class="whs8">vector = smesh.DirStruct(point)</p>
+
+<p class="whs8"> </p>
+
+<p class="whs8"># create a group to
+ be extruded</p>
+
+<p class="whs8">GroupTri = mesh.GroupOnGeom(face,
+ "Group of faces (extrusion)", smesh.FACE)</p>
+
+<p class="whs8"> </p>
+
+<p class="whs8"># perform extrusion
+ of the group</p>
+
+<p class="whs8">mesh.ExtrusionSweepObject(GroupTri,
+ vector, 5)</p>
+
+<p class="whs8"> </p>
+
+<p class="whs9"><span style="font-family: 'Lucida Console', monospace;">salome.sg.updateObjBrowser(1)</span>
+ </p>
+
+<h3><a name=bookmark10>Extrusion along a Path</a></h3>
+
+<p class="whs8">import math</p>
+
+<p class="whs8">import salome</p>
+
+<p class="whs8"> </p>
+
+<p class="whs9"># Geometry</p>
+
+<p class="whs8">import geompy</p>
+
+<p class="whs9"> </p>
+
+<p class="whs9"># 1. Create points</p>
+
+<p class="whs8">points = [[0, 0],
+ [50, 30], [50, 110], [0, 150], [-80, 150], [-130, 70], [-130, -20]]</p>
+
+<p class="whs8"> </p>
+
+<p class="whs8">iv = 1</p>
+
+<p class="whs8">vertices = []</p>
+
+<p class="whs8">for point in points:</p>
+
+<p class="whs10">vert
+ = geompy.MakeVertex(point[0], point[1], 0)</p>
+
+<p class="whs10">geompy.addToStudy(vert,
+ "Vertex_" + `iv`)</p>
+
+<p class="whs10">vertices.append(vert)</p>
+
+<p class="whs10">iv
+ += 1</p>
+
+<p class="whs10">pass</p>
+
+<p class="whs9"> </p>
+
+<p class="whs9"># 2. Create edges and
+ wires</p>
+
+<p class="whs8">Edge_straight = geompy.MakeEdge(vertices[0],
+ vertices[4])</p>
+
+<p class="whs8">Edge_bezierrr = geompy.MakeBezier(vertices)</p>
+
+<p class="whs8">Wire_polyline = geompy.MakePolyline(vertices)</p>
+
+<p class="whs8">Edge_Circle =
+ geompy.MakeCircleThreePnt(vertices[0], vertices[1], vertices[2])</p>
+
+<p class="whs8"> </p>
+
+<p class="whs8">geompy.addToStudy(Edge_straight,
+ "Edge_straight")</p>
+
+<p class="whs8">geompy.addToStudy(Edge_bezierrr,
+ "Edge_bezierrr")</p>
+
+<p class="whs8">geompy.addToStudy(Wire_polyline,
+ "Wire_polyline")</p>
+
+<p class="whs8">geompy.addToStudy(Edge_Circle
+ , "Edge_Circle")</p>
+
+<p class="whs9"> </p>
+
+<p class="whs9"># 3. Explode wire on
+ edges, as they will be used for mesh extrusion</p>
+
+<p class="whs8">Wire_polyline_edges
+ = geompy.SubShapeAll(Wire_polyline, geompy.ShapeType["EDGE"])</p>
+
+<p class="whs8">for ii in range(len(Wire_polyline_edges)):</p>
+
+<p class="whs10">geompy.addToStudyInFather(Wire_polyline,
+ Wire_polyline_edges[ii], "Edge_" + `ii + 1`)</p>
+
+<p class="whs10">pass</p>
+
+<p class="whs11"> </p>
+
+<p class="whs9"># Mesh</p>
+
+<p class="whs8">import smesh</p>
+
+<p class="whs8"> </p>
+
+<p class="whs12"># Mesh
+ the given shape with the given 1d hypothesis</p>
+
+<p class="whs13">def
+ Mesh1D(shape1d, nbSeg, name):</p>
+
+<p class="whs13"> mesh1d_tool
+ = smesh.Mesh(shape1d, name)</p>
+
+<p class="whs13"> algo
+ = mesh1d_tool.Segment()</p>
+
+<p class="whs13"> hyp
+ = algo.NumberOfSegments(nbSeg)</p>
+
+<p class="whs13"> isDone
+ = mesh1d_tool.Compute()</p>
+
+<p class="whs13"> if
+ not isDone: print 'Mesh ', name, ': computation failed'</p>
+
+<p class="whs13"> return
+ mesh1d_tool</p>
+
+<p class="whs12"># Create
+ a mesh with six nodes, seven edges and two quadrangle faces</p>
+
+<p class="whs13">def
+ MakeQuadMesh2(mesh_name):</p>
+
+<p class="whs13"> quad_1
+ = smesh.Mesh(name = mesh_name)</p>
+
+<p class="whs13"> </p>
+
+<p class="whs13"> <span
+ style="font-family: 'Times New Roman', serif;">#
+ six nodes</span></p>
+
+<p class="whs13"> n1
+ = quad_1.AddNode(0, 20, 10)</p>
+
+<p class="whs13"> n2
+ = quad_1.AddNode(0, 40, 10)</p>
+
+<p class="whs13"> n3
+ = quad_1.AddNode(0, 40, 30)</p>
+
+<p class="whs13"> n4
+ = quad_1.AddNode(0, 20, 30)</p>
+
+<p class="whs13"> n5
+ = quad_1.AddNode(0, 0,
+ 30)</p>
+
+<p class="whs13"> n6
+ = quad_1.AddNode(0, 0,
+ 10)</p>
+
+<p class="whs13"> </p>
+
+<p class="whs13"> <span
+ style="font-family: 'Times New Roman', serif;">#
+ seven edges</span></p>
+
+<p class="whs13"> quad_1.AddEdge([n1,
+ n2]) # 1</p>
+
+<p class="whs13"> quad_1.AddEdge([n2,
+ n3]) # 2</p>
+
+<p class="whs13"> quad_1.AddEdge([n3,
+ n4]) # 3</p>
+
+<p class="whs13"> quad_1.AddEdge([n4,
+ n1]) # 4</p>
+
+<p class="whs13"> quad_1.AddEdge([n4,
+ n5]) # 5</p>
+
+<p class="whs13"> quad_1.AddEdge([n5,
+ n6]) # 6</p>
+
+<p class="whs13"> quad_1.AddEdge([n6,
+ n1]) # 7</p>
+
+<p class="whs13"> </p>
+
+<p class="whs13"> <span
+ style="font-family: 'Times New Roman', serif;">#
+ two quadrangle faces</span></p>
+
+<p class="whs13"> quad_1.AddFace([n1,
+ n2, n3, n4]) # 8</p>
+
+<p class="whs13"> quad_1.AddFace([n1,
+ n4, n5, n6]) # 9</p>
+
+<p class="whs13"> return
+ [quad_1, [1,2,3,4,5,6,7], [8,9]]</p>
+
+<p class="whs13"> </p>
+
+<p class="whs12"># Path
+ meshes</p>
+
+<p class="whs13">Edge_straight_mesh
+ = Mesh1D(Edge_straight, 7, "Edge_straight")</p>
+
+<p class="whs13">Edge_bezierrr_mesh
+ = Mesh1D(Edge_bezierrr, 7, "Edge_bezierrr")</p>
+
+<p class="whs13">Wire_polyline_mesh
+ = Mesh1D(Wire_polyline, 3, "Wire_polyline")</p>
+
+<p class="whs13">Edge_Circle_mesh
+ =
+ Mesh1D(Edge_Circle ,
+ 8, "Edge_Circle")</p>
+
+<p class="whs13"> </p>
+
+<p class="whs12"># Initial
+ meshes (to be extruded)</p>
+
+<p class="whs13">[quad_1,
+ ee_1, ff_1] = MakeQuadMesh2("quad_1")</p>
+
+<p class="whs13">[quad_2,
+ ee_2, ff_2] = MakeQuadMesh2("quad_2")</p>
+
+<p class="whs13">[quad_3,
+ ee_3, ff_3] = MakeQuadMesh2("quad_3")</p>
+
+<p class="whs13">[quad_4,
+ ee_4, ff_4] = MakeQuadMesh2("quad_4")</p>
+
+<p class="whs13">[quad_5,
+ ee_5, ff_5] = MakeQuadMesh2("quad_5")</p>
+
+<p class="whs13">[quad_6,
+ ee_6, ff_6] = MakeQuadMesh2("quad_6")</p>
+
+<p class="whs13">[quad_7,
+ ee_7, ff_7] = MakeQuadMesh2("quad_7")</p>
+
+<p class="whs13"> </p>
+
+<p class="whs12"># ExtrusionAlongPath</p>
+
+<p class="whs12"># IDsOfElements,
+ PathMesh, PathShape, NodeStart,</p>
+
+<p class="whs12"># HasAngles,
+ Angles, HasRefPoint, RefPoint</p>
+
+<p class="whs13">refPoint
+ = smesh.PointStruct(0, 0, 0)</p>
+
+<p class="whs13">a10
+ = 10.0*math.pi/180.0</p>
+
+<p class="whs13">a45
+ = 45.0*math.pi/180.0</p>
+
+<p class="whs13"> </p>
+
+<p class="whs12"># 1.
+ Extrusion of two mesh edges along a straight path</p>
+
+<p class="whs13">error
+ = quad_1.ExtrusionAlongPath([1,2], Edge_straight_mesh, Edge_straight,
+ 1,</p>
+
+<p class="whs13"> 0,
+ [], 0, refPoint)</p>
+
+<p class="whs13"> </p>
+
+<p class="whs12"># 2.
+ Extrusion of one mesh edge along a curved path</p>
+
+<p class="whs13">error
+ = quad_2.ExtrusionAlongPath([2], Edge_bezierrr_mesh, Edge_bezierrr, 1,</p>
+
+<p class="whs13"> 0,
+ [], 0, refPoint)</p>
+
+<p class="whs13"> </p>
+
+<p class="whs12"># 3.
+ Extrusion of one mesh edge along a curved path with usage of angles</p>
+
+<p class="whs13">error
+ = quad_3.ExtrusionAlongPath([2], Edge_bezierrr_mesh, Edge_bezierrr, 1,</p>
+
+<p class="whs13"> 1,
+ [a45, a45, a45, 0, -a45, -a45, -a45], 0, refPoint)</p>
+
+<p class="whs13"> </p>
+
+<p class="whs12"># 4.
+ Extrusion of one mesh edge along the path, which is a part of a meshed
+ wire</p>
+
+<p class="whs13">error
+ = quad_4.ExtrusionAlongPath([4], Wire_polyline_mesh, Wire_polyline_edges[0],
+ 1,</p>
+
+<p class="whs13"> 1,
+ [a10, a10, a10], 0, refPoint)</p>
+
+<p class="whs13"> </p>
+
+<p class="whs12"># 5.
+ Extrusion of two mesh faces along the path, which is a part of a meshed
+ wire</p>
+
+<p class="whs13">error
+ = quad_5.ExtrusionAlongPath(ff_5 , Wire_polyline_mesh, Wire_polyline_edges[2],
+ 4,</p>
+
+<p class="whs13"> 0,
+ [], 0, refPoint)</p>
+
+<p class="whs13"> </p>
+
+<p class="whs12"># 6.
+ Extrusion of two mesh faces along a closed path</p>
+
+<p class="whs13">error
+ = quad_6.ExtrusionAlongPath(ff_6 , Edge_Circle_mesh, Edge_Circle, 1,</p>
+
+<p class="whs13"> 0,
+ [], 0, refPoint)</p>
+
+<p class="whs13"> </p>
+
+<p class="whs12"># 7.
+ Extrusion of two mesh faces along a closed path with usage of angles</p>
+
+<p class="whs13">error
+ = quad_7.ExtrusionAlongPath(ff_7, Edge_Circle_mesh, Edge_Circle, 1,</p>
+
+<p class="whs13"> 1,
+ [a45, -a45, a45, -a45, a45, -a45, a45, -a45], 0, refPoint)</p>
+
+<p class="whs13"> </p>
+
+<p class="whs13">salome.sg.updateObjBrowser(1)
+ </p>
+
+<h3><a name=bookmark11>Revolution</a></h3>
+
+<p class="whs8">import math</p>
+
+<p class="whs8"> </p>
+
+<p class="whs8">import SMESH_mechanic</p>
+
+<p class="whs8"> </p>
+
+<p class="whs8"><span style="font-family: 'Lucida Console', monospace;">mesh
+ = SMESH_mechanic.mesh</span></p>
+
+<p class="whs8">smesh = SMESH_mechanic.smesh</p>
+
+<p class="whs8"> </p>
+
+<p class="whs9"># create a group of faces
+ to be revolved</p>
+
+<p class="whs8">FacesRotate = [492,
+ 493, 502, 503]</p>
+
+<p class="whs8">GroupRotate = mesh.CreateGroup(SMESH.FACE,"Group
+ of faces (rotate)")</p>
+
+<p class="whs8">GroupRotate.Add(FacesRotate)</p>
+
+<p class="whs9"> </p>
+
+<p class="whs9"># define revolution angle
+ and axis</p>
+
+<p class="whs8">angle45 = 45 * math.pi
+ / 180</p>
+
+<p class="whs8">axisXYZ = SMESH.AxisStruct(-38.3128,
+ -73.3658, -23.321, -13.3402, -13.3265, 6.66632)</p>
+
+<p class="whs8"> </p>
+
+<p class="whs9"># perform revolution
+ of an object</p>
+
+<p class="whs8">mesh.RotationSweepObject(GroupRotate,
+ axisXYZ, angle45, 4, 1e-5) </p>
+
+<p class="whs9"> </p>
+
+<h3><a name=bookmark13>Pattern Mapping</a></h3>
+
+<p class="whs8">import geompy</p>
+
+<p class="whs8"> </p>
+
+<p class="whs8">import smesh</p>
+
+<p class="whs8"> </p>
+
+<p class="whs8"># define the geometry</p>
+
+<p class="whs9">Box_1 = geompy.MakeBoxDXDYDZ(200.,
+ 200., 200.)</p>
+
+<p class="whs9">geompy.addToStudy(Box_1,
+ "Box_1")</p>
+
+<p class="whs9"> </p>
+
+<p class="whs9">faces = geompy.SubShapeAll(Box_1,
+ geompy.ShapeType["FACE"])</p>
+
+<p class="whs9">Face_1 = faces[0]</p>
+
+<p class="whs9">Face_2 = faces[1]</p>
+
+<p class="whs9"> </p>
+
+<p class="whs9">geompy.addToStudyInFather(Box_1,
+ Face_1, "Face_1")</p>
+
+<p class="whs9">geompy.addToStudyInFather(Box_1,
+ Face_2, "Face_2")</p>
+
+<p class="whs9"> </p>
+
+<p class="whs8"># build a quadrangle
+ mesh 3x3 on Face_1</p>
+
+<p class="whs9">Mesh_1 = smesh.Mesh(Face_1)</p>
+
+<p class="whs9">algo1D = Mesh_1.Segment()</p>
+
+<p class="whs9">algo1D.NumberOfSegments(3)</p>
+
+<p class="whs9">Mesh_1.Quadrangle()</p>
+
+<p class="whs9"> </p>
+
+<p class="whs9">isDone = Mesh_1.Compute()</p>
+
+<p class="whs9">if not isDone: print
+ 'Mesh Mesh_1 : computation failed'</p>
+
+<p class="whs9"> </p>
+
+<p class="whs8"># build a triangle
+ mesh on Face_2</p>
+
+<p class="whs9">Mesh_2 = smesh.Mesh(Face_2)</p>
+
+<p class="whs9"> </p>
+
+<p class="whs9">algo1D = Mesh_2.Segment()</p>
+
+<p class="whs9">algo1D.NumberOfSegments(1)</p>
+
+<p class="whs9">algo2D = Mesh_2.Triangle()</p>
+
+<p class="whs9">algo2D.MaxElementArea(240)</p>
+
+<p class="whs9"> </p>
+
+<p class="whs9">isDone = Mesh_2.Compute()</p>
+
+<p class="whs9">if not isDone: print
+ 'Mesh Mesh_2 : computation failed'</p>
+
+<p class="whs9"> </p>
+
+<p class="whs8"># create a pattern</p>
+
+<p class="whs9">pattern = smesh.GetPattern()</p>
+
+<p class="whs9"> </p>
+
+<p class="whs9">isDone = pattern.LoadFromFace(Mesh_2.GetMesh(),
+ Face_2, 0)</p>
+
+<p class="whs9">if (isDone != 1): print
+ 'LoadFromFace :', pattern.GetErrorCode()</p>
+
+<p class="whs9"> </p>
+
+<p class="whs8"># apply the pattern
+ to a face of the first mesh</p>
+
+<p class="whs9">pattern.ApplyToMeshFaces(Mesh_1.GetMesh(),
+ [17], 0, 0)</p>
+
+<p class="whs9"> </p>
+
+<p class="whs9">isDone = pattern.MakeMesh(Mesh_1.GetMesh(),
+ 0, 0)</p>
+
+<p class="whs9">if (isDone != 1): print
+ 'MakeMesh :', pattern.GetErrorCode() </p>
+
+<script type="text/javascript" language="javascript1.2">
+<!--
+if (window.writeIntopicBar)
+ writeIntopicBar(0);
+//-->
+</script>
+</body>
+</html>