opencascade计算面惯性矩的函数

Opencascade提供了计算面惯性矩的函数，可以使用以下代码：

#include <gp_Pnt.hxx>
#include <gp_Vec.hxx>
#include <BRep_Tool.hxx>
#include <BRepGProp.hxx>
#include <GProp_GProps.hxx>

//获取面的惯性矩
void getInertia(const TopoDS_Face& face, gp_Pnt& centroid, gp_Vec& vect, double& ixx, double& iyy, double& izz)
{
    GProp_GProps props;
    BRepGProp::SurfaceProperties(face, props);

    //获取面的重心
    centroid = props.CentreOfMass();

    //获取面的法向量
    vect = props.Normal();

    //获取面的惯性矩
    ixx = props.Ix();
    iyy = props.Iy();
    izz = props.Iz();
}

该函数使用BRepGProp::SurfaceProperties计算面的惯性矩，参数face为面的TopoDS_Shape对象，centroid为计算得到的面的重心，vect为计算得到的面的法向量，ixx、iyy、izz为计算得到的面的惯性矩