1. 课程导入与工程问题拆解

在二维有限元分析中，求解器真正“看到”的并不是我们画在CAD里的那条光滑轮廓，而是一张由单元和节点组成的离散网格，以及施加在节点或边上的边界条件。几何模型只是分析的起点，网格划分把连续体变成有限自由度的代数系统，边界施加则决定了这个系统是否有唯一解、解是否符合真实约束。可以说，网格决定计算能不能算、算得准不准，边界决定算出来的解有没有物理意义。本课程围绕这两条主线展开，帮助你把“会点软件按钮”提升为“能判断结果是否可信”。

先看一个贯穿课程的典型二维案例：一块带中心圆孔的矩形钢板，尺寸取 200 mm × 100 mm，板厚 5 mm，材料为Q235钢，弹性模量取 206 GPa、泊松比 0.3，圆孔直径 20 mm，位于板中心。左端面完全固定，右端面承受 10 MPa 的均布拉伸载荷。这个案例看似简单，却同时包含了几何简化、应力集中、网格密度分布、约束施加和载荷换算等多个关键问题，非常适合作为二维有限元全流程的载体。

从图中可以看到，同一个物理对象在有限元中要经历“几何→离散→边界”的转换。左端固定意味着该边上所有节点的水平与竖直位移被约束为零；右端 10 MPa 的均布载荷不能直接“贴”在几何线上，而要换算为作用在右端各节点上的等效节点力。对于均布拉伸，若右端被划分为 n 段等长边，每段长度 L/n、板厚 t，则每个内部节点分担的等效拉力约为 σ·t·L/n，角节点减半。这一步换算如果做错，后续应力结果会整体偏移，而网格再漂亮也无法挽回。

网格划分的核心判断在于“把单元放在哪里”。带孔板的应力集中出现在孔边，理论应力集中系数约为 3，即远场 10 MPa 下孔边局部应力可达 30 MPa 量级（示例估算，实际值受几何比例与约束方式影响）。若孔边单元尺寸与远场一样粗，峰值应力会被严重低估。合理的做法是：远场采用较大单元控制规模，孔周设置 2～3 层渐进加密的环形区域，并保证孔边单元尺寸约为孔径的 1/10 以下。同时要检查单元质量，避免出现过大的长宽比、过小的内角或畸变单元，这些都会让刚度矩阵条件变差，甚至导致求解不收敛。

边界施加则要回答三个问题：约束哪些自由度、载荷加在哪里、参照方向如何定义。以本案例为例，坐标系取水平向右为 x 正方向、竖直向上为 y 正方向，则左端固定为 ux = uy = 0；右端载荷沿 x 正方向。若只约束左端一个点，结构会产生刚体转动，求解器可能报奇异矩阵；若把右端也错误地全约束，则载荷无处传递，结果接近零。这些错误在云图上往往表现得很“平静”，却完全失真，因此边界条件必须结合力学路径逐条核对。

本课程的学习路径按“对象—离散—边界—求解—判读”组织，但不会把每一节都做成同样的流水线。后续将依次展开：二维单元类型与形函数、网格尺寸与局部加密策略、网格质量指标、约束与载荷的等效换算、对称与周期性边界、接触与连接简化、求解设置与收敛判断、结果云图与路径提取、误差评估与网格无关性验证，最后以一个完整带孔板案例串联全流程。工具环境以通用有限元软件的操作逻辑为主线，参数与命令保持通用表达，便于你在Abaqus、ANSYS、COMSOL或自编程序中迁移。

需要特别提醒的是，二维分析中的“厚度”不是可有可无的输入。平面应力与平面应变的选择直接改变本构关系：薄板类问题取平面应力，厚壁或长柱类问题取平面应变。本案例板厚 5 mm、面内尺寸 200 mm，属于典型平面应力状态。若误选平面应变，孔边应力与变形都会出现系统性偏差。类似地，载荷是“每单位面积”还是“每单位长度”，约束是“边”还是“点”，都必须在建模前用一张简明的边界条件表固定下来，避免前后不一致。

把这张表填清楚，再动手划分网格和施加边界，你就已经避开了二维有限元中最常见的一批“隐形错误”。接下来的章节会逐项深入这些要素背后的原理与操作细节，让你不仅能做出云图，更能说清楚云图为什么长这样。