七七影院色七七_免费观看欧美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 n
clc;
c1 w5 K3 I9 v' u2 B0 M
L1=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:m
7 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 `: Q
XQM=[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! B
Q1=[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& {/ y
M2=[x1(i),Mmin];
( u- J! Y- x' t: o. [* E" k
QDX=[Q1;Q2]
% v9 O8 Q0 r4 ]( N' ~( y* Q6 l( z
MDX=[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; }, K
end
* 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# c
for j=1:n
0 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; x
end
" u+ u1 G9 q; d* L+ |% a, ^* P8 m
QQ(n1:n)=QQ(n1:n)-P;
* n/ _: t/ S0 D, u- B
MM(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/ U
function [QQ,MM]=QMQ(n,x1,c,d,q,QQ,MM)
) N: ^. w8 Y6 A3 R) p l$ Y) v8 M
for 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)>d
1 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% v
function TU_QM(x1,QQ,MM)
9 C5 b5 e0 J( f$ T% b) v0 h
subplot(2,1,1);
# D) S% ?9 j8 A+ b0 a
plot(x1,QQ);
. I5 r) e9 g' q( W# H' M+ A* O
grid;
2 |( ^1 h& Y# U1 m
title('剪力圖');
2 J4 g# i( y K! x6 P7 O
subplot(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; }' Z
grid;
& { 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