数学建模社区-数学中国

标题: 用MAPLE解BVP遇到newton iteration not converging [打印本页]

---

作者: lambeef    时间: 2010-2-7 13:13
标题: 用MAPLE解BVP遇到newton iteration not converging
具体代码如下：
5 e' f2 H& q, z* a9 t" k3 Jl := 1000*10^(-9);
w := 50*10^(-9);
. {3 b9 j4 h: Bk := (1/12)*w^3;
& E7 C/ K. j+ y3 @* }" wE := 150*10^9;
2 k, |. M" }: Rd := 20*10^(-9);
+ z# S! t$ N) _+ D v := 2.3;
electrostatic := 8.854*10^(-12)*v^2*(1/2);
' ]8 ~" A2 J: V' Gode := diff(y(x), x$4) = electrostatic/((d-y(x))^2/E/k);
bcs1 := y(0) = 0, (D(y))(0) = 0, (D(D(y)))(l) = 0, (D(D(D(y))))(l) = 0;
; U" Y4 `/ F( K; Q- S) [ansl1 := dsolve({ode, bcs1}, numeric, method = bvp, abserr = 28)
问题是求在静电场作用下，梁的弯曲问题。v是电压,当电压很小时，可以算出结果，但是电压大于1V时出现‘newton iteration not converging’. 想了很长时间做不出来，请各位帮个忙！谢谢！

---

欢迎光临 数学建模社区-数学中国 (http://www.madio.net/)

Powered by Discuz! X2.5