volume_polyhedron
- compas.geometry.volume_polyhedron(polyhedron)[source]
Compute the volume of a polyhedron represented by a closed mesh.
- Parameters:
- polyhedrontuple[sequence[[float, float, float] |
compas.geometry.Point], sequence[sequence[int]]] The vertices and faces of the polyhedron.
- polyhedrontuple[sequence[[float, float, float] |
- Returns:
- float
The volume of the polyhedron.
Warning
The volume computed by this funtion is only correct if the polyhedron is convex, has planar faces, and is positively oriented (all face normals point outwards).
Notes
This implementation is based on the divergence theorem, the fact that the area vector is constant for each face, and the fact that the area of each face can be computed as half the length of the cross product of two adjacent edge vectors [1].
\begin{align} V = \int_{P} 1 &= \frac{1}{3} \int_{\partial P} \mathbf{x} \cdot \mathbf{n} \\ &= \frac{1}{3} \sum_{i=0}^{N-1} \int{A_{i}} a_{i} \cdot n_{i} \\ &= \frac{1}{6} \sum_{i=0}^{N-1} a_{i} \cdot \hat n_{i} \end{align}References
[1]Nurnberg, R. Calculating the area and centroid of a polygon in 2d. Available at: http://wwwf.imperial.ac.uk/~rn/centroid.pdf