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

HEATS

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

9 u: B' w; v& g7 }' a0 Z

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

* t( n9 Y$ N: x# @

　　(1)前处理

( B, T0 T7 M4 b/ L- ] `+ N

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

+ \9 k$ ]* z& [5 I e

　　(2)后处理

^! y' U8 [; M

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

( ]4 d- F7 ^8 m+ l" | 6 J/ M: i O: p: z7 W: m

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

　　(defun c:seepmesh()

　　(if fn

7 M" I2 }8 K5 V* w $ Z% U$ ?) V" K( S" @. j

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

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

6 k7 f# S: s+ B+ d

　　(setq f (open fn "r"))

2 K8 }$ S; {$ N

　　(setq md (read-line f))

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

: d, R* e! P! V' F7 n: t + B ~* T8 ?7 y/ T/ |2 C

　　(setq md (read-line f))

0 W6 o" O5 V$ w$ y% v6 [- w- O: d

　　)；定位单元总数

( Z# x1 T$ }: Y$ O) ?( R5 F

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

6 | }, Q. u! b# {8 {& a& w+ Z! k

　　(setq md (read-line f))

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

　　(setq md (read-line f))

4 M! f3 G8 F, X, x1 W) l

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

　　(setq m 0)

; o. }2 x$ M% k $ u/ X( T& z9 q4 v. @

　　(repeat clh

1 ~; C8 p' ?! X, u2 d& s' h& W

　　(setq m (+ 1 m))

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

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

. k3 T3 N |+ Q' |; k( P4 k: ]/ ^ Q

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

( [7 l) C% v: g3 g: z$ Q6 i% C + i5 Z) X2 u5 R+ b

　　(setq m 0 n 0)

& w1 t2 Q( B! D' N

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

$ n7 d1 r! N5 F# x9 `- B( Z ( c7 V6 h0 f0 u2 k8 U

　　(setq m (+ m 1))

- |- l" U5 Y+ K* Q; k* a

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

1 P9 Z) [+ i8 K! j) T3 ]8 ?9 j

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

0 R# c* _# ], |$ I7 q7 m; V& Z

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

, F: ]! C! Q: g! L0 P+ @! W o

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

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

' d A0 ^# V5 c" x/ Q

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

& a& h9 Y: W, ?' _

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

/ ]; C/ r5 e$ x

　　(nth jdh1 jdzbb)

9 w, s- y1 {7 n1 x+ U , c2 b4 E; ^" Q

　　(nth jdh2 jdzbb)

　　(nth jdh3 jdzbb)

$ O: C6 T6 Q7 P% ~( m" J

　　(nth jdh4 jdzbb) "c")

: N2 k3 q6 n( ]3 x5 ~' p$ d( i( Q

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

7 o `& J/ `$ w

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

5 H- `# P, U0 {- |7 _4 ^7 Y( ]

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

2 J: |* a4 j" Q7 X) C* }! |( ]% O : D8 X/ t1 @) b# d# C

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

* ^, [. s$ e( j. R$ c2 d% q' s

　　；标注单元号

. F( O/ @( e7 A) x

　　(setq m 0)

　　(repeat jdzs

' X, z+ g1 N: T* p' W . n. b9 W) m5 ]$ F5 x1 l2 k

　　(setq m (+ m 1))

9 M6 ^. Y' k5 a: ^& I

　　(setq bzdzb (nth m jdzbb))

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

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

4 H! ^$ R' L" @% S- z: h5 S

　　))；标注节点号

　　(close f)

　　)

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

　　3结语

- y* u; E o( l# K

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

. k5 s$ m9 u! v6 f & v& D2 O8 _" C- n' G! k( T+ i# H