Python CSGCGALDomain3D示例

示例#1

文件： meshes.py 项目： janblechschmidt/NonvarFEM

def mesh_Cube(n=10):
    dom = mshr.Box(Point(0., 0., 0.), Point(1., 1., 1.))
    gen = mshr.CSGCGALMeshGenerator3D()
    #gen.parameters["facet_angle"] = 30.0
    #gen.parameters["facet_size"] = 0.5
    #gen.parameters["edge_size"] = 0.5
    gen.parameters["mesh_resolution"] = 0.01
    gen.parameters["odt_optimize"] = True
    mesh = gen.generate(mshr.CSGCGALDomain3D(dom))
    return mesh

示例#2

imp = implicitModeller(
    mergeActors(tf, bx),
    distance=0.04,
    outer=True,
    #                       res=[50,20,10],
    res=[110, 40, 20],
    bounds=[-1.0, 10.0, -1.0, 2.0, -.5, .5],
    maxdist=0.25,
)

imp.write(mytext + '.stl', binary=False)
show(imp)

#### tetralize
import mshr
import dolfin

surface = mshr.Surface3D(mytext + '.stl')
polyhedral_domain = mshr.CSGCGALDomain3D(surface)
dolfin.info(polyhedral_domain, True)

generator = mshr.CSGCGALMeshGenerator3D()
generator.parameters["mesh_resolution"] = 40.0
mesh = generator.generate(polyhedral_domain)

dolfin.File(mytext + ".xml") << mesh
printc('saved', mytext + ".xml")

#from vtkplotter.dolfin import plot
#plot(mesh)

示例#3

文件： ex02_tetralize-mesh.py 项目： zlinzju/vtkplotter

> vtkplotter shuttle.xml
"""
import mshr
import dolfin
from vtkplotter.dolfin import datadir, plot

fname = "shuttle.stl"

surface = mshr.Surface3D(datadir + fname)

# add a cylinder
cyl = mshr.Cylinder(dolfin.Point(-1.4, 0, 0), dolfin.Point(-1.0, 0, 0), 0.5,
                    0.7)
totdomain = surface + cyl

polyhedral_domain = mshr.CSGCGALDomain3D(totdomain)
dolfin.info(polyhedral_domain, True)

generator = mshr.CSGCGALMeshGenerator3D()

#### Try automatic
generator.parameters["mesh_resolution"] = 35.0
mesh = generator.generate(polyhedral_domain)

xmlname = fname.split(".")[0] + ".xml"
dolfin.File(xmlname) << mesh
print(mesh, "saved to " + xmlname)

##################################
plot(mesh, text=__doc__)
##################################

示例#4

文件： remove_degenerate.py 项目： Nasrollah/mshr

from __future__ import print_function
import mshr
import dolfin

TOLERANCE = 1e-10

# This geometry generates a lot of degenerate facets
cone = mshr.Cylinder(dolfin.Point(-1.0, 1.0, 1.0), dolfin.Point(1.0, -1.0, -1.0), .5, .5)
cyl = mshr.Cone(dolfin.Point(1.0, -1.0, 1.0), dolfin.Point(-1.0, 1.0, -1.0), .5)
geometry = cone + cyl;

polyhedral_domain = mshr.CSGCGALDomain3D(geometry)

print("Degenerate facets after boolean operation: {0}".format(polyhedral_domain.num_degenerate_facets(TOLERANCE)))
polyhedral_domain.remove_degenerate_facets(TOLERANCE)

dolfin.info(polyhedral_domain, True)

示例#5

文件： mwe_3d_mesh_refinement.py 项目： janblechschmidt/NonvarFEM

#     print('Dim V: {}'.format(V.dim()))
#     if k < K:
#         m = refine(m)

dom = mshr.Box(Point(0., 0., 0.), Point(1., 1., 1.))
# m_orig = mshr.generate_mesh(dom,0.01)

gen = mshr.CSGCGALMeshGenerator3D()

#gen.parameters["facet_angle"] = 30.0
#gen.parameters["facet_size"] = 0.5
#gen.parameters["edge_size"] = 0.5
gen.parameters["mesh_resolution"] = 0.01
gen.parameters["exude_optimize"] = True

m_orig = gen.generate(mshr.CSGCGALDomain3D(dom))

xf = File('mesh_orig.pvd')
xf << m_orig

V = FunctionSpace(m_orig, 'DG', 0)
w = TestFunction(V)
s = assemble(w * dx)
print(s.get_local())

# m = refine(m_orig)
#
# V = FunctionSpace(m, 'Lagrange', 2)
# u = TrialFunction(V)
# v = TestFunction(V)
#

示例#6

  vertex 0 1 0
  vertex 0 1 1
endloop
endfacet
facet normal 0 0 0
outer loop
  vertex 0 0 1
  vertex 0 0 0
  vertex 1 0 1
endloop
endfacet
facet normal 0 0 0
outer loop
  vertex 0 0 0
  vertex 1 0 0
  vertex 1 0 1
endloop
endfacet
endsolid
"""

fd, filename = tempfile.mkstemp(suffix=".stl")
os.write(fd, cube)

s = mshr.Surface3D(filename)
d = mshr.CSGCGALDomain3D(s)
os.close(fd)
os.remove(filename)

assert d.num_holes() == 0

示例#7

文件： concentric_cylinders.py 项目： zzyatlantise/fenics-tutorial

from dolfin import *
import mshr

# Parameters for geometry
a = 0.04
b = a + 0.004
c = a + 0.01
L = 0.5

# Define cylinders
inner = mshr.CSGCGALDomain3D(
    mshr.Cylinder(Point(0, 0, 0), Point(0, 0, L), a, a))
mid = mshr.CSGCGALDomain3D(mshr.Cylinder(Point(0, 0, 0), Point(0, 0, L), b, b))
outer = mshr.CSGCGALDomain3D(
    mshr.Cylinder(Point(0, 0, 0), Point(0, 0, L), c, c))

generator = mshr.TetgenMeshGenerator3D()
generator.parameters["preserve_surface"] = True
generator.parameters["mesh_resolution"] = 16.

# Mesh inner cylinder
inner_mesh = generator.generate(inner)

# Mesh mid part
mid_mesh = generator.generate(mshr.CSGCGALDomain3D(mid - inner))

# Mesh outer part
outer_mesh = generator.generate(mshr.CSGCGALDomain3D(outer - mid))

# Glue together inner and mid mesh
inner_mid_mesh = mshr.DolfinMeshUtils.merge_meshes(inner_mesh, mid_mesh)