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

模具设计培训zb

1、整体刚度和质量矩阵的提取。
该功能需要进行二次开发，由ansys形成
+ K! O, O7 k( X+ ^/ M的二进制文件.full提取整体刚度和质量
, r$ Y' b  t# z# ~; `$ O! c矩阵。

/ F( B. ~' s3 ^0 Z基于ansys的一个用户开发程序例子编
了一个程序（附件中）。

开发环境：compaq fortran 6.5
. K: v9 q- Q! X; \/ t1 S' b运行环境：win2000。
; t, Q% h/ E' o* q: s一个主文件：self.for，
另一个文件matrixout.f90用于矩阵输出
binlib.lib为ansys提供的库文件，将其
5 M: z8 G$ O) \0 s' k引入项目中（也可直接扔进debug目录）
3 A$ \4 \9 N$ Y' s, M- X! K
- x& a9 h8 `$ |0 {+ z& }.full文件由子空间迭代模态分析获得。
( t7 E& S. o1 ^+ x. D运行编译后的可执行文件.exe

2、单元刚度和质量矩阵的提取。
# d  O0 B: z; M  J( S/ n/DEBUG命令。详细说明可由以下转载文章

finish
2 f" N' o, H) M/clear

' K( w7 r  V( V$ a; nPI=3.1415926
; B( \; F- ?5 l' {. S2 M, A
2 Q5 I  k; E/ K# @; c# z4 ?5 lw1=3
6 B9 |9 O, C8 V3 d8 J2 pw2=10
w3=6
w4=1.2
3 ?! ~. q- _$ tr=.8
t=0.08

' H4 c0 z' C6 A/PREP7
4 b' s% t. ?& _) y) g% N$ g!*
ET,1,SHELL63
5 V/ t5 V: C$ z( L+ I$ uR,1,t

  Q2 d7 E, X  ^. {. L" PET,2,MASS21
, O  Z$ n/ @. X1 a$ x# lR,2,500,500,500,2000,2000,2000,
!*
" O1 \" z& A, T4 t; I8 b" hUIMP,1,EX, , ,2e11
UIMP,1,NUXY, , ,0.3,
/ k+ U/ R& ^* _* SUIMP,1,DAMP, , ,0.2,
; O: d3 O# B, @0 pUIMP,1,DENS, , ,7800,
$ X& m! ?$ {9 Y) u
  e' x6 u1 A7 H  {, }% uBLC4,0,0,w2,w1
# Y  @7 `+ K9 P* `; {ESIZE,1.5,0,
0 C2 }. d& v# G9 W& X+ kAMESH,all

NSEL,S,LOC,X,0.0
; R2 M8 P+ P' a4 y2 hD,all, , , , , ,ALL, , , , ,

, h+ v- n5 W0 w# P- Sallsel,all
. K/ u$ M# `" Z; tSFA,all,1,PRES,12
7 v2 E& g) c" h: kFINISH
9 o- D- _8 o) j/ p8 C7 ^+ \( W
) ^, A/ `/ s& G  b$ b* ]/OUTPUT,cp,out,, ! 将输出信息送到cp.out文件
0 W' m" _  @, O0 R/debug,-1,,,1 ! 指定输出单元矩阵
& P8 C& R9 l$ B6 ?% s# ^. G9 l% _$ {8 ?
/SOLU
: a0 e. N$ r6 jSOLVE
finish
$ Y$ b. p, T: |5 M
( q  W( p  r1 q, Z* `9 }! Z8 D/OUTPUT, TERM ! 将输出信息送到output windows中
1 k. u' U% k8 |7 [. f
5 ^" c! |5 a2 R4 I! 这时用编辑器打开cp.out文件，可以看到按单元写出的质量、刚度等矩阵