//--------------------------------------------------------------------------------------------------
bool _IsFaceOfBendSection(TopoDS_Face face, BendParameter bendParams)
{
if (bendParams == null)
{
bendParams = new BendParameter(); // Temporary usage
}
// Each face of a bend section has two circular edges
var edges = face.Edges();
var foundCircularEdges = 0;
foreach (var edge in edges)
{
var edgeAdaptor = new BRepAdaptor_Curve(edge);
if (edgeAdaptor.GetGeomType() == GeomAbs_CurveType.GeomAbs_Circle)
{
foundCircularEdges++;
if (bendParams.Edges[0] == null)
{
bendParams.Edges[0] = edge;
bendParams.Axes[0] = edgeAdaptor.Circle().Position();
bendParams.Radii[0] = edgeAdaptor.Circle().Radius();
bendParams.AngleRad = edgeAdaptor.LastParameter() - edgeAdaptor.FirstParameter();
}
else
{
if (!bendParams.Axes[0].Axis.IsCoaxial(edgeAdaptor.Circle().Axis(), 0.01f, 0.001f) ||
!bendParams.AngleRad.IsEqual(edgeAdaptor.LastParameter() - edgeAdaptor.FirstParameter(), 0.01))
{
return(false); // Circular edge with unproper parameters detected
}
if (bendParams.Edges[1] == null)
{
// Second edge of bend surface
bendParams.Edges[1] = edge;
bendParams.Axes[1] = edgeAdaptor.Circle().Position();
break;
}
// Additional edges, find free place
for (int edgeIndex = 0; edgeIndex < bendParams.Edges.Length; edgeIndex++)
{
if (bendParams.Edges[edgeIndex] == null)
{
bendParams.Edges[edgeIndex] = edge;
break;
}
if (bendParams.Edges[edgeIndex].IsSame(edge))
{
break;
}
}
if (!edgeAdaptor.Circle().Radius().IsEqual(bendParams.Radii[0], 0.000001))
{
bendParams.Radii[1] = edgeAdaptor.Circle().Radius();
}
}
}
}
// Coaxial circle edges not found, this is not a bend section
if (foundCircularEdges != 2)
{
return(false);
}
return(true);
}
/// <summary>
/// Searches for the longest edge of a shape, and find the longest parallel edge
/// </summary>
public static (TopoDS_Edge edge, Ax1?axis, TopoDS_Edge opEdge, Ax1?opAxis) FindLongestEdge(TopoDS_Shape shape)
{
TopoDS_Edge foundEdge = null;
Ax1? foundAxis = null;
TopoDS_Edge opEdge = null;
Ax1? opAxis = null;
// Find the longest edge to revolve around
double edgeLen = -1;
var edgeInfos = new List <(TopoDS_Edge edge, double len, Ax1 axis)>();
var edges = shape.Edges();
foreach (var edge in edges)
{
var brepAdaptor = new BRepAdaptor_Curve(edge);
if (brepAdaptor.GetGeomType() != GeomAbs_CurveType.GeomAbs_Line)
{
break;
}
var v1 = brepAdaptor.Value(brepAdaptor.FirstParameter());
var v2 = brepAdaptor.Value(brepAdaptor.LastParameter());
var len = v1.Distance(v2);
var axis = new Ax1(v1, new Vec(v1, v2).ToDir());
if (edge.Orientation() == TopAbs_Orientation.TopAbs_REVERSED)
{
axis.Reverse();
}
edgeInfos.Add((edge, len, axis));
if (foundEdge != null && foundEdge.IsSame(edge))
{
// Same edge with another orientation
if (foundEdge.Orientation() == TopAbs_Orientation.TopAbs_REVERSED &&
edge.Orientation() == TopAbs_Orientation.TopAbs_FORWARD)
{
// Prefer forward edges
foundEdge = edge;
}
continue;
}
if (len > edgeLen)
{
foundEdge = edge;
edgeLen = len;
foundAxis = axis;
}
}
if (foundAxis != null)
{
// Find opposite edge
edgeLen = -1;
foreach (var ei in edgeInfos)
{
if (ei.edge == foundEdge)
{
continue;
}
if (!ei.axis.IsParallel(foundAxis.Value, 0.00001))
{
continue;
}
if (ei.len > edgeLen)
{
opEdge = ei.edge;
edgeLen = ei.len;
opAxis = ei.axis;
}
}
}
return(foundEdge, foundAxis, opEdge, opAxis);
}
