src/SMESH_SWIG/ex01_cube2build.py

   1 #  -*- coding: iso-8859-1 -*-
   2 # Copyright (C) 2007-2011  CEA/DEN, EDF R&D, OPEN CASCADE
   3 #
   4 # Copyright (C) 2003-2007  OPEN CASCADE, EADS/CCR, LIP6, CEA/DEN,
   5 # CEDRAT, EDF R&D, LEG, PRINCIPIA R&D, BUREAU VERITAS
   6 #
   7 # This library is free software; you can redistribute it and/or
   8 # modify it under the terms of the GNU Lesser General Public
   9 # License as published by the Free Software Foundation; either
  10 # version 2.1 of the License.
  11 #
  12 # This library is distributed in the hope that it will be useful,
  13 # but WITHOUT ANY WARRANTY; without even the implied warranty of
  14 # MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
  15 # Lesser General Public License for more details.
  16 #
  17 # You should have received a copy of the GNU Lesser General Public
  18 # License along with this library; if not, write to the Free Software
  19 # Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA  02111-1307 USA
  20 #
  21 # See http://www.salome-platform.org/ or email : webmaster.salome@opencascade.com
  22 #
  23
  24 # =======================================
  25 #
  26 from geompy import *
  27
  28 import smesh
  29
  30 # Geometry
  31 # ========
  32
  33 # A small cube centered and put on a great cube build by points, edges, faces and solids
  34
  35 # Points
  36 # ------
  37
  38 greatPoint111 = MakeVertex( 0,  0,  0)
  39 greatPoint211 = MakeVertex(10,  0,  0)
  40 greatPoint311 = MakeVertex(20,  0,  0)
  41 greatPoint411 = MakeVertex(30,  0,  0)
  42
  43 greatPoint121 = MakeVertex( 0, 10,  0)
  44 greatPoint221 = MakeVertex(10, 10,  0)
  45 greatPoint321 = MakeVertex(20, 10,  0)
  46 greatPoint421 = MakeVertex(30, 10,  0)
  47
  48 greatPoint112 = MakeVertex( 0,  0, 10)
  49 greatPoint212 = MakeVertex(10,  0, 10)
  50 greatPoint312 = MakeVertex(20,  0, 10)
  51 greatPoint412 = MakeVertex(30,  0, 10)
  52
  53 greatPoint122 = MakeVertex( 0, 10, 10)
  54 greatPoint222 = MakeVertex(10, 10, 10)
  55 greatPoint322 = MakeVertex(20, 10, 10)
  56 greatPoint422 = MakeVertex(30, 10, 10)
  57
  58 greatPoint113 = MakeVertex( 0,  0, 20)
  59 greatPoint213 = MakeVertex(10,  0, 20)
  60 greatPoint313 = MakeVertex(20,  0, 20)
  61 greatPoint413 = MakeVertex(30,  0, 20)
  62
  63 greatPoint123 = MakeVertex( 0, 10, 20)
  64 greatPoint223 = MakeVertex(10, 10, 20)
  65 greatPoint323 = MakeVertex(20, 10, 20)
  66 greatPoint423 = MakeVertex(30, 10, 20)
  67
  68 greatPoint114 = MakeVertex( 0,  0, 30)
  69 greatPoint214 = MakeVertex(10,  0, 30)
  70 greatPoint314 = MakeVertex(20,  0, 30)
  71 greatPoint414 = MakeVertex(30,  0, 30)
  72
  73 greatPoint124 = MakeVertex( 0, 10, 30)
  74 greatPoint224 = MakeVertex(10, 10, 30)
  75 greatPoint324 = MakeVertex(20, 10, 30)
  76 greatPoint424 = MakeVertex(30, 10, 30)
  77
  78
  79 smallPoint111 = greatPoint222
  80 smallPoint211 = greatPoint322
  81 smallPoint121 = MakeVertex(10, 20, 10)
  82 smallPoint221 = MakeVertex(20, 20, 10)
  83
  84 smallPoint112 = greatPoint223
  85 smallPoint212 = greatPoint323
  86 smallPoint122 = MakeVertex(10, 20, 20)
  87 smallPoint222 = MakeVertex(20, 20, 20)
  88
  89 # Edges
  90 # -----
  91
  92 smallEdgeX11 = MakeEdge(smallPoint111, smallPoint211)
  93 smallEdgeX21 = MakeEdge(smallPoint121, smallPoint221)
  94 smallEdgeX12 = MakeEdge(smallPoint112, smallPoint212)
  95 smallEdgeX22 = MakeEdge(smallPoint122, smallPoint222)
  96
  97 smallEdgeY11 = MakeEdge(smallPoint111, smallPoint121)
  98 smallEdgeY21 = MakeEdge(smallPoint211, smallPoint221)
  99 smallEdgeY12 = MakeEdge(smallPoint112, smallPoint122)
 100 smallEdgeY22 = MakeEdge(smallPoint212, smallPoint222)
 101
 102 smallEdgeZ11 = MakeEdge(smallPoint111, smallPoint112)
 103 smallEdgeZ21 = MakeEdge(smallPoint211, smallPoint212)
 104 smallEdgeZ12 = MakeEdge(smallPoint121, smallPoint122)
 105 smallEdgeZ22 = MakeEdge(smallPoint221, smallPoint222)
 106
 107
 108 greatEdgeX111 = MakeEdge(greatPoint111, greatPoint211)
 109 greatEdgeX211 = MakeEdge(greatPoint211, greatPoint311)
 110 greatEdgeX311 = MakeEdge(greatPoint311, greatPoint411)
 111 greatEdgeX121 = MakeEdge(greatPoint121, greatPoint221)
 112 greatEdgeX221 = MakeEdge(greatPoint221, greatPoint321)
 113 greatEdgeX321 = MakeEdge(greatPoint321, greatPoint421)
 114
 115 greatEdgeX112 = MakeEdge(greatPoint112, greatPoint212)
 116 greatEdgeX212 = MakeEdge(greatPoint212, greatPoint312)
 117 greatEdgeX312 = MakeEdge(greatPoint312, greatPoint412)
 118 greatEdgeX122 = MakeEdge(greatPoint122, greatPoint222)
 119 greatEdgeX222 = smallEdgeX11
 120 greatEdgeX322 = MakeEdge(greatPoint322, greatPoint422)
 121
 122 greatEdgeX113 = MakeEdge(greatPoint113, greatPoint213)
 123 greatEdgeX213 = MakeEdge(greatPoint213, greatPoint313)
 124 greatEdgeX313 = MakeEdge(greatPoint313, greatPoint413)
 125 greatEdgeX123 = MakeEdge(greatPoint123, greatPoint223)
 126 greatEdgeX223 = smallEdgeX12
 127 greatEdgeX323 = MakeEdge(greatPoint323, greatPoint423)
 128
 129 greatEdgeX114 = MakeEdge(greatPoint114, greatPoint214)
 130 greatEdgeX214 = MakeEdge(greatPoint214, greatPoint314)
 131 greatEdgeX314 = MakeEdge(greatPoint314, greatPoint414)
 132 greatEdgeX124 = MakeEdge(greatPoint124, greatPoint224)
 133 greatEdgeX224 = MakeEdge(greatPoint224, greatPoint324)
 134 greatEdgeX324 = MakeEdge(greatPoint324, greatPoint424)
 135
 136 greatEdgeY11 = MakeEdge(greatPoint111, greatPoint121)
 137 greatEdgeY21 = MakeEdge(greatPoint211, greatPoint221)
 138 greatEdgeY31 = MakeEdge(greatPoint311, greatPoint321)
 139 greatEdgeY41 = MakeEdge(greatPoint411, greatPoint421)
 140
 141 greatEdgeY12 = MakeEdge(greatPoint112, greatPoint122)
 142 greatEdgeY22 = MakeEdge(greatPoint212, greatPoint222)
 143 greatEdgeY32 = MakeEdge(greatPoint312, greatPoint322)
 144 greatEdgeY42 = MakeEdge(greatPoint412, greatPoint422)
 145
 146 greatEdgeY13 = MakeEdge(greatPoint113, greatPoint123)
 147 greatEdgeY23 = MakeEdge(greatPoint213, greatPoint223)
 148 greatEdgeY33 = MakeEdge(greatPoint313, greatPoint323)
 149 greatEdgeY43 = MakeEdge(greatPoint413, greatPoint423)
 150
 151 greatEdgeY14 = MakeEdge(greatPoint114, greatPoint124)
 152 greatEdgeY24 = MakeEdge(greatPoint214, greatPoint224)
 153 greatEdgeY34 = MakeEdge(greatPoint314, greatPoint324)
 154 greatEdgeY44 = MakeEdge(greatPoint414, greatPoint424)
 155
 156 greatEdgeZ111 = MakeEdge(greatPoint111, greatPoint112)
 157 greatEdgeZ211 = MakeEdge(greatPoint211, greatPoint212)
 158 greatEdgeZ311 = MakeEdge(greatPoint311, greatPoint312)
 159 greatEdgeZ411 = MakeEdge(greatPoint411, greatPoint412)
 160
 161 greatEdgeZ121 = MakeEdge(greatPoint121, greatPoint122)
 162 greatEdgeZ221 = MakeEdge(greatPoint221, greatPoint222)
 163 greatEdgeZ321 = MakeEdge(greatPoint321, greatPoint322)
 164 greatEdgeZ421 = MakeEdge(greatPoint421, greatPoint422)
 165
 166 greatEdgeZ112 = MakeEdge(greatPoint112, greatPoint113)
 167 greatEdgeZ212 = MakeEdge(greatPoint212, greatPoint213)
 168 greatEdgeZ312 = MakeEdge(greatPoint312, greatPoint313)
 169 greatEdgeZ412 = MakeEdge(greatPoint412, greatPoint413)
 170
 171 greatEdgeZ122 = MakeEdge(greatPoint122, greatPoint123)
 172 greatEdgeZ222 = smallEdgeZ11
 173 greatEdgeZ322 = smallEdgeZ21
 174 greatEdgeZ422 = MakeEdge(greatPoint422, greatPoint423)
 175
 176 greatEdgeZ113 = MakeEdge(greatPoint113, greatPoint114)
 177 greatEdgeZ213 = MakeEdge(greatPoint213, greatPoint214)
 178 greatEdgeZ313 = MakeEdge(greatPoint313, greatPoint314)
 179 greatEdgeZ413 = MakeEdge(greatPoint413, greatPoint414)
 180
 181 greatEdgeZ123 = MakeEdge(greatPoint123, greatPoint124)
 182 greatEdgeZ223 = MakeEdge(greatPoint223, greatPoint224)
 183 greatEdgeZ323 = MakeEdge(greatPoint323, greatPoint324)
 184 greatEdgeZ423 = MakeEdge(greatPoint423, greatPoint424)
 185
 186 # Faces
 187 # -----
 188
 189 smallFaceX1 = MakeQuad(smallEdgeY11, smallEdgeZ11, smallEdgeY12, smallEdgeZ12)
 190 smallFaceX2 = MakeQuad(smallEdgeY21, smallEdgeZ21, smallEdgeY22, smallEdgeZ22)
 191 smallFaceY1 = MakeQuad(smallEdgeX11, smallEdgeZ11, smallEdgeX12, smallEdgeZ21)
 192 smallFaceY2 = MakeQuad(smallEdgeX21, smallEdgeZ12, smallEdgeX22, smallEdgeZ22)
 193 smallFaceZ1 = MakeQuad(smallEdgeX11, smallEdgeY11, smallEdgeX21, smallEdgeY21)
 194 smallFaceZ2 = MakeQuad(smallEdgeX12, smallEdgeY12, smallEdgeX22, smallEdgeY22)
 195
 196
 197 greatFaceX11 = MakeQuad(greatEdgeY11, greatEdgeZ111, greatEdgeY12, greatEdgeZ121)
 198 greatFaceX21 = MakeQuad(greatEdgeY21, greatEdgeZ211, greatEdgeY22, greatEdgeZ221)
 199 greatFaceX31 = MakeQuad(greatEdgeY31, greatEdgeZ311, greatEdgeY32, greatEdgeZ321)
 200 greatFaceX41 = MakeQuad(greatEdgeY41, greatEdgeZ411, greatEdgeY42, greatEdgeZ421)
 201
 202 greatFaceX12 = MakeQuad(greatEdgeY12, greatEdgeZ112, greatEdgeY13, greatEdgeZ122)
 203 greatFaceX22 = MakeQuad(greatEdgeY22, greatEdgeZ212, greatEdgeY23, greatEdgeZ222)
 204 greatFaceX32 = MakeQuad(greatEdgeY32, greatEdgeZ312, greatEdgeY33, greatEdgeZ322)
 205 greatFaceX42 = MakeQuad(greatEdgeY42, greatEdgeZ412, greatEdgeY43, greatEdgeZ422)
 206
 207 greatFaceX13 = MakeQuad(greatEdgeY13, greatEdgeZ113, greatEdgeY14, greatEdgeZ123)
 208 greatFaceX23 = MakeQuad(greatEdgeY23, greatEdgeZ213, greatEdgeY24, greatEdgeZ223)
 209 greatFaceX33 = MakeQuad(greatEdgeY33, greatEdgeZ313, greatEdgeY34, greatEdgeZ323)
 210 greatFaceX43 = MakeQuad(greatEdgeY43, greatEdgeZ413, greatEdgeY44, greatEdgeZ423)
 211
 212 greatFaceY111 = MakeQuad(greatEdgeX111, greatEdgeZ111, greatEdgeX112, greatEdgeZ211)
 213 greatFaceY211 = MakeQuad(greatEdgeX211, greatEdgeZ211, greatEdgeX212, greatEdgeZ311)
 214 greatFaceY311 = MakeQuad(greatEdgeX311, greatEdgeZ311, greatEdgeX312, greatEdgeZ411)
 215 greatFaceY121 = MakeQuad(greatEdgeX121, greatEdgeZ121, greatEdgeX122, greatEdgeZ221)
 216 greatFaceY221 = MakeQuad(greatEdgeX221, greatEdgeZ221, greatEdgeX222, greatEdgeZ321)
 217 greatFaceY321 = MakeQuad(greatEdgeX321, greatEdgeZ321, greatEdgeX322, greatEdgeZ421)
 218
 219 greatFaceY112 = MakeQuad(greatEdgeX112, greatEdgeZ112, greatEdgeX113, greatEdgeZ212)
 220 greatFaceY212 = MakeQuad(greatEdgeX212, greatEdgeZ212, greatEdgeX213, greatEdgeZ312)
 221 greatFaceY312 = MakeQuad(greatEdgeX312, greatEdgeZ312, greatEdgeX313, greatEdgeZ412)
 222 greatFaceY122 = MakeQuad(greatEdgeX122, greatEdgeZ122, greatEdgeX123, greatEdgeZ222)
 223 greatFaceY222 = smallFaceY1
 224 greatFaceY322 = MakeQuad(greatEdgeX322, greatEdgeZ322, greatEdgeX323, greatEdgeZ422)
 225
 226 greatFaceY113 = MakeQuad(greatEdgeX113, greatEdgeZ113, greatEdgeX114, greatEdgeZ213)
 227 greatFaceY213 = MakeQuad(greatEdgeX213, greatEdgeZ213, greatEdgeX214, greatEdgeZ313)
 228 greatFaceY313 = MakeQuad(greatEdgeX313, greatEdgeZ313, greatEdgeX314, greatEdgeZ413)
 229 greatFaceY123 = MakeQuad(greatEdgeX123, greatEdgeZ123, greatEdgeX124, greatEdgeZ223)
 230 greatFaceY223 = MakeQuad(greatEdgeX223, greatEdgeZ223, greatEdgeX224, greatEdgeZ323)
 231 greatFaceY323 = MakeQuad(greatEdgeX323, greatEdgeZ323, greatEdgeX324, greatEdgeZ423)
 232
 233 greatFaceZ11 = MakeQuad(greatEdgeX111, greatEdgeY11, greatEdgeX121, greatEdgeY21)
 234 greatFaceZ21 = MakeQuad(greatEdgeX211, greatEdgeY21, greatEdgeX221, greatEdgeY31)
 235 greatFaceZ31 = MakeQuad(greatEdgeX311, greatEdgeY31, greatEdgeX321, greatEdgeY41)
 236
 237 greatFaceZ12 = MakeQuad(greatEdgeX112, greatEdgeY12, greatEdgeX122, greatEdgeY22)
 238 greatFaceZ22 = MakeQuad(greatEdgeX212, greatEdgeY22, greatEdgeX222, greatEdgeY32)
 239 greatFaceZ32 = MakeQuad(greatEdgeX312, greatEdgeY32, greatEdgeX322, greatEdgeY42)
 240
 241 greatFaceZ13 = MakeQuad(greatEdgeX113, greatEdgeY13, greatEdgeX123, greatEdgeY23)
 242 greatFaceZ23 = MakeQuad(greatEdgeX213, greatEdgeY23, greatEdgeX223, greatEdgeY33)
 243 greatFaceZ33 = MakeQuad(greatEdgeX313, greatEdgeY33, greatEdgeX323, greatEdgeY43)
 244
 245 greatFaceZ14 = MakeQuad(greatEdgeX114, greatEdgeY14, greatEdgeX124, greatEdgeY24)
 246 greatFaceZ24 = MakeQuad(greatEdgeX214, greatEdgeY24, greatEdgeX224, greatEdgeY34)
 247 greatFaceZ34 = MakeQuad(greatEdgeX314, greatEdgeY34, greatEdgeX324, greatEdgeY44)
 248
 249 # Solids
 250 # ------
 251
 252 smallBlock   = MakeHexa(smallFaceX1, smallFaceX2, smallFaceY1, smallFaceY2, smallFaceZ1, smallFaceZ2)
 253
 254 greatBlock11 = MakeHexa(greatFaceX11, greatFaceX21, greatFaceY111, greatFaceY121, greatFaceZ11, greatFaceZ12)
 255 greatBlock21 = MakeHexa(greatFaceX21, greatFaceX31, greatFaceY211, greatFaceY221, greatFaceZ21, greatFaceZ22)
 256 greatBlock31 = MakeHexa(greatFaceX31, greatFaceX41, greatFaceY311, greatFaceY321, greatFaceZ31, greatFaceZ32)
 257
 258 greatBlock12 = MakeHexa(greatFaceX12, greatFaceX22, greatFaceY112, greatFaceY122, greatFaceZ12, greatFaceZ13)
 259 greatBlock22 = MakeHexa(greatFaceX22, greatFaceX32, greatFaceY212, greatFaceY222, greatFaceZ22, greatFaceZ23)
 260 greatBlock32 = MakeHexa(greatFaceX32, greatFaceX42, greatFaceY312, greatFaceY322, greatFaceZ32, greatFaceZ33)
 261
 262 greatBlock13 = MakeHexa(greatFaceX13, greatFaceX23, greatFaceY113, greatFaceY123, greatFaceZ13, greatFaceZ14)
 263 greatBlock23 = MakeHexa(greatFaceX23, greatFaceX33, greatFaceY213, greatFaceY223, greatFaceZ23, greatFaceZ24)
 264 greatBlock33 = MakeHexa(greatFaceX33, greatFaceX43, greatFaceY313, greatFaceY323, greatFaceZ33, greatFaceZ34)
 265
 266 # Compound
 267 # --------
 268
 269 c_l = []
 270 c_l.append(smallBlock)
 271 c_l.append(greatBlock11)
 272 c_l.append(greatBlock21)
 273 c_l.append(greatBlock31)
 274 c_l.append(greatBlock12)
 275 c_l.append(greatBlock22)
 276 c_l.append(greatBlock32)
 277 c_l.append(greatBlock13)
 278 c_l.append(greatBlock23)
 279 c_l.append(greatBlock33)
 280
 281 piece = MakeCompound(c_l)
 282
 283 # Add in study
 284 # ------------
 285
 286 piece_id = addToStudy(piece, "ex01_cube2build")
 287
 288 # Meshing
 289 # =======
 290
 291 smesh.SetCurrentStudy(salome.myStudy)
 292
 293 # Create hexahedrical mesh on piece
 294 # ---------------------------------
 295
 296 hexa = smesh.Mesh(piece, "ex01_cube2build:hexa")
 297
 298 algo = hexa.Segment()
 299 algo.NumberOfSegments(4)
 300
 301 hexa.Quadrangle()
 302
 303 hexa.Hexahedron()
 304
 305 # Create local hypothesis
 306 # -----------------------
 307
 308 algo = hexa.Segment(greatEdgeX111)
 309
 310 algo.Arithmetic1D(1, 4)
 311
 312 algo.Propagation()
 313
 314 # Compute the mesh
 315 # ----------------
 316
 317 hexa.Compute()