OpenCASCADE C++ 代码计算面惯性矩

以下是使用 OpenCASCADE 计算面惯性矩的 C++ 代码示例：

#include <stdio.h>
#include <math.h>
#include <BRepGProp.hxx>
#include <BRepPrimAPI_MakeBox.hxx>

int main()
{
    // 创建一个长方体
    BRepPrimAPI_MakeBox box(10.0, 20.0, 30.0);

    // 计算面的重心和惯性矩
    GProp_GProps props;
    BRepGProp::SurfaceProperties(box.Shape(), props);

    // 输出结果
    printf("面的重心：(%f, %f, %f)\n", props.CentreOfMass().X(), props.CentreOfMass().Y(), props.CentreOfMass().Z());
    printf("面的质量：%f\n", props.Mass());
    printf("面的惯性矩Ixx：%f\n", props.MatrixOfInertia().Value(1, 1));
    printf("面的惯性矩Iyy：%f\n", props.MatrixOfInertia().Value(2, 2));
    printf("面的惯性矩Izz：%f\n", props.MatrixOfInertia().Value(3, 3));

    return 0;
}

在此示例中，我们创建了一个长方体并使用 BRepGProp::SurfaceProperties 函数计算了其面的重心和惯性矩。最后，我们输出了结果。请注意，props.MatrixOfInertia() 返回的是一个 3x3 的矩阵，我们使用 Value() 函数获取其对应位置上的值。