有限元试卷答案.doc

1、1. 写出线弹性平面问题三类基本方程和二类边界条件（分量或指标形式），并指出相应的自变量。答：三个基本方程 平衡方程 本构方程 平面应力平面应变 应变协调方程 几何方程 二类边界条件 力的边界条件 位移边界条件 如今给定的位移边界为，则有（在），其中分别为边界上x，y方向上的位移分量 2. 简述有限元法分析的基本步骤和相对应的基本表达式。答：江见鲸的教材P26页（4步的这是），如果想要写完成7步的话那就是P8步骤： 将结构离散化 单元分析，求得单元节点位移与节点力的关系，计算单元刚度矩阵 以节点为隔离体，建立平衡方程 施加荷载 引入边界条件 求解方程，求得节点位移 对每一单元循环，由单元节点位

2、移通过单元刚度矩阵求得单元应力或杆件内力 表达式： 位移模式 几何矩阵 B= 弹性矩阵 应力矩阵 3. 单元分析中,假设的位移模式应满足哪些条件,为什么?平面三结点三角形单元中，能否构造如下的位移模式(说明原因)答：这类问题参照江见鲸的教材P13，要求满足的三个条件不能。因为不满足完备性，缺少表示刚体位移的常数项和表示应变是位移一阶导数的常应变项不能保证解的收敛性。4. 简述加权余量法、半解析法、样条有限元法、边界单元法的特点。 答：加权余量法：当n有限时，定解方程存在偏差（余量），取权函数，强迫余量在某种平均意义上均为采用使余量的加权积分为0的等效积分以“弱”形式来求得微分方程近似解的方法。

3、半解析法：离散与解析相结合的方法，减少计算工作量，降低费用。样条有限元法：具有紧凑型及良好的光滑性，明确的表达式的优点，所得到的结果均在单元节点上，在数据的后处理方面更为方便和精确。边界单元法：将所研究问题的偏微分方程，设法转换为在边界上定义的边界积分方程，然后将边界积分方程离散化为只含有边界结点未知量的代数方程组，解此方程组可得边界节点上的未知量并可由此进一步求得所研究区域中的未知量，它除了能处理有限元方法所适应的大部分问题外，还能处理有限元法不易解决的无限域问题。5. 验证3结点三角形单元的位移插值函数满足证明：由原题目所知得又因为故，得结论。这个题参照教材12页到14页6. 推导一维杆单

4、元的形函数、几何矩阵、应力矩阵、刚度矩阵答：形函数参照教材11到14页，几何矩阵在15页，应力矩阵在17页，刚度矩阵在19页。具体不是很会，大家自行解答吧7. 就下列言论写出自己的看法：（有限元分析之大腕版）一定要选最变态的题目，什么材料非线性啊，几何非线性啊，接触非线性啊，多物理场耦合啊，都给他弄进去。是个模型就几百万个单元，上千万个节点，画个剖面图就要十几个小时。再整一并行机群，TOP500的，张口就是 High Performance Computation，一口地道的伦敦腔，倍儿有面子。题目一扔进去就跑个把月，你要是一个星期以内出结果，你都不好意思和别人打招呼。你说这样一趟算下来要发多

5、少Paper? 10篇？10篇？就1篇！你还别嫌少，说不定人家还发在会议上。你得琢磨牛人的心理啊，有能耐算这样题目的人，根本就不在乎多发一篇两篇文章。什么叫大牛知道么？大牛就是不求灌水，但求经典。答：这道题。谁知道想问什么。8. 采用矩形薄板单元（4个结点，12个结点位移分量，1个挠度独立变量）计算受中心集中力的四边支承板，位移函数：计算结果如下表(边长为1，厚度为0.01，弹模为1，波松比为0.3)单元数（1/4板）四边固定板中心挠度wD/PL2边中点弯矩M/P220.00614-0.1178440.00580-0.1233660.00571-0.1245理论解0.00560-0.1257试

6、分析本题中有限元解位移大于解析解、弯矩小于解析解的原因。答：有限元解位移大于解析解的原因是单元为非完全协调单元。 挠度w是弯曲问题中的基本未知函数且由于忽略了Z方向的变化，因此它只是x,y的函数：，若w已知，则唯一、内力、应力均可按上述相应公式求出。在经典解析法中，w（x，y）常设为三角级数形式。9. 三角形单元的结点坐标如图所示，设单元中一点A的坐标为（0.6，0.3）。已知三角形三结点单元的i 结点位移为（2.0,1.0），j 结点位移为(2.1,1.1)，m 结点位移为(2.15,1.05)。1）写出单元的位移函数；2）求A 点的位移分量。答：设位移函数为，将I,j,m点坐标代入u,v,

7、最终得： 将A点坐标代入位移函数得 即 A（2.095，1.045）。10、已知外力P=50103N，确定位移场、应力、应变和支座反力。取E=2.1104N/mm2。答：结构刚度矩阵K为：整体荷载列阵为：引入边界条件Q1=0和Q3=1.2mm，修正方程如下：解得位移场：单元应力、应变为：支座反力为：另一套的部分题（重复和不会的就被我忽略了。）1.各学科对有限元法的定义略不相同,但有一条相同的,即都认为有限元法是有效的_数值分析_方法. 2. 求解区域的离散化:即把实际问题,的整体求解域剖分为有限个单元,这可以通过有限个_离散_的三角形,四边形或多边形来分割求解区域. 3. 在划分单元时,任意一

8、个三角形单元的_节点_,必须也是相邻单元的_节点_,而不是相邻单元边上的_点_. 二、简答题（40分）3. 单元分析中,假设的位移模式应满足哪些条件,为什么? 答：见教材P134. 试阐述总体刚度矩阵的特点.(这里问总刚在教材P23，如果问你单刚那就是在P20)答：整体刚度矩阵是对称矩阵 整体刚度矩阵的主对角线上的元素总是正的 整体刚度矩阵是一个稀疏阵 施加荷载没有支承的整体刚度矩阵是一个奇异阵答： 理由：单元刚度矩阵不随单元（或坐标轴）的平行移动或作n元（n为整数）角度的移动而改变三、（15分）如下常微分方程du/dx+u=0 x=(0,1),边界条件为:u|x=0=1,设满足边界条件的试函数u(x)=1+cx,试用Galerkin法确定c值.答：解：利用Galerkin法。把试函数取为权函数残值，又由，得。