七七影院色七七_免费观看欧美a一级黄片_亚洲综合久久久久久中文字幕_国产999999在线视频免费观看,国产小视频无码,国产精品亚洲日日摸夜夜添,女人高潮潮叫免费网站,久久影院国产精品,日韩成人在线影院,欧美囗交XX×BBB视频,色在线综合高清

機械社區(qū)

 找回密碼
 注冊會員

QQ登錄

只需一步,快速開始

搜索
查看: 4827|回復: 0
打印 上一主題 下一主題

[matlab] 求助:關于材料力學課程設計matlab程序

[復制鏈接]
跳轉到指定樓層
1#
發(fā)表于 2013-5-21 11:11:37 | 只看該作者 回帖獎勵 |倒序瀏覽 |閱讀模式
本帖最后由 獨孤峰yi 于 2013-5-21 11:13 編輯
4 A4 B' U/ r- |, `. S: L0 |0 u/ k/ f' l" u4 ^
材料力學課程設計《左端外伸梁的剪力與彎矩圖》,,由于沒有matlab基礎,早上七點半糾結到現(xiàn)在實在是hold不住了。我在書上找到了例題有正確答案卻不知如和用程序表達,。以下是題和我寫的程序,,不知道哪里錯了,希望了解的能幫幫我,,謝謝�,。�, o' y" N6 D0 z. O  w
題:. n# f! G& `4 z* {+ |5 c/ \
如圖所示F=20KN,Me=160KNm,q=20KN/m,試做梁的剪力圖和彎矩圖
1 j# A  k( Z! ~9 q主程序:
' w- {  R0 [6 R8 b0 A' a) T# j/ Aclear;
( V# M: y! k0 \: F% ?* @4 ~clc;
& N1 A  U) b) ]L1=2;
6 J+ {2 B3 F3 l7 z+ X* \8 {! zx=[0,2,10,12];
& P* P, o$ C- L9 G8 {9 H; AMPQ = [2  20  0 0
3 e' N5 g& D& h. @+ S7 T       3  20  0 10
7 [7 Q( K% S( p; e8 C       1  160 10 0];
0 {: C" ^, q6 m# ]" r/ p' UXQM=QMZW(x,L1,MPQ);3 b( [7 a, {! D
第一子程序:
; A! {* j4 U) @; S  dfunction XQM=QMZW(x,L1,MPQ)
" F, F1 H0 e; D+ |6 X! B: I3 g7 Y   [n,m]=size(x);/ w1 L3 [) f' H( v( R7 S
   L=x(m);
2 y; C/ p8 N: R5 u' b! R   L1=2;
9 q9 {* |" V. C) e! ^1 _, {   x1=[];1 b& S  a$ C8 B% c$ m3 @
   for i = 1: m-1
9 z" T: f" Q6 d# K         x1 = [x1,linspace(x(i),x(i+1),50)];! k- j) L6 ^3 n& \8 g
end                         6 L, h, G6 M, `/ H9 R) z
   MM = zeros(size(x1));8 y3 m  ~% a5 v- O" ?& v
   QQ = zeros(size(x1));. }" O4 @% O% a4 w4 u) q) h
   [m,t] = size(MPQ);
/ f+ r: M8 N& Z" i6 e; O7 V+ {2 i   [t,n] = size(x1);
, R' s$ F8 H3 m) ^; c9 t   for i=1:m
9 I+ M, b) d$ D2 I& Q" J         switch MPQ(i,1)1 P9 d2 d, l% _3 ~4 |+ t; _
case 11 I# M# Z( O7 z* ?' K; o; u+ K
                        a=MPQ(i,3);
4 {# V3 D8 u. P) ~! ]( @; g                        M=MPQ(i,2);
. Z& ~2 `: _+ ^  u* [! ^                        RA=M/(L-L1);! ]! {2 y/ G  N9 k3 d
                        RB=-RA;
( Y! s1 r9 S# n5 {5 w6 N3 Y0 V                        QQ=QQ+RA;' v3 N* e; \! a. t
                        MM=MM+RA*x1;1 q) V; @' E. F# P) `6 J6 q0 ]0 |* B
                        if a>= 0&&a<L
3 s: D: _; o4 {                            MM=QMM(n,x1,a,M,MM);                                % ]2 z* T* d. Q% Y" T1 t
                        end& z2 w( b! B9 ^
                        if a==0
' Q& H% ]9 ?- C7 e% a- f( L/ Z3 l2 m                           MM=MM+M;
( |7 w. r+ G1 x                         end : M% L; l, u% m' I& g
                           [QQ,MM]=QMP(n,x1,L-L1,-RB,QQ,MM);               
$ _( ?. ^/ R3 g  case 28 b8 U& Y; c% i, F. f! G
                        P=MPQ(i,2);
  z9 E2 y. a6 \                        b=MPQ(i,3);
) T( p8 W2 ^+ Z* Y- D+ S( Q9 s) v" g- A0 ]                        RA=P*(L-b)/(L-L1);
. y5 h3 N5 r( Y' U5 G4 W/ I. i                        RB=P+P*(L-b)/(L-L1);; Q. j2 K+ @$ G2 K
                        QQ=QQ+RA;
" Q  w3 F; k& L                        MM=MM+RA*x1;3 q6 R1 i. M' _' a/ m1 \
                        if b>0&&b<L
8 P8 l1 Z& Z2 }0 c                           [QQ,MM]=QMP(n,x1,b,P,QQ,MM);   
8 G0 \1 d7 P& m/ `/ k                        end% r& W4 [# a' Z% I$ s% q: m
                       [QQ,MM]=QMP(n,x1,L-L1,RB,QQ,MM);' q% J, f/ [% b9 W2 O5 K) `: G5 k$ [# n7 A
case 3% V- M% e$ Y/ S' C5 l. n6 j
                        q=MPQ(i,2);
5 f4 i" y, i6 Y                        c=MPQ(i,3);
6 f0 X" s& K: w                        d=MPQ(i,4);7 K; n! j, S+ ?1 ~  z5 D# l4 Z
                        RA=q*(d-c)*(L-0.5*(d+c))/(L-L1);
3 d) i3 N: s) b$ A. C: [- L. y                        RB=q*(d-c)-q*(d-c)*(L-0.5*(d+c))/(L-L1);
1 i0 }# v! G( @6 q                        QQ=QQ+RA;; b( N' I& Y: A* D
                        MM=MM+RA*x1;: D8 ?+ Z' F$ X. H
                        [QQ,MM]=QMQ(n,x1,c,d,q,QQ,MM); ! s8 v( h. U  p9 M
                        [QQ,MM]=QMP(n,x1,L-L1,RB,QQ,MM);% v5 o' b4 T1 Q7 D5 s6 U; @9 ?
         end
1 c0 L) J$ b9 _    end4 R: ?7 P- Y  Q
    [QDX,MDX,XQM]=MAX_MIN(x1,QQ,MM);        
. S( H+ \$ ^7 c    TU_QM(x1,QQ,MM);                           # Q) u& i+ U2 R# ?
第二子程序:+ Q5 H% @5 s9 {+ ~% p
function [QDX,MDX,XQM]=MAX_MIN(x1,QQ,MM)
- p, X; L9 A8 r" ]XQM=[x1',QQ',MM'];
. q3 p' A- d' a5 R! x6 \3 ]. `[Qmax,i]=max(QQ);3 Q9 [6 q$ L8 h! B& u; M- c
Q1=[x1(i),Qmax];+ ?! U: _3 ?5 \( ?
[Qmin,i]=min(QQ);2 g5 }1 V0 j! H7 o# T3 r
Q2=[x1(i),Qmin];' F, W7 J# t! P8 S$ O
    [Mmax,i]=max(MM);
8 x: W+ ]9 m, U; E# ^2 _    M1=[x1(i),Mmax];
5 s& q6 ~$ X7 _7 g- Q  T# L[Mmin,i]=min(MM);
& S! p2 f6 i3 ~$ s9 c1 A+ {  nM2=[x1(i),Mmin];2 H; `% l* _0 l" S9 A1 U
QDX=[Q1;Q2]. E- {) H) D% J- }
MDX=[M1;M2]
/ q  l* [0 r' O8 M1 I8 }! {( v第三子程序:
6 e' \; ?8 a/ ]0 d( y) Afunction MM=QMM(n,x1,a,M,MM)
, |( F7 m! A. p/ n' {0 d" ]) X( bfor j=1:n# o8 R# p6 d( P: o6 O4 E, \
    if x1(j)==a
$ g: E6 W% t) q1 _4 \9 o        n1=j;0 O6 F' j. u; l2 I: @: s
    end
/ M5 C" k% m* z9 o' eend* D6 K, N' [. W5 d6 N5 r# g- m
    MM(n1:n)=MM(n1:n)+M;1 C7 O/ R' ]# j- t. b& S  h
第四子程序:
) ]1 |3 b6 o/ Y  B. v% c+ M5 Ffunction [QQ,MM]=QMP(n,x1,b,P,QQ,MM)4 K' e0 t- r6 I- ?) Y% j
for j=1:n
" x/ e4 P0 \. C  y0 @    if x1(j)==b* n  R9 g1 C9 j% o
        n1=j;
: U! d8 V6 Z( t7 S2 v+ ]: j2 b* C    end
& I) ~' V. I/ s" _end
" G/ ?, l3 ~8 |. |2 E5 mQQ(n1:n)=QQ(n1:n)-P;
. C9 A4 R2 w1 I0 k0 P& h" L4 fMM(n1:n)=MM(n1:n)-P*(x1(n1:n)-b);
' V! R% \- V2 V第五子程序:( I, W. J) q0 ]6 l" @* R
function [QQ,MM]=QMQ(n,x1,c,d,q,QQ,MM)
/ w7 ~1 ]# ]8 q6 B! Pfor j=1:n
. Z: t* u$ E1 G0 p& N    if x1(j)>c# L) x4 g& G. _9 a. m9 G
        QQ(j)=QQ(j)-q*(x1(j)-c);
( m+ x; R: U% e        MM(j)= MM(j)-0.5*q*(x1(j)-c)^2;1 z( |( z! b# Z( [% ?2 T6 {4 ~; j
    end
- v( w* O) ?) l& D( Y) P+ b      if x1(j)>d
3 U# ?; @' ^6 W          QQ(j)=QQ(j)+q*(x1(j)-d);7 i" I0 a, ^. x4 {; Q, j4 \
           MM(j)= MM(j)+0.5*q*(x1(j)-d)^2;
% Z: L& p3 `, i5 ^  c9 a$ Q      end6 E& {5 o+ s2 w) V" t
end
8 a1 \$ I% P  D6 b第六子程序:
2 ?1 r' E% R) p. q, S) pfunction TU_QM(x1,QQ,MM)$ f, q0 s7 N: \& b
subplot(2,1,1);- k0 p3 M4 m" b0 m; e
plot(x1,QQ);
" d: B9 i9 @* f' X. z" w: B0 Mgrid;( N: h4 {. _4 k2 D9 N9 h' q8 e
title('剪力圖');
* v% C" u3 Z6 i9 F  w  Q* S8 usubplot(2,1,2);! h& Q& o$ ^  W4 d0 J# M, q: g# ^) Y
plot(x1,MM);
3 u5 t- D$ t9 _( Q( Ogrid;
! @- T' Y3 S4 }) o( Xtitle('彎矩圖');
4 X8 U, Z0 N5 X) u
5 T9 V! \0 `" I; m% G

本帖子中包含更多資源

您需要 登錄 才可以下載或查看,,沒有帳號,?注冊會員

x
您需要登錄后才可以回帖 登錄 | 注冊會員

本版積分規(guī)則

小黑屋|手機版|Archiver|機械社區(qū) ( 京ICP備10217105號-1,京ICP證050210號,,浙公網(wǎng)安備33038202004372號 )

GMT+8, 2025-4-26 13:04 , Processed in 0.056700 second(s), 16 queries , Gzip On.

Powered by Discuz! X3.4 Licensed

© 2001-2017 Comsenz Inc.

快速回復 返回頂部 返回列表