OilfieldToolsNet/ResInsight
4
0
Fork 0
mirror of https://github.com/OPM/ResInsight.git synced 2025-02-25 18:55:39 -06:00
Code Issues Packages Projects Releases Wiki Activity

Flip Y to get local coordinate element nodes in correct order for abaqus.

Browse Source 
Bump version to dev.03
This commit is contained in:
jonjenssen
2023-12-18 00:48:25 +01:00
committed by jonjenssen
parent e67f09a4ff
commit 4aa5a38d63
3 changed files with 10 additions and 3 deletions
Download Patch File Download Diff File Expand all files Collapse all files

2
ApplicationLibCode/ReservoirDataModel/RigFaultReactivationModelGenerator.cpp
View File

@@ -154,7 +154,7 @@ void RigFaultReactivationModelGenerator::setupLocalCoordinateTransform()
{
auto [xNormal, yNormal] = modelLocalNormalsXY();
yNormal = -1.0 * yNormal;
// yNormal = yNormal * -1.0;
m_localCoordTransform = cvf::Mat4d::fromCoordSystemAxes( &xNormal, &yNormal, &cvf::Vec3d::Z_AXIS );
cvf::Vec3d center = m_startPosition * -1.0;

9
ApplicationLibCode/ReservoirDataModel/RigGriddedPart3d.cpp
View File

@@ -687,8 +687,15 @@ void RigGriddedPart3d::generateLocalNodes( const cvf::Mat4d transform )
{
m_localNodes.clear();
// need to flip the Y axis for the element corners to be in an acceptable order for abaqus and the IJK numbering algorithm in resinsight
cvf::Vec3d xAxis = { 1.0, 0.0, 0.0 };
cvf::Vec3d yAxis = { 0.0, -1.0, 0.0 };
cvf::Vec3d zAxis = { 0.0, 0.0, 1.0 };
cvf::Mat4d flipY = cvf::Mat4d::fromCoordSystemAxes( &xAxis, &yAxis, &zAxis );
for ( auto& node : m_nodes )
{
m_localNodes.push_back( node.getTransformedPoint( transform ) );
auto tn = node.getTransformedPoint( transform );
m_localNodes.push_back( tn.getTransformedPoint( flipY ) );
}
}