本教程详细介绍了如何在cadquery或build123d环境中，通过访问底层open cascade technology (occt)几何数据来识别模型中的圆角边缘。文章将阐述传统方法在检测特定边缘特征时的局限性，并提供一个基于`ocp.brep.brep_tool`的实用解决方案，通过分析边缘的动态几何类型（如`geom_trimmedcurve`）来精确识别圆角，并附带示例代码及注意事项。
引言
在使用CadQuery或其现代封装build123d进行CAD模型处理时，我们经常需要从现有的STEP文件中提取几何信息，例如不同类型边的数量、面的类型、包围盒尺寸以及零件体积等。然而，当涉及到识别模型中特定特征，如圆角（fillet）或倒角（chamfer）时，CadQuery的默认API往往侧重于对象的创建而非现有对象的特征分析，这给用户带来了挑战。尽管PythonOCC提供了更底层的访问能力，但许多用户希望在CadQuery/build123d的框架内解决问题，以保持代码的一致性和简洁性。
挑战与传统方法局限
CadQuery提供了强大的选择器，例如file.edges('%PARABOLA')可以根据几何类型选择边。然而，对于圆角和倒角这类由特定操作生成的特征，它们在底层几何上可能表现为更通用的曲线类型（例如圆弧、椭圆弧或B样条曲线的片段），而非一个独立的、可直接识别的“圆角”类型。这意味着仅仅依靠edges('%ARC')或edges('%SPLINE')等选择器不足以精确识别所有圆角，并且无法区分普通圆弧与圆角操作产生的圆弧。因此，我们需要一种更深入的方法来探究边缘的几何属性。
基于底层几何类型识别的解决方案
解决此问题的关键在于利用CadQuery/build123d所基于的Open CASCADE Technology (OCCT)的底层几何工具。build123d库（作为CadQuery的现代封装）允许我们访问其包装的OCCT对象，进而使用OCP.BRep模块中的工具来分析边缘的几何特性。具体来说，我们可以通过检查边缘底层曲线的动态类型来判断其是否为圆角。
核心原理：BRep_Tool与几何类型分析
OCCT中的BRep_Tool类提供了一系列静态方法，用于从拓扑实体（如边）中提取其几何表示。其中，BRep_Tool.Curve_s(edge.wrapped, 0.0, 1.0)方法能够从给定的拓扑边（TopoDS_Edge，通过edge.wrapped获取）中提取其底层几何曲线。0.0和1.0参数代表曲线的起始和结束参数，这里用于获取整条曲线。
获取到曲线句柄后，我们可以调用其DynamicType().Name()方法来获取该几何曲线的类型名称字符串。对于由圆角操作生成的曲线，它们通常是“修剪过的曲线”（Trimmed Curve），其动态类型名称通常为Geom_TrimmedCurve。这种曲线表示一个完整曲线（如圆、椭圆或B样条）的一个片段，非常符合圆角（通常是圆弧或复杂曲线的片段）的几何特征。
示例代码与详细解析
以下代码演示了如何加载一个模型（或创建一个带有圆角的示例模型），然后遍历所有边，利用BRep_Tool识别并收集圆角边缘：
from build123d import *
from ocp_vscode import *
from OCP.BRep import BRep_Tool
# 1. 创建一个带有圆角的示例几何体
# 这里创建一个立方体并对其一个边缘进行圆角操作
b = Box(1, 1, 1)
# 对Z轴最高的一组边进行圆角，通常是顶部的四条边
b = fillet(b.edges().group_by(Axis.Z)[-1], 0.2)
# 2. 或者，如果你有一个STEP文件，可以这样导入：
# b = import_step("your_step_file.step")
# 3. 初始化一个列表用于存储检测到的圆角边缘
fillets = []
# 4. 遍历模型中的所有边缘
for edge in b.edges():
# 获取边缘的底层OCCT TopoDS_Edge对象
# 然后使用BRep_Tool.Curve_s提取其几何曲线
# 0.0和1.0表示获取整条曲线的参数范围
curve_handle = BRep_Tool.Curve_s(edge.wrapped, 0.0, 1.0)
# 获取几何曲线的动态类型名称
type_name = curve_handle.DynamicType().Name()
# 检查类型名称是否为 "Geom_TrimmedCurve"
# Geom_TrimmedCurve通常表示由圆角操作产生的曲线段
if type_name == "Geom_TrimmedCurve":
fillets.append(edge) # 如果是，则将其添加到圆角列表中
# 5. 可视化原始模型和检测到的圆角边缘
# 检测到的圆角边缘会以不同颜色显示
show(b, fillets)代码解析：
- *`from build123d import **: 导入build123d`库的所有功能，它是CadQuery的现代化封装，提供了更简洁的API。
- *`from ocp_vscode import **: 导入ocp_vscode`用于在VS Code中进行可视化，方便调试和结果展示。
- from OCP.BRep import BRep_Tool: 从OCP模块中导入BRep_Tool，这是访问OCCT底层几何工具的关键。
- b = Box(1, 1, 1) 和 b = fillet(...): 这部分代码用于创建一个带有圆角的示例模型。fillet函数对指定边缘进行圆角操作。在实际应用中，这通常会被import_step("your_step_file.step")替代。
- for edge in b.edges():: 迭代模型中的所有拓扑边缘。
- curve_handle = BRep_Tool.Curve_s(edge.wrapped, 0.0, 1.0): 这是核心步骤。edge.wrapped获取build123d边缘对象所包装的底层TopoDS_Edge对象。BRep_Tool.Curve_s则从这个拓扑边中提取出其几何表示（Geom_Curve类型）。
- type_name = curve_handle.DynamicType().Name(): 获取几何曲线的运行时类型名称。
- if type_name == "Geom_TrimmedCurve":: 判断曲线类型。Geom_TrimmedCurve是识别圆角的关键。
- fillets.append(edge): 将识别出的圆角边缘添加到列表中。
- show(b, fillets): 可视化原始模型和高亮显示的圆角边缘，以便直观验证结果。
注意事项与扩展
-
Geom_TrimmedCurve作为圆角指示器：
- Geom_TrimmedCurve是一个非常好的圆角指示器，因为它表示一个曲线的片段。圆角通常是圆弧、椭圆弧或B样条曲线的片段。
- 然而，并非所有Geom_TrimmedCurve都是圆角。例如，一个被修剪的直线段也可能是Geom_TrimmedCurve。但在CAD模型中，由圆角操作产生的曲线片段是Geom_TrimmedCurve的常见且重要用例。
- 对于更复杂的场景，可能需要结合边缘的拓扑邻接信息（例如，检查边缘所连接的两个面的夹角或曲率连续性）来进一步确认其是否为圆角。
-
倒角（Chamfers）的识别：
- 此方法主要针对圆角（fillets）。倒角（chamfers）通常是将一个尖锐的边替换为一条新的直线边。因此，倒角的边缘在几何上通常是Geom_Line，而不是Geom_TrimmedCurve。
- 要识别倒角，可能需要不同的策略，例如：
- 查找Geom_Line类型的边缘。
- 检查这些直线边缘与其相邻面的法线方向，看它们是否形成特定的角度（例如45度）。
- 分析边缘的拓扑上下文，判断其是否连接了两个通过倒角操作创建的平面。
-
性能考虑：
- 对于包含大量边缘的复杂STEP文件，遍历所有边缘并进行底层几何类型查询可能会有一定的性能开销。
- 如果需要优化，可以考虑在导入时进行一次性分析，并将结果缓存起来。
总结
通过利用build123d（或CadQuery）与OCCT底层库的结合，我们可以有效地识别CAD模型中的圆角边缘。BRep_Tool.Curve_s方法结合DynamicType().Name()提供了一种强大而灵活的方式来探究边缘的几何本质。虽然此方法主要适用于识别圆角，但理解其原理也为开发识别其他复杂几何特征的自定义分析工具奠定了基础。
