均布荷载作用下的简支梁结构有限元分析1.docxVIP

哈工程有限元大作业 均布荷载作用下简支梁结构分析 院（系）名称： 船舶工程学院 专业名称： 港口航道与海岸工程 学生姓名： 白天华 学 号: 03  摘要 本文利用ANSY湫件中的BEAM系列单元建立简支梁有限元模型，对其进行 静力分析与模态分析，得出梁的结构变形，分析梁的受力情况。并用有限元刚度 矩阵知识求解简支梁端点处得位移和旋度。在此基础上，利用经典力学对以上所 得的结果进行梁的有关计算，并将结果与有限元刚度矩阵和 ANSYSB件所得结 果进行比较。通过比较得出不同方法在简支梁求解过程中自己的优势和缺点。 1.问题求解 问题描述 利用力学方法求解 运用力学方法将上述结构求解，易得A、B支座反力相等为500N,该简支梁 的计算简图、弯矩图以及剪力图如下图所示 _L1000N/m _L 1000mm 图2简支梁计算简图 跨中弯矩：125N . 图3简支梁弯矩图 利用ANSY散件建立模型与求解 通过关键点创建实体模型，然后定义材料及单元届性，然后划分网格，建宜 有限元模型。 具体步骤包括：添加标题、定义关键点、定义直线、选择单元，定义实常数、 定义材料届性、设定网格尺寸、划分网格、施加荷载求解（选择分析类型、定义 约束、施加荷载）查看分析结果。 图5简支梁变形前后的情况 2计算结果对比 简支梁内力分析结果比较 节点应力有下面公式计算求得: My ? . Iz 有限元计算所得结果与力学的计算结果对比如下表所示: 单位(N/ m2) 下点应力 下点 应力 1 0 2 270 3 480 4 630 5 720 6 750 7 720 8 630 9 480 10 270 结构力学计算结果 PRINT ELEMENT TABLE ITEMS PER ELEMENT ***** P0ST1 ELEMENT TABLE LISTING ***** STAT CUBFEH7 EL￡M SMAXI 0.18639E-11 270,60 4G0.B0 TOC \o 1-5 \h \z 630.丽 720. 7S^.0S ? 720.60 630.晌 4朋.丽 10 270.63 MIHIMUM VALUES ELEH 1 UALUE B.18629E-11 NAM HUM UALUES ELEM 6 UftLUE X队前 ANSYS1态结果 简支梁竖向位移分析结果比较 结构力学计算求得的简支梁最大位移由下面图乘法求得: 实际荷载作用下梁弯矩表达式: M(x)=500x-500/ 单位荷载作用下梁弯矩表达式： Mp= J (1-a)x (0xa) 1 a(1-x) (ax1) 则在梁上任意点的竖向位移f: f=500j. 11 心+500 ” 77. dx V JZi-1 6 LJ. =?????)? 分别代入分段点的a的数值得各点的位移如下表： a 位移 0 有限元计算所得简支梁y方向位移如下图8所示： LOAD STFP= 1 SIIRSTFP= 1 IIML- 1 .UldVld LOAD GnSE= ? THF FOL LOU [NG DEGRFF OF FPFFDOM RFSIJLTS AR NODF a. uuidu H?Q0Q6 -24-S2K 4 5 -63.525 b -V4.4MU 7 -78.12E R -74.4Of? -4.d00 -24.52S MfiWIMUM ABSOLUTE UALUES NODE 7 UAWE -7B.125 端点旋度分析结果比较 利用结构力学图乘法求得端点处得旋度 旋度：仲Lx打|金 利用有限元刚度矩阵求得端点位移与旋度为: 假设梁的两端固定，并计算等价的节点荷载用以表示均匀变化的荷载力 M1R1『 M1 R1 『-1/2qL) -1/12qL2 -1/2qL J =EI/L3 I 1/12qL2 「12 6L -12 6L 八6L 4L2 -6L 2L2 -12L -6L 12 -6L , 6L 2L2 -6L 4L2 方程(a)是固定的精确模型，因为如果从中解出的所有位移和旋度，它们 的计算值都将为零。利用边界条件，得到矩阵方程： r -? L2/30 ]=EI/L?「4L2 2L2「？1] -?L2/20」 I 2L2 4L2 1 ?/ (b) 解方程组(b),得每个点处得旋度大小为： ?1 =?2=qL3/24EI (c) r R1)r r R1) r 12 6L -12 6R (0 ] 「1/2qL) M1 =EI/L3 6L 4L2 -6L2 2L2 -qL3/24EI + 1/12qL2 4 R2 -12 -6L 12 -6L 0 qL/2 * I M2) 6L 2L2 -6L2 4L2 i qL3/24EIj 、-1/12qL2 力，计算出最后的反作用力: 求解矩阵