东北大学成型力学博士考题

第一题：泛函可变化的函数称为可变函数，随可变函数而变的量，称为泛函。常应变单元B矩阵中的元素都是常量，单元中各点的应变分量也都是常量这种单元就是常应变单元。形状因子(,)(,,)kNxykijm是节点坐标的函数，只决定于单元的形状、节点配置和插值方式，所以把(,)kNxy称为形状函数或者形状因子。虚位移场：在uS上满足0iv，在U内满足,,1()2ijijjiuu的位移场iu成为虚位移场。[D]矩阵：与材料有关的弹性矩阵，对平面应力状态2101011002ED[B]矩阵：单元应变矩阵，它可以写成分块形式ijmBBBB，而子矩阵1(,,)2iiiibBcijmAcbi由于,,,,,ijmijmbbbccc都是常量，所以矩阵[B]中的元素都是常量，因而单元中各点的应变速率分量也都是常量，故通常称为这种单元为常应变单元。等参单元：具有相同的形式，用同样个数的相应结点值做参数，并具有完全相同的形状函数作为这些结点值前面的系数，具有这样特征的单元称为等参单元。第二题：采用能量法解平辊轧制矩形件：1、采用hill速度场，并设xcb则得速度场如下：xxx''xxyxxxxx'xzxxU1v=ChbbbU1v=y=vyChbbbhv=vzh式中000UvhbC，为秒流量；()xxbbx为轧件宽度的一半。由于'[]xyybxxvbv，则满足稳定轧制过程条件。2、由于'',xxzxxyxxzRybvvhvvb则该速度场满足表面的边界条件。而'''',,xxxxxxyxzxxxxxbhvvvbhbhbh由于0xyz满足体积不变条件。3、按塑性变分原理，泛函dfsNNN应取极小值。内部变形功率014()32lsdNUpxdx摩擦功率fN04()3lsUmQxdxsN为出入口的剪切功率。式中：''''()(,,,)()(,,,)xxxxxnpxbbbxQxbbxx则''''''0(,,,)(,,,,)lxxxxxxnbbbxFbbbxxdx使泛函求极值的函数xb可由如下欧拉方程确定：2'''0xxxFdFdFbdxbdxb而中性面坐标由0nx确定。4、用里兹法求的近似极值系数()xbx设202021101223322()()()xbabaaabxbaxxllll把上式带入泛函中，有12(,,)naax可通过12000naax解出12naax，求出()xbx。5、由于上述方程组是非线性的，使用Newton-Laphson法线性化，令三个方程为：()0(,1,2,3)ijfxij线性化后有：00,()0jjiijijjxxffxfdx从而求出123dxdxdx然后迭代计算，直到得出符合要求的解。6、求出()xbx之后，便可求出功和功率等。第三题：刚塑性有限元的解法有拉格朗日乘子法、体积可压缩法和罚函数法。拉格朗日乘子法能量泛函表达式：1、把刚塑性材料不完全变分原理写成矩阵形式：23TTTsvspvdvvpdscdv2、将连续体离散化成共有m个单元的n个节点，则对第e个单元，单元能量泛函：1/22([][])[N]{}[]{}3TTTeepeeeeeTTsvvsuBBudvupdsuBCdv其中：123121(,,......,,......)menmeuuu3、求泛函的极值条件：1()0eemijeijvv由于iv和j是任意的独立变量，所以有1()0emieivv（i=1,2,3,4……3n）1()0emjej(j=1,2,3,4……m)由以上m+3n个方程联立求出iv和；先求{}eev和ee1/2[][]2[N]{}[]{}{}3([][])TTTeepeTeeeeTseeeeTvvsBBvdvvpdsvBCdvvvBBv{}[]{}eeTeeevCBvdv然后对所有单元组合，使得到的未知数个数相同的联立方程组。但此方程组关于{v}是非线性的，所以采用摄动法进行线性化。线性化之后，得到线性方程组，便可以得出{v}的收敛解。拉格朗日法与罚函数相比，由于前者引入了拉格朗日乘子，若单元数目很大，这个附加的未知量将会增加联立方程数目或者增加系数矩阵的带宽，从而增加了计算时间，体积可压缩法和罚函数法虽然避免了附加未知量，但若初始速度场设定不好，会导致vm非常大，难以得出正确解。四弹性有限元单元总势能表达式：弹塑性有限元单元总势能表达式：11/21([][][])[N]{}2eeTeeTeTepevsuBDBudvuqds[D]为单元弹性矩阵[]epD为弹塑性矩阵。异同点：当物体产生塑性变形后，塑性变形区内的几何方程和平衡方程与弹性变形时相同，但是，其应力应变关系，则由线性变成非线性的，应力和应变也不在一一对应。弹塑性有限元不能一步得出结果，可以使用变刚度法或初载荷法求解。[]epD只与加载前的应力水平有关，而与增量无关，又由于位移插值函数、位移－应变的几何关系与弹性变形时相同，所以对于每次加载，解弹塑性问题都可用弹性有限元法完全相同的计算格式相同。四．给出弹、塑性有限元单元总势能泛函的矩阵表达式，概述两种解法的主要异同。弹性有限元：e12TTepTTeeeevsuBDBudVuNPdS弹塑性：e12TTpeTeeeeeevuBDBudVuF相同点：计算格式相同，几何方程、平衡方程相同不同点：1.弹性有限元应力应变关系是线性的，弹塑性有限元是非线性的，应力和应变不算一一对应。2.刚度矩阵不同。弹性：eTeevKBDBdV弹塑性：eTeeepvKBDBdV3.使用范围不同。弹性有限元适用于未屈服的弹性状态弹塑性有限元适用于塑性变形状态。第五题：1、616264KKK为0矩阵，[K]矩阵时一个稀疏阵，这种稀疏阵的形成是由于整体刚度矩阵是把全部单元刚度矩阵按节点编号叠加所致。对于[Krs]，只有当下下标r=s，或者r和s同属一个单元节点的号码时，才不为零，其它都为0。整体刚度矩阵[K]，的非零元素成带状分布在主元素附近。2、整体刚度矩阵[K]的主要性质为：1)[K]中每列元素的物理意义是某一节点的坐标轴方向发生单位位移，而其它节点位移都约束为0时，在所有节点上坐标轴方向需施加的节点力；2)[K]的主元素是正定的；[K]是一个对称矩阵；[K]是一个稀疏阵；[K]是一个奇异阵，当排除刚性位移后是正定阵。3、自由度为2；半带宽为8;4、下半带宽的行为第5行、半带宽内元素的列为4列，元素个数为8个，子矩阵个数为4个。五单元类型：四边形环等参单元；速度差值公式ezreVNvvV}]{[}{}{；形函数计算通式4)1)(1(kkKN（k=1,2,34）;形参矩阵]N0N0N0N00N0N0N0N[[N]43214321;自由度f=2;2应变速度列阵erzzre}]{[}{][VBT。。。。。；[B]由几何方程确定，4行8列,T]0111[[C]}C{T;3表面力列阵}0{}{fP；4速度边界条件：r=0时，1rv=3rv=5rv=0;初始速度场vr2=vr4=vr6=v0r0/2h0z=0时，1zv=2zv=0;vz3=vz4=–v0/2z=h0时，vz5=vz6=-v0单元①的速度摄动量列阵：Tzrzri]vvvv[}v{14432,与单元②叠加7：1；5摄动量组合后联立方程组的个数与特点：理论上2n+m=2ⅹ6+2=14个方程组，实际上只有3个方程组，因为7个速度已知，考虑到速度边界条件以摄动法解析时每个单元方程个数是5，结合后的方程组含有7个方程组；拉氏乘子数：2个，叠加结点数2个（结点3、4），摄动后矩阵Tks1}{的结数：7:7；收敛判据：}{}v{v某个小数，式中niiv12}v{，niiv12}v{一名词解释虚位移场：在uS上满足0iv，在U内满足,,1()2ijijjiuu的位移场iu成为虚位移场。[D]矩阵：与材料有关的弹性矩阵，对平面应力状态2101011002ED[B]矩阵：单元应变矩阵，它可以写成分块形式ijmBBBB，而子矩阵1(,,)2iiiibBcijmAcbi由于,,,,,ijmijmbbbccc都是常量，所以矩阵[B]中的元素都是常量，因而单元中各点的应变速率分量也都是常量，故通常称为这种单元为常应变单元。等参单元：具有相同的形式，用同样个数的相应结点值做参数，并具有完全相同的形状函数作为这些结点值前面的系数，具有这样特征的单元称为等参单元。二结合实例说明能量法（变分法）解析的主要步骤：小林史郎三维轧制（1）Hill速度场：''2211,(),()xxyzxxxxxhUydyZdZvvUUvUUhhdxhdxhh，（2）取xcb，得1,xxxUvchb'1xxyxxxxxdbbUyvvychbdxbb'1xxzxxxxxdbbUzvvzchbdxbb（3）由此速度场确定应变速率场为：''xxxxxyxxvbhvvbh，'yxyxxvvbyb，'xzxxvvzhzh'''''21()[2()]22xxxxxzzxzxxxxvhhbhvUxzchxbxhhbh'''''21()[2()]22yyxxxxxxyxxxxvUvbhbbxychxbxbhbb（4）按刚—塑性材料马尔可夫变分原理，范函dfsNNN取极小值内部变形功率：2222000243xxhbldVxxUNdVgNIZdxdydzcbh222222(22,3xyzxyxzN='''''21[2()]2xxxxxxxxbhbbbhbb''''222[()()()(),xxxxxxxxbhbhgbhbhI='''''21[2()]2xxxxxxxxhhbhhhbh摩擦功率：'200004413xxlblbfsfffxvNFvdSmUhdxdyU切向速度不连续量：22''221xxRfxxxxhbyVvUchbUchb在中性面处沿圆周方向相对滑动速度为零，2'210xRxxhVUchb出口断面处的剪切功率0012'20100014{[]()}3bhbyssyxxxlvNUvhdzdyhdyUU（5）里兹法求范函的近似极值函数()xbx，取c=1，()xbx,()xbx满足稳定轧制过程条件，速度边界条件，体积不变条件：a)23020211012222322()xbabaaabxbaxxxllllb)12(,,)naaxc)120,0,0,naax解得12,,,naax从而解得()xbx，展宽022[()]xxBbbxb（6）minmindfsNNN，确定单棍纯轧制力矩上界解：min,2vMWWR(7)120,0,0,naax为(