有限元分析期末考试题目

1、试说明有限单元法解题的主要步骤网格划分：先将弹性体划分为有限个单元组成的离散体，单元之间通过单元节点相连接；单元分析：建立各个单元的节点位移和节点力之间的关系式；整体分析：对由各个单元组成的整体进行分析，建立节点外载荷与结点位移之间的关系式，以求解结点位移。2、单刚和总刚各有什么特征单刚ek][是对称的，反映出单元抵抗这种变形的能力，里面的每一个元素ijk表示的含义为：当j号节点位移分量为1，且其他节点位移分量皆为零时，对应的i号节点力分量。总刚meekK1][][是各单元刚度的总和（叠加），也是对称的，并且是稀疏的，呈带状分布。即整体刚阵内有很多的零元素，且非零元素都集中在对角线附近。它反映出整个结构抵抗这种变形的能力。3、刚阵中，每一项元素的物理意义是什么单刚ek][里面的每一个元素ijk表示的含义为：当j号节点位移分量为1，且其他节点位移分量皆为零时，对应的i号节点力分量。4、三角形三节点形函数[N]的特征（1）形函数在各单元节点上的值具有“本点为1，它点为零”的性质，即kikiyxNkki01),(),,(轮换mji（2）在单元内任意点上，三个形函数之和等于1，即1mjiNNN（3）三角形任意单元一条边上的形函数，仅与该边的两端节点坐标有关，在ij边上：ijiixxxxyxN1),(，ijijxxxxyxN),(，0),(yxNm5、为了满足收敛性条件，位移模式满足哪些条件（1）相容性：所选定的函数在整个求解域内有一定的连续性，即形状函数在单元内都是连续的，因为连续性要求只反映在单元之间；（2）完备性：为能实现求解函数的的任意可变性，选定的试验函数在在整个求解域内应能表现出任意可能的变化形式，即要求试验函数是完备的。（3）几何对称性：单元内的插值函数选定是应满足函数形式上的几何对称性。6、位移模式中哪些项反映单元的刚体位移和常应变（1）平面问题：位移模式：yaxaavyaxaau654321,中的常数项反映单元的刚体位移，一次项反映刚体的常应变。（2）薄板弯曲问题：位移模式31231131029283726524321xyayxayaxyayxaxayaxyaxayaxaaw多项式中的前三项反映出中面平板无弯曲的刚体位移。三个二次项经二阶微分后反映出中面变形的3中应变形式。7、试述边界条件处理的大数法原理和步骤以平面问题为例（1）固定约束：n2个位移分量中第i号位移分量i为0，则对应的方程被改为0)...(222211nniiiiiiKKKMK其中M为一个相当大的数（如3010M），则解出有0)(iiiMK相对小量（2）强迫约束：n2个位移分量中第i号位移分量ii，载荷列阵}{Q中对应项的未知约束改为iiMKQ，则第i行方程式改为iiinniiiiiiMKKKKMK)...(222211同理，M为一个相当大的数（如3010M），则解出有iiiiiiiiiMKMKMK)(相对小量步骤:引入约束条件后，结构的整体刚阵][K和载荷列阵}{Q都在相应节点的元素进行改变乘以一个相当大的数3010M，但其体积及编号不变，仍可记为}{}]{[QK。8、用最小势能原理推导}{}]{[QK单元应变能：ehdDUTe}]{[}{21，而eB}]{[}{，代入有eeTeeTTeekhdBDBUe}{][}{21}]{][[][}{21，其中hdBDBkeTe]][[][][为单元刚阵。而结构总应变能}]{[}{21}{][}{2111KkUUTemeeTemee单元外力功：meTTeedTudpuW1}{}{}{}{，将eNu}]{[}{代入有eTTeepTeTTeTTeQQdTNdpNWee}{}{}{}{}{][}{}{][}{其中eedTNQdpNQTeTTep}{][}{,}{][}{，则整个结构的外力功}{}{}{}{QQWTeTe.结构总势能：}{}{}]{[}{21}{}{}]{[}{211QKQkWUTTmeeTeeTe由最小势能原理0i),...,2,1(Ni可得节点位移方程}{}]{[QK9、平面三节点三角形单元中，能否选取如下的位移模式？为什么？yaaxayxu3120),(.12654),(yaxaayxv23221),(.2yaxyaxayxu26524),(yaxyaxayxv均不能选取上述方程为位移模式，因为均不满足位移模式的收敛性条件，这样无法得出正确的结果。10、已知直梁单元单刚为ek][（课本第7页）。试推导局部坐标下平面刚架的单刚2222346266126122646612612][lllllllllllllEJke平面刚架i节点的载荷列阵可表示为TiiiiMYXQ][}{，而e单元的6个节点位移分量和6个节点力分量均以列阵表示为Tjjjiiieff}{'TjjjiiiemqTmqTp}{'有材料力学可知，轴向位移只与轴向力T有关，弯曲位移f、只与弯曲力q、m有关，由直梁的弯曲变形关系，弯曲部分的节点位移与节点力关系为：jjiiebjjiiffkmqmq][式中的eebkk][][，对于轴向变形jiesjikTT][，由刚度矩阵的物理意义，当0,1ji时，lEAkTlEAkTsjijsiii,，同理可求得][][sjjsijkk和，则轴向变形刚度矩阵：1111][lEAks则局部做标下平面刚架的单刚为：jjjiiijjjiiifflEJlEJlEJlEJlEJlEJlEJlEJlEJlEJlEAlEJlEAlEJlEJlEJlEJlEAlEJlEAmqTmqT466122661226046061200060120022322322232311、平面三节点三角形单元应力应变和位移有什么特征？为什么？平面三节点三角形单元的位移是连续线性变化的，应变是常应变，如果单元是均质的，应力也是常应力。因为位移模式的选取为yaxaavyaxaau654321,，而单元应变的求解方程为：eNxyyx}]{[00}{，应力求解方程为}]{[}{D，而eNvu]][[，三节点三角形单元的位移模式中最高项为x、y的一次项，进行一次偏导后，应变值为常数，而如果单元是均质的，应力也将是常应力。12、平面三节点三角形单元，一次多项式位移模式yaxaau321yaxaav654试写出][N][B][Sek][的表达式nnmmllNNNNNNN000000其中222ycxbaNycxbaNycxbaNnnnnmmmmllll，为三角形单元的面积；nnnnmmmmllllbccbbccbbccbB00000021iiiiiiibcccbbES2121)1(2][2),,(nmli]][[][]][[][][BDBhdVBDBkTVTeennmmlyxyxanmlyyb11nmlxxc11nnllmyxyxanlmyyb11nlmxxc11mmllnyxyxamlnyyb11mlnxxc1113、上题，试分析其单元收敛性（1）相容性：1）单元内的形状函数都是x、y的连续函数有12题中的形状函数方程可以看出，lN、mN、nN均是x、y的一次函数，于是可得形状函数都是x、y的连续函数2）相邻单元交界面上，两单元应有相同的位移形状函数是连续的，设相邻单元（1）和（2）的节点号分别为a、b、c和a、d、b，则ab为交界边，变形前ab为一直线。变形后，移至'a、'b、'c、'd。由于单元内位移采用线性插值，单元（1）变形后仍是一条直线''ba；单元（2）变形后也是一条直线''ba，由于同一节点的位移式一样的，有两点确定一条直线可知两个单元的''ba边是重合的。（2）完备性：该位移模式能反映出此单元的三种独立的刚体位移及三种常应变形式，因此这种单元满足完备性要求。综上所述，该单元是收敛的。14、三角形单元等效节点力的计算公式，算一个具体数字的题目（体积力或表面力）15、哪个节点编号好？为什么？只需要满足半宽带越小越好16、轴对称问题怎样由三维问题转化为二维问题P684.3轴对称问题：几何形状对称于中心轴，载荷对称于中心轴，变形也对称于中心轴17、轴对称问题应力应变分量有几项应力应变分量各有4项，应力分量：Trzzr}{应变分量：rzzr}{18、轴对称问题的简单三角形单元是否是常应力，常应变？为什么？不是常应力和常应变。因为应变与位移分量的关系式为：wurzzrrrwzuzwrururzzr0010}{，这里面除含有微分算符外，还包含了r的倒数项r1，则即使位移模式为线性的，但由于该项的存在，使得应变与坐标有关，即不会是常应变。应力应变的物理关系为D，由于应变不是常应变，则所求得的应力也不会是常应力。19、轴对称问题等效节点力计算中的积分域有什么要注意的？积分域为子午截面上单元三角形面积域20、复杂函数的积分计算如何近似处理？常须计算三角形域内的面积积分，其形式为eSdrdzzrfF),(，一般情况下，上式应采用适当的数值积分。对简单的三角形单元，上述积分可直接近似求得。如l、m、n为单元的3个节点，当单元较小时，可取)(3),(nmlSfffdrdzzrfFelf、mf、nf为被积函数),(zrf在l、m、n3点的值。有时也有)(31nmlrrrr)(31nmlzzzz则eScfdrdzzrfF),(21、四节点的矩形薄板单元位移模式，并讨论其连续性31231131029283726524321xyayxayaxyayxaxayaxyaxayaxaaw22、等参单元的优点参数单元可以具有曲面形状，以便于拼成复杂的实际结构，也可以方便地构造多节点，高精度的参数单元，有利于复杂三维结构的有限元分析。而等参单元可以通过坐标变换将实际单元“映射”为一个正方形，拥有相同的形状函数，只是等参单元的形状函数不表示为直角坐标x、y的显函数，而表示为自然坐标、的函数。23、等参单元为何引入自然坐标等参单元引入自然坐标是为了把实际单元通过坐标变换“映射”为一个正方形24、求形函数表达式的方法（划线法）25、雅可比矩阵在等参元中的意义雅克比矩阵yxyxJ][为坐标变换矩阵，可以求出有整体坐标x、y到局部坐标的、关系26、采用何种方法进行积分（数值积分）采用高斯积分进行数值积分：111)()(niiiwxfdxxfI式中)(ixf为积分点ix处被积函数值，iw为权重，N为积分点的数目（P113）27、等参的含义用同样的形状函数和同样的节点插值表示出单元的几何坐标