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

模具设计培训zb

1、整体刚度和质量矩阵的提取。
& S& ^) o3 u3 |% K3 D. I) Q该功能需要进行二次开发，由ansys形成
" A& s# |+ m" }4 Q; H的二进制文件.full提取整体刚度和质量
矩阵。

基于ansys的一个用户开发程序例子编
. L5 P: X- N: k. J" V: E+ ?; _了一个程序（附件中）。

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

: [! w9 Q, V% P  l  Y.full文件由子空间迭代模态分析获得。
运行编译后的可执行文件.exe
* }# H2 n+ E* A0 Z, w! ]( K# u' _
2、单元刚度和质量矩阵的提取。
( ^* W, \% S# B6 {. R  Y# H/DEBUG命令。详细说明可由以下转载文章

+ J4 W5 E, ]5 ?3 S- k8 j; C1 C. yfinish
/clear
' \+ u% |$ D7 J! e  l; o
PI=3.1415926
2 R" ^9 [: }- X' Y# T
w1=3
& }3 |# @# R8 v$ c7 d, Uw2=10
, W4 T0 r8 \# v- q6 v1 C3 Jw3=6
* O' f6 g/ U- K0 h4 Cw4=1.2
r=.8
t=0.08
+ [# z+ r' v# Y1 e5 J1 h$ q6 H9 U
  C# [; x' s; Z. _* R* s/PREP7
!*
0 ]# a% d  g; IET,1,SHELL63
R,1,t

7 k$ G* P7 K" E% f- RET,2,MASS21
R,2,500,500,500,2000,2000,2000,
!*
* _$ N5 p: d( K8 ?7 A& {3 j1 F, XUIMP,1,EX, , ,2e11
* @( j$ z* p, ~9 }UIMP,1,NUXY, , ,0.3,
UIMP,1,DAMP, , ,0.2,
UIMP,1,DENS, , ,7800,
3 w; ?& @6 j6 E& g8 h- X
2 f: s1 P# l( V' U3 \BLC4,0,0,w2,w1
ESIZE,1.5,0,
6 v/ g; G# }' ]2 MAMESH,all
2 {2 y0 V: x" s2 ?& {) ?' E
' v0 M( N) Y/ x; JNSEL,S,LOC,X,0.0
& J# ^1 s' v  M  Y* a$ aD,all, , , , , ,ALL, , , , ,
: b; L) x7 \- C- e( |' `# w
allsel,all
SFA,all,1,PRES,12
FINISH
  h0 b& {9 Z7 y) p8 p
/OUTPUT,cp,out,, ! 将输出信息送到cp.out文件
0 m, b6 Z7 k. N, }- i/debug,-1,,,1 ! 指定输出单元矩阵
* {( H; [! ^* |9 i% c5 a
8 y. x) R$ K3 I# h/SOLU
7 k' A$ g0 K& ?SOLVE
# k8 E6 K# X) {8 H3 Q. c& A8 T/ Yfinish
! [6 j" e4 I2 T) r1 B- ^2 }" e) b5 B
/OUTPUT, TERM ! 将输出信息送到output windows中
; s+ |3 G: F; p  L$ r
! 这时用编辑器打开cp.out文件，可以看到按单元写出的质量、刚度等矩阵