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

模具设计培训zb

1、整体刚度和质量矩阵的提取。
+ {" n- R8 E6 W+ C  e该功能需要进行二次开发，由ansys形成
2 ^2 R" a9 e1 M, [的二进制文件.full提取整体刚度和质量
6 t  }2 @  {; r9 w6 K3 {  _( M矩阵。
5 B& a' M) I  D1 H4 u! L9 x% q. _
, h. ?$ N4 |% r1 M' c3 D基于ansys的一个用户开发程序例子编
: Z. t: m, I- t" W( s了一个程序（附件中）。

开发环境：compaq fortran 6.5
运行环境：win2000。
一个主文件：self.for，
另一个文件matrixout.f90用于矩阵输出
binlib.lib为ansys提供的库文件，将其
2 x' |. P, m5 r引入项目中（也可直接扔进debug目录）

.full文件由子空间迭代模态分析获得。
3 ?% o3 L( N4 l' a, l. F* H9 |运行编译后的可执行文件.exe
4 K- P- m! u5 k. [
2、单元刚度和质量矩阵的提取。
- q) M3 W% q( F8 y9 Q( J' l' p/DEBUG命令。详细说明可由以下转载文章

; X- r2 W6 J# f  Z  n. Ifinish
/clear

- c# b0 D/ S+ A4 [PI=3.1415926

w1=3
: S8 \( l, n2 @' q: U% N8 D' i4 lw2=10
; f4 u0 X# c7 o/ Bw3=6
0 S! I7 l  u, e8 a; [! E/ C  Q, lw4=1.2
r=.8
t=0.08

3 x6 I: t6 {( p$ v& n6 Q; X5 b/PREP7
' [+ G; B7 e+ v2 J!*
ET,1,SHELL63
! j0 x. l$ t' p0 mR,1,t

4 q$ G0 f6 \  {4 m' o  WET,2,MASS21
R,2,500,500,500,2000,2000,2000,
!*
UIMP,1,EX, , ,2e11
& b" {/ l# J. g, G3 X" IUIMP,1,NUXY, , ,0.3,
UIMP,1,DAMP, , ,0.2,
UIMP,1,DENS, , ,7800,

  Y! [8 q! ~0 z: y, vBLC4,0,0,w2,w1
ESIZE,1.5,0,
" S* j6 [" Q8 S; _" AAMESH,all

NSEL,S,LOC,X,0.0
D,all, , , , , ,ALL, , , , ,

allsel,all
4 [" y' L( u. P6 H5 I$ JSFA,all,1,PRES,12
- ~5 V9 K* F. T; l! y8 j1 wFINISH
! m! o  X; L" B2 @% X
/OUTPUT,cp,out,, ! 将输出信息送到cp.out文件
. u! _3 f! `% j2 I7 {6 u/debug,-1,,,1 ! 指定输出单元矩阵

/SOLU
0 E! Z% ^5 X* [SOLVE
finish

/OUTPUT, TERM ! 将输出信息送到output windows中
+ K" B- S# r6 y
7 P! k7 m, V  {% Z5 Y! 这时用编辑器打开cp.out文件，可以看到按单元写出的质量、刚度等矩阵