板壳理论课程设计.docx

1、板壳理论课程设计变 z 可以不板壳理论课程设计第一部分 综述这学期我们学习了板壳理论，也就是弹性力学的下册。经过这学期的学习，我对弹性力学的概念也越发的清晰， 也认识到自己在公式推导方面的不足。 不过作为一名力学专业的学生， 这是最基本的数学素质要求。 通过学习板壳理论， 我对弹性力学问题的分析思路更加清晰， 尤其是对薄板问题的处理， 有了更好的认识和理解。弹性力学的研究对象是完全弹性体，根据问题的性质，忽略一些很小的次要因素，对物体的材料性质采用了一些基本假定， 即弹性力学的基本假定， 主要有连续性、完全弹性、均匀性、各向同性，符合以上假定的物体，就称为理想弹性体。弹性体是变形体的一种， 在

2、外力作用下物体变形， 当外力不超过某一限度时，出去外力后，除去外力后物体即恢复原状。在板壳理论中，这一原则的应用则更加广泛。在薄板的小挠度弯曲理论中，垂直于中面方向的正应计 ；应力分量 zx, zy , z 所引起的变形可以不计； 薄板中面内的各点都没有平行于中面的位移， 这是以三个计算假定为基础的。 例如在薄板弯曲问题中， 一定载荷引起的弯应力和扭应力， 在数值上最大， 因而是主要的应力； 横向切应力在数值上较小，是次要的应力；挤压应力在数值上更小，是更次要的应力。因此，在计算薄板的内力时，主要是计算弯矩和扭矩，横向剪力一般都无须计算。这学期板壳理论的学习，让我在上学期弹性力学上册的基础上有

3、了新的收获。上学期弹性力学的学习， 我感觉整本书就讲了十五个控制方程解十五个未知数。而剩下的问题就是如何求解这些方程的问题， 这也是数学和力学结合最紧密的地方。而求解的方法无外乎两种 :基于位移的求解和基于应力的求解，而前人的研究大部分都是如何使这些方程求解起来更方便。 弹性力学思路清晰， 但是方程和公式复杂。例如 ,应力函数的引入就是因为同时满足平衡方程和应力表达的相容方程是很难找到的。 这学期，我们学习了伽辽金法的应用及其举例， 认识了伽辽金位移函数 它使得原本要求的方程 (非齐次微分方程 )转化为求拉普拉期方程 ,而拉普拉斯方程在数学上（复变函数）已经研究的很透彻，因而大大简化了求解的难

4、度。而近代即二十世纪以来发展起来的能量法更是如此:对位移的变分方程代替了以位移表达的平衡方程及应力边界条件， 对应力的变分代替了相容方程及位移边界条件 这无疑都大大简化了弹性力学基本方程的求解过程。随着计算机的发展，各类软件也应用在了各行各业，这学期我们还学习了ABAQUS和 ANSYS两类有限元软件。 此外，通过数学软件 Matlab 和 Mathematica与有限元软件的结合应用， 使得求解速度大大加快， 便于方便实验， 这也使得许多从前很难解决的问题基本上都能获得满足工程精度的解答。 在传统理论解和有限元方法共同合作下， 弹性力学的发展会更加迅速， 它的应用范围更加广泛， 前景是非常可

5、观的。第二部分 解答题目：四边为夹支边的正方形薄板受均布荷载集度的解法1.里茨法设有一个正方形薄板，边长为 a=1m,厚度 =0.015m 如图所示，四边均为夹支边，在薄板受有均布荷载的作用， q0 100000N / m2 。取坐标轴如图所示，则有位移边界条件为( w )( w )( w )( w )x0xay0y a0,(w)x0,(w)x0,(w)y0,(w)yx0xay0y a0;0;0;0;在薄板的弯曲问题中， 一定荷载引起的弯应力和扭应力， 在数值上最大，因而是主要的应力；横向切应力在数值上较小，是次要的应力。因此在计算薄板的内力时 i ，主要是计算弯矩和扭矩，横向剪力一般都无需计

