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

#914 Created stress tensor rotation method, including tests

Browse Source
This commit is contained in:
Jacob Støren 2016-10-17 15:14:45 +02:00
parent 323f7893cb
commit 941eed6625
3 changed files with 140 additions and 5 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

115
ApplicationCode/UnitTests/ScalarMapper-Test.cpp
View File

@ -121,3 +121,118 @@ TEST(OffshoreSphericalCoords, OffshoreSphericalCoords)
}
#include "../cafTensor/cafTensor3.h"
#include "cvfMath.h"
//--------------------------------------------------------------------------------------------------
///
//--------------------------------------------------------------------------------------------------
TEST(TensorRotation, TensorRotation)
{
{
caf::Ten3f orgT(1.0f, 0.5f, 0.2f, 0, 0, 0);
cvf::Mat3f rotMx = cvf::Mat3f::fromRotation(cvf::Vec3f(1, 0, 0), cvf::PI_F);
caf::Ten3f rotT = orgT.rotate(rotMx);
EXPECT_NEAR(1.0f, rotT[caf::Ten3f::SXX], 1e-4);
EXPECT_NEAR(0.5f, rotT[caf::Ten3f::SYY], 1e-4);
EXPECT_NEAR(0.2f, rotT[caf::Ten3f::SZZ], 1e-4);
EXPECT_NEAR(0.0f, rotT[caf::Ten3f::SXY], 1e-4);
EXPECT_NEAR(0.0f, rotT[caf::Ten3f::SZX], 1e-4);
EXPECT_NEAR(0.0f, rotT[caf::Ten3f::SYZ], 1e-4);
}
{
caf::Ten3f orgT(1.0f, 0.5f, 0.2f, 0, 0, 0);
cvf::Mat3f rotMx = cvf::Mat3f::fromRotation(cvf::Vec3f(1, 0, 0), 0.5*cvf::PI_F);
caf::Ten3f rotT = orgT.rotate(rotMx);
EXPECT_NEAR(1.0f, rotT[caf::Ten3f::SXX], 1e-4);
EXPECT_NEAR(0.2f, rotT[caf::Ten3f::SYY], 1e-4);
EXPECT_NEAR(0.5f, rotT[caf::Ten3f::SZZ], 1e-4);
EXPECT_NEAR(0.0f, rotT[caf::Ten3f::SXY], 1e-4);
EXPECT_NEAR(0.0f, rotT[caf::Ten3f::SZX], 1e-4);
EXPECT_NEAR(0.0f, rotT[caf::Ten3f::SYZ], 1e-4);
}
{
caf::Ten3f orgT(1.0f, 0.5f, 0.2f, 0, 0, 0);
cvf::Mat3f rotMx = cvf::Mat3f::fromRotation(cvf::Vec3f(0, 0, 1), 0.5*cvf::PI_F);
caf::Ten3f rotT = orgT.rotate(rotMx);
EXPECT_NEAR(0.5f, rotT[caf::Ten3f::SXX], 1e-4);
EXPECT_NEAR(1.0f, rotT[caf::Ten3f::SYY], 1e-4);
EXPECT_NEAR(0.2f, rotT[caf::Ten3f::SZZ], 1e-4);
EXPECT_NEAR(0.0f, rotT[caf::Ten3f::SXY], 1e-4);
EXPECT_NEAR(0.0f, rotT[caf::Ten3f::SZX], 1e-4);
EXPECT_NEAR(0.0f, rotT[caf::Ten3f::SYZ], 1e-4);
}
{
caf::Ten3f orgT(1.0f, 0.5f, 0.2f, 0, 0, 0);
cvf::Mat3f rotMx = cvf::Mat3f::fromRotation(cvf::Vec3f(0, 0, 1), 0.25*cvf::PI_F);
caf::Ten3f rotT = orgT.rotate(rotMx);
EXPECT_NEAR(0.75f,rotT[caf::Ten3f::SXX], 1e-4);
EXPECT_NEAR(0.75f,rotT[caf::Ten3f::SYY], 1e-4);
EXPECT_NEAR(0.2f, rotT[caf::Ten3f::SZZ], 1e-4);
EXPECT_NEAR(0.25f,rotT[caf::Ten3f::SXY], 1e-4);
EXPECT_NEAR(0.0f, rotT[caf::Ten3f::SZX], 1e-4);
EXPECT_NEAR(0.0f, rotT[caf::Ten3f::SYZ], 1e-4);
}
{
caf::Ten3f orgT(0.75f, 0.75f, 0.2f, 0.25, 0, 0);
cvf::Mat3f rotMx = cvf::Mat3f::fromRotation(cvf::Vec3f(0, 0, 1), -0.25*cvf::PI_F);
caf::Ten3f rotT = orgT.rotate(rotMx);
EXPECT_NEAR(1.0f, rotT[caf::Ten3f::SXX], 1e-4);
EXPECT_NEAR(0.5f, rotT[caf::Ten3f::SYY], 1e-4);
EXPECT_NEAR(0.2f, rotT[caf::Ten3f::SZZ], 1e-4);
EXPECT_NEAR(0.0f, rotT[caf::Ten3f::SXY], 1e-4);
EXPECT_NEAR(0.0f, rotT[caf::Ten3f::SZX], 1e-4);
EXPECT_NEAR(0.0f, rotT[caf::Ten3f::SYZ], 1e-4);
}
{
caf::Ten3f orgT(1.0f, 0.5f, 0.2f, 0, 0, 0);
cvf::Mat3f rotMx = cvf::Mat3f::fromRotation(cvf::Vec3f(1, 1, 1), 0.2*cvf::PI_F);
caf::Ten3f rotT = orgT.rotate(rotMx);
EXPECT_NEAR(0.8320561f, rotT[caf::Ten3f::SXX], 1e-4);
EXPECT_NEAR(0.5584094f, rotT[caf::Ten3f::SYY], 1e-4);
EXPECT_NEAR(0.3095343f, rotT[caf::Ten3f::SZZ], 1e-4);
EXPECT_NEAR(0.2091861f, rotT[caf::Ten3f::SXY], 1e-4);
EXPECT_NEAR(-0.2258091f, rotT[caf::Ten3f::SZX], 1e-4);
EXPECT_NEAR(0.0166228f, rotT[caf::Ten3f::SYZ], 1e-4);
}
{
caf::Ten3f orgT(0.8320561f, 0.5584094f, 0.3095343f, 0.2091861f, 0.0166228f, -0.2258091f);
cvf::Mat3f rotMx = cvf::Mat3f::fromRotation(cvf::Vec3f(1, 1, 1), -0.2*cvf::PI_F);
caf::Ten3f rotT = orgT.rotate(rotMx);
EXPECT_NEAR(1.0f, rotT[caf::Ten3f::SXX], 1e-4);
EXPECT_NEAR(0.5f, rotT[caf::Ten3f::SYY], 1e-4);
EXPECT_NEAR(0.2f, rotT[caf::Ten3f::SZZ], 1e-4);
EXPECT_NEAR(0.0f, rotT[caf::Ten3f::SXY], 1e-4);
EXPECT_NEAR(0.0f, rotT[caf::Ten3f::SZX], 1e-4);
EXPECT_NEAR(0.0f, rotT[caf::Ten3f::SYZ], 1e-4);
}
}

