[求助]时变动力微分方程该如何求解？

xinquan · 发表于 2006-4-4 13:04

马上注册，结交更多好友，享用更多功能，让你轻松玩转社区。

您需要登录才可以下载或查看，没有账号？我要加入

x

假如一个多自由度结构，其质量，阻尼为常数，而刚度是随时间而变化的，请问一下，其动力方程的求解应该怎么实现？ 谢谢！十分着急！SOS!

siyanger · 发表于 2006-4-4 15:43

说清楚点吧，多少个自由度啊？我现在也正在解这类微分方程，15个自由度，用RK法没解出来呢。已换算法了！

xinquan · 发表于 2006-4-4 18:31

我现在只做简单的取三个自由度,如果找到方法就要根据实际情况改变自由度了,不知楼上的找到了什么好的方法啊? 还请指教啊!

siyanger · 发表于 2006-4-4 19:14

还在看资料哪！但算法一换，又要忙半天，还可能什么都出不来啊！我的是机械振动问题，自由度我觉得不算很多，不知道为什么算不出来。（是算出来，结果不对，也许和初值有关！）

happy · 发表于 2006-4-4 21:26

自由度比较少用ode45完全是可以的

siyanger · 发表于 2006-4-5 09:28

还有一说法，说微分方程组右端如果显含时间t，将时间t作为附加的状态变量，t=t,t'=1.这样相当于多加一个自由度。这知道这样是不是会好点？

xinquan · 发表于 2006-4-6 10:57

请问HAPPY: ode45可以用来解动力微分方程，那方程右边随时间变化的量该怎么处理呢？

yjh2911 · 发表于 2006-4-12 19:55

我做的是汽车振动方面的，6个自由度的二阶动力微分方程 能不能用RK法解阿？ 我用ODE45解总提示初始值少

happy · 发表于 2006-4-14 20:48

<DIV class=quote>以下是引用yjh2911在2006-4-12 19:55:07的发言： 
我做的是汽车振动方面的，6个自由度的二阶动力微分方程 能不能用RK法解阿？ 我用ODE45解总提示初始值少 </DIV>
最好把代码给出来，问题应该是你给的初值少于12个量

[此贴子已经被作者于2006-4-14 20:49:44编辑过]

xinquan · 发表于 2006-4-8 15:24

自己编的动力方程求解程序，结果能运行出来，但是得到的结果无论是有阻尼还是无阻尼，都是一个衰减信号,查不出问题啊，请高手帮忙看看! function vtb4(m,c,k,x0,v0,tf,w,f0,delt) fid1=fopen('e.txt','wt'); wn=sqrt(k/m); for t=0:delt:tf; xdd=(f0*sin(w*t)-k*x0-c*v0)/m; x3d=(f0*w*cos(w*t)-k*v0-c*xdd)/m; xd=v0+xdd*delt+x3d*delt^2/2; x=x0+xd*delt+xdd*delt^2/2; fprintf(fid1,'%10.4f',x); x0=x;v0=xd; t end fid2=fopen('e.txt','rt'); n=tf/delt; x=fscanf(fid2,'%f',[1,n]); t=1:n; plot(t,x),grid; xlabel('time(s)'); ylabel('displacement'); title('time and displacement relation'); 急啊!帮忙指点一下!不甚感激!

happy · 发表于 2006-4-9 15:47

把输入参数给出来

xinquan · 发表于 2006-4-14 15:16

输入参数如下:m=1e7,c=0,k=4e8,x0=1,v0=0,tf=10,w=2*pi,f0=1,delt=0.01 改变c得值，取c不等于0得时候，其输出图形基本没什么变化！

账号		自动登录	找回密码
密码			我要加入

[求助]时变动力微分方程该如何求解？

马上注册，结交更多好友，享用更多功能，让你轻松玩转社区。

回复：（xinquan）[求助]时变动力微分方程该如何求解？...

请问

回复：（yjh2911）请问

[求助]动力时程分析程序

回复：（xinquan）[求助]动力时程分析程序