topic Re: Calculate volume using 3d Polyine in Visual LISP, AutoLISP and General Customization Forum

Calculate volume using 3d Polyine

cadgeek33 — Sun, 13 Jun 2021 16:42:32 GMT

Hi All,

I want to calculate cut volume for a pond. I have 2 Surfaces as 3d polylines with elevation in each vertices. Is there any lisp or routine I can follow rather than getting the volume using Civil 3D, because we have more than 500 similar ponds in this project. Sample dwg. and Image Attached for reference.

TIA.

Re: Calculate volume using 3d Polyine

Sea-Haven — Mon, 14 Jun 2021 01:50:50 GMT

I think it can be done found a volume of a 3dface formula, in the Civ3d build is example code about creating surfaces so pick 2 - 3d plines make 2 surfaces convert to 3d faces and read these calculating volume, I think you will still need a Vol1-vol2 calc. If it was possible to somehow run the CIV3D volume calc from a command point of view or reset the surface names or better still add surface names so do 500 shapes then export a results file.

Whilst you can get into CIV3D sometimes using VL no idea where the volume stuff is saved.

There is surface to 3d solid problem is it makes a shell of your surface.

Re: Calculate volume using 3d Polyine

john.uhden — Mon, 14 Jun 2021 20:20:01 GMT

I've never messed with solids, but I was thinking that a rough calculation could be made by using an average elevation for the top and and average elevation for the bottom, get each 2D area, and use the truncated prism method for computing the volume. I thinks it's V= (A1 +A2 + (A1*A2)^0.5)*(Z2-Z1)/3.

Re: Calculate volume using 3d Polyine

stevor — Mon, 14 Jun 2021 23:09:37 GMT

One way may be to make them into one 3D Solid,

and the command massprop .

Re: Calculate volume using 3d Polyine

hosneyalaa — Tue, 15 Jun 2021 05:31:27 GMT

see this link

https://forums.autodesk.com/t5/visual-lisp-autolisp-and-general/area-capacity-table/td-p/7417797?fbclid=IwAR10033g2rxoLzKN2WzDKFlGuvsF-jy7Mpq-S-9Cui3t1JyTGHLpJawhM84

Re: Calculate volume using 3d Polyine

Sea-Haven — Wed, 16 Jun 2021 01:03:57 GMT

Acpac2 looks interesting I think Johns idea of vol by prisms may be the way to go, I know I did intersect 3d faces years ago so could do create the cross section area of a slice, there is triangles by YMG at Cadtutor and pretty sure has cross section option, this would mean a program option not using CIV3D. The issue I see with Acpac is it needs contours but I did not test.

No matter what somehow must have a 3d mesh. Still thinking about it. Will try to find the civ3d make surface example.

Have a look at the readme.

Re: Calculate volume using 3d Polyine

Sea-Haven — Wed, 16 Jun 2021 01:06:39 GMT

I could be wrong maybe missing something but it makes a shell not a filled in solid as I suggested earlier there is a direct civ3d command to convert surface to solid.

Re: Calculate volume using 3d Polyine

cadgeek33 — Wed, 16 Jun 2021 04:53:53 GMT

Had tried Civil 3d, but its taking longer than usual. Any ways thanks alot.

Re: Calculate volume using 3d Polyine

cadgeek33 — Wed, 16 Jun 2021 04:55:09 GMT

hmmm.. seems like this will work for me. I would try. Thanks alot

Re: Calculate volume using 3d Polyine

cadgeek33 — Wed, 16 Jun 2021 04:58:14 GMT

Ya i am following the same concept until i get a Lisp. Thank You

Re: Calculate volume using 3d Polyine

cadgeek33 — Wed, 16 Jun 2021 04:59:57 GMT

This works with the 2D surfaces only. I am getting null value for 3D surfaces. If someone could modify the coding for the 3D surfaces that would be much appreciable. Thank you.

Re: Calculate volume using 3d Polyine

cadgeek33 — Wed, 16 Jun 2021 05:01:44 GMT

ya its making a shell only.

Re: Calculate volume using 3d Polyine

cadgeek33 — Wed, 16 Jun 2021 05:02:42 GMT

thanks alot, I would try this as well.

Re: Calculate volume using 3d Polyine

Sea-Haven — Wed, 16 Jun 2021 07:05:15 GMT

Really have no idea if this is correct

