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