FEBio/FESlidingInterface.h

// FESlidingInterface.h: interface for the FESlidingInterface class.
//

#if !defined(AFX_FESLIDINGINTERFACE_H__742CFDED_B4BF_47AF_B673_1D26FE03F934__INCLUDED_)
#define AFX_FESLIDINGINTERFACE_H__742CFDED_B4BF_47AF_B673_1D26FE03F934__INCLUDED_

#if _MSC_VER > 1000
#pragma once
#endif // _MSC_VER > 1000

#include "FEContactInterface.h"
#include "FEContactSurface.h"

//-----------------------------------------------------------------------------
class FESlidingSurface : public FEContactSurface
{
public:
        FESlidingSurface(FEMesh* pm=0) : FEContactSurface(pm) { m_NQ.Attach(this); }

        bool Init();

        void ShallowCopy(FESlidingSurface& s)
        {
                Lm  = s.Lm;
                gap = s.gap;
                zero(m_pme);
                Lt  = s.Lt;
                m_Ln = s.m_Ln;
        }

        void Update() {}

        FEElement* FindMasterSegment(vec3d& x, vec3d& q, vec2d& r, bool& binit_nq, double tol);

        vec3d traction(int inode);

        void UpdateNormals();

        void Serialize(DumpFile& ar);

public:
        vector<double>                          gap;    
        vector<vec3d>                           nu;             
        vector<FESurfaceElement*>       m_pme;  
        vector<vec2d>                           rs;             
        vector<vec2d>                           rsp;    
        vector<double>                          Lm;             
        vector<mat2d>                           M;              
        vector<vec2d>                           Lt;             
        vector<double>                          off;    
        vector<double>                          eps;    
        vector<double>                          m_Ln;   

        FENNQuery               m_NQ;           
};

//-----------------------------------------------------------------------------


class FESlidingInterface : public FEContactInterface 
{
public:
        FESlidingInterface(FEModel* pfem);

        virtual ~FESlidingInterface(){}

        bool Init();

        void Update();

        void ProjectSurface(FESlidingSurface& ss, FESlidingSurface& ms, bool bupseg, bool bmove = false);

        void ShallowCopy(FEContactInterface& ci)
        {
                FESlidingInterface& si = dynamic_cast<FESlidingInterface&>(ci);
                m_ss.ShallowCopy(si.m_ss);
                m_ms.ShallowCopy(si.m_ms);
        }

        double Penalty() { return (m_pplc?m_eps*m_pplc->Value():m_eps);}

        virtual void ContactForces(vector<double>& F);

        virtual void ContactStiffness();

        virtual bool Augment(int naug);

        void Serialize(DumpFile& ar);

        void UpdateContactPressures();

protected:
        void CalcAutoPenalty(FESlidingSurface& s);

        void ContactNodalForce(int m, FESlidingSurface& ss, FESurfaceElement& mel, vector<double>& fe);

        void ContactNodalStiffness(int m, FESlidingSurface& ss, FESurfaceElement& mel, matrix& ke);

        void MapFrictionData(int inode, FESlidingSurface& ss, FESlidingSurface& ms, FESurfaceElement& sn, FESurfaceElement& so, vec3d& q);

private:
        FESlidingInterface(FESlidingInterface& si){}

public:
        FESlidingSurface        m_ss;   
        FESlidingSurface        m_ms;   

        bool                    m_btwo_pass;    

        int                             m_naugmax;      
        int                             m_naugmin;      
        double                  m_gtol;         
        double                  m_atol;         
        double                  m_ktmult;       
        double                  m_knmult;       

        double                  m_stol;         

        int                             m_nautopen;     
        double                  m_eps;          

        bool                    m_breloc;       

        double                  m_mu;           
        double                  m_epsf;         

        int                             m_nsegup;       
        int                             m_nplc;         
        FELoadCurve*    m_pplc;         

        DECLARE_PARAMETER_LIST();
};

#endif // !defined(AFX_FESLIDINGINTERFACE_H__742CFDED_B4BF_47AF_B673_1D26FE03F934__INCLUDED_)