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

HEATS

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

$ ~0 d5 L7 u$ Z `% D 5 n8 p! o" f1 `) o3 m

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

( j4 `: K5 s4 f6 E; F

　　(1)前处理

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

8 T& g% p* m. |

　　(2)后处理

& B4 `, [$ J* G& A# e

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

# ^3 U9 u6 \9 m' D. b9 P0 V' A 6 j) q Q8 h' F, o6 m, |; l& J

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

, |- z) p( [5 j ( j' O# k; U4 D/ E1 n

　　(defun c:seepmesh()

5 d# o' O( |2 u( P0 t* {

　　(if fn

$ h6 t( t" o- ?. ]- C4 |8 T2 r ( F* M( l5 V# u9 ^+ Y

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

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

7 L- W+ B/ A7 f0 p! k1 q

　　(setq f (open fn "r"))

* y8 D1 }; F$ x" x1 K' v1 E : U0 `, k+ O$ l8 b5 P/ @

　　(setq md (read-line f))

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

9 v0 b' ]# e5 b, Z. U

　　(setq md (read-line f))

　　)；定位单元总数

2 ?2 ?. e9 L+ f/ v

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

! |8 G" D) J: b! b

　　(setq md (read-line f))

2 c# c5 ?9 h. P: V" | * H" U2 v/ @$ D5 r W

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

　　(setq md (read-line f))

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

( P( _7 R) `* r! v5 O" t

　　(setq m 0)

6 L$ k" y' O/ y0 @: n0 }7 y. }

　　(repeat clh

　　(setq m (+ 1 m))

. }9 c7 x# `% u9 Y , ]1 x, Z; t, n& @) N& N

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

/ f# H! i! J% `% g! u/ P

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

# ~2 m. L2 O/ [" t4 \! a! v6 R/ g

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

: m/ f2 m4 Y3 k! I1 _1 S

　　(setq m 0 n 0)

3 W: C( d; l( E M; F3 T 7 P& n6 J1 z9 |+ \. `, U( e

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

. s1 f% U' h9 c, C/ g

　　(setq m (+ m 1))

n! \. o) N: b : C3 i9 Z+ ^+ v' {6 o6 b9 T

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

8 J2 h7 Q0 o B( d

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

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

8 `, x8 B! R% b" p: h# F

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

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

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

1 U" L) p1 I3 _9 w7 V; G I8 A, u

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

7 o- J; t& S, p* D1 | ) N& y/ A& i1 Y8 s+ O7 M ]4 m

　　(nth jdh1 jdzbb)

/ [+ w# S$ I1 c) o6 \ & j& o# w6 ?: D' W

　　(nth jdh2 jdzbb)

' o4 X: n5 z; ~7 s 8 [. I$ ]- D/ U8 x5 U6 V

　　(nth jdh3 jdzbb)

$ x" A) \1 @4 o" Z$ C

　　(nth jdh4 jdzbb) "c")

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

: \% \- B8 u) @* _, }1 I3 r

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

! z* L! j) r# T/ k 0 d3 \; e$ C: h6 w! g

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

& V4 r8 J, J2 ~3 Q @) U @2 |

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

3 A% ]! k6 ~4 D# |+ N; B1 s

　　；标注单元号

* |! a: C' N2 b# s& s$ U! o1 J7 {

　　(setq m 0)

1 S0 C! b% X2 u! h% N5 ^

　　(repeat jdzs

) g' w! V, @ L7 K

　　(setq m (+ m 1))

! Q& |( t3 e7 o# R( F7 u 5 d" e6 r/ u5 c9 f" @; x

　　(setq bzdzb (nth m jdzbb))

# U9 v3 m3 k) Q) L/ k. M e2 K

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

2 a( F7 R$ x+ J1 }" a2 c

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

1 I% D4 S3 {' n8 c6 P# o ' R) J6 C! _; t

　　))；标注节点号

3 P# c: v3 T% C. A7 A 1 h$ t j4 [' E

　　(close f)

7 S* m8 q% E; n6 K

　　)

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

* w( q# k' B& |! u

　　3结语

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

0 }* J4 J- b0 i8 T! [ 4 n5 m/ E- k0 Z: |: ?/ ^