『求助』一个关于动力学方程求解的问题
浏览：532356 回答：1
动力学方程
M*X''+C*X'+K*X=R(X,X')
M为质量矩阵，C为阻尼矩阵，K为刚度矩阵，X为广义坐标列阵，X'为广义速度列阵，X''为广义加速度列阵，R(X,X')为外激励列阵。M和K都是正定矩阵。X、X'、X''都是时间t的函数，但是X、X'、X''与t之间的表达关系式未知。
如果想要求从时间t1变化到时间t2时，X的变化，Matlab的程序该怎么编？假设t1＝0时，X、X'、X''都为0列阵，
原来想用NewMark法求解，但是由于R(X,X')中含有关于时间t的未知函数关系，就不知道该怎么写了。
M*X''+C*X'+K*X=R(X,X')
M为质量矩阵，C为阻尼矩阵，K为刚度矩阵，X为广义坐标列阵，X'为广义速度列阵，X''为广义加速度列阵，R(X,X')为外激励列阵。M和K都是正定矩阵。X、X'、X''都是时间t的函数，但是X、X'、X''与t之间的表达关系式未知。
如果想要求从时间t1变化到时间t2时，X的变化，Matlab的程序该怎么编？假设t1＝0时，X、X'、X''都为0列阵，
原来想用NewMark法求解，但是由于R(X,X')中含有关于时间t的未知函数关系，就不知道该怎么写了。
| R(X,X')是滑动轴承的非线性油膜力。如果用符号函数求解，求解的时间太长了。 请各位兄弟姐妹帮帮忙了。 |
![[免费案例]Ensight案例教程分享](https://img.jishulink.com/cimage/245b3ca9e2c939e40491a25edae94515.jpeg?image_process=resize,fw_576,fh_320,)
楼主会了吗