FEBio/FEPressureLoad.h

#pragma once
#include "FESurfaceLoad.h"

//-----------------------------------------------------------------------------
class FEPressureLoad : public FESurfaceLoad
{
public:
        struct LOAD
        {
                double  s[4];           // nodal scale factors
                int             lc;                     // load curve

                LOAD() { s[0] = s[1] = s[2] = s[3] = 1.0; }
        };

public:
        FEPressureLoad(FESurface* ps, bool blinear = false) : FESurfaceLoad(ps) { m_blinear = blinear; }

        void create(int n) { m_PC.resize(n); }

/*      FEDomain* Clone()
        {
                FEPressureLoad* ps = new FEPressureLoad(m_surf.GetMesh());
                ps->m_PC = m_PC;
                return ps;
        }
*/

        LOAD& PressureLoad(int n) { return m_PC[n]; }

        void StiffnessMatrix(FESolver* psolver);

        void Residual(FESolver* psolver, vector<double>& R);

        void Serialize(DumpFile& ar);

protected:
        void PressureStiffness(FESurfaceElement& el, matrix& ke, vector<double>& tn);

        bool PressureForce(FESurfaceElement& el, vector<double>& fe, vector<double>& tn);

        bool LinearPressureForce(FESurfaceElement& el, vector<double>& fe, vector<double>& tn);

protected:
        bool                    m_blinear;      
        vector<LOAD>    m_PC;           
};