| 查看: 1562 | 回复: 3 | ||
zhuang1017铁虫 (小有名气)
|
[求助]
间隙碰撞振动matlab编程求助 已有1人参与
|
|
各位求助，关于两个自由度间隙碰撞的程序只能通过ode45算出没有发生碰撞时两个物体位移和速度，考虑间隙大小，发生碰撞后速度和位移都会发生变化，该怎样编程，自己编了一个程序，把两位移之差满足间隙作为条件语句，然后速度发生变化，但是位移任为未发生碰撞的位移，明显错误，现向各位求教，下面我贴出我的程序，请指正 function dy=test_fun(t,y) g(1)=0.02; g(2)=0.02; u=0.02; % f(t)=1; % fe=1; dy=zeros(4,1); % y=zeros(4,1); dy(1)=y(3); dy(2)=y(4); dy(3)= -2*g(1)*y(3)-y(1)+u*2*g(2)*y(4)+u*y(2)+cos(t);% c-b*x(1)-w*x(1)*x(2); dy(4)= -2*(1+u)*g(2)*y(4)-(1+u)*y(2)+2*g(1)*y(3)+y(1)-cos(t); % w*x(1)*x(2)-(b+m)*x(2); end y0=[0,0,0,0]; [t,y]=ode45('test_fun',[0,30],y0); T=30; y0=[0,0,0,0]; [t,y]=ode45('test_fun',T,y0); T=[0:0.1:30]; y0=[0,0,0,0]; [t,y]=ode45('test_fun',T,y0); u=0.02; r=0.2; for i=1:length(T) if abs(y(i,1)-y(i,2))>=0.8 Y(i,3)=(1-u*r)./(1+u)*y(i,3)+u*(1+r)./(1+u)*y(i,4); Y(i,4)=(1+r)./(1+u)*y(i,3)+(u-r)./(1+u)*y(i,4); else Y(i,3)=y(i,3); Y(i,4)=y(i,4); end Y(i,1)=y(i,1); Y(i,2)=y(i,2); end  碰撞.png  模型.jpg |
回复此楼» 猜你喜欢
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有6人回复
-
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有3人回复
-
为什么资助数各大高校都创新高，自己申请怎么就这么难
已经有15人回复
-
科研人应该花精力去思考如何解决问题，而不是去凝练问题
已经有14人回复
-
要骂人了，新模版改版就是要淡化问题凝练这种虚的东西，结果有个评委还在说凝练得不够
已经有17人回复
-
基金系统什么内容也没有
已经有10人回复
-
学科评审组评审是指会评吗？
已经有5人回复
-
面上合作单位盖章
已经有9人回复
-
面上函评意见出来了，像什么等级？
已经有17人回复
-
小白求助 投论文要求的highlights应该如何写
已经有4人回复
ftc_简单
银虫 (正式写手)
- 应助: 99 (初中生)
- 金币: 1276.3
- 散金: 151
- 红花: 29
- 帖子: 436
- 在线: 129.4小时
- 虫号: 4907455
- 注册: 2016-08-09
- 性别: GG
- 专业: 动力学与控制
2楼2017-06-08 15:36:23
zhuang1017
铁虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 758
- 散金: 50
- 帖子: 63
- 在线: 38.6小时
- 虫号: 4078324
- 注册: 2015-09-17
- 性别: GG
- 专业: 机械动力学
3楼2017-06-08 15:42:13
zhuang1017
铁虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 758
- 散金: 50
- 帖子: 63
- 在线: 38.6小时
- 虫号: 4078324
- 注册: 2015-09-17
- 性别: GG
- 专业: 机械动力学
4楼2017-06-08 15:42:31
30