《有限元习题与答案》由会员分享，可在线阅读，更多相关《有限元习题与答案（14页珍藏版）》请在金锄头文库上搜索。
1、习题2.1 解释如下的概念：应力、应变，几何方程、物理方程、虚位移原理。解 应力是某截面上的应力在该处的集度。 应变是指单元体在某一个方向上有一个U的伸长量，其相对变化量就是应变。表示在x轴的方向上的正应变，其包括正应变和剪应变。 几何方程是表示弹性体内节点的应变分量与位移分量之间的关系，其完整表示如下： 物理方程：表示应力和应变关系的方程某一点应力分量与应变分量之间的关系如下：虚位移原理：在弹性有一虚位移情况下，由于作用在每个质点上的力系，在相应的虚位移上虚功总和为零，即为：若弹性体在已知的面力和体力的作用下处于平衡状态，那么使弹性体产生虚位移，所有作用在弹性体上的体力在虚位移上所做的工就等
2、于弹性体所具有的虚位能。2.2说明弹性体力学中的几个基本假设。 连续性假设：就是假定整个物体的体积都被组成该物体的介质所填满，不存在任何间隙。 完全弹性假设：就是假定物体服从虎克定律。 各向同性假设：就是假定整个物体是由同意材料组成的。 小变形和小位移假设：就是指物体各点的位移都远远小于物体原来的尺寸，并且其应变和转角都小于1。2.3简述线应变与剪应变的几何含义。线应变：应变和刚体转动与位移导数的关系，剪应变表示单元体棱边之间夹角的变化。2.4 推到平面应变平衡微分方程。解：对于单元体而言其平衡方程：在平面中有 代入上式的 2.5 如题图2.1所示，被三个表面隔离出来平面应力状态中的一点，求和
3、的值。解：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是常量。求讨论位移场
4、的相容性解： 因为 所以满足相容性条件有广义胡克定律得又则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 根据弹性力学平面问题的几何方程，证明应变分量满足下列方程，并解释该方程的意义。证
5、明：弹性力学平面问题的几何方程为：将方程，分别对y和x求二阶偏导并相加得：等式右端项，该方程为相容方程中的第一式，其意义为弹性体内任一点都有确定的位移，且同一点不可能有连个不同的位移，应变分量应满足相容方程，否则，变形后的微元体之间有可能出现开裂与重叠。2.14 假设Airy应力函数为，其中为常数，求，并求这些变量间的约束关系。解：由，对该应力函数求偏导得;对以上两式的偏导可求得：考虑相容性条件,将上式代入可得各常量间的关系如下：2.15 对给定的应力矩阵，求最大Tresca和Von.Mises应力。将Von Mises应力和Tresca应力 20 10 10进行比较，= 10 20 10 M
6、pa。10 10 20 z xy xz解：由Tresca准则：= y yz 故有s=20Mpa，max=s/2=10Mpa z1=（x+y）/2=30Mpa 2=10Mpa由Von Mises准则：2s2=6（xy2+yz2+yz2）解得s=30Mpa 30 -15 202.16 一点出的应力状态由应力矩阵给出，即= -15 -25 10 Mpa，若E=70Gpa， 20 10 40 =0.33，求单位体积的应变能。解：单位体积应变能：=1/2Ex2+y2+z2-2u(xy+yz+zz)+2(1+u)(xy2 +xz2+yz2)u=（E-2）/2 =0.33带入可得：=420.75J3.11
7、如图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=-4 b1=x2-x3=-8 b1=y3-y1=4 b1=x3-x1=0 b1=y1-y2=0 b1=x1-x2=8 b1 0 b2
8、 0 b3 0 -4 0 4 0 0B=1/2* 0 c1 0 c2 0 c3 =1/32* 0 -8 0 0 8 c1 b1 c2 b2 c3 b3 -8 4 0 8 0 1 0 1 0.3 0D=E/(1-2)* 1 0 =E/0.91* 0.3 1 0 0 0 (1-)/2 0 0 0.35 1 0.3 0 -0.125 0 0.125 0 0S=D*B=E/0.91* 0.3 1 0 * 0 -0.25 0 0 0.25 0 0 0.35 -0.25 0.125 0 0.25 0 1.4 0 -1.4 -0.7 0 0.7 0 4 -0.6 -4 0 0K=BT*D*B*t*=E/36
9、.4* -1.4 -0.6 2.4 1.3 0.6 0.7 -0.7 -4 1.3 -0.6 -1 0.35 0 0 0.6 -1 -0.6 0 0.7 0 0.7 -0.35 0 0 1 0 0 0.6 -1 -0.6 0 0.35 0.7 0 -0.7 -0.35 0 0.7 1.4 0 -1.4 -0.7K=BT*D*B*t*=E/36.4* 0.6 0 0 4 -0.6 -4 1 -0.7 -1.4 -0.6 2.4 1.3 0.6 -0.35 -1.4 -4 1.3 3.53.12 求下图中所示的三角形的单元插值函数矩阵和应变矩阵，u1=2.0mm，v1=1.2mm，u2=2.4mm
10、，v2=1.2mm，u3=2.1mm，v3=1.4mm，求单元内的应变和应力，求出主应力和方向。若在单元jm边作用有线性分布面载荷（x轴），求结点的的载荷分量。解：如图2=64/3，解得以下参数：a1=19 a2=-2 a3=6； b1=-3 b2=4 b3=-1；c1=-1 c2=-3 c3=4；N1=64/3*(19-3x-y) N2=64/3*(-2-3x-3y)N3=64/3*(6-x+4y)故N= Ni 0 Nj 0 Nm 0 0 Ni 0 Nj 0 Nm 1 0 1 0 1 0 = 0 1 0 1 0 1 bi 0 bj 0 bm 0B=1/2* 0 ci 0 cj 0 cm ci bi cj bj cm bm -3 0 4 0 -1 0=64/3* 0 -1 0 -3 0 4 -1 -3 -3 4 4 -1 1 0D=E/(1-2)* 1 0 0 0 (1-)/2 1 0 -3 0 4 0 -1 0单元应力矩阵S=D*B= E/13(1-2)* 1 0 * 0 -1 0 -3 0 4 0 0 (1-)/2 -1 -3 -3 4 4 -1 2 1.1 -3 -u 4 3u -1 4u 2.4单元应力=S*q= E/13
