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

HEATS

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

- k" u$ n. a2 C! ~

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

) t' _5 `! \0 b0 S8 i

　　(1)前处理

3 D9 p' K$ k* S$ c: T

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

8 N5 [$ q, p, U9 }1 r$ ]( F8 u % W. `8 e" w: J! U4 W N3 _

　　(2)后处理

+ T& C( D. K: b6 l- ~ ! ]0 k& x% c/ L

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

7 R' }5 v* A* G& s4 Y. m' F4 c( x+ U : Q2 _6 ?9 C' G3 e

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

　　(defun c:seepmesh()

0 A- _& |1 q) y3 g9 \

　　(if fn

3 r. U7 p* s9 ~, t( Y

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

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

　　(setq f (open fn "r"))

　　(setq md (read-line f))

$ B& l4 H3 l# g4 b7 J- l

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

　　(setq md (read-line f))

　　)；定位单元总数

; d0 N( c6 S* u d1 j3 n, Q5 \9 b

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

1 \7 o) c7 Z Z, T, X2 P

　　(setq md (read-line f))

E6 y% {8 t9 Y; F+ n. ?6 }

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

( v8 P A. L& V& _; ~/ H/ { ) G4 B9 v! U$ j3 ^" W( [% i

　　(setq md (read-line f))

* i/ a$ g% e( K ! E3 s. l! h6 W* }1 j

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

4 l9 ]2 u$ s8 L0 W0 p

　　(setq m 0)

3 |' h$ h4 V5 ]/ I7 B' H

　　(repeat clh

　　(setq m (+ 1 m))

: i! w/ h9 y/ Z* r: Z " ^6 K/ T5 `( D- T; w, E

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

/ T9 v# @/ P$ j . g' f* J/ M- c% D& w) w4 }

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

0 ~' A; r) a( l: h$ U0 E% [' e! t

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

/ d: d+ N6 s$ _6 z2 V s/ J

　　(setq m 0 n 0)

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

　　(setq m (+ m 1))

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

c x- f1 X: X9 G

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

4 S# I1 m; f1 C# S& \ 4 Z, \* S% U, \" P" Q

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

$ o* p g1 @- s2 ?# y

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

: Z- u4 w, s" S+ l' v8 K4 n

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

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

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

: A$ O$ o) L+ p) o

　　(nth jdh1 jdzbb)

3 {! k7 x7 Z' m* W* y

　　(nth jdh2 jdzbb)

2 f, M+ O# h1 q& j1 t4 n/ t8 H: @ # L( L, U$ x1 ^' O' T! _

　　(nth jdh3 jdzbb)

~3 _$ e8 u' [

　　(nth jdh4 jdzbb) "c")

" b8 H1 V1 n" E1 _2 }

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

: o. v6 l4 y8 F7 A

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

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

+ @* h3 ?- s; H+ \$ Z

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

　　；标注单元号

+ X0 Z/ ^- A8 ?/ Z- D1 [

　　(setq m 0)

, R% k: T% d! e + E" I4 y$ A$ N1 @. S7 U4 x

　　(repeat jdzs

; H& C9 P5 I9 x+ c8 G$ x* V+ ?2 `* A " ~# j6 U% U, w

　　(setq m (+ m 1))

　　(setq bzdzb (nth m jdzbb))

O" h9 o2 t6 L0 U) Z' B/ N% A

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

- p8 u, I8 _% G2 Q6 w. V, u: p6 A 3 D% Q% o0 g* \/ T" {" u) k/ O/ ?! @

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

2 B+ o& ]3 l5 B+ T( B8 s; C + P- F$ p9 Q9 s: d( K8 m; y

　　))；标注节点号

5 E5 L" A. t" S" a 7 O ?4 c& b1 w, G( l

　　(close f)

, @- `, f: `: c+ D2 I! Y : C. F$ t+ i b' d X2 |

　　)

4 }. a5 ^+ n1 A$ F4 L! |$ Y* X* x+ R

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

2 I; a; C/ K7 I6 R- ^1 S' T

　　3结语

|* i8 o, m, R% Z* k

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

' M! K' R& K4 c$ O) x) l" ` # x$ G, i- k4 W; c/ l) h9 h