质量阵和刚度阵如何提取？

lisen832000

请问beam188梁的质量阵和刚度阵如何提取？谢谢

anybody_szf

1、整体刚度和质量矩阵的提取。  
该功能需要进行二次开发，由ansys形成  
的二进制文件.full提取整体刚度和质量  
矩阵。  
  
基于ansys的一个用户开发程序例子编  
了一个程序（附件中）。  
  
开发环境：compaq fortran 6.5  
运行环境：win2000。  
一个主文件：self.for，  
另一个文件matrixout.f90用于矩阵输出  
binlib.lib为ansys提供的库文件，将其  
引入项目中（也可直接扔进debug目录）  
  
.full文件由子空间迭代模态分析获得。  
运行编译后的可执行文件.exe  
  
2、单元刚度和质量矩阵的提取。  
/DEBUG命令。详细说明可由以下转载文章  
  
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文件，可以看到按单元写出的质量、刚度等矩阵

lisen832000

整体刚度和质量矩阵的提取能不能再说的详细点?