有限单元法基础教学课件U4.pptxVIP

第4章 单元插值函数的构造4.1引 言 1 单元类型平面应力、平面应变问题二维实体单元实体单元轴对称问题轴对称单元三维实体问题三维实体单元一维结构单元杆单元梁单元索单元结构单元壳单元二维结构单元板单元膜单元2 单元形状二维实体单元：平面应力单元平面应力（应变）单元矩形四边形三角形三角形单元离散平面二维结构轴对称实体单元：轴对称环状单元）轴对称实体矩形四边形三角形对称面有限元离散三维实体单元：六面体四面体三棱柱单元四面体单元六面体单元一维结构单元：杆单元、梁单元、索单元空间直线空间曲线曲线单元（索单元）空间杆单元空间梁单元直线单元（杆、梁单元）二维结构单元平板单元壳单元空间矩形空间三角形平板单元壳单元3 单元阶次一次单元：线性单元，只有角结点。二次单元：在角结点间的边界上配置一个边内结点。三次单元：边界上配置二个内结点，有时含内部结点。4 特殊单元弹簧单元、阻尼单元、间隙单元、界面单元、刚体单元、集中质量单元等。无限元模拟裂纹的奇异单元4.2 一维单元插值函数 1 拉格朗日（Lagrange）一维单元1)插值函数已知函数f(x)在[a,b]上一系列点上的值用一较简单的函数代替 f(x)，使称 为插值函数。若 为多项式，称为多项式插值函数。可由 确定[a,b]上其它点的函数值。例子：2点插值线性组合2) Lagrange插值公式设函数f(x)在[a,b]上一系列互不相同的点上的值构造一个（n-1）次的多项式讨论：具有形状函数Ni(x)的性质，可以作为有限单元的插值函数。ξξ1xn10-1x13) Lagrange单元自然坐标作变换：x2一维 n结点单元插值0可得：可得：作变换：N11N2121ξ2(+1)1(-1)21例子：一次单元插值等参变换ξ2(1)1(-1)例子：一维3结点Lagrange插值函数二次单元插值3(0)ξ1-12 埃尔米特（Hermite）一维单元 若场方程的最高阶导数为m，要求场函数及其m-1阶导数在结点上均连续。如梁、板、壳结点的挠度和转角均要连续，转角即为挠度的一阶导数关系。2结点单元：假设Hi为三次多项式（i=1~4）类似地，可写出其它三组方程，解方程得到所有待定参数（i=1~4） Hermite插值函数H1在结点1处为1.0，在结点2处为0.0；H3在结点1处为0.0，在结点2处为1.0；H2的导数在结点1处为1.0，在结点2处为0.0；H4的导数在结点1处为0.0，在结点2处为1.0。4.3 二维单元插值函数 1 三角形单元1(1,0,0)1）面积坐标L2=0L3=0A3A2P特点：A1① 在本点为1，其它点为03(0,0,1)L1=02(0,1,0)②③ 与三角形的具体形状及其在总体坐标系中的位置无关。自然坐标2）与直角坐标间的转换关系结论： 三角形域内点的坐标和位移可以通过相同的插值函数分别由结点的坐标和位移得到。等参变换3）面积坐标的微分运算4）面积坐标表示的插值函数插值函数的构造式通过除结点 i 以外所有结点的直线方程的左端项划线法直线方程在结点 i 的取值一次单元3 (0,0,1)L1=0L2=02(0,1,0)L3=01 (1,0,0)二次单元3 (0,0,1)L3-1/2=0L1=0L2=05(0,1/2,1/2)6 (1/2,0,1/2)L2-1/2=0L1-1/2=02(0,1,0)L3=04 (1/2,1/2,0)1 (1,0,0)2 拉格朗日（Lagrange ）矩形单元一次单元二次单元三次单元插值函数变化规律缺点：随着插值函数方次的增高而增加内结点，从而增加自由度。3 Serendipity矩形单元划线法1）4结点单元同理：统一形式：双一次单元与Lagrange一次单元相同2）变结点单元可以构造变结点的过渡单元。可以构造5-8结点单元。3）8结点单元4.4 三维单元 1 四面体单元利用体积坐标：4一次单元：4结点1038976二次单元：10结点152角结点：棱内结点：2 拉格朗日（Lagrange）六面体单元沿三个方向分别有r、p、q个结点，可以分别构造（r-1）阶、（p-1）阶、（q-1）阶的Lagrange插值基函数：典型的三维Lagrange单元3 Serendipity六面体单元一次单元：8结点ζη81972018561715161314次单元：20结点129角结点：典型的棱内结点4 三棱柱单元Serendipity类单元插值函数由三角形单元插值函数与Serendipity单元插值函数相乘得到一次单元Lagrange插值单元插值函数由三角形单元插值函数与Lagrange插值函数相乘得到Serendipity类一次单元二次单元Li为面积坐标三棱柱Serendipity单元4.5 阶谱单元 自适应分析：p方案：网格不变，提高单元的阶次后重