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

模具设计培训zb

1、整体刚度和质量矩阵的提取。
. ]0 S: e: i& w' s  x% {该功能需要进行二次开发，由ansys形成
  g8 i9 K" s& x9 A; R; _9 d的二进制文件.full提取整体刚度和质量
" X, p8 K+ y- Y4 h) I; W% o9 s矩阵。
& T" S! c8 c0 y' Z
9 W2 d9 J; p: [9 K基于ansys的一个用户开发程序例子编
1 t% v) Y3 C( c8 A2 \了一个程序（附件中）。

4 T# a9 J+ K8 ]4 v0 P0 `4 [开发环境：compaq fortran 6.5
3 r  X+ |/ I6 H6 z. k8 F6 N运行环境：win2000。
6 A' z, h- J% O9 ~6 R/ K一个主文件：self.for，
0 K1 n* G; P' @2 }另一个文件matrixout.f90用于矩阵输出
+ b6 D6 U, b0 [) Q5 u. y6 fbinlib.lib为ansys提供的库文件，将其
引入项目中（也可直接扔进debug目录）

* S$ `: Q, w! \! q. T) v.full文件由子空间迭代模态分析获得。
$ }$ H' ]; r3 u3 F2 \1 |- E/ Y9 l0 ^运行编译后的可执行文件.exe

2、单元刚度和质量矩阵的提取。
/DEBUG命令。详细说明可由以下转载文章

finish
9 V* p( J( M2 a* o6 E* k: \- }/clear

PI=3.1415926

- ]2 u' `. S% p4 [9 s: @! {w1=3
w2=10
3 t% X# Y! I% F! _6 b$ |3 ~2 nw3=6
w4=1.2
+ h4 @' S) X2 br=.8
t=0.08
( E4 l/ W$ Z- `8 Y
/PREP7
!*
5 u, k$ D' ^  N. C( ^' T  I  g  c4 PET,1,SHELL63
, t7 O" P, g' g6 |, t$ @5 D5 i# MR,1,t
* ^. _5 c# E( g2 S
ET,2,MASS21
R,2,500,500,500,2000,2000,2000,
7 Q6 O' x3 P7 l, q/ I  V1 O/ \+ V) t!*
UIMP,1,EX, , ,2e11
. _' e2 _. k  t3 fUIMP,1,NUXY, , ,0.3,
& U( t; R5 L) t: N. _2 R; I1 I1 \UIMP,1,DAMP, , ,0.2,
- V  X1 w' R  C' O& ]: mUIMP,1,DENS, , ,7800,

$ R  F( e+ _! \* hBLC4,0,0,w2,w1
* S( g8 D3 L! d3 T& SESIZE,1.5,0,
AMESH,all
1 o7 p, |# j) Y8 F2 h9 H
0 [9 V6 c9 Z+ h/ R" WNSEL,S,LOC,X,0.0
D,all, , , , , ,ALL, , , , ,

allsel,all
$ z3 k* b+ r2 R8 P$ t) @SFA,all,1,PRES,12
FINISH

/OUTPUT,cp,out,, ! 将输出信息送到cp.out文件
  [- X6 a4 F8 g/debug,-1,,,1 ! 指定输出单元矩阵
& w6 ]2 v" b% U: k8 ~+ [6 d
- z$ b  A, B8 e, W. \- K" e" ?4 g/SOLU
2 z" I" ?9 x& ]SOLVE
6 M" K" ~/ `# afinish

/OUTPUT, TERM ! 将输出信息送到output windows中
1 _7 E1 S* Y! Z( n; s
! 这时用编辑器打开cp.out文件，可以看到按单元写出的质量、刚度等矩阵