怎么从ADINA中导出刚度矩阵
怎么从ADINA中导出刚度矩阵,哪位高手给说下啊谢谢 引用一下。呵呵Q:怎样将总体刚度矩阵和质量矩阵保存到文件当中?
A:可以用命令DISK-STORAGE GLOBAL-MATRIX=SAVE 总体刚度矩阵以及隐式动力分析中的质量矩阵。在用户界面中,该选项是Control > Miscellaneous File I/O… ,然后在对话框中选择Save the Global Matrix 。在Save Step 中定义非线性分析中需要保存矩阵的时间步。对于非线性分析,所保存的是开始迭代之前的矩阵。
在Unix 系统中,数据保存在<problem_name>.mtx 文件中。对于Windows 系统,文件名为fort.70 。在下一个版本中将改为<problem_name>.mtx 。
保存数据的Fortran 格式如下所示:
WRITE (70,*) ' MAXA ARRAY'
WRITE (70,1000) ((I,MAXA(I)),I=1,NEQ+1)
1000 ?FORMAT(1X,I8,5X,I8)
WRITE (70,*) ' GLOBAL STIFFNESS MATRIX'
WRITE (70,1001) ((K,A(K)),I=1,NWK)
1001 ?FORMAT(1X,I8,5X,E20.13)
WRITE (70,*) ' GLOBAL MASS MATRIX'
然后保存质量矩阵。对于一致质量矩阵,格式为:
WRITE (70,1002) ((K,AM(K)),I=1,NWK)
1002 ?FORMAT(1X,I8,5X,E20.13)
对于集中质量矩阵格式为:
WRITE (70,1002) ((K,AM(K)),I=1,NWK)
1002 ?FORMAT(1X,I8,5X,E20.13)
参数:
NEQ-方程的总数
NWK-矩阵元素的总数
可以在<problem_name>.out 文件中找到这两个参数。
注意:只有使用DIRECT 求解器(使用高斯消去法),命令行为MASTER SOLVER=DIRECT,这时的矩阵才有意义。
如果在模型中使用了子结构,对每个子结构都重复上面的步骤,最后是整个结构。 楼上的是牛人 fort.70 中
ADINA 的刚度矩阵和质量矩阵是一维变带宽存储的,数据前面有一个辅助数组记录数据的位置而且数组中最后一个数据是数组长度,把它看明白了 就可以编制个小程序把矩阵转换为原来形式 2楼是牛人,不过我还没看懂
页:
[1]