数学建模社区-数学中国
标题:
求助：关于‘newton iteration not converging’问题
[打印本页]
作者:
lambeef
时间:
2010-2-7 13:17
标题:
求助：关于‘newton iteration not converging’问题
具体代码如下：
w8 m* S9 r0 L2 ?5 Y
l := 1000*10^(-9);
# a8 x5 x0 c1 o- W& `! P/ @3 |
w := 50*10^(-9);
r% u% F; z* `
k := (1/12)*w^3;
5 m# L- b( Q0 l; a7 e& s* `
E := 150*10^9;
9 S# B7 ?- l: r% m! ?0 d
d := 20*10^(-9);
6 m6 J: _- @- ]1 A9 z1 H% M
v := 2.3;
% J" x, v. }4 B P' E. m- y8 }: ^
electrostatic := 8.854*10^(-12)*v^2*(1/2);
7 `* U8 i8 N1 {; ^: R% M9 l. k) x( b
ode := diff(y(x), x$4) = electrostatic/((d-y(x))^2/E/k);
* ?2 Q! D6 V" V& k) `, s* \, E$ n
bcs1 := y(0) = 0, (D(y))(0) = 0, (D(D(y)))(l) = 0, (D(D(D(y))))(l) = 0;
1 }$ R3 u7 z: P. Y; z7 P
ansl1 := dsolve({ode, bcs1}, numeric, method = bvp, abserr = 28)
: t4 R% ?0 E9 e+ B$ K( N3 I
问题是求在静电场作用下，梁的弯曲问题。v是电压,当电压很小时，可以算出结果，但是电压大于1V时出现‘newton iteration not converging’. 想了很长时间做不出来，请各位帮个忙！谢谢！
作者:
eigen
时间:
2010-3-23 11:54
应该是迭代的初值没有选取好，或者这个问题本身不稳定。建议用Matlab软件再算一算。
欢迎光临 数学建模社区-数学中国 (http://www.madio.net/)
Powered by Discuz! X2.5