doc/salome/gui/SMESH/transforming_meshes.htm

   1 <!doctype HTML public "-//W3C//DTD HTML 4.0 Frameset//EN">\r
   2 \r
   3 <html>\r
   4 \r
   5 <head>\r
   6 <title>Transforming Meshes</title>\r
   7 <meta http-equiv="content-type" content="text/html; charset=windows-1252">\r
   8 <meta name="generator" content="RoboHelp by eHelp Corporation www.ehelp.com"><style type="text/css">\r
   9 <!--\r
  10 h3.whs1 { margin-top:0pt; margin-bottom:0pt; }\r
  11 h4.whs2 { margin-top:0pt; margin-bottom:0pt; }\r
  12 p.whs3 { font-family:'Lucida Console' , monospace; margin-top:0pt; margin-bottom:0pt; }\r
  13 p.whs4 { margin-top:0pt; margin-bottom:0pt; }\r
  14 p.whs5 { margin-top:0pt; margin-bottom:0pt; font-family:'Lucida Console' , monospace; }\r
  15 -->\r
  16 </style><script type="text/javascript" language="JavaScript">\r
  17 <!--\r
  18 if ((navigator.appName == "Netscape") && (parseInt(navigator.appVersion) == 4))\r
  19 {\r
  20   var strNSS = "<style type='text/css'>";\r
  21   strNSS += "h3.whs1 {margin-top:1pt;margin-bottom:1pt; }";\r
  22   strNSS += "h4.whs2 {margin-top:1pt;margin-bottom:1pt; }";\r
  23   strNSS += "p.whs3 {margin-top:1pt;margin-bottom:1pt; }";\r
  24   strNSS += "p.whs4 {margin-top:1pt;margin-bottom:1pt; }";\r
  25   strNSS += "p.whs5 {margin-top:1pt;margin-bottom:1pt; }";\r
  26   strNSS +="</style>";\r
  27   document.write(strNSS);\r
  28 }\r
  29 //-->\r
  30 </script>\r
  31 <script type="text/javascript" language="JavaScript" title="WebHelpInlineScript">\r
  32 <!--\r
  33 function reDo() {\r
  34   if (innerWidth != origWidth || innerHeight != origHeight)\r
  35      location.reload();\r
  36 }\r
  37 if ((parseInt(navigator.appVersion) == 4) && (navigator.appName == "Netscape")) {\r
  38         origWidth = innerWidth;\r
  39         origHeight = innerHeight;\r
  40         onresize = reDo;\r
  41 }\r
  42 onerror = null; \r
  43 //-->\r
  44 </script>\r
  45 <style type="text/css">\r
  46 <!--\r
  47 div.WebHelpPopupMenu { position:absolute; left:0px; top:0px; z-index:4; visibility:hidden; }\r
  48 -->\r
  49 </style><script type="text/javascript" language="javascript1.2" src="whmsg.js"></script>\r
  50 <script type="text/javascript" language="javascript" src="whver.js"></script>\r
  51 <script type="text/javascript" language="javascript1.2" src="whproxy.js"></script>\r
  52 <script type="text/javascript" language="javascript1.2" src="whutils.js"></script>\r
  53 <script type="text/javascript" language="javascript1.2" src="whtopic.js"></script>\r
  54 </head>\r
  55 <body><script type="text/javascript" language="javascript1.2">\r
  56 <!--\r
  57 if (window.gbWhTopic)\r
  58 {\r
  59         if (window.addTocInfo)\r
  60         {\r
  61         addTocInfo("MESH module\nTUI Scripts\nTransforming Meshes");\r
  62 addButton("show",BTN_TEXT,"Show","","","","",0,0,"","","");\r
  63 \r
  64         }\r
  65         if (window.writeBtnStyle)\r
  66                 writeBtnStyle();\r
  67 \r
  68         if (window.writeIntopicBar)\r
  69                 writeIntopicBar(1);\r
  70 \r
  71         if (window.setRelStartPage)\r
  72         {\r
  73         setRelStartPage("smesh.htm");\r
  74 \r
  75                 autoSync(1);\r
  76                 sendSyncInfo();\r
  77                 sendAveInfoOut();\r
  78         }\r
  79 }\r
  80 else\r
  81         if (window.gbIE4)\r
  82                 document.location.reload();\r
  83 //-->\r
  84 </script>\r
  85 <h1>Transforming Meshes</h1>\r
  86 \r
  87 <h3 class="whs1">Transforming Meshes</h3>\r
  88 \r
  89 <h4 class="whs2">&nbsp;</h4>\r
  90 \r
  91 <h4 class="whs2"><a name=bookmark>Translation</a></h4>\r
  92 \r
  93 <p class="whs3">&nbsp;</p>\r
  94 \r
  95 <p class="whs3">import SMESH</p>\r
  96 \r
  97 <p class="whs3">import SMESH_mechanic</p>\r
  98 \r
  99 <p class="whs3">&nbsp;</p>\r
 100 \r
 101 <p class="whs3">mesh &nbsp;&nbsp;= \r
 102  SMESH_mechanic.mesh</p>\r
 103 \r
 104 <p class="whs3">salome = SMESH_mechanic.salome</p>\r
 105 \r
 106 <p class="whs4">&nbsp;</p>\r
 107 \r
 108 <p class="whs4"># translate a mesh</p>\r
 109 \r
 110 <p class="whs3">anEditor = mesh.GetMeshEditor()</p>\r
 111 \r
 112 <p class="whs3">point = SMESH.PointStruct(-150, -150, \r
 113  0)</p>\r
 114 \r
 115 <p class="whs3">vector = SMESH.DirStruct(point) </p>\r
 116 \r
 117 <p class="whs3">anEditor.TranslateObject(mesh, vector, \r
 118  1)</p>\r
 119 \r
 120 <p class="whs4">&nbsp;</p>\r
 121 \r
 122 <p class="whs4"><span style="font-family: 'Lucida Console', monospace;">salome.sg.updateObjBrowser(1)</span> \r
 123  </p>\r
 124 \r
 125 <p class="whs4">&nbsp;</p>\r
 126 \r
 127 <h4 class="whs2"><a name=bookmark1>Rotation</a></h4>\r
 128 \r
 129 <p class="whs3">&nbsp;</p>\r
 130 \r
 131 <p class="whs3">import SMESH</p>\r
 132 \r
 133 <p class="whs3">import SMESH_mechanic</p>\r
 134 \r
 135 <p class="whs3">&nbsp;</p>\r
 136 \r
 137 <p class="whs3">mport math</p>\r
 138 \r
 139 <p class="whs3">mesh &nbsp;&nbsp;= \r
 140  SMESH_mechanic.mesh</p>\r
 141 \r
 142 <p class="whs3">salome = SMESH_mechanic.salome</p>\r
 143 \r
 144 <p class="whs4">&nbsp;</p>\r
 145 \r
 146 <p class="whs4"># rotate a mesh </p>\r
 147 \r
 148 <p class="whs3">anEditor = mesh.GetMeshEditor()</p>\r
 149 \r
 150 <p class="whs3">axisXYZ = SMESH.AxisStruct(0, 0, 0, 5, \r
 151  5, 20)</p>\r
 152 \r
 153 <p class="whs3">angle180 = &nbsp;1.5*math.pi</p>\r
 154 \r
 155 <p class="whs3">anEditor.RotateObject(mesh, axisXYZ, angle180, \r
 156  1)</p>\r
 157 \r
 158 <p class="whs3">&nbsp;</p>\r
 159 \r
 160 <p class="whs5">salome.sg.updateObjBrowser(1)</p>\r
 161 \r
 162 <p class="whs4">&nbsp;</p>\r
 163 \r
 164 <p class="whs4">&nbsp;</p>\r
 165 \r
 166 <h4 class="whs2"><a name=bookmark2>Symmetry</a></h4>\r
 167 \r
 168 <p class="whs3">&nbsp;</p>\r
 169 \r
 170 <p class="whs3">import SMESH</p>\r
 171 \r
 172 <p class="whs3">import SMESH_mechanic</p>\r
 173 \r
 174 <p class="whs3">&nbsp;</p>\r
 175 \r
 176 <p class="whs3">import math</p>\r
 177 \r
 178 <p class="whs3">mesh &nbsp;&nbsp;= \r
 179  SMESH_mechanic.mesh</p>\r
 180 \r
 181 <p class="whs3">salome = SMESH_mechanic.salome</p>\r
 182 \r
 183 <p class="whs4">&nbsp;</p>\r
 184 \r
 185 <p class="whs4"># create a symmetrical \r
 186  copy of the mesh mirrored through a point</p>\r
 187 \r
 188 <p class="whs3">anEditor = mesh.GetMeshEditor()</p>\r
 189 \r
 190 <p class="whs3">anEditor.MirrorObject(mesh, SMESH.AxisStruct(0, \r
 191  0, 0, 0, 0, 0), SMESH.SMESH_MeshEditor.POINT, 1) </p>\r
 192 \r
 193 <p class="whs3">&nbsp;</p>\r
 194 \r
 195 <p class="whs4"><span style="font-family: 'Lucida Console', monospace;">salome.sg.updateObjBrowser(1)</span> \r
 196  </p>\r
 197 \r
 198 <p class="whs4">&nbsp;</p>\r
 199 \r
 200 <p class="whs4">&nbsp;</p>\r
 201 \r
 202 <h4 class="whs2"><a name=bookmark3>Merging \r
 203  Nodes</a></h4>\r
 204 \r
 205 <p class="whs3">&nbsp;</p>\r
 206 \r
 207 <p class="whs3">import SMESH</p>\r
 208 \r
 209 <p class="whs3">import SMESH_mechanic</p>\r
 210 \r
 211 <p class="whs3">&nbsp;</p>\r
 212 \r
 213 <p class="whs3">smesh &nbsp;= \r
 214  SMESH_mechanic.smesh</p>\r
 215 \r
 216 <p class="whs3">mesh &nbsp;&nbsp;= \r
 217  SMESH_mechanic.mesh</p>\r
 218 \r
 219 <p class="whs3">salome = SMESH_mechanic.salome</p>\r
 220 \r
 221 <p class="whs4">&nbsp;</p>\r
 222 \r
 223 <p class="whs4"># merge nodes</p>\r
 224 \r
 225 <p class="whs3">aMeshEditor = mesh.GetMeshEditor()</p>\r
 226 \r
 227 <p class="whs3">Tolerance = 25.0</p>\r
 228 \r
 229 <p class="whs3">&nbsp;</p>\r
 230 \r
 231 <p class="whs3">GroupsOfNodes = aMeshEditor.FindCoincidentNodes(Tolerance)</p>\r
 232 \r
 233 <p class="whs3">aMeshEditor.MergeNodes(GroupsOfNodes)</p>\r
 234 \r
 235 <p class="whs3">&nbsp;&nbsp;&nbsp;&nbsp;</p>\r
 236 \r
 237 <p class="whs4"><span style="font-family: 'Lucida Console', monospace;">salome.sg.updateObjBrowser(1)</span> \r
 238  </p>\r
 239 \r
 240 <p class="whs4">&nbsp;</p>\r
 241 \r
 242 <h3>Sewing Meshes</h3>\r
 243 \r
 244 <h4><a name=bookmark4>Sew Meshes Border to Border</a></h4>\r
 245 \r
 246 <p class="whs5">import salome</p>\r
 247 \r
 248 <p class="whs5">import geompy</p>\r
 249 \r
 250 <p class="whs5">import smesh</p>\r
 251 \r
 252 <p class="whs4">&nbsp;</p>\r
 253 \r
 254 <p class="whs4"># create two faces of a \r
 255  box</p>\r
 256 \r
 257 <p class="whs5">box1 = geompy.MakeBox(0., \r
 258  0., -10., 30., 20., 25.)</p>\r
 259 \r
 260 <p class="whs5">subShapeList1 = \r
 261  geompy.SubShapeAll(box1, geompy.ShapeType[&quot;FACE&quot;])</p>\r
 262 \r
 263 <p class="whs5">&nbsp;</p>\r
 264 \r
 265 <p class="whs5">box2 = geompy.MakeBox(0., \r
 266  5., 0., 20., 20., 15.)</p>\r
 267 \r
 268 <p class="whs5">subShapeList2 = \r
 269  geompy.SubShapeAll(box2, geompy.ShapeType[&quot;FACE&quot;])</p>\r
 270 \r
 271 <p class="whs5">EdgesList = geompy.SubShapeAll(subShapeList2[ \r
 272  1 ], geompy.ShapeType[&quot;EDGE&quot;])</p>\r
 273 \r
 274 <p class="whs5">&nbsp;</p>\r
 275 \r
 276 <p class="whs5">aComp = geompy.MakeCompound( \r
 277  [subShapeList1[ 2 ], subShapeList2[ 1 ]] )</p>\r
 278 \r
 279 <p class="whs5">idComp = geompy.addToStudy( \r
 280  aComp, &quot;Two faces&quot; )</p>\r
 281 \r
 282 <p class="whs5">aCompobj &nbsp;= \r
 283  salome.IDToObject( idComp )</p>\r
 284 \r
 285 <p class="whs5">&nbsp;</p>\r
 286 \r
 287 <p class="whs4"># create a mesh on two \r
 288  faces</p>\r
 289 \r
 290 <p class="whs4">&nbsp;</p>\r
 291 \r
 292 <p class="whs5">mesh = smesh.Mesh(aCompobj, \r
 293  &quot;Two faces : quadrangle mesh&quot;)</p>\r
 294 \r
 295 <p class="whs5">algo = mesh.Segment()</p>\r
 296 \r
 297 <p class="whs5">&nbsp;</p>\r
 298 \r
 299 <p class="whs4"># define &quot;NumberOfSegments&quot; \r
 300  hypothesis to cut an edge in a fixed number of segments</p>\r
 301 \r
 302 <p class="whs5">algo.NumberOfSegments(9)</p>\r
 303 \r
 304 <p class="whs4">&nbsp;</p>\r
 305 \r
 306 <p class="whs4"># create a quadrangle 2D \r
 307  algorithm for faces</p>\r
 308 \r
 309 <p class="whs5">mesh.Quadrangle()</p>\r
 310 \r
 311 <p class="whs4">&nbsp;</p>\r
 312 \r
 313 <p class="whs4"># create a local hypothesis</p>\r
 314 \r
 315 <p class="whs5">algo = mesh.Segment(EdgesList[2])</p>\r
 316 \r
 317 <p class="whs4">&nbsp;</p>\r
 318 \r
 319 <p class="whs4"># define &quot;Arithmetic1D&quot; \r
 320  hypothesis to cut an edge in several segments with increasing arithmetic \r
 321  length </p>\r
 322 \r
 323 <p class="whs5">algo.Arithmetic1D(1, \r
 324  4)</p>\r
 325 \r
 326 <p class="whs4">&nbsp;</p>\r
 327 \r
 328 <p class="whs4"># define &quot;Propagation&quot; \r
 329  hypothesis that propagates all other hypothesis on all edges on the opposite \r
 330  side in case of quadrangular faces</p>\r
 331 \r
 332 <p class="whs5">algo.Propagation()</p>\r
 333 \r
 334 <p class="whs5">mesh.Compute()</p>\r
 335 \r
 336 <p class="whs4">&nbsp;</p>\r
 337 \r
 338 <p class="whs4"># sew free borders</p>\r
 339 \r
 340 <p class="whs5">anEditor = mesh.GetMesh().GetMeshEditor()</p>\r
 341 \r
 342 <p class="whs5">anEditor.SewBorderToSide(5, \r
 343  45, 6, 113, 109, 0, 0)</p>\r
 344 \r
 345 <p class="whs5">&nbsp;</p>\r
 346 \r
 347 <p class="whs4"><span style="font-family: 'Lucida Console', monospace;">salome.sg.updateObjBrowser(1)</span> \r
 348  </p>\r
 349 \r
 350 <p class="whs4">&nbsp;</p>\r
 351 \r
 352 <h4><a name=bookmark5>Sew Conform Free Borders</a></h4>\r
 353 \r
 354 <p class="whs5">import salome</p>\r
 355 \r
 356 <p class="whs5">import geompy</p>\r
 357 \r
 358 <p class="whs5">import smesh</p>\r
 359 \r
 360 <p class="whs4">&nbsp;</p>\r
 361 \r
 362 <p class="whs4"># create two faces of the \r
 363  box</p>\r
 364 \r
 365 <p class="whs5">box1 = geompy.MakeBox(0., \r
 366  0., -10., 20., 20., 15.)</p>\r
 367 \r
 368 <p class="whs5">subShapeList1 = \r
 369  geompy.SubShapeAll(box1, geompy.ShapeType[&quot;FACE&quot;])</p>\r
 370 \r
 371 <p class="whs5">&nbsp;</p>\r
 372 \r
 373 <p class="whs5">box2 = geompy.MakeBox(0., \r
 374  5., 0., 20., 20., 15.)</p>\r
 375 \r
 376 <p class="whs5">subShapeList2 = \r
 377  geompy.SubShapeAll(box2, geompy.ShapeType[&quot;FACE&quot;])</p>\r
 378 \r
 379 <p class="whs5">EdgesList = geompy.SubShapeAll(subShapeList2[ \r
 380  1 ], geompy.ShapeType[&quot;EDGE&quot;])</p>\r
 381 \r
 382 <p class="whs5">&nbsp;</p>\r
 383 \r
 384 <p class="whs5">aComp = geompy.MakeCompound( \r
 385  [subShapeList1[ 2 ], subShapeList2[ 1 ]] )</p>\r
 386 \r
 387 <p class="whs5">idComp = geompy.addToStudy( \r
 388  aComp, &quot;Two faces&quot; )</p>\r
 389 \r
 390 <p class="whs5">aCompobj &nbsp;= \r
 391  salome.IDToObject( idComp )</p>\r
 392 \r
 393 <p class="whs4">&nbsp;</p>\r
 394 \r
 395 <p class="whs4"># create a mesh on two \r
 396  faces</p>\r
 397 \r
 398 <p class="whs5">mesh = smesh.Mesh(aCompobj, \r
 399  &quot;Two faces : quadrangle mesh&quot;)</p>\r
 400 \r
 401 <p class="whs5">algo = mesh.Segment()</p>\r
 402 \r
 403 <p class="whs4">&nbsp;</p>\r
 404 \r
 405 <p class="whs4"># define &quot;NumberOfSegments&quot; \r
 406  hypothesis to cut an edge in a fixed number of segments</p>\r
 407 \r
 408 <p class="whs5">algo.NumberOfSegments(9)</p>\r
 409 \r
 410 <p class="whs4">&nbsp;</p>\r
 411 \r
 412 <p class="whs4"># create a quadrangle 2D \r
 413  algorithm for faces</p>\r
 414 \r
 415 <p class="whs5">mesh.Quadrangle()</p>\r
 416 \r
 417 <p class="whs4">&nbsp;</p>\r
 418 \r
 419 <p class="whs4"># create a local hypothesis</p>\r
 420 \r
 421 <p class="whs5">algo = mesh.Segment(EdgesList[2])</p>\r
 422 \r
 423 <p class="whs4">&nbsp;</p>\r
 424 \r
 425 <p class="whs4"># define &quot;Arithmetic1D&quot; \r
 426  hypothesis to cut an edge in several segments with increasing arithmetic \r
 427  length </p>\r
 428 \r
 429 <p class="whs5">algo.Arithmetic1D(1, \r
 430  4)</p>\r
 431 \r
 432 <p class="whs4">&nbsp;</p>\r
 433 \r
 434 <p class="whs4"># define &quot;Propagation&quot; \r
 435  hypothesis that propagate all other hypothesis on all edges on<span style="margin-top: 0pt; margin-bottom: 0pt;"> \r
 436  the opposite side in case of quadrangular faces</span></p>\r
 437 \r
 438 <p class="whs5">algo.Propagation()</p>\r
 439 \r
 440 <p class="whs5">mesh.Compute()</p>\r
 441 \r
 442 <p class="whs4">&nbsp;</p>\r
 443 \r
 444 <p class="whs4"># sew free borders</p>\r
 445 \r
 446 <p class="whs5">anEditor = mesh.GetMesh().GetMeshEditor()</p>\r
 447 \r
 448 <p class="whs5">anEditor.SewConformFreeBorders(5, \r
 449  45, 6, 3, 24)</p>\r
 450 \r
 451 <p class="whs4"><span style="font-family: 'Lucida Console', monospace;">salome.sg.updateObjBrowser(1)</span> \r
 452  </p>\r
 453 \r
 454 <p class="whs4">&nbsp;</p>\r
 455 \r
 456 <h4><a name=bookmark6>Sew Free Borders</a></h4>\r
 457 \r
 458 <p class="whs5">import salome</p>\r
 459 \r
 460 <p class="whs5">import geompy</p>\r
 461 \r
 462 <p class="whs5">import smesh</p>\r
 463 \r
 464 <p class="whs4">&nbsp;</p>\r
 465 \r
 466 <p class="whs4"># create two faces of the \r
 467  box</p>\r
 468 \r
 469 <p class="whs5">box1 = geompy.MakeBox(0., \r
 470  0., 0., 20., 20., 15.)</p>\r
 471 \r
 472 <p class="whs5">subShapeList1 = \r
 473  geompy.SubShapeAll(box1, geompy.ShapeType[&quot;FACE&quot;])</p>\r
 474 \r
 475 <p class="whs4">&nbsp;</p>\r
 476 \r
 477 <p class="whs5">box2 = geompy.MakeBox(0., \r
 478  5., 0., 20., 20., 15.)</p>\r
 479 \r
 480 <p class="whs5">subShapeList2 = \r
 481  geompy.SubShapeAll(box2, geompy.ShapeType[&quot;FACE&quot;])</p>\r
 482 \r
 483 <p class="whs5">EdgesList = geompy.SubShapeAll(subShapeList2[ \r
 484  1 ], geompy.ShapeType[&quot;EDGE&quot;])</p>\r
 485 \r
 486 <p class="whs5">&nbsp;</p>\r
 487 \r
 488 <p class="whs5">aComp = geompy.MakeCompound( \r
 489  [subShapeList1[ 2 ], subShapeList2[ 1 ]] )</p>\r
 490 \r
 491 <p class="whs5">idComp = geompy.addToStudy( \r
 492  aComp, &quot;Two faces&quot; )</p>\r
 493 \r
 494 <p class="whs5">aCompobj &nbsp;= \r
 495  salome.IDToObject( idComp )</p>\r
 496 \r
 497 <p class="whs5">&nbsp;</p>\r
 498 \r
 499 <p class="whs4"># create a mesh on two \r
 500  faces</p>\r
 501 \r
 502 <p class="whs5">mesh = smesh.Mesh(aCompobj, \r
 503  &quot;Two faces : quadrangle mesh&quot;)</p>\r
 504 \r
 505 <p class="whs5">algo = mesh.Segment()</p>\r
 506 \r
 507 <p class="whs4">&nbsp;</p>\r
 508 \r
 509 <p class="whs4"># define &quot;NumberOfSegments&quot; \r
 510  hypothesis to cut an edge in a fixed number of segments</p>\r
 511 \r
 512 <p class="whs5">algo.NumberOfSegments(4)</p>\r
 513 \r
 514 <p class="whs4">&nbsp;</p>\r
 515 \r
 516 <p class="whs4"># creates a quadrangle \r
 517  2D algorithm for faces</p>\r
 518 \r
 519 <p class="whs5">mesh.Quadrangle()</p>\r
 520 \r
 521 <p class="whs4">&nbsp;</p>\r
 522 \r
 523 <p class="whs4"># create a local hypothesis</p>\r
 524 \r
 525 <p class="whs5">algo = mesh.Segment(EdgesList[2])</p>\r
 526 \r
 527 <p class="whs4">&nbsp;</p>\r
 528 \r
 529 <p class="whs4"># define &quot;Arithmetic1D&quot; \r
 530  hypothesis to cut an edge in several segments with &nbsp;increasing \r
 531  arithmetic length</p>\r
 532 \r
 533 <p class="whs5">algo.Arithmetic1D(1, \r
 534  4)</p>\r
 535 \r
 536 <p class="whs4">&nbsp;</p>\r
 537 \r
 538 <p class="whs4"># define &quot;Propagation&quot; \r
 539  hypothesis that propagate all other hypothesis on all edges on<span style="margin-top: 0pt; margin-bottom: 0pt;"> \r
 540  the opposite side in case of quadrangular faces</span></p>\r
 541 \r
 542 <p class="whs5">algo.Propagation()</p>\r
 543 \r
 544 <p class="whs5">mesh.Compute()</p>\r
 545 \r
 546 <p class="whs4">&nbsp;</p>\r
 547 \r
 548 <p class="whs4"># sew free borders</p>\r
 549 \r
 550 <p class="whs5">anEditor = mesh.GetMesh().GetMeshEditor()</p>\r
 551 \r
 552 <p class="whs5">anEditor.SewFreeBorders(6, \r
 553  21, 5, 1, 12, 3, 0, 0)</p>\r
 554 \r
 555 <p class="whs5">&nbsp;</p>\r
 556 \r
 557 <p class="whs4"><span style="font-family: 'Lucida Console', monospace;">salome.sg.updateObjBrowser(1)</span> \r
 558  </p>\r
 559 \r
 560 <p class="whs4">&nbsp;</p>\r
 561 \r
 562 <h4><a name=bookmark7>Sew Side Elements</a></h4>\r
 563 \r
 564 <p class="whs5">import salome</p>\r
 565 \r
 566 <p class="whs5">import geompy</p>\r
 567 \r
 568 <p class="whs5">import smesh</p>\r
 569 \r
 570 <p class="whs5">&nbsp;</p>\r
 571 \r
 572 <p class="whs4"># create two faces of the \r
 573  box</p>\r
 574 \r
 575 <p class="whs5">box1 = geompy.MakeBox(0., \r
 576  &nbsp;0., 0., \r
 577  10., 10., 10.)</p>\r
 578 \r
 579 <p class="whs5">box2 = geompy.MakeBox(0., \r
 580  15., 0., 20., 25., 10.)</p>\r
 581 \r
 582 <p class="whs5">EdgesList = geompy.SubShapeAll(box2, \r
 583  geompy.ShapeType[&quot;EDGE&quot;])</p>\r
 584 \r
 585 <p class="whs4">&nbsp;</p>\r
 586 \r
 587 <p class="whs5">aComp = geompy.MakeCompound( \r
 588  [box1, box2] )</p>\r
 589 \r
 590 <p class="whs5">idComp = geompy.addToStudy( \r
 591  aComp, &quot;Two faces&quot; )</p>\r
 592 \r
 593 <p class="whs5">aCompobj &nbsp;= \r
 594  salome.IDToObject( idComp )</p>\r
 595 \r
 596 <p class="whs5">&nbsp;</p>\r
 597 \r
 598 <p class="whs4"># create a mesh on two \r
 599  faces</p>\r
 600 \r
 601 <p class="whs5">mesh = smesh.Mesh(aCompobj, \r
 602  &quot;Two faces : quadrangle mesh&quot;)</p>\r
 603 \r
 604 <p class="whs5">algo = mesh.Segment()</p>\r
 605 \r
 606 <p class="whs4">&nbsp;</p>\r
 607 \r
 608 <p class="whs4"># define &quot;NumberOfSegments&quot; \r
 609  hypothesis to cut an edge in a fixed number of segments</p>\r
 610 \r
 611 <p class="whs5">algo.NumberOfSegments(2)</p>\r
 612 \r
 613 <p class="whs4">&nbsp;</p>\r
 614 \r
 615 <p class="whs4"># create a quadrangle 2D \r
 616  algorithm for faces</p>\r
 617 \r
 618 <p class="whs5">mesh.Quadrangle()</p>\r
 619 \r
 620 <p class="whs4">&nbsp;</p>\r
 621 \r
 622 <p class="whs4"># create a local hypothesis</p>\r
 623 \r
 624 <p class="whs5">algo = mesh.Segment(EdgesList[8])</p>\r
 625 \r
 626 <p class="whs4">&nbsp;</p>\r
 627 \r
 628 <p class="whs4"># define &quot;Arithmetic1D&quot; \r
 629  hypothesis to cut an edge in several segments with increasing arithmetic \r
 630  length </p>\r
 631 \r
 632 <p class="whs5">algo.NumberOfSegments(4)</p>\r
 633 \r
 634 <p class="whs4">&nbsp;</p>\r
 635 \r
 636 <p class="whs4"># define &quot;Propagation&quot; \r
 637  hypothesis that propagates all other hypothesis on all edges on <span \r
 638  style="margin-top: 0pt; margin-bottom: 0pt;">the opposite side in case \r
 639  of quadrangular faces</span></p>\r
 640 \r
 641 <p class="whs5">algo.Propagation()</p>\r
 642 \r
 643 <p class="whs5">mesh.Compute()</p>\r
 644 \r
 645 <p class="whs4">&nbsp;</p>\r
 646 \r
 647 <p class="whs4"># sew free borders</p>\r
 648 \r
 649 <p class="whs5">anEditor = mesh.GetMesh().GetMeshEditor()</p>\r
 650 \r
 651 <p class="whs5">anEditor.SewSideElements([69, \r
 652  70, 71, 72], [91, 92, 89, 90], 8, 38, 23, 58)</p>\r
 653 \r
 654 <p class="whs4"><span style="font-family: 'Lucida Console', monospace;">salome.sg.updateObjBrowser(1)</span> \r
 655  </p>\r
 656 \r
 657 <script type="text/javascript" language="javascript1.2">\r
 658 <!--\r
 659 if (window.writeIntopicBar)\r
 660         writeIntopicBar(0);\r
 661 //-->\r
 662 </script>\r
 663 </body>\r
 664 </html>\r