ANSYS中整体单元刚度和质量矩阵的提取方法详解

1、整体刚度和质量矩阵的提取。
该功能需要进行二次开发，由ansys形成
的二进制文件.full提取整体刚度和质量
矩阵。
基于ansys的一个用户开发程序例子编
了一个程序（附件中）。
开发环境：compaq fortran 6.5
运行环境：win2000。
一个主文件：self.for，
另一个文件matrixout.f90用于矩阵输出
binlib.lib为ansys提供的库文件，将其
引入项目中（也可直接扔进debug目录）
.full文件由子空间迭代模态分析获得。
运行编译后的可执行文件.exe
2、单元刚度和质量矩阵的提取。
/DEBUG命令。详细说明可由以下转载文章
获得：
发信人: tomath1 (雨天), 信区: FEA
标 题: Re: ANSYS中命令/DEBUG那儿有
发信站: BBS 水木清华站 (Thu Jun 27 14:28:51 2002)
finish
/clear
PI=3.1415926
w1=3
w2=10
w3=6
w4=1.2
r=.8
t=0.08
/PREP7
!*
ET,1,SHELL63
R,1,t
ET,2,MASS21
R,2,500,500,500,2000,2000,2000,
!*
UIMP,1,EX, , ,2e11
UIMP,1,NUXY, , ,0.3,
UIMP,1,DAMP, , ,0.2,
UIMP,1,DENS, , ,7800,
BLC4,0,0,w2,w1
ESIZE,1.5,0,
AMESH,all
NSEL,S,LOC,X,0.0
D,all, , , , , ,ALL, , , , ,
allsel,all
SFA,all,1,PRES,12
FINISH
/OUTPUT,cp,out,, ! 将输出信息送到cp.out文件
/debug,-1,,,1 ! 指定输出单元矩阵
/SOLU
SOLVE
finish
/OUTPUT, TERM ! 将输出信息送到output windows中
! 这时用编辑器打开cp.out文件，可以看到按单元写出的质量、刚度等矩
阵