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

HEATS

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

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

# {3 {0 B( h9 f

　　(1)前处理

- H: E% T' N% T: G; Y* B t0 T

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

　　(2)后处理

3 ]1 G q& z# L6 `: R

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

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

" ]& b3 P5 k2 A4 M B

　　(defun c:seepmesh()

$ X4 N. A9 Z3 g' ?! h/ @8 n ) ^, t) W8 y9 `& ~8 Z6 N# S9 f" q

　　(if fn

U' R) I/ V: z! ] $ W) `, ~! R& f, s n$ \9 z2 x8 v

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

3 S# v! w* l! i$ n% j% ~( @

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

　　(setq f (open fn "r"))

3 _; j5 M7 p* K$ e5 F

　　(setq md (read-line f))

# m6 P+ b1 k' t5 u6 T* y

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

　　(setq md (read-line f))

, R$ j5 ~$ s6 `; h) ^0 T( _

　　)；定位单元总数

6 {1 L" Q. }+ Z / l' I% j+ ]2 ~# l1 A( r! }

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

, s, N. n$ g3 ^* F 5 O% p9 k% J9 h% l

　　(setq md (read-line f))

3 x D5 L; t- @+ C% p . s6 R' g" T$ e# _: C4 s$ K# S

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

0 d& F' N/ _( ~) A/ p' j5 R8 w 5 ?# U2 |8 {3 j8 N# P, `+ O

　　(setq md (read-line f))

9 r6 s* k: c) d2 E1 _( f+ e * ?. H3 p( G+ v8 k1 r

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

　　(setq m 0)

2 H8 O6 ?2 D' ] K) W) x3 q

　　(repeat clh

7 Z) E" Y% S) w( F; S0 A/ Z: h- a7 T 5 l3 M+ k& U* e1 A5 w: O+ C9 H% f

　　(setq m (+ 1 m))

$ O7 F3 j5 k, T- q; Z9 e

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

) _7 _) \+ \- G

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

/ |8 a0 L6 f2 B4 `9 ?/ ^

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

8 y, M, z6 m" Y1 p

　　(setq m 0 n 0)

7 o; N2 h5 I6 U& e W4 i2 h6 g

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

+ `9 R/ g2 J: w( x

　　(setq m (+ m 1))

* Y* F2 p t, ~" ? z, a8 ]$ X # s- Q5 T! ?* r2 B

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

6 s. }$ J; r7 n$ V- ]) Q

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

- ] x, v. V3 ^/ r2 G z# B7 L

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

7 ]$ X% J1 h6 N& @0 z8 y. N$ a8 m

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

+ C% |; P! u p ) J: H% T+ b8 D: O

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

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

% K! `2 ?( ]1 w

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

\+ r- Z- [% D2 t2 ~ + `, _3 a1 }! A0 A1 b6 S

　　(nth jdh1 jdzbb)

% m9 ^; U# O4 @1 e. w , a( U% A( j& o1 ^8 W( R2 B$ w

　　(nth jdh2 jdzbb)

& `8 |; G8 p5 N2 `) R

　　(nth jdh3 jdzbb)

8 s2 z9 ]( S3 Q }( s 0 c6 e) J9 \# w$ s( q% ?

　　(nth jdh4 jdzbb) "c")

7 n& h' B6 b& Q1 \' S0 u # e& B; Y6 a: F3 ?* t( i2 ?. a

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

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

* Q% e u0 [, } d% J

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

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

! I r x: u' y" I/ Q% D" y

　　；标注单元号

　　(setq m 0)

　　(repeat jdzs

* ]4 Z' ]) g0 O2 I1 D6 `+ l

　　(setq m (+ m 1))

- p: H% W- f, B3 D

　　(setq bzdzb (nth m jdzbb))

/ f, E2 n+ t: Q5 o$ N# Y" U

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

+ J) L( f t# ]: v7 R' z! G

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

. R8 I5 b$ ?& K; @4 _$ o) u

　　))；标注节点号

# N+ M: Y9 H' D# O

　　(close f)

! A* \1 j* v% t; E 1 Z9 S8 K/ }% N5 l( ]# C- E

　　)

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

　　3结语

6 _7 Y9 ^2 X& |! t7 t5 j! x

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

: N c$ z8 Z0 C& z