计算力学复习题答案

计算力学试题答案

1. 有限单元法和经典Ritz 法的主要区别是什么？

答：经典Ritz 法是在整个区域内假设未知函数，适用于边界几何形状简单的情形；有限单元法是将整个区域离散，分散成若干个单元，在单元上假设未知函数。有限单元法是单元一级的Ritz 法。 2、单元刚度矩阵和整体刚度矩阵各有什么特征？刚度矩阵[K ]奇异有何物理意义？在 求解问题时如何消除奇异性？

答：单元刚度矩阵的特征：⑴对称性⑵奇异性⑶主元恒正⑷平面图形相似、弹性矩阵D 、厚度t 相同的单元，e K 相同⑸e K 的分块子矩阵按结点号排列，每一子矩阵代表一个结点，占两行两列，其位置与结点位置对应。

整体刚度矩阵的特征：⑴对称性⑵奇异性⑶主元恒正⑷稀疏性⑸非零元素呈带状分布。

[]K 的物理意义是任意给定结构的结点位移所得到的结构结点力总体上满足力和力矩的平衡。 为消除[]K 的奇异性，需要引入边界条件，至少需给出能限制刚体位移的约束条件。

3. 列式说明乘大数法引入给定位移边界条件的原理？

答：设：j j a a =，则将 jj jj k k α=

j jj j P k a α=

即：

修改后的第j 个方程为

112222j j jj j j n n jj j k a k a k a k a k a αα++++

+=

由于

得 jj j jj j k a k a αα≈ 所以 j j a a ≈

对于多个给定位移()12,,

,l j c c c =时，则按序将每个给定位移都作上述修正，得到全部进行修正

后的K 和P ，然后解方程即可得到包括给定位移在内的全部结点位移值。 4. 何为等参数单元？为什么要引入等参数单元？

答：等参变换是对单元的几何形状和单元内的场函数采用相同数目的结点参数及相同的插值函数进行变换，采用等参变换的单元称之为等参数单元。

借助于等参数单元可以对于一般的任意几何形状的工程问题和物理问题方便地进行有限元离散，其优点有：对单元形状的适应性强；单元特性矩阵的积分求解方便（积分限标准化）；便于编制通用化程序。

5、对于平面4节点（线性）和8节点（二次）矩形单元，为了得到精确的刚度矩阵，

需要多少个Gauss 积分点？说明理由。

1112

1121121222222212

2212222222j n j n

j j jj j n j jj j n n nj n n n n k k k k a P k k k k a P k k k k a k a k k k k a P αα⎡⎤⎧⎫⎧⎫⎢⎥⎪⎪⎪⎪⎢⎥⎪⎪⎪⎪⎢⎥⎪⎪⎪⎪⎪⎪⎪⎪

=⎢

⎥⎨⎬⎨⎬

⎢⎥⎪⎪⎪⎪

⎢⎥⎪⎪⎪⎪⎢

⎥⎪⎪⎪⎪⎪⎪⎪⎪⎢⎥⎩⎭⎩⎭

⎣

⎦15

10α≈0 ()

ij

jj

k i j k α≈≠ ()

jj ij k k i j α>>≠

答：对于平面4节点（线性）矩形单元：

(,)i N ξη∝1,,,ξηξη

T B DB 221,,,,,ξηξηξη∝ =J 常数 所以2m =

因而积分点数为：22⨯矩阵

对于平面8节点（二次）矩形单元：

(,)i N ξη∝22221,,,,,,ξηξηξηξη

T B DB 221341,,,,,,ξηξηξηη∝ =J 常数 所以4m = 因而积分点数为：33⨯矩阵

⑴矩形、正方形、平行四边形=J 常数

1、有限单元法的解题步骤如何？它与经典Ritz 法有何区别？

答：⑴划分单元，输入结点和单元信息； ⑵单元分析：e

e

N K P 、、

⑶整体分析：

引入位移边界条件得到：=Ka P

⑷求解方程得到解a

⑸对位移a 结果进行有关整理、计算单元或结点的应力、应变

2、总刚度矩阵[K]的任一元素k ij 的物理意义是什么？如何解释总刚度矩阵的奇异性和

带状稀疏性？

答：K 中元素的ij K 物理意义：当结构的第j 个结点位移方向上发生单位位移，而其它结点位移方向上位移为零时，需在第i 个结点位移方向上施加的结点力大小。

奇异性：K

=0,力学意义是对任意给定的结点位移所得到的结构结点力总体上是满足力和力矩的平

衡。反之，给定任意满足力和力矩平衡的结点载荷P ，由于K 的奇异性却不能解得结构的位移a ,因

而结构仍可能发生任意的刚体位移。为消除[]K 的奇异性，结构至少需给出能限制刚体位移的约束条件。

带状稀疏性：由于连续体离散为有限个单元体时，每个结点的相关单元只是围绕在该结点周围为数甚少的几个，一个结点通过相关单元与之发生关系的相关结点也只是它周围的少数几个，因此虽然总体单元数和结点数很多，结构刚度矩阵的阶数很高，但刚度系数中非零系数却很少，即为总刚度矩阵的稀疏性。另外，只要结点编号是合理的，这些稀疏的非零元素将集中在以主对角线为中心的一条带状区域内，即为总刚度矩阵的带状分布特性。

1 1.5

2

m n +≥=141

2.522

m n ++≥==1

,

e

n e e e e T ==∑K G K G 1

e

n e e e T ==∑P G P

3、以3节点三角形单元为例证明插值函数特性

11

i =∑=n

i

N

，n 为节点数。

答： 图形见课本P105图3.6 由面积坐标：

插值函数：i i N L =

()i j m P L ，L ，L

所以

4、什么是等参单元？等参单元的收敛性如何？

答：等参变换是对单元的几何形状和单元内的场函数采用相同数目的结点参数及相同的插值函数进行变换，采用等参变换的单元称之为等参元。

等参单元满足收敛性需满足两个条件：即单元必须是协调的和完备的。完备性条件：要求插值函数中包含完全的线性项（包含常数项和一次项）。协调性条件：单元边界上位移连续，相邻单元边界具有相同的结点，每一单元沿边界的坐标和未知函数采用相同的插值函数。

5、对于空间8节点（线性）和20节点（二次）六面体单元，为了得到精确的刚度矩

阵，需要多少个Gauss 积分点？说明理由。 答：

对于空间8节点（线性）六面体单元：

(,)i N ξη∝1,,,,,,,x y z xy yz zx xyz

T B DB 221,,,,,x xy xz x y ∝ =J 常数 所以2m = 因而积分点数为：222⨯⨯矩阵

对于空间20节点（二次）六面体单元：

(,)i N ξη∝2223332222221,,,,,,,,,,,,,,,,,,,x y z x y z x y z xy yz zx x y xy x z xz y z yz xyz

T B DB 41,,,,,x xy xz x ∝ =J 常数

所以4m =

因而积分点数为：333⨯⨯矩阵

i

i A L A =

(i,j,m)i =n

j i j m i m i i 1+1A A A A A A A N A A A =+=++==∑(i,j,m)

i =1

1.5

2m n +≥=1 2.5

2

m n +≥=