有限元大作业-孙博伟(1).docxVIP

SHANGHAI 《有限元分析》大作业报告 学生姓名: 孙博伟 学生学号: 5120109145 专 业: 船舶与海洋工程 日 期: 2015/6/20 1加筋板横向分布载荷分析 板单元参数：长1000mm，宽500mm，厚度15mm；E=200GPa，μ=0.3 加强筋单元：长500mm，宽15mm，高20mm；E=200GPa，μ=0.3 约束条件：四边固定 外加载荷：均布载荷0.001MPa 里兹法计算 板的刚度： 断面惯性矩： 将板的挠曲面表示成如下级数形式： 取其第一项，既 力函数： 板的应变能： 梁的应变能： 系统总位能： 又根据可得： 当时，取得最大值0.002708mm  有限元板单元 图1：加筋板板单元挠度分析 由图1知板单元加筋板受横向均布载荷的最大位移为 图2：加筋板板单元应力分析 由图2知板单元加筋板受横向均布载荷的最大应力为 有限元体单元 图3：加筋板体单元挠度分析 由图3知体单元加筋板受横向均布载荷的最大位移为 图4：加筋板体单元应力分析 由图4知体单元加筋板受横向均布载荷的最大应力为 有限元板+梁单元 图5：加筋板板+梁单元挠度分析 由图5知板+梁单元加筋板受横向均布载荷的最大挠度为 图6：加筋板板+梁单元应力分析 由图6知板+梁单元加筋板受横向均布载荷的最大应力为 对比总结 位移（mm） 误差 里兹法 0.00271 \ 板单元 0.00249 8.12% 体单元 0.00257 5.17% 板+梁单元 0.00222 18.1% 表1：加筋板最大挠度分析结果对比 四种方法得到的最大挠度结果分析总结在表1中显示，发现板+梁单元的误差最大。可能的原因是有限元方法所取的网格数目不足，而且里兹法中只取了一阶近似，在这里也存在误差。总的来说，三种有限元方法的结果都有不同程度的偏差，说明不同的单元取法得到的结果会有不同误差。 2方板圆孔受面内拉伸分析 方板参数：长、宽：600mm，厚15mm，中心圆孔直径：100mm 板材料参数：E=200GPa，μ=0.3 外加载荷：两边受拉力600N 2.1 有限元三节点划分 图7：三节点网格X-component分析结果 2.2 有限元六节点分析 图8：六节点网格X-component分析结果  2.3 有限元四节点分析 图9：四节点网格X-component分析结果 2.4 有限元八节点分析 图10：八节点网格X-component分析结果 2.5 查阅解析解 查阅Heywood经验公式有 又有，由此可得 2.6 对比总结 总结图像中的数据可得到下表： 最大应力（MPa） 平均应力（MPa） 应力集中系数 三节点单元 0.150 0.0667 2.249 六节点单元 0.203 0.0667 3.043 四节点单元 0.174 0.0667 2.609 八节点单元 0.188 0.0667 2.819 表2：圆孔方板应力集中系数结果比较 其中用到的公式有： 平均应力： 应力集中系数： 得到的几个系数均与经验公式得到的结果相近，其中三节点的差别最多，说明节点数目过少会导致计算的精准度降低。 3 个人感想 虽然早就在以前的力学课程中对有限元方法有所耳闻，但是这却是第一次真正地动手操作。这也再一次让我看到了将理论运用于实践的重要性。虽然完成本次大作业的过程并不是一帆风顺，在刚开始的阶段遇到了很多棘手的问题，但是好在有老师、助教与同学们的帮助，在学习的过程中一点点探索，一点点改进。在最初阶段为了作出一个图像绞尽脑汁，到最后熟练地完成作业任务，这中间的经历实在让人收获颇多。我也认识到了作为交大学子，在学习基础知识的同时也要注重实践，多多培养自己运用各种专业软件的能力。最后也感谢各位老师助教在完成本项大作业过程中对我提供的帮助与指导！