[教育]有限元 第三章.ppt

3.5 阶谱单元（简介） 3.5.1 一维阶谱单元 1. 标准式的C0型单元 —— Lagrange单元 ——未知场函数 （1） ?i ——各结点的参数都具有相同的物理意义 （结点的函数值）； （2） （3） —— 构造插值函数 Ni 的依据（方法） 缺点： （1） 由低阶单元升为高阶单元时，原低阶单元的各插值函数随之变化。如 x l 1 2 一次单元： 二次单元： 3 （2） 不便于编制通用计算程序（需重新构造Ni 、Ke、Pe等） （3） 不适合于自适应有限元分析。 2. 阶谱单元的概念 功用： 当发现用较低阶单元对某一问题分析的精度不够时，可在单元网格不变的条件下，提高单元的阶次（增加结点）以提高单元的精度（ANSYS — P方法），同时低阶单元的刚度矩阵又能被利用，提高计算效率。 缺点： （1） 由低阶单元升为高阶单元时，原低阶单元的各插值函数随之变化。如 x l 1 2 一次单元： 二次单元： 3 * 第 3 章 单元和插值函数的构造 要点： —— 插值函数 Ni 的构造 插值函数 Ni 的构造的基本理论与方法 重点： —— C0类单元的插值函数 （Ni ）构造 主 要 内 容 3.1 引 言 3.2 一维单元 3.3 二维单元 3.4 三维单元 3.6 小 结 3.5 阶谱单元 3.1 引 言 一、有限元法的求解步骤 （1）划分单元； （2）构造单元位移，求插值函数Ni （3）求应变矩阵B、应力矩阵S； （4）求 K e、P e； （5）集成 K 、P ； （6）引入位移边界条件 ； 修正：K 、P （7）求解方程 及其它计算。 二、影响插值函数 Ni 的因素 （1）单元形状； 直线 （2）单元的结点数目； 二、影响插值函数 Ni 的因素 （1）单元形状； （3）单元的结点数目； （2）结点参数的类型： （a）未知场函数的结点值； （b）未知场函数导数的结点值； 如：梁单元、板单元等。 高次单元的优点 （a）提高精度； （b）模拟一些曲边单元，更好地适应边界形状。 三、常用插值函数 Ni 的形式 （1）Lagrange插值函数； （2）Hermite 插值函数； 适用于只需对未知场函数的结点值插值的情形。 适用于除未知场函数的结点值外，还需对未知函数导数值插值的情形。 —— C0类插值 —— C1类插值 3.2 一维单元 3.2.1 Lagrange单元 1、插值函数 Ni 的一般性质 x ? 1 n 2 3 i-1 i i+1 l —— 插值函数 Ni 的一般形式 性质： （1） （2） 2、插值函数的Lagrange构造方法 （1）一般坐标形式： （3.2.1） （3.2.2） （3.2.3） x ? 1 n 2 3 i-1 i i+1 l 2、插值函数的Lagrange构造方法 （1）一般坐标形式： 多项式次数 显然满足： （2） （1） 如： x =0 =l x1 x2 l 1 2 （2）自然坐标形式： x =0 =l x1 x2 l 1 2 令： 则 变为 （3.2.5） （3.2.6） 如： 则： x 1 2 3 x1 x2 x2 l 有 （3.2.7） x 1 2 3 x1 x2 x2 l 则有 如： （3.2.8） （3.2.9） 自然坐标的另一种形式： x l x1 x2 1 2 xc 如： （3.2.10） x 1 2 3 x1 x2 x2 l 如： （3.2.11） （3）Lagrange单元插值函数 Ni 的一般形式： （3.2.12） 其中： —— 分子因子 —— 分母因子 说明： 显然，有 （1） —— ? 的 n 次完全多项式 也是 x 的 n 次完全多项式 —— 满足完备性准则。 （2） 说明： （1） —— ? 的 n 次完全多项式 也是 x 的 n 次完全多项式 —— 满足完备性准则。 —— 当通过坐标变换将直线单元变为曲线单元时，单元的场函数仍满足收敛性准则。 3.2.2 Hermite 单元 适用于结点参数除含未知场函数的结点值外，还含未知函数导数的结点值。 如： 梁单元、板壳单元 1. 2 结点 Hermit 插值函数 （3.2.13） 或： （3.2.14） 1. 2 结点 Hermit 插值函数 （3.2.13） 或： （3.2.14） 其中： 一阶Hermit多项式 一阶Hermit多项式的性质： （3.2.15） 一阶Hermit多项式的性质： （3.2.15） ——是构造一阶Hermit函数的准则 对于2结点情形： 取： —— 三次多项式 1 2 1 1 1 1 （3.2.16） （3.2.17） 2阶Hermit多项式（2结点）： 或： （3.2