薄板弯曲问题的有限元分析报告.pdfVIP

变分原理与有限元素法课程报告 报告名称：薄板弯曲问题的有限元分析 姓名： 学号： 导师： 专业： 2015.5.15 1 目录 1. 问题描述3 2. 理论基础3 2.1 矩形薄板弯曲单元 3 2.1.1 挠度函数 3 2.1.2 单元刚度矩阵 5 2.2 四边简支矩形板的纳维叶解法 5 3. 有限元模型6 4. 结果与分析7 4.1 均布载荷作用下四边简支板 7 4.2 集中载荷作用下四边简支板 8 4.2 均布载荷作用下四边固支板 9 4.2 集中载荷作用下四边固支板 10 4.5 总结11 2 1.问题描述 1 1 一块方板，边长为L，厚度为t （ L/ t ），材料为铝，分别用不同密 80 5 度的四节点12个自由度的矩形单元来划分网格。 要求：考虑四边简支和四边固支两种边界情况，分别计算受均匀载荷q和在 板中心处受集中载荷P两种载荷情况下，板的中心挠度 (不超过板厚t 的 max 1/5)，进而计算出不同情况下的方板的中心挠度系数；将计算出的系数与精确解 进行比较，通过比较发现不同有限元网格密度对薄板弯曲问题计算结果的影响。 本例中，方板边长L=40mm，厚度t=1mm，铝的弹性模量E=70GPa，泊松比  0.3 ，粗略计算当q=0.1MPa或者P=50N时，板中心挠度小于板厚的1/5，属 于小挠度弯曲，因此载荷可取这两个值。 2.理论基础 2.1矩形薄板弯曲单元 2.1.1 挠度函数 薄板弯曲单元中比较简单的是四节点12个自由度的矩形单元，将矩形薄板 沿坐标方向划分为若干矩形单元，如图1所示，每个单元设有四个节点，每个节 点位移有三个分量：挠度w，绕x轴的转角 w/ yx ，绕y轴的转角 w/ xy ， 即   wi  wi      { }    ( w / y) ( i  1,2 ,3 ,4 )     i xi i     yi  ( w / x) i    图1 3 单元的节点位移为 e T T T T T {} [{ } { } { } {