fct_taille_fem_solution_fonction_isotrope.cpp

Aller à la documentation de ce fichier.

//####//------------------------------------------------------------
//####//------------------------------------------------------------
//####//  MAGiC
//####//  Jean Christophe Cuilliere et Vincent FRANCOIS
//####//  Departement de Genie Mecanique - UQTR
//####//------------------------------------------------------------
//####//  MAGIC est un projet de recherche de l equipe ERICCA
//####//  du departement de genie mecanique de l Universite du Quebec a Trois Rivieres
//####//  http://www.uqtr.ca/ericca
//####//  http://www.uqtr.ca/
//####//------------------------------------------------------------
//####//------------------------------------------------------------
//####//
//####//       fct_taille_fem_solution_fonction_isotrope.cpp
//####//
//####//------------------------------------------------------------
//####//------------------------------------------------------------
//####//    COPYRIGHT 2000-2024
//####//  jeu 13 jun 2024 11:58:52 EDT
//####//------------------------------------------------------------
//####//------------------------------------------------------------
#include "gestionversion.h"
#include "mg_gestionnaire.h"
#include <math.h>
#include <fem_solution.h>
#include "fct_taille_fem_solution_fonction_isotrope.h"
 
 
double FCT_TAILLE_FEM_SOLUTION_DEFINITION_FONCTION_ISOTROPE::constante[100]={0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.};
 
 
double FCT_TAILLE_FEM_SOLUTION_DEFINITION_FONCTION_ISOTROPE::exponentielle(double x)
{
double coeff=FCT_TAILLE_FEM_SOLUTION_DEFINITION_FONCTION_ISOTROPE::constante[1];
double res = (exp(coeff*x))*FCT_TAILLE_FEM_SOLUTION_DEFINITION_FONCTION_ISOTROPE::constante[0];
return res;
}
double FCT_TAILLE_FEM_SOLUTION_DEFINITION_FONCTION_ISOTROPE::lineaire(double x)
{
double coeff=FCT_TAILLE_FEM_SOLUTION_DEFINITION_FONCTION_ISOTROPE::constante[1];
double res = (coeff*x) + FCT_TAILLE_FEM_SOLUTION_DEFINITION_FONCTION_ISOTROPE::constante[0];
return res;
}
 
 
FCT_TAILLE_FEM_SOLUTION_FONCTION_ISOTROPE::FCT_TAILLE_FEM_SOLUTION_FONCTION_ISOTROPE(MG_GESTIONNAIRE *gesttmp,FEM_MAILLAGE* femtmp,char *nomsolution):FCT_TAILLE_FEM_SOLUTION(gesttmp,femtmp)
{
strcpy(nom,nomsolution);
isotrope=true;
}
 
FCT_TAILLE_FEM_SOLUTION_FONCTION_ISOTROPE::FCT_TAILLE_FEM_SOLUTION_FONCTION_ISOTROPE(FCT_TAILLE_FEM_SOLUTION_FONCTION_ISOTROPE& mdd):FCT_TAILLE_FEM_SOLUTION(mdd)
{
strcpy(nom,mdd.nom);
isotrope=true;
}
 
FCT_TAILLE_FEM_SOLUTION_FONCTION_ISOTROPE::~FCT_TAILLE_FEM_SOLUTION_FONCTION_ISOTROPE()
{
  
}
 
 
void FCT_TAILLE_FEM_SOLUTION_FONCTION_ISOTROPE::construit_depuis_scalaire(FEM_SOLUTION *solbase,int champbase,double (*f)(double x))
{
sol=new FEM_SOLUTION(fem,1,nom,1,(char*)"Carte de taille",MAGIC::ENTITE_SOLUTION::ENTITE_NOEUD,MAGIC::TYPE_SOLUTION::TENSEUR);
gest->ajouter_fem_solution(sol);
sol->change_legende(0,"metrique");
LISTE_FEM_NOEUD::iterator it;
int i=0;
for (FEM_NOEUD* no=fem->get_premier_noeud(it);no!=NULL;no=fem->get_suivant_noeud(it))
    {
    double val;
    val=solbase->lire(i,champbase);
    double ecart_nodal=(*f)(val);
    sol->ecrire(1./ecart_nodal/ecart_nodal,i,0,0);
    sol->ecrire(0.,i,0,1);
    sol->ecrire(0.,i,0,2);
    sol->ecrire(0.,i,0,3);
    sol->ecrire(1./ecart_nodal/ecart_nodal,i,0,4);
    sol->ecrire(0.,i,0,5);
    sol->ecrire(0.,i,0,6);
    sol->ecrire(0.,i,0,7);
    sol->ecrire(1./ecart_nodal/ecart_nodal,i,0,8);
    i++;  
    }
isotrope=true;
numchamp=0;
sol->active_solution(numchamp);
}