东南大学计算力学习题及答案汇总(版).doc

Four short words sum up what has lifted most successful individuals above the crowd: a little bit more.-author-date东南大学计算力学习题及答案汇总(2011版)东南大学计算力学习题及答案汇总(2011版)第三章1如图所示一三角形钢板，两个结点固定，对第三个结点施以单位水平位移，测出所施加的力，从而得出相应的刚度系数。其他点依此类推，这样测得的刚度系数所组成的刚度矩阵，是否与按照常规三角形单元刚度矩阵计算公式所得结果一样？用这样实测所得的刚度矩阵能否进行有限元分析？为什么？解：不一样。单元刚度矩阵中每个元素的物理意义：表示单元第个自由度产生单位位移，其它自由度固定时，第个自由度产生的节点力。单元刚度矩阵是在单元处于平衡状态的前提下得出的，单元作为分离体看待，作用在它上面的外力（单元力）必是平衡力系，然而研究单元平衡时没有引入约束承受平衡力系作用的无约束单元，其变形是确定，但位移是不能确定的，即单元可发生任意的刚体位移。不能。因为与有限元中单元与单元之间的约束情况不一样，不能进行有限元分析。2以位移为基本未知量的有限元法其解具有下限性质，试证明之。解：系统总位能的离散形式将求解的方程带入可得在平衡情况下，系统总位能等于负的应变能。在有限元解中，由于假定的近似位移模式一般来说总与精确解有差别的。设近似解为、，真实解为、且根据最小势能原理，得到的系统的总位能总会比真正的总位能要大，故则则近似解的位移总体上小于精确解的位移解释如下：单元原是连续体的一部分，具有无限多个自由度，在假定了单元的位移函数后，自由度限制为只有以结点位移表示的有限自由度，引入了更多的约束和限制，使得单元刚度较实际连续体加强了，连续体的整体刚度随之增加，所以有限元解整体上较真实解偏小。3 请分别阐述单元刚度矩阵和整体刚度矩阵中任一元素的物理意义。解：在单刚中，表示单元第个位移产生一单位位移，其它位移为零时，第个位移方向上引起的节点力。在整体刚度中，表示第个自由度产生一单位位移，其它自由度为零时，第个自由度上引起的节点力。4 简述虚功原理，且使用虚功原理导出外荷载与节点荷载的等效关系式。解:虚功原理：变形体中任意满足平衡的力系在任意满足协调条件的变形状态上作的虚功等于零，即体系外力的虚功与内力的虚功之和等于零。设为外荷载（此处为体力），为节点荷载，为单元内位移场，为结点位移场根据虚功原理由于故则5试述弹性力学中按位移求解与有限单元法中按位移求解之间的异同点。解：弹性力学有限单元法物理模型连续体离散化结构基本方程几何方程物理方程平衡微分方程几何方程物理方程结点平衡方程解法解微分方程解代数方程解答形式用函数表示用数值表示解答精度精确解近似解6 如果三节点三角形单元绕其中某一个节点作小的刚体转动，其转角为，证明单元内所有的应力均为零。解：在三角形单元中由于三角形单元绕其中某一个节点作小的刚体转动，各节点的位移可表示为：则可知节点位移向量故应变由于弹性矩阵为常量矩阵，应变向量为零向量，故为零向量，即单元内所有的应力为零。7二维单元在坐标内平面平移到不同位置，单元刚度矩阵相同吗？在平面内旋转时又怎样？试证明之。解：二维单元在坐标内平面移到不同位置时，刚度矩阵相同。在平面内旋转时，刚度矩阵也相同。刚度矩阵单元平移或旋转时，不变，故单元刚度矩阵不变。8 判断有限元网格离散合理性 a) 对图1(a)所示的有限元网格，评论网格的优劣性，指出模型中的错误，并加以改正。b) 评论图1(b)的网格划分合理吗?为什么?请加以改正。图1解：（a）网格划分不合理。 1）无过渡单元 2）无边界条件 3）夹角区应力集中，应适当加密风格 4）对称结构网格应对称划分（b）不合理。 1）左部网格应适当加密2）由于三角形单元会造成局部精度不够，过渡区可采用其它单元划分3）右部单元的长宽比较大，就进行适当调整。9 如图2所示，平面三角形构件以x-y坐标系表示的刚度矩阵方程如下：试建立以,(与图中同向的位移)及,来表示的刚度矩阵方程。解：用坐标变换则其中 ， 由10 某平面结构采用四节点矩形单元和三节点三角形单元建立有限元计算模型，其如图3所示。试求结点2的等效荷载列阵。解：单元，荷载作用于边上，故等效节点力只与号节点有关形函数，在边上，则线性分布面力则单元，形函数 在1-2边上，故节点2的等效荷载列阵11 试求如图4所示的有限元网格的整体刚度矩阵，假设每个节点的自由度数为1，且设表示第e个单元的单元刚度矩阵（注意：结果应该用表示）。图5 图4 解：单元刚度矩阵，整体刚度矩阵：12 图5中两个三角形单元组成平行四边形，已知单元按局部编码的单元刚度矩阵K和应力矩阵S是按图5示单元的局部编码写出K，S。解：由图可知则由得到13如图6所示8结点矩形单元(每边中点为结点), 3点为坐标原点，a=b=2，单元厚为t。求该单元的位移函数和形函数和并检验其是否满足收敛性条件。求在2-6-3边作用均布水平荷载q时的等效结点荷载。解：(1)位移函数：引入无量纲的局部坐标则故则时，则角节点的形函数为边中节点的形函数为证明收敛性：位移函数中表示刚体位移，和表示常应变，故位移函数具有完备性设相邻单元公共边界上的直线方程是（或），代入位移函数中为（或）的2次函数，而边界上三点确定的位移函数为也为二次曲线，故单元在公共边界连续，故位移函数收敛(2)荷载作用在边上，故等效节点力只与号节点有关在边上计算 第四章1经典梁理论和Timoshenko梁理论有哪些相同点和哪些不同点?基于以上两种理论的梁单元各有何特性？解：经典梁理论Timoshenko梁理论相同点Kirchhoff假设不同点型单元弯曲梁单元截面转动是挠度的一阶导数，只有挠度是独立的采用Hermite插值型单元考虑剪切变形影响挠度和截面转动各自独立插值采用拉格朗日插值特性梁的高度远小于跨度梁很薄时，会造成剪切锁死现象2 写出杆件的应变能计算公式，并给出推导过程。解：将只考虑轴向变形的杆件划分成个单元，节点坐标为单元的位移函数 （）用形函数近似位移函数得，其中单元的应变单元的应力单元应变能其中3在杆系系统中,除了采用凝聚自由度的方法实现铰接端条件， 还有什么方法可以实现以上条件，并比较这几种方法的优缺点。解： 优点缺点凝聚自由度法4利用最小势能原理，推导图1所示弹性基础上梁单元方程，其中该梁的势能为：Lw(x)xkf图1解：根据最小势能原理可知故有对第一项分部积分则引入强制边界条件和自然边界条件使由于的任意性故控制微分方程为此梁的位移函数，则由于物理关系可知则由得则梁单元刚度方程为其中 5 图2所示刚架1) 如何进行节点编号使整体刚度矩阵K的带宽最小？2) 刚架的整体刚度矩阵中a节点的总刚度矩阵Kaa和的总刚度矩阵Kbc各由哪些单元的哪些分块矩阵叠加组成（自行确定单元局部坐标方向）3) 试按照二维等带宽存储和一维变带宽存储方式确定Kaa中对角元素的在相应存储数组中的位置。图2 有铰点的刚架解：1）考虑每个节点有两个自由度由于半带宽d=（相邻结点码的最大差值+1）*2故节点编号如图所示可使单元内节点编码相差3，使得带宽d=82） 3）考虑单元节点1自由度的凝聚可知中对角线元素在原整体刚度矩阵中第6行第7列和第7行第8列则用二维等带宽存储后在矩阵中的第6行第2列和第7行第2列用一维变带宽存储后在中对角线元素在数组中的位置为9和10 第五章1根据以下形函数表达式画出形函数N1和N3以及导数(dN2/dx)和(dN4/dx),它们代表梁单元整个长度上形状变化。2、对于图1中给出的四节点二次应变一维等参单元，试确定：a) 形函数N1，N2,N3，N4；b) 单元刚度矩阵k。图1解：（a）由拉格朗日插值函数可知， （b）一维问题中，单元刚度矩阵，3试利用变节点数法构造插值函数的，构造出图2所示的三次三角形单元的形函数及相应的位移函数。图2解：位移函数： （1）构造不考虑边节点和内部节点的角节点的插值函数：（2）构造不考虑内部节点的边节点的插值函数： ，（3）内部节点插值函数：（4）修正边中点的插值函数：同理得，（4）修正角节点的插值函数：4试构造如图3所示的15结点三棱柱体单元的插值函数，并判断其构造的位移函数是否收敛。解：（1）不考虑边中点构造三角形角节点的插值函数：，（2）构造三角形边中点的插值函数：，（3）构造四边形边中点的插值函数：，（4）修正角节点的插值函数：第六章1 等参元的收敛性证明。证明：（1）协调性：考察单元之间的公共边，为了保证协调性，相邻单元在这些公共边（或面）上应有完全相同的结点，同时每一单元沿这些边的坐标和未知函数应采用相同的插值函数加以确定。（2）完备性：三维等参元中， 有限元中，将场函数离散为各个单元局部场函数的集合体单元内场函数为则当时，表明单元能够表示线性变化的场函数，满足了完备性的要求。2 等参元的优点是什么？解：1）等参单元为协调元，满足有限元解收敛的充要条件2）将不规则单元转换为规则母单元后，容易构造位移函数和形函数3）当单元边界呈二次以上的曲线时，容易用很少的单元去逼近曲线边界3什么是位移的零能模式，在什么条件下会发生？如何检验它是否存在和如何防止它的出现。解：（1）由于采用减缩积分方案导致其应变能为零，而自身有别于刚体运动的位移模式称为位移的零能模式。（2）通过检查的非奇异性条件是否得到满足来验证是否存在零能模式。（3）高斯积分点提供应变分量的数目大于系统独立自由度数目，是保证系统刚度矩阵非奇异性的必要条件。系统不出现对应于除刚体运动以外位移模式的零特征值，是保证系统刚度矩阵非奇异性的充分条件。4 请阐说减缩积分概念,并分析其优缺点. 解:在数值积分中，能够保证不降低收敛速度的条件下求解各种条件有限元问题的最小阶次，比精确积分低阶的积分可称为减缩积分。一维问题刚度矩阵的积分中，如果插值函数中的多项式阶数为，微分算子中的导数的阶次是，则有限元得到的被积函数是次多项式。为了保证原积分的精度，选择高斯积分的阶次，可精确积分至次多项式，可达到精确积分刚度矩阵的要求。在二维单元和三维单元中仍按来确定积分阶次，即高斯积分阶数低于被积函数所有项次精确积分所需要阶数的积分方案，称为减缩积分。优缺点：（1）精确积分是由插值函数中非完全项的最高方次所要求，而决定有限元精度的通常是完全多项式的方次。这些非完全的最高方次项往往不能提高精度，反而带来不好影响。取较低阶的高斯积分，使积分精度正好保证完全多项式方次的要求，而不包括更高次的非完全多项式的要求，在一定情况下改善了单元的精度。（2）在最小位能原理基础上建立的位移有限元，位移解具有下限性质。有限元的计算模型具有较实际结构偏大的整体刚度。选取减缩积分方案使有限元计算模型的刚度有所降低，有助于提高计算精度。（3）采用减缩积分可能使系统刚度矩阵奇异，出现有别于刚体运动的位移零能模式。5 如需要对二维三次Serendipity单元进行精确积分，试讨论所需的Gauss积分的阶次（假定为常数）。解：插值函数中的多项式阶数为，微分算子中的导数的阶次是被积函数是非完全次项的最高次为次多项式，完全项的最高次为次多项式若为精确积分，需要高斯积分点故积分点数目为若为减缩积分，需要高斯积分点故积分点数目为6 求图1所示单元的节点等效荷载；图1解：在上则 7如图2所示12节点正方形单元，求其Jacobi行列式;图2解：求形函数：（1） 构造角节点形函数：（2） 构造边节点的形函数：，（3） 修正角节点形函数：8试构造如图3所示的6结点斜三棱柱体等参单元的插值函数，并证明其合理性。图3解：对图中6结点斜三棱柱体进行等参变换9空间八结点等参数单元各边与坐标轴x,y,z平行，在y方向作用有线性变化体力，若用高斯积分法分析结点荷载的精确值，试求所需要的最少积分点数。解：空间八结点等参单元中形函数的阶次为，且方向作用有阶次为的线性变化体力，由于单元各边与坐标轴x,y,z平行，故为常数，故被积函数的阶次为精确积分所需要高斯积分点，积分点数为第七章1 采用矩形薄板单元计算薄壳问题时，其单刚方程有何特点？解：采用矩形薄板单元计算薄壳时，为了简单计算，平板的面内变形与弯曲变形可认为是互不影响的，即板内变形和受力可看成是平面应力和平板弯曲两状态的迭加，结点未知数为单刚方程中，其中，特点：（1）基于（克希霍夫假设）中面无伸缩假设，可知与无关。（2）由于平行于中面的各层相互不挤压，不拉伸，沿方向不会引起翘曲，故与无关（3）对结点力不起作用，但为了计算不共面的相邻单元的弯扭应力，必须考虑。（4）和对应的刚度系数设定为零。2设薄板矩形单元，节点的位移未知数为：若位移模式取试判断该位移模式是否收敛？ 解：位移函数：完备性：挠度位移曲线中代表薄板的刚体位移，其中代表薄板在方向的移动，和分别代表薄板单元绕轴和轴的刚体转动。代表薄板弯曲的常应变（常曲率和常扭率）则挠度位移曲线满足完备性要求连续性：当x=常数(或y=常数)的边界上，挠度位移曲线是三次变化的曲线，可由两端节点的挠度值和转角值可唯一确定，故在单元交界面上是连续的。和也分别是和的三次曲线，由两端节点的转角和扭率可唯一确定，则也可唯一确定，故转角位移函数连续。所以，该位移模式收敛。3 论证矩形4节点12自由度薄板单元是完备的非协调单元。解：由于薄板弯曲变形时，可由中面的挠度表示，故4节点12自由度薄板单元的位移函数完备性：挠度位移曲线中代表薄板的刚体位移，其中代表薄板在方向的移动，和分别代表薄板单元绕轴和轴的刚体转动。代表薄板弯曲的常应变（常曲率和常扭率）则挠度位移曲线满足完备性要求连续性：当x=常数或y=常数的边界上，挠度位移曲线是三次变化的曲线，可由两端节点的挠度值和转角值唯一确定，故在单元交界面上是连续的。和分别是和的三次曲线，但是由两端节点的转角不能唯一确定，故转角位移函数不连续。所以，薄板矩形单元是非协调单元。4. 四边固定的正方形薄板，边长为4m，板厚为0.1m，弹性模量E为常量，。在板中心还联结有4根弹性杆件支承，杆长均为4m，所有杆与地面的夹角都为45o,弹性模量也为E，截面积为0.04m2。图4(1)为结构示意图，其中A-A,B-B 剖面见图4(2)，当板上受均匀分布荷载q0作用时，试求单元中点1的挠度（见图4(1)） 图4(1) 图4 (2) A-A 、B-B剖面解：分析单元，（1） 等效节点荷载:（2） 引入边界条件后（3） 考虑弹性支承杆的刚度（4） 在单刚方程中引入边界条件，并叠加弹性支承杆的刚度代入数据得-