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

模具设计培训zb

1、整体刚度和质量矩阵的提取。
该功能需要进行二次开发，由ansys形成
# h' \1 y4 U4 R的二进制文件.full提取整体刚度和质量
  u6 Y: D) e3 T: q6 t' `. q$ K矩阵。

基于ansys的一个用户开发程序例子编
- S) B, [6 l, t. l: e了一个程序（附件中）。
& H9 }# Q/ G2 F; d' C! K0 E) C
/ z  U8 c, p# J" q- n开发环境：compaq fortran 6.5
: c' o. {7 u4 j$ f运行环境：win2000。
* m( Y, u* Z$ i一个主文件：self.for，
; a! E( k5 p2 f7 Y5 U9 {9 Z另一个文件matrixout.f90用于矩阵输出
- _( R0 x! A) e+ _! }+ p: fbinlib.lib为ansys提供的库文件，将其
引入项目中（也可直接扔进debug目录）

.full文件由子空间迭代模态分析获得。
运行编译后的可执行文件.exe
. M4 i* q3 _# l) O3 j# L* ?; x
1 q3 t) G& Z9 W0 f% _2、单元刚度和质量矩阵的提取。
$ i" y. X; v& k) S5 H" j6 m" J6 M/DEBUG命令。详细说明可由以下转载文章
/ r; g) W. @! {9 n4 e) |
finish
/clear

PI=3.1415926
: X. p' g0 ~" C* b# ?
# {7 K; d. ?) l6 [" z0 Ow1=3
w2=10
w3=6
8 X/ Q4 V2 _! S6 c2 Hw4=1.2
7 x0 ~# S( [) k' z7 Z! Yr=.8
t=0.08

/PREP7
: F6 r5 V9 h) L% X  W# ^6 q1 Y4 ]/ }0 C!*
, A$ D+ D$ U1 n% D. b; kET,1,SHELL63
- k2 V* k( }. m& f: AR,1,t

' t6 i# N  h7 D: QET,2,MASS21
& A" V9 |# G2 rR,2,500,500,500,2000,2000,2000,
3 g, s& F: `( M8 J, p3 M!*
' v% T6 |! r* r6 Q8 UUIMP,1,EX, , ,2e11
UIMP,1,NUXY, , ,0.3,
UIMP,1,DAMP, , ,0.2,
UIMP,1,DENS, , ,7800,
$ R( F' S/ y, V7 n! ]2 D
+ l7 y3 P6 E; G! n2 }BLC4,0,0,w2,w1
! k: X% o5 [; K( y% LESIZE,1.5,0,
- k2 Y6 z( _1 G% \' i. fAMESH,all
) K8 g5 B; N$ ]; h- j2 U
( c- f: ^) u% ^4 Z5 r/ uNSEL,S,LOC,X,0.0
3 R  v9 F) G7 I2 G  u4 s8 w6 oD,all, , , , , ,ALL, , , , ,

1 r+ g# I5 X3 Qallsel,all
' T- M/ d+ @8 `SFA,all,1,PRES,12
FINISH
2 S1 u* t* ]' Y6 S& q1 n
4 p9 X. k/ l  G& G/OUTPUT,cp,out,, ! 将输出信息送到cp.out文件
. e( N; ^4 Z' a& v$ h( }1 p/debug,-1,,,1 ! 指定输出单元矩阵
" O  @; G" _! |8 n: r2 m9 `
/SOLU
SOLVE
8 y: S1 w! I9 S- h# Ffinish
! v8 |, _( }5 S
/OUTPUT, TERM ! 将输出信息送到output windows中
. t/ {: R' l& x
/ o; [1 k4 N9 |5 t: t3 J! 这时用编辑器打开cp.out文件，可以看到按单元写出的质量、刚度等矩阵