数学建模社区-数学中国
标题:
求助：关于‘newton iteration not converging’问题
[打印本页]
作者:
lambeef
时间:
2010-2-7 13:17
标题:
求助：关于‘newton iteration not converging’问题
具体代码如下：
& F4 A* w1 m- n, a1 N
l := 1000*10^(-9);
3 c* F- h# \+ ^3 P9 [" K
w := 50*10^(-9);
/ ~: {7 y2 g5 V- d
k := (1/12)*w^3;
( E6 ^' Y' o: M& ~; f
E := 150*10^9;
; g4 M* U9 R& ^3 B* j
d := 20*10^(-9);
6 i. h4 P5 y% j1 f& N5 m( E* q
v := 2.3;
2 o& P: a& X/ _# x) `7 K2 w. h% `3 l% \
electrostatic := 8.854*10^(-12)*v^2*(1/2);
! e3 y% w6 R5 e+ s
ode := diff(y(x), x$4) = electrostatic/((d-y(x))^2/E/k);
( l7 |6 P( e' G, V; ]7 j
bcs1 := y(0) = 0, (D(y))(0) = 0, (D(D(y)))(l) = 0, (D(D(D(y))))(l) = 0;
* i* s1 A' I) p7 `9 T& e5 F4 y
ansl1 := dsolve({ode, bcs1}, numeric, method = bvp, abserr = 28)
9 F0 i- y p$ H p K4 f
问题是求在静电场作用下，梁的弯曲问题。v是电压,当电压很小时，可以算出结果，但是电压大于1V时出现‘newton iteration not converging’. 想了很长时间做不出来，请各位帮个忙！谢谢！
作者:
eigen
时间:
2010-3-23 11:54
应该是迭代的初值没有选取好，或者这个问题本身不稳定。建议用Matlab软件再算一算。
欢迎光临 数学建模社区-数学中国 (http://www.madio.net/)
Powered by Discuz! X2.5