简介：本资源是一份高质量的三维CAD课程设计源码，面向计算机、自动化等专业本科生及三维建模初学者，聚焦几何建模核心能力训练——基于半边数据结构实现欧拉操作与扫掠建模，并通过OpenGL完成实体可视化。项目完整实现5种欧拉操作（如MEV、KEMR等）及基于其构建的扫掠操作，支持带孔多边形底面沿指定向量生成三维实体，图形界面可交互旋转视角并实时渲染。压缩包共446个文件，含272个hpp头文件（主体逻辑与数据结构定义）、136个inl内联实现、20个h接口声明、4个cpp核心算法实现（EulerOperation/Sweep等），以及GLFW/Glad/OpenGL相关库与着色器文件，总大小816KB，结构清晰、模块解耦。已有281人学习下载，代码经VS2019严格调试，评审分达95分，附详细readme、输入样例（in.txt）及底面示意图，适合作为期末大作业、毕设参考或半边结构进阶实践范例。
1. 半边结构不是“画线工具”，而是三维CAD建模的底层契约
很多初学几何建模的同学一看到“半边数据结构”就下意识认为这是个“高级链表”或“带方向的边”，结果在实现欧拉操作时反复崩塌——顶点连不上、面法向翻转、扫掠后出现自交洞。其实，半边（Half-Edge）根本不是为“画图”服务的，它是对流形曲面拓扑关系的精确编码协议：每条物理边被拆成两条有向半边，每条半边明确绑定一个面、一个顶点、一个下一跳半边、一个对偶半边。这种设计让“添加一个孔”“拉伸一个面”这类操作不再依赖全局遍历或启发式猜测，而是通过固定指针跳转完成原子更新。本项目正是用 C++ 在 OpenGL 3.3 上落地这套协议：从HalfEdgeDataStructure.h中 7 个核心指针成员开始，到EulerOperation.cpp里makeFace()、killFace()等 5 个欧拉操作的指针重连逻辑，再到Sweep.cpp中将底面环沿向量平移生成侧壁并自动缝合——所有操作都严格满足欧拉公式 V − E + F = 2（对球面拓扑）。它不追求炫酷渲染，但每个旋转中的实体模型背后，都是半边指针在内存中无声而精准的握手。适合计算机图形学课程设计、CAD底层原理实践者，以及想摆脱“只会调库建模”困境的 C++ 开发者。
2. 半边数据结构的设计动机与内存布局实现
2.1 为什么必须用半边？对比面片、翼边、邻接表的失效场景
在三维 CAD 建模中，单纯存储顶点坐标（如std::vector<glm::vec3>）或三角面片（std::vector<std::array<int,3>>）无法支撑欧拉操作。例如执行MkFace（创建新面）时，若仅靠顶点索引，系统无法判断某条边属于哪个面、该边逆时针绕行时下一个顶点是谁、该边对面是否已存在——这会导致面法向混乱、孔洞无法识别。翼边结构（Winged-Edge）虽记录左右面，但未强制方向性，处理带孔多边形时需额外标记内外环；邻接表则完全丢失环序信息。半边结构通过强制方向性+双向绑定解决这些问题：每条半边HalfEdge明确指向其起点顶点origin、所属面face、下一跳半边next、对偶半边twin。这种设计使getOuterLoop()（获取外环）、isHole()（判断内环）等查询可在 O(1) 时间完成，且所有欧拉操作均只修改局部指针，不触发全局重索引。
提示：本项目中
HalfEdgeDataStructure.h的HalfEdge结构体定义为：struct HalfEdge { Vertex* origin = nullptr; Face* face = nullptr; HalfEdge* next = nullptr; HalfEdge* prev = nullptr; // 由 next 反推，但缓存提升效率 HalfEdge* twin = nullptr; Edge* edge = nullptr; // 指向共享边对象，用于属性管理 };注意
prev成员非必需，但项目中显式缓存以避免每次遍历next链反推，这是典型的空间换时间优化。
2.2 顶点、边、面三类实体的内存组织与生命周期管理
半边结构的健壮性依赖于三类实体的协同管理。本项目采用手动内存池+RAII 封装策略，而非裸指针堆分配：
Vertex类仅存储坐标 (glm::vec3 pos) 和一条关联半边 (HalfEdge* incident_edge)，后者指向以该顶点为起点的任意半边，用于快速进入环遍历；Edge类作为物理边容器，持有两条半边指针 (HalfEdge* he1,HalfEdge* he2) 和几何属性（如是否为边界边），避免半边重复计算长度；Face类存储面 ID、法向量 (glm::vec3 normal) 和一条起始半边 (HalfEdge* outer_component)，该半边必须属于外环（逆时针方向）。
关键约束在HalfEdgeDataStructure构造函数中强制执行：
// 初始化时预分配内存池，避免频繁 new/delete vertex_pool = std::make_unique<std::vector<Vertex>>(initial_capacity); edge_pool = std::make_unique<std::vector<Edge>>(initial_capacity * 2); face_pool = std::make_unique<std::vector<Face>>(initial_capacity);所有实体通过createVertex()、createEdge()等工厂方法从池中获取，析构时统一归还。这种设计杜绝了悬空指针——当killFace()删除面时，其所有半边的face指针被置为nullptr，但半边本身仍在池中待复用，后续makeFace()可直接重用内存地址。
2.3 输入解析与半边网构建：从 in.txt 到拓扑连接
main.cpp中parseInput()函数将in.txt转为半边网，流程分三步：
- 顶点批量注册：读取所有环的所有点，调用
hed.createVertex(pos)生成顶点，并存入std::vector<Vertex*> vertices； - 外环半边链构建：对每个环（按输入顺序），依次创建半边并链接
next指针，首尾相接形成闭环；同时设置origin和face（初始面 ID 递增）； - 内环与对偶关系建立：对内环（非首个环），先构建同向半边链，再调用
hed.connectHoleToOuter()—— 该函数在HalfEdgeDataStructure.h中实现，核心逻辑是：找到外环上距离内环某顶点最近的边，将其拆分为两条半边，插入内环半边作为twin，确保内环半边face指向同一面但next方向与外环相反。
此过程严格保证：每个面的外环半边face->outer_component指向逆时针环，内环半边通过twin与外环边关联，且所有twin指针双向可查。调试时可通过hed.validateTopology()检查V−E+F是否恒等于 2（对单连通物体）。
| 验证项 | 检查方式 | 失败示例 |
|---|---|---|
| 半边配对完整性 | 遍历所有半边，确认he->twin != nullptr && he->twin->twin == he | 某条半边twin为空，导致扫掠时侧壁缺失 |
| 面环方向一致性 | 对每个面，遍历outer_component链，计算多边形有向面积符号 | 内环被误判为外环，扫掠后出现面翻转 |
| 边界边识别 | 统计he->face == nullptr的半边数，应等于孔洞数×2 | 孔洞未正确连接，Sweep生成非流形几何 |
3. 五个欧拉操作的指针重连逻辑与边界条件处理
3.1 欧拉操作的数学本质：保持 V−E+F 不变量的拓扑变换
欧拉操作（Euler Operations）并非任意编辑，而是满足欧拉示性数χ = V − E + F守恒的原子操作。本项目实现的五个操作对应经典 CAD 建模原语：
MkFace：创建新面（V−E+F → (V)−(E+1)+(F+1)，χ 不变）KillFace：删除面（V−E+F → (V)−(E−1)+(F−1)，χ 不变）MkEdge：分割边（V−E+F → (V+1)−(E+1)+F，χ 不变）KillEdge：合并边（V−E+F → (V−1)−(E−1)+F，χ 不变）MkVertex：在边上插入顶点（V−E+F → (V+1)−(E+1)+F，χ 不变）
所有操作均通过修改半边指针实现，不改变顶点坐标。例如MkEdge(v1, v2, f)并非“画一条线”，而是：找到v1和v2所在面f的公共半边链，插入新半边并重连next/twin，使v1→v2成为新面边界。这种纯拓扑操作是扫掠（Sweep）能正确缝合侧壁的基础。
3.2MkFace与KillFace的实现细节：面创建与孔洞管理
MkFace是扫掠操作的前置依赖，其核心是构建新面并正确关联半边。EulerOperation.cpp中mkFace(std::vector<Vertex*> loop)实现如下：
Face* f = hed.createFace(); // 分配新面 HalfEdge* first_he = hed.createHalfEdge(); first_he->origin = loop[0]; first_he->face = f; first_he->next = nullptr; // 临时置空 HalfEdge* curr = first_he; for (size_t i = 1; i < loop.size(); ++i) { HalfEdge* next_he = hed.createHalfEdge(); next_he->origin = loop[i]; next_he->face = f; curr->next = next_he; curr = next_he; } curr->next = first_he; // 闭环 f->outer_component = first_he; // 关键：为每条新边创建对偶半边（若对面存在） for (HalfEdge* he = first_he; ; he = he->next) { Edge* e = hed.findOrCreateEdge(he->origin, he->next->origin); if (e->he1 == nullptr) { e->he1 = he; he->edge = e; } else if (e->he2 == nullptr) { e->he2 = he; he->twin = e->he1; e->he1->twin = he; } if (he->next == first_he) break; }注意：
findOrCreateEdge()通过顶点对哈希查找边，避免重复创建。twin指针在he1/he2分配后立即建立，确保后续KillFace可安全解除绑定。
KillFace则需谨慎处理孔洞：若被删面含内环，其内环半边twin必须重定向至相邻面。killFace(Face* f)中关键步骤：
// 1. 标记面内所有半边 face = nullptr for (HalfEdge* he = f->outer_component; ; he = he->next) { he->face = nullptr; if (he->next == f->outer_component) break; } // 2. 对每个内环半边，找到相邻面并重连 twin for (auto& hole_he : f->inner_components) { HalfEdge* adj_he = hole_he->twin; if (adj_he && adj_he->face) { // 相邻面存在 adj_he->face->addInnerComponent(hole_he); // 将 hole_he 归入相邻面 hole_he->twin = nullptr; // 断开旧 twin } }3.3MkEdge的边界处理：如何避免非流形几何
MkEdge(v1, v2, f)在面f内部添加对角线，但必须满足：v1和v2均在f的外环上，且不相邻（否则退化为边）。项目中通过getRingVertices(f)获取环顶点列表，再检查v1/v2索引差是否 ≥2。若满足，插入逻辑为：
- 创建新半边
he_new1（v1→v2）和he_new2（v2→v1）； - 找到
v1在环中的前驱prev_v1和后继next_v1，断开prev_v1→next_v1链，插入he_new1； - 同理处理
v2，插入he_new2； - 设置
he_new1->twin = he_new2，he_new2->twin = he_new1。
失败场景常因顶点不在同一面：MkEdge会返回nullptr并输出错误日志，而非崩溃。这种防御性编程是课程设计高分的关键——评审者看重鲁棒性，而非仅功能实现。
4. 扫掠操作的几何生成与 OpenGL 渲染管线适配
4.1 扫掠（Sweep）的拓扑-几何双阶段实现
扫掠操作sweep(HalfEdgeDataStructure& hed, const glm::vec3& direction)并非简单平移顶点，而是分两阶段：
- 拓扑阶段：基于现有底面半边网，生成侧壁半边结构。对底面每个半边
he_bottom，创建两条新半边he_side1（he_bottom->origin → he_bottom->next->origin平移后）、he_side2（反向），并链接成四边形环；同时为顶面生成新半边链，方向与底面相反（保证法向一致）。 - 几何阶段：计算所有新顶点坐标。底面顶点
v平移得v_top = v + direction，侧壁顶点由v和v_top线性插值得到（实际渲染用v和v_top构成矩形）。
Sweep.cpp中generateSideWalls()函数核心逻辑：
for (Face* f : hed.getFaces()) { if (f->isBottom()) { // 仅处理底面 HalfEdge* he = f->outer_component; do { Vertex* v1 = he->origin; Vertex* v2 = he->next->origin; // 创建侧壁四边形：v1→v2→v2_top→v1_top Vertex* v1_top = hed.createVertex(v1->pos + direction); Vertex* v2_top = hed.createVertex(v2->pos + direction); // 构建半边链：he1(v1→v2), he2(v2→v2_top), he3(v2_top→v1_top), he4(v1_top→v1) HalfEdge* he1 = hed.createHalfEdge(v1, v2, side_face); HalfEdge* he2 = hed.createHalfEdge(v2, v2_top, side_face); HalfEdge* he3 = hed.createHalfEdge(v2_top, v1_top, side_face); HalfEdge* he4 = hed.createHalfEdge(v1_top, v1, side_face); // 链接 next he1->next = he2; he2->next = he3; he3->next = he4; he4->next = he1; // 设置 twin：he1 与顶面边 twin，he2/he4 与相邻侧壁 twin he1->twin = findTopEdge(v1, v2); // 从顶面环查找 he2->twin = findAdjacentSideEdge(v2, v2_top); // 需遍历邻接面 he = he->next; } while (he != f->outer_component); } }参数说明：
direction是in.txt末尾输入的扫掠向量，单位为模型空间坐标系。项目中未做单位归一化，故输入0.0 0.0 20.0即沿 Z 轴移动 20 单位，符合 CAD 作业要求。
4.2 OpenGL 渲染管线的半边网映射：从拓扑到顶点缓冲区
Draw.h负责将半边网转换为 OpenGL 可绘制的std::vector<glm::vec3>。由于半边结构天然支持面遍历，drawSolid()函数流程为：
- 遍历所有
Face* f，跳过底面/顶面（因带孔多边形渲染未实现）； - 对每个侧面
f，提取其半边环顶点：std::vector<glm::vec3> vertices; - 每个四边形面拆为两个三角形：
(v0,v1,v2)和(v0,v2,v3)； - 将三角形顶点压入
std::vector<glm::vec3> positions，并同步填充normals（面法向）和uvs（简易纹理坐标）； - 绑定 VAO/VBO，调用
glDrawArrays(GL_TRIANGLES, 0, positions.size())。
关键优化在calculateFaceNormal()：
glm::vec3 normal(0); HalfEdge* he = f->outer_component; do { glm::vec3 e1 = he->next->origin->pos - he->origin->pos; glm::vec3 e2 = he->next->next->origin->pos - he->next->origin->pos; normal += glm::cross(e1, e2); // 累加叉积，避免单三角形误差 he = he->next; } while (he != f->outer_component); f->normal = glm::normalize(normal);4.3 交互控制与视角变换：Camera.h 的增量式更新
Camera.h实现第一人称视角，核心是processKeyboard()和processMouseScroll()。main.cpp中每帧调用：
camera.ProcessKeyboard(FORWARD, deltaTime); // W camera.ProcessKeyboard(BACKWARD, deltaTime); // S camera.ProcessKeyboard(LEFT, deltaTime); // A camera.ProcessKeyboard(RIGHT, deltaTime); // D其中ProcessKeyboard()更新cameraPos向量，再通过glm::lookAt()生成视图矩阵：
glm::mat4 getViewMatrix() { glm::vec3 front; front.x = cos(glm::radians(yaw)) * cos(glm::radians(pitch)); front.y = sin(glm::radians(pitch)); front.z = sin(glm::radians(yaw)) * cos(glm::radians(pitch)); front = glm::normalize(front); glm::vec3 right = glm::normalize(glm::cross(front, worldUp)); glm::vec3 up = glm::normalize(glm::cross(right, front)); return glm::lookAt(cameraPos, cameraPos + front, up); }注意：
yaw/pitch由鼠标移动更新，worldUp固定为(0,1,0)，确保 Y 轴为“上”。旋转实体效果由main.cpp中model = glm::rotate(model, (float)glfwGetTime(), glm::vec3(0.0f, 1.0f, 0.0f));实现，与相机解耦。
5. 编译配置、调试技巧与常见运行时问题排查
5.1 Visual Studio 2019 环境配置要点
项目使用.vcxproj文件，需确保以下依赖正确链接：
- GLFW：下载预编译二进制（
glfw-3.3.8.bin.WIN64），将include目录加入附加包含目录，lib-vc2019中glfw3.lib加入附加依赖项，DLL 放入hello_opengl.exe同目录； - GLAD：
glad.c需添加到项目源文件，glad.h路径加入包含目录，必须在glfwInit()后调用gladLoadGLLoader((GLADloadproc)glfwGetProcAddress)； - GLM：头文件库，无需编译，
#include <glm/glm.hpp>即可； - 运行库：项目属性 → C/C++ → 代码生成 → 运行库设为
/MD（动态链接），避免与 GLFW/GLAD 的 CRT 版本冲突。
提示：若出现
LNK2019: unresolved external symbol __imp__gladLoadGLLoader，检查glad.c是否在项目中编译（右键 → 属性 → 常规 → 项类型 = C++ 源文件），且glad.h路径无拼写错误。
5.2in.txt输入格式验证与调试输出
parseInput()函数内置严格校验：
- 第一行
n必须 ≥1，否则报错"Invalid number of loops"； - 每个环顶点数
m必须 ≥3，否则"Loop must have at least 3 vertices"； - 扫掠向量不能为零向量，否则
"Sweep direction vector is zero"。
调试时启用#define DEBUG_PRINT（在main.cpp顶部），程序启动后输出：
Parsed 3 loops: Loop 0 (outer): 4 vertices Loop 1 (hole): 3 vertices Loop 2 (hole): 4 vertices Sweep direction: (0.000, 0.000, 20.000) Generated 12 side faces, 0 top/bottom faces此输出可快速定位输入解析阶段问题。
5.3 图形显示异常的三层排查法
当 OpenGL 窗口黑屏或模型错乱时，按以下顺序排查：
| 层级 | 检查点 | 命令/操作 | 预期结果 |
|---|---|---|---|
| 着色器层 | shader.vs和shader.fs是否编译成功 | 在Shader::Shader()中添加glGetShaderiv(shader, GL_COMPILE_STATUS, &success) | success == GL_TRUE，否则打印glGetShaderInfoLog() |
| VAO/VBO层 | 顶点数据是否正确上传 | 在Draw.h的drawSolid()中，glBufferData()后调用glGetBufferParameteriv(GL_ARRAY_BUFFER, GL_BUFFER_SIZE, &size) | size等于positions.size() * sizeof(glm::vec3) |
| 拓扑层 | 半边网是否有效 | 在main.cpprender()前插入hed.validateTopology() | 输出"Topology valid: V=XX, E=YY, F=ZZ, χ=2" |
最常见问题是glVertexAttribPointer()的stride参数错误：本项目中顶点数据为std::vector<glm::vec3>，故stride = sizeof(glm::vec3)，若误设为0或sizeof(float)*3会导致三角形错位。
5.4 性能优化建议：从课程设计到工业级的演进路径
本项目为课程设计，但可扩展为工业级 CAD 内核：
- 内存优化：当前使用
std::vector存储实体，改为std::pmr::vector+ 自定义内存池，减少碎片； - 并行扫掠：
generateSideWalls()中每个环独立，可用std::execution::par_unseq并行化； - GPU 加速拓扑：将半边指针数组上传至 SSBO，用 Compute Shader 执行
MkFace，降低 CPU-GPU 数据拷贝； - 持久化支持：在
HalfEdgeDataStructure中添加saveToFile(const std::string& path)，序列化为.off或.obj格式。
这些优化不改变核心算法，但让代码从“能跑”迈向“可工程化”。例如，将glad.c替换为glad2并启用GLAD_GLAPI_EXPORT，即可支持 OpenGL 4.6 的glCreateBuffers()，为后续 GPU 加速铺路。
本文还有配套的精品资源，点击获取
