有限元实例分析大作业二.docxVIP

周边固支受均布载荷作用的圆板： 问题描述： 几何参数为R=40cm,,板厚h=1cm的一周边固支圆板，材料参数为E=2.0x10，v=0.3，=17.32kg/cm，均布载荷强度q=-10.0kg/cm,画出圆板的变形图并计算圆心处的挠度大小。 求解分析： 对于该问题采用了patran-nastran建模求解，首先建立模型，利用patran中的P-shape功能创建平面，然后利用meshseed撒种生成mesh网格，网格采用quad4格式，周向圆盘边界种子数80，径向半径范围内种子数40，由于网格划分过程中产生重节点，所以划分完成后要利用patran的Equvalece功能完成重节点的删除，具体的网格划分情况如下图： 由于该圆板的周边属于固支，所以边界条件就是将圆盘周边的各点的沿x、y、z方向的位移都限制为0，绕x轴、y轴、z轴的转动也分别约束为0。 对于外载荷的施加，选择pressure，施加区域为圆盘的上表面，方向向下，如此就形成了均布载荷的作用。 对于材料的定义，由于该变形过程中涉及弹性和塑性变形，所以本解法在材料的属性中分别定义了其弹性参数E、v和弹塑性参数，硬化系数默认为0。 对于模型的属性选择2D /Shell，选择上步定义的材料并定义板的厚度为h=1cm。 最后是关于模型的求解参数的定义，由于该模型属于大变形，所以涉及了有限元的集合非线性，而材料方面又存在塑性变形，因此还涉及到材料的非线性，因此求解类型需要选择非线性求解，在Solution Type选择Nonlinear，由于该问题涉及到了材料非线性及几何非线性，较为复杂，为了保证分析的顺利进行及结果的精度，选择了弧长法进行求解，弧长法对步长的要求较为严格，为了保证收敛及计算效率，此处设置载荷迭代步数为50，弧长法的参数设置如下图： 设置完成后提交运算即可。 结果分析： 计算完成后，查看分析结果，在Patran/Result/Quick plot功能下绘制位移云图及变形情况如下： 可见圆盘变形呈向下凸起的形式，圆心处挠度最大，沿半径方向向周边走挠度逐渐减小，最后在固支处变为0，径向各点的挠度变化情况如下图： 随着载荷的增加圆心处的挠度的变化曲线如下： 初始载荷为总始载荷的2%，最后结束时为总载荷的100%，在此时圆心处的挠度最大，为1.20cm，该问题的理论解为1.229cm，相对误差为2.36%，计算结果与理论结果符合较好，误差在允许范围内。