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

模具设计培训zb

1、整体刚度和质量矩阵的提取。
6 x3 ~6 h9 B8 Q该功能需要进行二次开发，由ansys形成
的二进制文件.full提取整体刚度和质量
' M& u+ C) t" `矩阵。
( E) o+ [7 G  Z! }( h
基于ansys的一个用户开发程序例子编
# r; z* i* J' K& v  p. k了一个程序（附件中）。

. X, d7 L- M: x7 d; f开发环境：compaq fortran 6.5
运行环境：win2000。
7 s: _: _9 J( k1 ?一个主文件：self.for，
另一个文件matrixout.f90用于矩阵输出
binlib.lib为ansys提供的库文件，将其
引入项目中（也可直接扔进debug目录）
' p! x* e5 z' M0 h( `
' X2 r! y% X' C$ i8 p.full文件由子空间迭代模态分析获得。
5 f" j/ y9 }9 @. n' G运行编译后的可执行文件.exe
. K' Q' R3 ]" R+ }3 T
2、单元刚度和质量矩阵的提取。
9 _2 g( P) e/ K! s8 _- g/DEBUG命令。详细说明可由以下转载文章

finish
/clear

PI=3.1415926

2 B1 V& y. D- n5 `w1=3
* x: ~2 p% G8 @4 Y- D/ {w2=10
- r* _! L/ r% K  t  mw3=6
w4=1.2
r=.8
+ E, Z% Q) j& w" w5 ~t=0.08

/PREP7
9 R" c1 n6 n4 W& J4 ]2 w9 t* S!*
7 i8 u7 r' x/ b' aET,1,SHELL63
R,1,t
" C/ U* ~6 I& M- z5 ^% X
' A3 B! f& z" T) }+ o, qET,2,MASS21
7 u- \' q* `2 ^3 {  S' IR,2,500,500,500,2000,2000,2000,
6 V/ I( \4 S  z5 s# L!*
UIMP,1,EX, , ,2e11
UIMP,1,NUXY, , ,0.3,
4 z: c! @  M  z: H/ f1 P, D5 PUIMP,1,DAMP, , ,0.2,
( ~1 _/ ?+ k6 @) j6 ^9 m  x. tUIMP,1,DENS, , ,7800,
) E8 T( C2 Q7 h9 b, z& z7 X* i
BLC4,0,0,w2,w1
ESIZE,1.5,0,
AMESH,all

NSEL,S,LOC,X,0.0
$ C" n2 j1 E+ o- k/ b4 r' o) @D,all, , , , , ,ALL, , , , ,
3 M4 w  G$ \; `6 G' L* `
allsel,all
+ g" K; s5 {9 X3 dSFA,all,1,PRES,12
FINISH
" M' e9 O8 ~6 J& t. ]1 A
0 l- E0 v0 i5 X) O: p' L6 D/OUTPUT,cp,out,, ! 将输出信息送到cp.out文件
. c2 b. H/ Z. {" ~8 @$ l/ k$ E/debug,-1,,,1 ! 指定输出单元矩阵

/SOLU
- {5 |/ Q% C2 F3 t5 c  y: iSOLVE
" v) F0 {$ v/ T  Xfinish
/ q" j  C  O: \- d5 i3 c' y
/OUTPUT, TERM ! 将输出信息送到output windows中

. r; ^- j8 O& r  i! 这时用编辑器打开cp.out文件，可以看到按单元写出的质量、刚度等矩阵