空间问题的有限单元法.pptVIP

所谓等参单元：即以规则形状单元（如正四边形、正六面体单元等）的位移函数相同阶次函数为单元几何边界的变换函数，进行坐标变换所获得的单元。由于单元几何边界的变换式与规则单元的位移函数有相同的节点参数，故称由此获得的单元为等参单元。借助于等参单元可以对一般任意形状的求解域方便地进行有限元离散。 进行有限元分析时，单元离散化会带来计算误差，主要采用两种方法来降低单元离散化产生的误差： 1）提供单元划分的密度，称为h方法； 2）提高单元位移函数多项式的阶次，称为p方法。 从前可知，矩形单元比三角形有更高的精度，而三角形有较矩形单元更好的边界适应性。实际工程中，往往更希望有单元精度高、边界适应性好的单元。 任意四边形单元网格划分不受边界形状限制，单元大小可以不相等。但直接对其进行单元分析及其困难，因为它的几何形状不规则，没有统一的形状，对各个单元按不同公式计算，工作量过大难以进行。 (7-18) 在自然坐标系中，各结点的形状函数可写成如下形式 对于8个顶角结点（ i＝1，2，……，8） 式中 xi、yi、zi——结点i的坐标； ui、vi、wi——结点i沿x、y、z方向的位移； Ni——对应于i结点的形状函数。 对于 的边上点（i＝9，11，13，15） 对于 的边上点（i＝10，12，14，17） 第7章 空间问题的有限单元法 对于 的边上点（i＝17，18，19，20） (7-19) 7.4.2 单元刚度矩阵 三维变形状态下，一点的应变与位移的几何关系为 (7-20) 第7章 空间问题的有限单元法 式中 为单元结点位移列向量，而单元应变转换矩阵 [B]可按结点分块表示为 其中每个子矩阵又可分为上下两块，有 （i＝1，2，…20） (7-21) 单元的刚度矩阵为 第7章 空间问题的有限单元法 (7-22) 为便于以下计算，弹性矩阵[D]可分块写为 (7-23) 令 , 则 ， 为60×60的方阵，可按结点写为子块形式 第7章 空间问题的有限单元法 式中第i行j列的子矩阵为 (7-24) 将(7-20)、(7-22)分块式代入(7-23)，其被积函数可写为 (7-25) 式中 按坐标变换式(7-17)，应有 第7章 空间问题的有限单元法 同样可有 (7-27) 三维7面体的雅可比矩阵为 (7-27) 同理可采用三维高斯求积公式计算单元刚度矩阵。即 第7章 空间问题的有限单元法 式中，L，M，N为沿 、 、 方向的积分点数目，而积分点坐标 及权重 可由高斯积分表查得。 7.5 ANSYS桁架结构计算示例 7.5.1 问题描述 如图7-4所示，一个圆柱实体。柱高0.2m，圆柱横截面直径为0.1m。约束方式：底面全约束。承受载荷： 第7章 空间问题的有限单元法 A点承受Z方向集中载荷Fz=5000N和Y方向集中载荷 Fy=-5000N；B点承受X方向集中载荷Fx=5000N；C点承受Z方向集中载荷Fz=-5000N；D点承受X方向集中载荷 Fx=-5000N。弹性模量为EX=210GP，泊松比μ=0.3。 图7-4 圆柱实体示意图 7.5.2 ANSYS求解操作过程 打开Ansys软件，在Ansys环境下做如下操作: 第7章 空间问题的有限单元法 （1）选择单元类型 运行PreprocessorElement TypeAdd/Edit/Delete，弹出Element Types对话框，如图7-5所示。然后单击 Add，弹出Library of Element Types窗口，如图7-7所 示，选择SOLID45单元，单击OK。 图7-5 单元类型对话框 图7-7 单元类型库对话框 第7章 空间问题的有限单元法 （2）设置材料属性 运行PreprocessorMaterial PropsMaterial Models，弹出如图7-7所示对话框。双击Isotropic，弹出Linear Isotropic Properties for Material Number1对话框，如图7-8所示，在EX选项