如何推导粘滞阻尼器单元的刚度阵/阻尼阵？

xiaoqiushui318

如何推导粘滞阻尼器单元的刚度阵/阻尼阵或者说如何在有限元程序里实现粘滞阻尼器单元？
Sap2000里有一种粘滞阻尼器单元，用Maxwell模型描述，单元由一个弹簧串联一个阻尼器得到。其内力由下式描述
        f=K*dk=c*dvc^beta
其中K为弹簧弹性常数，c、beta为阻尼器常数，dk为弹簧的变形，dvc为阻尼器的速度。
Sap2000里关于此单元的说明见附件。
请问这样一个单元如何能到处其单元刚度矩阵或阻尼矩阵呢？