|
本帖最后由 獨孤峰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
|