有限元法基础课件作者赵维涛陈孝珍5空间问题与板壳单元.ppt

矩形薄板单元的位移模式可取为 将薄板单元的坐标和节点位移代入上式，求解得到 ，并代回上式，则上式表示为 式中： 式中： 5.3.2 单元位移模式 式中： （i、j、m、n轮换） 式中： 5.3.2 单元位移模式 将位移代入应变表达式，可得到单元应变 式中： （i、j、m、n轮换） 5.3.3 单元应变与应力 化简后得到单元应变矩阵，如下 则单元应力可由物理方程求得，如下 5.3.3 单元应变与应力 单元刚度矩阵用分块形式表示为 将应变矩阵 、弹性矩阵 代入单元刚度一般表达式 中，得单元刚度矩阵为 式中子矩阵为3×3的矩阵，如下 5.3.4 单元刚度矩阵 式中： 5.3.4 单元刚度矩阵 式中： 至于单元节点等效载荷的求法以及总体刚阵的组装均与平面问题是类似的，这里不再给出。把四节点矩形平面应力单元和该节讨论的四节点矩形弯曲薄板单元结合起来，可以构造同时考虑中面的伸展和板弯曲的四节点矩形板单元，进而分析壳体单元。 5.3.4 单元刚度矩阵 如图（a）表示一个边长为1m的正方形薄板，四边固支，中心受垂直板面的集中载荷P作用，计算薄板中心挠度。边长l=1m，厚度h=0.01m，弹性模量E=210GPa，泊松比 ，载荷 kN。 （a） （b） 5.4 算 例 解：由于是对称结构，取1/4板进行分析，共划分4个单元，9个节点，网格划分如图（b）所示。由于对称结构，则作用于节点1上的载荷为P/4=0.25kN，其中 轴上的节点有 ，在 轴上的节点有 ，计算结果见表6.1。 表6.1计算结果 划分单元数（1/4板） 2×2 4×4 6×6 8×8 10×10 解析解 板中心挠度（10-3m） 0.31094 0.29873 0.29534 0.29398 0.29329 0.29185 5.4 算 例 尚辅网 / 第五章 空间问题与板壳单元 第五章 空间问题与板壳单元 § 5.1 空间轴对称问题 § 5.2 常用的体单元 § 5.3 板壳问题 § 5.4 算例 5.1 空间轴对称问题 在用有限元法求解轴对称问题时，采用的单元一般为整圆环，如图所示。它们和子午面rz面相交的截面，可以是直边三角形或矩形，也可以是任意四边形、曲边三角形、曲边四边形等。各个单元之间以圆环形的铰相互连接，而每一个铰与子午面rz面的交点就称为节点。如图上的i、j、m等。所有单元将在子午面rz面上形成有限元网格，与在平面问题中形成的网格一样。因为在轴对称问题中采用的单元是一个整圆环，所以在计算单元的体积时要注意到这一点。下面以三角形截面单元为例，说明如何求解轴对称问题。 如图所示子午面rz面上的一个三角形单元，设单元上任一点的径向位移（沿r向位移）分量为 ，轴向位移（沿z向位移）为量为 。由于这两个位移分量仅是r和z的函数，故令 与平面问题一样，可将位移用形函数及节点位移表示为 5.1.1 位移模式 即 式中： ——二阶单位矩阵； 、 、 ——形函数矩阵，如下 （i、j、m 轮换） 式中： 5.1.1 位移模式 由弹性力学知，轴对称问题中除平面内的应变分量 外，还有环向应变 。 因此其几何方程为 将位移代入上式，得 5.1.2 单元应变 引入 则 5.1.2 单元应变 应变矩阵 中的元素不全是常量，因此单元内的应变也不是常量，这是因为轴对称问题中采用的单元是圆环，径向的位移 必引起环向应变 ，而此应变的大小与点的位置有关。另外，由于 中含有 项，使单元的应变、应力及单元刚度矩阵的计算比平面问题复杂得多。 根据弹性力学可知 ，对于轴对称问题，有 5.1.3 单元应力 根据虚功原理或用最小位能原理，可以和平面问题一样推得其单元刚度矩阵的表达式为 ，在轴对称问题中，由于单元是一圆环，上述积分式中的微分体积 可取为微分圆环的体积，即 ，故单元刚度矩阵为 与平面问题一样，单元刚度矩阵 是一个 阶的方阵，矩阵 可分成三块，故 也可分成 个子矩阵，每个子矩阵为 阶