An Object Finite Element Library

EC2D1T3.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 class EC2D1T3
     for Eddy Current Problems in Two-Dimensions with a scalar magnetic field
                           using the 3-Node triangle

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


#ifndef __EC2D1T3_H
#define __EC2D1T3_H


#include "Equa_Electromagnetics.h"
#include "Triang3.h"
#include "Line2.h"


namespace OFELI {

class EC2D1T3 : public Equa_Electromagnetics<complex_t,3,3,2,2>
{

 public :

    EC2D1T3() { }

    EC2D1T3(const Element *el);

    EC2D1T3(const Side *side);

    EC2D1T3(const Element *el, const Vect<complex_t> &u, const double &time=0.);

    EC2D1T3(const Side *sd, const Vect<complex_t> &u, const double &time=0.);

    ~EC2D1T3();

    void Magnetic(double omega, double coef=1.);

    void Electric(double coef=1.);

    double Joule();

    complex_t IntegMF();

    complex_t IntegND(const Vect<complex_t> &h, int opt=GLOBAL_ARRAY);

    double VacuumArea();

 protected :

    void ElementSet(const Element *el);
    void SideSet(const Side *el);

 private :
   Triang3 *_tr;
   Line2   *_ln;
};

} /* namespace OFELI */

#endif

Copyright © 1998-2007 Rachid Touzani