图1-11(a NDIV 为50的网格划分及约束受载图
图1-11(b NDIV 为50的位移分布图
图1-11(c NDIV 为50的应力分布图
由以上不同单元数目的位移应力分布图可以看出，大坝截面所受位移和应力的变化趋势是相同的，最大应力都发生在坝底和水的交界点附近，最小应力发生在大坝顶端；最大变形位移也是发生在坝顶。不同单元数目下计算的数据如表1-2所示。
（4）结果分析
由以上分析结果可知：
（1）随着单元数目的增加，最大位移变化不大，应力变化范围逐步增大； （2）随着单元数目的增加，即网格划分越密，分析的结果准确度将会提高；但是单元数目的增加和节点数目的增加都会造成计算量的增加和计算速度的下降的问题。
（3）对于本次计算结果，仍可能存在虚假应力，应力的准确值无法准确得出，只是网格划分越密，计算结果越精确。所以减少虚假应力影响的措施之一就是增加单元的数目，提高网格划分的密度。
五、第3问的有限元建模及计算结果
由图1-1所示的划分方案可知，需采用手动划分网格：首先创建6个节点，然后采用不同的方式连接节点创建单元，从而分别得到两种不同的网格划分方式，见下图1-12所示。
对底边的三个节点施加全约束；载荷建立方程式并创建table ；其他的处理方式与第1小题相同。