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

QQ375930253

　　2 .7应用软件前、后处理
　　由于软件升级滞后，一些常用软件后处理功能很弱或没有，可以利用LSP程序结合AutoCAD增加或简化前后处理功能。下面结合平面渗流分析程序STSE软件对其前、后处理功能的实现进行介绍：
　　(1)前处理
　　STSE为平面渗流有限元程序，单元划分和单元、节点编号工作量最大，可以借助其它通用有限元软件的前处理功能如Ansis、Algor、SAP84等进行初步处理，生成单元、节点编号和节点坐标，而后按照STSE数据文件的格式要求对数据文件进行编辑。
　　(2)后处理
　　后处理成果主要为：单元网格图，浸润线和等势线。单元网格图中有节点和单元编号，不同渗透系数的单元采用不同的颜色，以便根据网格图直观地判断数据文件中几何参数和材料特性的正误；
% E$ H1 N  C$ ^' v: s　　首先分析STSE的输出结果文件的格式，搜索并筛选其中主要参数如单元总数、节点总数和材料总数，把节点坐标和单元信息分别存储，通过对单元循环生成单元网格图。
5 u* ^2 i$ S! q9 ]  G* a　　(defun c:seepmesh()
* ^0 n" R$ T: y5 j4 S( B! R* J　　(if fn
4 P5 B4 t8 n! n- N3 j! h2 L　　(setq fn (getfiled "渗流结果文件名" fn "" 2))
. a1 x" C" @: T* k9 @" K　　(setq fn (getfiled "渗流结果文件名" "" "" 2)))
" I8 s$ F/ ~+ U9 `  P　　(setq f (open fn "r"))
5 L7 G3 u5 H) ]: [3 X　　(setq md (read-line f))
　　(while (/= (substr md 25 8) "单元总数")
/ J4 q$ d# x2 Y4 c$ q1 k/ \# ~" i) {3 J: l9 c　　(setq md (read-line f))
3 j" J# H1 C7 l+ W　　)；定位单元总数
  c7 Q3 M- f  T; D　　(setq dyzs (atoi (substr md 58 5)))；读单元总数
8 }8 r/ f. |) B3 J, i　　(setq md (read-line f))
% \+ ]/ @& F7 _9 j4 V) ?　　(setq jdzs (atoi (substr md 58 5)))；读节点总数
" A0 B/ G5 S1 d" |$ }  z$ Z: A　　(setq md (read-line f))
　　(setq clh (atoi (substr md 58 5)))；读材料总数
1 K0 [# Y) r" k. k, b9 B7 V　　(setq m 0)
0 f9 u/ p7 N6 }: ]* {　　(repeat clh
2 O* d2 Q; G+ L1 Q5 w　　(setq m (+ 1 m))
　　(setq tcm (strcat "zclh" (itoa m)))
6 m. G, Z  J: t% }7 P4 [) `　　(command "layer" "m" tcm "c" (itoa m) tcm "")
+ Y5 \: E, c1 I) @1 ], Y3 F　　)；按材料种类生成图层名称
　　(setq m 0 n 0)
　　(repeat dyzs；对单元循环，dycfb中存放单元信息
. q; k, @  e$ G: k9 }0 p# l　　(setq m (+ m 1))
( k2 ?' y% ]; l' w8 S　　(setq clh (nth 0 (nth m dycfb)))
　　(setq jdh1 (nth 1 (nth m dycfb)))
　　(setq jdh2 (nth 2 (nth m dycfb)))
　　(setq jdh3 (nth 3 (nth m dycfb)))
　　(setq jdh4 (nth 4 (nth m dycfb)))
　　(command "layer" "s" (strcat "zclh" (itoa clh)) "")
% {$ s/ j5 [6 ]: h* E( W2 W　　(command "pline"；绘制单元网格
( E' X  K8 d4 n9 D- Q　　(nth jdh1 jdzbb)
: e# O2 V1 a* W4 }　　(nth jdh2 jdzbb)
4 r# }8 h6 m5 Y9 D+ e% A. C$ Z' K　　(nth jdh3 jdzbb)
　　(nth jdh4 jdzbb) "c")
  `$ T+ g( [' k6 L2 l0 v　　(setq bzdzb (mapcar '+ (nth jdh1 jdzbb) (nth jdh2 jdzbb)
$ c: m1 _* i$ y5 d　　(nth jdh3 jdzbb) (nth jdh4 jdzbb)))
; W8 G" n3 E: n9 F6 X　　(setq bzdzb (list (/ (car bzdzb) 4.0) (/ (cadr bzdzb) 4.0)))
　　(command "text" "m" bzdzb "0.5" "0" (itoa m))
　　；标注单元号
　　(setq m 0)
+ }. d8 b- S/ V# v: c9 f( Y+ v4 Z　　(repeat jdzs
6 N+ e& o. q6 C) K" X　　(setq m (+ m 1))
　　(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 V; l1 y$ Q8 V# j% P* t　　))；标注节点号
　　(close f)
( R: r, ]1 e/ ?　　)
! e0 }" R/ I  [7 y5 {. M2 n. g8 `　　限于篇幅，等势线和浸润线生成程序不再赘述。同样，对工程设计中遇到的其它软件如：STAB边坡稳定分析、SAP84的输出文件均可进行类似后处理工作。
　　3结语
0 v% a' d" }7 G2 `　　AutoCAD软件应用已经在设计单位普及，VisualLisp的出现对于提高Lisp语言的编程效率和AutoCAD应用水平具有重要意义。本文涉及的几个实例已经在实际工作中发挥很大作用，充分说明在工程设计当中结合本专业要求开发一些简单实用的Lisp小程序可以有效减轻设计人员的劳动强度，使AutoCAD不仅仅是绘图工具，而真正成为设计人员得心应手的高效设计平台。
5 O& T6 r' {' Q& J- y文章关键词：