有限元法基础-5等参元与数值积分说课讲解.pptVIP

5.4 数值积分方法 有限元法基础 例：两点Gauss积分 积分点位置： i=0 i=1 5.4 数值积分方法 有限元法基础 得到 求解高阶积分点坐标和权系数，一般利用Legendre多项式来进行。 5.4 数值积分方法 有限元法基础 5.4 数值积分方法 有限元法基础 n=2 Newton-Cotes Gauss 5.4 数值积分方法 有限元法基础 二维和三维Gauss积分 对二维积分 首先令 为常数，对 积分 再对 积分，得到 5.4 数值积分方法 有限元法基础 类似地三维积分为 注：每个方向可以选取不同的积分点数。 5.4 数值积分方法 有限元法基础 3）Irons积分方案 对三维六面体积分 每个方向使用n点Newton-Cotes积分，需 n3个点，在 每个方向的精度为n-1次。 每个方向使用m点Gauss积分，需 m3个点，在每个方 向的精度为2m-1次。 Irons积分方案通过三个方向优化节点位置，提高积分精度。 5.4 数值积分方法 有限元法基础 5.4 数值积分方法 有限元法基础 5.4 数值积分方法 有限元法基础 4）Hammer积分方案 讨论对象为面积坐标和体积坐标的积分 5.4 数值积分方法 有限元法基础 5.1等参变换的概念 有限元法基础 5）面积坐标 直边三角形时： 5.1等参变换的概念 有限元法基础 6）体积坐标 5.2 等参变换的条件与收敛性 有限元法基础 等参变换的条件 等参变换中，需计算Jacobi矩阵的逆 是否存在？ 存在的条件是 这是两个坐标系间一对一变换的条件 5.2 等参变换的条件与收敛性 有限元法基础 以二维情况为例说明 1）子单元与母单元的单元节点编号顺序相反， ，顺序相同 2） 若子单元与母单元同样是凸的，即各节点处 5.2 等参变换的条件与收敛性 有限元法基础 畸变单元举例 节点1 节点2 节点3 由于 是连续函数，故在1-2边至到2-3边时 必有一点 ，不具备等参变换条件。 5.2 等参变换的条件与收敛性 有限元法基础 畸变单元举例 边1-2 退化为一个节点 在该点处 ，也不具备 等参变换条件。 实际计算单元刚度矩阵是用数值积分， 并不会出现奇异性，应用中仍可使用； 四边形退化为三角形单元的积分精度较差。 5.2 等参变换的条件与收敛性 有限元法基础 等参单元的收敛性 弹性力学问题的收敛性包括完备性和协调性： 完备性：场插值至少一阶完备，能正确反映刚体位移和常应变。 协调性：单元内部位移连续且满足几何方程，单元间的位移场是连续的。 5.2 等参变换的条件与收敛性 有限元法基础 完备性 设单元内任一点i的位移场为 代入位移插值函数 5.2 等参变换的条件与收敛性 有限元法基础 注意到等参变换 5.2 等参变换的条件与收敛性 有限元法基础 只要 Ni 满足形函数性质，完备性就得到满足， 插值函数能够反映刚体位移和常应变。 5.2 等参变换的条件与收敛性 有限元法基础 协调性 单元间边界上的位移场： 具有相同的节点和相同的节点数 插值函数相同，有连续的位移场 插值函数满足 5. 等参元与数值积分 有限元法基础 练习题： 什么是等参元满足有限元收敛准则的条件？同样条件可否适用于次参和超参单元？ 证明边界为直线的三角形和平行四边形的二维单元的Jacobi矩阵是常数矩阵。 证明面积坐标的幂函数的积分公式。 （提示：利用面积坐标之和等于1的关系消去被积函数中的一个坐标，并注意积分上下限设置。） 5. 3 弹性力学中等参单元的一般列式方法 有限元法基础 有限元方程为 单元刚度矩阵为 5. 3 弹性力学中等参单元的一般列式方法 有限元法基础