Open CASCADE学习|曲面上一点的曲率及切平面

曲率（Curvature）是一个几何学的概念，用于描述一个物体的形状在某一点上的弯曲程度。在我们日常生活中，曲率与我们的生活息息相关，如道路的弯道、建筑物的拱形结构、自然界的山脉等等。了解曲率的概念和计算方法，对于工程设计、自然科学研究以及艺术创作等领域具有重要的实际意义。

曲率是曲线和曲面的重要性质，表征了几何的“弯曲程度”。曲率越大，越弯曲；曲率越小，越趋近于直线、平面。直线、平面的曲率为零。

在一定条件下，过曲面Σ上的某一点M的曲线有无数多条，每一条曲线在点M处有一条切线，在一定的条件下这些切线位于同一平面，称这个平面为曲面Σ在点M处的切平面（tangent plane）。点M叫做切点。

#include <Geom_ToroidalSurface.hxx>
#include <BRepBuilderAPI_MakeFace.hxx>
#include <BRepAdaptor_Surface.hxx>
#include <GeomLProp_SLProps.hxx>
#include <Geom_Line.hxx>
#include "Viewer.h"
​
int main(int argc, char* argv[])
{gp_Torus aTorus(gp_Ax3(), 8.0, 4.0);Geom_ToroidalSurface  torSur(aTorus);TopoDS_Face aTourF = BRepBuilderAPI_MakeFace(aTorus);gp_Vec DirTanV;gp_Vec DirTanU;gp_Vec DirZMain;gp_Pnt p0;torSur.D1(1.0, 1.0, p0, DirTanU, DirTanV);DirZMain = DirTanU;DirZMain.Cross(DirTanV);gp_Pln TPlane(p0, DirZMain); //求切平面TopoDS_Face tanF = BRepBuilderAPI_MakeFace(TPlane,-10,10,-10,10);//求曲率BRepAdaptor_Surface adaptor(aTourF);GeomLProp_SLProps aProp(adaptor.Surface().Surface(), 2, 0.000001);aProp.SetParameters(1.0, 1.0);
​gp_Dir Max; //主曲率方向gp_Dir Min; //主曲率方向Standard_Real maxCurvature = aProp.MaxCurvature();Standard_Real minCurvature = aProp.MinCurvature();aProp.CurvatureDirections(Max, Min);
​Geom_Line aMaxVec(p0, Max);Geom_Line aMinVec(p0, Min);//aProp.GaussianCurvature();//aProp.MeanCurvature();Viewer vout(50, 50, 500, 500);vout << aTourF;vout << tanF;vout.StartMessageLoop();return 0;
}
​