哈工大材料力学上机程序.doc

《哈工大材料力学上机程序.doc》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。

1、材料力学电算大作业题目名称：有限差分法求梁的弯曲变形通用程序作者班号：XXXXXXXXX作者学号：XXXXXXXXX作者姓名：XXXXXX指导教师：完成时间：2013年6月27日评语：成绩（满分10分:签名:时间:有限差分法求梁的弯曲变形通用程序设计说明书作品简介本作品采用强大的数据处理软件Matlab编程，可以计算简支梁和悬臂梁在复杂受力载荷(弯矩，集中力，分布力)作用下的弯矩，并画出图形。本作品有良好的图形交互GUI界面，使用方便。同时本作品已经编译成独立运行的可安装・exe软件，只要目标计算机上安装有MCR编译器，本程序可以完美运行。1确定求解流程:选择梁的类型2数学模型建立和

2、对应matlab程序文件复杂受力情况下，采用叠加法求梁的任一点的弯矩。2.1简支梁的弯矩的数学模型图一复杂受力载荷F的简支梁2.1叠加法求任意一点弯矩：M(x)M_M(X)=-M/L*x(0<=x

3、2*x-d-c)(d<=x

4、f)-F*(x(n_b+l:nf)-b);M_F=[M_F1,M_F2];%分布力产生的弯矩n_c=round(c/dx)+l;n__d=round(d/dx)+l;M_ql=(L-0.5*(c+d))/L*q*(d-c)*x(l:n_c);M_q2=(L-0.5*(c+d))/L*q*(d-c)*x(n_c+l:n_d)-0.5*q*(x(n_c+l:n_d)-c).A2;M_q3=(L-0.5*(c+d))/L*q*(d-c)*x(n_d+l:nf)-0.5*q*(d-c)*(2*x(n__d+l:nf)-d-c);M_q=[M_ql,M_q2,M_q3];%总的弯矩M_x=M_

5、M+M_F+M_q;end3创建GUI图形交互界面GUI图形交互界面特点：1给出操作提示，方便参数输入；2同时界面绘出弯矩图；GUI部分代码(绘弯矩图代码)：%简支梁读入参数L=str2num(get(handles.edit_L,'string'))；M=str2num(get(handles.edit_M,'string1));a=str2num(get(handles.edit_a/string1));F=str2num(get(handles.edit^B'string'))；b=str2num(get(handles.edit_b,'string'))；q=str2num(

6、get(handles・edit_q/string'));c=str2num(get(handles.edit_c/string'))；d=str2num(get(handles.edit_d,'string'));El=str2num(get(handles.edit_EI/string));%求弯矩，画弯矩图axes(handles.axeslB);cla[Mx]=M_x(L,M/a/F/b/q/c/d);line(linspace(0,L/3001),Mx);gridon;title(•简支梁的弯矩图J;xlabelCx/m*)；ylabel(,M_x/kN*m,);GUI界面

7、展示：4算例4』简支梁的算例:例题1（材料力学书上例7・5改编，求弯矩图。此处给出具体数值）L=2a=2;q=10KN/m/M=qaA2=10KN*m/a=l,F=0,b=0/c=0,d=l.;力加载方式与程序界面提示方向相同。程序运行结果:*选择简支梁入区.1/F简支梁的弯矩图x/mQL(m)M(KN*m)a(m)F(KN)b(m)q(KN/m)10简支0的W距c(m)d(m)WtErttttAE简支梁的弯矩图L(m)M(KNwm)a(m)21°F(KN