ForkMaster/smesh
1
0
Fork 0
mirror of https://git.salome-platform.org/gitpub/modules/smesh.git synced 2025-02-06 02:34:17 +05:00
Code Issues Packages Projects Releases Wiki Activity

PAL13615 (EDF PAL 315/31 GEOM SMESH : meshing of a "5 edges quadrangle")

Browse Source 
add GetLength(const SMESH_Mesh* aMesh, const double edgeLength)
This commit is contained in:
eap 2007-02-20 07:10:31 +00:00
parent a7051ca684
commit 1863381f98
2 changed files with 129 additions and 79 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

103
src/StdMeshers/StdMeshers_AutomaticLength.cxx
View File

@ -97,29 +97,51 @@ void StdMeshers_AutomaticLength::SetFineness(double theFineness)
} }
} }
//================================================================================ namespace {
/*!
//================================================================================
/*!
* \brief Return pointer to TopoDS_TShape * \brief Return pointer to TopoDS_TShape
* \param theShape - The TopoDS_Shape * \param theShape - The TopoDS_Shape
* \retval inline const TopoDS_TShape* - result * \retval inline const TopoDS_TShape* - result
*/ */
//================================================================================ //================================================================================
inline const TopoDS_TShape* getTShape(const TopoDS_Shape& theShape) inline const TopoDS_TShape* getTShape(const TopoDS_Shape& theShape)
{ {
return theShape.TShape().operator->(); return theShape.TShape().operator->();
} }
//================================================================================
/*! //================================================================================
/*!
* \brief computes segment length by S0 and edge length
*/
//================================================================================
const double a14divPI = 14. / PI;
inline double segLength(double S0, double edgeLen, double minLen )
{
// PAL10237
// S = S0 * f(L/Lmin) where f(x) = 1 + (2/Pi * 7 * atan(x/5) )
// =>
// S = S0 * ( 1 + 14/PI * atan( L / ( 5 * Lmin )))
return S0 * ( 1. + a14divPI * atan( edgeLen / ( 5 * minLen )));
}
//================================================================================
/*!
* \brief Compute segment length for all edges * \brief Compute segment length for all edges
* \param theMesh - The mesh * \param theMesh - The mesh
* \param theTShapeToLengthMap - The map of edge to segment length * \param theTShapeToLengthMap - The map of edge to segment length
*/ */
//================================================================================ //================================================================================
static void computeLengths( SMESHDS_Mesh* aMesh, void computeLengths( SMESHDS_Mesh* aMesh,
map<const TopoDS_TShape*, double> & theTShapeToLengthMap) map<const TopoDS_TShape*, double> & theTShapeToLengthMap,
{ double & theS0,
double & theMinLen)
{
theTShapeToLengthMap.clear(); theTShapeToLengthMap.clear();
TopoDS_Shape aMainShape = aMesh->ShapeToMesh(); TopoDS_Shape aMainShape = aMesh->ShapeToMesh();
@ -147,17 +169,17 @@ static void computeLengths( SMESHDS_Mesh* aMesh,
// image attached to PAL10237 // image attached to PAL10237
// NbSeg // NbSeg
// ^ // ^
// | // |
// 10|\ // 10|\
// | \ // | \
// | \ // | \
// | \ // | \
// 5| -------- // 5| --------
// | // |
// +------------> // +------------>
// 1 10 Lmax/Lmin // 1 10 Lmax/Lmin
const int NbSegMin = 5, NbSegMax = 10; // on axis NbSeg const int NbSegMin = 5, NbSegMax = 10; // on axis NbSeg
const double Lrat1 = 1., Lrat2 = 10.; // on axis Lmax/Lmin const double Lrat1 = 1., Lrat2 = 10.; // on axis Lmax/Lmin
@ -171,18 +193,37 @@ static void computeLengths( SMESHDS_Mesh* aMesh,
MESSAGE( "S0 = " << S0 << ", Lmin = " << Lmin << ", Nbseg = " << (int) NbSeg); MESSAGE( "S0 = " << S0 << ", Lmin = " << Lmin << ", Nbseg = " << (int) NbSeg);
// Compute segments length for all edges // Compute segments length for all edges
// S = S0 * f(L/Lmin) where f(x) = 1 + (2/Pi * 7 * atan(x/5) )
// =>
// S = S0 * ( 1 + 14/PI * atan( L / ( 5 * Lmin )))
const double a14divPI = 14. / PI, a5xLmin = 5 * Lmin;
map<const TopoDS_TShape*, double>::iterator tshape_length = theTShapeToLengthMap.begin(); map<const TopoDS_TShape*, double>::iterator tshape_length = theTShapeToLengthMap.begin();
for ( ; tshape_length != theTShapeToLengthMap.end(); ++tshape_length ) for ( ; tshape_length != theTShapeToLengthMap.end(); ++tshape_length )
{ {
double & L = tshape_length->second; double & L = tshape_length->second;
L = S0 * ( 1. + a14divPI * atan( L / a5xLmin )); L = segLength( S0, L, Lmin );
} }
theS0 = S0;
theMinLen = Lmin;
}
}
//=============================================================================
/*!
* \brief Computes segment length for an edge of given length
*/
//=============================================================================
double StdMeshers_AutomaticLength::GetLength(const SMESH_Mesh* theMesh,
const double theEdgeLength)
throw(SALOME_Exception)
{
if ( !theMesh ) throw SALOME_Exception(LOCALIZED("NULL Mesh"));
SMESHDS_Mesh* aMeshDS = const_cast< SMESH_Mesh* > ( theMesh )->GetMeshDS();
if ( theMesh != _mesh )
{
computeLengths( aMeshDS, _TShapeToLength, _S0, _minLen );
_mesh = theMesh;
}
double L = segLength( _S0, theEdgeLength, _minLen );
return L / (theCoarseConst + theFineConst * _fineness);
} }
//============================================================================= //=============================================================================
@ -203,7 +244,7 @@ double StdMeshers_AutomaticLength::GetLength(const SMESH_Mesh* theMesh,
if ( theMesh != _mesh ) if ( theMesh != _mesh )
{ {
SMESHDS_Mesh* aMeshDS = const_cast< SMESH_Mesh* > ( theMesh )->GetMeshDS(); SMESHDS_Mesh* aMeshDS = const_cast< SMESH_Mesh* > ( theMesh )->GetMeshDS();
computeLengths( aMeshDS, _TShapeToLength ); computeLengths( aMeshDS, _TShapeToLength, _S0, _minLen );
_mesh = theMesh; _mesh = theMesh;
} }

11
src/StdMeshers/StdMeshers_AutomaticLength.hxx
View File

@ -51,9 +51,18 @@ public:
StdMeshers_AutomaticLength(int hypId, int studyId, SMESH_Gen * gen); StdMeshers_AutomaticLength(int hypId, int studyId, SMESH_Gen * gen);
virtual ~ StdMeshers_AutomaticLength(); virtual ~ StdMeshers_AutomaticLength();
/*!
* \brief Computes segment for a given edge
*/
double GetLength(const SMESH_Mesh* aMesh, const TopoDS_Shape& anEdge) double GetLength(const SMESH_Mesh* aMesh, const TopoDS_Shape& anEdge)
throw(SALOME_Exception); throw(SALOME_Exception);
/*!
* \brief Computes segment length for an edge of given length
*/
double GetLength(const SMESH_Mesh* aMesh, const double edgeLength)
throw(SALOME_Exception);
/*! /*!
* \brief Set Fineness * \brief Set Fineness
* \param theFineness - The Fineness value [0.0-1.0], * \param theFineness - The Fineness value [0.0-1.0],
@ -89,7 +98,7 @@ public:
protected: protected:
std::map<const TopoDS_TShape*, double> _TShapeToLength; std::map<const TopoDS_TShape*, double> _TShapeToLength;
const SMESH_Mesh* _mesh; const SMESH_Mesh* _mesh;
double _fineness; double _fineness, _S0, _minLen;
}; };
#endif #endif