数学建模社区-数学中国
标题:
用MAPLE解BVP遇到newton iteration not converging
[打印本页]
作者:
lambeef
时间:
2010-2-7 13:13
标题:
用MAPLE解BVP遇到newton iteration not converging
具体代码如下：
9 o% M p+ _9 E* S8 M+ _# k
l := 1000*10^(-9);
2 u, Q8 b' f4 O- `* ` U) q$ U, J: V
w := 50*10^(-9);
0 P. N; p0 @ {, y8 v
k := (1/12)*w^3;
, g( e. e$ h" p+ q1 q6 }
E := 150*10^9;
7 c# j: Z5 C% v* f: g
d := 20*10^(-9);
( X$ g- t* w/ Y1 }& K N; I1 z
v := 2.3;
6 [6 c2 ~5 }' O. {3 F
electrostatic := 8.854*10^(-12)*v^2*(1/2);
Q9 I$ g$ T) ?% c4 {- F
ode := diff(y(x), x$4) = electrostatic/((d-y(x))^2/E/k);
: v; t3 G! _# @" R9 ^
bcs1 := y(0) = 0, (D(y))(0) = 0, (D(D(y)))(l) = 0, (D(D(D(y))))(l) = 0;
# B" y$ Q# G' `4 b8 D
ansl1 := dsolve({ode, bcs1}, numeric, method = bvp, abserr = 28)
0 d1 T2 l2 ~1 d
问题是求在静电场作用下，梁的弯曲问题。v是电压,当电压很小时，可以算出结果，但是电压大于1V时出现‘newton iteration not converging’. 想了很长时间做不出来，请各位帮个忙！谢谢！
欢迎光临 数学建模社区-数学中国 (http://www.madio.net/)
Powered by Discuz! X2.5