求问导出刚度矩阵

liulenju

我用pde做了一个流体中声波的传播（频域），中心有个point source，study step选的是stationary
现在想把系统矩阵提出来，总是出线问题。
4.3a的版本介绍里说可以提出Modal Solver的系统矩阵，这个Modal Solver是那几种？

mxio

执行：

1. model = mphload('point_source');
   
2. str = mphmatrix(model,'sol1','out',{'K','L','M','N'},'initmethod','init');
   

复制代码

在 str 这个结构体中，就有你想要的刚度矩阵
调用解释：
mphmatrix(模型对象名,解标签,'out',{刚度矩阵们},线性化方法，解序号)
线性化方法 即 initmethod， 有两种选择 init 和 sol
如果选择 init模式 那么 需要联用 Initsol 它有 string 和 zero两种选择
至于输出 'out'
K    刚度矩阵
L    力矢
M    约束矢量
N    约束Jacobian
D    阻尼系数矩阵
E    质量矩阵
此外还有些其他的，不多讲了，具体情况具体分析。。。。

SMT_Captain

学习了，楼上的回答很有启发。

xuxumi

赞！！！！非常好~~