无网格法大作业.docxVIP

二、利用伽辽金无网格法实现悬臂梁弯曲问题 算例描述：悬臂梁自由度受集中力作用如图1所示，编写2D FEGM程序计算位移、 应变、应力。 梁的尺寸：L = 42m, D = 12m,厚度默认为lm; 弹性常数：E = 30MPa,v = 0.3; 不考虑梁的自重，在自由端施加的集中载荷为P = -1000 N; 图1自由端受集中力作用的悬臂梁 弹性力学解析解 TOC \o 1-5 \h \z pv n2 d x 方向位移：以兀,y) = 一[(6L 一 3x)x + (2 + v)(^2 + —) -D(y + —] 6 El 4 2 P ￡)2 y 方向位移：v(x,y) = -—[3v^2(￡-x) + (4 + 5v) —— + (3￡-x)x2] 6EI 4 梁截面的法向应力：cyxx^y) = -P(L~X)y P n2 y方向正应力:crvv = 0梁截而剪应力：(x,y) = — (— -,2) J ￡ t 应变能：￡ = -[己 DzdG u 4.45 1、节点间距收敛性分析 为研究节点间距对计算结果的影响，考虑如下5种等间距节点分布情况： 3X3(轴向X竖向)、6X3、11X3、11X5 11X7。依次增加轴向和竖向的节点数, 考察应变能残量和悬臂梁中心轴y=o的计算位移的收敛性，以及某一梁截面 (x=23.667)的计算应力收敛性。此部分程序与相应结果见附录lo