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

HEATS

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

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

　　(1)前处理

" L' M' G& N0 ~9 i

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

: v3 X+ M- N3 G( t4 s

　　(2)后处理

9 ?6 [( e3 e, Z! C

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

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

　　(defun c:seepmesh()

6 ^# j: @8 G8 g0 s1 |

　　(if fn

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

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

　　(setq f (open fn "r"))

( S1 E# u5 _9 d* i. P. U, @9 X" b

　　(setq md (read-line f))

" A+ M9 |/ I1 v& n" ?* A

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

. N9 r# A1 u Y1 h7 ?0 }

　　(setq md (read-line f))

$ a0 [8 T7 T k4 ^$ Q3 ^ # \- j8 T4 O7 p, ]0 e# w. m

　　)；定位单元总数

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

　　(setq md (read-line f))

: Q! W, z, g8 g4 v" q * N; O+ z1 X$ \* u) d9 a( z- z

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

J4 _: a. s2 c2 M 0 z9 U, X. X% s8 t U7 _% H

　　(setq md (read-line f))

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

　　(setq m 0)

9 r9 U. r& N- S" s

　　(repeat clh

4 S3 Z) Z8 Q7 ?4 ?: p; ^' d3 Y

　　(setq m (+ 1 m))

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

' b8 S ^1 ]8 F0 m9 I

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

3 {# ?. P: g5 G) d! N( l

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

% q. z! E4 v7 A& k4 D

　　(setq m 0 n 0)

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

$ G6 w( o8 s* U" |+ m

　　(setq m (+ m 1))

0 U/ P4 q0 Z1 t. |7 N* q

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

) l$ J9 a1 G- v# b4 s3 ~4 m/ ~

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

8 O8 c5 O8 s; v7 K

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

3 M' f; j( J4 Z% p ) Z; e5 h7 B( F5 r' \8 h

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

5 N, u. \ c- [8 `4 ~2 G, d8 b 1 v) \2 U9 T) z% b* s8 Z" N

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

! |0 @( N8 Q, ?

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

$ m1 p, z2 \8 D2 i

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

　　(nth jdh1 jdzbb)

　　(nth jdh2 jdzbb)

6 z p6 I, k" g; q. `- h v8 b, n9 ^1 R

　　(nth jdh3 jdzbb)

　　(nth jdh4 jdzbb) "c")

, s/ Y" k6 M6 i9 K. J7 e1 G

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

& e8 X6 x4 ^; [8 ]! C1 H# }* X$ X4 z 3 k) F* z( l/ \

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

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

' q* j/ {* h X* }$ T9 } 9 i: F% t. j8 x# i! `" z

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

　　；标注单元号

$ H2 a9 v9 }$ G5 u! \# Q 3 I7 d/ Z) P! R2 F/ P ?1 Z2 Z/ C

　　(setq m 0)

, R" q& v0 l3 Q5 q1 o

　　(repeat jdzs

& Z. Z; s8 r( T/ f& _# z

　　(setq m (+ m 1))

' c( R; }% }6 x) A, C$ {

　　(setq bzdzb (nth m jdzbb))

) J- k9 w6 p. K8 q' q; f

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

$ Q& _# j$ x" c

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

, K5 \" ]* A% x9 D4 w( t$ O

　　))；标注节点号

0 W* \/ T$ q- s* J- n- l

　　(close f)

　　)

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

+ A3 t j5 C; a: Z( t

　　3结语

1 }5 P! `( _( h5 d* I" `

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

" L J; g: y& m9 |7 Y8 M; { $ D' j0 [, x( ?( ?3 U8 J6 L