mirror of
https://git.salome-platform.org/gitpub/modules/smesh.git
synced 2024-12-30 19:30:35 +05:00
2023 lines
52 KiB
HTML
Executable File
2023 lines
52 KiB
HTML
Executable File
<!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-weight:bold; color:#ff0000; }
|
|
p.whs2 { margin-top:0pt; margin-bottom:0pt; }
|
|
p.whs3 { margin-top:0pt; margin-bottom:0pt; font-family:'Lucida Console' , monospace; }
|
|
p.whs4 { margin-top:0pt; margin-bottom:0pt; font-family:'Times New Roman' , serif; }
|
|
h4.whs5 { margin-top:0pt; margin-bottom:0pt; }
|
|
p.whs6 { margin-top:0.5pt; margin-bottom:0pt; font-family:'Lucida Console' , monospace; }
|
|
p.whs7 { margin-top:0.5pt; margin-bottom:0pt; }
|
|
p.whs8 { margin-top:0.5pt; margin-bottom:0pt; font-family:'Times New Roman' , serif; }
|
|
p.whs9 { font-family:'Lucida Console' , monospace; }
|
|
p.whs10 { font-family:'Times New Roman' , serif; }
|
|
p.whs11 { margin-left:36px; font-family:'Lucida Console' , monospace; }
|
|
p.whs12 { margin-left:36px; }
|
|
p.whs13 { margin-left:192px; font-family:'Lucida Console' , monospace; }
|
|
p.whs14 { 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 += "p.whs4 {margin-top:1pt;margin-bottom:1pt; }";
|
|
strNSS += "h4.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.whs8 {margin-top:1pt;margin-bottom:1pt; }";
|
|
strNSS += "p.whs14 {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>
|
|
|
|
<p class="whs1"># Attention! The scripts for Adding nodes
|
|
and Elements have
|
|
been written using the old approach basing on direct usage of SMESH idl
|
|
interface.</p>
|
|
|
|
<p class="whs1"># For the moment smesh package doesn't provide
|
|
methods to add nodes and elements. </p>
|
|
|
|
<p class="whs2"><span style="font-weight: bold; color: #ff0000;"><font color=red ><B>#
|
|
In the next SALOME version the scripts will be updated to use only the
|
|
commands from smesh package.</B></font></span></p>
|
|
|
|
<h4>Add Node</h4>
|
|
|
|
<p class="whs3"><span style="font-family: 'Lucida Console', monospace;">import
|
|
SMESH</span></p>
|
|
|
|
<p class="whs3">import SMESH_mechanic</p>
|
|
|
|
<p class="whs3"> </p>
|
|
|
|
<p class="whs3">mesh = SMESH_mechanic.mesh</p>
|
|
|
|
<p class="whs3">aMeshEditor = mesh.GetMeshEditor()</p>
|
|
|
|
<p class="whs3"> </p>
|
|
|
|
<p class="whs4"># add node</p>
|
|
|
|
<p class="whs3">new_id = aMeshEditor.AddNode(50,
|
|
10, 0)</p>
|
|
|
|
<p class="whs3">print ""</p>
|
|
|
|
<p class="whs3">if new_id == 0:
|
|
print "KO node addition."</p>
|
|
|
|
<p class="whs3">else: print
|
|
"New Node has been added with ID ", new_id </p>
|
|
|
|
<p class="whs3"> </p>
|
|
|
|
<h4>Add Edge</h4>
|
|
|
|
<p class="whs3"><span style="font-family: 'Lucida Console', monospace;">import
|
|
SMESH</span></p>
|
|
|
|
<p class="whs3">import SMESH_mechanic</p>
|
|
|
|
<p class="whs3"> </p>
|
|
|
|
<p class="whs3">mesh = SMESH_mechanic.mesh</p>
|
|
|
|
<p class="whs3">aMeshEditor = mesh.GetMeshEditor()</p>
|
|
|
|
<p class="whs3">print ""</p>
|
|
|
|
<p class="whs3"> </p>
|
|
|
|
<p class="whs4"># add node</p>
|
|
|
|
<p class="whs3">n1 = aMeshEditor.AddNode(50,
|
|
10, 0)</p>
|
|
|
|
<p class="whs3">if n1 == 0: print
|
|
"KO node addition."</p>
|
|
|
|
<p class="whs3"> </p>
|
|
|
|
<p class="whs4"># add edge</p>
|
|
|
|
<p class="whs3">e1 = aMeshEditor.AddEdge([n1,
|
|
38])</p>
|
|
|
|
<p class="whs3">if e1 == 0: print
|
|
"KO edge addition."</p>
|
|
|
|
<p class="whs3">else: print
|
|
"New Edge has been added with ID ", e1 </p>
|
|
|
|
<p class="whs3"> </p>
|
|
|
|
<h4>Add Triangle</h4>
|
|
|
|
<p class="whs3"><span style="font-family: 'Lucida Console', monospace;">import
|
|
SMESH</span></p>
|
|
|
|
<p class="whs3">import SMESH_mechanic</p>
|
|
|
|
<p class="whs3"> </p>
|
|
|
|
<p class="whs3">mesh = SMESH_mechanic.mesh</p>
|
|
|
|
<p class="whs3">aMeshEditor = mesh.GetMeshEditor()</p>
|
|
|
|
<p class="whs3">print ""</p>
|
|
|
|
<p class="whs3"> </p>
|
|
|
|
<p class="whs4"># add node</p>
|
|
|
|
<p class="whs3">n1 = aMeshEditor.AddNode(50,
|
|
10, 0)</p>
|
|
|
|
<p class="whs3">if n1 == 0: print
|
|
"KO node addition."</p>
|
|
|
|
<p class="whs3"> </p>
|
|
|
|
<p class="whs4"># add triangle</p>
|
|
|
|
<p class="whs3">t1 = aMeshEditor.AddFace([n1,
|
|
38, 39])</p>
|
|
|
|
<p class="whs3">if t1 == 0: print
|
|
"KO triangle addition."</p>
|
|
|
|
<p class="whs3">else: print
|
|
"New Triangle has been added with ID ", t1 </p>
|
|
|
|
<h4>Add Quadrangle</h4>
|
|
|
|
<p class="whs3">import SMESH</p>
|
|
|
|
<p class="whs3">import SMESH_mechanic</p>
|
|
|
|
<p class="whs3"> </p>
|
|
|
|
<p class="whs3">mesh = SMESH_mechanic.mesh</p>
|
|
|
|
<p class="whs3">aMeshEditor = mesh.GetMeshEditor()</p>
|
|
|
|
<p class="whs3">print ""</p>
|
|
|
|
<p class="whs2"> </p>
|
|
|
|
<p class="whs2"># add node</p>
|
|
|
|
<p class="whs3">n1 = aMeshEditor.AddNode(50,
|
|
10, 0)</p>
|
|
|
|
<p class="whs3">if n1 == 0: print
|
|
"KO node addition."</p>
|
|
|
|
<p class="whs3"> </p>
|
|
|
|
<p class="whs3">n2 = aMeshEditor.AddNode(40,
|
|
20, 0)</p>
|
|
|
|
<p class="whs3">if n2 == 0: print
|
|
"KO node addition."</p>
|
|
|
|
<p class="whs3"> </p>
|
|
|
|
<p class="whs2"># add quadrangle</p>
|
|
|
|
<p class="whs3">q1 = aMeshEditor.AddFace([n2,
|
|
n1, 38, 39])</p>
|
|
|
|
<p class="whs3">if q1 == 0: print
|
|
"KO quadrangle addition."</p>
|
|
|
|
<p class="whs2"><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="whs3">import SMESH </p>
|
|
|
|
<p class="whs3">import SMESH_mechanic</p>
|
|
|
|
<p class="whs3"> </p>
|
|
|
|
<p class="whs3">mesh = SMESH_mechanic.mesh</p>
|
|
|
|
<p class="whs3">aMeshEditor = mesh.GetMeshEditor()</p>
|
|
|
|
<p class="whs3">print ""</p>
|
|
|
|
<p class="whs3"> </p>
|
|
|
|
<p class="whs4"># add node</p>
|
|
|
|
<p class="whs3">n1 = aMeshEditor.AddNode(50,
|
|
10, 0)</p>
|
|
|
|
<p class="whs3">if n1 == 0: print
|
|
"KO node addition."</p>
|
|
|
|
<p class="whs3"> </p>
|
|
|
|
<p class="whs4"># add tetrahedron</p>
|
|
|
|
<p class="whs3">t1 = aMeshEditor.AddVolume([n1,
|
|
38, 39, 246])</p>
|
|
|
|
<p class="whs3">if t1 == 0: print
|
|
"KO tetrahedron addition."</p>
|
|
|
|
<p class="whs2"><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="whs3"><span style="font-family: 'Lucida Console', monospace;">import
|
|
SMESH</span></p>
|
|
|
|
<p class="whs3">import SMESH_mechanic</p>
|
|
|
|
<p class="whs3"> </p>
|
|
|
|
<p class="whs3">mesh = SMESH_mechanic.mesh</p>
|
|
|
|
<p class="whs3">aMeshEditor = mesh.GetMeshEditor()</p>
|
|
|
|
<p class="whs3">print ""</p>
|
|
|
|
<p class="whs3"> </p>
|
|
|
|
<p class="whs4"># add nodes</p>
|
|
|
|
<p class="whs3">nId1 = aMeshEditor.AddNode(50,
|
|
10, 0)</p>
|
|
|
|
<p class="whs3">nId2 = aMeshEditor.AddNode(47,
|
|
12, 0)</p>
|
|
|
|
<p class="whs3">nId3 = aMeshEditor.AddNode(50,
|
|
10, 10)</p>
|
|
|
|
<p class="whs3">nId4 = aMeshEditor.AddNode(47,
|
|
12, 10)</p>
|
|
|
|
<p class="whs3"> </p>
|
|
|
|
<p class="whs3">if nId1 == 0 or
|
|
nId2 == 0 or nId3 == 0 or nId4 == 0: print "KO node addition."</p>
|
|
|
|
<p class="whs3"> </p>
|
|
|
|
<p class="whs4"># add hexahedron</p>
|
|
|
|
<p class="whs3">vId = aMeshEditor.AddVolume([nId2,
|
|
nId1, 38, 39, nId4, nId3, 245, 246])</p>
|
|
|
|
<p class="whs3">if vId == 0: print
|
|
"KO Hexahedron addition."</p>
|
|
|
|
<p class="whs3">else: print
|
|
"New Hexahedron has been added with ID ", vId </p>
|
|
|
|
<p class="whs3"> </p>
|
|
|
|
<h4>Add Polygon</h4>
|
|
|
|
<p class="whs3"><span style="font-family: 'Lucida Console', monospace;">import
|
|
math</span></p>
|
|
|
|
<p class="whs3">import salome</p>
|
|
|
|
<p class="whs3">import smesh</p>
|
|
|
|
<p class="whs3"> </p>
|
|
|
|
<p class="whs4"># create an empty mesh
|
|
structure</p>
|
|
|
|
<p class="whs3">gen = smesh.smesh</p>
|
|
|
|
<p class="whs3">mesh = gen.CreateEmptyMesh()</p>
|
|
|
|
<p class="whs3">MeshEditor = mesh.GetMeshEditor()</p>
|
|
|
|
<p class="whs3"> </p>
|
|
|
|
<p class="whs4"># a method to build
|
|
a polygonal mesh element with <nb_vert> angles:</p>
|
|
|
|
<p class="whs3">def MakePolygon
|
|
(a_mesh, x0, y0, z0, radius, nb_vert):</p>
|
|
|
|
<p class="whs3"> al
|
|
= 2.0 * math.pi / nb_vert</p>
|
|
|
|
<p class="whs3"> node_ids
|
|
= []</p>
|
|
|
|
<p class="whs3"> </p>
|
|
|
|
<p class="whs4"> #
|
|
Create nodes for a polygon</p>
|
|
|
|
<p class="whs3"> for
|
|
ii in range(nb_vert):</p>
|
|
|
|
<p class="whs3"> nid
|
|
= MeshEditor.AddNode(x0 + radius * math.cos(ii*al),</p>
|
|
|
|
<p class="whs3"> y0
|
|
+ radius * math.sin(ii*al),</p>
|
|
|
|
<p class="whs3"> z0)</p>
|
|
|
|
<p class="whs3"> node_ids.append(nid)</p>
|
|
|
|
<p class="whs3"> pass</p>
|
|
|
|
<p class="whs3"> </p>
|
|
|
|
<p class="whs4"> #
|
|
Create a polygon</p>
|
|
|
|
<p class="whs3"> return
|
|
MeshEditor.AddPolygonalFace(node_ids)</p>
|
|
|
|
<p class="whs3"> </p>
|
|
|
|
<p class="whs4"># Create three polygons</p>
|
|
|
|
<p class="whs3">f1 = MakePolygon(mesh,
|
|
0, 0, 0,
|
|
30, 13)</p>
|
|
|
|
<p class="whs3">f2 = MakePolygon(mesh,
|
|
0, 0, 10, 21, 9)</p>
|
|
|
|
<p class="whs3">f3 = MakePolygon(mesh,
|
|
0, 0, 20, 13, 6)</p>
|
|
|
|
<p class="whs3"> </p>
|
|
|
|
<p class="whs3">salome.sg.updateObjBrowser(1)
|
|
</p>
|
|
|
|
<h4>Add Polyhedron</h4>
|
|
|
|
<p class="whs3"><span style="font-family: 'Lucida Console', monospace;">import
|
|
salome</span></p>
|
|
|
|
<p class="whs3">import math</p>
|
|
|
|
<p class="whs3">import smesh</p>
|
|
|
|
<p class="whs3"> </p>
|
|
|
|
<p class="whs4"># create an empty mesh
|
|
structure</p>
|
|
|
|
<p class="whs3">gen = smesh.smesh</p>
|
|
|
|
<p class="whs3">mesh = gen.CreateEmptyMesh()</p>
|
|
|
|
<p class="whs3">MeshEditor = mesh.GetMeshEditor()</p>
|
|
|
|
<p class="whs3"> </p>
|
|
|
|
<p class="whs4"># Create nodes for
|
|
12-hedron with pentagonal faces</p>
|
|
|
|
<p class="whs3">al = 2 * math.pi
|
|
/ 5.0</p>
|
|
|
|
<p class="whs3">cosal = math.cos(al)</p>
|
|
|
|
<p class="whs3">aa = 13</p>
|
|
|
|
<p class="whs3">rr = aa / (2.0 *
|
|
math.sin(al/2.0))</p>
|
|
|
|
<p class="whs3">dr = 2.0 * rr *
|
|
cosal</p>
|
|
|
|
<p class="whs3">r1 = rr + dr</p>
|
|
|
|
<p class="whs3">dh = rr * math.sqrt(2.0
|
|
* (1.0 - cosal * (1.0 + 2.0 * cosal)))</p>
|
|
|
|
<p class="whs3">hh = 2.0 * dh -
|
|
dr * (rr*(cosal - 1) + (rr + dr)*(math.cos(al/2) - 1)) / dh</p>
|
|
|
|
<p class="whs3"> </p>
|
|
|
|
<p class="whs3">dd = [] <span style="font-family: 'Times New Roman', serif;">#
|
|
top</span></p>
|
|
|
|
<p class="whs3">cc = [] <span style="font-family: 'Times New Roman', serif;">#
|
|
below top</span></p>
|
|
|
|
<p class="whs3">bb = [] <span style="font-family: 'Times New Roman', serif;">#
|
|
above bottom</span></p>
|
|
|
|
<p class="whs3">aa = [] <span style="font-family: 'Times New Roman', serif;">#
|
|
bottom</span></p>
|
|
|
|
<p class="whs3"> </p>
|
|
|
|
<p class="whs3">for i in range(5):</p>
|
|
|
|
<p class="whs3"> cos_bot
|
|
= math.cos(i*al)</p>
|
|
|
|
<p class="whs3"> sin_bot
|
|
= math.sin(i*al)</p>
|
|
|
|
<p class="whs3"> </p>
|
|
|
|
<p class="whs3"> <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="whs3"> sin_top
|
|
= math.sin(i*al + al/2.0)</p>
|
|
|
|
<p class="whs3"> </p>
|
|
|
|
<p class="whs3"> nd
|
|
= MeshEditor.AddNode(rr * cos_top, rr * sin_top, hh )
|
|
<span style="font-family: 'Times New Roman', serif;"># top</span></p>
|
|
|
|
<p class="whs3"> nc
|
|
= MeshEditor.AddNode(r1 * cos_top, r1 * sin_top, hh - dh) <span style="font-family: 'Times New Roman', serif;">#
|
|
below top</span></p>
|
|
|
|
<p class="whs3"> nb
|
|
= MeshEditor.AddNode(r1 * cos_bot, r1 * sin_bot, dh)
|
|
<span style="font-family: 'Times New Roman', serif;"># above bottom</span></p>
|
|
|
|
<p class="whs3"> na
|
|
= MeshEditor.AddNode(rr * cos_bot, rr * sin_bot, 0)
|
|
<span style="font-family: 'Times New Roman', serif;"># bottom</span></p>
|
|
|
|
<p class="whs3"> dd.append(nd)
|
|
<span style="font-family: 'Times New Roman', serif;"># top</span></p>
|
|
|
|
<p class="whs3"> cc.append(nc)
|
|
<span style="font-family: 'Times New Roman', serif;"># below top</span></p>
|
|
|
|
<p class="whs3"> bb.append(nb)
|
|
<span style="font-family: 'Times New Roman', serif;"># above bottom</span></p>
|
|
|
|
<p class="whs3"> aa.append(na)
|
|
<span style="font-family: 'Times New Roman', serif;"># bottom</span></p>
|
|
|
|
<p class="whs3"> </p>
|
|
|
|
<p class="whs3"> pass</p>
|
|
|
|
<p class="whs3"> </p>
|
|
|
|
<p class="whs4"># Create a polyhedral
|
|
volume (12-hedron with pentagonal faces)</p>
|
|
|
|
<p class="whs3">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="whs3"> dd[0],
|
|
cc[0], bb[1], cc[1], dd[1], #
|
|
-</p>
|
|
|
|
<p class="whs3"> dd[1],
|
|
cc[1], bb[2], cc[2], dd[2], #
|
|
-</p>
|
|
|
|
<p class="whs3"> dd[2],
|
|
cc[2], bb[3], cc[3], dd[3], <span style="font-family: 'Times New Roman', serif;">#
|
|
- below top</span></p>
|
|
|
|
<p class="whs3"> dd[3],
|
|
cc[3], bb[4], cc[4], dd[4], #
|
|
-</p>
|
|
|
|
<p class="whs3"> dd[4],
|
|
cc[4], bb[0], cc[0], dd[0], #
|
|
-</p>
|
|
|
|
<p class="whs3"> aa[4],
|
|
bb[4], cc[4], bb[0], aa[0], #
|
|
.</p>
|
|
|
|
<p class="whs3"> aa[3],
|
|
bb[3], cc[3], bb[4], aa[4], #
|
|
.</p>
|
|
|
|
<p class="whs3"> aa[2],
|
|
bb[2], cc[2], bb[3], aa[3], <span style="font-family: 'Times New Roman', serif;">#
|
|
. above bottom</span></p>
|
|
|
|
<p class="whs3"> aa[1],
|
|
bb[1], cc[1], bb[2], aa[2], #
|
|
.</p>
|
|
|
|
<p class="whs3"> aa[0],
|
|
bb[0], cc[0], bb[1], aa[1], #
|
|
.</p>
|
|
|
|
<p class="whs3"> aa[0],
|
|
aa[1], aa[2], aa[3], aa[4]], <span style="font-family: 'Times New Roman', serif;">#
|
|
bottom</span></p>
|
|
|
|
<p class="whs3"> [5,5,5,5,5,5,5,5,5,5,5,5])</p>
|
|
|
|
<p class="whs3">salome.sg.updateObjBrowser(1)
|
|
</p>
|
|
|
|
<p class="whs3"> </p>
|
|
|
|
<h3><a name=bookmark1>Removing Nodes and Elements</a></h3>
|
|
|
|
<p class="whs1"><span style="font-weight: bold; color: #ff0000;"><font color=red ><B>#
|
|
Attention! This script has been written using the old approach basing
|
|
on direct usage of SMESH idl interface.</B></font></span></p>
|
|
|
|
<p class="whs1"># For the moment smesh package doesn't provide
|
|
methods to remove nodes and elements. </p>
|
|
|
|
<p class="whs1"># In the next SALOME version the scripts will
|
|
be updated to use only the commands from smesh package.</p>
|
|
|
|
<h4 class="whs5"> </h4>
|
|
|
|
<h4 class="whs5">Removing Nodes</h4>
|
|
|
|
<p class="whs2"> </p>
|
|
|
|
<p class="whs3"><span style="font-family: 'Lucida Console', monospace;">import
|
|
SMESH_mechanic</span></p>
|
|
|
|
<p class="whs3"> </p>
|
|
|
|
<p class="whs3">mesh = SMESH_mechanic.mesh</p>
|
|
|
|
<p class="whs3">aMeshEditor = mesh.GetMeshEditor()</p>
|
|
|
|
<p class="whs3"> </p>
|
|
|
|
<p class="whs4"># remove nodes #246
|
|
and #255</p>
|
|
|
|
<p class="whs3">res = aMeshEditor.RemoveNodes([246,
|
|
255])</p>
|
|
|
|
<p class="whs3">if res == 1: print
|
|
"Nodes removing is OK!"</p>
|
|
|
|
<p class="whs3">else: print
|
|
"KO nodes removing." </p>
|
|
|
|
<p class="whs2"> </p>
|
|
|
|
<h4>Removing Elements</h4>
|
|
|
|
<p class="whs3"><span style="font-family: 'Lucida Console', monospace;">import
|
|
SMESH_mechanic</span></p>
|
|
|
|
<p class="whs3">mesh = SMESH_mechanic.mesh</p>
|
|
|
|
<p class="whs3"> </p>
|
|
|
|
<p class="whs3">anEditor = mesh.GetMeshEditor()</p>
|
|
|
|
<p class="whs3"> </p>
|
|
|
|
<p class="whs4"># remove three elements:
|
|
#850, #859 and #814</p>
|
|
|
|
<p class="whs3">res = anEditor.RemoveElements([850,
|
|
859, 814])</p>
|
|
|
|
<p class="whs3">if res == 1: print
|
|
"Elements removing is OK!"</p>
|
|
|
|
<p class="whs3">else: print
|
|
"KO Elements removing." </p>
|
|
|
|
<p class="whs3"> </p>
|
|
|
|
<h3><a name=bookmark2>Renumbering Nodes and Elements</a></h3>
|
|
|
|
<p class="whs1"><span style="font-weight: bold; color: #ff0000;"><font color=red ><B>#
|
|
Attention! This script has been written using the old approach basing
|
|
on direct usage of SMESH idl interface.</B></font></span></p>
|
|
|
|
<p class="whs1"># For the moment smesh package doesn't provide
|
|
methods to renumber nodes and elements. </p>
|
|
|
|
<p class="whs1"># In the next SALOME version the scripts will
|
|
be updated to use only the commands from smesh package.</p>
|
|
|
|
<p class="whs6"> </p>
|
|
|
|
<p class="whs6">import SMESH_mechanic</p>
|
|
|
|
<p class="whs6"> </p>
|
|
|
|
<p class="whs6">mesh = SMESH_mechanic.mesh</p>
|
|
|
|
<p class="whs6">anEditor = mesh.GetMeshEditor()</p>
|
|
|
|
<p class="whs7"> </p>
|
|
|
|
<p class="whs7"><span style="font-family: 'Lucida Console', monospace;">anEditor.RenumberNodes()</span>
|
|
</p>
|
|
|
|
<h3><a name=bookmark3>Moving Nodes</a></h3>
|
|
|
|
<p class="whs1"><span style="font-weight: bold; color: #ff0000;"><font color=red ><B>#
|
|
Attention! This script has been written using the old approach basing
|
|
on direct usage of SMESH idl interface.</B></font></span></p>
|
|
|
|
<p class="whs1"># For the moment smesh package doesn't provide
|
|
methods to move nodes. </p>
|
|
|
|
<p class="whs1"># In the next SALOME version the scripts will
|
|
be updated to use only the commands from smesh package.</p>
|
|
|
|
<p class="whs3"> </p>
|
|
|
|
<p class="whs3">import SMESH</p>
|
|
|
|
<p class="whs3">import SMESH_mechanic</p>
|
|
|
|
<p class="whs3"> </p>
|
|
|
|
<p class="whs3">mesh = SMESH_mechanic.mesh</p>
|
|
|
|
<p class="whs3">aMeshEditor = mesh.GetMeshEditor()</p>
|
|
|
|
<p class="whs2"> </p>
|
|
|
|
<p class="whs2"># move node #38</p>
|
|
|
|
<p class="whs3">aMeshEditor.MoveNode(38,
|
|
20., 10., 0.) </p>
|
|
|
|
<p class="whs2"> </p>
|
|
|
|
<h3><a name=bookmark4>Diagonal Inversion</a></h3>
|
|
|
|
<p class="whs1"><span style="font-weight: bold; color: #ff0000;"><font color=red ><B>#
|
|
Attention! This script has been written using the old approach basing
|
|
on direct usage of SMESH idl interface.</B></font></span></p>
|
|
|
|
<p class="whs1"># For the moment smesh package doesn't provide
|
|
methods to produce a diagonal inversion. </p>
|
|
|
|
<p class="whs1"># In the next SALOME version the scripts will
|
|
be updated to use only the commands from smesh package.</p>
|
|
|
|
<p class="whs7"> </p>
|
|
|
|
<p class="whs6">import salome</p>
|
|
|
|
<p class="whs6">import smesh</p>
|
|
|
|
<p class="whs6"> </p>
|
|
|
|
<p class="whs8"># create an empty mesh
|
|
structure</p>
|
|
|
|
<p class="whs6">gen = smesh.smesh</p>
|
|
|
|
<p class="whs6">mesh = gen.CreateEmptyMesh()</p>
|
|
|
|
<p class="whs6">aMeshEditor = mesh.GetMeshEditor()</p>
|
|
|
|
<p class="whs6"> </p>
|
|
|
|
<p class="whs8"># create the following
|
|
mesh:</p>
|
|
|
|
<p class="whs6"># .----.----.----.</p>
|
|
|
|
<p class="whs6"># | /|
|
|
/|
|
|
/|</p>
|
|
|
|
<p class="whs6"># | /
|
|
| / | /
|
|
|</p>
|
|
|
|
<p class="whs6"># | / |
|
|
/ | / |</p>
|
|
|
|
<p class="whs6"># |/ |/
|
|
|/
|
|
|</p>
|
|
|
|
<p class="whs6"># .----.----.----.</p>
|
|
|
|
<p class="whs6"> </p>
|
|
|
|
<p class="whs6">bb = [0, 0, 0, 0]</p>
|
|
|
|
<p class="whs6">tt = [0, 0, 0, 0]</p>
|
|
|
|
<p class="whs6">ff = [0, 0, 0, 0,
|
|
0, 0]</p>
|
|
|
|
<p class="whs6"> </p>
|
|
|
|
<p class="whs6">bb[0] = aMeshEditor.AddNode(
|
|
0., 0., 0.)</p>
|
|
|
|
<p class="whs6">bb[1] = aMeshEditor.AddNode(10.,
|
|
0., 0.)</p>
|
|
|
|
<p class="whs6">bb[2] = aMeshEditor.AddNode(20.,
|
|
0., 0.)</p>
|
|
|
|
<p class="whs6">bb[3] = aMeshEditor.AddNode(30.,
|
|
0., 0.)</p>
|
|
|
|
<p class="whs6"> </p>
|
|
|
|
<p class="whs6">tt[0] = aMeshEditor.AddNode(
|
|
0., 15., 0.)</p>
|
|
|
|
<p class="whs6">tt[1] = aMeshEditor.AddNode(10.,
|
|
15., 0.)</p>
|
|
|
|
<p class="whs6">tt[2] = aMeshEditor.AddNode(20.,
|
|
15., 0.)</p>
|
|
|
|
<p class="whs6">tt[3] = aMeshEditor.AddNode(30.,
|
|
15., 0.)</p>
|
|
|
|
<p class="whs6"> </p>
|
|
|
|
<p class="whs6">ff[0] = aMeshEditor.AddFace([bb[0],
|
|
bb[1], tt[1]])</p>
|
|
|
|
<p class="whs6">ff[1] = aMeshEditor.AddFace([bb[0],
|
|
tt[1], tt[0]])</p>
|
|
|
|
<p class="whs6">ff[2] = aMeshEditor.AddFace([bb[1],
|
|
bb[2], tt[2]])</p>
|
|
|
|
<p class="whs6">ff[3] = aMeshEditor.AddFace([bb[1],
|
|
tt[2], tt[1]])</p>
|
|
|
|
<p class="whs6">ff[4] = aMeshEditor.AddFace([bb[2],
|
|
bb[3], tt[3]])</p>
|
|
|
|
<p class="whs6">ff[5] = aMeshEditor.AddFace([bb[2],
|
|
tt[3], tt[2]])</p>
|
|
|
|
<p class="whs6"> </p>
|
|
|
|
<p class="whs8"># inverse the diagonal
|
|
bb[1] - tt[2]</p>
|
|
|
|
<p class="whs6">print "\nDiagonal
|
|
inversion ... ",</p>
|
|
|
|
<p class="whs6">res = aMeshEditor.InverseDiag(bb[1],
|
|
tt[2])</p>
|
|
|
|
<p class="whs6">if not res: print
|
|
"failed!"</p>
|
|
|
|
<p class="whs6">else: print
|
|
"done."</p>
|
|
|
|
<p class="whs7"> </p>
|
|
|
|
<p class="whs7"><span style="font-family: 'Lucida Console', monospace;">salome.sg.updateObjBrowser(1)</span>
|
|
</p>
|
|
|
|
<h3><a name=bookmark5>Uniting two Triangles</a></h3>
|
|
|
|
<p class="whs1"><span style="font-weight: bold; color: #ff0000;"><font color=red ><B>#
|
|
Attention! This script has been written using the old approach basing
|
|
on direct usage of SMESH idl interface.</B></font></span></p>
|
|
|
|
<p class="whs1"># For the moment smesh package doesn't provide
|
|
methods to unite two triangles. </p>
|
|
|
|
<p class="whs1"># In the next SALOME version the scripts will
|
|
be updated to use only the commands from smesh package.</p>
|
|
|
|
<p> </p>
|
|
|
|
<p class="whs9">import salome</p>
|
|
|
|
<p class="whs9">import smesh</p>
|
|
|
|
<p> </p>
|
|
|
|
<p># create an empty mesh structure</p>
|
|
|
|
<p class="whs9">gen = smesh.smesh</p>
|
|
|
|
<p class="whs9">mesh = gen.CreateEmptyMesh()</p>
|
|
|
|
<p class="whs9">aMeshEditor = mesh.GetMeshEditor()</p>
|
|
|
|
<p> </p>
|
|
|
|
<p># create the following mesh:</p>
|
|
|
|
<p class="whs9"># .----.----.----.</p>
|
|
|
|
<p class="whs9"># | /|
|
|
/|
|
|
/|</p>
|
|
|
|
<p class="whs9"># | /
|
|
| / | /
|
|
|</p>
|
|
|
|
<p class="whs9"># | / |
|
|
/ | / |</p>
|
|
|
|
<p class="whs9"># |/ |/
|
|
|/
|
|
|</p>
|
|
|
|
<p class="whs9"># .----.----.----.</p>
|
|
|
|
<p> </p>
|
|
|
|
<p class="whs9">bb = [0, 0, 0, 0]</p>
|
|
|
|
<p class="whs9">tt = [0, 0, 0, 0]</p>
|
|
|
|
<p class="whs9">ff = [0, 0, 0, 0,
|
|
0, 0]</p>
|
|
|
|
<p class="whs9"> </p>
|
|
|
|
<p class="whs9">bb[0] = aMeshEditor.AddNode(
|
|
0., 0., 0.)</p>
|
|
|
|
<p class="whs9">bb[1] = aMeshEditor.AddNode(10.,
|
|
0., 0.)</p>
|
|
|
|
<p class="whs9">bb[2] = aMeshEditor.AddNode(20.,
|
|
0., 0.)</p>
|
|
|
|
<p class="whs9">bb[3] = aMeshEditor.AddNode(30.,
|
|
0., 0.)</p>
|
|
|
|
<p class="whs9"> </p>
|
|
|
|
<p class="whs9">tt[0] = aMeshEditor.AddNode(
|
|
0., 15., 0.)</p>
|
|
|
|
<p class="whs9">tt[1] = aMeshEditor.AddNode(10.,
|
|
15., 0.)</p>
|
|
|
|
<p class="whs9">tt[2] = aMeshEditor.AddNode(20.,
|
|
15., 0.)</p>
|
|
|
|
<p class="whs9">tt[3] = aMeshEditor.AddNode(30.,
|
|
15., 0.)</p>
|
|
|
|
<p class="whs9"> </p>
|
|
|
|
<p class="whs9">ff[0] = aMeshEditor.AddFace([bb[0],
|
|
bb[1], tt[1]])</p>
|
|
|
|
<p class="whs9">ff[1] = aMeshEditor.AddFace([bb[0],
|
|
tt[1], tt[0]])</p>
|
|
|
|
<p class="whs9">ff[2] = aMeshEditor.AddFace([bb[1],
|
|
bb[2], tt[2]])</p>
|
|
|
|
<p class="whs9">ff[3] = aMeshEditor.AddFace([bb[1],
|
|
tt[2], tt[1]])</p>
|
|
|
|
<p class="whs9">ff[4] = aMeshEditor.AddFace([bb[2],
|
|
bb[3], tt[3]])</p>
|
|
|
|
<p class="whs9">ff[5] = aMeshEditor.AddFace([bb[2],
|
|
tt[3], tt[2]])</p>
|
|
|
|
<p> </p>
|
|
|
|
<p># delete the diagonal bb[1] - tt[2]</p>
|
|
|
|
<p class="whs9">print "\nUnite
|
|
two triangles ... ",</p>
|
|
|
|
<p class="whs9">res = aMeshEditor.DeleteDiag(bb[1],
|
|
tt[2])</p>
|
|
|
|
<p class="whs9">if not res: print
|
|
"failed!"</p>
|
|
|
|
<p class="whs9">else: print
|
|
"done."</p>
|
|
|
|
<p class="whs9"> </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="whs1"><span style="font-weight: bold; color: #ff0000;"><font color=red ><B>#
|
|
Attention! This script has been written using the old approach basing
|
|
on direct usage of SMESH idl interface.</B></font></span></p>
|
|
|
|
<p class="whs1"># For the moment smesh package doesn't provide
|
|
methods to unite a set of triangles. </p>
|
|
|
|
<p class="whs1"># In the next SALOME version the scripts will
|
|
be updated to use only the commands from smesh package.</p>
|
|
|
|
<p class="whs10"> </p>
|
|
|
|
<p class="whs9">import salome</p>
|
|
|
|
<p class="whs9">import smesh</p>
|
|
|
|
<p class="whs9"> </p>
|
|
|
|
<p class="whs10"># create an empty mesh
|
|
structure</p>
|
|
|
|
<p class="whs9">gen = smesh.smesh</p>
|
|
|
|
<p class="whs9">mesh = gen.CreateEmptyMesh()</p>
|
|
|
|
<p class="whs9">aMeshEditor = mesh.GetMeshEditor()</p>
|
|
|
|
<p class="whs9"> </p>
|
|
|
|
<p class="whs10"># create the following
|
|
mesh:</p>
|
|
|
|
<p class="whs9"># .----.----.----.</p>
|
|
|
|
<p class="whs9"># | /|
|
|
/|
|
|
/|</p>
|
|
|
|
<p class="whs9"># | /
|
|
| / | /
|
|
|</p>
|
|
|
|
<p class="whs9"># | / |
|
|
/ | / |</p>
|
|
|
|
<p class="whs9"># |/ |/
|
|
|/
|
|
|</p>
|
|
|
|
<p class="whs9"># .----.----.----.</p>
|
|
|
|
<p class="whs9"> </p>
|
|
|
|
<p class="whs9">bb = [0, 0, 0, 0]</p>
|
|
|
|
<p class="whs9">tt = [0, 0, 0, 0]</p>
|
|
|
|
<p class="whs9">ff = [0, 0, 0, 0,
|
|
0, 0]</p>
|
|
|
|
<p class="whs9"> </p>
|
|
|
|
<p class="whs9">bb[0] = aMeshEditor.AddNode(
|
|
0., 0., 0.)</p>
|
|
|
|
<p class="whs9">bb[1] = aMeshEditor.AddNode(10.,
|
|
0., 0.)</p>
|
|
|
|
<p class="whs9">bb[2] = aMeshEditor.AddNode(20.,
|
|
0., 0.)</p>
|
|
|
|
<p class="whs9">bb[3] = aMeshEditor.AddNode(30.,
|
|
0., 0.)</p>
|
|
|
|
<p class="whs9"> </p>
|
|
|
|
<p class="whs9">tt[0] = aMeshEditor.AddNode(
|
|
0., 15., 0.)</p>
|
|
|
|
<p class="whs9">tt[1] = aMeshEditor.AddNode(10.,
|
|
15., 0.)</p>
|
|
|
|
<p class="whs9">tt[2] = aMeshEditor.AddNode(20.,
|
|
15., 0.)</p>
|
|
|
|
<p class="whs9">tt[3] = aMeshEditor.AddNode(30.,
|
|
15., 0.)</p>
|
|
|
|
<p class="whs9"> </p>
|
|
|
|
<p class="whs9">ff[0] = aMeshEditor.AddFace([bb[0],
|
|
bb[1], tt[1]])</p>
|
|
|
|
<p class="whs9">ff[1] = aMeshEditor.AddFace([bb[0],
|
|
tt[1], tt[0]])</p>
|
|
|
|
<p class="whs9">ff[2] = aMeshEditor.AddFace([bb[1],
|
|
bb[2], tt[2]])</p>
|
|
|
|
<p class="whs9">ff[3] = aMeshEditor.AddFace([bb[1],
|
|
tt[2], tt[1]])</p>
|
|
|
|
<p class="whs9">ff[4] = aMeshEditor.AddFace([bb[2],
|
|
bb[3], tt[3]])</p>
|
|
|
|
<p class="whs9">ff[5] = aMeshEditor.AddFace([bb[2],
|
|
tt[3], tt[2]])</p>
|
|
|
|
<p class="whs9"> </p>
|
|
|
|
<p class="whs10"># unite a set of triangles</p>
|
|
|
|
<p class="whs9">aFilterMgr = smesh.smesh.CreateFilterManager()</p>
|
|
|
|
<p class="whs9">aFunctor = aFilterMgr.CreateMinimumAngle()</p>
|
|
|
|
<p class="whs9">print "\nUnite
|
|
a set of triangles ... ",</p>
|
|
|
|
<p class="whs9">res = aMeshEditor.TriToQuad([ff[2],
|
|
ff[3], ff[4], ff[5]], aFunctor, 60.)</p>
|
|
|
|
<p class="whs9">if not res: print
|
|
"failed!"</p>
|
|
|
|
<p class="whs9">else: print
|
|
"done."</p>
|
|
|
|
<p class="whs9"> </p>
|
|
|
|
<p class="whs9">salome.sg.updateObjBrowser(1)
|
|
</p>
|
|
|
|
<h3><a name=bookmark12>Orientation</a></h3>
|
|
|
|
<p class="whs1"><span style="font-weight: bold; color: #ff0000;"><font color=red ><B>#
|
|
Attention! This script has been written using the old approach basing
|
|
on direct usage of SMESH idl interface.</B></font></span></p>
|
|
|
|
<p class="whs1"># For the moment smesh package doesn't provide
|
|
methods to change orientation. </p>
|
|
|
|
<p class="whs1"># In the next SALOME version the scripts will
|
|
be updated to use only the commands from smesh package.</p>
|
|
|
|
<p class="whs9"> </p>
|
|
|
|
<p class="whs9">import salome</p>
|
|
|
|
<p class="whs9">import smesh</p>
|
|
|
|
<p class="whs9"> </p>
|
|
|
|
<p class="whs10"># SMESH module</p>
|
|
|
|
<p class="whs9">gen = smesh.smesh</p>
|
|
|
|
<p class="whs9">mesh = gen.CreateEmptyMesh()</p>
|
|
|
|
<p class="whs9">MeshEditor = mesh.GetMeshEditor()</p>
|
|
|
|
<p class="whs9"> </p>
|
|
|
|
<p class="whs10"># build five quadrangles:</p>
|
|
|
|
<p class="whs9">dx = 10</p>
|
|
|
|
<p class="whs9">dy = 20</p>
|
|
|
|
<p class="whs9"> </p>
|
|
|
|
<p class="whs9">n1 =
|
|
MeshEditor.AddNode(0.0 * dx, 0, 0)</p>
|
|
|
|
<p class="whs9">n2 =
|
|
MeshEditor.AddNode(1.0 * dx, 0, 0)</p>
|
|
|
|
<p class="whs9">n3 =
|
|
MeshEditor.AddNode(2.0 * dx, 0, 0)</p>
|
|
|
|
<p class="whs9">n4 =
|
|
MeshEditor.AddNode(3.0 * dx, 0, 0)</p>
|
|
|
|
<p class="whs9">n5 =
|
|
MeshEditor.AddNode(4.0 * dx, 0, 0)</p>
|
|
|
|
<p class="whs9">n6 =
|
|
MeshEditor.AddNode(5.0 * dx, 0, 0)</p>
|
|
|
|
<p class="whs9">n7 =
|
|
MeshEditor.AddNode(0.0 * dx, dy, 0)</p>
|
|
|
|
<p class="whs9">n8 =
|
|
MeshEditor.AddNode(1.0 * dx, dy, 0)</p>
|
|
|
|
<p class="whs9">n9 =
|
|
MeshEditor.AddNode(2.0 * dx, dy, 0)</p>
|
|
|
|
<p class="whs9">n10 = MeshEditor.AddNode(3.0
|
|
* dx, dy, 0)</p>
|
|
|
|
<p class="whs9">n11 = MeshEditor.AddNode(4.0
|
|
* dx, dy, 0)</p>
|
|
|
|
<p class="whs9">n12 = MeshEditor.AddNode(5.0
|
|
* dx, dy, 0)</p>
|
|
|
|
<p class="whs9"> </p>
|
|
|
|
<p class="whs9">f1 = MeshEditor.AddFace([n1,
|
|
n2, n8 , n7 ])</p>
|
|
|
|
<p class="whs9">f2 = MeshEditor.AddFace([n2,
|
|
n3, n9 , n8 ])</p>
|
|
|
|
<p class="whs9">f3 = MeshEditor.AddFace([n3,
|
|
n4, n10, n9 ])</p>
|
|
|
|
<p class="whs9">f4 = MeshEditor.AddFace([n4,
|
|
n5, n11, n10])</p>
|
|
|
|
<p class="whs9">f5 = MeshEditor.AddFace([n5,
|
|
n6, n12, n11])</p>
|
|
|
|
<p class="whs9"> </p>
|
|
|
|
<p class="whs10"># Change the orientation
|
|
of the second and the fourth faces.</p>
|
|
|
|
<p class="whs9">MeshEditor.Reorient([2,
|
|
4])</p>
|
|
|
|
<p class="whs9"> </p>
|
|
|
|
<p class="whs9">salome.sg.updateObjBrowser(1)
|
|
</p>
|
|
|
|
<p class="whs9"> </p>
|
|
|
|
<h3><a name=bookmark7>Cutting Quadrangles</a></h3>
|
|
|
|
<p class="whs1"><span style="font-weight: bold; color: #ff0000;"><font color=red ><B>#
|
|
Attention! This script has been written using the old approach basing
|
|
on direct usage of SMESH idl interface.</B></font></span></p>
|
|
|
|
<p class="whs1"># For the moment smesh package doesn't provide
|
|
methods to cut quadrangles. </p>
|
|
|
|
<p class="whs1"># In the next SALOME version the scripts will
|
|
be updated to use only the commands from smesh package.</p>
|
|
|
|
<p class="whs9"> </p>
|
|
|
|
<p class="whs9">import SMESH</p>
|
|
|
|
<p class="whs9">import SMESH_mechanic</p>
|
|
|
|
<p class="whs9"> </p>
|
|
|
|
<p class="whs9">smesh = SMESH_mechanic.smesh</p>
|
|
|
|
<p class="whs9">mesh =
|
|
SMESH_mechanic.mesh</p>
|
|
|
|
<p class="whs9"> </p>
|
|
|
|
<p class="whs10"># cut two quadrangles:
|
|
405 and 406</p>
|
|
|
|
<p class="whs9">aFilterMgr = smesh.CreateFilterManager()</p>
|
|
|
|
<p class="whs9">aFunctor = aFilterMgr.CreateMinimumAngle()</p>
|
|
|
|
<p class="whs9"> </p>
|
|
|
|
<p class="whs9">aMeshEditor = mesh.GetMeshEditor()</p>
|
|
|
|
<p class="whs9">aMeshEditor.QuadToTri([405,
|
|
406], aFunctor) </p>
|
|
|
|
<p class="whs10"> </p>
|
|
|
|
<h3><a name=bookmark8>Smoothing</a></h3>
|
|
|
|
<p class="whs1"><span style="font-weight: bold; color: #ff0000;"><font color=red ><B>#
|
|
Attention! This script has been written using the old approach basing
|
|
on direct usage of SMESH idl interface.</B></font></span></p>
|
|
|
|
<p class="whs1"># For the moment smesh package doesn't provide
|
|
methods to perform smoothing. </p>
|
|
|
|
<p class="whs1"># In the next SALOME version the scripts will
|
|
be updated to use only the commands from smesh package.</p>
|
|
|
|
<p class="whs1"> </p>
|
|
|
|
<p class="whs9">import salome</p>
|
|
|
|
<p class="whs9">import geompy</p>
|
|
|
|
<p class="whs9"> </p>
|
|
|
|
<p class="whs9">import SMESH</p>
|
|
|
|
<p class="whs9">import SMESH_mechanic</p>
|
|
|
|
<p class="whs9"> </p>
|
|
|
|
<p class="whs9">mesh = SMESH_mechanic.mesh</p>
|
|
|
|
<p class="whs9">aMeshEditor = mesh.GetMeshEditor()</p>
|
|
|
|
<p class="whs9"> </p>
|
|
|
|
<p class="whs10"># select the top face</p>
|
|
|
|
<p class="whs9">faces = geompy.SubShapeAllSorted(SMESH_mechanic.shape_mesh,
|
|
geompy.ShapeType["FACE"])</p>
|
|
|
|
<p class="whs9">face = faces[3]</p>
|
|
|
|
<p class="whs9">geompy.addToStudyInFather(SMESH_mechanic.shape_mesh,
|
|
face, "face planar with hole")</p>
|
|
|
|
<p class="whs9"> </p>
|
|
|
|
<p class="whs10"># create a group of faces
|
|
to be smoothed</p>
|
|
|
|
<p class="whs9">GroupSmooth = mesh.CreateGroupFromGEOM(SMESH.FACE,
|
|
"Group of faces (smooth)", face)</p>
|
|
|
|
<p class="whs9"> </p>
|
|
|
|
<p class="whs10"># perform smoothing</p>
|
|
|
|
<p class="whs10"># boolean SmoothObject(Object,
|
|
IDsOfFixedNodes, MaxNbOfIterations, MaxAspectRatio, Method)</p>
|
|
|
|
<p class="whs9">res = aMeshEditor.SmoothObject(GroupSmooth,
|
|
[], 20, 2., SMESH.SMESH_MeshEditor.CENTROIDAL_SMOOTH)</p>
|
|
|
|
<p class="whs9"> </p>
|
|
|
|
<p class="whs9">print "\nSmoothing
|
|
... ",</p>
|
|
|
|
<p class="whs9">if not res: print
|
|
"failed!"</p>
|
|
|
|
<p class="whs9">else: print
|
|
"done."</p>
|
|
|
|
<p class="whs9"> </p>
|
|
|
|
<p class="whs9">salome.sg.updateObjBrowser(1)
|
|
</p>
|
|
|
|
<p class="whs9"> </p>
|
|
|
|
<h3><a name=bookmark9>Extrusion</a></h3>
|
|
|
|
<p class="whs1"><span style="font-weight: bold; color: #ff0000;"><font color=red ><B>#
|
|
Attention! This script has been written using the old approach basing
|
|
on direct usage of SMESH idl interface.</B></font></span></p>
|
|
|
|
<p class="whs1"># For the moment smesh package doesn't provide
|
|
methods to perform extrusion. </p>
|
|
|
|
<p class="whs1"># In the next SALOME version the scripts will
|
|
be updated to use only the commands from smesh package.</p>
|
|
|
|
<p class="whs9"> </p>
|
|
|
|
<p class="whs9">import salome</p>
|
|
|
|
<p class="whs9">import geompy</p>
|
|
|
|
<p class="whs9"> </p>
|
|
|
|
<p class="whs9">import SMESH</p>
|
|
|
|
<p class="whs9">import SMESH_mechanic</p>
|
|
|
|
<p class="whs9"> </p>
|
|
|
|
<p class="whs9">mesh = SMESH_mechanic.mesh</p>
|
|
|
|
<p class="whs9">aMeshEditor = mesh.GetMeshEditor()</p>
|
|
|
|
<p class="whs9"> </p>
|
|
|
|
<p class="whs10"># select the top face</p>
|
|
|
|
<p class="whs9">faces = geompy.SubShapeAllSorted(SMESH_mechanic.shape_mesh,
|
|
geompy.ShapeType["FACE"])</p>
|
|
|
|
<p class="whs9">face = faces[7]</p>
|
|
|
|
<p class="whs9">geompy.addToStudyInFather(SMESH_mechanic.shape_mesh,
|
|
face, "face circular top")</p>
|
|
|
|
<p class="whs9"> </p>
|
|
|
|
<p class="whs10"># create a vector for
|
|
extrusion</p>
|
|
|
|
<p class="whs9">point = SMESH.PointStruct(0.,
|
|
0., 5.)</p>
|
|
|
|
<p class="whs9">vector = SMESH.DirStruct(point)</p>
|
|
|
|
<p class="whs9"> </p>
|
|
|
|
<p class="whs10"># create a group to be
|
|
extruded</p>
|
|
|
|
<p class="whs9">GroupTri = mesh.CreateGroupFromGEOM(SMESH.FACE,
|
|
"Group of faces (extrusion)", face)</p>
|
|
|
|
<p class="whs9"> </p>
|
|
|
|
<p class="whs10"># perform extrusion of
|
|
the group</p>
|
|
|
|
<p class="whs9">aMeshEditor.ExtrusionSweepObject(GroupTri,
|
|
vector, 5)</p>
|
|
|
|
<p class="whs10"> </p>
|
|
|
|
<p class="whs10"><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="whs1"><span style="font-weight: bold; color: #ff0000;"><font color=red ><B>#
|
|
Attention! This script has been written using the old approach basing
|
|
on direct usage of SMESH idl interface.</B></font></span></p>
|
|
|
|
<p class="whs1"># For the moment smesh package doesn't provide
|
|
methods to perform extrusion along a path. </p>
|
|
|
|
<p class="whs1"># In the next SALOME version the scripts will
|
|
be updated to use only the commands from smesh package.</p>
|
|
|
|
<p class="whs9"> </p>
|
|
|
|
<p class="whs9">import math</p>
|
|
|
|
<p class="whs9">import salome</p>
|
|
|
|
<p class="whs9"> </p>
|
|
|
|
<p class="whs10"># Geometry</p>
|
|
|
|
<p class="whs9">import geompy</p>
|
|
|
|
<p class="whs10"> </p>
|
|
|
|
<p class="whs10"># 1. Create points</p>
|
|
|
|
<p class="whs9">points = [[0, 0],
|
|
[50, 30], [50, 110], [0, 150], [-80, 150], [-130, 70], [-130, -20]]</p>
|
|
|
|
<p class="whs9"> </p>
|
|
|
|
<p class="whs9">iv = 1</p>
|
|
|
|
<p class="whs9">vertices = []</p>
|
|
|
|
<p class="whs9">for point in points:</p>
|
|
|
|
<p class="whs11">vert
|
|
= geompy.MakeVertex(point[0], point[1], 0)</p>
|
|
|
|
<p class="whs11">geompy.addToStudy(vert,
|
|
"Vertex_" + `iv`)</p>
|
|
|
|
<p class="whs11">vertices.append(vert)</p>
|
|
|
|
<p class="whs11">iv
|
|
+= 1</p>
|
|
|
|
<p class="whs11">pass</p>
|
|
|
|
<p class="whs10"> </p>
|
|
|
|
<p class="whs10"># 2. Create edges and
|
|
wires</p>
|
|
|
|
<p class="whs9">Edge_straight = geompy.MakeEdge(vertices[0],
|
|
vertices[4])</p>
|
|
|
|
<p class="whs9">Edge_bezierrr = geompy.MakeBezier(vertices)</p>
|
|
|
|
<p class="whs9">Wire_polyline = geompy.MakePolyline(vertices)</p>
|
|
|
|
<p class="whs9">Edge_Circle =
|
|
geompy.MakeCircleThreePnt(vertices[0], vertices[1], vertices[2])</p>
|
|
|
|
<p class="whs9"> </p>
|
|
|
|
<p class="whs9">geompy.addToStudy(Edge_straight,
|
|
"Edge_straight")</p>
|
|
|
|
<p class="whs9">geompy.addToStudy(Edge_bezierrr,
|
|
"Edge_bezierrr")</p>
|
|
|
|
<p class="whs9">geompy.addToStudy(Wire_polyline,
|
|
"Wire_polyline")</p>
|
|
|
|
<p class="whs9">geompy.addToStudy(Edge_Circle
|
|
, "Edge_Circle")</p>
|
|
|
|
<p class="whs10"> </p>
|
|
|
|
<p class="whs10"># 3. Explode wire on
|
|
edges, as they will be used for mesh extrusion</p>
|
|
|
|
<p class="whs9">Wire_polyline_edges
|
|
= geompy.SubShapeAll(Wire_polyline, geompy.ShapeType["EDGE"])</p>
|
|
|
|
<p class="whs9">for ii in range(len(Wire_polyline_edges)):</p>
|
|
|
|
<p class="whs11">geompy.addToStudyInFather(Wire_polyline,
|
|
Wire_polyline_edges[ii], "Edge_" + `ii + 1`)</p>
|
|
|
|
<p class="whs11">pass</p>
|
|
|
|
<p class="whs12"> </p>
|
|
|
|
<p class="whs10"># Mesh</p>
|
|
|
|
<p class="whs9">import smesh</p>
|
|
|
|
<p class="whs9">import SMESH</p>
|
|
|
|
<p class="whs9"> </p>
|
|
|
|
<p class="whs9">gen = smesh.smesh</p>
|
|
|
|
<p class="whs9"> </p>
|
|
|
|
<p class="whs9">smeshgui = salome.ImportComponentGUI("SMESH")</p>
|
|
|
|
<p class="whs9">smeshgui.Init(salome.myStudyId)</p>
|
|
|
|
<p class="whs10"> </p>
|
|
|
|
<p class="whs10"># 1D algorithm and three
|
|
1D hypotheses</p>
|
|
|
|
<p class="whs9">Wire_discretisation
|
|
= gen.CreateHypothesis('Regular_1D', 'libStdMeshersEngine.so')</p>
|
|
|
|
<p class="whs9"> </p>
|
|
|
|
<p class="whs9">Nb_Segments_3 = gen.CreateHypothesis('NumberOfSegments',
|
|
'libStdMeshersEngine.so')</p>
|
|
|
|
<p class="whs9">Nb_Segments_7 = gen.CreateHypothesis('NumberOfSegments',
|
|
'libStdMeshersEngine.so')</p>
|
|
|
|
<p class="whs9">Nb_Segments_8 = gen.CreateHypothesis('NumberOfSegments',
|
|
'libStdMeshersEngine.so')</p>
|
|
|
|
<p class="whs9"> </p>
|
|
|
|
<p class="whs9">Nb_Segments_3.SetNumberOfSegments(3)</p>
|
|
|
|
<p class="whs9">Nb_Segments_7.SetNumberOfSegments(7)</p>
|
|
|
|
<p class="whs9">Nb_Segments_8.SetNumberOfSegments(8)</p>
|
|
|
|
<p class="whs10"> </p>
|
|
|
|
<p class="whs10"># Mesh given shape with
|
|
given 1d hypothesis</p>
|
|
|
|
<p class="whs9">def Mesh1D(shape1d,
|
|
hyp1d, name):</p>
|
|
|
|
<p class="whs11">mesh1d_tool
|
|
= smesh.Mesh(shape1d)</p>
|
|
|
|
<p class="whs11">mesh1d
|
|
= mesh1d_tool.GetMesh()</p>
|
|
|
|
<p class="whs11">status
|
|
= mesh1d.AddHypothesis(shape1d, hyp1d)</p>
|
|
|
|
<p class="whs11">status
|
|
= mesh1d.AddHypothesis(shape1d, Wire_discretisation)</p>
|
|
|
|
<p class="whs11">isDone
|
|
= mesh1d_tool.Compute()</p>
|
|
|
|
<p class="whs11">if
|
|
not isDone: print 'Mesh ', name, ': computation failed'</p>
|
|
|
|
<p class="whs11">return
|
|
mesh1d</p>
|
|
|
|
<p class="whs9"> </p>
|
|
|
|
<p class="whs10"># Create a mesh with
|
|
six nodes, seven edges and two quadrangle faces</p>
|
|
|
|
<p class="whs9">def MakeQuadMesh2(mesh_name):</p>
|
|
|
|
<p class="whs11">quad_1
|
|
= gen.CreateEmptyMesh()</p>
|
|
|
|
<p class="whs11">smeshgui.SetName(salome.ObjectToID(quad_1),
|
|
mesh_name)</p>
|
|
|
|
<p class="whs11">editor_1
|
|
= quad_1.GetMeshEditor()</p>
|
|
|
|
<p class="whs12"> </p>
|
|
|
|
<p class="whs12"># six nodes</p>
|
|
|
|
<p class="whs11">n1
|
|
= editor_1.AddNode(0, 20, 10)</p>
|
|
|
|
<p class="whs11">n2
|
|
= editor_1.AddNode(0, 40, 10)</p>
|
|
|
|
<p class="whs11">n3
|
|
= editor_1.AddNode(0, 40, 30)</p>
|
|
|
|
<p class="whs11">n4
|
|
= editor_1.AddNode(0, 20, 30)</p>
|
|
|
|
<p class="whs11">n5
|
|
= editor_1.AddNode(0, 0,
|
|
30)</p>
|
|
|
|
<p class="whs11">n6
|
|
= editor_1.AddNode(0, 0,
|
|
10)</p>
|
|
|
|
<p class="whs12"> </p>
|
|
|
|
<p class="whs12"># seven edges</p>
|
|
|
|
<p class="whs11">editor_1.AddEdge([n1,
|
|
n2]) # 1</p>
|
|
|
|
<p class="whs11">editor_1.AddEdge([n2,
|
|
n3]) # 2</p>
|
|
|
|
<p class="whs11">editor_1.AddEdge([n3,
|
|
n4]) # 3</p>
|
|
|
|
<p class="whs11">editor_1.AddEdge([n4,
|
|
n1]) # 4</p>
|
|
|
|
<p class="whs11">editor_1.AddEdge([n4,
|
|
n5]) # 5</p>
|
|
|
|
<p class="whs11">editor_1.AddEdge([n5,
|
|
n6]) # 6</p>
|
|
|
|
<p class="whs11">editor_1.AddEdge([n6,
|
|
n1]) # 7</p>
|
|
|
|
<p class="whs12"> </p>
|
|
|
|
<p class="whs12"># two quadrangle faces</p>
|
|
|
|
<p class="whs11">editor_1.AddFace([n1,
|
|
n2, n3, n4]) # 8</p>
|
|
|
|
<p class="whs11">editor_1.AddFace([n1,
|
|
n4, n5, n6]) # 9</p>
|
|
|
|
<p class="whs11">return
|
|
[quad_1, editor_1, [1,2,3,4,5,6,7], [8,9]]</p>
|
|
|
|
<p class="whs10"> </p>
|
|
|
|
<p class="whs10"># Path meshes</p>
|
|
|
|
<p class="whs9">Edge_straight_mesh
|
|
= Mesh1D(Edge_straight, Nb_Segments_7, "Edge_straight")</p>
|
|
|
|
<p class="whs9">Edge_bezierrr_mesh
|
|
= Mesh1D(Edge_bezierrr, Nb_Segments_7, "Edge_bezierrr")</p>
|
|
|
|
<p class="whs9">Wire_polyline_mesh
|
|
= Mesh1D(Wire_polyline, Nb_Segments_3, "Wire_polyline")</p>
|
|
|
|
<p class="whs9">Edge_Circle_mesh =
|
|
Mesh1D(Edge_Circle ,
|
|
Nb_Segments_8, "Edge_Circle")</p>
|
|
|
|
<p class="whs10"> </p>
|
|
|
|
<p class="whs10"># Initial meshes (to
|
|
be extruded)</p>
|
|
|
|
<p class="whs9">[quad_1, editor_1,
|
|
ee_1, ff_1] = MakeQuadMesh2("quad_1")</p>
|
|
|
|
<p class="whs9">[quad_2, editor_2,
|
|
ee_2, ff_2] = MakeQuadMesh2("quad_2")</p>
|
|
|
|
<p class="whs9">[quad_3, editor_3,
|
|
ee_3, ff_3] = MakeQuadMesh2("quad_3")</p>
|
|
|
|
<p class="whs9">[quad_4, editor_4,
|
|
ee_4, ff_4] = MakeQuadMesh2("quad_4")</p>
|
|
|
|
<p class="whs9">[quad_5, editor_5,
|
|
ee_5, ff_5] = MakeQuadMesh2("quad_5")</p>
|
|
|
|
<p class="whs9">[quad_6, editor_6,
|
|
ee_6, ff_6] = MakeQuadMesh2("quad_6")</p>
|
|
|
|
<p class="whs9">[quad_7, editor_7,
|
|
ee_7, ff_7] = MakeQuadMesh2("quad_7")</p>
|
|
|
|
<p class="whs10"> </p>
|
|
|
|
<p class="whs10"># ExtrusionAlongPath</p>
|
|
|
|
<p class="whs10"># IDsOfElements, PathMesh,
|
|
PathShape, NodeStart,</p>
|
|
|
|
<p class="whs10"># HasAngles, Angles,
|
|
HasRefPoint, RefPoint</p>
|
|
|
|
<p class="whs10"> </p>
|
|
|
|
<p class="whs9">refPoint = SMESH.PointStruct(0,
|
|
0, 0)</p>
|
|
|
|
<p class="whs9"> </p>
|
|
|
|
<p class="whs9">a10 = 10.0*math.pi/180.0</p>
|
|
|
|
<p class="whs9">a45 = 45.0*math.pi/180.0</p>
|
|
|
|
<p class="whs10"> </p>
|
|
|
|
<p class="whs10"># 1. Extrusion of two
|
|
mesh edges along a straight path</p>
|
|
|
|
<p class="whs9">error = editor_1.ExtrusionAlongPath([1,2],
|
|
Edge_straight_mesh, Edge_straight, 1,</p>
|
|
|
|
<p class="whs13"> 0,
|
|
[], 0, refPoint)</p>
|
|
|
|
<p class="whs10"> </p>
|
|
|
|
<p class="whs10"># 2. Extrusion of one
|
|
mesh edge along a curved path</p>
|
|
|
|
<p class="whs9">error = editor_2.ExtrusionAlongPath([2],
|
|
Edge_bezierrr_mesh, Edge_bezierrr, 1,</p>
|
|
|
|
<p class="whs13"> 0,
|
|
[], 0, refPoint)</p>
|
|
|
|
<p class="whs10"> </p>
|
|
|
|
<p class="whs10"># 3. Extrusion of one
|
|
mesh edge along a curved path with usage of angles </p>
|
|
|
|
<p class="whs9">error = editor_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="whs10"> </p>
|
|
|
|
<p class="whs10"># 4. Extrusion of one
|
|
mesh edge along the path, which is a part of a meshed wire</p>
|
|
|
|
<p class="whs9">error = editor_4.ExtrusionAlongPath([4],
|
|
Wire_polyline_mesh, Wire_polyline_edges[0], 1,</p>
|
|
|
|
<p class="whs13"> 1,
|
|
[a10, a10, a10], 0, refPoint)</p>
|
|
|
|
<p class="whs10"> </p>
|
|
|
|
<p class="whs10"># 5. Extrusion of two
|
|
mesh faces along the path, which is a part of a meshed wire</p>
|
|
|
|
<p class="whs9">error = editor_5.ExtrusionAlongPath(ff_5
|
|
, Wire_polyline_mesh, Wire_polyline_edges[2], 4,</p>
|
|
|
|
<p class="whs13"> 0,
|
|
[], 0, refPoint)</p>
|
|
|
|
<p class="whs10"> </p>
|
|
|
|
<p class="whs10"># 6. Extrusion of two
|
|
mesh faces along a closed path</p>
|
|
|
|
<p class="whs9">error = editor_6.ExtrusionAlongPath(ff_6
|
|
, Edge_Circle_mesh, Edge_Circle, 1,</p>
|
|
|
|
<p class="whs13"> 0,
|
|
[], 0, refPoint)</p>
|
|
|
|
<p class="whs10"> </p>
|
|
|
|
<p class="whs10"># 7. Extrusion of two
|
|
mesh faces along a closed path with usage of angles </p>
|
|
|
|
<p class="whs9">error = editor_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="whs14">salome.sg.updateObjBrowser(1)</p>
|
|
|
|
<h3><a name=bookmark11>Revolution</a></h3>
|
|
|
|
<p class="whs1"><span style="font-weight: bold; color: #ff0000;"><font color=red ><B>#
|
|
Attention! This script has been written using the old approach basing
|
|
on direct usage of SMESH idl interface.</B></font></span></p>
|
|
|
|
<p class="whs1"># For the moment smesh package doesn't provide
|
|
methods to perform revolution. </p>
|
|
|
|
<p class="whs1"># In the next SALOME version the scripts will
|
|
be updated to use only the commands from smesh package.</p>
|
|
|
|
<p class="whs1"> </p>
|
|
|
|
<p class="whs9">import math</p>
|
|
|
|
<p class="whs9">import SMESH</p>
|
|
|
|
<p class="whs9">import SMESH_mechanic</p>
|
|
|
|
<p class="whs9"> </p>
|
|
|
|
<p class="whs9">mesh = SMESH_mechanic.mesh</p>
|
|
|
|
<p class="whs9">aMeshEditor = mesh.GetMeshEditor()</p>
|
|
|
|
<p class="whs9"> </p>
|
|
|
|
<p class="whs10"># create a group of faces
|
|
to be revolved</p>
|
|
|
|
<p class="whs9">FacesRotate = [492,
|
|
493, 502, 503]</p>
|
|
|
|
<p class="whs9">GroupRotate = mesh.CreateGroup(SMESH.FACE,"Group
|
|
of faces (rotate)")</p>
|
|
|
|
<p class="whs9">GroupRotate.Add(FacesRotate)</p>
|
|
|
|
<p class="whs10"> </p>
|
|
|
|
<p class="whs10"># define revolution angle
|
|
and axis</p>
|
|
|
|
<p class="whs9">angle45 = 45 * math.pi
|
|
/ 180</p>
|
|
|
|
<p class="whs9">axisXYZ = SMESH.AxisStruct(-38.3128,
|
|
-73.3658, -23.321, -13.3402, -13.3265, 6.66632)</p>
|
|
|
|
<p class="whs9"> </p>
|
|
|
|
<p class="whs10"># perform revolution
|
|
of an object</p>
|
|
|
|
<p class="whs9">aMeshEditor.RotationSweepObject(GroupRotate,
|
|
axisXYZ, angle45, 4, 1e-5) </p>
|
|
|
|
<p class="whs10"> </p>
|
|
|
|
<h3><a name=bookmark13>Pattern Mapping</a></h3>
|
|
|
|
<p class="whs1"><span style="font-weight: bold; color: #ff0000;"><font color=red ><B>#
|
|
Attention! This script has been written using the old approach basing
|
|
on direct usage of SMESH idl interface.</B></font></span></p>
|
|
|
|
<p class="whs1"># For the moment smesh package doesn't provide
|
|
methods to perform pattern mapping. </p>
|
|
|
|
<p class="whs1"># In the next SALOME version the scripts will
|
|
be updated to use only the commands from smesh package.</p>
|
|
|
|
<p class="whs9"> </p>
|
|
|
|
<p class="whs9">import geompy</p>
|
|
|
|
<p class="whs9">import smesh</p>
|
|
|
|
<p class="whs9"> </p>
|
|
|
|
<p class="whs10"># 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">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="whs10"># 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="whs10"># build a triangle mesh
|
|
on Face_2</p>
|
|
|
|
<p class="whs9">Mesh_2 = smesh.Mesh(Face_2)</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">isDone = Mesh_2.Compute()</p>
|
|
|
|
<p class="whs9">if not isDone: print
|
|
'Mesh Mesh_2 : computation failed'</p>
|
|
|
|
<p class="whs9"> </p>
|
|
|
|
<p class="whs10"># create a pattern</p>
|
|
|
|
<p class="whs9">pattern = smesh.smesh.GetPattern()</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="whs10"># 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">isDone = pattern.MakeMesh(Mesh_1.GetMesh(),
|
|
0, 0)</p>
|
|
|
|
<p class="whs9">if (isDone != 1):
|
|
print 'MakeMesh :', pattern.GetErrorCode() </p>
|
|
|
|
<p class="whs9"> </p>
|
|
|
|
<p class="whs10"> </p>
|
|
|
|
<script type="text/javascript" language="javascript1.2">
|
|
<!--
|
|
if (window.writeIntopicBar)
|
|
writeIntopicBar(0);
|
|
//-->
|
|
</script>
|
|
</body>
|
|
</html>
|