- 积分
- 1
- 注册时间
- 2007-1-26
- 仿真币
-
- 最后登录
- 1970-1-1
|
发表于 2008-8-24 17:17:12
|
显示全部楼层
来自 河南郑州
引用一下。呵呵
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,这时的矩阵才有意义。
如果在模型中使用了子结构,对每个子结构都重复上面的步骤,最后是整个结构。 |
|