悬臂梁有限元方法采用龙哥库塔RK方法求解响应，解发散？
用四阶龙哥库塔方法求解一个悬臂梁的振动响应，搞不明白，为什么求解出的结果发散呢？有质量、刚度和阻尼矩阵，再加上一个脉冲响应力，按理说结构的自由度也不多，才20多个，而且求解时的步长也很小，我试着都搞到2Mhz的采样频率了，怎么还是发散解呢？有没有哪位高人遇到过此类的问题，求解释！
本帖最后由 欧阳中华 于 2013-10-6 19:32 编辑
.
这个问题可以这样测试：
1.将计算的时域值做FFT变换，看看得到的频域里峰值是不是对应问题的固有频率，不是就算法有问题；
2.时间取的长一些，建议是最低固有周期的10倍以上，再看看是否发散；
3.将阻尼系数设定为0，看看计算结果。
这些方法都有问题，就找个已知方程作为考题，与解析解比比. . .
欧阳中华 发表于 2013-10-6 19:30 static/image/common/back.gif
.
这个问题可以这样测试：
1 计算出的时域值做FFT的结果就是，一条零线加上最后的地方是个单点峰值
2 时间取的已经足够长了，5秒钟的计算结果还是这样，到后边就是完全是Nan了
3 这个我还没有试过，可以试试
我想问问的是，这个算法有没有人也出现过这种问题，是我自己编程的问题，还是说这是算法本身有问题，然后怎么设置参数才能规避这些问题？
1.先找一个简单的问题调试龙格库塔程序
2.对比其他软件查看质量矩阵，刚度矩阵以及阻尼矩阵
3.你的结果好奇怪，没有频率主峰啊。
页:
[1]