CAD 几何内核底层数学分类

以下对OCCT (OpenCASCADE)、Parasolid和ACIS三大几何内核所依赖的底层数学进行系统分类。一、基础数学1.1 线性代数内容应用场景向量运算点积、叉积法线计算、方向判定、共面/共线检测矩阵运算3×3, 4×4 齐次矩阵几何变换平移、旋转、缩放、镜像仿射变换坐标系转换、几何映射射影变换NURBS 齐次坐标表示、有理曲线/曲面的投影空间操作特征值 / 特征向量惯性张量分析、主轴方向提取奇异值分解 (SVD)伪逆 (Moore-Penrose) 计算、退化判断、鲁棒拟合1.2 解析几何点、线、面的参数方程与隐式方程距离计算点到线、点到面、线到线二次曲面 (Quadric) 方程球、柱、锥、椭球、双曲面、抛物面四次曲面环面 Torus 的隐式方程最高次项为 4 次1.3 微积分内容应用场景偏导数 / 梯度曲面法线、切平面计算曲线弧长积分曲线测量、等弧长参数化曲面面积 / 体积积分质量属性计算散度定理 (Divergence Theorem)将体积分转化为面积分用于体积和惯性矩计算雅可比矩阵参数空间到物理空间的映射、变量替换质量属性的面积分方法三大内核通用Volume13∯Sr⋅n^,dA\text{Volume} \frac{1}{3} \oiint_S \mathbf{r} \cdot \hat{\mathbf{n}} , dAVolume31​∬​S​r⋅n^,dA更高阶矩重心、惯性矩通过类似的面积分公式计算避免了直接体积分。1.4 逼近论 (Approximation Theory)内容应用场景Weierstrass 逼近定理多项式逼近的理论基础最佳一致逼近 / Chebyshev 逼近曲线/曲面的最优多项式逼近最小二乘逼近离散点集的曲线/曲面拟合插值理论 (Lagrange / Hermite)过已知点的曲线构造二、曲线与曲面几何核心2.1 参数曲线按数学表示分类: ├── 解析曲线: │ ├── 直线 (Line) │ ├── 圆 / 椭圆 (Circle / Ellipse) │ └── 双曲线 / 抛物线 (Hyperbola / Parabola) └── 自由曲线: ├── Bézier 曲线 (Bernstein 基函数) ├── B-Spline 曲线 (Cox-de Boor 递推) ├── NURBS 曲线 (有理 B-Spline) └── 插值样条 (Hermite / 三次样条)2.2 参数曲面按数学表示分类: ├── 解析曲面: │ ├── 平面 (Plane) │ ├── 二次曲面 (柱面 / 锥面 / 球面 / 椭球面 / 双曲面 / 抛物面) │ └── 环面 (Torus, 隐式四次, 但内核中作为解析曲面处理) └── 自由曲面: ├── Bézier 曲面 (张量积) ├── B-Spline 曲面 └── NURBS 曲面 按构造方式分类: ├── 旋转面 (Surface of Revolution) ├── 扫掠面 (Swept Surface) ├── 放样面 (Lofted Surface) ├── 偏置曲面 (Offset Surface) └── 裁剪曲面 (Trimmed Surface)注两种分类维度正交。例如一个旋转面可能被存储为解析曲面如球面也可能被存储为 NURBS 曲面。2.3 PCurve参数空间曲线概念说明定义三维边 (Edge) 在其所属面的参数空间 (u, v) 中的二维曲线表示数学内容二维参数空间中的曲线表示、三维曲线到曲面参数空间的投影核心作用连接拓扑与几何的关键要素裁剪曲面 (Trimmed Surface) 的数学基础每条边在其相邻的两个面上各有一条 PCurve是 BRep 数据结构中不可缺少的几何信息。2.4 NURBS 数学三大内核共同核心数学内容说明Bernstein 多项式Bézier 曲线/曲面的基函数Cox-de Boor 递推公式B-Spline 基函数的递归求值节点向量 (Knot Vector)控制基函数的支撑域与连续性齐次坐标与有理化精确表示圆锥曲线/二次曲面节点插入 (Knot Insertion)Boehm 算法、Oslo 算法节点删除 (Knot Removal)曲线/曲面简化升阶 (Degree Elevation)增加多项式次数而不改变形状de Casteljau 算法Bézier 曲线/曲面的稳定求值与细分曲线/曲面拟合最小二乘逼近、插值三、微分几何3.1 曲线微分几何概念说明切向量T曲线方向主法向量N、副法向量B构成 Frenet-Serret 标架曲率 κ弯曲程度的标量度量挠率 τ空间扭转程度的标量度量弧长参数化均匀采样、几何不变性Frenet 标架的局限性退化情况问题直线段曲率为零法向量N未定义标架退化拐点 (Inflection Point)曲率过零点法向量翻转截面突变近直线段数值不稳定替代方案方法说明校正标架 (Corrected Frenet)对直线段做插值修正平行传输标架 (Bishop Frame)无突变最为稳定固定参考方向投影将参考向量投影到法平面3.2 曲面微分几何概念说明第一基本形式 (I): E, F, G曲面上的内蕴度量长度、角度、面积第二基本形式 (II): L, M, N曲面相对于切平面的弯曲度量高斯曲率 K κ₁κ₂内蕴曲率凸/鞍/平判断平均曲率 H (κ₁κ₂)/2曲面光顺性判断主曲率 κ₁, κ₂最大/最小法曲率的数值标量主方向e₁,e₂主曲率对应的切方向向量法曲率给定切方向上的截面曲率测地线曲面上的最短路径3.3 连续性理论连续性定义应用C⁰ / C¹ / C²参数连续性依赖参数化B-Spline 节点处的数学连续性G⁰ / G¹ / G² / G³几何连续性不依赖参数化不同曲面片间的光顺拼接关键区别C¹ 连续要求参数化一致方向上的切向量相等G¹ 只要求切平面连续允许参数化不同。在 NURBS 多片拼接和圆角曲面质量控制中这个区分至关重要。3.4 曲面光顺 (Fairing / Smoothing)内容数学基础曲率梳 (Curvature Comb)曲率连续性的可视化评估能量最小化薄板样条能量、弯曲能量的泛函最小化反射线分析 (Reflection Lines)曲面质量的视觉评估方法四、拓扑学4.1 点集拓扑基础概念在 BRep 中的作用开集、闭集、边界面开区域与边界闭曲线的严格定义连通性壳体连通性判断、孔洞识别流形 / 非流形实体合法性判断可定向性 (Orientability)法线一致性排除 Möbius 带等不可定向曲面4.2 BRep 拓扑结构典型 BRep 拓扑层次: Vertex ──── Edge ──── Wire/Loop ──── Face ──── Shell ──── Solid (顶点) (边) (环) (面) (壳) (实体)三大内核的关键拓扑区别内核边-面关联机制说明OCCTWinged-Edge 变体7 层拓扑结构支持非流形ParasolidFin功能类似有向半边记录边在面上的使用信息ACISCoedge共边承载 PCurve 与方向信息注三者均非教科书标准的 Half-Edge 数据结构但都实现了类似的有向边使用机制。ACIS 的完整拓扑层次为 Body → Lump → Shell → Face → Loop →Coedge→ Edge → Vertex。4.3 欧拉-庞加莱公式简化形式单壳体、无内环V−EF2(1−g)V - E F 2(1 - g)V−EF2(1−g)完整形式含内环、多壳体V−EF−(L−F)−2(S−G)0V - E F - (L - F) - 2(S - G) 0V−EF−(L−F)−2(S−G)0即V−EF−L2S−2G0V - E F - L 2S - 2G 0V−EF−L2S−2G0其中V 顶点数E 边数F 面数L 环数包括外环和内环S 壳体数G 穿透孔数 (genus / through-holes)用于拓扑有效性验证和欧拉操作 (MEV, MEF 等) 的增量构造。4.4 图论概念应用邻接图面-面、边-面、顶点-边的连接关系表示图遍历 (BFS / DFS)拓扑遍历、连通性分析对偶图Shell 的连通性分析五、布尔运算数学5.1 曲面求交 (Surface-Surface Intersection, SSI)方法说明解析求交平面-平面、平面-二次曲面等有封闭解的情况隐式化方法将参数曲面转为隐式方程利用代数方法求解行进法 (Marching Method)从种子点出发沿交线步进追踪三大内核主流方法递归细分法 (Subdivision)递归细分包围盒逐步逼近交线5.2 曲线求交 (Curve-Curve / Curve-Surface)方法适用场景Newton-Raphson 迭代光滑曲线/曲面给定良好初值Bézier ClippingBézier 曲线间的鲁棒求交隐式化求交低次曲线间的精确求交区间方法保证找到所有根的鲁棒方法5.3 代数几何与消元理论方法应用结式 (Resultant)两多项式曲线求交的消元Gröbner 基多项式方程组的系统求解隐式化 (Implicitization)将参数曲线/曲面转化为隐式方程5.4 布尔运算流程布尔运算 (Union / Intersection / Subtraction): ├── 1. 求交 (Intersection) → SSI / CSI 计算交线 ├── 2. 边分裂 (Edge Split) → 在交点处分割边 ├── 3. 面分裂 (Face Split) → 交线将面分为子面 ├── 4. 分类 (Classification) → 判断各子面位于另一实体的内/外/上 ├── 5. 选择与拓扑重建 (Selection Rebuild) → 根据布尔类型选取保留的子面 └── 6. 结果验证与修复 (Validation Healing) → 拓扑一致性检查六、特征操作的数学基础6.1 偏置与抽壳运算数学基础曲面偏置 (Offset)沿法线方向等距推移核心难点自交检测与消除抽壳 (Shell)多曲面联合偏置需处理相邻偏置面的相交与修剪Minkowski 和 (理论基础)实体与闭球体的 Minkowski 和 外偏置膨胀内偏置需要 Minkowski 差。实际内核中较少直接计算更多作为理论分析框架6.2 倒角与圆角运算数学基础恒定半径圆角 (Constant Fillet)滚球法球沿边滚动生成的包络面变半径圆角 (Variable Fillet)变半径滚球的包络面倒角 (Chamfer)沿边的平面/规则面截切管道曲面 (Canal Surface)沿曲线移动的球族的包络面6.3 拔模 (Draft)内容说明基本原理相对于中性面 (Neutral Plane)沿拔模方向 (Pull Direction) 将面倾斜指定角度数学本质被拔模面的锥面化变换中性线面与中性面的交线保持不动其余部分按角度倾斜难点相邻面的连续性维护、自交检测6.4 扫掠与放样运算数学基础平移扫掠 (Extrude)截面沿方向向量的直线运动生成柱面旋转扫掠 (Revolve)截面沿轴旋转生成旋转面路径扫掠 (Sweep)截面沿路径曲线运动涉及截面定位标架的选择见 3.1 节 Frenet 标架及其替代方案放样 (Loft)多截面的光滑插值构造 B-Spline/NURBS 蒙皮面七、计算几何与算法7.1 空间索引与加速结构算法/结构应用k-d Tree最近点查询、范围搜索BVH (层次包围体)碰撞检测、求交加速包围盒 (AABB / OBB)快速相交预判、空间剪裁均匀网格 / 八叉树空间划分7.2 经典计算几何算法算法应用Delaunay 三角剖分曲面网格化Voronoi 图中轴变换 (Medial Axis)、距离场计算凸包算法包围体构造、干涉检测点在多边形/多面体判断布尔运算的内外分类7.3 网格化 (Tessellation / Meshing)方法说明参数空间映射网格化将 UV 空间的网格映射到三维空间Advancing Front (前沿推进法)从边界向内推进生成三角形Delaunay 网格生成满足 Delaunay 准则的三角/四面体网格自适应网格细分基于曲率/弦高误差的局部加密八、数值方法8.1 非线性方程求解方法应用Newton-Raphson 法曲线/曲面求交、最近点投影二分法 / Brent 法区间内保证收敛的根搜索不动点迭代某些特殊构造的迭代求解8.2 线性方程组求解方法应用高斯消元 / LU 分解一般线性方程组QR 分解最小二乘问题、数值稳定性Cholesky 分解对称正定系统拟合问题的法方程8.3 数值积分方法应用Gauss-Legendre 积分弧长、面积、体积计算Gauss-Kronrod 积分自适应积分与误差估计Simpson / Romberg 积分质量属性计算8.4 优化方法方法应用最小二乘法线性/非线性曲线/曲面拟合约束优化 (Lagrange 乘子法)最近点、切点、极值点计算Levenberg-Marquardt 法非线性最小二乘曲面拟合8.5 区间算术 (Interval Arithmetic)内容应用区间运算鲁棒的几何谓词判断点在线的哪侧根隔离保证找到所有交点包含函数 (Inclusion Function)Bézier/B-Spline 的值域估计8.6 容差与数值稳定性内容说明几何容差 (Geometric Tolerance)点重合、曲线/曲面相切的判定阈值拓扑容差 (Topological Tolerance)拓扑实体的合并/分离判断浮点误差分析IEEE 754 双精度下的误差传播与累积控制容差一致性几何容差与拓扑容差的协调——工程核心难题九、延伸隐式几何方法虽然三大传统内核以参数化 BRep 为核心但以下隐式方法在现代几何计算中日益重要方法数学基础特点SDF (Signed Distance Field)标量场、水平集 (Level Set)布尔运算简化为 min/max 操作F-Rep (Function Representation)实值函数的集合运算表达能力强适合增材制造十、总结数学层级架构┌──────────────────────────────────────────────────┐ │ 应用层 (Feature Operations) │ │ 布尔运算 / 圆角 / 抽壳 / 扫掠 / 放样 / 拔模 │ ├──────────────────────────────────────────────────┤ │ 算法层 (Algorithms) │ │ SSI / 网格化 / 内外分类 / 偏置 / 容差管理 │ ├──────────┬───────────────────────────────────────┤ │ 拓扑层 │ 几何层 │ │ BRep 结构 │ NURBS / 微分几何 / 解析曲面 / 参数化 │ │ 欧拉操作 │ 曲线曲面求值与微分 │ │ 方向与连通 │ 连续性理论 / 逼近论 / 光顺 │ ├──────────┴───────────────────────────────────────┤ │ 数值计算层 (Numerics) │ │ Newton法 / 数值积分 / 最小二乘 / 区间算术 / 容差 │ ├──────────────────────────────────────────────────┤ │ 基础数学层 (Foundations) │ │ 线性代数 / 微积分 / 解析几何 / 代数几何 / 图论 │ └──────────────────────────────────────────────────┘架构说明拓扑层与几何层设为并列关系。在 BRep 架构中拓扑提供结构与连接关系“什么连着什么”几何提供形状与度量信息“长什么样”二者互补共生。核心结论三大内核的数学基础高度相似——都以NURBS 几何 BRep 拓扑 布尔运算为核心。差异主要体现在算法实现的鲁棒性行进法、布尔引擎等核心算法的工程优化程度不同容差管理策略各内核采用不同的容差体系来协调几何精度与拓扑一致性拓扑遍历机制OCCTWinged-Edge 变体、ParasolidFin、ACISCoedge实现不同目标一致