| 查看: 511 | 回复: 0 | ||
[求助]
行星齿轮均载特性研究matlab编程
|
|
有会的前辈可以帮我看下们，为啥解不出来，公式没问题，编程估计不对，刚刚学matlab t=0:0.001:0.05; rs=19*8.4667; w=2; ka=0.25; km=2.5*1.0e8; rp=33*8.4667; E=10e-6; a1=19/(1+(86/19))*t+pi/2-pi/9; a2=19/(1+(86/19))*t+pi/2-pi/9+2*pi/3; a3=19/(1+(86/19))*t+pi/2-pi/9+4*pi/3; b1=19/(1+(86/19))*t+pi/2+pi/9; b2=19/(1+(86/19))*t+pi/2+pi/9+2*pi/3; b3=19/(1+(86/19))*t+pi/2+pi/9+4*pi/3; syms x1 y1 x2 y2 fs1 fs2 fs3 fr1 fr2 fr3 v v1 v2 v3 S=solve([(fs1.*cos(a1)+fs2.*cos(a2)+fs3.*cos(a3))-2.*2.7530.*10^9.*x1==0,(fs1.*sin(a1)+fs2.*sin(a2)+fs3.*sin(a3))-2.*2.7530.*10^9.*y1==0,(fr1.*cos(b1)+fr2.*cos(b2)+fr3.*cos(b3))-2.*3.7061.*10^9.*x2==0,(fr1.*sin(b1)+fr2.*sin(b2)+fr3.*sin(b3))-2.*3.7061.*10^9.*y2==0,fs1==fr1,fs2==fr2,fs3==fr3,41379-rs.*(fs1+fs2+fs3)==0,fr1==km.*(1+ka.*cos(w.*t)).*(rp.*v1+E.*cos(w.*t)),fr2==km.*(1+ka.*cos(w.*t)).*(rp.*v2+E.*cos(w.*t)),fr3==km.*(1+ka.*cos(w.*t)).*(rp.*v3+E.*cos(w.*t)),fs1==km.*(1+ka.*cos(w.*t)).*(rs.*v-rp.*v1+E.*cos(w.*t)),fs2==km.*(1+ka.*cos(w.*t)).*(rs.*v-rp.*v2+E.*cos(w.*t)),fs3==km.*(1+ka.*cos(w.*t)).*(rs.*v-rp.*v3+E.*cos(w.*t))],[x1,y1,x2,y2,v,v1,v2,v3,fs1,fs2,fs3,fr1,fr2,fr3]); @ctgu_zheng 发自小木虫IOS客户端 |
回复此楼» 猜你喜欢
北京理工大学 空天学院-招收发动机相关专业（力学、燃烧传热等方向）博士生
已经有2人回复
-
陶瓷销盘实验销材料求助
已经有2人回复
-
机械工程论文润色/翻译怎么收费?
已经有166人回复
-
C/C-SiC复合材料材料性能问题
已经有5人回复
-
定向工程博士如何申请
已经有2人回复
-
紧急招收2026年秋季入学博士生1名（湘潭大学 固体废弃物低碳利用湖南省工程研究中心）
已经有20人回复
-
【急招】“双一流”高校-新能源材料课题组招收2026年秋季入学博士生1名(湘潭大学)
已经有22人回复
-
祈祷面上基金能中
已经有129人回复
-
力学研究所某专项团队支撑岗位和特别研究助理岗位招聘启事(长期有效)
已经有10人回复
-
力学研究所某专项团队支撑岗位和特别研究助理岗位招聘启事
已经有10人回复
-
自用分享｜桌面深层文件夹分组映射工具 HyperDesk 免费体验版
已经有2人回复