南京理工大学研究生-有限元方法理论及应用考试-个人答案.docVIP

目录

TOC\o1-3\h\u1 等参单元及其应用 1

1.1 概述 1

等参单元的概念、原理 1

等参单元对有限元法工程应用的意义 1

1.2 等参单元的数值积分方法 1

等参单元刚度矩阵的数值积分方法 1

确定积分阶的原理 2

全积分单元与减缩积分单元讨论和评价 3

1.3 线性等参单元 3

全积分、减缩积分线性等参单元有关问题的分析讨论 3

1.4 等参单元的应用 5

2 分析与计算 6

2.1 四节点平面等参单元的收敛协调性 6

2.2 八节点平面等参单元 8

2.3 3节点平面三角形单元 9

2.4 20节点六面体等参单元 10

2.5 20节点六面体等参单元 11

3 上机实验 15

3.1 实验一 15

实验题目 15

实验目的 15

建模概述 15

计算结果分析与结论 16

实验体会与总结 32

3.2 实验二 33

实验题目 33

实验目的 33

建模概述 33

计算结果分析与讨论 34

实验体会与总结 36

3.3 实验三 36

实验题目 36

实验目的 36

建模概述 37

计算结果分析与结论 37

实验体会与总结 44

等参单元及其应用

概述

等参单元的概念、原理

普通单元受到两个方面的限制：（1）单元的精度。单元的节点数越多，单元精度越高；（2）单元几何上的限制。普通矩形和六面体单元都不能模拟任意形状几何体，所有几种普通单元都是直线边界，处理曲边界几何体误差较大。

为了解决上述矛盾，方法就是突破矩形单元和六面体单元几何方面的限制，使其成为任意四边形和任意六面体单元，这类单元位移模式和形函数的构造和单元列式的导出不能沿用构造简单单元的方法，必须引入等参变换，采用相同的插值函数对单元的节点坐标和节点位移在单元上进行插值。这种单元叫做等参单元。所引入的等参变换意义为在把母单元映射为实际单元的几何变换过程中采用了与位移模式相同的参数（插值函数）。

等参单元对有限元法工程应用的意义

等参单元的原理对有限元法工程应用的意义有以下几点：（1）等参单元形状、方位任意，容易构造高阶单元，适应性好，精度高；（2）等参单元列式具有统一的形式，规律性强，采用数值积分计算，程序处理方便；（3）由于等参单元涉及单元几何形状的变换，对实际单元的形态有一定要求。单元形态好坏影响计算结果的精度。单元形态应满足：A单元各方向的尺寸尽量接近；B单元边界不能过于曲折，不能有拐点和折点，尽量接近直线或抛物线；C边之间夹角接近直角。（4）高阶等参元精度高，描述复杂边界和形状的能力强，所需单元少，在结构应力分析中应用最广泛；（5）借助于等参元可以对一般的任意几何形状的工程问题和物理问题方便地进行有限元离散。因此，等参元的提出为有限元法成为现代工程实际领域最有效的数值分析方法迈出了重要的一步；（6）由于等参变换的采用使实际问题物理坐标系内的等参单元的单元刚度、质量、阻尼、载荷等特性矩阵的计算仍能在局部的自然坐标表示的规则域内进行，实现了矩阵积分式内被积函数中的导数、体积微元、面积微元、线段微元的变换及积分限的置换。因此不管各个积分式的被积函数如何复杂，都可以方便地采用标准化的数值积分方法进行计算，从而使各类不同工程实际问题的有限元分析纳入了统一的通用化程序。

等参单元的数值积分方法

等参单元刚度矩阵的数值积分方法

计算单元特性矩阵的方法有Newton-Cotes积分、高斯积分、Irons积分、Hammer积分等。为了减少积分点的数目和便于程序编制，一般采用高斯数值积分方法。

一维高斯数值积分：

结论：n阶高斯积分公式对2n-1次多项式被积函数可求得精确积分。

二位高斯数值积分：

积分公式对ξ，η方向最高方次为2n-1的多项式可求得精确值。

确定积分阶的原理

1保证刚度矩阵积分精度的积分阶选择

（1）完全积分方案

选择积分阶的基本考虑是保证被积函数所有项次精确积分，这种积分方案称为完全积分方案。

以二维单元刚度矩阵高斯数值积分为例：

单元刚度矩阵为

根据这个原理，具有规则形状的单元（|J|=常数）完全积分方案如下：

二维4节点、三维8节点等参元分别是2×2、2×2×2积分；

二维8节点、三维20节点等参元分别是3×3、3×3×3积分。

对于|J|≠常数，需要增加积分阶数。

（2）减缩积分方案

实际应用中选取的积分阶往往可以低于被积函数所有项次精确积分所需要的阶数，这种积分方案称为减缩积分。

对二、三维连续体单元，通常按形函数中完全多项式阶数p