7
Fwk/AppFwk/cafTensor/cafTensor3.h
View File

@ -46,10 +46,11 @@ public:
void setFromInternalLayout(S* tensorData);
void setFromAbaqusLayout(S* tensorData);
cvf::Vec3f calculatePrincipals(cvf::Vec3f principalDirections[3]);
float calculateVonMises();
};
cvf::Vec3f calculatePrincipals(cvf::Vec3f principalDirections[3]) const;
float calculateVonMises() const;
Tensor3 rotate(const cvf::Matrix3<S>& rotMx) const;
};
typedef Tensor3<float> Ten3f;

23
Fwk/AppFwk/cafTensor/cafTensor3.inl
View File

@ -168,7 +168,7 @@ cvf::Vec3d eigenVector3(const cvf::Mat3d& mx, double eigenValue, bool* computedO
/// The tensor must be laid out as follows: SXX, SYY, SZZ, SXY, SYZ, SZX
//--------------------------------------------------------------------------------------------------
template<typename S>
cvf::Vec3f Tensor3<S>::calculatePrincipals( cvf::Vec3f principalDirections[3])
cvf::Vec3f Tensor3<S>::calculatePrincipals( cvf::Vec3f principalDirections[3]) const
{
CVF_TIGHT_ASSERT(m_tensor);
@ -314,12 +314,31 @@ cvf::Vec3f Tensor3<S>::calculatePrincipals( cvf::Vec3f principalDirections[3])
///
//--------------------------------------------------------------------------------------------------
template< typename S>
float caf::Tensor3<S>::calculateVonMises()
float caf::Tensor3<S>::calculateVonMises() const
{
return (float) sqrt( ( (m_tensor[0]*m_tensor[0] + m_tensor[1]*m_tensor[1] + m_tensor[2]*m_tensor[2]) ) +
( -(m_tensor[0]*m_tensor[1] + m_tensor[1]*m_tensor[2] + m_tensor[0]*m_tensor[2]) ) +
( 3*(m_tensor[3]*m_tensor[3] + m_tensor[4]*m_tensor[4] + m_tensor[5]*m_tensor[5]) ) );
}
//--------------------------------------------------------------------------------------------------
/// Calculates Trot = rotMx*T*transpose(rotMx)
//--------------------------------------------------------------------------------------------------
template< typename S>
Tensor3<S> caf::Tensor3<S>::rotate(const cvf::Matrix3<S>& rotMx) const
{
cvf::Matrix3<S> tensor(m_tensor[SXX], m_tensor[SXY], m_tensor[SZX],
m_tensor[SXY], m_tensor[SYY], m_tensor[SYZ],
m_tensor[SZX], m_tensor[SYZ], m_tensor[SZZ]);
cvf::Matrix3<S> transposedRotMx = rotMx;
transposedRotMx.transpose();
cvf::Matrix3<S> rotatedTensor = rotMx * tensor * transposedRotMx;
return Tensor3(rotatedTensor(0,0), rotatedTensor(1,1), rotatedTensor(2,2),
rotatedTensor(1,0), rotatedTensor(1,2), rotatedTensor(0,2));
}
}