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

模具设计培训zb

1、整体刚度和质量矩阵的提取。
该功能需要进行二次开发，由ansys形成
的二进制文件.full提取整体刚度和质量
矩阵。
5 b- a! T2 W  |% x" Y9 Q2 b
基于ansys的一个用户开发程序例子编
了一个程序（附件中）。

开发环境：compaq fortran 6.5
运行环境：win2000。
一个主文件：self.for，
6 }5 w  P' A" q0 k另一个文件matrixout.f90用于矩阵输出
: g  S0 {. p6 y! W4 nbinlib.lib为ansys提供的库文件，将其
5 A3 k  g3 {  Y引入项目中（也可直接扔进debug目录）

.full文件由子空间迭代模态分析获得。
运行编译后的可执行文件.exe

2、单元刚度和质量矩阵的提取。
/DEBUG命令。详细说明可由以下转载文章
; e! @& V2 z; g, C! K' X
finish
/clear
7 v8 f2 \" k; r$ t, ]+ w
PI=3.1415926

9 Z3 |8 ?5 |# k$ \& G/ Uw1=3
w2=10
w3=6
/ e$ W2 C9 H# b" R# iw4=1.2
r=.8
* Y" ?, l$ ]+ ~6 R/ nt=0.08

/PREP7
!*
2 {( S  U8 i& T8 N, LET,1,SHELL63
0 \- z: I4 S4 o3 C# DR,1,t
1 x! g, I% l# J) {
7 U% G1 ~7 i% B  _) M" o! n6 B& ?ET,2,MASS21
: E3 _/ a3 R, k6 J7 M4 p; mR,2,500,500,500,2000,2000,2000,
!*
UIMP,1,EX, , ,2e11
6 O6 E4 U" Q' ]7 e$ i+ kUIMP,1,NUXY, , ,0.3,
% B1 V9 ^- k" n( l% U- u0 XUIMP,1,DAMP, , ,0.2,
# l$ b- ~! _; nUIMP,1,DENS, , ,7800,
" |- G" N% j4 s$ j
) M/ `# U( }( A/ P1 Y. y3 N, ~BLC4,0,0,w2,w1
ESIZE,1.5,0,
  `3 Q, T. I3 Z9 k& W3 p$ g4 q. yAMESH,all
. E% A) a# k6 i: U/ }; d
$ ~  _) z) D: ]% K5 n) T# BNSEL,S,LOC,X,0.0
D,all, , , , , ,ALL, , , , ,
' t  \! ^. i2 S  Z7 [; q* O! t
/ x; @; {" B! uallsel,all
SFA,all,1,PRES,12
0 K3 t1 N2 S" i8 X5 A( |8 y( X  @FINISH

/ g; P9 r3 O3 {) M: y& `1 S/OUTPUT,cp,out,, ! 将输出信息送到cp.out文件
- X0 j  O& T' A2 Y4 ^) E  e/debug,-1,,,1 ! 指定输出单元矩阵
# D" z, N) E0 n. Q
' E2 Z( T% ^5 m& f, q0 U) {6 ?/SOLU
+ c' R1 {$ ~" `. q3 X( vSOLVE
finish
4 G  q& z$ f9 F: q
4 a& G: a% ~8 E7 \1 ^+ i$ w/OUTPUT, TERM ! 将输出信息送到output windows中

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