MATLAB 2020a有限元分析实战：圆孔应力集中问题保姆级教程（附完整代码）

MATLAB 2020a有限元分析实战：圆孔应力集中问题保姆级教程（附完整代码）

有限元分析，听起来像是大型工业软件和高端科研的专属领域，但你知道吗？其实在MATLAB这个我们熟悉的计算环境中，就藏着一套相当强大的内置工具，足以让你从零开始，亲手完成一次完整的结构力学仿真。对于很多机械、土木、航空专业的工程师和学生来说，第一次接触有限元往往被复杂的商业软件界面和晦涩的理论吓退。今天，我们就换个思路，用MATLAB自带的PDE工具箱，来啃下“应力集中”这块经典的硬骨头。

为什么选择圆孔应力集中问题？因为它太典型了。一块平板上开个孔，在承受拉伸时，孔边局部区域的应力会急剧升高，远大于远处的名义应力。这个“应力集中系数”是工程设计中的关键参数，直接关系到结构的疲劳寿命和安全性。通过MATLAB来复现和分析这个问题，不仅能让你直观理解应力集中的现象，更能让你掌握一套从几何建模到结果可视化的完整仿真流程。无论你是想验证理论解、完成课程作业，还是为更复杂的项目打基础，这篇教程都将手把手带你走完全程，并提供可直接运行的完整代码，让你避开我当初踩过的那些坑。

1. 从零搭建你的第一个有限元模型：几何与网格

万事开头难，有限元分析的第一步——几何建模，往往就卡住了不少人。商业软件里复杂的草图绘制和特征操作，在MATLAB里可以用一种更“数学化”的方式来实现。MATLAB的PDE工具箱采用构造实体几何（CSG）的方法，通过基本图形的布尔运算（并、交、差）来构建复杂形状。对于我们的带孔平板，思路非常直接：先画一个正方形，再画一个圆，然后从正方形里“挖掉”这个圆。

1.1 使用decsg函数进行几何建模

核心函数是 decsg。它需要三个输入：几何矩阵 gm、集合公式 sf 和命名空间 ns。这听起来有点抽象，我们拆开看。

首先，几何矩阵 gm 是一个列向量，用来描述基本图形的尺寸和位置。MATLAB内置了四种基本图形，用数字1到4表示：

注意，这里的矩形是用四个顶点来定义的。对于我们的1x1正方形（左下角在原点），和半径为0.1、圆心在原点的圆，定义如下：

% 定义正方形：顶点顺序为左下、右下、右上、左上
square = [3, 4, 0, 1, 1, 0, 0, 0, 1, 1]';

% 定义圆：圆心(0,0)，半径0.1
circle = [1, 0, 0, 0.1]';

% 关键一步：对齐列向量长度
% decsg要求所有描述列向量长度一致，短的需补零
circle = [circle; zeros(length(square)-length(circle), 1)];

% 组合成几何矩阵gm，每列代表一个图形
gm = [square, circle];

注意：补零操作至关重要，否则 decsg 函数会报错。这是新手常忽略的一个细节。

其次，命名空间 ns 给每个基本图形一个标签，用于后续的布尔运算。