高等有限元课后题答案
- 格式:doc
- 大小:370.00 KB
- 文档页数:19
课程名称有限元单元法学院机械工程学院专业机械制造及其自动化学生姓名周祥态学号 110010021 教师陈爱军理论部分习题及解答:1.已知平面应变下三节点三角形单元的节点坐标()0,0i 、()1,6j 和()4,3m ,单元的节点位移分量1==i i v u 、0====m m j j v u v u ,材料的弹性模量E ,泊松比为ν。
试求单元的应变和应力分量?解:单元面积:1100112111612221134i i j j mmx y A x y x y === 23j m m j ai x y x y =-= 0j a = 0m a =3i j m b y y =-=- 4j m i b y y =-= 1m i j b y y =-=-3i m j c x x =-=- 3j i m c x x =-=- 6m j i c x x =-=1(133)211(43)211(6)21i j m N x y N x y N x y ⎧=--⎪⎪⎪∴=-⎨⎪⎪=-+⎪⎩[]304010103030621333461B --⎡⎤⎢⎥∴=--⎢⎥⎢⎥----⎣⎦∴单元应变{}[]{}113040103011030306302121333461600e eB εδ⎡⎤⎢⎥⎢⎥---⎡⎤⎡⎤⎢⎥⎢⎥⎢⎥==--=-⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥----⎣⎦⎣⎦⎢⎥⎢⎥⎣⎦平面应变下:[]101(1)10(1)(12)112002(1)E D μμμμμμμμμ⎡⎤⎢⎥-⎢⎥⎢⎥-=⎢⎥+--⎢⎥⎢⎥-⎢⎥-⎣⎦{}[]{}311(1)3121(1)(12)17(1)(12)1231ee E E D μμσεμμμμμμμ⎧⎫-⎪⎪-⎪⎪⎧⎫⎪⎪---⎪⎪∴===⎨⎬⎨⎬+--+-⎪⎪⎪⎪-⎩⎭⎪⎪-⎪⎪-⎩⎭2.图示为一个三个节点的杆单元,O 为坐标原点,其位移模式取为2321x C x C C U ⋅+⋅+=。
有限元课后习题答案1.1有限元法的基本思想和基本步骤是什么首先,将表示结构的连续离散为若干个子域,单元之间通过其边界上的节点连接成组合体。
其次,用每个单元内所假设的近似函数分片地表示求解域内待求的未知厂变量。
步骤:结构的离散化,单元分析,单元集成,引入约束条件,求解线性方程组,得出节点位移。
1.2有限元法有哪些优点和缺点优点:有限元法可以模拟各种几何形状复杂的结构,得出其近似解;通过计算机程序,可以广泛地应用于各种场合;可以从其他CAD软件中导入建好的模型;数学处理比较方便,对复杂形状的结构也能适用;有限元法和优化设计方法相结合,以便发挥各自的优点。
缺点:有限元计算,尤其是复杂问题的分析计算,所耗费的计算时间、内存和磁盘空间等计算资源是相当惊人的。
对无限求解域问题没有较好的处理办法。
1.3有限元法在机械工程中有哪些具体的应用静力学分析模态分析动力学分析热应力分析其他分析2.1杆件结构划分单元的原则是什么?1)杆件的交点一定要取为节点2)阶梯形杆截面变化处一定要取为节点3)支撑点和自由端要取为节点4)集中载荷作用处要取为节点5)欲求位移的点要取为节点6)单元长度不要相差太多2.2简述单元刚度矩阵的性质。
单元刚度矩阵是描述单元节点力与节点位移之间关系的矩阵。
2.3有限元法基本方程中每一项的意义是什么?{Q}---整个结构的节点载荷列阵(包括外载荷、约束力);{}---整个结构的节点位移列阵;[K]---结构的整体刚度矩阵,又称总刚度矩阵。
2.4简述整体刚度矩阵的性质和特点。
对称性奇异性稀疏性主对角上的元素恒为正2.5位移边界条件和载荷边界条件的意义是什么由于刚度矩阵的线性相关性不能得到解,从而引入边界条件。
2.6写出平面刚架问题中单元刚度矩阵的坐标变换式2.7推导平面刚架局部坐标系下的单元刚度矩阵。
2.8简述整体坐标的概念。
单元刚度矩阵的坐标变换式把平面刚架的所有单元在局部坐标系X’O’Y’下的单元刚度矩阵变换到一个统一的坐标系xOy下,这个统一的坐标系xOy称为整体坐标系。
有限元试题及答案 有限元试题及答案 一 判断题(20分)(×)1. 节点的位置依赖于形态,而并不依赖于载荷的位置(√)2. 对于高压电线的铁塔那样的框架结构的模型化处理使用梁单元 (×)3. 不能把梁单元、壳单元和实体单元混合在一起作成模型 (√)4. 四边形的平面单元尽可能作成接近正方形形状的单元(×)5. 平面应变单元也好,平面应力单元也好,如果以单位厚来作模型化 处理的话会得到一样的答案(×)6. 用有限元法不可以对运动的物体的结构进行静力分析 (√)7. 一般应力变化大的地方单元尺寸要划的小才好(×)8. 所谓全约束只要将位移自由度约束住,而不必约束转动自由度 (√)9. 同一载荷作用下的结构,所给材料的弹性模量越大则变形值越小 (√)10一维变带宽存储通常比二维等带宽存储更节省存储量。
二、填空(20分)1.平面应力问题与薄板弯曲问题的弹性体几何形状都是 薄板 ,但前者受力特点是: 平行于板面且沿厚度均布载荷作用 ,变形发生在板面内;后者受力特点是: 垂直于板面 的力的作用,板将变成有弯有扭的曲面。
2.平面应力问题与平面应变问题都具有三个独立的应力分量: σx ,σy ,τxy ,三个独立的应变分量:εx ,εy ,γxy ,但对应的弹性体几何形状前者为 薄板 ,后者为 长柱体 。
3.位移模式需反映 刚体位移 ,反映 常变形 ,满足 单元边界上位移连续 。
4.单元刚度矩阵的特点有:对称性 , 奇异性 ,还可按节点分块。
5.轴对称问题单元形状为:三角形或四边形截面的空间环形单元 ,由于轴对称的特性,任意一点变形只发生在子午面上,因此可以作为 二 维问题处理。
6.等参数单元指的是:描述位移和描述坐标采用相同的形函数形式。
等参数单元优点是:可以采用高阶次位移模式,能够模拟复杂几何边界,方便单元刚度矩阵和等效节点载荷的积分运算。
7.有限单元法首先求出的解是 节点位移 ,单元应力可由它求得,其计算公式为{}{}[][]eD B σδ=。
有限元习题及答案有限元习题及答案有限元方法是一种常用的数值计算方法,用于求解各种工程和科学问题。
在学习有限元方法的过程中,练习习题是非常重要的,可以帮助学生巩固所学的知识,并提高解决实际问题的能力。
本文将介绍一些有限元习题及其答案,希望对学习有限元方法的同学有所帮助。
习题一:一维热传导问题考虑一个长度为L的一维杆,其两端固定,杆上的温度满足以下热传导方程:∂²T/∂x² = 0,其中T为温度,x为位置。
已知杆的两端温度分别为T1和T2,求解杆上的温度分布。
解答一:根据热传导方程,可以得到温度分布的一般解为T(x) = Ax + B,其中A和B为常数。
根据边界条件,可以得到方程组:T(0) = B = T1T(L) = AL + B = T2解方程组可得A = (T2 - T1) / L,B = T1。
因此,温度分布为T(x) = ((T2 - T1) / L) * x + T1。
习题二:二维弹性问题考虑一个矩形薄板,其长为L,宽为W,材料的弹性模量为E,泊松比为ν。
已知薄板的边界上施加了一定的边界条件,求解薄板上的位移场。
解答二:对于二维弹性问题,可以使用平面应力假设,即假设薄板内部的应力只有两个分量σx和σy,并且与z轴无关。
根据平面应力假设和胡克定律,可以得到位移场的偏微分方程:∂²u/∂x² + ν * (∂²u/∂y²) + (1 - ν) * (∂²v/∂x∂y) = 0∂²v/∂y² + ν * (∂²v/∂x²) + (1 - ν) * (∂²u/∂x∂y) = 0其中u和v分别为位移场在x和y方向上的分量。
边界条件根据具体情况给定。
通过数值方法,如有限元方法,可以求解位移场的近似解。
习题三:三维流体力学问题考虑一个三维流体力学问题,流体在一个封闭容器内流动,容器的形状为一个长方体,已知流体的速度场和压力场的初始条件,求解流体的运动状态。
有限元试题及答案一、选择题1. 有限元法是一种数值方法,主要用于求解什么类型的数学问题?A. 线性代数方程B. 微分方程C. 积分方程D. 偏微分方程答案:D2. 在有限元分析中,以下哪项不是网格划分的基本原则?A. 网格应尽量均匀B. 网格应避免交叉C. 网格应尽量小D. 网格应适应几何形状答案:C3. 有限元方法中,单元的局部刚度矩阵可以通过以下哪种方式获得?A. 直接积分B. 矩阵乘法C. 线性插值D. 经验公式答案:A二、填空题1. 有限元方法中,______ 是指将连续的域离散化成有限数量的小单元。
答案:离散化2. 在进行有限元分析时,______ 是指在单元内部使用插值函数来近似求解场变量。
答案:近似3. 有限元法中,______ 是指在单元边界上满足的连续性条件。
答案:边界条件三、简答题1. 简述有限元法的基本步骤。
答案:有限元法的基本步骤包括:(1)定义问题域;(2)离散化问题域,生成网格;(3)为每个单元定义局部坐标系和形状函数;(4)组装全局刚度矩阵和载荷向量;(5)施加边界条件;(6)求解线性代数方程;(7)提取结果并进行后处理。
2. 描述有限元分析中的单元类型有哪些,并简述每种单元的特点。
答案:常见的单元类型包括:(1)一维单元,如杆单元和梁单元,特点是沿一个方向传递力;(2)二维单元,如三角形和四边形单元,特点是在平面内传递力;(3)三维单元,如四面体和六面体单元,特点是在空间内传递力。
每种单元都有其特定的形状函数和刚度矩阵。
四、计算题1. 给定一个简单的一维弹性杆问题,其长度为L,两端固定,中间施加集中力P。
使用有限元法求解该杆的位移和应力分布。
答案:首先,将杆离散化为一个单元。
使用一维杆单元的局部刚度矩阵和形状函数,可以推导出全局刚度矩阵。
然后,施加边界条件,即杆的两端位移为零。
最后,将集中力P转换为等效节点载荷,求解线性代数方程,得到节点位移。
应力可以通过位移和杆的截面特性计算得出。
有限元课后第三章习题答案有限元课后第三章习题答案第一题:根据题目给出的信息,我们可以得出以下结论:1. 题目中提到了一个平面问题,即只考虑二维情况。
2. 材料的弹性模量为E = 210 GPa。
3. 材料的泊松比为ν = 0.3。
4. 材料的厚度为t = 10 mm。
5. 材料的长度为L = 100 mm。
6. 材料的宽度为W = 50 mm。
7. 材料的边界条件为固定边界。
根据以上信息,我们可以开始解题。
首先,我们需要确定有限元模型的几何形状和单元类型。
由于题目给出的是一个平面问题,我们可以选择使用二维平面应力单元来建模。
根据题目给出的材料尺寸,我们可以选择一个矩形区域作为有限元模型的几何形状。
接下来,我们需要确定有限元模型的单元划分。
由于题目没有给出具体的单元划分要求,我们可以根据经验选择适当的单元尺寸和划分密度。
在这里,我们可以将矩形区域划分为若干个等大小的四边形单元。
然后,我们需要确定有限元模型的边界条件。
根据题目给出的信息,材料的边界条件为固定边界。
这意味着模型的边界上的节点在计算过程中将保持固定位置,不发生位移。
因此,我们需要将边界上的节点固定。
接下来,我们可以开始进行有限元计算。
首先,我们需要确定有限元模型的节点和单元编号。
然后,我们可以根据材料的弹性模量和泊松比,以及节点和单元的位置信息,计算出每个节点和单元的刚度矩阵。
然后,我们可以根据边界条件,将固定边界上的节点的位移设置为0。
这样,我们就可以得到一个由位移未知数构成的线性方程组。
通过求解这个线性方程组,我们可以得到模型中每个节点的位移。
最后,我们可以根据节点的位移和单元的刚度矩阵,计算出每个单元的应力和应变。
根据题目给出的材料厚度,我们可以得到每个单元的应力和应变的平均值。
综上所述,根据题目给出的信息,我们可以使用有限元方法来求解这个平面问题。
通过建立有限元模型,确定边界条件,进行有限元计算,我们可以得到模型中每个节点的位移和每个单元的应力和应变。
1.1 有限单元法中“离散”的含义是什么?有限单元法是如何将具有无限自由度的连续介质问题转变成有限自由度问题的?位移有限元法的标准化程式是怎样的?(1)离散的含义即将结构离散化,即用假想的线或面将连续体分割成数目有限的单元,并在其上设定有限个节点;用这些单元组成的单元集合体代替原来的连续体,而场函数的节点值将成为问题的基本未知量。
(2)给每个单元选择合适的位移函数或称位移模式来近似地表示单元内位移分布规律,即通过插值以单元节点位移表示单元内任意点的位移。
因节点位移个数是有限的,故无限自由度问题被转变成了有限自由度问题。
(3)有限元法的标准化程式:结构或区域离散,单元分析,整体分析,数值求解。
1.3 单元刚度矩阵和整体刚度矩阵各有哪些性质?各自的物理意义是什么?两者有何区别?单元刚度矩阵的性质:对称性、奇异性(单元刚度矩阵的行列式为零)。
整体刚度矩阵的性质:对称性、奇异性、稀疏性。
单元 Kij 物理意义 Kij 即单元节点位移向量中第 j 个自由度发生单位位移而其他位移分量为零时,在第 j 个自由度方向引起的节点力。
整体刚度矩阵 K 中每一列元素的物理意义是:要迫使结构的某节点位移自由度发生单位位移,而其他节点位移都保持为零的变形状态,在所有个节点上需要施加的节点荷载。
2.2 什么叫应变能?什么叫外力势能?试叙述势能变分原理和最小势能原理,并回答下述问题:势能变分原理代表什么控制方程和边界条件?其中附加了哪些条件?(1)在外力作用下,物体内部将产生应力σ和应变ε,外力所做的功将以变形能的形式储存起来,这种能量称为应变能。
(2)外力势能就是外力功的负值。
(3)势能变分原理可叙述如下:在所有满足边界条件的协调位移中,那些满足静力平衡条件的位移使物体势能泛函取驻值,即势能的变分为零δ∏p=δ Uε+δV=0此即变分方程。
对于线性弹性体,势能取最小值,即δ2∏P=δ2Uε+δ2V≥0此时的势能变分原理就是著名的最小势能原理。
2 弹性力学问题的有限单元法思考题2.1 有限元法离散结构时为什么要在应力变化复杂的地方采用较密网格,而在其他地方采用较稀疏网格?答:在应力变化复杂的地方每一结点与相邻结点的应力都变化较大,若网格划分较稀疏,则在应力突变处没有设置结点,而使得所求解的误差很大,若网格划分较密时,则应力变化复杂的地方可以设置更多的结点,从而使得所求解的精度更高一些。
2.2 因为应力边界条件就是边界上的平衡方程,所以引用虚功原理必然满足应力边界条件,对吗?答:对。
2.3 为什么有限元只能求解位移边值问题和混合边值问题?弹性力学中受内压和外压作用的圆环能用有限元方法求解吗?为什么?答:有限元法是一种位移解法,故只能求解位移边值问题和混合边值问题。
而应力边值问题没有确定的位移约束,不能用位移法求解,所以也不能用有限元法求解。
2.4 矩形单元旋转一个角度后还能够保持在单元边界上的位移协调吗?答:能。
矩形单元的插值函数满足单元内部和单元边界上的连续性要求,是一个协调元。
矩形的插值函数只与坐标差有关,旋转一个角度后各个结点的坐标差保持不变,所以插值函数保持不变。
因此矩形单元旋转一个角度后还能够保持在单元边界上的位移协调。
2.5 总体刚度矩阵呈带状分布,与哪些因素有关?如何计算半带宽? 答:因素:总体刚度矩阵呈带状分布与单元内最大结点号与最小结点号的差有关。
计算:设半带宽为B ,每个结点的自由度为n ,各单元中结点整体码的最大差值为D ,则B=n(D+1),在平面问题中n=2。
2.6 为什么单元尺寸不要相差太大,如果这样,会导致什么结果? 答:由于实际工程是一个二维或三维的连续体,将其分为具有简单而规则的几何单元,这样便于网格计算,还可以通过增加结点数提高单元精度。
在几何形状上等于或近似与原来形状,减小由于形状差异过大带来的误差。
若形状相差过大,使结构应力分析困难加大,误差同时也加大。
2.7 剖分网格时,在边界出现突变和有集中力作用的地方要设置结点或单元边界,试说明理由。
答:有限元处于弹性力学问题的方法是离散法。
它将一个受外力作用的连续弹性体离散成一定数量的有限小的单元集合体,单元之间只在结点上相互联系,即只有结点才能传递力。
所以在边界出现突变和有集中力作用的地方要设置结点和单元边界。
2.8 为什么说三角形三结点单元是常应变单元,如果在每边中点增加一个结点,那么单元内应力如何分布?答:(1)应变矩阵[B]中的参数m j i m j i c c c b b b 、、、、、由坐标变量x 、y 之差确定。
当单元的坐标差确定之后,这些参数与坐标变量x 、y 无关,因此[B]为常量阵。
当单元的结点位移{a}确定后,由[B]转换求得的单元应变都是常量,也就是说在荷载作用下单元中各点具有统一的xy y x γεε、、值。
因此三结点三角形单元称为常应变单元。
(2)如果在每边中点增加一个结点,单元内的应力为线性分布。
习 题2.1试证明x 、y 与面积坐标的关系 证明:设P 点坐标为(x,y ) 同理可求得:由面积坐标定义得:由此推出坐标y x 、与面积坐标的函数关系:()()22j i i j j i i j i j j i m j j m m j m j m j j m A c L c L a c a c x b c b c A b L b L a b b a y b c b c ⎧-+-=⎪-⎪⎨-+-⎪=⎪-⎩式(2.1)面积:代入式(2.1)有: 其中形状参数由下式确定: 代入上式(2.1)可转化为:再加上 m j i L L L ++=1所以用面积坐标表示直角坐标矩阵形式如下: 2.2 试证明两相似三角形的单元刚度矩阵相同。
证明:由于两个三角形相似,故设h A A =21, h 为一常数。
三角形:()111121i j ji i c b c b A -=参数Λj i j i c c b b 、、、,只与坐标差有关,所以单元刚度矩阵通式为: 故[][]21rs rs K K =所以因此两相似三角形的单元刚度矩阵相同。
2.3 直角三角形固定在刚性基础上,受齐顶的水压力和自重作用,如图2.14所示。
若按一个单元计算,水的容重g γ,三角形平面构件容重g ρ,取泊松比v =1/6,试求顶点位移和固定面上的反力。
解:按逆时针编码,局部编码与整体编码相同:1-2-3建立坐标())0,0(3)3,0(20,21:a a xoy(1) 求形函数矩阵: 图(2.14) 形函数: 所以: 形函数的矩阵为:(2) 刚度矩阵 可得:(3)位移列向量和右端项 由边界条件可确定:{}{}Teu a 000022υ=水压力和构件厚分别为: 自重为W 与支座反力:所以:{}Ty x y x eW R h q R W h q W R R R ⎭⎬⎫⎩⎨⎧-+--=33363303011由[]{}{}eeeR a K =得到下列矩阵方程组:化简得:可得:将⎭⎬⎫⎩⎨⎧22υu 代入下式: 固定面上的反力:a h ga gh q 330===γγ从而可得支座反力为:2.4 试从式(2.69)说明对角线元素改1法只能用于给定零位移的情形,而对角元素充大数法可以适用任意指定位移的情形。
解:(1)在式(2.69)中, 采用对角线元素改一法,则:所以在式(2.69a )中可改为:需满足:{}{}a b a a ''=因为 :ab K ]'[为非零矩阵 所以 :{}b a '应需为向量则这种方法只能用于给定零位移。
(2)同上在式(2.69a )中{}{}a K R aa a''='α (其中α可取2010左右甚至更大的量级)根据主元充大数:{}a aa a K ']'[α》》{}b ab a K ']'[由于所以右边近似等于{}a aa a K ']'[α 所以左边等于右边即主元充大数法满足任何给定位移(零值和非零值)2.5 仿照例2.3,试证明矩形单元是完备协调元。
..证明:在平面问题中,场函数为位移场函数,4结点矩形单元的位移插值函数为:平面问题的泛函数中所包含的物理量的最高导数为2阶,而位移插值函数中包含的位移场函数(考虑到对称性,位移插值函数中只取了xy 项)及其直至2阶的完全多项式,所以是完备的。
要证明协调性,只需考虑单元边界上的连续性。
为此,参考上图,对于j-m 边的公共结点j ,m 的位移j u ,m u 完全确定,所以在边界上是协调的。
2.6 仿照例2.3,试证明图2.15所示任意四边形单元在选取位移模式 是完备的但不是协调的。
证明:完备性证明同2.5;协调性:对于边i-j ,直线方程为Bx A y +=所以代入上式,得:公共边i-j ,i ,j 两点的位移不能完全确定A ',B ',C '三个为未知量故任意四边形单元是完备的但不是协调的。
2.7上端受均布荷载作用的墙体视为平面构件,如图2.16(a)所示,划分成为8个单元如图2.16(b),试在图2.16(b)上画出边界条件和等效结点荷载;倘若利用对称性取一半,如图2.16(c),试在图2.16(c)上画出边界条件和等效结点荷载。
设E=30Gpa ,v =1/6,均布荷载q=20KN/m ,试按图2.16(c)求构件的结点位移。
3单元和插值函数思 考 题3.1什么是面积坐标?如何计算三角形内某点的面积坐标?答:(1)如(a )图所示,三角形内任一点P (x ,y ),将P 与三角形三个顶点i ,j ,m 连成3个三角形。
令i A 为i 点所对应三角形pjm 的面积,j A 为j 点所对应的三角形pmi 的面积,m A 为m 点所对应的三角形pij 的面积,面积坐标定义为:r L =r A /A (i ,j ,m ),其中A 为三角形ijm 的面积,点p (x ,y )用面积坐标可以写为P (i L ,j L ,m L ),且i L +j L +m L =1。
(2)求某点面积坐标除用定义外,还可用如图(b )所示的方法,即三角形内某点的面积坐标可通过同底三角形的高度比来计算。
如图(b )中的i L =h i /i H 。
(a ) (b) 图3.3 面积坐标3.2 什么是划线法?如何用划线法形成单元的插值函数?答:(1)划线法是根据形函数的0-1特性,将需要等于零的各结点用直线连接起来(划线);(2)在该直线上为零,则在该直线上的各结点的值也为零,为此形函数一定包含了该直线方程的因子,将需要等于零的各个因子乗起来即得到该单元的行函数。
3.3 下列平面单元的位移具有连续性吗? (1)平面三角形二次单元;连续(2)平面三角形三次单元; 连续(3)8结点矩形单元; 连续(4)8结点任意四边形单元。
连续3.4 下列单元满足收敛的充分必要条件∑Ni=1吗? (1)平面三角形三次单元; 满足 (2)变结点单元; 满足(3)长方体20结点单元。
满足3.5 对于非协调的薄板单元如何进行分片检验?答:当赋予单元片各个结点以与常应变状态相应的位移值和载荷值时,校验0)(m1=-∑=ei je e ij P a K 是否满足,如能满足则认为通过分片检验。
3.6在平面壳单元中如何判别共面点?可用什么方法进行处理?答:(1)在平面壳体单元中,如果某一点的各个单元面法向不同,经局部坐标转化到整体坐标后,该点的总体位移有6 个,若方向相同,常称此点为共面点。
(2)处理方法有两种:i 、在局部坐标系内建立结点平衡方程,并删去zi θ方向的平衡方程,于是剩下的方程满足唯一解的条件。
ii 、在此结点上,给一任意的的刚度系数z k θ,这时在局部坐标系中,此结点在zi θ方向的平衡方程经变换后,总体坐标中的系统方程满足唯一条件,它不影响单元应力。
习题3.1试利用面积坐标构造10结点三角形单元(图3.22)的9、10结点的插值函数。
解:利用划线法可得:)31(13199-=L L L N α 即 )3132(323219-⨯⨯⨯=α所以 )13(29)31(2271311319-=-=L L L L L L N 所以 3211027L L L N = 3.2利用构造变结点数单元插值函数的方法,构造图3.22所示三次三角形单元的插值函数,并和式(3.5)的结果进行比较。
解:由划线法可得设1N 为原三结点三角形的形函数,即11L N = 其余点在4结点的形函数均为0,32)4(1=L 所以 13204⨯'+=α 324-='α0)5(1=N ,15=N , 其余点在5结点的形函数均为0,31)5(1=L 1(8)0N =,81N =,其余结点在8点的形函数为0,1(8)13L =1(9)0N =,91N =,其余点在9结点的形函数均为0,1(9)23L =∴923α'=- 同理 1013α'=- 故 109854113132313132N N N N N L N -----= 以此类推经与式(3. 5)比较,所得结果相同。