FEBio/FETriphasicDomain.h

#pragma once
#include "FEElasticSolidDomain.h"

//-----------------------------------------------------------------------------
class FETriphasicDomain : public FEElasticSolidDomain
{
public:
        FETriphasicDomain(FEMesh* pm, FEMaterial* pmat) : FEElasticSolidDomain(pm, pmat) { m_ntype = FE_TRIPHASIC_DOMAIN; }
        
        FEDomain* Clone()
        {
                FETriphasicDomain* pd = new FETriphasicDomain(m_pMesh, m_pMat);
                pd->m_Elem = m_Elem; pd->m_pMesh = m_pMesh; pd->m_Node = m_Node;
                return pd;
        }
        
        void InitElements();
        
        void StiffnessMatrix(FESolidSolver* psolver);
        
        void Residual(FESolidSolver* psolver, vector<double>& R);
        
        bool Initialize(FEModel& fem);
        
        // update stresses
        void UpdateStresses(FEModel& fem);
        
        void ElementStiffness(FEM& fem, int iel, matrix& ke) {
                FESolidElement& el = Element(iel);
                ElementTriphasicStiffness(fem, el, ke);
        }
        
protected:
        bool InternalFluidWork(FEM& fem, FESolidElement& elem, vector<double>& fe);
        
        bool InternalFluidWorkSS(FEM& fem, FESolidElement& elem, vector<double>& fe);
        
        bool InternalSoluteWork(FEM& fem, FESolidElement& elem, vector<double>& fe, const int ion);
        
        bool InternalSoluteWorkSS(FEM& fem, FESolidElement& elem, vector<double>& fe, const int ion);
        
        bool ElementTriphasicStiffness(FEM& fem, FESolidElement& el, matrix& ke);
        
        bool ElementTriphasicStiffnessSS(FEM& fem, FESolidElement& el, matrix& ke);
        
        void SolidElementStiffness(FEM& fem, FESolidElement& el, matrix& ke);
        
        void TriphasicMaterialStiffness(FEM& fem, FESolidElement& el, matrix& ke);
        
};