高等有限元课后题答案分析.doc

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]中的参数由坐标变量x、y之差确定。当单元的坐标差确定之后，这些参数与坐标变量x、y无关，因此[B]为常量阵。当单元的结点位移{a}确定后，由[B]转换求得的单元应变都是常量，也就是说在荷载作用下单元中各点具有统一的值。因此三结点三角形单元称为常应变单元。 (2)如果在每边中点增加一个结点，单元内的应力为线性分布。 什么是面积坐标？如何计算三角形内某点的面积坐标？ 答：(1)如（a）图所示，三角形内任一点P（x，y），将P与三角形三个顶点i，j，m连成3个三角形。令为i点所对应三角形pjm的面积，为j点所对应的三角形pmi的面积，为m点所对应的三角形pij的面积，面积坐标定义为：= /A（i，j，m），其中A为三角形ijm的面积，点p（x，y）用面积坐标可以写为P（，，），且++ =1。 （2）求某点面积坐标除用定义外，还可用如图（b）所示的方法，即三角形内某点的面积坐标可通过同底三角形的高度比来计算。如图（b）中的= /。 （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 对于非协调的薄板单元如何进行分片检验？ 答：当赋予单元片各个结点以与常应变状态相应的位移值和载荷值时，校验是否满足，如能满足则认为通过分片检验。 在平面壳单元中如何判别共面点？可用什么方法进行处理？ 答：（1）在平面壳体单元中，如果某一点的各个单元面法向不同，经局部坐标转化到整体坐标后，该点的总体位移有6 个，若方向相同，常称此点为共面点