An Object Finite Element Library

Material.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 Material to select materials

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

#ifndef __MATERIAL_H
#define __MATERIAL_H


#include <string>

#include <iostream>
using std::ostream;
using std::endl;
using std::cerr;

#include <iomanip>
using std::setw;

#include "OFELI_Config.h"
#include "Point.h"
#include "FFI.h"
#include "fparser.h"

using std::string;

extern FunctionParser _theParser;

namespace OFELI {

enum {
   BY_VALUE,
   BY_POSITION,
   BY_FIELD
};

#ifndef DOXYGEN_SHOULD_SKIP_THIS
struct Prop {
   bool    exist;
   int     type;
   double  value;
   string  fp_xyzt;
};
#endif /* DOXYGEN_SHOULD_SKIP_THIS */


#ifndef DOXYGEN_SHOULD_SKIP_THIS
struct ErrorInMaterial {
   void Message(const char *file, size_t line, int code, int p1, const string &p2);
};
#endif /* DOXYGEN_SHOULD_SKIP_THIS */

class Material
{

 public:

    Material();

    Material(const Material &m);

    ~Material();

    int set(int m, const string &name, int type);

#ifndef DOXYGEN_SHOULD_SKIP_THIS
    void scan();
    void scanXML();
#endif /* DOXYGEN_SHOULD_SKIP_THIS */

    string getName(int m) const;

    int getCode(size_t i) const;

    size_t getNbMat() const;

    void setCode(int m);

    int check(int c);

    double Density();

    double Density(const Point<double> &x, double t);

    double SpecificHeat();

    double SpecificHeat(const Point<double> &x, double t);

    double ThermalConductivity();

    double ThermalConductivity(const Point<double> &x, double t);

    double MeltingTemperature();

    double MeltingTemperature(const Point<double> &x, double t);

    double EvaporationTemperature();

    double EvaporationTemperature(const Point<double> &x, double t);

    double ThermalExpansion();

    double ThermalExpansion(const Point<double> &x, double t);

    double LatentHeatForMelting();

    double LatentHeatForMelting(const Point<double> &x, double t);

    double LatentHeatForEvaporation();

    double LatentHeatForEvaporation(const Point<double> &x, double t);

    double DielectricConstant();

    double DielectricConstant(const Point<double> &x, double t);

    double ElectricConductivity();

    double ElectricConductivity(const Point<double> &x, double t);

    double ElectricResistivity();

    double ElectricResistivity(const Point<double> &x, double t);

    double MagneticPermeability();

    double MagneticPermeability(const Point<double> &x, double t);

    double Viscosity();

    double Viscosity(const Point<double> &x, double t);

    double YoungModulus();

    double YoungModulus(const Point<double> &x, double t);

    double PoissonRatio();

    double PoissonRatio(const Point<double> &x, double t);

    double Property(int i);

    double Property(int i, const Point<double> &x, double t);

    Material & operator=(const Material &m)
    {
       _nb_mat = m._nb_mat;
       _path = m._path;
       for (size_t i=0; i<_nb_mat; i++) {
          _code[i] = m._code[i];
          _mat[i] = m._mat[i];
       }
       _ff = m._ff;
       set_kw();
       return *this;
    }

    friend class XMLParser;

 private:

   int _state, _type;
   size_t _nb_mat, _index_mat;
   string _mat[MAX_NB_MATERIALS], _name, _path;
   int _code[MAX_NB_MATERIALS];
   FFI _ff;
   valarray<string> _kw;
   void getValue(size_t i);
   void getFunction(size_t i);
   void getTable() { }
   double _data[10];
   int read(size_t);
   ErrorInMaterial _e;

   void set_kw();
   double getProperty(Prop &prop);
   double getProperty(Prop &prop, const Point<double> &x, const double &t);
   Prop _density[MAX_NB_MATERIALS];
   Prop _specific_heat[MAX_NB_MATERIALS];
   Prop _thermal_conductivity[MAX_NB_MATERIALS];
   Prop _melting_temperature[MAX_NB_MATERIALS];
   Prop _evaporation_temperature[MAX_NB_MATERIALS];
   Prop _thermal_expansion[MAX_NB_MATERIALS];
   Prop _latent_heat_melting[MAX_NB_MATERIALS];
   Prop _latent_heat_evaporation[MAX_NB_MATERIALS];
   Prop _dielectric_constant[MAX_NB_MATERIALS];
   Prop _electric_conductivity[MAX_NB_MATERIALS];
   Prop _electric_resistivity[MAX_NB_MATERIALS];
   Prop _magnetic_permeability[MAX_NB_MATERIALS];
   Prop _viscosity[MAX_NB_MATERIALS];
   Prop _young_modulus[MAX_NB_MATERIALS];
   Prop _poisson_ratio[MAX_NB_MATERIALS];
   Prop _prop[MAX_NB_MATERIALS][10];
};


    ostream& operator<<(ostream& s, const Material &m);

} /* namespace OFELI */

#endif

Copyright © 1998-2007 Rachid Touzani