QQ登录

只需要一步,快速开始

 注册地址  找回密码
查看: 3831|回复: 1
打印 上一主题 下一主题

求助:关于‘newton iteration not converging’问题

[复制链接]
字体大小: 正常 放大
lambeef 实名认证       

2

主题

0

听众

8

积分

升级  3.16%

该用户从未签到

自我介绍
200 字节以内

不支持自定义 Discuz! 代码
跳转到指定楼层
1#
发表于 2010-2-7 13:17 |只看该作者 |倒序浏览
|招呼Ta 关注Ta
具体代码如下:
9 E& X7 q) U+ R. jl := 1000*10^(-9); : F! S5 l' @4 z1 N. F
w := 50*10^(-9);
. j; f' P: h# l8 Dk := (1/12)*w^3; 1 _) K/ p( U% ]3 m( N& }: U
E := 150*10^9; - N  B! ?- }4 X' U, d
d := 20*10^(-9);
( k. i1 v2 m; k! p1 w# c v := 2.3; & ~8 z2 c- v- t* P8 Y! l# Y% t' e
electrostatic := 8.854*10^(-12)*v^2*(1/2);
1 X8 H) h4 {  S1 [* P; _$ u( code := diff(y(x), x$4) = electrostatic/((d-y(x))^2/E/k); ' J0 w" Y2 f  O$ y: Z9 v) G
bcs1 := y(0) = 0, (D(y))(0) = 0, (D(D(y)))(l) = 0, (D(D(D(y))))(l) = 0;
, X* ?9 W; ]9 Lansl1 := dsolve({ode, bcs1}, numeric, method = bvp, abserr = 28)6 L6 Q% C- B6 `3 S) I- e3 X8 S: ?
问题是求在静电场作用下,梁的弯曲问题。v是电压,当电压很小时,可以算出结果,但是电压大于1V时出现‘newton iteration not converging’. 想了很长时间做不出来,请各位帮个忙!谢谢!
zan
转播转播0 分享淘帖0 分享分享0 收藏收藏0 支持支持0 反对反对0 微信微信
eigen        

6

主题

4

听众

924

积分

升级  81%

  • TA的每日心情
    开心
    2021-11-8 14:12
  • 签到天数: 58 天

    [LV.5]常住居民I

    自我介绍
    eigen

    新人进步奖

    群组Mathematica研究小组

    群组数学建摸协会

    群组数学建模培训课堂1

    应该是迭代的初值没有选取好,或者这个问题本身不稳定。建议用Matlab软件再算一算。
    回复

    使用道具 举报

    您需要登录后才可以回帖 登录 | 注册地址

    qq
    收缩
    • 电话咨询

    • 04714969085
    fastpost

    关于我们| 联系我们| 诚征英才| 对外合作| 产品服务| QQ

    手机版|Archiver| |繁體中文 手机客户端  

    蒙公网安备 15010502000194号

    Powered by Discuz! X2.5   © 2001-2013 数学建模网-数学中国 ( 蒙ICP备14002410号-3 蒙BBS备-0002号 )     论坛法律顾问:王兆丰

    GMT+8, 2026-4-20 08:57 , Processed in 0.454517 second(s), 63 queries .

    回顶部