| 查看: 1001 | 回复: 6 | |||
| 本帖产生 1 个 机械EPI ，点击这里进行查看 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
[交流]
Matlab编程 已有1人参与
|
|||
a1=0:0.05:0.5;a=-a1:0.1:0.0826;请用matlab求关于a和kb的曲线，谢谢大神 4(]Q1R(C~B}URL4YH212T~B.jpg |
回复此楼» 猜你喜欢
TC4钛合金激光焊接金相组织分析
已经有2人回复
-
读博求助-本硕专业和工作行业不一样，还可不可以读博，可以的话，读哪个呢
已经有3人回复
-
金属材料论文润色/翻译怎么收费?
已经有221人回复
-
26申博-西工大材料专业-电化学方向-情况有点特殊
已经有6人回复
-
AI赋能复合材料力学：机器学习、PINN与多尺度仿真实战
已经有1人回复
-
合金高温性能
已经有0人回复
-
AI+ CFD！20 个实战案例拆解 fluent、pinn 与智能流体应用
已经有1人回复
-
为什么复合材料+机器学习那么好发文章？
已经有0人回复
-
纳米颗粒分散度差，容易团聚怎么办
已经有1人回复
-
求一篇论文，出自Molecular Simulation，
已经有0人回复
送红花一朵 5楼2013-07-17 19:46:48
Remantic
木虫 (正式写手)
元始天尊
- 机械EPI: 1
- 应助: 19 (小学生)
- 贵宾: 0.01
- 金币: 4190.6
- 散金: 1002
- 红花: 3
- 沙发: 1
- 帖子: 828
- 在线: 1134.9小时
- 虫号: 966445
- 注册: 2010-03-10
- 专业: 机器人学及机器人技术
2楼2013-07-16 22:01:20
3楼2013-07-17 17:33:53
Remantic
木虫 (正式写手)
元始天尊
- 机械EPI: 1
- 应助: 19 (小学生)
- 贵宾: 0.01
- 金币: 4190.6
- 散金: 1002
- 红花: 3
- 沙发: 1
- 帖子: 828
- 在线: 1134.9小时
- 虫号: 966445
- 注册: 2010-03-10
- 专业: 机器人学及机器人技术
★ ★ ★ ★ ★ ★
小木虫: 金币+0.5, 给个红包，谢谢回帖
茉莉新西兰: 金币+5, 机械EPI+1, 非常专业的应助，辛苦了，奖励EPI一枚 2013-07-18 22:43:20
小木虫: 金币+0.5, 给个红包，谢谢回帖
茉莉新西兰: 金币+5, 机械EPI+1, 非常专业的应助，辛苦了，奖励EPI一枚 2013-07-18 22:43:20
|
实际上应是a1与kb的关系曲线，因为kb只是a1的函数，属于变下限积分问题，偶编了个MATLAB程序，你看看对不对： % 小木虫 % http://muchong.com/bbs/viewthread.php?tid=6130700&pid=1#pid1 % date 2013-07-18 clc;clear; kb = zeros(size(0:0.05:0.5)); a2 = 0.0826; E =2.1e11;L=90; i=0; for a1 =0:0.05:0.5; i=i+1; FH = @(a)3*(1+cos(a1).*((a2-a)).*sin(a)-cos(a)).^2.*(a2-a).*cos(a)./(2*E*L*(sin(a)+(a2-a).*cos(a)).^3); kb(i) = quadl(FH,-a1,a2,1e-6); end plot(0:0.05:0.5,1./kb,'r','linewidth',3); xlabel('\alpha ');ylabel('Kb');  untitled.jpg |
» 本帖已获得的红花（最新10朵）
4楼2013-07-17 19:29:39
