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

HEATS

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

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

4 U( C4 e+ a) I3 W5 Z " R* Z$ }* B- b' U: l; z; j) R

　　(1)前处理

8 E3 E$ w9 c/ B3 l' x& N6 X# q$ B

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

　　(2)后处理

/ Q; g/ }" N) l1 _: c% ~. O. u8 D

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

) D8 h9 P' C$ Y) `) c, k5 q! B4 r1 p) B

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

S' d* y8 B, ^/ F; w2 B

　　(defun c:seepmesh()

　　(if fn

7 r4 ^ o: d8 D8 h* a! h

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

% b- n- `7 B2 {6 c7 C

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

" y2 T9 b" F% M e9 T4 G5 G0 |+ k

　　(setq f (open fn "r"))

# n9 [! [% n; e/ ~2 S5 J7 ^

　　(setq md (read-line f))

( P3 P- U/ S. b) R 9 u1 `6 T5 v3 u4 m

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

　　(setq md (read-line f))

) d2 P4 a: G& Z6 N

　　)；定位单元总数

0 \4 d% r: L5 ^4 K

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

& R% \, R8 G+ d( A3 B$ i" v" g, M: F

　　(setq md (read-line f))

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

% M% c. [9 n- R- I

　　(setq md (read-line f))

9 f0 D" Q7 c4 V; } 9 ^, ^ t$ A3 {0 U, O+ {: I1 T

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

/ G/ H) \# u' A( e7 r0 S 3 ?) U/ P/ \/ s* @

　　(setq m 0)

, l6 S* j" k, ~- ]5 Z0 K% S, U

　　(repeat clh

3 F7 Q/ o! E/ F2 q G 2 U2 ^/ H* I' E$ n9 d4 G$ C

　　(setq m (+ 1 m))

. e& Q/ l* a: N$ @0 _: h

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

1 I6 {. s# V; ^4 X, w7 q ' C G- Z1 P" T9 L+ `' A

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

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

u7 x/ l7 R" C/ x; P

　　(setq m 0 n 0)

5 @+ S1 U1 k& v) g

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

7 H1 x# J" A9 R' Q) V( L

　　(setq m (+ m 1))

R0 O$ }' J' ^. ^# E7 u- D

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

/ M. J2 C( L, b: U ' J$ q: S- Z$ a- H* U

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

1 I5 J. }$ U3 |# v3 e% S

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

+ u2 a4 ~& Q2 U7 a1 x

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

0 C" ]: E' G- k ' V9 n4 W- Y1 \7 m; ~& T

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

2 r& B* ? x" Z8 C. H' V( o

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

3 g$ ~/ z2 q2 G/ ?9 n: _1 L; p

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

　　(nth jdh1 jdzbb)

1 ~+ i- R- Z" Y7 }; E

　　(nth jdh2 jdzbb)

0 Q4 C& Q" U j* R

　　(nth jdh3 jdzbb)

; Z9 j# R0 G2 b! o 6 L5 d" g6 f0 j

　　(nth jdh4 jdzbb) "c")

2 K7 S) v% e/ Z% S0 Y' W+ o - @; ?, c1 P7 e

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

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

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

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

6 w c$ ]2 C! \

　　；标注单元号

, y; D$ V8 i6 e! r

　　(setq m 0)

3 q$ n4 _4 M R- A4 f7 p) o

　　(repeat jdzs

/ P/ A8 v6 D. B- J( k $ L {0 U0 @; Z+ K) r

　　(setq m (+ m 1))

) W1 [. X2 Q% a7 t7 Y, f( L0 ^

　　(setq bzdzb (nth m jdzbb))

8 G- v( a- [) O2 o; P

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

3 ~1 |# O* Q0 S# {" h/ Q. N; s! E

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

　　))；标注节点号

" k. l/ [! k. G. R

　　(close f)

　　)

+ n+ j6 B1 O0 q- |9 f

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

$ V# p* g& r9 ~6 e# q

　　3结语

3 S9 r1 y4 Z( p ^& N5 ]9 a8 T

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

% e* v; R/ a2 e8 ] T1 N+ @( i" n! N