def fillet(event=None): Box = BRepPrimAPI_MakeBox(gp_Pnt(-400, 0, 0), 200, 230, 180).Shape() fillet = BRepFilletAPI_MakeFillet(Box) # Add fillet on each edge for e in Topo(Box).edges(): fillet.Add(20, e) blendedBox = fillet.Shape() P1 = gp_Pnt(250, 150, 75) S1 = BRepPrimAPI_MakeBox(300, 200, 200).Shape() S2 = BRepPrimAPI_MakeBox(P1, 120, 180, 70).Shape() Fuse = BRepAlgoAPI_Fuse(S1, S2) FusedShape = Fuse.Shape() fill = BRepFilletAPI_MakeFillet(FusedShape) for e in Topo(FusedShape).edges(): fill.Add(e) for i in range(1, fill.NbContours() + 1): length = fill.Length(i) Rad = 0.15 * length fill.SetRadius(Rad, i, 1) blendedFusedSolids = fill.Shape() display.EraseAll() display.DisplayShape(blendedBox) display.DisplayShape(blendedFusedSolids) display.FitAll()
def clone_tooth(base_shape): clone = gp_Trsf() grouped_shape = base_shape # Find a divisor, between 1 and 8, for the number_of teeth multiplier = 1 max_multiplier = 1 for i in range(0, 8): if num_teeth % multiplier == 0: max_multiplier = i + 1 multiplier = max_multiplier for i in range(1, multiplier): clone.SetRotation(gp_OZ(), -i * tooth_angle) rotated_shape = BRepBuilderAPI_Transform(base_shape, clone, True).Shape() grouped_shape = BRepAlgoAPI_Fuse(grouped_shape, rotated_shape).Shape() # Rotate the basic tooth and fuse together aggregated_shape = grouped_shape for i in range(1, int(num_teeth / multiplier)): clone.SetRotation(gp_OZ(), -i * multiplier * tooth_angle) rotated_shape = BRepBuilderAPI_Transform(grouped_shape, clone, True).Shape() aggregated_shape = BRepAlgoAPI_Fuse(aggregated_shape, rotated_shape).Shape() cylinder = BRepPrimAPI_MakeCylinder(gp_XOY(), top_radius - roller_diameter, thickness) aggregated_shape = BRepAlgoAPI_Fuse(aggregated_shape, cylinder.Shape()).Shape() return aggregated_shape
def boolean_fuse(shapeToCutFrom, joiningShape): from OCC.BRepAlgoAPI import BRepAlgoAPI_Fuse join = BRepAlgoAPI_Fuse(shapeToCutFrom, joiningShape) join.RefineEdges() join.FuseEdges() shape = join.Shape() #join.Destroy() return shape
def pipe_3(radius_pipe, radius_sphere): pipe1 = pipe(gp_Pnt(0, 0, 0), gp_Pnt(0, 0, 1), radius_pipe) sphere1 = sphere(gp_Pnt(0, 0, 1), radius_sphere) pipe2 = pipe(gp_Pnt(0, 0, 1), gp_Pnt(0, 1, 2), radius_pipe) sphere2 = sphere(gp_Pnt(0, 1, 2), radius_sphere) pipe3 = pipe(gp_Pnt(0, 1, 2), gp_Pnt(0, 2, 2), radius_pipe) glued1 = BRepAlgoAPI_Fuse(pipe1, sphere1).Shape() glued2 = BRepAlgoAPI_Fuse(glued1, pipe2).Shape() glued3 = BRepAlgoAPI_Fuse(glued2, sphere2).Shape() glued4 = BRepAlgoAPI_Fuse(glued3, pipe3).Shape() return glued4
def fuse(event=None): display.EraseAll() box1 = BRepPrimAPI_MakeBox(2, 1, 1).Shape() box2 = BRepPrimAPI_MakeBox(2, 1, 1).Shape() box1 = translate_topods_from_vector(box1, gp_Vec(.5, .5, 0)) fuse_shp = BRepAlgoAPI_Fuse(box1, box2).Shape() display.DisplayShape(fuse_shp) display.FitAll()
def createModel(self): Pnt = getGpPt(self.secOrigin) edges = makeEdgesFromPoints(self.points) wire = makeWireFromEdges(edges) aFace = makeFaceFromWire(wire) extrudeDir = self.L * (self.wDir) # extrudeDir is a numpy array prism = makePrismFromFace(aFace, extrudeDir) my_sphere = BRepPrimAPI_MakeSphere(Pnt, 5.0).Shape() spherebody = BRepAlgoAPI_Fuse(prism, my_sphere).Shape() return prism
def createModel(self): edges = makeEdgesFromPoints(self.points) wire = makeWireFromEdges(edges) aFace = makeFaceFromWire(wire) extrudeDir = -self.T * self.shaftDir # extrudeDir is a numpy array boltHead = makePrismFromFace(aFace, extrudeDir) cylOrigin = self.origin boltCylinder = BRepPrimAPI_MakeCylinder(gp_Ax2(getGpPt(cylOrigin), getGpDir(self.shaftDir)), self.r, self.H).Shape() whole_Bolt = BRepAlgoAPI_Fuse(boltHead, boltCylinder).Shape() return whole_Bolt
def create_model(self): edges = make_edges_from_points(self.points) wire = make_wire_from_edges(edges) aFace = make_face_from_wire(wire) extrudeDir = -self.T * self.shaft_dir # extrudeDir is a numpy array boltHead = make_prism_from_face(aFace, extrudeDir) cylOrigin = self.origin boltCylinder = BRepPrimAPI_MakeCylinder( gp_Ax2(get_gp_pt(cylOrigin), get_gp_dir(self.shaft_dir)), self.r, self.H).Shape() whole_Bolt = BRepAlgoAPI_Fuse(boltHead, boltCylinder).Shape() return whole_Bolt
def fuse(self, toFuse): """ Fuse shapes together """ fuse_op = BRepAlgoAPI_Fuse(self.wrapped, toFuse.wrapped) fuse_op.RefineEdges() fuse_op.FuseEdges() # fuse_op.SetFuzzyValue(TOLERANCE) fuse_op.Build() return Shape.cast(fuse_op.Shape())
def simple_mesh(): # # Create the shape # shape = BRepPrimAPI_MakeBox(200, 200, 200).Shape() theBox = BRepPrimAPI_MakeBox(200, 60, 60).Shape() theSphere = BRepPrimAPI_MakeSphere(gp_Pnt(100, 20, 20), 80).Shape() shape = BRepAlgoAPI_Fuse(theSphere, theBox).Shape() # # Mesh the shape # BRepMesh_IncrementalMesh(shape, 0.8) builder = BRep_Builder() comp = TopoDS_Compound() builder.MakeCompound(comp) bt = BRep_Tool() ex = TopExp_Explorer(shape, TopAbs_FACE) while ex.More(): face = topods_Face(ex.Current()) location = TopLoc_Location() facing = (bt.Triangulation(face, location)).GetObject() tab = facing.Nodes() tri = facing.Triangles() for i in range(1, facing.NbTriangles()+1): trian = tri.Value(i) index1, index2, index3 = trian.Get() for j in range(1, 4): if j == 1: m = index1 n = index2 elif j == 2: n = index3 elif j == 3: m = index2 me = BRepBuilderAPI_MakeEdge(tab.Value(m), tab.Value(n)) if me.IsDone(): builder.Add(comp, me.Edge()) ex.Next() display.EraseAll() display.DisplayShape(shape) display.DisplayShape(comp, update=True)
def createModel(self): edges = makeEdgesFromPoints(self.points) wire = makeWireFromEdges(edges) aFace = makeFaceFromWire(wire) extrudeDir = -self.T * self.shaftDir # extrudeDir is a numpy array boltHead = makePrismFromFace(aFace, extrudeDir) mkFillet = BRepFilletAPI_MakeFillet(boltHead) anEdgeExplorer = TopExp_Explorer(boltHead, TopAbs_EDGE) while anEdgeExplorer.More(): aEdge = topods.Edge(anEdgeExplorer.Current()) mkFillet.Add(self.T / 17. , aEdge) anEdgeExplorer.Next() boltHead = mkFillet.Shape() cylOrigin = self.origin boltCylinder = BRepPrimAPI_MakeCylinder(gp_Ax2(getGpPt(cylOrigin), getGpDir(self.shaftDir)), self.r, self.H).Shape() whole_Bolt = BRepAlgoAPI_Fuse(boltHead,boltCylinder).Shape() mkFillet = BRepFilletAPI_MakeFillet(whole_Bolt) return whole_Bolt
def create_single_assembly(self, tire_radius, wheel_radius, struct_length, scale=1): tire_radius = tire_radius wheel_radius = wheel_radius struct_length = struct_length straight_arm_height = tire_radius * 1.5 ring_radius = wheel_radius torus_radius = tire_radius torus = BRepPrimAPI_MakeTorus(torus_radius, ring_radius).Shape() straight_arm_radius = ring_radius / 2 straight_arm = BRepPrimAPI_MakeCylinder(straight_arm_radius, straight_arm_height).Shape() angled_arm_radius = straight_arm_radius angled_arm_height = tire_radius * 1.5 angle = -40 angled_arm = BRepPrimAPI_MakeCylinder(angled_arm_radius, angled_arm_height).Shape() traf = CTiglTransformation() traf.add_rotation_x(angle) angled_arm = traf.transform(angled_arm) traf = CTiglTransformation() traf.add_translation(0, angled_arm_height * sin(angle), 0) angled_arm = traf.transform(angled_arm) traf = CTiglTransformation() traf.add_translation(0, angled_arm_height * cos(angle), angled_arm_height * sin(angle)) straight_arm = traf.transform(straight_arm) arm_assembly = BRepAlgoAPI_Fuse(straight_arm, angled_arm).Shape() traf = CTiglTransformation() traf.add_rotation_x(90) torus = traf.transform(torus) pin_radius = wheel_radius pin_height = tire_radius pin = BRepPrimAPI_MakeCylinder(pin_radius, pin_height).Shape() traf = CTiglTransformation() traf.add_translation(0, 0, -straight_arm_height + torus_radius / 2) torus = traf.transform(torus) traf = CTiglTransformation() traf.add_rotation_x(90) pin = traf.transform(pin) traf = CTiglTransformation() traf.add_translation(0, 0, -straight_arm_height + pin_radius) pin = traf.transform(pin) traf = CTiglTransformation() traf.add_mirroring_at_xzplane() mirrored_arm_assembly = traf.transform(arm_assembly) traf = CTiglTransformation() traf.add_mirroring_at_xzplane() mirrored_pin = traf.transform(pin) arm_assembly = BRepAlgoAPI_Fuse(mirrored_arm_assembly, arm_assembly).Shape() struct_radius = wheel_radius - 0.5 * wheel_radius struct_height = struct_length struct = BRepPrimAPI_MakeCylinder(struct_radius, struct_height).Shape() traf = CTiglTransformation() traf.add_translation(0, 0, torus_radius) struct = traf.transform(struct) pin = BRepAlgoAPI_Fuse(mirrored_pin, pin).Shape() part_1 = BRepAlgoAPI_Fuse(arm_assembly, pin).Shape() part_2 = BRepAlgoAPI_Fuse(part_1, struct).Shape() part_3 = BRepAlgoAPI_Fuse(part_2, pin).Shape() return BRepAlgoAPI_Fuse(part_3, torus).Shape()
def boolean_fuse(base): ring_radius = 0.25 torus_radius = 1.0 - ring_radius torus = BRepPrimAPI_MakeTorus(torus_radius, ring_radius).Shape() fuse = BRepAlgoAPI_Fuse(base, torus).Shape() return fuse
def brep_feat_extrusion_protrusion(event=None): #Extrusion S = BRepPrimAPI_MakeBox(400.,250.,300.).Shape() faces = Topo(S).faces() F = faces.next() surf1 = BRep_Tool_Surface(F) Pl1 = Handle_Geom_Plane_DownCast(surf1).GetObject() D1 = Pl1.Pln().Axis().Direction().Reversed() MW = BRepBuilderAPI_MakeWire() p1,p2 = gp_Pnt2d(200.,-100.), gp_Pnt2d(100.,-100.) aline = GCE2d_MakeLine(p1,p2).Value() MW.Add(BRepBuilderAPI_MakeEdge(aline,surf1,0.,p1.Distance(p2)).Edge()) p1,p2 = gp_Pnt2d(100.,-100.), gp_Pnt2d(100.,-200.) aline = GCE2d_MakeLine(p1,p2).Value() MW.Add(BRepBuilderAPI_MakeEdge(aline,surf1,0.,p1.Distance(p2)).Edge()) p1,p2 = gp_Pnt2d(100.,-200.), gp_Pnt2d(200.,-200.) aline = GCE2d_MakeLine(p1,p2).Value() MW.Add(BRepBuilderAPI_MakeEdge(aline,surf1,0.,p1.Distance(p2)).Edge()) p1,p2 = gp_Pnt2d(200.,-200.), gp_Pnt2d(200.,-100.) aline = GCE2d_MakeLine(p1,p2).Value() MW.Add(BRepBuilderAPI_MakeEdge(aline,surf1,0.,p1.Distance(p2)).Edge()) MKF = BRepBuilderAPI_MakeFace() MKF.Init(surf1,False,1e-6) MKF.Add(MW.Wire()) FP = MKF.Face() BRepLib_BuildCurves3d(FP) # MKP = BRepFeat_MakePrism(S,FP,F,D1,0,True) # MKP.Perform(-200) # print 'depth 200' # res1 = MKP.Shape() # display.DisplayShape(res1) # time.sleep(1) display.EraseAll() MKP = BRepFeat_MakePrism(S,FP,F,D1,0,True) MKP.PerformThruAll() print 'depth thru all' res1 = MKP.Shape() # display.DisplayShape(res1) # Protrusion faces.next() F2 = faces.next() surf2 = BRep_Tool_Surface(F2) Pl2 = Handle_Geom_Plane_DownCast(surf2).GetObject() D2 = Pl2.Pln().Axis().Direction().Reversed() MW2 = BRepBuilderAPI_MakeWire() p1, p2 = gp_Pnt2d(100.,100.), gp_Pnt2d(200.,100.) # p1, p2 = gp_Pnt2d(100.,100.), gp_Pnt2d(150.,100.) aline = GCE2d_MakeLine(p1,p2).Value() MW2.Add(BRepBuilderAPI_MakeEdge(aline,surf2,0.,p1.Distance(p2)).Edge()) p1, p2 = gp_Pnt2d(200.,100.), gp_Pnt2d(150.,200.) aline = GCE2d_MakeLine(p1,p2).Value() MW2.Add(BRepBuilderAPI_MakeEdge(aline,surf2,0.,p1.Distance(p2)).Edge()) p1, p2 = gp_Pnt2d(150.,200.), gp_Pnt2d(100.,100.) aline = GCE2d_MakeLine(p1,p2).Value() MW2.Add(BRepBuilderAPI_MakeEdge(aline,surf2,0.,p1.Distance(p2)).Edge()) MKF2 = BRepBuilderAPI_MakeFace() MKF2.Init(surf2,False,1e-6) MKF2.Add(MW2.Wire()) MKF2.Build() # display.DisplayShape(MW2.Wire()) FP = MKF2.Face() BRepLib_BuildCurves3d(FP) MKP2 = BRepFeat_MakePrism(res1,FP,F2,D2,0,True) MKP2.PerformThruAll() display.EraseAll() # display.DisplayShape(MKP2.Shape()) trf = gp_Trsf() trf.SetTranslation(gp_Vec(0,0,300)) gtrf = gp_GTrsf() gtrf.SetTrsf(trf) tr = BRepBuilderAPI_GTransform(MKP2.Shape(), gtrf, True) from OCC.BRepAlgoAPI import BRepAlgoAPI_Fuse fused = BRepAlgoAPI_Fuse(tr.Shape(), MKP2.Shape()) fused.RefineEdges() fused.Build() print 'boolean operation error status:', fused.ErrorStatus() display.DisplayShape(fused.Shape())
def union(shp1: TopoDS_Shape, shp2: TopoDS_Shape, merge=False) -> TopoDS_Shape: """ shp1 + shape2 """ intrs = BRepAlgoAPI_Fuse(shp1, shp2) return _bool_op(intrs)
def _do_operation(self,shape1,shape2): d = self.declaration args = [shape1,shape2] if d.pave_filler: args.append(d.pave_filler) return BRepAlgoAPI_Fuse(*args)
BIHSm = 0.5 # height of small CI BIRLg = 4 # radius of large CI BIHLg = float(BIRLg*BIHSm)/BIRSm # height of large CI (proportional) BIS = 0.5 # seperation between CIs BIHPack = BIHSm+BIHLg+2*BIS # pack length BIno = int(BCL//BIHPack) # number of packs fitted into core length BIrem = BCL % BIHPack # remainder length # building bushing shape Btap = BRepPrimAPI_MakeCylinder(BOR,BTL).Shape() p = gp_Ax2(gp_Pnt(0,0,BTL),gp_DZ()) # placement Bcore = BRepPrimAPI_MakeCylinder(p,BCR,BCL).Shape() p = gp_Ax2(gp_Pnt(0,0,BTL+BCL),gp_DZ()) # placement Bhead = BRepPrimAPI_MakeCylinder(p,BOR,BHL).Shape() bushing = BRepAlgoAPI_Fuse(Btap, Bcore).Shape() bushing = BRepAlgoAPI_Fuse(bushing, Bhead).Shape() # adding composite insulators to shape BIinit = BTL+float(BIrem)/2+BIS # initial height for first cone for i in range(0,BIno): Bconez = BIinit + i*BIHPack # local height for cones p = gp_Ax2(gp_Pnt(0,0,Bconez),gp_DZ()) # placement BconeLg = BRepPrimAPI_MakeCone(p,BIRLg,BCR,BIHLg).Shape() bushing = BRepAlgoAPI_Fuse(bushing, BconeLg).Shape() p = gp_Ax2(gp_Pnt(0,0,Bconez+BIHLg+BIS),gp_DZ()) # placement BconeSm = BRepPrimAPI_MakeCone(p,BIRSm,BCR,BIHSm).Shape() bushing = BRepAlgoAPI_Fuse(bushing, BconeSm).Shape() BSIN = BOR*math.sin(math.radians(BA)) ; BCOS = BOR*math.cos(math.radians(BA))
def char_to_solid(self, glyph): layer = glyph.layers['Fore'] bodies = [] for contour in layer: i = 0 total = len(contour) curve_points = [] if total <= 2: # Can't make solid out of 1 or 2 points continue wire = BRepBuilderAPI_MakeWire() for point in contour: if point.on_curve: if i > 0: # Complete old curve curve_points.append(gp_Pnt(point.x, point.y, 0)) self.add_to_wire(curve_points, wire) if i < total: # Start new curve curve_points = [gp_Pnt(point.x, point.y, 0)] if i == total-1: first = contour[0] curve_points.append(gp_Pnt(first.x, first.y, 0)) self.add_to_wire(curve_points, wire) else: curve_points.append(gp_Pnt(point.x, point.y, 0)) if i == total-1: first = contour[0] curve_points.append(gp_Pnt(first.x, first.y, 0)) self.add_to_wire(curve_points, wire) i += 1 face = BRepBuilderAPI_MakeFace(wire.Wire()) extrusion_vector = gp_Vec(0, 0, self.thickness) prism = BRepPrimAPI_MakePrism(face.Shape(), extrusion_vector) bodies.append(dict( prism = prism, isClockwise = contour.isClockwise(), )) if len(bodies) > 0: if len(bodies) == 1: return bodies[0]['prism'].Shape() elif len(bodies) > 1: final = None positive_union = None for body in bodies: if body['isClockwise'] == 1: if positive_union: positive_union = BRepAlgoAPI_Fuse( positive_union.Shape(), body['prism'].Shape()) else: positive_union = body['prism'] negative_union = None for body in bodies: if body['isClockwise'] == 0: if negative_union: negative_union = BRepAlgoAPI_Fuse( negative_union.Shape(), body['prism'].Shape()) else: negative_union = body['prism'] if positive_union and negative_union: final = BRepAlgoAPI_Cut( positive_union.Shape(), negative_union.Shape()) elif positive_union: final = positive_union elif negative_union: final = negative_union return final.Shape()
anEdgeExplorer.Next() myBody = mkFillet # Create the neck of the bottle neckLocation = gp_Pnt(0, 0, height) neckAxis = gp_DZ() neckAx2 = gp_Ax2(neckLocation, neckAxis) myNeckRadius = thickness / 4.0 myNeckHeight = height / 10.0 mkCylinder = BRepPrimAPI_MakeCylinder(neckAx2, myNeckRadius, myNeckHeight) myBody = BRepAlgoAPI_Fuse(myBody.Shape(), mkCylinder.Shape()) # Our goal is to find the highest Z face and remove it faceToRemove = None zMax = -1 # We have to work our way through all the faces to find # the highest Z face so we can remove it for the shell aFaceExplorer = TopExp_Explorer(myBody.Shape(), TopAbs_FACE) while aFaceExplorer.More(): aFace = topods.Face(aFaceExplorer.Current()) if face_is_plane(aFace): aPlane = geom_plane_from_face(aFace) # We want the highest Z face, so compare this to the previous faces