[matlab] 求助：關(guān)于材料力學(xué)課程設(shè)計(jì)matlab程序

[復(fù)制鏈接]

獨(dú)孤峰yi

電梯直達(dá)

1^#

發(fā)表于 2013-5-21 11:11:37 | 只看該作者 |倒序?yàn)g覽 |閱讀模式

本帖最后由獨(dú)孤峰yi 于 2013-5-21 11:13 編輯

材料力學(xué)課程設(shè)計(jì)《左端外伸梁的剪力與彎矩圖》，由于沒有matlab基礎(chǔ)，早上七點(diǎn)半糾結(jié)到現(xiàn)在實(shí)在是hold不住了。我在書上找到了例題有正確答案卻不知如和用程序表達(dá)。以下是題和我寫的程序，不知道哪里錯(cuò)了，希望了解的能幫幫我，謝謝啊！
題：
如圖所示F=20KN,Me=160KNm,q=20KN/m,試做梁的剪力圖和彎矩圖
主程序：
clear;
clc;
L1=2;
x=[0,2,10,12];
MPQ = [2  20  0 0
   3  20  0 10
   1  160 10 0];
XQM=QMZW(x,L1,MPQ);
第一子程序：
function XQM=QMZW(x,L1,MPQ)
[n,m]=size(x);
L=x(m);
L1=2;
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-L1);
                     RB=-RA;
                     QQ=QQ+RA;
                     MM=MM+RA*x1;
                     if a>= 0&&a<L
                        MM=QMM(n,x1,a,M,MM);
                     end
                     if a==0
                        MM=MM+M;
                     end
                        [QQ,MM]=QMP(n,x1,L-L1,-RB,QQ,MM);
  case 2
                     P=MPQ(i,2);
                     b=MPQ(i,3);
                     RA=P*(L-b)/(L-L1);
                     RB=P+P*(L-b)/(L-L1);
                     QQ=QQ+RA;
                     MM=MM+RA*x1;
                     if b>0&&b<L
                        [QQ,MM]=QMP(n,x1,b,P,QQ,MM);
                     end
                     [QQ,MM]=QMP(n,x1,L-L1,RB,QQ,MM);
case 3
                     q=MPQ(i,2);
                     c=MPQ(i,3);
                     d=MPQ(i,4);
                     RA=q*(d-c)*(L-0.5*(d+c))/(L-L1);
                     RB=q*(d-c)-q*(d-c)*(L-0.5*(d+c))/(L-L1);
                     QQ=QQ+RA;
                     MM=MM+RA*x1;
                     [QQ,MM]=QMQ(n,x1,c,d,q,QQ,MM);
                     [QQ,MM]=QMP(n,x1,L-L1,RB,QQ,MM);
      end
end
[QDX,MDX,XQM]=MAX_MIN(x1,QQ,MM);
TU_QM(x1,QQ,MM);
第二子程序：
function [QDX,MDX,XQM]=MAX_MIN(x1,QQ,MM)
XQM=[x1',QQ',MM'];
[Qmax,i]=max(QQ);
Q1=[x1(i),Qmax];
[Qmin,i]=min(QQ);
Q2=[x1(i),Qmin];
[Mmax,i]=max(MM);
M1=[x1(i),Mmax];
[Mmin,i]=min(MM);
M2=[x1(i),Mmin];
QDX=[Q1;Q2]
MDX=[M1;M2]
第三子程序：
function MM=QMM(n,x1,a,M,MM)
for j=1:n
if x1(j)==a
      n1=j;
end
end
MM(n1:n)=MM(n1:n)+M;
第四子程序：
function [QQ,MM]=QMP(n,x1,b,P,QQ,MM)
for j=1:n
if x1(j)==b
      n1=j;
end
end
QQ(n1:n)=QQ(n1:n)-P;
MM(n1:n)=MM(n1:n)-P*(x1(n1:n)-b);
第五子程序：
function [QQ,MM]=QMQ(n,x1,c,d,q,QQ,MM)
for j=1:n
if x1(j)>c
      QQ(j)=QQ(j)-q*(x1(j)-c);
      MM(j)= MM(j)-0.5*q*(x1(j)-c)^2;
end
   if x1(j)>d
      QQ(j)=QQ(j)+q*(x1(j)-d);
         MM(j)= MM(j)+0.5*q*(x1(j)-d)^2;
   end
end
第六子程序：
function TU_QM(x1,QQ,MM)
subplot(2,1,1);
plot(x1,QQ);
grid;
title('剪力圖');
subplot(2,1,2);
plot(x1,MM);
grid;
title('彎矩圖');

本帖子中包含更多資源

您需要登錄才可以下載或查看，沒有賬號(hào)？注冊(cè)會(huì)員

分享到: QQ好友和群 QQ空間

收藏0 轉(zhuǎn)播 分享 淘帖0 問題專業(yè)，描述清楚0 伸手黨/灌水/看不懂0

回復(fù)

使用道具舉報(bào)

返回列表發(fā)新帖

		自動(dòng)登錄	找回密碼
密碼			注冊(cè)會(huì)員