- 在线时间
- 0 小时
- 最后登录
- 2010-2-8
- 注册时间
- 2010-2-7
- 听众数
- 0
- 收听数
- 0
- 能力
- 0 分
- 体力
- 22 点
- 威望
- 0 点
- 阅读权限
- 20
- 积分
- 8
- 相册
- 0
- 日志
- 0
- 记录
- 0
- 帖子
- 3
- 主题
- 2
- 精华
- 0
- 分享
- 0
- 好友
- 0
升级   3.16% 该用户从未签到 - 自我介绍
- 200 字节以内
不支持自定义 Discuz! 代码
 |
具体代码如下:
0 P5 I- Z1 n" s9 rl := 1000*10^(-9); . Y8 a# w: t! H. F- e1 U" T0 j
w := 50*10^(-9);
. O: I- R0 b2 I9 G2 L* z: Bk := (1/12)*w^3; 1 x, Y: x. s# N# B
E := 150*10^9;
: P6 A+ k2 \7 wd := 20*10^(-9);
7 D* ~4 z: `: s% I v := 2.3;
# W8 l% j; V1 g3 ielectrostatic := 8.854*10^(-12)*v^2*(1/2); ( p g7 E2 g( }$ W
ode := diff(y(x), x$4) = electrostatic/((d-y(x))^2/E/k);
6 i2 D: H; O1 j" R* Q5 `, Tbcs1 := y(0) = 0, (D(y))(0) = 0, (D(D(y)))(l) = 0, (D(D(D(y))))(l) = 0; / v% g, y: s v: i$ t
ansl1 := dsolve({ode, bcs1}, numeric, method = bvp, abserr = 28)
( H+ n/ Y% B5 z1 f问题是求在静电场作用下,梁的弯曲问题。v是电压,当电压很小时,可以算出结果,但是电压大于1V时出现‘newton iteration not converging’. 想了很长时间做不出来,请各位帮个忙!谢谢! |
zan
|