之前使用matlab模拟本征值的链接如下： https://sunwaybits.tech/eigenvalue-analysis-of-a-beam.html 现在使用bout3中的slepc求解器尝试对这一问题重复操作： int rhs(BoutReal t) { mesh->communicate(u, v); // First-ord…
eigenbox代码 原eigenbox的代码如下： int init(bool restarting) { solver->add(f, "f"); solver->add(g, "g"); return 0; } int rhs(BoutReal t) { mesh->communicate(f); ddt(g) = D2DX…
以一个Euler-Bernoulli Beam模型为例，解释本征模分析，包括理论和代码。下面是本征值问题的最初形态与求解方式。一般先求本征值，再求解本征向量…