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

模具设计培训zb

1、整体刚度和质量矩阵的提取。
! U6 w3 H0 x( s2 {! M9 ]9 i该功能需要进行二次开发，由ansys形成
的二进制文件.full提取整体刚度和质量
矩阵。
; T, D8 H0 b$ y$ y- A9 Q
" O" W( N1 O" F0 m6 |$ [基于ansys的一个用户开发程序例子编
了一个程序（附件中）。

1 F, Y) _5 ?# E' k3 D开发环境：compaq fortran 6.5
  T- w& J- [: K4 g/ d( D% R( j运行环境：win2000。
- `" J1 a2 b; i. R$ U一个主文件：self.for，
另一个文件matrixout.f90用于矩阵输出
binlib.lib为ansys提供的库文件，将其
引入项目中（也可直接扔进debug目录）

; S& F1 I3 [+ k4 K  J# v.full文件由子空间迭代模态分析获得。
/ Y3 ?# |8 H- C( c0 w; H7 |* c2 [运行编译后的可执行文件.exe

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

finish
/clear
2 r1 d' D) ?# f
1 G9 B8 z% |) x1 c, F/ IPI=3.1415926

w1=3
w2=10
w3=6
w4=1.2
( s$ y; e7 r/ z- }& w$ e' Mr=.8
! g; r' n% k1 ?; z6 J7 L; yt=0.08

5 F- M# e4 V& r. u( ?1 ~! A, x4 l/PREP7
!*
ET,1,SHELL63
. Y; v2 f. o3 w; R6 {; jR,1,t
: B( U- I" B' i' V; r& Z
7 D1 X& A" k3 M( G" i" N6 {ET,2,MASS21
4 I6 [2 }6 n, X; [5 H% C8 vR,2,500,500,500,2000,2000,2000,
!*
+ [6 T2 ^; R7 c4 m/ V" @% KUIMP,1,EX, , ,2e11
UIMP,1,NUXY, , ,0.3,
UIMP,1,DAMP, , ,0.2,
* G' Y; K: ?# s1 n6 N; I. G  rUIMP,1,DENS, , ,7800,

BLC4,0,0,w2,w1
- n$ P$ M' Y/ t2 K4 J' {3 _ESIZE,1.5,0,
+ c) ?1 K9 b' U9 ?# BAMESH,all

NSEL,S,LOC,X,0.0
D,all, , , , , ,ALL, , , , ,
9 O& Z* o% m" E- l
allsel,all
- D. c, l  L9 [# N# Z& kSFA,all,1,PRES,12
FINISH

% m% t  @9 T) ^5 ?- _4 ^# R/OUTPUT,cp,out,, ! 将输出信息送到cp.out文件
/debug,-1,,,1 ! 指定输出单元矩阵
8 Z! ?$ q0 x( [8 N, l% g
# F% U  d+ X9 x" o  u; o/SOLU
9 A9 N% r5 E' J" cSOLVE
! z& h. O+ ]; i6 C  Zfinish

/OUTPUT, TERM ! 将输出信息送到output windows中

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