def UpdateZone(self, Surfaces, display):
print(self.Surfaces)
sewing = BRepBuilderAPI_Sewing()
for i in self.Surfaces:
sewing.Add(Surfaces[i].GHShape)
sewing.Perform()
sewed_shape = sewing.SewedShape()
tds = topods()
# FixedShape = ShapeFix_Shell(tds.Shell(sewed_shape))
# FixedShape.Perform()
solid = BRepBuilderAPI_MakeSolid(tds.Shell(sewed_shape))
# Test = BRepClass3d_SolidClassifier()
# Test.Perform()
display.DisplayShape(solid.Shape())
#print('done')
Props = GProp_GProps()
brepgprop.VolumeProperties(solid.Shape(), Props)
self.Vol = Props.Mass()
print(self.Vol)
return display
def make_sewed(*faces) -> TopoDS_Shell:
bld = BRepBuilderAPI_Sewing()
for f in faces:
bld.Add(f)
bld.Perform()
result = bld.SewedShape()
return ShapeToTopology()(result)
def sew_shapes(shapes, tolerance=0.001):
sew = BRepBuilderAPI_Sewing(tolerance)
for shp in shapes:
if isinstance(shp, list):
for i in shp:
sew.Add(i)
else:
sew.Add(shp)
sew.Perform()
print("n degenerated shapes", sew.NbDegeneratedShapes())
print("n deleted faces:", sew.NbDeletedFaces())
print("n free edges", sew.NbFreeEdges())
print("n multiple edges:", sew.NbMultipleEdges())
result = ShapeToTopology()(sew.SewedShape())
return result
def merge_faces(occ_face_list, tolerance=1e-06):
sew = BRepBuilderAPI_Sewing(tolerance)
for shp in occ_face_list:
if isinstance(shp, list):
for i in shp:
sew.Add(i)
else:
sew.Add(shp)
sew.Perform()
nfreeedge = sew.NbFreeEdges()
free_edges = []
for fe_cnt in range(nfreeedge):
free_edges.append(sew.FreeEdge(fe_cnt + 1))
face_list = wire_frm_loose_edges(free_edges)
return face_list
def make_shell_frm_faces(occ_face_list, tolerance=1e-06):
#make shell
sewing = BRepBuilderAPI_Sewing()
sewing.SetTolerance(tolerance)
for f in occ_face_list:
sewing.Add(f)
sewing.Perform()
sewing_shape = fetch.shape2shapetype(sewing.SewedShape())
#topo_dict = fetch.topos_frm_compound(sewing_shape)
#shell_list = topo_dict["shell"]
shell_list = fetch.geom_explorer(sewing_shape, "shell")
return shell_list
def combine_faces(compshape, sew_tolerance):
"""
Method to combine faces in a shell by adding connectivity and continuity
:param compshape: TopoDS_Shape
:param sew_tolerance: tolerance for sewing
:return: Topo_Shell
"""
offsew = BRepOffsetAPI_FindContigousEdges(sew_tolerance)
sew = BRepBuilderAPI_Sewing(sew_tolerance)
face_explorers = TopExp_Explorer(compshape, TopAbs_FACE)
n_faces = 0
# cycle on Faces
while face_explorers.More():
tface = topods.Face(face_explorers.Current())
sew.Add(tface)
offsew.Add(tface)
n_faces += 1
face_explorers.Next()
offsew.Perform()
offsew.Dump()
sew.Perform()
shell = sew.SewedShape()
sew.Dump()
shell = topods.Shell(shell)
shell_fixer = ShapeFix_Shell()
shell_fixer.FixFaceOrientation(shell)
if shell_fixer.Perform():
print("{} shells fixed! ".format(shell_fixer.NbShells()))
else:
print "Shells not fixed! "
new_shell = shell_fixer.Shell()
if brepalgo_IsValid(new_shell):
print "Shell valid! "
else:
print "Shell failed! "
return new_shell
e4 = BRepBuilderAPI_MakeEdge(v3, v2).Edge() e5 = BRepBuilderAPI_MakeEdge(v3, v1).Edge() # create wires w0 = BRepBuilderAPI_MakeWire(e5, e3, e4).Wire() w1 = BRepBuilderAPI_MakeWire(e1, e3, e0).Wire() w2 = BRepBuilderAPI_MakeWire(e0, e5, e2).Wire() w3 = BRepBuilderAPI_MakeWire(e2, e4, e1).Wire() # then create faces f0 = BRepBuilderAPI_MakeFace(w0).Face() f1 = BRepBuilderAPI_MakeFace(w1).Face() f2 = BRepBuilderAPI_MakeFace(w2).Face() f3 = BRepBuilderAPI_MakeFace(w3).Face() # sew the faces together to create a shell sew = BRepBuilderAPI_Sewing() sew.Add(f0) sew.Add(f1) sew.Add(f2) sew.Add(f3) sew.Perform() tetrahedron_shell = sew.SewedShape() # display the result from OCC.Display.SimpleGui import init_display display, start_display, add_menu, add_function_to_menu = init_display() display.DisplayShape(tetrahedron_shell, update=True) start_display()
def get_current_loft(self):
wires = []
curves = []
propeller_number_, xz_mirror_, xy_mirror_, yz_mirror_ \
, rot_x_, hub_length, pitch_angle, root_le_pos_x_, \
root_le_pos_y_, root_le_pos_z_, section_1_length_, section_2_length_, \
section_3_length_, section_4_length_, section_5_length_, section_1_profile_, \
section_2_profile_, section_3_profile_, section_4_profile_, section_5_profile_, \
section_1_z_, section_2_z_, section_3_z_, section_4_z_, section_5_z_ \
, section_1_chord_, section_2_chord_, section_3_chord_, section_4_chord_, \
section_5_chord_, section_1_pitch_angle_, section_2_pitch_angle_, \
section_3_pitch_angle_, section_4_pitch_angle_, section_5_pitch_angle_ = read_propeller_parameters(
name=self.name)
interval = 5
chords = [
section_1_chord_, section_1_chord_, section_2_chord_,
section_3_chord_, section_4_chord_, section_5_chord_, 0.001
]
profile = [
section_1_profile_, section_1_profile_, section_2_profile_,
section_3_profile_, section_4_profile_, section_5_profile_,
section_5_profile_
]
length = [
0, section_1_length_, section_2_length_, section_3_length_,
section_4_length_, section_5_length_ / 2, section_5_length_ / 2
]
z = [
section_1_z_, section_1_z_, section_2_z_, section_3_z_,
section_4_z_, section_5_z_, section_5_z_
]
pitch = [
0, section_1_pitch_angle_, section_2_pitch_angle_,
section_3_pitch_angle_, section_4_pitch_angle_,
section_5_pitch_angle_, section_5_pitch_angle_
]
n = random.random()
lifting_surface = self.config.get_wings().create_wing(
f"{n}", len(chords), f"naca4412")
sections = [-hub_length / 2, 0.1 * hub_length, hub_length / 2]
radius = [
0.00,
rot_x_ / 2,
rot_x_ / 2,
]
x_ = []
x_.extend(np.linspace(sections[0], sections[1], num=interval))
x_.append(0.9 * hub_length)
radii = []
radii.extend(np.linspace(radius[0], radius[1], num=interval))
radii.append((rot_x_ / 2))
print(x_, radii)
n = random.random()
hub = self.config.get_fuselages().create_fuselage(
f"hub{n}", len(x_), "circularProfile")
for (x, rad, index) in zip(x_, radii, range(1, len(x_) + 1)):
section = hub.get_section(index)
sectionElement = section.get_section_element(1)
sectionElementCenter = sectionElement.get_ctigl_section_element()
sectionElementCenter.set_center(tigl3.geometry.CTiglPoint(x, 0, 0))
sectionElementCenter.set_area(getArea(rad))
n_sections = lifting_surface.get_section_count()
lifting_surface.set_root_leposition(
tigl3.geometry.CTiglPoint(-chords[0] / 2, 0, 0))
print(n_sections)
y = 0.0
x = hub_length / 2
for (z_, chord, length_, pitch_, profile_,
idx) in zip(z, chords, length, pitch, profile,
range(1, n_sections + 1)):
profile__ = "naca" + profile_
constant = 0.0
nacanumber = profile__.split("naca")[1]
if nacanumber.isdigit():
if len(nacanumber) == 4:
constant = int(nacanumber[2:]) * 0.01
s = lifting_surface.get_section(idx)
e = s.get_section_element(1)
ce = e.get_ctigl_section_element()
ce.set_width(chord)
ce.set_height(chord * constant)
center = ce.get_center()
y += length_
center.x = x
center.y = y
center.z = 0.0
ce.set_center(center)
ce.set_profile_uid(f"{profile_}")
e.set_rotation(tigl3.geometry.CTiglPoint(0, pitch_, 0))
lifting_surface.set_rotation(
tigl3.geometry.CTiglPoint(0, pitch_angle, 0))
for idx in range(1, n_sections + 1):
s = lifting_surface.get_section(idx)
e = s.get_section_element(1)
ce = e.get_ctigl_section_element()
wires.append(ce.get_wire())
for l in wires:
adapt = BRepAdaptor_CompCurve(l)
curve = Handle_BRepAdaptor_HCompCurve(
BRepAdaptor_HCompCurve(adapt))
approx = Approx_Curve3d(curve, 0.001, GeomAbs_C2, 200, 12)
if (approx.IsDone() and approx.HasResult()):
curves.append(approx.Curve())
surface1 = tigl3.surface_factories.interpolate_curves(curves)
face1 = BRepBuilderAPI_MakeFace(surface1, 1e-6).Face()
face2 = BRepBuilderAPI_MakeFace(surface1, 1e-6).Face()
sew = BRepBuilderAPI_Sewing()
sew.Add(face1)
sew.Add(face2)
sew.Perform()
shape = sew.SewedShape()
print(shape)
tds = topods()
model = BRepBuilderAPI_MakeSolid()
model.Add(tds.Shell(shape))
solid = model.Solid()
print(solid)
rot_trafo = tigl3.geometry.CTiglTransformation()
rot_trafo.add_rotation_x(90)
loft = []
loft.append(hub.get_loft().shape())
loft.append(
tigl3.geometry.CNamedShape(rot_trafo.transform(solid),
"cut").shape())
if propeller_number_ == 2:
trafo = tigl3.geometry.CTiglTransformation()
trafo.add_mirroring_at_xzplane()
loft.append(
tigl3.geometry.CNamedShape(trafo.transform(loft[1]),
"cut").shape())
elif propeller_number_ >= 3:
delta = 360 / propeller_number_
for i in range(1, propeller_number_ + 1):
loft_copy = deepcopy(loft[1])
trafo = tigl3.geometry.CTiglTransformation()
trafo.add_rotation_x(delta * i)
print(delta * i)
loft.append(
tigl3.geometry.CNamedShape(trafo.transform(loft_copy),
"cut").shape())
builder = BRep_Builder()
assembly = TopoDS_Compound()
builder.MakeCompound(assembly)
for l in loft:
builder.Add(assembly, l)
trafo = tigl3.geometry.CTiglTransformation()
trafo.add_translation(root_le_pos_x_, root_le_pos_y_, root_le_pos_z_)
assembly = trafo.transform(assembly)
return [assembly]
def curved_surface(self):
self.read_parameters()
n = random.random()
self.lifting_surface = self.wings.create_wing(f"{n}", 5,
self.airfoil_type)
chords = [
self.section_1_chord_, self.section_2_chord_,
self.section_3_chord_, self.section_4_chord_, self.section_5_chord_
]
y = [
self.section_1_y_, self.section_2_y_, self.section_3_y_,
self.section_4_y_, self.section_5_y_
]
x = [
self.section_1_x_, self.section_2_x_, self.section_3_x_,
self.section_4_x_, self.section_5_x_
]
z = [
self.section_1_z_, self.section_2_z_, self.section_3_z_,
self.section_4_z_, self.section_5_z_
]
twist = [
self.section_1_twist_angle_, self.section_2_twist_angle_,
self.section_3_twist_angle_, self.section_4_twist_angle_,
self.section_5_twist_angle_
]
profile = self.airfoil_type
constant = 0.5
nacanumber = profile.split("naca")[1]
if nacanumber.isdigit():
if len(nacanumber) == 4:
constant = int(nacanumber[2:]) * 0.01
# decrease section size towards wing tips
wires = []
curves = []
n_sections = self.lifting_surface.get_section_count()
for idx in range(1, n_sections + 1):
s = self.lifting_surface.get_section(idx)
e = s.get_section_element(1)
ce = e.get_ctigl_section_element()
ce.set_width(chords[idx - 1])
ce.set_height(chords[idx - 1] * constant)
center = ce.get_center()
center.x = x[idx - 1]
center.y = y[idx - 1]
center.z = z[idx - 1]
ce.set_center(center)
self.lifting_surface.set_rotation(
tigl3.geometry.CTiglPoint(self.rot_x_, self.rot_y_, self.rot_z_))
self.lifting_surface.set_root_leposition(
tigl3.geometry.CTiglPoint(self.root_le_pos_x_, self.root_le_pos_y_,
self.root_le_pos_z_))
for idx in range(1, n_sections + 1):
s = self.lifting_surface.get_section(idx)
e = s.get_section_element(1)
ce = e.get_ctigl_section_element()
wires.append(ce.get_wire())
for l in wires:
adapt = BRepAdaptor_CompCurve(l)
curve = Handle_BRepAdaptor_HCompCurve(
BRepAdaptor_HCompCurve(adapt))
approx = Approx_Curve3d(curve, 0.001, GeomAbs_C2, 200, 12)
if (approx.IsDone() and approx.HasResult()):
curves.append(approx.Curve())
surface = tigl3.surface_factories.interpolate_curves(curves)
sew = BRepBuilderAPI_Sewing()
face1 = BRepBuilderAPI_MakeFace(surface, 1e-6).Face()
face2 = BRepBuilderAPI_MakeFace(surface, 1e-6).Face()
sew.Add(face1)
sew.Add(face2)
sew.Perform()
shape = sew.SewedShape()
tds = topods()
model = BRepBuilderAPI_MakeSolid()
model.Add(tds.Shell(shape))
solid = model.Solid()
loft = []
loft.append(solid)
if self.xy_mirror_:
trafo = tigl3.geometry.CTiglTransformation()
trafo.add_mirroring_at_xyplane()
loft.append(
tigl3.geometry.CNamedShape(trafo.transform(loft[0]),
"cut").shape())
elif self.xz_mirror_:
trafo = tigl3.geometry.CTiglTransformation()
trafo.add_mirroring_at_xzplane()
loft.append(
tigl3.geometry.CNamedShape(trafo.transform(loft[0]),
"cut").shape())
elif self.yz_mirror_:
trafo = tigl3.geometry.CTiglTransformation()
trafo.add_mirroring_at_yzplane()
loft.append(
tigl3.geometry.CNamedShape(trafo.transform(loft[0]),
"cut").shape())
self.old_profile = self.airfoil_type
return loft