基于组件的三维CAD系统开发的关键技术研究（四）

HEATS

　　2 .7应用软件前、后处理

Z# r/ | N2 U2 F0 j( r

　　由于软件升级滞后，一些常用软件后处理功能很弱或没有，可以利用LSP程序结合AutoCAD增加或简化前后处理功能。下面结合平面渗流分析程序STSE软件对其前、后处理功能的实现进行介绍：

& }; k+ `% Z I3 i" a" F

　　(1)前处理

) ]% F) [& L, [- s& [

　　STSE为平面渗流有限元程序，单元划分和单元、节点编号工作量最大，可以借助其它通用有限元软件的前处理功能如Ansis、Algor、SAP84等进行初步处理，生成单元、节点编号和节点坐标，而后按照STSE数据文件的格式要求对数据文件进行编辑。

　　(2)后处理

6 h5 r1 @0 T# p. }, @

　　后处理成果主要为：单元网格图，浸润线和等势线。单元网格图中有节点和单元编号，不同渗透系数的单元采用不同的颜色，以便根据网格图直观地判断数据文件中几何参数和材料特性的正误；

# o: W' N" Q! H& m. b ! j2 R9 I) Z( I- s9 u, e

　　首先分析STSE的输出结果文件的格式，搜索并筛选其中主要参数如单元总数、节点总数和材料总数，把节点坐标和单元信息分别存储，通过对单元循环生成单元网格图。

* i3 [- q; o* L- \ V 4 y2 O7 n2 a/ b9 q4 H

　　(defun c:seepmesh()

5 U3 G* B- o/ t5 p * R, E/ H- i5 X+ A5 ?* T( z

　　(if fn

2 ^/ R6 _9 a* d4 _0 a6 ~1 p

　　(setq fn (getfiled "渗流结果文件名" fn "" 2))

0 g8 B' b7 n. T

　　(setq fn (getfiled "渗流结果文件名" "" "" 2)))

7 S4 V0 R: f" I* x# L

　　(setq f (open fn "r"))

* j" O; z6 U `# Y6 _

　　(setq md (read-line f))

" _( G6 V) K( v1 _" b4 h/ B- e

　　(while (/= (substr md 25 8) "单元总数")

　　(setq md (read-line f))

- G+ {1 D( B8 K, T1 i8 P% A

　　)；定位单元总数

. @ p: G: K2 V. [1 t) a6 g# y , \; m4 {* I: x* x5 @3 g

　　(setq dyzs (atoi (substr md 58 5)))；读单元总数

; ~! n& o2 _: p' s" c& x

　　(setq md (read-line f))

0 Z4 _4 }6 w5 k6 J! K & [4 _& p" y0 `" G- _* W, l

　　(setq jdzs (atoi (substr md 58 5)))；读节点总数

" @9 ~3 |. G! k" M+ m

　　(setq md (read-line f))

　　(setq clh (atoi (substr md 58 5)))；读材料总数

: H( h |! {1 ~9 |

　　(setq m 0)

3 S: @+ H$ O' }

　　(repeat clh

0 z1 T& \" m3 _ ) ]( x: @ |+ `9 Y! A1 F4 z% ?

　　(setq m (+ 1 m))

2 e; j7 ?1 l! }$ ^) Z* {

　　(setq tcm (strcat "zclh" (itoa m)))

3 e4 A2 W7 B0 ~! v8 o2 Y

　　(command "layer" "m" tcm "c" (itoa m) tcm "")

. b6 C# i$ n; g3 E9 a. H ! L# x( r; g9 o) R! r2 _$ W* e

　　)；按材料种类生成图层名称

( F, ^9 C0 q& N4 O# T

　　(setq m 0 n 0)

e* v) Y/ p8 O4 Q - g/ H% T/ u# D0 H! {

　　(repeat dyzs；对单元循环，dycfb中存放单元信息

5 R5 ?* W/ _- x; ~7 x' L/ S6 W# y * V) J* l7 q2 A2 D# \

　　(setq m (+ m 1))

. p2 e) [8 `. |( o* F- p

　　(setq clh (nth 0 (nth m dycfb)))

' `4 H5 F) u: U5 {* |' i P( z1 i' F8 p2 }2 s9 l% P

　　(setq jdh1 (nth 1 (nth m dycfb)))

7 S$ j5 W: o) A7 U/ o) r

　　(setq jdh2 (nth 2 (nth m dycfb)))

/ n) q- M0 c, v7 _

