数学建模社区-数学中国
标题:
用MAPLE解BVP遇到newton iteration not converging
[打印本页]
作者:
lambeef
时间:
2010-2-7 13:13
标题:
用MAPLE解BVP遇到newton iteration not converging
具体代码如下：
8 E# Z; R& v" e0 H Z/ E( g( C
l := 1000*10^(-9);
# [! Y& _9 j0 ~4 x
w := 50*10^(-9);
$ M6 b3 g. _: d7 `* n+ C
k := (1/12)*w^3;
, N, k1 I5 J+ `( _0 a0 N& b
E := 150*10^9;
( ~: O, i: {: @( M0 ^
d := 20*10^(-9);
$ d; J Z. |8 s4 M$ A& N( o) y. V
v := 2.3;
" L+ L- q8 ~1 S9 ?/ t) @, U: \
electrostatic := 8.854*10^(-12)*v^2*(1/2);
0 e% E, n& N; s. m3 d: G6 ~, j
ode := diff(y(x), x$4) = electrostatic/((d-y(x))^2/E/k);
# _ R3 q3 s6 B+ Q: V4 P- O
bcs1 := y(0) = 0, (D(y))(0) = 0, (D(D(y)))(l) = 0, (D(D(D(y))))(l) = 0;
4 `3 |6 u$ l" C3 F# L* Z* S$ C
ansl1 := dsolve({ode, bcs1}, numeric, method = bvp, abserr = 28)
: B6 q" p( L( |! ^" E9 I
问题是求在静电场作用下，梁的弯曲问题。v是电压,当电压很小时，可以算出结果，但是电压大于1V时出现‘newton iteration not converging’. 想了很长时间做不出来，请各位帮个忙！谢谢！
欢迎光临 数学建模社区-数学中国 (http://www.madio.net/)
Powered by Discuz! X2.5