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

機(jī)械社區(qū)

標(biāo)題: 求助:關(guān)于材料力學(xué)課程設(shè)計(jì)matlab程序 [打印本頁(yè)]

作者: 獨(dú)孤峰yi    時(shí)間: 2013-5-21 11:11
標(biāo)題: 求助:關(guān)于材料力學(xué)課程設(shè)計(jì)matlab程序
本帖最后由 獨(dú)孤峰yi 于 2013-5-21 11:13 編輯 # o* ?: l% M& m, }- q: P$ ~% s( X7 Z

& r/ {1 O+ p5 H- U: u% q材料力學(xué)課程設(shè)計(jì)《左端外伸梁的剪力與彎矩圖》,,由于沒有matlab基礎(chǔ),,早上七點(diǎn)半糾結(jié)到現(xiàn)在實(shí)在是hold不住了。我在書上找到了例題有正確答案卻不知如和用程序表達(dá)。以下是題和我寫的程序,不知道哪里錯(cuò)了,希望了解的能幫幫我,謝謝啊,!
8 P/ ~$ q/ ], Y- L; v! p$ B題:6 u- r4 p1 D6 p& ~" y6 C8 S  h  A7 o
如圖所示F=20KN,Me=160KNm,q=20KN/m,試做梁的剪力圖和彎矩圖
/ s9 M6 E4 t5 @% j- s5 z# f主程序:  e) J* x$ {7 o- {# s7 W
clear;
) |( g3 F- V; j) g4 nclc;
  c1 w5 K3 I9 v' u2 B0 ML1=2;" H! h2 i$ A0 g: G: l5 i# X
x=[0,2,10,12];- b. U, L' S. j+ [' @( d
MPQ = [2  20  0 0
7 }, X" D  @& S3 {       3  20  0 10
- S4 F/ ?1 C& @: p2 h+ f       1  160 10 0];' U- r! m. k+ p# {* \" i
XQM=QMZW(x,L1,MPQ);
8 B8 c5 C/ a! O2 u& w* }第一子程序:/ b* l" n4 ?& f
function XQM=QMZW(x,L1,MPQ)
  h2 y4 h; ?# K/ C- _% ]! I* L5 b   [n,m]=size(x);
, M0 Z6 u3 r& Y& ^! ^5 R4 }: Y+ Y   L=x(m);8 }' J' n8 m2 v, S% u
   L1=2;
# [: D  u7 }* X: i4 O   x1=[];
6 H4 c+ _+ Z% v9 L1 z. D   for i = 1: m-1
7 t; Z; g2 O/ y! T6 |+ q+ B         x1 = [x1,linspace(x(i),x(i+1),50)];  f+ U* M4 {8 j. A8 @. J% O( B
end                        
: e- k* ^% Q! A" u   MM = zeros(size(x1));. d" i2 S. }1 U
   QQ = zeros(size(x1));3 z# y# f& ~4 B) P1 i% S
   [m,t] = size(MPQ);8 }) W5 J- P! @* e
   [t,n] = size(x1);1 x0 ?# c1 _/ @  y* W4 y" t" A
   for i=1:m7 B/ V1 h9 P' |6 V- }, [
         switch MPQ(i,1)7 L1 d1 y5 S: H- ?, A
case 1
  b0 E& A. E2 h) d( P9 {# H/ x                        a=MPQ(i,3);
9 S. R+ l: S4 u! `4 e# p                        M=MPQ(i,2);
; ~: G6 L9 }0 ?5 ~, @0 K                        RA=M/(L-L1);# u4 K  s6 h7 R" B$ U+ y$ U
                        RB=-RA;6 r5 G0 Q& d, f, y) m$ q( L' ]
                        QQ=QQ+RA;
  g$ j# b) l3 Z' e* @4 ?& P6 w                        MM=MM+RA*x1;; p. j9 m1 }% u, a  ]' l
                        if a>= 0&&a<L
- M; R# `$ j' w$ t3 {7 y* {: s9 Y                            MM=QMM(n,x1,a,M,MM);                                
7 L& T2 H9 H$ n$ Z, f# m* f( L                        end( X! r) ^# Y" W& t3 c- x8 r
                        if a==0
( B+ r9 ]7 z0 p3 k                           MM=MM+M;; d+ F: Y5 Y+ A  v2 N. ^
                         end   k5 Y+ L0 x# O) ]! L& `, {! K" o. Y0 e" _
                           [QQ,MM]=QMP(n,x1,L-L1,-RB,QQ,MM);               
  b5 t/ k! E& Z  case 2, N, a) O# \" G. H% y- S
                        P=MPQ(i,2);, S+ E& m+ s. m
                        b=MPQ(i,3);" O- o4 {  t, U! g) X; {+ H
                        RA=P*(L-b)/(L-L1);
9 r! i9 L3 F* R! l7 ~# q7 N                        RB=P+P*(L-b)/(L-L1);
3 t; u! g- x1 y                        QQ=QQ+RA;0 k2 a; w" M5 S) K- U
                        MM=MM+RA*x1;2 K; c/ }2 P& K+ N, T- m. s. O
                        if b>0&&b<L
, h, S# A- M2 c4 H$ E                           [QQ,MM]=QMP(n,x1,b,P,QQ,MM);    7 Y7 G; G; U1 x9 G9 h8 ~
                        end
# s5 S7 J% H8 |2 _                       [QQ,MM]=QMP(n,x1,L-L1,RB,QQ,MM);: T: f8 b8 c) g1 q+ n
case 3* @7 i( Q7 p( N9 v6 k
                        q=MPQ(i,2);
, K( Q% W4 Q* E                        c=MPQ(i,3);) a; s, T" D, ]! F% H
                        d=MPQ(i,4);
. m" j4 a) i" m  P) n% H                        RA=q*(d-c)*(L-0.5*(d+c))/(L-L1);6 v. w/ j4 F7 g) a5 v
                        RB=q*(d-c)-q*(d-c)*(L-0.5*(d+c))/(L-L1);
