有限元方法理论及应用

《有限元方法理论及应用》由会员上传分享，免费在线阅读，更多相关内容在学术论文-天天文库。

1、机械工程学院研究生研究型课程考试答卷课程名称：有限元方法理论及应用考试形式：□专题研究报告□论文□大作业□√综合考试学生姓名：学号：学生联系方式：导师：序号分项类别分值1理论知识探讨运用写作2理论分析与计算3上机实验报告4总分分1等参单元及其应用31.1概述等参单元的原理及其对有限元法工程应用的意义。31.2等参单元的数值积分方法41.3线性等参单元和非协调元101.4等参单元的应用132分析与计算152.1计算题一152.2计算题二162.3计算题三182.4计算题四223上机实验283.1第一题283.1.1实验题目283.1.2实验目的293.1.3建模概述293.1.4计算结果分析与

2、结论303.1.5实验体会与总结503.2实验二503.2.1实验题目503.2.2实验目的513.2.3建模概述513.2.4计算结果分析与结论533.2.5实验体会与总结553.3实验三553.3.1实验题目553.3.2实验目的553.3.3建模概述553.3.4计算结果分析与结论563.3.5实验体会与总结591等参单元及其应用1.1概述等参单元的原理及其对有限元法工程应用的意义。平面三角形单元、平面四面体单元、三维六面体单元这些单元受到两个方面的约束：其一是单元精度的约束，节点数越多，精度越高。其二是工程中的问题往往是复杂的几何体，规则的六面体和四面体不能准确地描述，且上述单元都是

3、直线边界，处理曲边界几何误差大。为了解决上述矛盾，可以使其成为任意四边形和任意六面体单元，显然，由于它已经不再是规则的四边形和六面体，所以它们的单元位移模式和形函数也不同于规则的四边形和六面体的形函数。为此必须引入所谓的等参变换。采用相同的插值函数对单元的节点坐标和节点位移在单元上进行插值，这种单元称为等参单元。等参单元的原理是通过等参变换，建立起局部（自然）坐标中几何形状规则的单元与总体（笛卡尔）坐标中几何形状扭曲的单元的一一对应的映射关系，以满足对一般形状求解域进行离散化的需要。为了得到上述映射的数学表达，引入对母单元节点上x,y,z坐标进行插值的思想，将母单元上每一点对应的x,y,z坐

4、标看成是对节点坐标的插值，插值函数与位移插值中的形函数相同：这样就得到了一个事实上的映射，n是节点总数，节点数越多，单元精度越高，是形状函数。通过上式建立起两个坐标系之间的变换，从而将自然坐标内的形状规则的单元变成为总体笛卡尔坐标内的形状扭曲的单元，通常称前者为母单元，后者为子单元。由于该几何变换式中采用了与位移模式相同的插值函数，因此称为等参变换。工程中一些结构的形状有的是比较复杂的、不规则的，有的具有曲边边界，如果用一般单元分析些类结构，需要划分大量的网格，取更多的节点，这样一来计算增大很多，而且处理曲边边界几何体误差也较大。对此，可以用等参单元来解决。等参单元具有曲面形状，可以用较少的

5、单元拼成复杂不规则的实际结构，大大减少计算量，同时也提高了计算精度。等参单元的优点如下：1、等参单元形状、方位任意，容易构造高阶单元，适应性好，精度高。2、等参单元列式具有统一的形式，规律性强，采用数值积分计算，程序处理方便。1.2等参单元的数值积分方法等参单元刚度矩阵的数值积分方法及确定积分阶的原理。全积分、减缩积分单元讨论和评价。对等参单元数值积分原理的讨论：在等参单元推求载荷向量或刚度矩阵是，需要进行如下形式的积分：其中被积分函数一般比较复杂，有的可以积分出结果，但式子很繁；有的甚至的不到它的显式。因此，一般都用数值积分代替函数积分，即，在单元内选出某些点，称为积分点，算出被积函数在这

6、些积分点处的函数值，然后用对应的加权系数，乘以这些函数值，在求出总和，将其作为近似的积分值。1.2.1一维数值积分设有积分式。现讨论它的积分。首先构造一个多项式上有，然后用近似函数的积分来近似原被积函数的积分，称为积分点或取样点。积分点的数目和位置决定了近似的程度，因而也就决定了数值积分的精度。（1）Newton-Cotes积分此种积分包括积分域端点在内的积分点按等间距分布。对于个积分点（或取样点），根据积分点上的被积函数值可以构造一个近似多项式，使在积分点上有：（1）这个近似多项式可以通过Lagrange多项式来表示。（2）其中是阶Lagrange插值函数，即有（3）由于Lagrange插

7、值函数有以下性质:（4）所以有（3）的插值函数是次多项式，因此近似函数也是次多项式，积分（5）令（6）则（3）可以写作（7）式中称为积分的权函数，可由（6）确定。可以看到加权系数与被积函数无关，只与积分点的个数和位置有关。Newton-Cotes积分中积分点的位置按等间距分布，即（=0,1,2,…）（8）其中是积分点间距。`现在用近似，根据式（7）可以写成（9）式中是余项。由于个积分点的Newton-Cote