Mod: simple structure compressed visually
This commit is contained in:
parent
cae96521b0
commit
eca4670b07
@ -7,38 +7,44 @@
|
||||
|
||||
from math import pi, sqrt
|
||||
|
||||
def simpleCubic(theta = 0.01, fillet = False, direction = [1, 0, 0]):
|
||||
class simple(object):
|
||||
def __init__(self, direction, theta, r0, L, radius, filletsEnabled, fillets):
|
||||
|
||||
self.direction = direction
|
||||
self.theta = theta
|
||||
self.r0 = r0
|
||||
self.L = L
|
||||
self.radius = radius
|
||||
self.filletsEnabled = filletsEnabled
|
||||
self.fillets = fillets
|
||||
|
||||
|
||||
def build(self):
|
||||
|
||||
###
|
||||
# Pore Cell
|
||||
##
|
||||
if self.direction in [[1, 0, 0], [0, 0, 1]]:
|
||||
###
|
||||
# Parameters
|
||||
##
|
||||
r0 = 1.0
|
||||
L = 2 * r0
|
||||
|
||||
radius = r0 / (1 - theta)
|
||||
xn, yn, zn = 3, 3, 3
|
||||
|
||||
length = L * sqrt(2)
|
||||
width = L * sqrt(2)
|
||||
height = L
|
||||
length = self.L * sqrt(2)
|
||||
width = self.L * sqrt(2)
|
||||
height = self.L
|
||||
|
||||
xl = sqrt(length ** 2 * 0.5)
|
||||
yw = xl
|
||||
zh = height
|
||||
|
||||
C1, C2 = 0.8, 0.5 #0.8, 0.05
|
||||
theta1, theta2 = 0.01, 0.28
|
||||
Cf = C1 + (C2 - C1) / (theta2 - theta1) * (theta - theta1)
|
||||
delta = 0.2
|
||||
filletradius = delta - Cf * (radius - r0)
|
||||
|
||||
scale = 100
|
||||
oo = geompy.MakeVertex(0, 0, 0)
|
||||
|
||||
###
|
||||
# Bounding box
|
||||
##
|
||||
if direction == [1, 0, 0]:
|
||||
if self.direction == [1, 0, 0]:
|
||||
sk = geompy.Sketcher3D()
|
||||
sk.addPointsAbsolute(xl, 0, 0)
|
||||
sk.addPointsAbsolute(0, yw, 0)
|
||||
@ -48,9 +54,9 @@ def simpleCubic(theta = 0.01, fillet = False, direction = [1, 0, 0]):
|
||||
|
||||
inletface = geompy.MakeFaceWires([sk.wire()], True)
|
||||
vecflow = geompy.GetNormal(inletface)
|
||||
cubic = geompy.MakePrismVecH(inletface, vecflow, width)
|
||||
poreCell = geompy.MakePrismVecH(inletface, vecflow, width)
|
||||
|
||||
elif direction == [0, 0, 1]:
|
||||
elif self.direction == [0, 0, 1]:
|
||||
sk = geompy.Sketcher3D()
|
||||
sk.addPointsAbsolute(0, yw, 0)
|
||||
sk.addPointsAbsolute(xl, 0, 0)
|
||||
@ -60,13 +66,13 @@ def simpleCubic(theta = 0.01, fillet = False, direction = [1, 0, 0]):
|
||||
|
||||
inletface = geompy.MakeFaceWires([sk.wire()], True)
|
||||
vecflow = geompy.GetNormal(inletface)
|
||||
cubic = geompy.MakePrismVecH(inletface, vecflow, height)
|
||||
poreCell = geompy.MakePrismVecH(inletface, vecflow, height)
|
||||
|
||||
|
||||
inletface = geompy.MakeScaleTransform(inletface, oo, scale)
|
||||
cubic = geompy.MakeScaleTransform(cubic, oo, scale)
|
||||
poreCell = geompy.MakeScaleTransform(poreCell, oo, scale)
|
||||
|
||||
faces = geompy.ExtractShapes(cubic, geompy.ShapeType["FACE"], False)
|
||||
faces = geompy.ExtractShapes(poreCell, geompy.ShapeType["FACE"], False)
|
||||
symetryface = []
|
||||
|
||||
for face in faces:
|
||||
@ -79,107 +85,25 @@ def simpleCubic(theta = 0.01, fillet = False, direction = [1, 0, 0]):
|
||||
else:
|
||||
symetryface.append(face)
|
||||
|
||||
###
|
||||
# Grains
|
||||
##
|
||||
ox = geompy.MakeVectorDXDYDZ(1, 0, 0)
|
||||
oy = geompy.MakeVectorDXDYDZ(0, 1, 0)
|
||||
oz = geompy.MakeVectorDXDYDZ(0, 0, 1)
|
||||
|
||||
grain = geompy.MakeSphereR(radius)
|
||||
lattice = geompy.MakeMultiTranslation2D(grain, ox, L, xn, oy, L, yn)
|
||||
lattice = geompy.MakeMultiTranslation1D(lattice, oz, L, zn)
|
||||
|
||||
grains = geompy.ExtractShapes(lattice, geompy.ShapeType["SOLID"], True)
|
||||
grains = geompy.MakeFuseList(grains, False, False)
|
||||
|
||||
grains = geompy.MakeScaleTransform(grains, oo, scale)
|
||||
grainsOrigin = None
|
||||
|
||||
if fillet:
|
||||
grainsOrigin = geompy.MakeScaleTransform(grains, oo, 1 / scale)
|
||||
grains = geompy.MakeFilletAll(grains, filletradius * scale)
|
||||
|
||||
###
|
||||
# Groups
|
||||
##
|
||||
shape = geompy.MakeCutList(cubic, [grains])
|
||||
shape = geompy.MakeScaleTransform(shape, oo, 1 / scale, theName = "simpleCubic")
|
||||
|
||||
sall = geompy.CreateGroup(shape, geompy.ShapeType["FACE"])
|
||||
geompy.UnionIDs(sall,
|
||||
geompy.SubShapeAllIDs(shape, geompy.ShapeType["FACE"]))
|
||||
|
||||
inlet = geompy.CreateGroup(shape, geompy.ShapeType["FACE"], theName = "inlet")
|
||||
inletshape = geompy.MakeCutList(inletface, [grains])
|
||||
inletshape = geompy.MakeScaleTransform(inletshape, oo, 1 / scale)
|
||||
geompy.UnionList(inlet, geompy.SubShapeAll(
|
||||
geompy.GetInPlace(shape, inletshape, True), geompy.ShapeType["FACE"]))
|
||||
|
||||
outlet = geompy.CreateGroup(shape, geompy.ShapeType["FACE"], theName = "outlet")
|
||||
outletshape = geompy.MakeCutList(outletface, [grains])
|
||||
outletshape = geompy.MakeScaleTransform(outletshape, oo, 1 / scale)
|
||||
geompy.UnionList(outlet, geompy.SubShapeAll(
|
||||
geompy.GetInPlace(shape, outletshape, True), geompy.ShapeType["FACE"]))
|
||||
|
||||
symetry = []
|
||||
for (n, face) in enumerate(symetryface):
|
||||
name = "symetry" + str(n)
|
||||
symetry.append(geompy.CreateGroup(shape, geompy.ShapeType["FACE"], theName = name))
|
||||
symetryshape = geompy.MakeCutList(face, [grains])
|
||||
symetryshape = geompy.MakeScaleTransform(symetryshape, oo, 1 / scale)
|
||||
geompy.UnionList(symetry[n], geompy.SubShapeAll(
|
||||
geompy.GetInPlace(shape, symetryshape, True), geompy.ShapeType["FACE"]))
|
||||
|
||||
groups = []
|
||||
groups.append(inlet)
|
||||
groups.append(outlet)
|
||||
groups.extend(symetry)
|
||||
|
||||
if fillet:
|
||||
strips = geompy.CreateGroup(shape, geompy.ShapeType["FACE"], theName = "strips")
|
||||
shapeShell = geompy.ExtractShapes(shape, geompy.ShapeType["SHELL"], True)
|
||||
stripsShape = geompy.MakeCutList(shapeShell[0], groups + [grainsOrigin])
|
||||
geompy.UnionList(strips, geompy.SubShapeAll(
|
||||
geompy.GetInPlace(shape, stripsShape, True), geompy.ShapeType["FACE"]))
|
||||
groups.append(strips)
|
||||
|
||||
wall = geompy.CutListOfGroups([sall], groups, theName = "wall")
|
||||
groups.append(wall)
|
||||
|
||||
return shape, groups
|
||||
|
||||
|
||||
def simpleHexagonalPrism(theta = 0.01, fillet = False, direction = [1, 1, 1]):
|
||||
|
||||
elif self.direction == [1, 1, 1]:
|
||||
###
|
||||
# Parameters
|
||||
##
|
||||
r0 = 1.0
|
||||
L = 2 * r0
|
||||
|
||||
radius = r0 / (1 - theta)
|
||||
xn, yn, zn = 3, 3, 3
|
||||
|
||||
length = L * sqrt(2)
|
||||
width = L * sqrt(2)
|
||||
height = L
|
||||
length = self.L * sqrt(2)
|
||||
width = self.L * sqrt(2)
|
||||
height = self.L
|
||||
|
||||
point = []
|
||||
xl, yw, zh = -L - L / 6, -L - L / 6, -L / 6
|
||||
point.append((L + xl, L + yw, L + zh))
|
||||
point.append((5 * L / 3 + xl, 2 * L / 3 + yw, 2 * L / 3 + zh))
|
||||
point.append((2 * L + xl, L + yw, 0 + zh))
|
||||
point.append((5 * L / 3 + xl, 5 * L / 3 + yw, -L / 3 + zh))
|
||||
point.append((L + xl, 2 * L + yw, 0 + zh))
|
||||
point.append((2 * L / 3 + xl, 5 * L / 3 + yw, 2 * L / 3 + zh))
|
||||
point.append((L + xl, L + yw, L + zh))
|
||||
|
||||
C1, C2 = 0.8, 0.5 # 0.8, 0.05
|
||||
theta1, theta2 = 0.01, 0.28
|
||||
Cf = C1 + (C2 - C1) / (theta2 - theta1) * (theta - theta1)
|
||||
delta = 0.2
|
||||
filletradius = delta - Cf * (radius - r0)
|
||||
xl, yw, zh = -self.L - self.L / 6, -self.L - self.L / 6, -self.L / 6
|
||||
point.append((L + xl, self.L + yw, self.L + zh))
|
||||
point.append((5 * self.L / 3 + xl, 2 * self.L / 3 + yw, 2 * self.L / 3 + zh))
|
||||
point.append((2 * self.L + xl, self.L + yw, 0 + zh))
|
||||
point.append((5 * self.L / 3 + xl, 5 * self.L / 3 + yw, -self.L / 3 + zh))
|
||||
point.append((L + xl, 2 * self.L + yw, 0 + zh))
|
||||
point.append((2 * self.L / 3 + xl, 5 * self.L / 3 + yw, 2 * self.L / 3 + zh))
|
||||
point.append((L + xl, self.L + yw, self.L + zh))
|
||||
|
||||
scale = 100
|
||||
oo = geompy.MakeVertex(0, 0, 0)
|
||||
@ -194,12 +118,12 @@ def simpleHexagonalPrism(theta = 0.01, fillet = False, direction = [1, 1, 1]):
|
||||
|
||||
inletface = geompy.MakeFaceWires([sk.wire()], False)
|
||||
vecflow = geompy.GetNormal(inletface)
|
||||
hexagonPrism = geompy.MakePrismVecH(inletface, vecflow, L * sqrt(3))
|
||||
poreCell = geompy.MakePrismVecH(inletface, vecflow, self.L * sqrt(3))
|
||||
|
||||
inletface = geompy.MakeScaleTransform(inletface, oo, scale)
|
||||
hexagonPrism = geompy.MakeScaleTransform(hexagonPrism, oo, scale)
|
||||
poreCell = geompy.MakeScaleTransform(poreCell, oo, scale)
|
||||
|
||||
faces = geompy.ExtractShapes(hexagonPrism, geompy.ShapeType["FACE"], False)
|
||||
faces = geompy.ExtractShapes(poreCell, geompy.ShapeType["FACE"], False)
|
||||
symetryface = []
|
||||
|
||||
for face in faces:
|
||||
@ -212,6 +136,9 @@ def simpleHexagonalPrism(theta = 0.01, fillet = False, direction = [1, 1, 1]):
|
||||
else:
|
||||
symetryface.append(face)
|
||||
|
||||
else:
|
||||
raise Exception(f"Direction { self.direction } is not implemented")
|
||||
|
||||
###
|
||||
# Grains
|
||||
##
|
||||
@ -219,9 +146,9 @@ def simpleHexagonalPrism(theta = 0.01, fillet = False, direction = [1, 1, 1]):
|
||||
oy = geompy.MakeVectorDXDYDZ(0, 1, 0)
|
||||
oz = geompy.MakeVectorDXDYDZ(0, 0, 1)
|
||||
|
||||
grain = geompy.MakeSphereR(radius)
|
||||
lattice = geompy.MakeMultiTranslation2D(grain, ox, L, xn, oy, L, yn)
|
||||
lattice = geompy.MakeMultiTranslation1D(lattice, oz, L, zn)
|
||||
grain = geompy.MakeSphereR(self.radius)
|
||||
lattice = geompy.MakeMultiTranslation2D(grain, ox, self.L, xn, oy, self.L, yn)
|
||||
lattice = geompy.MakeMultiTranslation1D(lattice, oz, self.L, zn)
|
||||
|
||||
grains = geompy.ExtractShapes(lattice, geompy.ShapeType["SOLID"], True)
|
||||
grains = geompy.MakeFuseList(grains, False, False)
|
||||
@ -229,14 +156,14 @@ def simpleHexagonalPrism(theta = 0.01, fillet = False, direction = [1, 1, 1]):
|
||||
grains = geompy.MakeScaleTransform(grains, oo, scale)
|
||||
grainsOrigin = None
|
||||
|
||||
if fillet:
|
||||
if self.filletsEnabled:
|
||||
grainsOrigin = geompy.MakeScaleTransform(grains, oo, 1 / scale)
|
||||
grains = geompy.MakeFilletAll(grains, filletradius * scale)
|
||||
grains = geompy.MakeFilletAll(grains, self.fillets * scale)
|
||||
|
||||
###
|
||||
# Groups
|
||||
##
|
||||
shape = geompy.MakeCutList(hexagonPrism, [grains])
|
||||
shape = geompy.MakeCutList(poreCell, [grains])
|
||||
shape = geompy.MakeScaleTransform(shape, oo, 1 / scale, theName = "simpleCubic")
|
||||
|
||||
sall = geompy.CreateGroup(shape, geompy.ShapeType["FACE"])
|
||||
@ -269,7 +196,7 @@ def simpleHexagonalPrism(theta = 0.01, fillet = False, direction = [1, 1, 1]):
|
||||
groups.append(outlet)
|
||||
groups.extend(symetry)
|
||||
|
||||
if fillet:
|
||||
if self.filletsEnabled:
|
||||
strips = geompy.CreateGroup(shape, geompy.ShapeType["FACE"], theName = "strips")
|
||||
shapeShell = geompy.ExtractShapes(shape, geompy.ShapeType["SHELL"], True)
|
||||
stripsShape = geompy.MakeCutList(shapeShell[0], groups + [grainsOrigin])
|
||||
@ -282,16 +209,3 @@ def simpleHexagonalPrism(theta = 0.01, fillet = False, direction = [1, 1, 1]):
|
||||
|
||||
return shape, groups
|
||||
|
||||
|
||||
def simple(theta, fillet, direction):
|
||||
#from salomepl import getGeom
|
||||
|
||||
if direction in [[1, 0, 0], [0, 0, 1]]:
|
||||
return simpleCubic(theta, fillet, direction)
|
||||
|
||||
elif direction == [1, 1, 1]:
|
||||
return simpleHexagonalPrism(theta, fillet, direction)
|
||||
|
||||
else:
|
||||
raise Exception("This direction is not implemented")
|
||||
|
||||
|
||||
Loading…
Reference in New Issue
Block a user