9.2.2　绘图中的细节考虑

9.2.2　绘图中的细节考虑

用上一节的方法计算出的圆盘的应力解存在一个问题，即图9-2中的圆盘之外的区域也存在应力值。为了消除这些区域的影响，可先判定圆盘的范围，然后将圆盘之外区域的数据点设为0。用一个语句可完成上述设置：

在应力分量数据矩阵生成后，就可以利用MATLAB绘图函数进行可视化了。图9-3是用contourf显示的σx的结果。可以发现，尽管将圆盘之外的区域设置为0，但这些设置却给绘图函数带来“麻烦”。在试件边界处，由于应力突跳，正式绘图区的数据很难显示出来（图9-3）。

图9-3　集中载荷下圆盘σx应力分布图

（a）20条等值线；（b）150条等值线

利用MATLAB的nan，可以完美解决上述问题。nan有两个特点，即参与任何运算时，其结果都是nan；绘图函数不处理nan数据点，在该数据点处显示背景色。

将圆盘区域之外的点设置为nan：

将其可视化后，可以明显地看出有用的绘图区域，但由于端部应力梯度较大，等值线非常集中，效果不够美观，如图9-4 （a）所示。可以利用shading flat命令去掉黑色的等值线，使图像更加好看些，如图9-4 （b）所示。

图9-4　巴西圆盘应力结果（σx）可视化

（a）未去掉等值线；（b）去掉等值线