Contenuto principale

volume

R2026b

Compute volume of 3-D mesh elements

Description

V = volume(mesh) returns the volume V of the entire mesh.

example

[V,VE] = volume(mesh) also returns a row vector VE containing volumes of each individual element of the mesh.

example

V = volume(mesh,elements) returns the combined volume of the specified elements of the mesh.

example

Examples

collapse all

Generate a 3-D mesh and find its volume.

Create a PDE model.

model = createpde;

Import and plot the geometry.

importGeometry(model,"BracketWithHole.stl");
pdegplot(model)

3-D bracket geometry with a circular hole

Generate a mesh and plot it.

mesh = generateMesh(model);
figure
pdemesh(model)

Tetrahedral mesh generated for the bracket geometry

Compute the volume of the entire mesh.

mv = volume(mesh)
mv = 
8.0244e-04

Generate a 3-D mesh and find the volume of each element.

Create a PDE model.

model = createpde;

Import and plot the geometry.

importGeometry(model,"BracketWithHole.stl");
pdegplot(model)

3-D bracket geometry with a circular hole

Generate a mesh and plot it.

mesh = generateMesh(model);
figure
pdemesh(model)

Tetrahedral mesh generated for the bracket geometry

Compute the volume of the entire mesh and the volume of each individual element of the mesh. Display the volumes of the first 5 elements.

[va,vi] = volume(mesh);
vi(1:5)
ans = 1×5
10-6 ×

    0.3801    0.4653    0.1624    0.4345    0.1729

Find the combined volume of a group of elements of a 3-D mesh.

Create a PDE model.

model = createpde;

Import and plot the geometry.

importGeometry(model,"BracketWithHole.stl");
pdegplot(model)

Figure contains an axes object. The axes object contains 6 objects of type quiver, text, patch, line.

Generate a mesh and plot it.

mesh = generateMesh(model);
figure
pdemesh(model)

Figure contains an axes object. The hidden axes object contains 5 objects of type quiver, text, patch.

Evaluate the shape quality of the mesh elements and find the elements with the quality values less than 0.5.

Q = meshQuality(mesh);
elemIDs = find(Q < 0.5);

Compute the total volume of these elements.

mv05 = volume(mesh,elemIDs)
mv05 = 
1.4865e-07

Find how much of the total mesh volume belongs to these elements. Return the result as a percentage.

mv05_percent = mv05/volume(mesh)*100
mv05_percent = 
0.0185

Input Arguments

collapse all

Mesh description, specified as an FEMesh object.

Element IDs, specified as a positive integer or a matrix of positive integers.

Example: [10 68 81 97 113 130 136 164]

Output Arguments

collapse all

Volume of the entire mesh or the combined volume of the specified elements of the mesh, returned as a positive number.

Volume of individual elements, returned as a row vector of positive numbers.

Version History

Introduced in R2018a