计算力学复习题答案

计算力学试题答案1.有限单元法和经典Ritz法的主要区别是什么？答：经典Ritz法是在整个区域内假设未知函数，适用于边界几何形状简单的情形；有限单元法是将整个区域离散，分散成若干个单元，在单元上假设未知函数。有限单元法是单元一级的Ritz法。2、单元刚度矩阵和整体刚度矩阵各有什么特征？刚度矩阵[K]奇异有何物理意义？在求解问题时如何消除奇异性？答：单元刚度矩阵的特征：⑴对称性⑵奇异性⑶主元恒正⑷平面图形相似、弹性矩阵D、厚度t相同的单元，eK相同⑸eK的分块子矩阵按结点号排列，每一子矩阵代表一个结点，占两行两列，其位置与结点位置对应。整体刚度矩阵的特征：⑴对称性⑵奇异性⑶主元恒正⑷稀疏性⑸非零元素呈带状分布。K的物理意义是任意给定结构的结点位移所得到的结构结点力总体上满足力和力矩的平衡。为消除K的奇异性，需要引入边界条件，至少需给出能限制刚体位移的约束条件。3.列式说明乘大数法引入给定位移边界条件的原理？答：设：jjaa，则将jjjjkkjjjjPka即：修改后的第j个方程为112222jjjjjjnnjjjkakakakaka由于得jjjjjjkaka所以jjaa对于多个给定位移12,,,ljccc时，则按序将每个给定位移都作上述修正，得到全部进行修正后的K和P，然后解方程即可得到包括给定位移在内的全部结点位移值。4.何为等参数单元？为什么要引入等参数单元？答：等参变换是对单元的几何形状和单元内的场函数采用相同数目的结点参数及相同的插值函数进行变换，采用等参变换的单元称之为等参数单元。借助于等参数单元可以对于一般的任意几何形状的工程问题和物理问题方便地进行有限元离散，其优点有：对单元形状的适应性强；单元特性矩阵的积分求解方便（积分限标准化）；便于编制通用化程序。5、对于平面4节点（线性）和8节点（二次）矩形单元，为了得到精确的刚度矩阵，需要多少个Gauss积分点？说明理由。111211211212222222122212222222jnjnjjjjjnjjjjnnnjnnnnkkkkaPkkkkaPkkkkakakkkkaP15100()ijjjkijk()jjijkkij答：对于平面4节点（线性）矩形单元：(,)iN1,,,BDB221,,,,,J常数所以2m因而积分点数为：22矩阵对于平面8节点（二次）矩形单元：(,)iN22221,,,,,,BDB221341,,,,,,J常数所以4m因而积分点数为：33矩阵⑴矩形、正方形、平行四边形J常数1、有限单元法的解题步骤如何？它与经典Ritz法有何区别？答：⑴划分单元，输入结点和单元信息；⑵单元分析：eeNKP、、⑶整体分析：引入位移边界条件得到：KaP⑷求解方程得到解a⑸对位移a结果进行有关整理、计算单元或结点的应力、应变2、总刚度矩阵[K]的任一元素kij的物理意义是什么？如何解释总刚度矩阵的奇异性和带状稀疏性？答：K中元素的ijK物理意义：当结构的第j个结点位移方向上发生单位位移，而其它结点位移方向上位移为零时，需在第i个结点位移方向上施加的结点力大小。奇异性：K=0,力学意义是对任意给定的结点位移所得到的结构结点力总体上是满足力和力矩的平衡。反之，给定任意满足力和力矩平衡的结点载荷P，由于K的奇异性却不能解得结构的位移a,因而结构仍可能发生任意的刚体位移。为消除K的奇异性，结构至少需给出能限制刚体位移的约束条件。带状稀疏性：由于连续体离散为有限个单元体时，每个结点的相关单元只是围绕在该结点周围为数甚少的几个，一个结点通过相关单元与之发生关系的相关结点也只是它周围的少数几个，因此虽然总体单元数和结点数很多，结构刚度矩阵的阶数很高，但刚度系数中非零系数却很少，即为总刚度矩阵的稀疏性。另外，只要结点编号是合理的，这些稀疏的非零元素将集中在以主对角线为中心的一条带状区域内，即为总刚度矩阵的带状分布特性。11.52mn1412.522mn1,eneeeeKGKG1eneeePGP3、以3节点三角形单元为例证明插值函数特性11iniN，n为节点数。答：图形见课本P105图3.6由面积坐标：插值函数：iiNLijmPL，L，L所以4、什么是等参单元？等参单元的收敛性如何？答：等参变换是对单元的几何形状和单元内的场函数采用相同数目的结点参数及相同的插值函数进行变换，采用等参变换的单元称之为等参元。等参单元满足收敛性需满足两个条件：即单元必须是协调的和完备的。完备性条件：要求插值函数中包含完全的线性项（包含常数项和一次项）。协调性条件：单元边界上位移连续，相邻单元边界具有相同的结点，每一单元沿边界的坐标和未知函数采用相同的插值函数。5、对于空间8节点（线性）和20节点（二次）六面体单元，为了得到精确的刚度矩阵，需要多少个Gauss积分点？说明理由。答：对于空间8节点（线性）六面体单元：(,)iN1,,,,,,,xyzxyyzzxxyzBDB221,,,,,xxyxzxyJ常数所以2m因而积分点数为：222矩阵对于空间20节点（二次）六面体单元：(,)iN2223332222221,,,,,,,,,,,,,,,,,,,xyzxyzxyzxyyzzxxyxyxzxzyzyzxyzBDB41,,,,,xxyxzxJ常数所以4m因而积分点数为：333矩阵iiALA(i,j,m)injijmimii1+1AAAAAAANAAA(i,j,m)i11.52mn12.52mn1、为什么说3节点三角形单元是常应变单元？答：常应变单元指的是在一个单元内的应变为常数，有限元中的常应变单元指的是线性三角形单元，线性三角形单元的位移场为线性的，应变为位移的一阶导数，故为常数，因此称为常应变单元。2.以平面4节点双线性单元为例，说明形函数的两个重要特性。答：图形见课件2.5矩形单元插值函数（形函数）的性质进而验证插值函数的性质：3、何为等参变换？等参元有那些优点？答：等参变换是对单元的几何形状和单元内的场函数采用相同数目的结点参数及相同的插值函数进行变换，采用等参变换的单元称之为等参元。借助于等参元可以对于一般的任意几何形状的工程问题和物理问题方便地进行有限元离散，其优点有：对单元形状的适应性强；单元特性矩阵的积分求解方便（积分限标准化）；便于编制通用化程序。二、如图所示平面问题有限元网格，每个单元4个节点，每个节点2个自由度，1.给出适当的节点编号，使总的系数矩阵的半带宽最小，并计算半带宽的值；2.在您的节点编号下，图中节点A的主对角线上的元素在总系数矩阵中的行号和列号如何？3.哪几个单元对节点A的主对角线上的系数有非零贡献？1、答：沿短边回头编号，存储量最小。带宽的计算：(1)D相邻结点号码的最大差值自由度2、答：由得节点A的主对角线上的元素11,11k、12,12k在总系数矩阵中的行号为11和12,列号为11,12。11114N21114N31114N41114N1()(,)0()ijjijijNij411iiN11,1111,126612,1112,12kkKkk3、答：2、3、4单元对A的主对角线上的系数有非零贡献。注意：杆件单元在每个节点上有1个自由度，带宽不用乘以2。三、图示6节点三角形单元的142边作用有均布侧压力q，单元厚度为t，求单元的等效节点载荷。见课件2.5.3—5采用面积坐标时，单元矩阵的计算中：（2）均布侧压力q和（4）x方向三角形分布力中的形函数。四、见课件3.3.31234567891011121314151234567A五、答：⑴观察插值函数，包含一次完全多项式，满足完备性要求，但其含有三次项37a与38a，此点与八点矩形单元四个边上均仅有三个点，最多满足二次项矛盾，不满足协调性，因而不满足收敛性。⑵观察差值函数，包含一次完全多项式，多项式最高次为2，明显既满足完备性和协调性，故收敛。第三页PDF：六、考虑等截面轴力杆单元，题图中分别示出2节点和3节点单元体，1.写出它们的位移插值函数；2.推导这两种单元体的刚度矩阵；3.对3节点单元体用静力凝聚法消去中间节点自由度，建立单元体刚度矩阵表达式。解：⑴图a，令1212212()2()xxxxxxxxL，则有121,1,11。故有，(1)211121()()(1)2Nl，(1)122211()()(1)2Nl。图b，令1212312()2()xxxxxxxxL，则有1231,0,1,11。故有，(2)23111213()()1()()(1)()()2Nl(2)213222123()()()()1()()Nl(2)12313132()()1()()(1)()()2Nl⑵图a：有应变矩阵12122dNdNdNdNddNBLNdxdxdddxdL。00011211TTLLLeTdNdNEAdNdNEAKBDBdxEAdxddxdxLddL。图b：有应变矩阵1111112dNdNdNdNdNdNddNBLNdxdxdxddddxdL。00014212216831814TTLLLeTdNdNEAdNdNEAKBDBdxEAdxddxdxLddL。⑶1122331421216831814uFEAuFLuF，有12322213332168(28)/16LLuuuFuFuuEAEA代入12313142LuuuFEA与12333814LuuuFEA，消去2u，得：121332171084132102FFuEAuLFF。二、用最小势能原理推导平面弹性问题的有限元方程。解：利用最小势能（位能）原理推导平面弹性问题有限元方程，见课本P61。第二章课件也有。三、如图所示的平面内部三角形单元网格每节点2个自由度，1.用图中所给节点编号计算总刚矩阵的半带宽；2.对节点重新编号，使结构总刚矩阵的半带宽最小，并说明此时中心节点主对角线上的元素在总刚矩阵中的行号和列号？3.是否所有单元对中心节点主对角线上的元素都有非零贡献？4.两种编号下，结构总刚矩阵中的非零元素是否相等？为什么？解：⑴半带宽（相邻结点号码的最大差值+1）自由度(91)216。⑵中心节点元素编号为5，其余元素编号顺时针依次编写为1-4，6-9。此时中心节点对应主对角线上的元素9,9K、10,10K在总刚度矩阵中的行号、列号分别为第9行和第9列，第10行和第10列⑶是。由于8个三角型单元均有一个节点是中心节点主对角线上元素，有123455555555555,,,,KKKKK678555555,,,KKK均不为零，故所有单元对主对角线上元素都有非零贡献。⑷两种编号方式，非零元素相等。编号的合理化只是将非零元素的位置集中在以主对角线为中心的一条带状区域内，但并不改变非零元素的个数。四、3节点三角形单元的jm边作用有分布侧压力，如图示，单元厚度为t，求单元的等效节点载荷。解：如图梯形分布力可以分解成为一