Abaqus-提取矩阵

        Abaqus中涉及的矩阵注意包括：质量矩阵，刚度矩阵，阻尼矩阵，载荷矩阵等。有时候，我们想提取模型的刚度矩阵，该如何实现呢？目前提取矩阵只能通过编辑关键字实现。

一、提取单元矩阵

在Step定义关键字后面添加如下关键字：

** STEP: Step-1***Step, name=Step-1, nlgeom=NO, perturbation*Frequency, eigensolver=Lanczos, acoustic coupling=on, normalization=displacement, , 10., , ,*Element Matrix Output,Elset=SET-2,File Name=filename,Frequency=1,Output File=User Defined,Stiffness=Yes,Mass=Yes

** STEP: Step-1
**
*Step, name=Step-1, nlgeom=NO, perturbation
*Frequency, eigensolver=Lanczos, acoustic coupling=on, normalization=displacement
, , 10., , ,
*Element Matrix Output,Elset=SET-2,
File Name=filename,Frequency=1,Output File=User Defined,Stiffness=Yes,Mass=Yes

** STEP: Step-1
**
*Step, name=Step-1, nlgeom=NO, perturbation
*Frequency, eigensolver=Lanczos, acoustic coupling=on, normalization=displacement
, , 10., , ,
*Element Matrix Output,Elset=SET-2,
File Name=filename,Frequency=1,Output File=User Defined,Stiffness=Yes,Mass=Yes

Elset, 是指需要提取矩阵的单元集 合. File Name, 为生成的矩阵文件自定义文件名. 在非线性分析中，可以使用Dload=Yes 来提取单元集 合的载荷分布矩阵。

二、提取全局矩阵

提取全局矩阵也是相当容易的，只需在inp最后添加如下关键字：

***STEP*MATRIX GENERATE, STIFFNESS, MASS*MATRIX OUTPUT, STIFFNESS, MASS, FORMAT=MATRIX INPUT*END STEP

计算完成后，在文件夹下会生成名为jobname_STIF2和jobname_MASS2的.mtx格式文件. 打开jobname_STIF2文件你会发现如下文字：

1,1, 1,1, 1.000000000000000e+361,5, 1,1, 9.283850055777577e-121,6, 1,1, -4.181074365176346e+04

其中，第一列为节点所在行编号，第二列为行节点自由度编号，第三列为节点所在列编号，第四列为列节点自由度编号，第五列为矩阵值。

*需注意的是，提取质量矩阵时，必须要定义密度，因为刚度矩阵Kx=F与质量无关。所有结构分析都涉及刚度矩阵的计算：

静力学分析: Kx = f
模态分析: Ma +Kx = 0
动力学分析: Ma + Cv + Kx = f(t)

*注意2：这些关键字不支持Explicit求解器，若添加，求解会报错。

而在Dynamic implicit则可以使用如下关键字：

** *STEP *MATRIX GENERATE, STIFFNESS, MASS, STRUCTURAL DAMPING, VISCOUS DAMPING*MATRIX OUTPUT, STIFFNESS, MASS, STRUCTURAL DAMPING, VISCOUS DAMPING, FORMAT=COORDINATE*END STEP