资源描述:
《用matlab解决简支梁问题.docx》由会员上传分享,免费在线阅读,更多相关内容在学术论文-天天文库。
1、如何用matlab绘制简支梁剪力图弯矩图说明:输入变量分段数组x分段点一般在集中力,集中力偶作用出和分布载荷的起末端。载荷数组MPQ若梁上的外载荷总数为PN,则用PN行四列的数组MPQ储存载荷,数组MPQ第一列代表载荷的类型:1为集中力偶,2为集中力,3为分布载荷,第二列代表载荷的大小,第三列代表集中力,集中力偶或者分布载荷左端与简支梁左端的距离,第四列代表均匀载荷右端与简支梁左端的距离,当载荷为集中力或者集中力偶时,第四列为0.符号规定集中力和均匀载荷向下为正,向上为负,集中力偶顺时针为正,逆时针为负。输出变量内力数组XQM如果梁被分为NN-1段,则内力数组XQM为NN行,三列的数组,第
2、一列代表梁的横截面的位置,第二列代表剪力,第三列代表弯矩。剪力极值及位置QDXQDX是一个二行二列的数组,第一列代表极值所在的位置,第二列代表极值弯矩极值及位置MDXMDX是一个二行二列的数组,第一列代表极值所在的位置,第二列代表极值子程序集中力偶对弯矩贡献的子函数QMM集中力对剪力和弯矩贡献的子函数QMP分布载荷对剪力和弯矩贡献的子函数QMQ求剪力和弯矩极值的子函数MAX_MIN绘制剪力图和弯矩图的子函数TU_QM计算分析程序简支梁QMDJ左端固定悬臂梁QMDXZ右端固定悬臂梁QMDXY左端外伸梁QMDWZ右端外伸梁QMDWY两端外伸梁QMDWLMatlab求解简支梁挠度曲线方程(201
3、4-01-1114:09:45)转载▼标签:截面gui挠度分类:Matlab项目经验载荷工程力学1. 项目概述使用Matlab求解简支梁挠度曲线方程,并将求解程序使用GUI图形界面实现。载荷分布可选:均布载荷、直线分布载荷、二次分布载荷。梁截面可选:矩形截面、环形界面、工字钢截面。 2. 关键词GUI;梁;挠度;载荷;截面;工程力学 3. 技术说明首先根据载荷分布曲线求解剪力分布方程,剪力分布积分得到弯矩分布方程,弯矩积分得到转角方程,转角积分得到挠度方程。然后根据简支梁的边界条件(简支梁两端支点挠度等于0)代入解出积分常数,即可得挠度方程。GUI运行模式:首
4、先接受用户输入参数,得到载荷分布参数、截面参数、材料弹性模量、梁长度等参数,用户点击Calculate(计算)按钮后,调用相应程序计算出挠度方程显示在相应文本框,并作出挠度曲线。 4. 程序源代码(1)均布载荷简支梁挠度曲线方程计算symsq0xlF=q0*l;%等效集中力Fa=F/2;%A处支反力Fb=F/2;%B处支反力qx=q0;%载荷密度Fx=Fa-int(qx,x);%剪力分布Mx=int(Fx,x);%弯矩分布symsC1C2%常数Sita=int(Mx,x)+C1;%转角wx=int(Sita,x)+C2;%挠度%根据边界条件确定常数eq1=subs(wx,x,0)
5、;%A处挠度为0eq2=subs(wx,x,l);%B处挠度为0g=solve(eq1,eq2,C1,C2);%解方程求C1C2C1=g.C1;C2=g.C2;wx=-subs(wx); (2)二次分布载荷简支梁挠度曲线方程计算symsq0xlF=0.5*q0*l;%等效集中力Fa=F/3;%A处支反力Fb=F*2/3;%B处支反力qx=q0*x/l;%载荷密度Fx=Fa-int(qx);%剪力分布Mx=int(Fx);%弯矩分布symsC1C2%常数Sita=int(Mx)+C1;%转角wx=int(Sita)+C2;%挠度%根据边界条件确定常数eq1=subs(wx,x,0);%A处挠
6、度为0eq2=subs(wx,x,l);%B处挠度为0g=solve(eq1,eq2,C1,C2);%解方程求C1C2C1=g.C1;C2=g.C2;wx=-subs(wx); (3)线性分布载荷简支梁挠度曲线方程计算symsq0xlqx=4*q0*(x/l-x^2/l^2);%载荷密度F=int(qx,x,0,l);%等效集中力Fa=F/2;%A处支反力Fb=F/2;%B处支反力Fx=Fa-int(qx);%剪力分布Mx=int(Fx);%弯矩分布symsC1C2%常数Sita=int(Mx)+C1;%转角wx=int(Sita)+C2;%挠度%根据边界条件确定常数eq1=subs(wx
7、,x,0);%A处挠度为0eq2=subs(wx,x,l);%B处挠度为0g=solve(eq1,eq2,C1,C2);%解方程求C1C2C1=g.C1;C2=g.C2;wx=-subs(wx); 5. 程序结果 6. 总结本项目的主要难点:一是挠度计算;二是绘图程序的编写,如铰支座、载荷分布曲线、梁截面的绘制。