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

move assignment of linStress into MechContainer

Browse Source
This commit is contained in:
Arne Morten Kvarving 2025-02-05 11:35:41 +01:00
parent 779a225c9d
commit 6809268d12
3 changed files with 19 additions and 8 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

13
opm/simulators/flow/MechContainer.cpp
View File

@ -117,6 +117,19 @@ assignDisplacement(const unsigned globalDofIdx,
this->disp_[2][globalDofIdx] = disp[2];
}
template<class Scalar>
void MechContainer<Scalar>::
assignLinStress(const unsigned globalDofIdx,
const Dune::FieldVector<Scalar,6>& linStress)
{
this->linstressXX_[globalDofIdx] = linStress[0];
this->linstressYY_[globalDofIdx] = linStress[1];
this->linstressZZ_[globalDofIdx] = linStress[2];
this->linstressYZ_[globalDofIdx] = linStress[3];
this->linstressXZ_[globalDofIdx] = linStress[4];
this->linstressXY_[globalDofIdx] = linStress[5];
}
template<class Scalar>
void MechContainer<Scalar>::
assignPotentialForces(const unsigned globalDofIdx,

3
opm/simulators/flow/MechContainer.hpp
View File

@ -60,6 +60,9 @@ public:
const Scalar pressForce,
const Scalar tempForce);
void assignLinStress(const unsigned globalDofIdx,
const Dune::FieldVector<Scalar,6>& linStress);
void assignStrain(const unsigned globalDofIdx,
const Dune::FieldVector<Scalar,6>& strain);

11
opm/simulators/flow/OutputBlackoilModule.hpp
View File

@ -215,6 +215,9 @@ public:
this->mech_.assignDisplacement(globalDofIdx,
model.disp(globalDofIdx, /*include_fracture*/true));
this->mech_.assignLinStress(globalDofIdx,
model.linstress(globalDofIdx));
this->mech_.assignPotentialForces(globalDofIdx,
model.mechPotentialForce(globalDofIdx),
model.mechPotentialPressForce(globalDofIdx),
@ -227,14 +230,6 @@ public:
this->mech_.assignStress(globalDofIdx,
model.stress(globalDofIdx, /*include_fracture*/true));;
const auto linstress = model.linstress(globalDofIdx);
this->mech_.linstressXX_[globalDofIdx] = linstress[Voigt::XX];
this->mech_.linstressYY_[globalDofIdx] = linstress[Voigt::YY];
this->mech_.linstressZZ_[globalDofIdx] = linstress[Voigt::ZZ];
this->mech_.linstressXY_[globalDofIdx] = linstress[Voigt::XY];
this->mech_.linstressXZ_[globalDofIdx] = linstress[Voigt::XZ];
this->mech_.linstressYZ_[globalDofIdx] = linstress[Voigt::YZ];
// is the tresagii stress which make rock fracture
const auto fracstress = model.fractureStress(globalDofIdx);
this->mech_.fracstressXX_[globalDofIdx] = fracstress[Voigt::XX];