如何不全部利用CASS算平场土方体积猜想

树櫴希德

先用CASS生成的三角网，记得三角网中任何三角形都不能穿越地形特征线，要符合地形地貌
   然后用下面这个小代码将三角网（三维多段线转化成三维面  3DFACE）

1. (defun vxs (e / i v lst)
   
2.   (setq i 0)
   
3.   (while
   
4.     (setq v (vlax-curve-getpointatparam e (setq i (1+ i))))
   
5.      (setq lst (cons v lst))
   
6.   )
   
7.   (reverse lst))
   
8. ;;;;;;;;;;;;;;;
   
9. 
   
10. (setq ssa (ssget '((0 . "POLYLINE") (8 . "sjw"))))
    
11.                 (setq ii   0
    
12.                       no  0
    
13.                   )
    
14.                   (repeat (sslength ssa)
    
15.                        (setq en (ssname ssa ii)
    
16.                             ptb (vxs en)
    
17.           pzx (append pzx (list ptb))
    
18.            ii  (1+ ii)               )
    
19.        ; (setq pzx (list (car pt) (cadr pt) (caddr pt)))
    
20.                        
    
21.                   )
    
22.   ;(write-line (strcat pzx) fff)
    
23.                  ; (close fff)
    
24.                   ;(princ (strcat "\n坐标已存入"" wjm ""中"))
    
25.    
    
26.     ;(setvar "cmdecho" cm)
    
27.     (princ)
    
28. 
    
29. 
    
30. (foreach n pzx
    
31. 
    
32. (entmake  (list '(0 . "3dface") (cons 10 (nth 0 n))
    
33.     (cons 11 (nth 1 n))
    
34.                 (cons 12 (nth 2 n))
    
35.                 (cons 13 (nth 0 n))
    
36.                 '(8 . "0")
    
37.     '(62 . 3)
    
38.     )  )
    
39. 
    
40. 
    
41.   )

(defun vxs (e / i v lst)<br />   (setq i 0)<br />   (while<br />     (setq v (vlax-curve-getpointatparam e (setq i (1+ i))))<br />      (setq lst (cons v lst))<br />   )<br />   (reverse lst))<br /> ;;;;;;;;;;;;;;;<br /> <br /> (setq ssa (ssget '((0 . "POLYLINE") (8 . "sjw"))))<br />                 (setq ii   0<br />                       no  0<br />                   )<br />                   (repeat (sslength ssa)<br />                        (setq en (ssname ssa ii)<br />                             ptb (vxs en)<br />           pzx (append pzx (list ptb))<br />            ii  (1+ ii)               )<br />        ; (setq pzx (list (car pt) (cadr pt) (caddr pt)))<br />                        <br />                   )<br />   ;(write-line (strcat pzx) fff) <br />                  ; (close fff)<br />                   ;(princ (strcat "\n坐标已存入"" wjm ""中"))<br />    <br />     ;(setvar "cmdecho" cm)<br />     (princ)<br /> <br /> <br /> (foreach n pzx<br /> <br /> (entmake  (list '(0 . "3dface") (cons 10 (nth 0 n))<br />     (cons 11 (nth 1 n))<br />                 (cons 12 (nth 2 n))<br />                 (cons 13 (nth 0 n))<br />                 '(8 . "0")<br />     '(62 . 3)<br />     )  )<br /> <br /> <br />   )

然后用 CONVTOSURFACE命令 将三维面3DFACE转换为曲面（平面，分散的曲面SURFACE PLANE），再用命令UNION（CAD2008不行，要高版本，2016使用可以）将分散的曲面SURFACE PLANE合并成完整的曲面SURFACE
   以标高最低点Z值画LWPOLYLINE多段线围住三角网一圈，并用拉伸命令EXTRUDE拉伸，高度为最高Z值减去最低点Z值形成三维实体3D solid
  用剖切命令SLICE将三维实体剖切，选合并成完整的曲面SURFACE为参照，删除不要的剖切后实体，剩余的就是要算体积的实体
用MASSPROP命令计算实体体积  完成
    这只是我的小小想法， 不足之处， 还望谅解、 指点，谢谢

gzxl

支持一下，也就是半自动化的思路咯。

Atsai

三角网=>取所有三角网的最低z值=>每个三角网跟最低z值面形成三角柱并求体积=>加总所有三角柱体积
这样也可以得到楼主用手动所求的资料。

树櫴希德

Atsai 发表于 2022-10-21 20:29
三角网=>取所有三角网的最低z值=>每个三角网跟最低z值面形成三角柱并求体积=>加总所有三角柱体积
这样也可 ...

但是你得不到一个完整的3D SOLID