带孔盘体应力分析

盘体的几何尺寸为：边长4 米，中心的空洞半径R为0.5米。

左右面分别承受了一个均匀的10kpa的载荷

在笛卡尔坐标系下面，关于第三个维度的行为存在两种假定：

（1）应力平面条件：适用于在这个2D 平面外的第三维度的应力是可忽略的情况；

（2）应变平面条件：适用于在这个2D 平面外的第三维度的应变是可忽略的情况。应力平面条件对于那些第三维度方向非常薄的算例是合适的。应变平面条件对于第三维方向比较厚的情况下适用。

网格划分

6 个patch 如下定义：平板的每个边为一个patch，空洞为一个patch，前后面各为一个patch。左边和下面的patch 是对称面。这个对称面是一个几何的对称面。因此它需要在网格里面进行定义而不是简单的定义一个初始边界条件。用户需要在blockMeshDict 里面使用symmetryPlane 关键词来定义。frontAndBack，这个patch 代表了这是一个2D 的算例，在计算中这个patch 将被忽略。再次提醒一下，这是一个几何方面的限制，因此需要在网格内进行定义。用户需要在blockMeshDict中将其定义为empty。

边界和初始条件

应力边界条件通过traction 字符来指定。它由以下两个关键词组合来构成：

（1）在关键词traction 下指定应力边界矢量和大小。

（2）在pressure关键词下，如果这个边界面法向应力指向表面之外，就被定义为负值，因为它的up 和hole 这两个面牵引力为0，因此边界牵引力和压力设定为0。对于right 边界牵引力应该为(1e4 0 0)Pa，pressure 应该为0Pa。

物性

控制

离散格式和求解器控制

fvSchemes文件

fvSolution文件