4 O4 E/ K0 S! a; |                        QQ=QQ+RA;
0 k& u$ ~' |+ f3 i5 V% w                        MM=MM+RA*x1;
5 V+ N/ }% U$ A                        [QQ,MM]=QMQ(n,x1,c,d,q,QQ,MM); 0 Z; Q' G" E- p0 m8 A
                        [QQ,MM]=QMP(n,x1,L-L1,RB,QQ,MM);1 O8 u/ I" |$ v. ^3 `; B
         end
1 o4 H1 S6 y' i; B    end
, U. {, V( ?- i  P    [QDX,MDX,XQM]=MAX_MIN(x1,QQ,MM);        
2 j' m, z& q0 M    TU_QM(x1,QQ,MM);                           . `7 n# S: C9 D( x
第二子程序:: K- p' |; E+ v9 b0 I, e' Y" E3 i2 _
function [QDX,MDX,XQM]=MAX_MIN(x1,QQ,MM)
( l- E7 R  m# u/ q7 `: QXQM=[x1',QQ',MM'];
) n; @3 k# \4 n9 G, P2 p% [[Qmax,i]=max(QQ);
7 q8 G5 p- ^0 D4 i- w; Q5 S$ I; {" T! BQ1=[x1(i),Qmax];& q% ^* J2 q" }1 l. N. s! l
[Qmin,i]=min(QQ);) j; H; \# y  B$ a" u
Q2=[x1(i),Qmin];! m, e4 D0 w5 P& G4 i8 O
    [Mmax,i]=max(MM);
; S" h4 k# c' f' f    M1=[x1(i),Mmax];
' o3 ~: `( ?5 W5 Z[Mmin,i]=min(MM);
7 l  d2 J2 t* m% x& {/ yM2=[x1(i),Mmin];
( u- J! Y- x' t: o. [* E" kQDX=[Q1;Q2]
% v9 O8 Q0 r4 ]( N' ~( y* Q6 l( zMDX=[M1;M2]
! N; H3 Q6 l9 }( K( y1 \6 C第三子程序:+ g2 p; B2 z: t/ W: V3 G5 u( v; m
function MM=QMM(n,x1,a,M,MM)2 o  W& D  o$ a1 x" y' \, b
for j=1:n! h" l+ W* s' U8 y1 v
    if x1(j)==a. X& d* `* L' s- C2 Q3 C: O! L: B6 W$ u
        n1=j;
" |* K  b4 G! ^% N    end
) o  A6 b3 l, Q, r6 M; }, Kend
* H) I' S4 n$ F& N9 V) T/ M    MM(n1:n)=MM(n1:n)+M;- v9 P) e7 J. h' ?1 i( N5 _0 D
第四子程序:8 w$ |! @6 i( n; X' a% F1 W
function [QQ,MM]=QMP(n,x1,b,P,QQ,MM)
9 W) r6 f5 e  D# cfor j=1:n0 l2 m/ L; L, X+ Q; i- }0 z
    if x1(j)==b
$ W! U& ?+ O* l5 Z        n1=j;6 O. d- M4 p  e3 q
    end
* e& F1 m- N, [3 I$ V4 o; xend
" u+ u1 G9 q; d* L+ |% a, ^* P8 mQQ(n1:n)=QQ(n1:n)-P;
* n/ _: t/ S0 D, u- BMM(n1:n)=MM(n1:n)-P*(x1(n1:n)-b);9 R* ^3 q" L4 m6 }4 \& Y
第五子程序:
) C  u# _6 E" o2 A3 q/ Ufunction [QQ,MM]=QMQ(n,x1,c,d,q,QQ,MM)
) N: ^. w8 Y6 A3 R) p  l$ Y) v8 Mfor j=1:n. M- M$ w5 Q- D2 R( r' B
    if x1(j)>c
) S0 H& {5 b+ g& n        QQ(j)=QQ(j)-q*(x1(j)-c);
$ R  U- z' |; V1 O        MM(j)= MM(j)-0.5*q*(x1(j)-c)^2;
2 X4 ^4 U8 O4 C) r( s2 ]    end$ J: B, n) W  f; R0 G. x/ W4 z' P
      if x1(j)>d1 R5 u7 w- _' e5 `* ~
          QQ(j)=QQ(j)+q*(x1(j)-d);7 n3 ^* R/ i" W# ~- M3 l
           MM(j)= MM(j)+0.5*q*(x1(j)-d)^2;4 D) S. V: x1 M
      end! M3 H! t3 ]& G% ?  v. P  T6 D
end! d- X  o& M2 C6 c1 ^. w4 p) {8 |* M
第六子程序:
; D) m! Z' e7 o: x% vfunction TU_QM(x1,QQ,MM)
9 C5 b5 e0 J( f$ T% b) v0 hsubplot(2,1,1);# D) S% ?9 j8 A+ b0 a
plot(x1,QQ);
. I5 r) e9 g' q( W# H' M+ A* Ogrid;2 |( ^1 h& Y# U1 m
title('剪力圖');
2 J4 g# i( y  K! x6 P7 Osubplot(2,1,2);* p. u) a& e1 u) R# \" o( e8 v4 c
plot(x1,MM);
& W7 N3 a, D8 J. F% h* Y& X, W; }' Zgrid;& {  y4 e9 A8 H
title('彎矩圖');! M( i: `: u* a) O/ N
# p& q8 F" Y2 ~6 E' a" O





歡迎光臨 機(jī)械社區(qū) (http://97307.cn/) Powered by Discuz! X3.4