梁计算机分析程序.doc

下载文档 降价啦

1
0
约4.66千字
约 9页
2017-02-13 发布于重庆
举报
版权申诉
保障服务

梁计算机分析程序.doc

1、本文档共9页，可阅读全部内容。
2、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。
3、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击这里二次下载。
4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。

梁计算机分析程序

简支梁计算机分析程序 function XQM = QMDJ(x,MPQ) [n,m] = size(x); L = x(m); x1 = []; for i = 1:m-1 x1 = [x1,linspace(x(i),x(i+1),50)]; end %离散简支梁 MM = zeros(size(x1)); QQ = zeros(size(x1)); %将离散点的剪力和弯矩数组QQ、MM赋零 [m,t]=size(MPQ); [t,n]=size(x1); for i = 1:m; switch MPQ(i,1) case 1 a=MPQ(i,3); M=MPQ(i,2); RA=-M/L; QQ=QQ+RA; MM=MM+RA*x1; if a0 aL MM=QMM(n,x1,a,M,MM); end if a==0 MM=MM+M; end case 2 P=MPQ(i,2); b=MPQ(i,3); RA=(L-b)*P/L; if b0 bL QQ=QQ+RA; MM=MM+RA*x1; [QQ,MM]=QMP(n,x1,b,P,QQ,MM); end case 3 q=MPQ(i,2); c=MPQ(i,3); d=MPQ(i,4); RA = (L-0.5*(c+d))*q*(d-c)/L; QQ=QQ+RA; MM=MM+RA*x1; [QQ,MM] = QMQ(n,x1,c,d,q,QQ,MM); end end [QDX,MDX,XQM] = MAX_MIN(x1,QQ,MM); TU_QM(x1,QQ,MM); 分析实例 clear; clc; x=[0,2,4] MPQ=[2,30,2,0] XQM=QMDJ(x,MPQ); 右端外伸梁计算机程序 function XQM=QMYW(x,L2,MPQ) [n,m]=size(x); L=x(m); x1=[]; for i=1:m-1 x1=[x1,linspace(x(i),x(i+1),50)]; end MM=zeros(size(x1)); QQ=zeros(size(x1)); [m,t]=size(MPQ); [t,n]=size(x1); for i=1:m switch MPQ(i,1); case 1 a=MPQ(i,3); M=MPQ(i,2); RA=-M/(L-L2); RB=RA; QQ=QQ+RA; MM=MM+RA*x1; MM=QMM(n,x1,a,M,MM); [QQ,MM]=QMP(n,x1,L-L2,RB,QQ,MM); end case 2 P=MPQ(i,2); b=MPQ(i,3); RA=P-Pb/(L-L2); RB=-PL/(L-L2); QQ=QQ+RA; MM=MM+RA*x1; [QQ,MM]=QMP(n,x1,b,P,MM); [QQ,MM]=QMP(n,x1,L-L2,RB,QQ,MM); case 3 q=MPQ(i,