opencascade计算面惯性矩的python代码

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

import math
from OCC.Core.BRepPrimAPI import BRepPrimAPI_MakeBox
from OCC.Core.GProp import GProp_GProps, GProp_PrincipalProps
from OCC.Core.BRepGProp import brepgprop_SurfaceProperties

# 创建一个长为10，宽为5，高为3的长方体
box = BRepPrimAPI_MakeBox(10, 5, 3).Shape()

# 计算面惯性矩
props = GProp_GProps()
brepgprop_SurfaceProperties(box, props)
principal_props = GProp_PrincipalProps(props)

# 输出结果
print("面积：", props.Mass())
print("质心：", props.CentreOfMass().Coord())
print("主惯性矩：", principal_props.I1(), principal_props.I2(), principal_props.I3())
print("主惯性轴：", principal_props.D1().Coord(), principal_props.D2().Coord(), principal_props.D3().Coord())

该代码使用BRepPrimAPI_MakeBox创建一个长方体，并使用brepgprop_SurfaceProperties计算其面惯性矩。结果包括面积，质心，主惯性矩和主惯性轴