QQ登录

只需要一步,快速开始

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

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

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

2

主题

0

听众

8

积分

升级  3.16%

该用户从未签到

自我介绍
200 字节以内

不支持自定义 Discuz! 代码
跳转到指定楼层
1#
发表于 2010-2-7 13:17 |只看该作者 |倒序浏览
|招呼Ta 关注Ta
具体代码如下:3 ]/ f  N6 A& z/ f" M+ a
l := 1000*10^(-9);
- D0 q: f% T- r6 Lw := 50*10^(-9);
# J. L8 M& h; zk := (1/12)*w^3;
& D9 Y3 z, |& o: ]. W8 _3 VE := 150*10^9; 0 j7 a. {4 u. {# ~+ O8 J
d := 20*10^(-9);
; n) ^; ?3 E) k7 w v := 2.3; % Z, t: t/ V+ }( t6 ?9 M) F' F
electrostatic := 8.854*10^(-12)*v^2*(1/2);
, b  I0 C/ H+ i+ \3 xode := diff(y(x), x$4) = electrostatic/((d-y(x))^2/E/k);
- \) j: d; n- W0 ubcs1 := y(0) = 0, (D(y))(0) = 0, (D(D(y)))(l) = 0, (D(D(D(y))))(l) = 0;
. d5 B0 ?: y7 V# s* G- N0 xansl1 := dsolve({ode, bcs1}, numeric, method = bvp, abserr = 28)
/ q5 j# u, ?$ ^6 Q; I; m0 M( N问题是求在静电场作用下,梁的弯曲问题。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-19 21:53 , Processed in 0.448524 second(s), 63 queries .

    回顶部