1、计算力学试题答案1. 有限单元法和经典 Ritz 法的主要区别是什么？答：经典 Ritz法是在整个区域内假设未知函数，适用于边界几何形状简单的情形；有限单元法是将整个区域离散，分散成若干个单元，在单元上假设未知函数。有限单元法是单元一级的 Ritz法。2、单元刚度矩阵和整体刚度矩阵各有什么特征？刚度矩阵K奇异有何物理意义？在求解问题时如何消除奇异性？答：单元刚度矩阵的特征：对称性奇异性主元恒正平面图形相似、弹性矩阵 D、厚度 t相同的单元， 相同 的分块子矩阵按结点号排列，每一子矩阵代表一个结点，占两行两列，eeK其位置与结点位置对应。整体刚度矩阵的特征：对称性奇异性主元恒正稀疏性非零元素呈带
2、状分布。的物理意义是任意给定结构的结点位移所得到的结构结点力总体上满足力和力矩的平衡。K为消除 的奇异性，需要引入边界条件，至少需给出能限制刚体位移的约束条件。3. 列式说明乘大数法引入给定位移边界条件的原理？答：设： ，则将 jjajjkPa即：修改后的第 个方程为j12 2jj jjnjkakaka 由于得 jj所以 ja对于多个给定位移 时，则按序将每个给定位移都作上述修正，得到全部进行修正12,lc后的 K和 P，然后解方程即可得到包括给定位移在内的全部结点位移值。4. 何为等参数单元？为什么要引入等参数单元？答：等参变换是对单元的几何形状和单元内的场函数采用相同数目的结点参数及相同的
3、插值函数进行变换，采用等参变换的单元称之为等参数单元。借助于等参数单元可以对于一般的任意几何形状的工程问题和物理问题方便地进行有限元离散，其优点有：对单元形状的适应性强；单元特性矩阵的积分求解方便（积分限标准化） ；便于编制通用化程序。5、对于平面4节点（线性）和8节点（二次）矩形单元，为了得到精确的刚度矩阵，1211212 21222 2jnjjjjnjjnnj nk PkakP 1500 ()iji ()jij需要多少个 Gauss 积分点？说明理由。答：对于平面 4节点（线性）矩形单元：(,)iN1,BD2,J常 数所以 m因而积分点数为： 矩阵2对于平面 8节点（二次）矩形单元：(,)
4、iN21,BD2134J常 数所以 4m因而积分点数为： 矩阵3矩形、正方形、平行四边形 J常 数1、有限单元法的解题步骤如何？它与经典 Ritz 法有何区别？答：划分单元，输入结点和单元信息；单元分析： eNKP、 、整体分析： 引入位移边界条件得到： KaP求解方程得到解 a对位移 结果进行有关整理、计算单元或结点的应力、应变2、总刚度矩阵K的任一元素kij 的物理意义是什么？如何解释总刚度矩阵的奇异性和带状稀疏性？答：K 中元素的 物理意义：当结构的第 个结点位移方向上发生单位位移，而其它结点位ij j移方向上位移为零时，需在第 个结点位移方向上施加的结点力大小。奇异性： =0,力学意义
5、是对任意给定的结点位移所得到的结构结点力总体上是满足力和力矩的平衡。反之，给定任意满足力和力矩平衡的结点载荷 P，由于 K的奇异性却不能解得结构的位移 ,因a而结构仍可能发生任意的刚体位移。为消除 的奇异性，结构至少需给出能限制刚体位移的约束条件。带状稀疏性：由于连续体离散为有限个单元体时，每个结点的相关单元只是围绕在该结点周围为数甚少的几个，一个结点通过相关单元与之发生关系的相关结点也只是它周围的少数几个，因此虽然总体单元数和结点数很多，结构刚度矩阵的阶数很高，但刚度系数中非零系数却很少，即为总刚度矩阵的稀疏性。另外，只要结点编号是合理的，这些稀疏的非零元素将集中在以主对角线为中心的一条带状
6、区域内，即为总刚度矩阵的带状分布特性。1.5n2.5n1,eneG1eneP3、以3节点三角形单元为例证明插值函数特性 ，n 为节点数。1iniN答： 图形见课本 P105图 3.6由面积坐标： 插值函数： iiNLijmPL， ，所以4、什么是等参单元？等参单元的收敛性如何？答：等参变换是对单元的几何形状和单元内的场函数采用相同数目的结点参数及相同的插值函数进行变换，采用等参变换的单元称之为等参元。等参单元满足收敛性需满足两个条件：即单元必须是协调的和完备的。完备性条件：要求插值函数中包含完全的线性项（包含常数项和一次项） 。协调性条件：单元边界上位移连续，相邻单元边界具有相同的结点，每一单
7、元沿边界的坐标和未知函数采用相同的插值函数。5、对于空间8节点（线性）和20节点（二次）六面体单元，为了得到精确的刚度矩阵，需要多少个 Gauss 积分点？说明理由。答：对于空间 8节点（线性）六面体单元：(,)iN1,xyzxyzBD2J常 数所以 2m因而积分点数为： 矩阵对于空间 20节点（二次）六面体单元：(,)iN232221,xyzxyzxyzxyzxBD4J常 数所以 4m因而积分点数为： 矩阵3iiA(i,j)njijmimii1 +1AA 1.5n12.5n1、为什么说 3 节点三角形单元是常应变单元？答：常应变单元指的是在一个单元内的应变为常数，有限元中的常应变单元指的是线
8、性三角形单元，线性三角形单元的位移场为线性的，应变为位移的一阶导数，故为常数，因此称为常应变单元。2.以平面 4 节点双线性单元为例，说明形函数的两个重要特性。答：图形见课件 2.5矩形单元插值函数（形函数）的性质进而验证插值函数的性质：3、何为等参变换？等参元有那些优点？答：等参变换是对单元的几何形状和单元内的场函数采用相同数目的结点参数及相同的插值函数进行变换，采用等参变换的单元称之为等参元。借助于等参元可以对于一般的任意几何形状的工程问题和物理问题方便地进行有限元离散，其优点有：对单元形状的适应性强；单元特性矩阵的积分求解方便（积分限标准化） ；便于编制通用化程序。2、如图所示平面问题有
9、限元网格，每个单元4 个节点，每个节点2 个自由度，1. 给出适当的节点编号，使总的系数矩阵的半带宽最小，并计算半带宽的值；2. 在您的节点编号下，图中节点A 的主对角线上的元素在总系数矩阵中的行号和列号如何？3. 哪几个单元对节点 A 的主对角线上的系数有非零贡献？1、答：沿短边回头编号，存储量最小。带宽的计算： ( 1)D相 邻 结 点 号 码 的 最 大 差 值 自 由 度2、答：由得节点 A的主对角线上的元素 、 在总系数矩阵中的行号为 11和 12,列号为 11,12。1,k2,11N24314N1 ()(,)0ijijijN41ii1,262kK3、答：2、3、4 单元对 A的主对
10、角线上的系数有非零贡献。注意：杆件单元在每个节点上有 1个自由度，带宽不用乘以 2。三、图示6 节点三角形单元的142 边作用有均布侧压力q，单元厚度为t，求单元的等效节点载荷。见课件 2.5.35采用面积坐标时，单元矩阵的计算中：（2）均布侧压力 q和（4）x 方向三角形分布力中的形函数。四、见课件 3.3.312345678910112131415234567A五、答：观察插值函数，包含一次完全多项式，满足完备性要求，但其含有三次项 与 ，此37a8点与八点矩形单元四个边上均仅有三个点，最多满足二次项矛盾，不满足协调性，因而不满足收敛性。观察差值函数，包含一次完全多项式，多项式最高次为 2
11、，明显既满足完备性和协调性，故收敛。第三页 PDF：六、考虑等截面轴力杆单元，题图中分别示出 2 节点和 3 节点单元体，1. 写出它们的位移插值函数；2. 推导这两种单元体的刚度矩阵；3. 对 3 节点单元体用静力凝聚法消去中间节点自由度，建立单元体刚度矩阵表达式。解：图 a， 令 ，则有 。1212()()xxL12,1故有， ， 。(1)211()Nl (1)122()Nl图 b，令 ，则有 。21231()xxL123,0,故有， (2)23111)()Nl(2) 232212)(l(2)313132)(1)Nl 图 a：有应变矩阵 。12dNdNdBLxxL。000 1T TL Le
12、T EAEAKDdxEA 图 b：有应变矩阵 。11112dNdNdNBxxxL。0001422683T TLLLeTdNEAdNEAKBDxEAxL ，有123142683uFEAL123221368()/6uFuuEA代入 与 ，消去 ，得：1214LuEA23384L2。12370321FEAL二、用最小势能原理推导平面弹性问题的有限元方程。解：利用最小势能（位能）原理推导平面弹性问题有限元方程，见课本 P61。第二章课件也有。三、如图所示的平面内部三角形单元网格每节点 2 个自由度，1. 用图中所给节点编号计算总刚矩阵的半带宽；2. 对节点重新编号，使结构总刚矩阵的半带宽最小，并说明此
13、时中心节点主对角线上的元素在总刚矩阵中的行号和列号？3. 是否所有单元对中心节点主对角线上的元素都有非零贡献？4. 两种编号下，结构总刚矩阵中的非零元素是否相等？为什么？解：半带宽 （相邻结点号码的最大差值+1） 自由度 。(91)26中心节点元素编号为 5，其余元素编号顺时针依次编写为 1-4，6-9。此时中心节点对应主对角线上的元素 、 在总刚度矩阵中的行号、列号分别为第 9行和第 9列，第 10行和第 10列9,K10,是。由于 8个三角型单元均有一个节点是中心节点主对角线上元素，有 1234555,K均不为零，故所有单元对主对角线上元素都有非零贡献。6755,两种编号方式，非零元素相等
14、。编号的合理化只是将非零元素的位置集中在以主对角线为中心的一条带状区域内，但并不改变非零元素的个数。四、3 节点三角形单元的jm 边作用有分布侧压力，如图示，单元厚度为 t，求单元的等效节点载荷。解：如图梯形分布力可以分解成为一个三角形分布力和一个均匀分布力，然后叠加。参加课本 P68页或者课件，易知：均匀分布部分， 10102TePqlt三角型分布部分： 12()3Telt两者相加， 21210006 Teltqq第四页 PDF五、块体棱柱单元的上下三角形表面 = 1，利用面积坐标L (i1, 2, 3) i 和等参坐标写出该线性单元的形函数。解： 11()()2LNL；22()() 33(
15、1)()2NL； ； 。411LN 55L 61()12NL六、对图示的矩形单元采用如下的插值函数1. 请分析该单元的协调性。2. 若将插值多项式写成解：如图所示的矩形单元若满足协调性，因为在 y方向有三个节点，x 方向有两个节点，故 y方向上，差值函数视 x为常数 y的次数至少为 3次，x 方向上，插值函数视 y为常数，x 的次数至少为 2次。对比给出的方程知，该单元若采用该插值函数满足协调性。插值函数改变后，增加了一个 ，而由图中矩形单元知道 x方向多项式最高次为 2次，显然不3x满足协调性。七、证明 4 节点平行四边形二维单元的雅可比矩阵是常数矩阵.对于二维平行四边形有，有 （7.1）
16、44 1312411 2 34(,)iiii xyNNxyxyJ等参变换下取， ，而在自然坐标下不妨设该平行四边形的四边的方程分别为：ii 从而知道：0abc0abc； ；1()(Nkabc2()(Nkabc3()(Nkabc。其中， ， ， ，4 1421ka144解得四个节点的坐标，一并代入方程 7.1得： 的四个元素均为常数，故有 4 节点平行四边形二J维单元的雅可比矩阵是常数矩阵。第五页 pdf该单元在结构中的位置使得总体节点编号分别为 19、20、30、31，回答：1. 在未引入边界条件以前，j 单元刚度矩阵K e 的系数将贡献给总体刚度矩阵K中的哪些行、列？2. 具体写出刚度矩阵K
17、 e 中的哪些元素对总体刚度矩阵K中的下列行和列有贡献，（1）59 行61 列；（2）38 行39 列；（3）59 行59 列；（4）37 行37 列。（题解：由节点位移矢量知该结构节点自由度为 2，设某一个节点的编号为 ，则该个节点将分别对i总体刚度矩阵的第 、 行与列，产生贡献。所以，总体编号为 19、20、30、31 的总体节点21ii将相应对总体刚度矩阵的 37、38 行和列，39、40 行和列，59、60 行和列，60、61 行和列产生贡献。对第 59行 61列产生贡献的是 ；对第 38行 39列产生贡献的是 ；对第 59行与 5935K27K列产生贡献的是 ；对 37行和 37列产生贡献的是 .3 1三、3 节点三角形单元的ij 边作用有三角形分布侧压力，单元厚度为t，求单元的等效节点载荷。 解：答案参见课本 P68页，或者课件。此处略。五、下面是3 节点轴力杆单元示意图及其刚度矩阵，试将刚度矩阵凝聚成只以u1 和u3表示。解： ，有12378166uFEAL12322136678(8)/7LLuFuuEA代入 与 ，消去 ，得：123178Lu233862。12337FEAL