; https://stackoverflow.com/questions/5695322/calculate-volume-of-3d-mesh (setq p1 (list 1 1 1)) (setq p2 (list 2 13 15)) (setq p3 (list 20 15 8)) (command "3dface" p1 p2 p3 "" "") ; var v321 = p3.X * p2.Y * p1.Z (setq v321 (*(nth 0 p3)(nth 1 p2)(nth 2 p1))) ; var v231 = p2.X * p3.Y * p1.Z; (setq v231 (*(nth 0 p2) (nth 1 p3)(nth 2 p1))) ;var v312 = p3.X * p1.Y * p2.Z; (setq v312 (* (nth 0 p3)(nth 1 p1)(nth 2 p2))) ;var v132 = p1.X * p3.Y * p2.Z; (setq v132 (* (nth 0 p1)(nth 1 p3)(nth 2 p2))) ;var v213 = p2.X * p1.Y * p3.Z; (setq v213 (* (nth 0 p2)(nth 1 p1)(nth 2 p3))) ;var v123 = p1.X * p2.Y * p3.Z; (setq v123 (* (nth 0 p1)(nth 1 p2)(nth 2 p3))) (princ (setq vol (* (/ 1.0 6.0) (+ (- v321) v231 v312 (- v132) (- v213) v123))))

But could loop through all 3dfaces.

Re: Calculate volume using 3d Polyine

hak_vz — Wed, 16 Jun 2021 08:18:12 GMT

@john.uhden wrote:
I've never messed with solids, but I was thinking that a rough calculation could be made by using an average elevation for the top and and average elevation for the bottom, get each 2D area, and use the truncated prism method for computing the volume. I thinks it's V= (A1 +A2 + (A1*A2)^0.5)*(Z2-Z1)/3.

That's the way to calculate volume when two bases are parallel horizontal planes. In case as presented in example, where we have two inclined bases there are various options how to calculate volume. One is use truncated prism with heights z1 being average z of lower base and z2 average z of upper base. Other method is to add volumes of truncated prism and two wedges. Third is to multiply volume of truncated prism with some coefficient representing difference in inclined bases i.e let say 1.05.

Best option is to use method of parallel vertical sections and calculate volume using truncated prism formula, as we elaborated in one of previous posts. I have to work on code for this option, since requests like this are frequent.

Re: Calculate volume using 3d Polyine

Sea-Haven — Thu, 17 Jun 2021 02:42:04 GMT

If I can find the maths for this would use base as all at Z=0

Say with z=0 for base vol = (z1+z2+z3)/3 * Area of base at z=0 seems to be answer. Please disregard if I am wrong.

So would get total vol= sum vol triangles1 - sum vol triangles2.

If all wrong disregard, for 3dfaces need lisp a few free make TIN out there.

Only really suitable for holes or pyramids no cut/fill

Re: Calculate volume using 3d Polyine

hak_vz — Thu, 17 Jun 2021 06:27:45 GMT

@Sea-Haven

It's truncated pyramid and not prism (actually in some combinations one can use both). Volume is calculated using parallel sections with vertical plane. Areas from neighboring sections are looked as a bases of truncated pyramid and used to calculate partial volumes. that summed up give volume. I have started to write some code for this.

Re: Calculate volume using 3d Polyine

john.uhden — Fri, 18 Jun 2021 11:35:42 GMT

@hak_vz , @cadgeek33

I have been working on the multiple plane theory with tessellation lines connecting bottom to top and connecting the points at given Zs to form the equivalent of contours, then using the truncated prism method for the incremental volume between levels and then a total. It's very rough if the 3DPolys don't have many vertices, so I am going to add a tessellation factor to add more vertices (to the data, not to the polyline).

The volume is computed from the minimum elevation of the bottom to the minimum elevation of the top on the premise that this is for a pond volume, which would overflow if the water rises above its lowest outlet.

This could also output a stage-storage table (if I knew how to make tables). Maybe someone else can add that when I'm done. I can create a list that looks like:

'((z depth area incvol cumvol) ;; repeated)

Re: Calculate volume using 3d Polyine

hak_vz — Fri, 18 Jun 2021 15:38:59 GMT

@john.uhden wrote:
@hak_vz , @cadgeek33
It's very rough if the 3DPolys don't have many vertices, so I am going to add a tessellation factor to add more vertices (to the data, not to the polyline).

@john.uhden I like your idea. You can try to use ruled surface (rulesurf) created between upper and lower base 3dpolylines to obtain contour line points at some level Z. If ruled surface is exploded it creates 3dfaces from whom is relatively ease to points at level Z and create contour.

Re: Calculate volume using 3d Polyine

john.uhden — Sat, 19 Jun 2021 00:54:29 GMT

Thanks, @hak_vz,
but I already have the interpolation function, and I don't like using
commands unless I have to. It's all about manipulating data.