茂之钳
73df04d5cb
M2.1 — G2/G3 连续性分析 (Agent #0): - surface_continuity.h/.cpp: G0/G1/G2/G3 curve/surface detection - Weingarten equation for curvature, Frénet frame, zebra stripe - surface_analysis.h/.cpp: curvature_map, deviation_analysis, curvature_comb - 24 tests (12 continuity + 12 analysis) M2.2 — 曲面延伸 + N边填充 (Agent #1): - surface_extension.h/.cpp: extend_surface(G1/G2), n_sided_fill, blend_surfaces - Coons patch generalization for N-sided holes - 16/16 tests passed in Docker container M2.3 — 高级过渡曲面 (Agent #2): - advanced_blend.h/.cpp: real implementations replacing stubs - variable_radius_blend, multi_face_blend, rolling_ball_blend, face_face_blend - Ball-rolling envelope + corner sphere filling - 15+ tests with validate() verification
117 lines
3.8 KiB
C++
117 lines
3.8 KiB
C++
#pragma once
|
||
#include "vde/curves/nurbs_surface.h"
|
||
#include "vde/curves/nurbs_curve.h"
|
||
#include "vde/core/point.h"
|
||
#include <vector>
|
||
#include <utility>
|
||
|
||
namespace vde::curves {
|
||
|
||
using core::Point3D;
|
||
using core::Vector3D;
|
||
|
||
/**
|
||
* @brief 连续性级别
|
||
* @ingroup curves
|
||
*/
|
||
enum class Continuity {
|
||
C0 = 0, ///< 仅位置连续
|
||
G1 = 1, ///< 切线/切平面连续
|
||
G2 = 2, ///< 曲率连续
|
||
};
|
||
|
||
/**
|
||
* @brief 曲面延伸
|
||
*
|
||
* 沿等参边界延伸曲面，支持 G1 切平面连续。
|
||
* 对于 G2 连续，在约束方向构造 3 排控制点以保持曲率。
|
||
*
|
||
* @param surface 输入 NURBS 曲面
|
||
* @param edge_id 边界边索引（0=umin, 1=umax, 2=vmin, 3=vmax）
|
||
* @param distance 延伸距离（参数空间）
|
||
* @param continuity 期望的连续性级别（默认 G1）
|
||
* @return 延伸后的 NURBS 曲面
|
||
*
|
||
* @note G1 模式：追加 1 排控制点，沿切向延伸
|
||
* @note G2 模式：追加 2 排控制点，第 2 排由曲率约束计算
|
||
* @note 延伸仅沿切向线性/二次外推，不保证与原曲面精确 G2 连续
|
||
*
|
||
* @code{.cpp}
|
||
* auto ext = extend_surface(surf, 1, 0.5, Continuity::G1);
|
||
* @endcode
|
||
* @ingroup curves
|
||
*/
|
||
[[nodiscard]] NurbsSurface extend_surface(const NurbsSurface& surface,
|
||
int edge_id, double distance,
|
||
Continuity continuity = Continuity::G1);
|
||
|
||
/**
|
||
* @brief N 边孔洞填充
|
||
*
|
||
* 用 N 条边界曲线构造 Coons patch 泛化的填充曲面。
|
||
* 支持三角形(N=3)、四边形(N=4)、五边形(N=5)等任意边数。
|
||
*
|
||
* 算法：将边界曲线映射到单位圆参数域，使用径向基函数插值
|
||
* 构造中心点，再以双线性 Coons 方式生成内部控制点。
|
||
*
|
||
* @param boundary_curves N 条边界曲线，需围成闭合环（首尾相连）
|
||
* @param continuity 期望的连续性级别（默认 G2）
|
||
* @return 填充 NURBS 曲面
|
||
*
|
||
* @note G1 模式：单层中间环，线性插值
|
||
* @note G2 模式：双层中间环，保持曲率衰减
|
||
* @note 至少需要 3 条边界曲线
|
||
*
|
||
* @code{.cpp}
|
||
* auto fill = n_sided_fill({c1, c2, c3}, Continuity::G2); // 三角形孔
|
||
* @endcode
|
||
* @ingroup curves
|
||
*/
|
||
[[nodiscard]] NurbsSurface n_sided_fill(const std::vector<NurbsCurve>& boundary_curves,
|
||
Continuity continuity = Continuity::G2);
|
||
|
||
/**
|
||
* @brief 两面间过渡曲面
|
||
*
|
||
* 检测两面公共边界，沿该边构造圆弧或 S 形过渡曲面。
|
||
* 与 nurbs_operations 中的 blend_surfaces 不同，此版本自动检测公共边。
|
||
*
|
||
* 算法：
|
||
* 1. 检测 surf_a 和 surf_b 的共享边界
|
||
* 2. 沿共享边界两侧各偏移 radius 提取曲线
|
||
* 3. 在两条偏移曲线间构造过渡曲面（直纹面或 Coons 面）
|
||
*
|
||
* @param surf_a 第一个曲面
|
||
* @param surf_b 第二个曲面
|
||
* @param radius 过渡半径
|
||
* @return 过渡 NURBS 曲面
|
||
*
|
||
* @note 如果两面无公共边界，返回退化的空曲面
|
||
* @note 过渡面在边界处与原曲面 G1 连续（等参方向）
|
||
*
|
||
* @code{.cpp}
|
||
* auto blend = blend_surfaces(surf_a, surf_b, 0.3);
|
||
* @endcode
|
||
* @ingroup curves
|
||
*/
|
||
[[nodiscard]] NurbsSurface blend_surfaces(const NurbsSurface& surf_a,
|
||
const NurbsSurface& surf_b,
|
||
double radius);
|
||
|
||
/**
|
||
* @brief 检查两曲面沿共享边界的 G1 连续性
|
||
*
|
||
* @param surf_a 第一个曲面
|
||
* @param surf_b 第二个曲面
|
||
* @param edge surf_a 的边界边索引（0=umin, 1=umax, 2=vmin, 3=vmax）
|
||
* @param samples 沿边界的采样点数（默认 10）
|
||
* @return 所有采样点均满足 G1 连续时返回 true
|
||
*
|
||
* @ingroup curves
|
||
*/
|
||
[[nodiscard]] bool is_g1_continuous(const NurbsSurface& surf_a,
|
||
const NurbsSurface& surf_b,
|
||
int edge, int samples = 10);
|
||
|
||
} // namespace vde::curves
|