得到的整体刚度矩阵怎么看

mlfeng

我得到了整体刚度矩阵，可是不知道要怎么看.
帮助里面说的含义如下:
The assembled sparse matrix operator data are written to the text file as a series of comma-separated lists. Each row in the file represents a single matrix entry; a row is written as a comma-separated list with the following elements:

• Row node label
• Degree of freedom for row node
• Column node label
• Degree of freedom for column node
• Matrix entry
  

得到的结果截取一部分如下:
1,1,1,1,1.262499989749938e-001
1,2,1,2,1.262499989749938e-001
1,3,1,3,1.262499989749938e-001
2,1,2,1,2.524999979499878e-001
2,2,2,2,2.524999979499878e-001
2,3,2,3,2.524999979499878e-001
3,1,3,1,2.525000168874888e-001
3,2,3,2,2.525000168874888e-001
3,3,3,3,2.525000168874888e-001
4,1,4,1,2.525000269874892e-001
4,2,4,2,2.525000269874892e-001
4,3,4,3,2.525000269874892e-001
5,1,5,1,2.525000168874888e-001
5,2,5,2,2.525000168874888e-001
5,3,5,3,2.525000168874888e-001
......
请问怎么样看这些数据，要怎么排列得到整体刚度矩阵？
请牛人帮帮忙看一下，麻烦啦~~~

mlfeng

唔，为什么没人回复？
有人帮忙看一下么？
谢谢

zsq-w

这个数据可以用matlab处理下，非常方便。

mlfeng

怎么做？能说具体些么？
我实在搞不懂这些数字的具体含义.....

woser

关注中，我也是
不懂

jacobi

1、这个是单元刚度阵吧？
2、参考ANSYS的话，应该是进行RCM或者AMD重排优化过的节点编号
pasuka 发表于 2009-9-11 20:44

从介绍来看，这个就是整体刚度矩阵，
Matrix generation:

• can output matrices representing the stiffness, mass, damping, and load vectors in a model;
• allows for the mathematical abstraction of model data such as mesh and material information;
• is a linear perturbation procedure;
• includes initial stress and load stiffness effects due to preloads and initial conditions if nonlinear geometric effects are included in the analysis;
• writes matrix data to a text file that can be read as input in other analyses; and
• creates matrices identical to those used in a subspace-based steady-state dynamic procedure (see “Subspace-based steady-state dynamic analysis,” Section 6.3.9).
  

其格式需要琢磨一下才可：

• Row node label
• Degree of freedom for row node
• Column node label
• Degree of freedom for column node
• Matrix entry
  

1 节点

lylguagn

文件中储存的为带状稀疏矩阵（总刚）的部分元素
第1，3为同一个单元中的两个节点号，直接按照mesh得到的号码，没经过优化。文件中每一对号码所包含的若干行代表总刚中的某个子矩阵
第2，4为两个节点交互影响的自由度，相当于子矩阵中的行列号
这就相当于如果将单元刚度矩阵中某两个节点对应的四个子块提取出来重新组合，由于对称性的存在，抛去第一象限部分，那么会留下两个三角和一个正方矩阵，为避免存储重复，排列规律为
按照第3由小到大的顺序，且1的数值大于3的数值
所以在文件中默认并不是以稀疏矩阵的排布规律进行存储，如果要得到带状稀疏的矩阵，还要将相应的元素还原到对应的位置，可以参考cae中的节点排布
如果不在乎计算成本，只想得到某种刚度矩阵，便可以直接以节点号乘以自由度作为其在总刚的初始位置来进行“集成”（此集成只需要找位置）

现在还有个问题是：存储的矩阵元素并不全，初步估计为带宽内计算为0的元素不予记录，由于单元刚阵的计算方程未知，所以还没经过完全证实，大家继续讨论。。。。

lylguagn

8# lylguagn
1、ABAQUS的这种存储，很有可能是分块压缩行存储BCSR或者北大袁明武称之为细胞索引的这类存储格式，可以把半带宽内的零元剔除；

2、半带宽内零元不存储的话，ABAQUS求解的时候势必对总体刚度阵和总 ...
pasuka 发表于 2009-9-13 12:42

aba采用什么技术存储是一个问题，我们得到什么样的刚阵是另外一个问题
现在行列号都有了，只要安放到相应位置，其他位置补0就行了，半带内有否0是结构本身的问题。
至于用什么语言处理，搞清楚了规则，就可以随便编写个脚本就行，不过一个矩阵而已
我想我们的重点还是文件中的这种表达方式的真正含义

ylj1986

我在这里问一下，利用ABAQUS建立的二维复合材料层合板材料，如何才能输出其对应的三维正交各向异性材料的材料属性呢？如果不能的话，可不可以输出其对应的刚度矩阵呢

Linda010

5,3,5,3,2.525000168874888e-001' ~~~~其中5,3,5,3,是什么意思呢

feilideliubo

请问你是怎么得到整体刚度矩阵的啊
我正在求结构的刚度，不知道怎么提取刚度矩阵
请高手指点

dolphin

我也想知道怎么得到，结构的刚度矩阵

shenchenghit

3# zsq-w

请问如何用matlab处理得到刚度矩阵？特别是刚度矩阵比较大的时候。

shenchenghit

12# dolphin
keywords里有："MATRIX GENERATE"，产生整体刚度矩阵

wangjihu220

lz怎么做能得到刚度矩阵啊，能说的详细点吗，

wangjihu220

总体刚度矩阵，需要在INP文件里填哪些内容啊？

wangjihu220

楼主，问题解决了吗？我现在也遇到这个问题，能和我讲解一下啊，谢谢啊！

leemu215

第一行为例
Row node label 1，节点1
Degree of freedom for row node，1节点1的第一个自由度 比如说X方向的位移
Column node label 1，节点1
Degree of freedom for column node 1节点1的第一个自由度 比如说X方向的位移
Matrix entry刚度矩阵的值

这个就应该是K11

tingliu

aba中总体刚度矩阵只贮存非零元素，似乎还是列优先的：即先从第一列开始，然后第二列，...

Maybn

abaqus能得到结构的整体刚度（质量）阵？