固体力学中的数值方法试题(答案).doc

有限单元法的解题步骤如何？有限单元法和经典Ritz 法的主要区别是什么？ 答：⑴ 结构离散：划分单元，结点编号和载荷移植； ⑵单元分析：单元体内假设位移模式，用节点的未知参数表征单元特性。 ⑶整体分析：将各单元关系式组合成包含节点未知数的整体方程组，引入位移边界条件得到：。 ⑷求解方程得到节点位移 ⑸对位移d结果进行有关整理、计算单元或结点的应力、应变 经典Ritz法是在整个区域内假设未知函数，适用于边界几何形状简单的情形；有限单元法是将整个区域离散，分散成若干个单元，在单元上假设未知函数。有限单元法是单元一级的Ritz法。 2、单元刚度矩阵和整体刚度矩阵各有什么特征？刚度矩阵[K]奇异有何物理意义？在 求解问题时如何消除奇异性？ 答：单元刚度矩阵的特征：⑴对称性（2）与单元位置无关（3）奇异性（4）主元恒正（4）平面图形相似、弹性矩阵D、厚度t相同的单元，相同⑸的分块子矩阵按结点号排列，每一子矩阵代表一个结点，占两行两列，其位置与结点位置对应。 整体刚度矩阵的特征：⑴对称性⑵奇异性⑶主元恒正⑷稀疏性⑸非零元素呈带状分布。 的物理意义是任意给定结构的结点位移所得到的结构结点力总体上满足力和力矩的平衡。 为消除的奇异性，需要引入边界条件，至少需给出能限制刚体位移的约束条件。 3. 列式说明乘大数法引入给定位移边界条件的原理？ 答：设：，则将 即： 修改后的第个方程为 由于 得 所以 对于多个给定位移时，则按序将每个给定位移都作上述修正，得到全部进行修正后的K和P，然后解方程即可得到包括给定位移在内的全部结点位移值。 何为等参数单元？为什么要引入等参数单元？ 答：等参变换是对单元的几何形状和单元内的场函数采用相同数目的结点参数及相同的插值函数进行变换，采用等参变换的单元称之为等参数单元。 借助于等参数单元可以对于一般的任意几何形状的工程问题和物理问题方便地进行有限元离散，其优点有：对单元形状的适应性强；单元特性矩阵的积分求解方便（积分限标准化）；便于编制通用化程序。 5、对于平面4节点（线性）和8节点（二次）矩形单元，为了得到精确的刚度矩阵， 需要多少个Gauss 积分点？说明理由。 答：对于平面4节点（线性）矩形单元： 所以 因而积分点数为：矩阵 对于平面8节点（二次）矩形单元： 所以 因而积分点数为：矩阵 ⑴矩形、正方形、平行四边形 6、总刚度矩阵[K]的任一元素kij 的物理意义是什么？如何解释总刚度矩阵的奇异性和 带状稀疏性？ 答：K 中元素的物理意义：当结构的第个结点位移方向上发生单位位移，而其它结点位移方向上位移为零时，需在第个结点位移方向上施加的结点力大小。 奇异性：=0,力学意义是对任意给定的结点位移所得到的结构结点力总体上是满足力和力矩的平衡。反之，给定任意满足力和力矩平衡的结点载荷P，由于K的奇异性却不能解得结构的位移,因而结构仍可能发生任意的刚体位移。为消除的奇异性，结构至少需给出能限制刚体位移的约束条件。 带状稀疏性：由于连续体离散为有限个单元体时，每个结点的相关单元只是围绕在该结点周围为数甚少的几个，一个结点通过相关单元与之发生关系的相关结点也只是它周围的少数几个，因此虽然总体单元数和结点数很多，结构刚度矩阵的阶数很高，但刚度系数中非零系数却很少，即为总刚度矩阵的稀疏性。另外，只要结点编号是合理的，这些稀疏的非零元素将集中在以主对角线为中心的一条带状区域内，即为总刚度矩阵的带状分布特性。 7、以3节点三角形单元为例证明插值函数特性，n 为节点数。 答： 图形见课本P105图3.6 由面积坐标： 插值函数： 所以 8、如图所示平面问题有限元网格，每个单元4 个节点，每个节点2 个自由度， 1. 给出适当的节点编号，使总的系数矩阵的半带宽最小，并计算半带宽的值； 2. 在您的节点编号下，图中节点A 的主对角线上的元素在总系数矩阵中的行号和列号如何？ 3. 哪几个单元对节点A 的主对角线上的系数有非零贡献？ 1、答：沿短边回头编号，存储量最小。 带宽的计算： 2、答：由 得节点A的主对角线上的元素、在总系数矩阵中的行号为11和12,列号为11,12。 3、答：2、3、4单元对A的主对角线上的系数有非零贡献。 注意：杆件单元在每个节点上有1个自由度，带宽不用乘以2。 答：⑴观察插值函数，包含一次完全多项式，满足完备性要求，但其含有三次项 与 ，此点与八点矩形单元四个边上均仅有三个点，最多满足二次项矛盾，不满足协调性，因而不满足收敛性。 ⑵观察差值函数，包含一次完全多项式，多项式最高次为2，明显既满足完备性和协调性，