ANSYS中整体、单元刚度和质量矩阵的提取

模具设计培训zb

1、整体刚度和质量矩阵的提取。
该功能需要进行二次开发，由ansys形成
8 _$ q" G: @+ y% |( n1 K3 h( r的二进制文件.full提取整体刚度和质量
矩阵。

基于ansys的一个用户开发程序例子编
/ X7 S! O% j" T2 e( z2 x了一个程序（附件中）。

3 e9 O7 N5 A4 B开发环境：compaq fortran 6.5
7 z2 e( r. e7 X! f! H( g; _6 u& {运行环境：win2000。
- p, e9 ?% I6 h9 d; q一个主文件：self.for，
另一个文件matrixout.f90用于矩阵输出
binlib.lib为ansys提供的库文件，将其
引入项目中（也可直接扔进debug目录）

.full文件由子空间迭代模态分析获得。
运行编译后的可执行文件.exe
" R/ a6 @( U# a6 c% P1 ^
2、单元刚度和质量矩阵的提取。
/DEBUG命令。详细说明可由以下转载文章

finish
/clear
9 Y6 P) O6 G" M
PI=3.1415926
" M, L5 a7 e: d
w1=3
w2=10
w3=6
$ F  F$ P  x' c9 T- u! L- X2 \w4=1.2
r=.8
t=0.08
* X9 J% Y* Z4 j5 e# b' |' B
/PREP7
!*
# O9 _- W+ k9 Q8 a% _" MET,1,SHELL63
R,1,t

ET,2,MASS21
R,2,500,500,500,2000,2000,2000,
!*
9 [; [7 X% S9 [6 g9 S; QUIMP,1,EX, , ,2e11
; x! R4 R- X5 S5 _+ J1 K* X0 QUIMP,1,NUXY, , ,0.3,
* ]% n; S3 K& ]7 }' i/ XUIMP,1,DAMP, , ,0.2,
& `) v& [3 }% v. gUIMP,1,DENS, , ,7800,

' `# @+ ~, ?+ e- F( N0 e6 q: M9 ]BLC4,0,0,w2,w1
0 N5 u- K, l- e3 G8 I% CESIZE,1.5,0,
AMESH,all
$ r' `8 r) ?' \7 M8 Y# ^9 d3 C
NSEL,S,LOC,X,0.0
D,all, , , , , ,ALL, , , , ,
/ e: C& ^! g7 ]% o; t' C) T! V
allsel,all
% j& {3 S& b2 l  u3 B- O% cSFA,all,1,PRES,12
FINISH

  O  y8 S9 W" j  {0 E) {- T: g4 w/OUTPUT,cp,out,, ! 将输出信息送到cp.out文件
/debug,-1,,,1 ! 指定输出单元矩阵
4 N1 P6 y% h/ x0 T; U
) z% a5 y5 U9 A) X  b( B9 T/SOLU
% l) M6 y, \- ^8 ISOLVE
finish

/OUTPUT, TERM ! 将输出信息送到output windows中
7 h( D1 s$ s( L( E2 Q
) O5 p0 G' @* E7 Y6 W+ z5 t! 这时用编辑器打开cp.out文件，可以看到按单元写出的质量、刚度等矩阵