def __init__(self, r=None, radius=None, phi=360, center=False):
if radius: r = radius
C = Geom2d_Circle( gp_OX2d(), r )
self.instance = BRepBuilderAPI_MakeEdge2d(C.Circ2d())
# self.instance = BRepPrimAPI_MakeSphere(r, phi*math.pi/180)
Solid.__init__(self, self.instance.Shape(), center=center)
def points_from_curve(event=None):
radius = 5.
abscissa = 3.
circle = Geom2d_Circle(gp_OX2d(), radius, True)
gac = Geom2dAdaptor_Curve(circle.GetHandle())
ua = GCPnts_UniformAbscissa(gac, abscissa)
aSequence = []
if ua.IsDone():
n = ua.NbPoints()
for count in range(1, n + 1):
p = gp_Pnt2d()
circle.D0(ua.Parameter(count), p)
aSequence.append(p)
# convert analytic to bspline
display.DisplayShape(circle, update=True)
i = 0
for p in aSequence:
i = i + 1
pstring = 'P%i : parameter %f' % (i, ua.Parameter(i))
pnt = gp_Pnt(p.X(), p.Y(), 0)
# display points
display.DisplayShape(pnt, update=True)
display.DisplayMessage(pnt, pstring)
def curves2d_from_curves(event=None):
major, minor = 12, 4
axis = gp_OX2d()
ellipse = GCE2d_MakeEllipse(axis, major, minor).Value()
trimmed_curve = Geom2d_TrimmedCurve(ellipse, -1, 2, True)
bspline = geom2dconvert_CurveToBSplineCurve(trimmed_curve.GetHandle(),
Convert_TgtThetaOver2)
display.DisplayShape(bspline, update=True)
def curves2d_from_curves():
major, minor = 12, 4
axis = gp_OX2d()
ellipse = GCE2d_MakeEllipse(axis, major, minor).Value()
trimmed_curve = Geom2d_TrimmedCurve(ellipse, -1, 2, True)
bspline = geom2dconvert_CurveToBSplineCurve(trimmed_curve.GetHandle(),
Convert_TgtThetaOver2)
display.DisplayShape(bspline, update=True)
def points_from_curve():
radius = 5.
abscissa = 3.
circle = Geom2d_Circle(gp_OX2d(), radius, True)
gac = Geom2dAdaptor_Curve(circle.GetHandle())
ua = GCPnts_UniformAbscissa(gac, abscissa)
a_sequence = []
if ua.IsDone():
n = ua.NbPoints()
for count in range(1, n + 1):
p = gp_Pnt2d()
circle.D0(ua.Parameter(count), p)
a_sequence.append(p)
# convert analytic to bspline
display.DisplayShape(circle, update=True)
i = 0
for p in a_sequence:
i += 1
pstring = 'P%i : parameter %f' % (i, ua.Parameter(i))
pnt = gp_Pnt(p.X(), p.Y(), 0)
# display points
display.DisplayShape(pnt, update=True)
display.DisplayMessage(pnt, pstring)
def cut_out(base):
outer = gp_Circ2d(gp_OX2d(), top_radius - 1.75 * roller_diameter)
inner = gp_Circ2d(gp_OX2d(), center_radius + 0.75 * roller_diameter)
geom_outer = GCE2d_MakeCircle(outer).Value()
geom_inner = GCE2d_MakeCircle(inner).Value()
Proxy(geom_inner).Reverse()
base_angle = (2. * M_PI) / mounting_hole_count
hole_angle = atan(hole_radius / mounting_radius)
correction_angle = 3 * hole_angle
left = gp_Lin2d(gp_Origin2d(), gp_DX2d())
right = gp_Lin2d(gp_Origin2d(), gp_DX2d())
left.Rotate(gp_Origin2d(), correction_angle)
right.Rotate(gp_Origin2d(), base_angle - correction_angle)
geom_left = GCE2d_MakeLine(left).Value()
geom_right = GCE2d_MakeLine(right).Value()
inter_1 = Geom2dAPI_InterCurveCurve(geom_outer, geom_left)
inter_2 = Geom2dAPI_InterCurveCurve(geom_outer, geom_right)
inter_3 = Geom2dAPI_InterCurveCurve(geom_inner, geom_right)
inter_4 = Geom2dAPI_InterCurveCurve(geom_inner, geom_left)
if inter_1.Point(1).X() > 0:
p1 = inter_1.Point(1)
else:
p1 = inter_1.Point(2)
if inter_2.Point(1).X() > 0:
p2 = inter_2.Point(1)
else:
p2 = inter_2.Point(2)
if inter_3.Point(1).X() > 0:
p3 = inter_3.Point(1)
else:
p3 = inter_3.Point(2)
if inter_4.Point(1).X() > 0:
p4 = inter_4.Point(1)
else:
p4 = inter_4.Point(2)
trimmed_outer = GCE2d_MakeArcOfCircle(outer, p1, p2).Value()
trimmed_inner = GCE2d_MakeArcOfCircle(inner, p4, p3).Value()
plane = gp_Pln(gp_Origin(), gp_DZ())
arc1 = BRepBuilderAPI_MakeEdge(geomapi_To3d(trimmed_outer, plane)).Edge()
lin1 = BRepBuilderAPI_MakeEdge(gp_Pnt(p2.X(), p2.Y(), 0),
gp_Pnt(p3.X(), p3.Y(), 0)).Edge()
arc2 = BRepBuilderAPI_MakeEdge(geomapi_To3d(trimmed_inner, plane)).Edge()
lin2 = BRepBuilderAPI_MakeEdge(gp_Pnt(p4.X(), p4.Y(), 0),
gp_Pnt(p1.X(), p1.Y(), 0)).Edge()
cutout_wire = BRepBuilderAPI_MakeWire(arc1)
cutout_wire.Add(lin1)
cutout_wire.Add(arc2)
cutout_wire.Add(lin2)
# Turn the wire into a face
cutout_face = BRepBuilderAPI_MakeFace(cutout_wire.Wire())
filleted_face = BRepFilletAPI_MakeFillet2d(cutout_face.Face())
explorer = BRepTools_WireExplorer(cutout_wire.Wire())
while explorer.More():
vertex = explorer.CurrentVertex()
filleted_face.AddFillet(vertex, roller_radius)
explorer.Next()
cutout = BRepPrimAPI_MakePrism(filleted_face.Shape(),
gp_Vec(0.0, 0.0, thickness)).Shape()
result = base
rotate = gp_Trsf()
for i in range(0, mounting_hole_count):
rotate.SetRotation(gp_OZ(), i * 2. * M_PI / mounting_hole_count)
rotated_cutout = BRepBuilderAPI_Transform(cutout, rotate, True)
result = BRepAlgoAPI_Cut(result, rotated_cutout.Shape()).Shape()
return result
