doc/examples/defining_hypotheses_adaptive1d.py

   1 # Usage of Adaptive 1D hypothesis
   2
   3 import math
   4
   5 import salome
   6 salome.salome_init_without_session()
   7
   8 from salome.geom import geomBuilder
   9 from salome.smesh import smeshBuilder
  10
  11 geom_builder = geomBuilder.New()
  12 smesh_builder = smeshBuilder.New()
  13
  14 box   = geom_builder.MakeBoxDXDYDZ( 100, 100, 100 )
  15 tool  = geom_builder.MakeTranslation( box, 50, 0, 10 )
  16 axis  = geom_builder.MakeVector( geom_builder.MakeVertex( 100, 0, 100 ),geom_builder.MakeVertex( 100, 10, 100 ),)
  17 tool  = geom_builder.Rotate( tool, axis, math.pi * 25 / 180. )
  18 shape = geom_builder.MakeCut( box, tool )
  19 cyl   = geom_builder.MakeCylinder( geom_builder.MakeVertex( -10,5, 95 ), geom_builder.MakeVectorDXDYDZ(1,0,0), 2, 90)
  20 shape = geom_builder.MakeCut( shape, cyl )
  21 tool  = geom_builder.MakeBoxTwoPnt( geom_builder.MakeVertex( -10, 2, 15 ), geom_builder.MakeVertex( 90, 5, 16 ))
  22 shape = geom_builder.MakeCut( shape, tool, theName="shape" )
  23
  24 # Parameters of Adaptive hypothesis. minSize and maxSize are such that they do not limit
  25 # size of segments because size of geometrical features lies within [2.-100.] range, hence
  26 # size of segments is defined by deflection parameter and size of geometrical features only.
  27 minSize = 0.1
  28 maxSize = 200
  29 deflection = 0.05
  30
  31 mesh = smesh_builder.Mesh( shape )
  32 mesh.Segment().Adaptive( minSize, maxSize, deflection )
  33 mesh.Triangle().MaxElementArea( 300 )
  34 if not mesh.Compute(): raise Exception("Error when computing Mesh")