我有一个具有以下运动方程的大型系统

$M\ddot z(t)+K z(t)=F(t)$

$M$是质量矩阵

$K$是刚度矩阵

我想用负载相关的 Ritz 分析求解位移的方程、特征值和特征向量（该算法在此处简要描述）。$z(t)$

我知道我可以使用APRACK 来解决大规模特征值问题，但是负载相关的 Ritz 分析呢？是否已经有一个求解器（商业的、开源的，一切都无关紧要）可以做这种事情？

为了清楚起见，我正在寻找我可以使用的编程库，而不是我无法以编程方式调用的商业软件。