数学建模社区-数学中国
标题:
用MAPLE解BVP遇到newton iteration not converging
[打印本页]
作者:
lambeef
时间:
2010-2-7 13:13
标题:
用MAPLE解BVP遇到newton iteration not converging
具体代码如下：
+ ], q+ a# H% B
l := 1000*10^(-9);
4 e, ?! ~* R/ Z# F
w := 50*10^(-9);
1 g+ g, D% M7 A3 V# \! x
k := (1/12)*w^3;
0 C; i5 {2 p' B& S
E := 150*10^9;
, E" `0 _1 g& E8 B
d := 20*10^(-9);
' F2 t8 I x8 R7 j1 F7 \7 K; Z% V
v := 2.3;
- o7 M8 v) ?% }% a3 ~1 u% @8 r
electrostatic := 8.854*10^(-12)*v^2*(1/2);
: B% v( I, s* R0 k9 _# k- W9 w
ode := diff(y(x), x$4) = electrostatic/((d-y(x))^2/E/k);
j- }9 ]2 x( _2 j w0 P5 ~( D
bcs1 := y(0) = 0, (D(y))(0) = 0, (D(D(y)))(l) = 0, (D(D(D(y))))(l) = 0;
) x, }+ l R% a2 q) p
ansl1 := dsolve({ode, bcs1}, numeric, method = bvp, abserr = 28)
# _. E+ c/ q3 A8 m" j7 x2 _
问题是求在静电场作用下，梁的弯曲问题。v是电压,当电压很小时，可以算出结果，但是电压大于1V时出现‘newton iteration not converging’. 想了很长时间做不出来，请各位帮个忙！谢谢！
欢迎光临 数学建模社区-数学中国 (http://www.madio.net/)
Powered by Discuz! X2.5