An Object Finite Element Library

Equa_Electromagnetics.h

Go to the documentation of this file.

/*==============================================================================

                                    O  F  E  L  I

                           Object  Finite  Element  Library

  ==============================================================================

   Copyright (C) 1998 - 2008 Rachid Touzani

   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; Version 2 of the License.

   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, write to the :

   Free Software Foundation
   Inc., 59 Temple Place - Suite 330
   Boston, MA  02111-1307, USA

  ==============================================================================

  Definition of abstract class 'FE_Electromagnetics' for Finite Element Equations
                         to solve electromagnetic problems

  ==============================================================================*/


#ifndef __FE_ELECTROMAGNETICS_H
#define __FE_ELECTROMAGNETICS_H

#include "Equation.h"
using std::complex;
using std::abs;

namespace OFELI {

template<class T_, size_t NEN_, size_t NEE_, size_t NSN_, size_t NSE_> class Equa_Electromagnetics;

template<class T_, size_t NEN_, size_t NEE_, size_t NSN_, size_t NSE_>
class Equa_Electromagnetics : virtual public Equation<T_,NEN_,NEE_,NSN_,NSE_>
{

 public:

   using Equation<T_,NEN_,NEE_,NSN_,NSE_>::setMaterialProperty;
   using Equation<T_,NEN_,NEE_,NSN_,NSE_>::_theElement;
   using Equation<T_,NEN_,NEE_,NSN_,NSE_>::_theSide;

 protected:

    void MagneticPermeability(const double &mu) { _mu = mu; }

    void MagneticPermeability(const char *exp) { setMaterialProperty(exp,"Magnetic Permeability"); }

    void ElectricConductivity(const double &sigma) { _sigma = sigma; }

    void ElectricConductivity(const char *exp) { setMaterialProperty(exp,"Electric Conductivity"); }

    void ElectricResistivity(const double &rho) { _rho = rho; }

    void ElectricResistivity(const char *exp) { setMaterialProperty(exp,"Electric Resistivity"); }

    void setMaterial()
    {
      _sigma = _theElement->getMaterial()->ElectricConductivity();
      _rho   = _theElement->getMaterial()->ElectricResistivity();
      _mu    = _theElement->getMaterial()->MagneticPermeability();
    }
#ifndef DOXYGEN_SHOULD_SKIP_THIS
   Point<double> _x[NEN_];
   double        _mu, _sigma, _rho, _body_source, _bound_source;
#endif /* DOXYGEN_SHOULD_SKIP_THIS */

};

} /* namespace OFELI */

#endif

Copyright © 1998-2007 Rachid Touzani