HELP:压电UEL单元的模态分析

你好,专家,

我正在尝试对用户自定义压电壳单元进行模态分析。我的自由振动的本征值问题在附件中给出。

我的问题是质量矩阵[M]不是正定的，因为电自由度没有与之相关的质量，并且给定的刚度矩阵[K]由于负介电项而是非正定的。我知道，有可能浓缩出上述方程并找到真实的特征值和特征向量，但如何在UEL子程序中使用Lanczos解算器做到这一点。在质量矩阵中引入任意/人工质量项，得到特征值和特征向量。还有其他合适的方法吗?/想法吗?

我需要尽快解决这个问题。我将非常感谢任何关于这方面的帮助。

提前感谢你，

Somu