/home/dorival/mechsys/lib/models/linflow.h

Go to the documentation of this file.

/************************************************************************
 * MechSys - Open Library for Mechanical Systems                        *
 * Copyright (C) 2005 Dorival M. Pedroso, Raul Durand                   *
 * Copyright (C) 2009 Sergio Galindo                                    *
 *                                                                      *
 * This program is free software: you can redistribute it and/or modify *
 * it under the terms of the GNU General Public License as published by *
 * the Free Software Foundation, either version 3 of the License, or    *
 * any later version.                                                   *
 *                                                                      *
 * This program is distributed in the hope that it will be useful,      *
 * but WITHOUT ANY WARRANTY; without even the implied warranty of       *
 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the         *
 * GNU General Public License for more details.                         *
 *                                                                      *
 * You should have received a copy of the GNU General Public License    *
 * along with this program. If not, see <http://www.gnu.org/licenses/>  *
 ************************************************************************/

#ifndef MECHSYS_LINFLOW_H
#define MECHSYS_LINFLOW_H

// MechSys
#include <mechsys/models/model.h>
#include <mechsys/models/flowstate.h>

class LinFlow : public Model
{
public:
    // Constructor
    LinFlow (int NDim, SDPair const & Prms);

    // Methods
    void InitIvs   (SDPair const & Ini, State * Sta)                             const;
    void TgIncs    (State const * Sta, Vec_t & DGra, Vec_t & DVel, Vec_t & DIvs) const { DVel = -1.0*D*DGra; }
    void Stiffness (State const * Sta, Mat_t & TheD)                             const { TheD = D; }

    // Data
    double kx;
    double ky;
    double kz;
    Mat_t  D;
};




inline LinFlow::LinFlow (int NDim, SDPair const & Prms)
    : Model (NDim,Prms,/*niv*/0,"LinFlow")
{
    // parameters
    if (Prms.HasKey("k"))
    {
        double k = Prms("k");
        kx = k;
        ky = k;
        kz = k;
    }
    else
    {
        kx = Prms("kx");
        ky = Prms("ky");  if (NDim==3)
        kz = Prms("kz");
    }

    // stiffness
    D.change_dim (NDim,NDim);
    if (NDim==2)
    {
        D =  kx, 0.0,
            0.0,  ky;
    }
    else if (NDim==3)
    {
        D =  kx, 0.0, 0.0,
            0.0,  ky, 0.0,
            0.0, 0.0,  kz;
    }
}

inline void LinFlow::InitIvs (SDPair const & Ini, State * Sta) const
{
    FlowState * sta = static_cast<FlowState*>(Sta);
    sta->Init (Ini);
}




Model * LinFlowMaker(int NDim, SDPair const & Prms, Model const * O) { return new LinFlow(NDim,Prms); }

int LinFlowRegister()
{
    ModelFactory   ["LinFlow"] = LinFlowMaker;
    MODEL.Set      ("LinFlow", (double)MODEL.Keys.Size());
    MODEL_PRM_NAMES["LinFlow"].Resize(3);
    MODEL_PRM_NAMES["LinFlow"] = "kx", "ky", "kz";
    MODEL_IVS_NAMES["LinFlow"].Resize(0);
    return 0;
}

int __LinFlow_dummy_int = LinFlowRegister();


#endif // MECHSYS_LINFLOW_H