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

均布荷载作用下的简支梁结构有限元分析1（总页）--本页仅作为文档封面，使用时请直接删除即可--元分析1（总页）--内页可以根据需求调整合适字体及大小--哈工程有限元大作业均布荷载作用下简支梁结构分析院（系）名称：船舶工程学院专业名称：港口航道与海岸工程学生姓名：白天华学号：03摘要本文利用ANSYS软件中的BEAM系列单元建立简支梁有限元模型，对其进行静力分析与模态分析，得出梁的结构变形，分析梁的受力情况。并用有限元刚度矩阵知识求解简支梁端点处得位移和旋度。在此基础上，利用经典力学对以上所得的结果进行梁的有关计算，并将结果与有限元刚度矩阵和ANSYS软件所得结果进行比较。通过比较得出不同方法在简支梁求解过程中自己的优势和缺点。1.问题求解问题描述钢制实心梁的截面尺寸为10mmX10mm（如图1所示），弹性模量为200GPa，均布荷载的大小及方向如图1所示。图1利用力学方法求解运用力学方法将上述结构求解，易得A、B支座反力相等为500N，该简支梁的计算简图、弯矩图以及剪力图如下图所示1000N/m.11IIn1000mm_

图2简支梁计算简图跨中弯矩：125N.皿图3简支梁弯矩图支座反力500N图4简支梁剪力图利用ANSYS软件建立模型与求解通过关键点创建实体模型，然后定义材料及单元属性，然后划分网格，建立有限元模型。具体步骤包括：添加标题、定义关键点、定义直线、选择单元，定义实常数、定义材料属性、设定网格尺寸、划分网格、施加荷载求解（选择分析类型、定义约束、施加荷载）查看分析结果。图5简支梁变形前后的情况

Hi日匕Load±n*j图6简支梁应力图图7简支梁剪力图2计算结果对比简支梁内力分析结果比较节点应力有下面公式计算求得：6=iIz有限元计算所得结果与力学的计算结果对比如下表所示:单位（N/皿）

PRINTELEMENTTABLEITEMSFERELEMENT*****POSTlELEMENTTAELELISTINGSTfiTCUREENIELEMSMAMI0.18629E-11270.60480.60S30.60720.6075^.60720.60630.60480.60270.60NIHinUMUALUESELEM1UALUE0.10629E-11HAMinUMUALUESELEM6节点应力102270348046305720675077208630948010270VALUE节点应力102270348046305720675077208630948010270ANSYS模态结果

结构力学计算结果简支梁竖向位移分析结果比较（1）结构力学计算求得的简支梁最大位移由下面图乘法求得:实际荷载作用下梁弯矩表达式:M(x)=500x-500x2单位荷载作用下梁弯矩表达式：Mp=(1-a)x(0<x<a)a(1-x)(a<x<1)则在梁上任意点的竖向位移f：f=500「"+500i心I」*dxU£2-13121=)分别代入分段点的a的数值得各点的位移如下表:有限元计算所得简支梁y方向位移如下图8所示:LOADSTEP-1SUBSTEP-1TIME=1LOADCfi£E=UTHEFOLLOWINGDEGREEOFFREEDOMRESULTSARNODE0-S000-24.525-46.400-63.525-7S.125-74=400-63.E2E10-4G.40011-24.525MAXIMUMABSOLUTEUALUESUALUE-78.125图8端点旋度分析结果比较（1）利用结构力学图乘法求得端点处得旋度旋度：如（；L；qJ-:）=顼|（2）利用有限元刚度矩阵求得端点位移与旋度为:假设梁的两端固定，并计算等价的节点荷载用以表示均匀变化的荷载力M1-M2cR2R1M1R2,-1/2qL、.-1/12qL2126L.-126Lv16L4L2-6L2L201-1/2qL=EI/L3-12L-6L12-6Lv2(a)1/12qL2J6L2L2-6L4L202方程（a）是固定的精确模型，因为如果从中解出的所有位移和旋度，它们的计算值都将为零。利用边界条件，得到矩阵方程:202-WL2/202L2-WL2/30=EI/L34L2202-WL2/202L2解方程组（b），得每个点处得旋度大小为:01=02=qL3/24EI（c）用实际节点荷载代替作用在梁上的荷载力力，计算出最后的反作用力：加上由节点旋度引起的反作用126lr-126L01/2qLM1fEI/L36一4L2-6L22L2-qL01=02=qL3/24EI（c）用实际节点荷载代替作用在梁上的荷载力力，计算出最后的反作用力：加上由节点旋度引起的反作用126lr-126L01/2qLM1fEI/L36一4L2-6L22L2-qL3/'24日1/12qL2R2-12-6L12-6L0qL/2（d）求解矩阵方程，得到最终结果:R1=qL/2R2=qL/2M1=M2=03结论（1）本文通过ANSYS有限元软件中BEAM3单元建立了简支梁模型，经过同种工况的力学静力分析，简支梁应力、位移结果相同。（2）用有限元刚度矩阵法求得的简支梁端点位移与旋度的结果和经典结构力学求得的结果一致。（3）对静定简支梁的分析，有限元软件ANSYS能直观的观察梁的各种物理变化，经典力学求解方法相对刚度矩阵法更加简洁方便，但刚度矩阵法对更加复杂结构的求解相对更方便。参考文献李家宝.结构