6、算。而由基尔霍夫指出，薄板任一边界上扭矩都可以变换为等效的横向剪力，和原来的横向剪力合并，内力表达式为MMMF2w2 w),xD (2y 2x2w2 w),yD (2x 2yxyM yxD (12 w),x ySxD4 w , F SyD4 w ,xy在里茨法中，内力边界条件可忽略，因此将挠度的表达式取为w C1w1C1 * x(xa) y( ya)sinxsinyaa（a）则上列位移边界条件都能满足，同时，式（a）在薄板的四边还满足了内力边界条件，即弯矩不等于零。薄板的形变势能表达式为：V1(x xy xyxyxy)dxdydz1(xy xyyzyzzxzx)dxdydz (b)22在薄板的

7、小挠度弯曲问题中，按照计算假定，z , yz , zx 形变分量不计，于是形变势能的表达式简化为2 w2w(2 w )2dxdy( w2 w )(w2 w ) dxdyx2y2x yx x y2y x x yV1xy )dxdydz(xxy xyxy(c)2xEz2(2 w2w),xz2w,1x2y22x根据物理方程Ez(2 w2 wz2w代入式（ c）y12y2x2 ),yy2xy1Ez2 w ,xy2z2wx yx y整理后得VE2(2 22(1)2w 2w2w)2dxdydz2)zw)22(2(1xyx y各项与 Z 无关V1D(2w)22(1)2w 2w(2w2dxdy222)（d）x

8、yx y在等厚度薄板中， D是常量，DE3式（ d）可以写为12(12)D2 22w 2w2w2V2(w)2(1 )x2y2( x y)dxdy（e）或写为 VD(22)D2w 2w(2w2dxdy （ f ）2w) dxdy (122)xyx y其中2 w 2w(2w)2dxdyw 2w(w2wx2y2x y(2) dxdyx x yy x x y由格林定理xP( x, y)Q( x, y)dxdyQ (x, y)dxP( x, y)dy .y得2 w 2 w(2wdxdyw2 ww 2 wdy,其中右边的积x2y2x)xdxxy2yx y分是沿薄板的边界进行的。本题中薄板全部边界条件都是夹

9、支边，有w0x式（ f ）可以简化为 VD( 2 w)2 dxdy（ g）2按式（ a）求挠度 w 对于坐标的二阶导数，得到2 w4C1 *sin2xy 2 (ay)22C1 sin 2 ( ax) y2 (ay)2x2a2a22 w4C1 *sin2x2 (ax) y22C1 sin 2 x2 (a x)(ay)2y2a2a2得 VD(2 w)2 dxdyD * 176a6C12 sin 44221575从而VD * C1 *176 a6sin44C11575（h）qw dxdy1a6 sin22m144(i)由 式 （ h ） 和 (i)求出C1， C12.98e 5代 入 式 （ a ）

10、 , 得w2.98e8* x( xa) y( ya)sinxsinyaa由题意知薄板中心a,a处挠度最大，为max1.8625e 6m22w2. 差分法2.1用 4*4网格求解ha。由于对称，只有3 个独立的未知值，即4w1 , w2 , w3 ，取坐标如下所示oX一般来说，夹支边外一行虚结点32处的挠度，就等于边界内一行相对3结点处的挠度，但计算精度较低，212有时引起很大误差。对夹支边来说，假定 w在边线上按三次式变化，323即 A B ( x ) C ( x ) 2D ( x )3hhh以边界结点4 为原点，则有边界条件Y图 1w(w)x00,(x )x00,(w)x hw3 ,(w)x

11、2hw2求得A 0,B 0,C2 w 3w 2 , Dw 3w 244从而得出 w(2 w 3w 2 )( x ) 2( w 3w 2 )( x ) 34h4h据此边界外虚结点 a, b, c,的挠度表达式分别为wa3w3w2 ；wb3w2w1 ； wc3w3w2222为结点 1,2,3建立差分方程如下：20 w18(4 w2 )2(4 w3)a)4q0,(D420 w28( w12w3 ) 2( w2 ) ( w2wb ) ( a)4 q0 ,4 D20 w38(2w2 )2w12w32( wc )( a )4 q0 ,4D带入得到20 w18(4 w2 )2(4 w3 )( a ) 4 q

12、0 ,4D20 w8( w2w )2(w )(w3ww1 )( a) 4 q0 ,21322224D20 w38(2 w2 )2w12w32(3w2w2 )( a ) 4 q0 ,24D整理得出关于 w 的线性方程组矩阵如下：20328w118.52816w2=2.60099e6 1216.528w1由此得到该3 个结点处的挠度为（单位m ）：w1 2.162e6，w21.378e6，w3 0.878e 6其中最大挠度为 wmaxw12.162e 62.2 用 8*8 的网格求解ha。由于对称，取1薄板为研究对象， 建立如84下坐标系，并标注结点如图所示o 123711245812356101

13、378109141112131415.同 4*4 网格差分法，边界外虚结点 a, b, c, d ,e 的挠度分别为：aw aw 2；.3 w3.bX2.cwb3w8w5 ； wc 3w10w6 ；.d22.e.Yw3ww10 ； w3ww13d9e1522为结点建立差分方程如下：图 220w8(4w )2(4w )4w( a)4 q0,12438D20w8(wwww )2(wwww )2www ( a)4q0 ,2434125525278D20w38(w22w5w7 )2(2 w42w8 )2w6w11w1( a)4q0 ,8D20w48(2w52w2 ) 2(w12w3w6 )2w42w8

14、( a )4 q0 ,8D20w58(w3w6w4w8 )2(w2w7w10w5 )w2w12w1020w68(2w52w10 )2(w4w8w8w9 ) 2w32w13a)4q0(D820w78(2w8w3w11 ) 2(2w52w2 )w22w10waa)4q0(D820w88(w5w7w10 )2(w3 w6 )w4w8wbw9( a)4 q0 ,8D20w98(w10w10 )2(w6 )2w82wd( a)4q0 ,8Dw15 ( a)4 q0 , 8 D,20w10 8(w8 w6 w9 ) 2(w5 w10 ) w5 wc w7( a) 4q0 .8D整理得出关于 w 的线性方程

15、组矩阵如下：2032 48000000w118 258 16601000w 211820416 28400w31216 42216 20200w 410388 2382803w50.162562 e10022162004216w661018.5 040231602w 7100218.5 282418w81000002022316.5w91000038.5 18825w101由此得到 10 个结点挠度如下 ( 单位 m ) ：w1 2.003e 6,w2 1.794e 6,w3 1.218e 6,w4 1.609e 6,w5 1.096e 6,w6 0.754e 6, 其中最大挠度为 w1 2.003e 6w7 0.466e 6,w8 0.422e 6,w9 0.128e 6,w10 0.298e 63有限元解法3.1 建立一个三维实体正方形薄板，边长为 1m， 厚度 0.015m3.1.1 建模分网（ 3D 实体）3.1.2 定义荷载及边界条件（四边夹支）3.1.2 求解查看结果结果：方板中心处挠度最大，最大挠度为 wmax 1.960e 6m3.2 建立一个二维壳，边长为 1m，厚度 0.015m ，与三维实体条件一样，求解并查看结果如下结果：方板中心处挠度最大，最大挠度为 wmax 2.006e 6m三种解法的解法比