高等有限元课后题答案.pdfVIP

2 弹性力学问题的有限单元法 思 考 题 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] 中的参数 bi 、 b j 、 bm 、 ci 、 c j 、 c m 由坐标变量 x、y 之差 确定。当单元的坐标差确定之后，这些参数与坐标变量 x、y 无关， 因此[B] 为常量阵。 当单元的结点位移 {a} 确定后， 由[B] 转换求得的单 元应变都是常量，也就是说在荷载作用下单元中各点具有统一的 x 、 y 、 xy 值。因此三结点三角形单元称为常应变单元。 (2)如果在每边中点增加一个结点，单元内的应力为线性分布。 习 题 2.1 试证明 x、y 与面积坐标的关系 证明 :设 P 点坐标为（ x,y ） 1 x y 1 A pij 1 x i y i 2 1 x j y j 1 x i y j y i x x j y x i y x j