tiange 发表于 2006-10-2 11:24

这个方程怎么编程序(关于悬臂梁受到一冲击)

在得到悬臂梁的运动方程,利用振型函数展开后,
如果给悬臂梁施加一个冲击力,怎么用MATLAB 编写程序求其广义坐标Y(t)?

dsolve('D2y+2*m*n*Dy+n*n*y=f(t)’)
但是这个f(t)是一个冲击函数 只是在t=0时刻的一个冲击力
在t=0时刻为一值,而在其他时刻均为零
怎么体现在函数里面的

笑石头 发表于 2006-10-2 20:18

在simulink下达建模型很简单的就解出来了。D2y+2*m*n*Dy+n*n*y=f(t)经拉普拉斯变换(假设初始条件为零)得
s^2*y+2*m*n*s*y+n^2*y=f(t),上式可以变换为y=f(t)/(s^2+2*m*n*s+n^2)。在simulink下信号源为f(t),即一个阶越信号(step time舍得很小,如0.01,initial value为t=0时刻的冲击力的值,final value为0),中间为传递函数1/(s^2+2*m*n*s+n^2),右边接一个示波器,示波器的输出就是所求的y了。

tiange 发表于 2006-10-4 10:11

谢谢
页: [1]
查看完整版本: 这个方程怎么编程序(关于悬臂梁受到一冲击)