七七影院色七七_免费观看欧美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 編輯
; o5 G( h3 {3 ~: j8 Q; W
/ A9 d, D7 q( P$ Y" I. I$ I6 L* N
材料力學(xué)課程設(shè)計(jì)《左端外伸梁的剪力與彎矩圖》,,由于沒有matlab基礎(chǔ),,早上七點(diǎn)半糾結(jié)到現(xiàn)在實(shí)在是hold不住了。我在書上找到了例題有正確答案卻不知如和用程序表達(dá),。以下是題和我寫的程序,,不知道哪里錯(cuò)了,,希望了解的能幫幫我,謝謝�,�,!
% G c) r6 M! A* n4 c) h
題:
' ]) }" v% a$ }0 O, R; Z) f
如圖所示F=20KN,Me=160KNm,q=20KN/m,試做梁的剪力圖和彎矩圖
* Z8 C5 f1 i' N7 B" s" v& B
主程序:
; v) y$ p+ j8 R9 K: T; J
clear;
, e0 C6 ^" ?6 ^" J$ N7 ~
clc;
: Q- A: a! E! Q5 I5 W; R% q5 m& t
L1=2;
2 x& U4 T) c& m7 J6 F5 U* u
x=[0,2,10,12];
0 g6 p" a& N8 T. B4 @% j
MPQ = [2 20 0 0
1 j2 E3 S, `; m
3 20 0 10
8 k j3 @0 ]; X, W
1 160 10 0];
. K9 F5 Y4 K X
XQM=QMZW(x,L1,MPQ);
% N% F( }9 G3 T: d
第一子程序:
1 {6 f s3 y% _: n
function XQM=QMZW(x,L1,MPQ)
7 O, R, v4 k% O5 c& J
[n,m]=size(x);
9 I2 a7 n( D$ j) V( j$ G
L=x(m);
& o! T7 v: v2 p
L1=2;
& Y0 |! u' S4 w8 p% P; Z
x1=[];
3 z6 R# j2 O% A2 K
for i = 1: m-1
- i3 w# B0 ^% J3 i! C& u
x1 = [x1,linspace(x(i),x(i+1),50)];
/ i5 h- e& ~3 {& h* l% j3 z
end
3 a0 A g( ^7 y% O, Z$ U) e
MM = zeros(size(x1));
3 j7 {: y& k( A- M5 N' U7 {. r
QQ = zeros(size(x1));
% K) z* w. V; M" y4 i
[m,t] = size(MPQ);
( \" g7 D2 d# g! p$ p) y- o( e: l
[t,n] = size(x1);
/ k1 M$ K& I4 x# G' Z9 ~2 `
for i=1:m
+ H1 s; u8 f' d* f7 a
switch MPQ(i,1)
8 Q) l5 l) u# u* }4 g* |# [ Q
case 1
2 T4 [/ x& o! j/ i0 A( n$ A& k. q* L
a=MPQ(i,3);
$ b' a& N/ L8 I, ~9 J
M=MPQ(i,2);
) c& B; |1 x# v; u- ?$ f, w, ^
RA=M/(L-L1);
) K1 Z" P+ ^+ w/ v
RB=-RA;
0 y2 ?: K3 G( K: l3 W
QQ=QQ+RA;
G0 _- I, z, c, x8 u( i
MM=MM+RA*x1;
& H, w/ l; M" i; h/ s( B
if a>= 0&&a<L
8 y+ P5 l- b1 x) ~. G
MM=QMM(n,x1,a,M,MM);
; l5 e) c! e2 U* S
end
5 Y$ j/ u: @; W2 u3 ~
if a==0
. b' f* p! n0 P% O
MM=MM+M;
4 ?4 C* O( ~5 g
end
) R' S. U' d+ K. b# E. n
[QQ,MM]=QMP(n,x1,L-L1,-RB,QQ,MM);
9 ^5 W( O( q* P0 O, ?: w
case 2
- y5 I, g5 L- X' h- k% ^
P=MPQ(i,2);
* x4 j* Q! p1 u# R9 S/ p
b=MPQ(i,3);
6 r) d0 e3 K2 a$ m4 h8 s
RA=P*(L-b)/(L-L1);
! k# z" V" s1 m f
RB=P+P*(L-b)/(L-L1);
" b8 L) W2 H/ K0 V# J2 a
QQ=QQ+RA;
v) w& b+ }, H/ `" k
MM=MM+RA*x1;
9 x+ }# R& N% B
if b>0&&b<L
# G/ E e$ P. K2 L7 B
[QQ,MM]=QMP(n,x1,b,P,QQ,MM);
/ |& P5 z& @% {+ \6 s: R( }
end
% r- [$ D) V( d2 M4 C3 f
[QQ,MM]=QMP(n,x1,L-L1,RB,QQ,MM);
$ B% y! y$ e3 s" b# ]
case 3
* I/ g- c( ^9 ?; ]5 @
q=MPQ(i,2);
1 i) H1 N t3 [1 h c2 J3 b" M
c=MPQ(i,3);
3 q% X1 O+ N3 d2 ]4 |7 r0 [
d=MPQ(i,4);
6 ?0 P' `2 w( c% j
RA=q*(d-c)*(L-0.5*(d+c))/(L-L1);
! T. W4 \% p1 ^7 z8 }# A- n5 Q! a
RB=q*(d-c)-q*(d-c)*(L-0.5*(d+c))/(L-L1);
, _' F' L7 G- ~2 K1 Z" c
QQ=QQ+RA;
' u7 w4 E7 n! L, E8 p7 y& g1 n
MM=MM+RA*x1;
9 R8 ^3 v3 V' m& F7 d7 w
[QQ,MM]=QMQ(n,x1,c,d,q,QQ,MM);
6 Q+ r# D6 K+ U. Y
[QQ,MM]=QMP(n,x1,L-L1,RB,QQ,MM);
- X& i3 L; b/ [! C: i! h
end
: l' g6 L9 m7 F; B5 |' F
end
$ O: C* C0 }5 W6 C+ y* P c
[QDX,MDX,XQM]=MAX_MIN(x1,QQ,MM);
3 E/ A- I% \5 j
TU_QM(x1,QQ,MM);
1 R" {( {& n% `( l. I
第二子程序:
& W, f. B% T+ `( Q
function [QDX,MDX,XQM]=MAX_MIN(x1,QQ,MM)
1 r& \" N" S+ Y
XQM=[x1',QQ',MM'];
/ T# f8 i) i+ T0 S0 w
[Qmax,i]=max(QQ);
; M/ u$ G/ ^9 G7 }5 }8 ^4 _
Q1=[x1(i),Qmax];
9 s d4 C3 ]2 R, k3 d( a
[Qmin,i]=min(QQ);
. i: B. L9 a; S/ z7 ^2 c: }
Q2=[x1(i),Qmin];
+ K/ c5 R9 \# X' U; a
[Mmax,i]=max(MM);
1 k# z( Y! f L. w
M1=[x1(i),Mmax];
# v# w, c% `5 ]7 Z& @6 W
[Mmin,i]=min(MM);
9 j; M9 A* ]( i3 x
M2=[x1(i),Mmin];
; j% V+ ?& E1 d2 Z `' R d% T: @
QDX=[Q1;Q2]
6 R# _ D- T' ~ @8 d% a
MDX=[M1;M2]
) c. j( @: v2 X
第三子程序:
* j! W6 h, E7 ]8 Y
function MM=QMM(n,x1,a,M,MM)
# A, v2 X k, K, o' E
for j=1:n
1 u# y% \+ h5 l: P4 {+ O0 }6 C
if x1(j)==a
, G; o7 c0 `, l# `% P3 ?
n1=j;
3 `( p/ B" J& \7 u4 H# v; r4 E9 A
end
0 K+ W! n6 A' V3 n2 N+ Q" B- L
end
. R1 B6 v2 q' {1 Z4 K% y( I
MM(n1:n)=MM(n1:n)+M;
- {* ]6 l6 y: q' P+ G) d; s
第四子程序:
' J3 u0 c! d" A( r+ \& r& t
function [QQ,MM]=QMP(n,x1,b,P,QQ,MM)
- v" d5 M; M t
for j=1:n
4 o7 d6 @3 h) A3 }; R# n5 Z# T
if x1(j)==b
4 |" w: u: K* D+ H# K' d
n1=j;
1 T& U0 \ P. r; d
end
+ X* y7 @; O/ M2 i0 g; Y) e: X
end
. H* `5 X/ E: F8 _0 w( {
QQ(n1:n)=QQ(n1:n)-P;
: v9 T- Q" t, H: R
MM(n1:n)=MM(n1:n)-P*(x1(n1:n)-b);
6 b1 e; ]$ y' a# }7 C- u+ ?
第五子程序:
) X( I& B! |/ [1 f" C3 H: ^' ^2 M
function [QQ,MM]=QMQ(n,x1,c,d,q,QQ,MM)
) Q6 ~( t1 K1 {9 {7 t
for j=1:n
3 s: ^' g2 p& N. n& Q
if x1(j)>c
, Y9 j' v6 W6 M9 g: l
QQ(j)=QQ(j)-q*(x1(j)-c);
2 P; [( v0 ~9 a' J
MM(j)= MM(j)-0.5*q*(x1(j)-c)^2;
6 U- Z# r" L9 b5 T+ R
end
, u6 h7 [( p. Z) _4 K8 S' r) G
if x1(j)>d
3 T$ |# |' g% ]# A. \
QQ(j)=QQ(j)+q*(x1(j)-d);
, p4 y/ J+ q" D
MM(j)= MM(j)+0.5*q*(x1(j)-d)^2;
& H. X2 z- {" |* a. K/ Q7 r
end
* Z3 T/ _) x, p1 k9 g
end
: C+ w& m: a2 F
第六子程序:
. S! [( k& B5 d; R: T
function TU_QM(x1,QQ,MM)
; u6 j- U2 x/ h
subplot(2,1,1);
# m* g+ r$ [5 a) O; n
plot(x1,QQ);
5 Q1 l: o! E! |, C! ~, H- @* S
grid;
( a& R- X) [- u0 [5 y% t; m
title('剪力圖');
8 @7 {8 j- ?+ x" W W$ f5 n' n. D/ A
subplot(2,1,2);
' d8 L: i$ t) _
plot(x1,MM);
& r, i6 t* {) t# D4 c( A
grid;
y# T: q! O+ r: @; X2 V$ U
title('彎矩圖');
( h/ @6 `5 \- T5 [. [( N
8 M3 H( s6 b5 ~/ E) G4 o
歡迎光臨 機(jī)械社區(qū) (http://97307.cn/)
Powered by Discuz! X3.4