07有限元Chapter7复杂单元及其实现.pptVIP

第七章 有限元分析中复杂单元及其实现 山东理工大学交通学院 李红艳 公共邮箱：yjs_fem @126.com 密码：123456 7.1 1D高阶单元 7.2 2D高阶单元 7.3 3D高阶单元 7.4 基于薄板理论的弯曲板单元 7.5 子结构和超级单元 7.6 特殊高精度单元 7.1 1D高阶单元 3节点杆单元 4节点杆单元 7.2 2D高阶单元 7.4基于薄板理论的弯曲板单元 在实际工程中，存在着大量的板壳构件，它的几何特点是其厚度远远小于其它两个方向的尺寸。因此，可以引入一定的假设来对厚度方向的受力特点进行简化，以充分体现薄板的力学特征，并且还可以减少力学变量的使用。 类似于细长梁的简化轻型，小挠度薄板理论中的简化假设就Kirchhoff假定，由于薄板中要保持转角的连续，可以承受弯矩，因此薄板问题是C1型问题。 基本变量与方程 三个基本假设（柯西霍夫假设）： 1、假设薄板中面的法线在变形后仍然为直法线，并且厚度方向的正应变很小。 2、假设薄板的中面内，没有横向位移。 3、假设z方向的正应力引起的变形很小，可以忽略。 基于以上基本假设，可以给出或者推导出三大类基本变量。 h为厚度 位移：采用中面位移为基本位移变量 应变： 应力：由前面的假设可知，剪应力为零而z方向的正应力引起的变形很小，因此可以采用平面范围内的应力分量来作为基本变量。 其中Mx，My分别为垂直于x和y轴截面上单位长度的弯矩， Mxy是垂直于y轴截面上单位长度的扭矩。 基于以上三个基本变量，可以推导出三大类基本方程。 平衡方程：外力和内力之间的平衡关系 几何方程：描述的是位移和应变之间关系 物理方程：应力和应变之间的关系 边界条件： 平衡方程 其中 薄板表面上所作用 的垂直分布载荷 平衡方程 几何方程 写成矩阵形式为 物理方程 边界条件： 1、位移和转角边界条件 2、位移和力矩边界条件 3、力矩和集中剪力边界条件 两大类非协调板单元 1、4节点矩形非协调板单元 2、3节点三角形非协调板单元 7.5 子结构和超级单元 结构的重复性表现在“几何空间” 和“计算时间” 上。 几何空间上的重复意味着结构形式在几何上是重复的；对于此可以用子结构的分析方法进行处理。 计算时将上的重复意味着结构的某一部分在多次计算中是重复的，多次计算一般出现在优化设计，非线性分析，动态分析中，需要重复进行反复的迭代。对于时间上的重复性，可以采用超级单元的方法来处理，从本质上讲，超级单元也是一种子结构。 子结构 系统中具有相同特征和性质的局部结构称为子结构。由于该结构具有几何上的多种重复性，可以将其划分为多级子结构。 取具有重复性的结构作为子结构 对最底层进行分析，形成刚度方程并进行缩聚。 对子结构进行装配形成上一级子结构（采用坐标变换并进行再缩聚） 对多级子结构全部处理完后，形成最后的整体刚度方程并求解 将结果进行回代，再求解各级子结构内部的节点位移和其它物理量 超级单元：它是一种广义的特定单元，在实际应用中可以根据需要具体产生，产生后的超级单元实际上是一个经缩聚内部节点自由度后的子结构，它表现为只有与外部有链接关系的节点位移自由度，构建超级单元的目的是减少计算量，特别是需要多次迭代的复杂计算过程可以充分体现它的优越性。由于超级单元的使用，可以大大减小每次生成刚度矩阵的计算量，同时也减小了计算规模，从而获得较高的计算效率。 一般说来，从节点看作内部节点，主节点看作主节点。相应的位移为主-从自由度。 7.6 特殊高精度单元 1、升阶谱单元 2、复合单元