机械有限元习题答案-哈工大

第二章习题2.1解释如下的概念：应力、应变，几何方程、物理方程、虚位移原理。解eq\o\ac(○,1)应力是某截面上的应力在该处的集度。eq\o\ac(○,2)应变是指单元体在某一个方向上有一个ΔU的伸长量，其相对变化量就是应变。表示在x轴的方向上的正应变，其包括正应变和剪应变。eq\o\ac(○,3)几何方程是表示弹性体内节点的应变分量与位移分量之间的关系，其完整表示如下：eq\o\ac(○,4)物理方程：表示应力和应变关系的方程某一点应力分量与应变分量之间的关系如下：eq\o\ac(○,5)虚位移原理：在弹性有一虚位移情况下，由于作用在每个质点上的力系，在相应的虚位移上虚功总和为零，即为：假设弹性体在的面力和体力的作用下处于平衡状态，那么使弹性体产生虚位移，所有作用在弹性体上的体力在虚位移上所做的工就等于弹性体所具有的虚位能。2.2说明弹性体力学中的几个根本假设。eq\o\ac(○,1)连续性假设：就是假定整个物体的体积都被组成该物体的介质所填满，不存在任何间隙。eq\o\ac(○,2)完全弹性假设：就是假定物体服从虎克定律。eq\o\ac(○,3)各向同性假设：就是假定整个物体是由同意材料组成的。eq\o\ac(○,4)小变形和小位移假设：就是指物体各点的位移都远远小于物体原来的尺寸，并且其应变和转角都小于1。2.3简述线应变与剪应变的几何含义。线应变：应变和刚体转动与位移导数的关系，剪应变表示单元体棱边之间夹角的变化。2.4推到平面应变平衡微分方程。解：对于单元体而言其平衡方程：在平面中有代入上式的2.5如题图2.1所示，被三个外表隔离出来平面应力状态中的一点，求和的值。解：x方向上：联立二式得：2.6相对于xyz坐标系，一点的应力如下某外表的外法线方向余弦值为，，求该外表的法相和切向应力。解：该平面的正应力全应力该平面的切应力2.7一点的应力如下MP求主应力和每一个主应力方向的方向余弦；球该店的最大剪应力。解：设主平面方向余弦为，由题知将代入得即，。最大剪应力〔1〕当时代入式〔2.21〕〔2〕当时代入式〔2.21〕且2.8一点P的位移场为，求该点p(1,0,2)的应变分量。解：p点沿坐标方向的位移分量为u,v,w点p(1,0,2)处线应变为，，剪应变为，，2.9一具有平面应力场的物体，材料参数为E、v。有如下位移场其中，a、b、c、d是常量。求讨论位移场的相容性解：因为所以满足相容性条件有广义胡克定律得又那么2.10一具有平面应力场的物体，材料性质是E=210GPa,v=0.3.并且有如下位移场当x=0.050m,y=0.020m时，求物体的应力和应变。位移场是否相容？解：由广义胡克定律，，满足相容性条件2.11对于一个没有任何体积力的圆盘，处于平面应力状态。其中a,b,c,d,e,f,g,h是常量。为了使应力满足平衡方程和相容方程，这些常量的约束条件是什么？解：由题意得：，，，代入平衡方程根据广义胡克定律：代入相容方程〔2〕代入〔1〕得其中2.13根据弹性力学平面问题的几何方程，证明应变分量满足以下方程，并解释该方程的意义。证明：弹性力学平面问题的几何方程为：①，②，③，将方程①，②分别对y和x求二阶偏导并相加得：等式右端项，该方程为相容方程中的第一式，其意义为弹性体内任一点都有确定的位移，且同一点不可能有连个不同的位移，应变分量应满足相容方程，否那么，变形后的微元体之间有可能出现开裂与重叠。2.14假设Airy应力函数为，其中为常数，求，并求这些变量间的约束关系。解：由，对该应力函数求偏导得;对以上两式的偏导可求得：考虑相容性条件,将上式代入可得各常量间的关系如下：2.15对给定的应力矩阵，求最大Tresca和Von.Mises应力。将VonMises应力和Tresca应力201010进行比拟，δ=102010Mpa。1020δzτxyτxz解：由Tresca准那么：δ=δyτyz故有δs=20Mpa，τmax=δs/2=10Mpaδzδ1=〔δx+δy〕/2=30Mpaδ2=10Mpa由VonMises准那么：2δs2=6〔τxy2+τyz2+τyz2〕解得δs=30Mpa30-15202.16一点出的应力状态由应力矩阵给出，即δ=-15-2510Mpa，假设E=70Gpa，γ201040=0.33，求单位体积的应变能。解：单位体积应变能：υ=1/2E{δx2+δy2+δz2-2u(δxδy+δyδz+δzδz)+2(1+u)(τxy2+τxz2+τyz2)}u=〔E-2γ〕/2γγ=0.33带入可得：υ=420.75J第三章3.1解释根本概念：位移插值函数、位移模式、单元刚度矩阵及其刚度系数、单元刚度矩阵的对称性和奇异性、结构刚度矩阵的集成、单元载荷向量、有限元解的收敛准那么、位移解的下限性质。答：位移插值函数：建立以单元结点位移表示的单元内各点位移的表达式，选择一个简单的单元位移模式，单元内各点的位移可按此位移模式由单元结点位移通过插值而获得。位移模式：单元位移模式的一般表达式为{d(x,y)}=[f(x,y)]a单元刚度矩阵及其刚度系数:Re=kδe是表示单元的结点力和结点位移之间关系的刚度方程，k就是单元刚度矩阵kBTDBtdxdy=∫∫单元刚度矩阵的对称性和奇异性：对称性是指第j个单位位移分量引起的第i个结点力分量等于由第i个单位位移分量引起的第j个结点力分量。奇异性：指单元刚度矩阵不存在逆矩阵。结构刚度矩阵的集成：对N个经推广的单元刚度矩阵进行求和，叠加，得到结构整体刚度矩阵。有限元解的收敛准那么：满足三个条件1、位移模式包含单元的刚体位移2、位移模式必须能包含单元的常应变3、位移模式在单元内要连续，且在相邻单元之间的位移必须协调位移解的下限性质：对于一个给定的位移模式，其刚度系数是数值比精确值要大。所以在给定的载荷之下，有限元计算模型的变形将比实际结构的变形小。因此细分单元网格，位移近似解将由下方收敛于精确解，即得到真实解的下界。3.11如图3.11所示的平面三角形单元，厚度t=1cm，弹性模量E=2.0*105mpa，泊松比γ=0.3，试求插值函数矩阵N，应变矩阵B，应力矩阵S，单元刚度矩阵Ke。解：此三角形单元可得：2△=〔10-2〕*4=32，故有a1=1/32*〔8u1-5u2-16u3〕a2=1/32*〔4u1-4u2〕a3=1/32*〔-8u1+8u3〕a4=1/32*〔56v1-8v2-16v3〕a5=1/32*〔-4v1+4v2〕a6=1/32*〔-8v1+8v3〕而b1=y2-y3=-4b1=x2-x3=-8b1=y3-y1=4b1=x3-x1=0b1=y1-y2=0b1=x1-x2=8b10b20b30-40400[B]=1/2△*0c10c20c3=1/32*0-8008c1b1c2b2c3b1γ010.30[D]=[E/(1-γ2)]*γ10=[E/0.91]*0.31000(1-γ)/2000.3510.30-0.12500.12500[S]=[D]*[B]={E/0.91}*0.310*0-0.25000.25000.35-0.250.12500.2501.40-1.4-0.700.704-0.6-400[K]①=BT*D*B①*t*△={E/36.4}*-1.4-0.62.41.30.60.7-0.7-41.3-0.6-10.35000.6-1-00.700.7-0.35001000.6-1-0.600.350.70-0.7-0.3500.71.40-1.4-0.7[K]②=BT*D*B②*t*△={E/36.4}*0.6004-0.6-41-0.7-1.4-0.62.41.30.6-0.35-1.4-41.33.53.12求以下图中所示的三角形的单元插值函数矩阵及应变矩阵，u1=2.0mm，v1=1.2mm，u2=2.4mm，v2=1.2mm，u3=2.1mm，v3=1.4mm，求单元内的应变和应力，求出主应力及方向。假设在单元jm边作用有线性分布面载荷〔x轴〕，求结点的的载荷分量。解：如图2△=64/3，解得以下参数：a1=19a2=-2a3=6；b1=-3b2=4b3=-1；c1=-1c2=-3c3=4；N1={64/3}*(19-3x-y)N2={64/3}*(-2-3x-3y)N3={64/3}*(6-x+4y)故N=Ni0Nj0Nm00Ni0Nj0Nm101010=010101bi0bj0bm0[B]={1/2△}*0ci0cj0ccibicjbjcmbm-3040-10={64/3}*0-10-304-1-3-344-11γ0[D]={E/(1-γ2)}*γ1000(1-γ)/21γ0-3040-10单元应力矩阵[S]=[D]*[B]={E/13(1-γ2)}*γ10*0-10-30400(1-γ)/2-1-3-344-121.1-3-u43u-14u2.4单元应力[δ]=[S]*[q]={E/13(1-γ2)}*-3u-14u-3-u4*1.2(u-1)/2(3u-3)/2(3u-3)/22-2u2-2u(u-1)/22.41.43.13解：二维单元在x,y坐标平面内平移到不同位置，单元刚度矩阵相同，在平面矩阵180°时变化，单元作上述变化时，应力矩阵不变化。〔0,1〕〔2,1〕3.14〔0,1〕〔2,1〕〔2,0〕〔0,0〕②①yx解：令，，而，，〔2,0〕〔0,0〕②①yx单元①单元②：由和扩充KZ〔总刚度阵〕而，其中，，化简得：那么，3.15如下图有限元网格，，单元厚度，弹性模量，泊松比。答复下述问题：〔1〕结点如何编号才能使结构刚度矩阵带宽最小？〔2〕如何设置位移边界条件才能约束结构的刚体移动？〔3〕形成单元刚度矩阵并集成结构刚度矩阵。〔4〕如果施加一定载荷，拟定求解步骤。(1)(2)〔3〕解：1、节点编号如图(2)所示；2、如图(3)设置位移边界条件才能约束结构的刚体移动；3、如图(2)所示各节点的坐标为(以m为单位)：1(0,0)，2(0.08,0)，3(0,0.04)，4(0.08,0.04)，5(0,0.08)，6(0.08,0.08)，7(0,0.12)，8(0.08,0.12)解：单元号 1 2 3 4 5 6 相邻结点 1 3 4 5 5 7 2 2 5 4 6 6 3 4 3 6 7 8对于单元号1：；；；；；；对于单元号2：；；；；；；对于单元号3：；；；；；；对于单元号4：；；；；；；对于单元号5：；；；；；；对于单元号6：；；；；；；平面三角形单元的面积均为弹性矩阵均为应变矩阵应力矩阵单元刚度矩阵结构刚度矩阵为：假设施加一定载荷，求解步骤为：1、对单元编号，并列出各单元三个结点的结点号；2、计算外载荷的等效结点力，列出结构结点载荷列阵；3、计算单元刚度矩阵，组集结构整体刚度矩阵4、引入边界条件，即根据约束情况修正结构有限元方程，特别是消除整体刚度矩阵的奇异性，得到考虑约束条件的可解的有限元方程。5、利用线性方程组的数值解法，对结构的有限元方程进行求解，得到所有各结点的位移向量。最后根据需要求解单元应力。3.16一长方形薄板如下图。其两端受均匀拉伸。板长12cm,宽4cm，厚1cm。材料，泊松比。均匀拉力。使用有限元法求解板的内应力，并和精确解比拟〔提示：可利用结构对称性，并用2个三角形单元对结构进行离散〕。解：解：结点编号1 2 3 4 单元号 1 2 X坐标 0 12 0 12 相邻结点 1 3 Y坐标 0 0 4 4 2 2 3 4平面三角形单元的面积均为应力矩阵为：单元1的应变距阵为：单元1的单元刚度矩阵为：单元2的应变距阵为：单元2的单元刚度矩阵为：总刚度矩阵为：位移分量为：载荷列阵为：因为可以得单元1的单元应力：单元2的单元应力：长方形薄板内应力的精确解为：拉应力，用有限元法求解出的结果与精确解大致相等。3.17验证三角形单元的位移差值函数满足及。解：平面三角形形函数为：，其中，，分别是行列式2A中的第一行，第二行和第三行各元素的代数余子式。行列式中，任一行的元素与其相应的代数余子式的乘积之和等于行列式的值，而任一行的元素与其它行对应元素的代数余子式乘积之和为零，故有：当，同时有，同理也有：，即。3.18推导如下图的9节点矩形单元的形函数。解:三维杆单元的形状函数，①在局部坐标系中令节点1,5,2所对应的带入①式得到节点1,5,2仅在x方向上的形函数：②同理可得：由，即节点2,6,3，可得到沿着全局坐标系y轴的形状函数(通过变量轮换)，节点1的形函数即x，y方向的乘积：由此可得：同理可整理得：，，，，，，3.19如下图为一个桁架单元，端点力为[U1，U2]，端点位移为[u1，u2]，设内部任一点的轴向位移u是坐标x的线性函数：推导其形函数矩阵N。解：轴向位移u是坐标x的线性函数，，写成向量形式为，设两个节点的坐标为，代入向量形式的位移函数解得：那么由位移函数可得形函数为：答：轴对称三角形环单元不是常应变单元，如果弹性体的几何形状、约束条件及载荷都对称于某一轴，那么所有的位移应变及应力也是对称于此轴，这样问题称为轴对称。轴对称三角形环单元与平面常应变单元是不同的，轴对称三角形环单元的应变不是常数矩阵，其应变矩阵B=[BBB],其中B=，〔i,j,m〕。应变分量，，都是常量，但环向应变不是常量，它与，，中的r和z有关。答：轴对称问题中，刚度自由度：环向位移，径向位移，轴向位移。以三角环单元平均半径、平均高度进行计算的单元刚度矩阵，配合以精确积分所得的等效结点载荷矩阵，计算的结果还是不错的！轴对称问题的两个单元a和b，设材料的弹性模量为E，泊松比为μ=0.15，试手算这两个单元的刚度矩阵。解：对于单元，由题可知：单元a的截面面积为单元a的刚度矩阵写成分块矩阵形式为：其中子矩阵可写为：所以的刚度矩阵为对于单元，由题可知单元的截面面积为单元的刚度矩阵写成分块矩阵形式为：其中子矩阵可写为：所以单元的刚度矩阵为5.1答：杆件受到纵向〔平行于杆轴〕载荷的作用，这样杆件的拉压问题；杆件受到横向〔垂直于杆轴〕载荷的作用，这是梁的弯曲问题。杆件受到力相似到薄板就有，薄板受到纵向载荷的作用，这是平面应力问题；薄板受到横向载荷的作用，这是薄板的弯曲问题。薄板的弯曲可以认为是梁弯曲的推广，是双向的弯曲问题，中面法线在变形后保持不伸缩，并且成为弹性曲面的法线，中面在变形后，其线段和面积的投影形状保持不变〔小挠度薄板〕。中面的挠度，而纵向位移、，主要应力分量，，。某一点的位移：，，。某一点的应力：,,弹性曲面微分方程,其中……板的抗挠刚度。5.2答：矩形薄板单元：薄板单元位移函数并不满足连续性或相容性要求，采用这种位移函数的单元是非协调单元，这种四节点矩形弯曲单元变形后，其挠度面在单元间虽然互相连续，但其法向导数并不连续，单元间在变形后是不连续光滑〔有棱〕的，当单元逐渐取小的时候，还能够收敛于精确解。三角形薄板单元：常使用面积坐标，分析说明，只以挠度及其一阶导数作为节点的位移函数用一般的形状函数是不可能构造满足相容性的薄板单元，需再加上二阶导数，就可以实现。在相邻单元之间，挠度是连续的，但法向的斜率是不连续的，这种位移模式是非协调单云，收敛不如矩形单元，单元足够小，节点增多，如六节点三角形，九节点三角形等。5.3谈论在平面应力和弯曲状态组合的情况下，三角形刚度矩阵的特点平面内的作用力产生的变形不影响弯曲变形，反之亦然节点把转向在两种应力状态下都不参加到变形中，相应的节点力也不存在，将平面应力状态和弯曲状态加以组合后，单元的每个节点的位移向量和节点力向量是要指出的是，在局部坐标系中，节点位移不包括，但为了下一步将局部坐标系的单元刚度阵换到总体坐标系下进行集成，由于平面应力状态下的节点力和平面应力状态下的节点位移互不影响，弯曲应力状态下的节点与平面应力状态下的节点位移互不影响，所以组合应力状态下的平板、薄板单元的单元刚度矩阵如下：，=其中矩阵和分别是平面应力问题和薄板弯曲问题的相应子矩阵，三角形单元的单元刚度矩阵是18×18矩阵。6.1结构的动态特性：结构的固有频率及其相应的模型，以及在随着时间而变形的外加激振力的鼓励下，机器或结构被激起的位移，应力或称被激起的动力响应，机械产品的动态性能是其重要的性能指标，尤其对现代复杂、高速、重载精密机械系统，动态性能是影响其工作性能及产品指标的关键技术指标，机械结构的动态特性问题早在上个世纪30年代就引起人们的重视，动态特性的开展为机械动态设计提供了坚实的根底。6.2结构离散后，在运动状态各节点的动力平衡为：其中,,分别以惯性力、阻尼力和动力载荷均为矢量，为弹性力，弹性力矢量可用