数学建模社区-数学中国
标题:
麻烦大家帮忙看下mathematical编的程序，不知道哪里出错了
[打印本页]
作者:
游龙游龙
时间:
2014-7-7 17:15
标题:
麻烦大家帮忙看下mathematical编的程序，不知道哪里出错了
本帖最后由 游龙游龙 于 2014-7-8 14:45 编辑
L" }1 z, h1 R
* Y& V) e6 x) y: R
麻烦大家帮忙看下这个程序，不知道哪里出错了，小弟是刚入门，请高手指教
; N6 O! w1 u, q2 D$ G
行列式求不出，画图找不到根
8 ]1 o; A7 t1 Y; f8 P) m( J* g! w
(*定义材料常数*)
: _0 s) U4 r; A% f% _
Subscript[\[Nu], TL] = 0.25;
- O6 V% X9 Q9 [; T5 C
Subscript[EE, L] = 53.8*10^9;(*L方向的弹性模量*)
3 j, ]& \% t7 u- p
Subscript[EE, T] = 17.93*10^9;(*T方向的弹性模量*)
8 @+ @& U4 f- K# K7 l; ~
Subscript[G, LT] = 8.96*10^9;(*LT方向的剪切模量*)
" E9 l/ w; E9 ?
Subscript[\[Nu], LT] =
2 |7 S( h) `9 |2 @
Subscript[\[Nu], TL]*Subscript[EE, L]/Subscript[EE, T];
+ } R; K1 T) V" R, F$ {
\[Mu] = 1 - Subscript[\[Nu], LT]*Subscript[\[Nu], TL];
& U" Z. [" f, c3 ] b4 J
Subscript[C, L] = Subscript[EE, L]/\[Mu];
* {1 U& X( C0 U. R Q0 o: k
Subscript[C, S] = Subscript[\[Nu], LT]*Subscript[EE, T]/\[Mu];
- j" O' M8 K) t. X" X0 o4 J8 H9 W
Subscript[C, T] = Subscript[EE, T]/\[Mu];
6 d" E8 i, a6 p$ F% m
Subscript[C, LT] = Subscript[G, LT];
* [$ t. x4 l9 R3 k
(*脱层几何参数*)
! S7 K: U. Y$ t: Z; b
\[Rho] = 1600;(*脱层密度*)
m9 g$ i' n2 P6 \% T* [( t: M
c = 0.06;(*阻尼*)
* s; o; P( v9 O/ M& n* m6 H
t = 0.01;(*t每层厚度*)
8 I; G& V0 R/ ^: z! g5 [, A6 L9 `
h = 0.05;(*h为梁总厚度*)
! q1 j: |- x, ~" _- T: c
l = 1;(*为梁总的长度*)
: s: d" J. ]- x, D* `* x. g+ H
Subscript[l, 1] = 0.3;(*为梁1区的长度*)
4 S! k6 D" s: z
Subscript[l, 2] = 0.3;(*为梁2区的长度*)
' U G( o8 p! y1 x/ e
Subscript[l, 3] = 0.4;(*为梁3区的长度*)
+ e. t9 V! n* `$ r0 S6 S B' ?! t
Subscript[l, 4] = 0.4;(*为梁4区的长度*)
) G) f6 k& I& r* m' V' u2 g
nn1 = h/t;(*nn1为总铺层厚度*)
- v7 A7 _# o' z* P) b) u H
nn2 = 2.0;(*为上子铺层层数*)
5 ~: Y/ _+ b5 X- |5 d7 X4 O. U
nn3 = h/t - nn2;(*下子铺层层数*)
: e* B# O- o0 x9 D; Y
Subscript[h, 3] = t*nn2;(*梁的上部子厚度*)
8 _3 C3 @8 y) p3 C8 f# V( `
Subscript[h, 4] = h - Subscript[h, 3];(*梁的下部子厚度*)
0 r6 C, X l: q3 c! Y2 P
Subscript[h, 1] = h;(*梁1区的厚度*)
8 I% Y6 |, g" F7 T& M
Subscript[h, 2] = h;(*梁2区的厚度*)
0 X3 a. X7 Q9 `6 C2 K3 ]& G
(*求各单层刚度*)
8 X* e1 @& v6 B- m3 ]3 V3 K- i
agk = {0 Pi/180, 90 Pi/180, 0 Pi/180, 90 Pi/180, 0 Pi/180};(*列表各铺层角度*)
5 Q6 s! s2 a: R2 W: `* a6 \
" r$ G- S( M( `* ]1 D2 q" I9 y& W
5 H$ B0 Z4 `- x0 f" W! S+ z
ck = Cos[agk];(*列表求各铺层角度余弦*)
( y- W6 H% D; i# D+ j! [2 q
sk = Sin[agk];(*列表求各铺层角度正弦*)
2 [: [) g {; f0 Q8 Z" i, }. P2 f
s11k = ck^4*Subscript[C, L] + 2 ck^2*sk^2*Subscript[C, S] +
, Q/ w/ ~; L& w
sk^4*Subscript[C, T] +
! `' L& \. {8 V4 F9 r8 I* x
4 ck^2*sk^2*Subscript[C, LT];(*求各区刚度*)
% w/ B" |) Z H; m3 O
A11 = \!\(
- w0 j, ^2 {! s2 x P$ x
\*UnderoverscriptBox[\(\[Sum]\), \(k =
; y+ `% C1 J8 h/ f$ Q4 W; u
1\), \(nn1\)]\((s11k[\([\)\(k\)\(]\)]*t)\)\)(*一区拉伸刚度*)
. S# L1 n5 c4 m4 f* j {
A21 = A11
/ ~4 Q1 m' w5 ?$ t" x1 `
A31 = (\!\(
0 ]! P/ W% B* L# x/ [' ~
\*UnderoverscriptBox[\(\[Sum]\), \(k =
( z6 Z* R$ L d8 V
1\), \(nn2\)]\(\((s11k[\([\)\(k\)\(]\)])\)*t\)\))(*三区拉伸刚度*)
& l) V0 h) k$ O* `
A41 = (\!\(
. ~) Z3 i. f) V1 j2 X; ?
\*UnderoverscriptBox[\(\[Sum]\), \(k =
$ ~6 P/ |" D4 }3 [9 a0 b
nn2 + 1\), \(nn1\)]\(\((s11k[\([\)\(k\)\(]\)])\)*t\)\))(*四区拉伸刚度*)
9 r, G9 h- n5 G9 z6 L. f
1 q* [/ K+ l! o0 }
; F; R# j* Z; e, [" T
Print["A21=A11;D21=D11"];(*二区拉伸刚度、弯曲刚度*)
. q/ M+ d- A% \+ E' `/ I, z+ l8 r
D11 = \!\(
7 @' R' c% k, b m3 T+ w- Y4 {
\*UnderoverscriptBox[\(\[Sum]\), \(k =
~5 Z) y2 M) Z! j
1\), \(nn1\)]\((s11k[\([\)\(k\)\(]\)]*
1 E, |7 B, R9 Z% c O
\*FractionBox[\(t^3\), \(12\)])\)\)
2 U! T0 O" V, A8 S( ^) a; d' A, f
D21 = D11
* k( p, P4 L* k, s C
D31 = \!\(
* m3 w2 S& f; d0 @( C. O1 @& F
\*UnderoverscriptBox[\(\[Sum]\), \(k =
+ [- z) |0 q7 a% C5 y3 l
1\), \(nn2\)]\((s11k[\([\)\(k\)\(]\)]*
' O2 u6 B' g2 y5 {8 L0 H
\*FractionBox[\(t^3\), \(12\)])\)\)(*三区弯曲刚度*)
0 _2 Q/ T, ]5 p7 y- R. K( P. o6 o( \
D41 = \!\(
: j' z0 m3 D& X
\*UnderoverscriptBox[\(\[Sum]\), \(k =
* v! q s# h2 i$ t. V! f
nn2 + 1\), \(nn1\)]\((s11k[\([\)\(k\)\(]\)]*
0 t( u! q, j5 G& L V4 j4 U
\*FractionBox[\(t^3\), \(12\)])\)\)(*四区弯曲刚度*)
7 U. }& U2 x3 [ o- t
k1 = Subscript[l, 1]/l; k2 = Subscript[l, 2]/l; k3 = Subscript[l, \
4 T. {' n5 ~/ u! i5 w# r% ^9 t) G
3]/l; k4 = Subscript[l, 4]/l; \[Alpha] = h/l;
g2 p7 U4 L- ] r
a1 = (\[Alpha]*k1)/(DD1*(\[Beta]1)^3); a2 = (\[Alpha]*k2)/(
, f9 y( E8 d; m$ d% E# f4 d/ w
DD2*(\[Beta]2)^3); a3 = (\[Alpha]*k3)/(
" W8 P' }' q$ U0 j2 o7 g
DD3*(\[Beta]3)^3); a4 = (\[Alpha]*k4)/(DD4*(\[Beta]4)^3);
* h4 L/ \: ]6 G
\[Beta]1 = h/Subscript[l, 1]; \[Beta]2 = h/Subscript[l, 2]; \[Beta]3 \
* f8 `0 K/ }; Z* I6 H# V- ^* q: p
= h/Subscript[l, 3]; \[Beta]4 = h/Subscript[l, 4];
& ], X/ z1 P/ \2 Z
b1 = Subscript[l, 1]/Subscript[l, 3]; b2 = Subscript[l, \
3 X. z4 R. T' J& a
1]/Subscript[l, 4]; b3 = Subscript[l, 2]/Subscript[l, 3]; b4 = \
* P- K+ E3 Q- R3 T- ]+ Z" `! r, F
Subscript[l, 2]/Subscript[l, 4];
& x0 \1 {6 H. {, U- a8 |. @, F# y
DD1 = D11/(A11*h^2); DD2 = D21/(A11*h^2); DD3 = D31/(A11*h^2);
6 W7 \- c/ {% U( v6 u) e
DD4 = D41/(A11*h^2);
9 c, N5 ?! Q% F8 R; @7 b. j
r11 = Power[a1*\[Omega]^2, (4)^-1]; r12 = Power[
1 K- i3 D6 P# Z& J6 J- T7 K. i% y
a2*\[Omega]^2, (4)^-1]; r13 = Power[
7 a" Z; L! A$ E" [* [' r: e
a3*\[Omega]^2, (4)^-1]; r14 = Power[a4*\[Omega]^2, (4)^-1];
& ~1 c/ P9 `5 P' D% M
BB = {{1, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0},
5 q' X" ~) P: e
{0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0},
# y6 i7 }8 b* E% m4 O
{0, Cosh[r12], 0, 0, 0, Sinh[r12], 0, 0, 0, Cos[r12], 0, 0, 0,
% t6 O7 `! a4 y+ p& i+ H3 }
Sin[r12], 0, 0},
( E0 b( p7 P+ `
{0, Sinh[r12], 0, 0, 0, Cosh[r12], 0, 0, 0, -Sin[r12], 0, 0, 0,
7 F% P3 H( V* L8 k3 |8 R8 \
Cos[r12], 0, 0},
% f: O" E; U3 J2 Q$ F* x
{-(r11)^2*DD1*(\[Beta]1)^2*Cosh[r11],
: a3 L/ b( s$ `& g, Y
0, (r13)^2*DD3*(\[Beta]3)^2, (r14)^2*DD4*(\[Beta]4)^2,
+ y- s# k) ?$ f6 t4 h
-(r11)^2*DD1*(\[Beta]1)^2*Sinh[r11], 0, 0,
3 n# _" n* l, Z, c, z2 B0 A
0, (r11)^2*DD1*(\[Beta]1)^2*Cos[r11],
t- ~" n j/ B, i/ D ^+ l4 i% u
0, -(r13)^2*DD3*(\[Beta]3)^2, -(r14)^2*DD4*(\[Beta]4)^2, (r11)^2*
$ g3 ~' D) W8 l) ~2 {) l" A# X
DD1*(\[Beta]1)^2*Sin[r11], 0, 0, 0},
1 ? L$ k7 m% A) i- q7 A% m( i( K
{-(r11)^3*DD1*(\[Beta]1)^3*Sinh[r11], 0, 0,
& e# z/ D) f4 `) T9 e& m7 j# m
0, -(r11)^3*DD1*(\[Beta]1)^3*Cosh[r11],
* \ u; C) P% r' l5 H0 }
0, (r13)^3*DD3*(\[Beta]3)^3, (r14)^3*DD4*(\[Beta]4)^3, -(r11)^3*
& j+ [: o V" T8 ?% B
DD1*(\[Beta]1)^3*Sin[r11], 0, 0, 0,
4 ?- s2 }! c0 c. q) P! [8 e- k6 D& m
(r11)^3*DD1*(\[Beta]1)^3*Cos[r11],
6 x- [: S9 \' v! r. ~! Z
0, -(r13)^3*DD3*(\[Beta]3)^3, -(r14)^3*DD4*(\[Beta]4)^3},
2 v$ O6 E( {+ T/ f( U( R5 s# U
{0, -(r12)^2*DD2*(\[Beta]2)^2, (r13)^2*DD3*(\[Beta]3)^2*
1 t) a" M3 E1 q8 p \( I: a
Cosh[r13], (r14)^2*DD4*(\[Beta]4)^2*Cosh[r14],
8 F# @$ R& S: F; s7 _! w0 _
0, 0, (r13)^2*DD3*(\[Beta]3)^2*Sinh[r13], (r14)^2*DD4*(\[Beta]4)^2*
0 w7 ^4 m6 c+ }+ o1 R
Sinh[r14], 0, (r12)^2*DD2*(\[Beta]2)^2,
$ q/ W L1 h$ ?: A+ i [; U; Q
-(r13)^2*DD3*(\[Beta]3)^2*Cos[r13], -(r14)^2*DD4*(\[Beta]4)^2*
3 J/ j9 z) T1 {8 P/ D) Q4 z
Cos[r14], 0, 0,
* Q% Y" i/ R0 G
-(r13)^2*DD3*(\[Beta]3)^2*Sin[r13], -(r14)^2*DD4*(\[Beta]4)^2*
( ]/ K R9 k+ |4 Y7 A4 i
Sin[r14]},
' i* w1 N" _8 e- S! O, {0 H8 i `
{0, 0, (r13)^3*DD3*(\[Beta]3)^3*Sinh[r13], (r14)^3*DD4*(\[Beta]4)^3*
+ e. V- C2 W2 q) _; W7 Q+ K! R
Sinh[r14], 0, -(r12)^3*DD2*(\[Beta]2)^3,
" f& I' z" ~: ]4 R( A6 x s. s: I" L
(r13)^3*DD3*(\[Beta]3)^3*Cosh[r13], (r14)^3*DD4*(\[Beta]4)^3*
( N. G0 x, b0 O2 k
Cosh[r14], 0, 0, (r13)^3*DD3*(\[Beta]3)^3*Sin[r13],
/ c3 Y: Y! E6 ~$ b$ t& P$ L% @
(r14)^3*DD4*(\[Beta]4)^3*Sin[r14],
0 u# e2 p9 R" s
0, (r12)^3*DD2*(\[Beta]2)^3, -(r13)^3*DD3*(\[Beta]3)^3*Cos[r13],
5 y% B4 y+ M2 E
-(r14)^3*DD4*(\[Beta]4)^3*Cos[r14]},
( a6 K6 b" E+ X& c. @
{Cosh[r11], 0, -1, 0, Sinh[r11], 0, 0, 0, Cos[r11], 0, -1, 0,
c% u5 r! `# Z
Sin[r11], 0, 0, 0},
4 X$ \# E* n7 R/ {! g
{r11*Sinh[r11], 0, 0, 0, r11*Cosh[r11], 0, -r13*b1,
2 p) Q4 Y0 \5 Q2 O- E" ?% D
0, -r11*Sin[r11], 0, 0, 0, r11*Cos[r11],
9 A- o+ ~! K/ N0 \
0, -r13*b1},
3 d+ C# m- O) h! ^
{Cosh[r11], 0, 0, -1, Sinh[r11], 0, 0, 0, Cos[r11], 0, 0, -1,
: o, }; a3 H6 y" K# ?! a# E: d
Sin[r11], 0, 0, 0},
- @' o/ r3 u9 r1 G( t0 w% g+ W
{r11*Sinh[r11], 0, 0, 0, r11*Cosh[r11], 0,
* j! o5 h! f* N7 X- B
0, -r14*b2, -r11*Sin[r11], 0, 0, 0, r11*Cos[r11],
4 V3 J4 P; s& H4 I8 q8 O5 g
0, 0, -r14*b2},
* A: P% O, ^, l
{0, 1, -Cosh[r13], 0, 0, 0, -Sinh[r13], 0, 0, 1, -Cos[r13], 0, 0, 0,
% \: g' n3 X4 j+ [' _: y c& \
0, -Sin[r13]},
; S# k; y. {8 }
{0, 0, -r13*b3*Sinh[r13], 0, 0, r12, -r13*b3*Cosh[r13], 0, 0, 0,
: ^- e' U6 f( G0 H
r13*b3*Sin[r13],
$ u8 m) d l0 \3 S* j
0, 0, r12, -r13*b3*Cos[r13], 0},
+ N$ n! ^- _* O' y( r4 G& N. }
{0, 1, 0, -Cosh[r14], 0, 0, 0, -Sinh[r14], 0, 1, 0, -Cos[r14], 0, 0,
4 u8 [: A+ p H, K
0, -Sin[r14]},
! w2 h6 W5 ^+ i6 M$ t6 S; h! D
{0, 0, 0, -r14*b4*Sinh[r14], 0, r12, 0, -r14*b4*Cosh[r14], 0, 0, 0,
' }8 H( `( |& S: c; Y
r14*b4*Sin[r14],
& u$ Y; Z; X& H- x0 d) `7 m' E
0, r12, 0, -r14*b4*Cos[r14]}}
0 {# p, E9 q0 y N2 I; r* J
aa1 = Det[BB]
1 I+ e# c% D1 r) Q
Plot[aa1, {\[Omega], -10, 10}]
( T2 U7 H0 _* J: g$ c# q$ {) W
ss = FindRoot[aa1 == 0, {\[Omega], 2}][[1]]
' p* O: ~; _6 {9 o9 n( L) r0 p
/ j2 f; ^) ], g. M1 z1 D
作者:
百年孤独
时间:
2014-7-8 09:30
这是什么语言编程？
5 |3 H: A* z- l& L% x; l, ~ d" y
作者:
游龙游龙
时间:
2014-7-8 14:44
百年孤独 发表于 2014-7-8 09:30
2 L/ }; t j" A5 d- ~3 C8 [2 a
这是什么语言编程？
( q% c* a) ~) |) G z4 Z6 y
mathematical编的
欢迎光临 数学建模社区-数学中国 (http://www.madio.net/)
Powered by Discuz! X2.5