　　(setq jdh3 (nth 3 (nth m dycfb)))

0 y9 ]: b, [4 L( o2 a / F" u. e7 U; U8 c

　　(setq jdh4 (nth 4 (nth m dycfb)))

, B+ C; S7 }5 K

　　(command "layer" "s" (strcat "zclh" (itoa clh)) "")

6 U$ _: O% i9 o- ]0 A) J. ^ . [( H! `- Q; y2 W

　　(command "pline"；绘制单元网格

6 E u( E1 \/ m+ Y0 } ) F; J6 g6 N4 E% [4 R

　　(nth jdh1 jdzbb)

　　(nth jdh2 jdzbb)

　　(nth jdh3 jdzbb)

4 G; c& o+ H& e* H

　　(nth jdh4 jdzbb) "c")

% q( J7 \( _1 n1 k$ }- o. x 7 _4 z) J8 o2 M- s$ Q( e) ^ d

　　(setq bzdzb (mapcar '+ (nth jdh1 jdzbb) (nth jdh2 jdzbb)

P- h# N; f1 u& l

　　(nth jdh3 jdzbb) (nth jdh4 jdzbb)))

9 F* {5 p4 s5 o3 `# w! v) x$ j . _9 e% A l( N H

　　(setq bzdzb (list (/ (car bzdzb) 4.0) (/ (cadr bzdzb) 4.0)))

4 K4 n% \6 Y% Y+ W- `$ x- ` 0 D. K' Z) t8 t- v+ L7 X) m; c1 g9 P

　　(command "text" "m" bzdzb "0.5" "0" (itoa m))

P' o" @# ~* h

　　；标注单元号

; z: _5 e# u9 A! ` {5 \ 0 _9 ? c. j+ W7 P* z" G6 }

　　(setq m 0)

( v9 n( s" @! {9 F' k, W1 B) \

　　(repeat jdzs

% b! c+ b0 h* e6 c0 ~. x. R+ n9 ? " R' P' s" m; O9 j

　　(setq m (+ m 1))

" `; v! T4 ~( p. c' f8 ?

　　(setq bzdzb (nth m jdzbb))

: X% b$ L& }% ^* \

　　(setq bzdzb (list (+ (car bzdzb) 0.0) (+ (cadr bzdzb) 0.0)))

6 s7 R$ y) d( c9 d" E5 S0 _. Z! p

　　(command "text" "m" bzdzb "0.5" "0" (itoa m))

　　))；标注节点号

9 B# R- r1 `- `3 W* N4 \ g( J

　　(close f)

　　)

6 k! _4 {, u5 L, c) K I ( E( O0 x: I- k: u! e* U- A4 {2 B

　　限于篇幅，等势线和浸润线生成程序不再赘述。同样，对工程设计中遇到的其它软件如：STAB边坡稳定分析、SAP84的输出文件均可进行类似后处理工作。

　　3结语

; ^1 S ~" M/ |) s0 f; I( C: B

　　AutoCAD软件应用已经在设计单位普及，VisualLisp的出现对于提高Lisp语言的编程效率和AutoCAD应用水平具有重要意义。本文涉及的几个实例已经在实际工作中发挥很大作用，充分说明在工程设计当中结合本专业要求开发一些简单实用的Lisp小程序可以有效减轻设计人员的劳动强度，使AutoCAD不仅仅是绘图工具，而真正成为设计人员得心应手的高效设计平台。

4 L& z2 I; Y* b! t* n( K1 e 6 X& F t& ?- S: U# ~" Q6 ~6 \