geometrie/src/fem_maillage.cpp

//------------------------------------------------------------
//------------------------------------------------------------
//  MAGiC
//  Jean Christophe Cuilli�e et Vincent FRANCOIS
//  D�artement de G�ie M�anique - UQTR
//------------------------------------------------------------
//  Le projet MAGIC est un projet de recherche du d�artement
//  de g�ie m�anique de l'Universit�du Qu�ec �
//  Trois Rivi�es
//  Les librairies ne peuvent �re utilis�s sans l'accord
//  des auteurs (contact : francois@uqtr.ca)
//------------------------------------------------------------
//------------------------------------------------------------
//
//       mg_maillage.cpp
//
//------------------------------------------------------------
//------------------------------------------------------------
//  COPYRIGHT 2000
//  Version du 02/03/2006 �11H22
//------------------------------------------------------------
//------------------------------------------------------------


#include "gestionversion.h"
#include <fstream>
#include <iostream>
#include <string.h>
#include <math.h>


#include "fem_maillage.h"
#include "fem_graph_noeud.h"
#include "mg_gestionnaire.h"
#include "mt_gestionnaire.h"
#include "mt_propriete.h"
#include "mt_materiau.h"
#include "ot_decalage_parametre.h"
#include "tpl_map_entite.h"
#include "tpl_liste_entite.h"
//#include "message.h"


FEM_MAILLAGE::FEM_MAILLAGE(MG_GEOMETRIE* geo,MG_MAILLAGE* mai,int deg):MG_IDENTIFICATEUR(),mg_geometrie(geo),maillage(mai),degre(deg),deforme(0)
{
}

FEM_MAILLAGE::FEM_MAILLAGE(unsigned long num,MG_GEOMETRIE* geo,MG_MAILLAGE* mai,int deg):MG_IDENTIFICATEUR(num),mg_geometrie(geo),maillage(mai),degre(deg),deforme(0)
{
}

FEM_MAILLAGE::~FEM_MAILLAGE()
{
supprimer_tout_fem_tetra();
supprimer_tout_fem_triangle();
supprimer_tout_fem_segment();
supprimer_tout_fem_noeud();
}

FEM_MAILLAGE* FEM_MAILLAGE::dupliquer(MG_GESTIONNAIRE* gest)
{
FEM_MAILLAGE* mai=new FEM_MAILLAGE(mg_geometrie,maillage,degre);
gest->ajouter_fem_maillage(mai);
long decalage=mai->get_id()-this->get_id();
LISTE_ENTITE::const_iterator i;
for (i=lst_entite.begin();i!=lst_entite.end();i++)
     {
        FEM_ELEMENT_MAILLAGE& ele=(FEM_ELEMENT_MAILLAGE&)(*(*i));
        ele.dupliquer(mai,decalage);
     }
return mai;
}

MG_GEOMETRIE* FEM_MAILLAGE::get_mg_geometrie(void)
{
return mg_geometrie;
}

MG_MAILLAGE* FEM_MAILLAGE::get_mg_maillage(void)
{
return maillage;
}

int FEM_MAILLAGE::get_degre(void)
{
return degre;
}

bool FEM_MAILLAGE::contient ( MG_IDENTIFICATEUR * id )
{
        return ( lst_entite.find ( id ) != lst_entite.end() );
}

void FEM_MAILLAGE::change_gestionnaire(MG_GESTIONNAIRE* mggest)
{
gest=mggest;
}

// ENTITE FEM_NOEUD
FEM_NOEUD*  FEM_MAILLAGE::ajouter_fem_noeud(MG_NOEUD* mgnoeud,unsigned long num)
{
FEM_NOEUD* femnoeud;
if (num==0) femnoeud=new FEM_NOEUD(mgnoeud); else femnoeud=new FEM_NOEUD(num,mgnoeud);
int resultat = ajouter_fem_noeud(femnoeud);
if (resultat==FAIL)
   {
      delete femnoeud;
      return NULL;
   }
return femnoeud;
}

int FEM_MAILLAGE::ajouter_fem_noeud(FEM_NOEUD *femnoeud)
{
gest->recherche_bonid(*femnoeud);
MG_IDENTIFICATEUR *id=femnoeud;
std::pair<LISTE_ENTITE::iterator,bool> p=lst_entite.insert(id);
if (!p.second)
   {
//      afficheur << ENTITEID << id->get_id() << DUPLIQUE << enderr;
      return FAIL;
   }

std::pair<const unsigned long,FEM_NOEUD*> tmp(femnoeud->get_id(),femnoeud);
lst_fem_noeud.insert(tmp);
return OK;
}


FEM_NOEUD* FEM_MAILLAGE::get_fem_noeudid(unsigned long num)
{
LISTE_FEM_NOEUD::iterator i=lst_fem_noeud.find(num);
if (i==lst_fem_noeud.end())
   {
//      afficheur << INEXISTE << enderr;
      return NULL;
   }
return ((*i).second);
}


FEM_NOEUD* FEM_MAILLAGE::get_fem_noeud(unsigned int num)
{
if (!(num<lst_fem_noeud.size()))
   {
//      afficheur << INEXISTE << enderr;
      return NULL;
   }
LISTE_FEM_NOEUD::iterator i = lst_fem_noeud.begin();
std::advance(i, num);

return ((*i).second);
}

FEM_NOEUD* FEM_MAILLAGE::get_premier_noeud(LISTE_FEM_NOEUD::iterator &  it)
{
it = lst_fem_noeud.begin();
if (it == lst_fem_noeud.end())
    return NULL;
return it->second ;
}

FEM_NOEUD* FEM_MAILLAGE::get_suivant_noeud(LISTE_FEM_NOEUD::iterator &  it)
{
it++;
if(it==lst_fem_noeud.end())
    return NULL;
return it->second;
}

unsigned int FEM_MAILLAGE::get_nb_fem_noeud(void)
{
return lst_fem_noeud.size();
}


int FEM_MAILLAGE::supprimer_fem_noeudid(unsigned long num)
{
FEM_NOEUD* femnoeud=get_fem_noeudid(num);
if (femnoeud==NULL)
   {
//      afficheur << INEXISTE2 << enderr;
      return FAIL;
   }
MG_IDENTIFICATEUR* id=femnoeud;
LISTE_ENTITE::iterator i=lst_entite.find(id);
lst_entite.erase(i);
LISTE_FEM_NOEUD::iterator j=lst_fem_noeud.find(num);
lst_fem_noeud.erase(j);
delete femnoeud;
return OK;
}


int FEM_MAILLAGE::supprimer_fem_noeud(unsigned int num)
{
FEM_NOEUD* femnoeud=get_fem_noeud(num);
if (femnoeud==NULL)
   {
//      afficheur << INEXISTE2 << enderr;
      return FAIL;
   }
MG_IDENTIFICATEUR* id=femnoeud;
LISTE_ENTITE::iterator i=lst_entite.find(id);
lst_entite.erase(i);
LISTE_FEM_NOEUD::iterator j=lst_fem_noeud.begin();
for (unsigned int k=0;k<num;k++) j++;
lst_fem_noeud.erase(j);
delete femnoeud;
return OK;
}


void FEM_MAILLAGE::supprimer_tout_fem_noeud(void)
{
while (get_nb_fem_noeud()!=0)
 {
    LISTE_FEM_NOEUD::iterator j=lst_fem_noeud.begin();
    FEM_NOEUD* femnoeud=(*j).second;
    MG_IDENTIFICATEUR* id=femnoeud;
    LISTE_ENTITE::iterator i=lst_entite.find(id);
    lst_entite.erase(i);
    lst_fem_noeud.erase(j);
    delete femnoeud;
 }
}


// ENTITE FEM_SEGMENT


int FEM_MAILLAGE::ajouter_fem_segment(FEM_SEGMENT *femsegment)
{
gest->recherche_bonid(*femsegment);
MG_IDENTIFICATEUR *id=femsegment;
std::pair<LISTE_ENTITE::iterator,bool> p=lst_entite.insert(id);
if (!p.second)
   {
//      afficheur << ENTITEID << id->get_id() << DUPLIQUE << enderr;
      return FAIL;
   }

std::pair<const unsigned long,FEM_SEGMENT*> tmp(femsegment->get_id(),femsegment);
lst_fem_segment.insert(tmp);
return OK;
}

FEM_SEGMENT* FEM_MAILLAGE::get_fem_segmentid(unsigned long num)
{
LISTE_FEM_SEGMENT::iterator i=lst_fem_segment.find(num);
if (i==lst_fem_segment.end())
   {
//      afficheur << INEXISTE << enderr;
      return NULL;
   }
return ((*i).second);
}

unsigned int FEM_MAILLAGE::get_nb_fem_segment(void)
{
return lst_fem_segment.size();
}


FEM_SEGMENT* FEM_MAILLAGE::get_fem_segment(unsigned int num)
{
if (!(num<lst_fem_segment.size()))
   {
//      afficheur << INEXISTE << enderr;
      return NULL;
   }
LISTE_FEM_SEGMENT::iterator i=lst_fem_segment.begin();
for (unsigned long j=0;j<num;j++) i++;
return ((*i).second);
}


int FEM_MAILLAGE::supprimer_fem_segmentid(unsigned long num)
{
FEM_SEGMENT* femsegment=get_fem_segmentid(num);
if (femsegment==NULL)
   {
//      afficheur << INEXISTE2 << enderr;
      return FAIL;
   }
MG_IDENTIFICATEUR* id=femsegment;
LISTE_ENTITE::iterator i=lst_entite.find(id);
lst_entite.erase(i);
LISTE_FEM_SEGMENT::iterator j=lst_fem_segment.find(num);
lst_fem_segment.erase(j);
delete femsegment;
return OK;
}

int FEM_MAILLAGE::supprimer_fem_segment(unsigned int num)
{
FEM_SEGMENT* femsegment=get_fem_segment(num);
if (femsegment==NULL)
   {
//      afficheur << INEXISTE2 << enderr;
      return FAIL;
   }
MG_IDENTIFICATEUR* id=femsegment;
LISTE_ENTITE::iterator i=lst_entite.find(id);
lst_entite.erase(i);
LISTE_FEM_SEGMENT::iterator j=lst_fem_segment.find(num);
lst_fem_segment.erase(j);
delete femsegment;
return OK;
}

void FEM_MAILLAGE::supprimer_tout_fem_segment(void)
{
while (get_nb_fem_segment()!=0)
 {
    LISTE_FEM_SEGMENT::iterator j=lst_fem_segment.begin();
    FEM_SEGMENT* femsegment=(*j).second;
    MG_IDENTIFICATEUR* id=femsegment;
    LISTE_ENTITE::iterator i=lst_entite.find(id);
    lst_entite.erase(i);
    lst_fem_segment.erase(j);
    delete femsegment;
 }
}


FEM_SEGMENT* FEM_MAILLAGE::get_premier_segment(LISTE_FEM_SEGMENT::iterator &  it)
{
it = lst_fem_segment.begin();
if (it == lst_fem_segment.end())
    return NULL;
return it->second ;
}

FEM_SEGMENT* FEM_MAILLAGE::get_suivant_segment(LISTE_FEM_SEGMENT::iterator &  it)
{
it++;
if(it==lst_fem_segment.end())
    return NULL;
return it->second;
}


// ENTITE FEM_TRIANGLE


int FEM_MAILLAGE::ajouter_fem_triangle(FEM_TRIANGLE *femtriangle)
{
gest->recherche_bonid(*femtriangle);
MG_IDENTIFICATEUR *id=femtriangle;
std::pair<LISTE_ENTITE::iterator,bool> p=lst_entite.insert(id);
if (!p.second)
   {
//      afficheur << ENTITEID << id->get_id() << DUPLIQUE << enderr;
      return FAIL;
   }

std::pair<const unsigned long,FEM_TRIANGLE*> tmp(femtriangle->get_id(),femtriangle);
lst_fem_triangle.insert(tmp);
return OK;
}

FEM_TRIANGLE* FEM_MAILLAGE::get_fem_triangleid(unsigned long num)
{
LISTE_FEM_TRIANGLE::iterator i=lst_fem_triangle.find(num);
if (i==lst_fem_triangle.end())
   {
//      afficheur << INEXISTE << enderr;
      return NULL;
   }
return ((*i).second);
}

unsigned int FEM_MAILLAGE::get_nb_fem_triangle(void)
{
return lst_fem_triangle.size();
}


FEM_TRIANGLE* FEM_MAILLAGE::get_fem_triangle(unsigned int num)
{
if (!(num<lst_fem_triangle.size()))
   {
//      afficheur << INEXISTE << enderr;
      return NULL;
   }
LISTE_FEM_TRIANGLE::iterator i=lst_fem_triangle.begin();
for (unsigned long j=0;j<num;j++) i++;
return ((*i).second);
}


int FEM_MAILLAGE::supprimer_fem_triangleid(unsigned long num)
{
FEM_TRIANGLE* femtriangle=get_fem_triangleid(num);
if (femtriangle==NULL)
   {
//      afficheur << INEXISTE2 << enderr;
      return FAIL;
   }
MG_IDENTIFICATEUR* id=femtriangle;
LISTE_ENTITE::iterator i=lst_entite.find(id);
lst_entite.erase(i);
LISTE_FEM_TRIANGLE::iterator j=lst_fem_triangle.find(num);
lst_fem_triangle.erase(j);
delete femtriangle;
return OK;
}

int FEM_MAILLAGE::supprimer_fem_triangle(unsigned int num)
{
FEM_TRIANGLE* femtriangle=get_fem_triangle(num);
if (femtriangle==NULL)
   {
//      afficheur << INEXISTE2 << enderr;
      return FAIL;
   }
MG_IDENTIFICATEUR* id=femtriangle;
LISTE_ENTITE::iterator i=lst_entite.find(id);
lst_entite.erase(i);
LISTE_FEM_TRIANGLE::iterator j=lst_fem_triangle.find(num);
lst_fem_triangle.erase(j);
delete femtriangle;
return OK;
}

void FEM_MAILLAGE::supprimer_tout_fem_triangle(void)
{
while (get_nb_fem_triangle()!=0)
 {
    LISTE_FEM_TRIANGLE::iterator j=lst_fem_triangle.begin();
    FEM_TRIANGLE* femtriangle=(*j).second;
    MG_IDENTIFICATEUR* id=femtriangle;
    LISTE_ENTITE::iterator i=lst_entite.find(id);
    lst_entite.erase(i);
    lst_fem_triangle.erase(j);
    delete femtriangle;
 }
}


FEM_TRIANGLE* FEM_MAILLAGE::get_premier_triangle(LISTE_FEM_TRIANGLE::iterator &  it)
{
it = lst_fem_triangle.begin();
if (it == lst_fem_triangle.end())
    return NULL;
return it->second ;
}

FEM_TRIANGLE* FEM_MAILLAGE::get_suivant_triangle(LISTE_FEM_TRIANGLE::iterator &  it)
{
it++;
if(it==lst_fem_triangle.end())
    return NULL;
return it->second;
}


// ENTITE FEM_TETRA


int FEM_MAILLAGE::ajouter_fem_tetra(FEM_TETRA *femtetra)
{
gest->recherche_bonid(*femtetra);
MG_IDENTIFICATEUR *id=femtetra;
std::pair<LISTE_ENTITE::iterator,bool> p=lst_entite.insert(id);
if (!p.second)
   {
//      afficheur << ENTITEID << id->get_id() << DUPLIQUE << enderr;
      return FAIL;
   }

std::pair<const unsigned long,FEM_TETRA*> tmp(femtetra->get_id(),femtetra);
lst_fem_tetra.insert(tmp);
return OK;
}

FEM_TETRA* FEM_MAILLAGE::get_fem_tetraid(unsigned long num)
{
LISTE_FEM_TETRA::iterator i=lst_fem_tetra.find(num);
if (i==lst_fem_tetra.end())
   {
//      afficheur << INEXISTE << enderr;
      return NULL;
   }
return ((*i).second);
}

unsigned int FEM_MAILLAGE::get_nb_fem_tetra(void)
{
return lst_fem_tetra.size();
}


FEM_TETRA* FEM_MAILLAGE::get_fem_tetra(unsigned int num)
{
if (!(num<lst_fem_tetra.size()))
   {
//      afficheur << INEXISTE << enderr;
      return NULL;
   }
LISTE_FEM_TETRA::iterator i=lst_fem_tetra.begin();
for (unsigned long j=0;j<num;j++) i++;
return ((*i).second);
}


int FEM_MAILLAGE::supprimer_fem_tetraid(unsigned long num)
{
FEM_TETRA* femtetra=get_fem_tetraid(num);
if (femtetra==NULL)
   {
//      afficheur << INEXISTE2 << enderr;
      return FAIL;
   }
MG_IDENTIFICATEUR* id=femtetra;
LISTE_ENTITE::iterator i=lst_entite.find(id);
lst_entite.erase(i);
LISTE_FEM_TETRA::iterator j=lst_fem_tetra.find(num);
lst_fem_tetra.erase(j);
delete femtetra;
return OK;
}

int FEM_MAILLAGE::supprimer_fem_tetra(unsigned int num)
{
FEM_TETRA* femtetra=get_fem_tetra(num);
if (femtetra==NULL)
   {
//      afficheur << INEXISTE2 << enderr;
      return FAIL;
   }
MG_IDENTIFICATEUR* id=femtetra;
LISTE_ENTITE::iterator i=lst_entite.find(id);
lst_entite.erase(i);
LISTE_FEM_TETRA::iterator j=lst_fem_tetra.find(num);
lst_fem_tetra.erase(j);
delete femtetra;
return OK;
}

void FEM_MAILLAGE::supprimer_tout_fem_tetra(void)
{
while (get_nb_fem_tetra()!=0)
 {
    LISTE_FEM_TETRA::iterator j=lst_fem_tetra.begin();
    FEM_TETRA* femtetra=(*j).second;
    MG_IDENTIFICATEUR* id=femtetra;
    LISTE_ENTITE::iterator i=lst_entite.find(id);
    lst_entite.erase(i);
    lst_fem_tetra.erase(j);
    delete femtetra;
 }
}


FEM_TETRA* FEM_MAILLAGE::get_premier_tetra(LISTE_FEM_TETRA::iterator &  it)
{
it = lst_fem_tetra.begin();
if (it == lst_fem_tetra.end())
    return NULL;
return it->second ;
}

FEM_TETRA* FEM_MAILLAGE::get_suivant_tetra(LISTE_FEM_TETRA::iterator &  it)
{
it++;
if(it==lst_fem_tetra.end())
    return NULL;
return it->second;
}


void FEM_MAILLAGE::construire(int num)
{
if (degre==1) construire_lineraire();
if (degre==2) construire_quadratique();
if (num>0) optimise_numerotation();
if (num>1) echange_numerotation();

}


void FEM_MAILLAGE::construire_lineraire(void)
{
LISTE_MG_NOEUD::iterator it_noeud;
MG_NOEUD* mgnoeud=maillage->get_premier_noeud(it_noeud);
vector<FEM_NOEUD*> lstnoeuddirect;
unsigned int i=0;
while(mgnoeud)
        {
        mgnoeud->change_nouveau_numero(i);
        FEM_NOEUD* femnoeud=new FEM_NOEUD(mgnoeud);
        femnoeud->change_numero(i+1);
        i++;
        ajouter_fem_noeud(femnoeud);
        mgnoeud=maillage->get_suivant_noeud(it_noeud);
        lstnoeuddirect.insert(lstnoeuddirect.end(),femnoeud);
        }
int dimsansgeo=0;
if (get_mg_geometrie()==NULL)
        {
        if (maillage->get_nb_mg_tetra()!=0) dimsansgeo=3;
        else if (maillage->get_nb_mg_triangle()!=0) dimsansgeo=2;
        else if (maillage->get_nb_mg_segment()!=0) dimsansgeo=1;
        }
if ((get_mg_geometrie()!=NULL) || ( (get_mg_geometrie()==NULL) && (dimsansgeo==1)))
{
LISTE_MG_SEGMENT::iterator it_seg;
MG_SEGMENT * mgseg = maillage->get_premier_segment(it_seg);
while(mgseg)
        {
        if (mgseg->get_lien_topologie()!=NULL) if (mgseg->get_lien_topologie()->get_dimension()!=1) { mgseg=maillage->get_suivant_segment(it_seg);continue;}
                {
                FEM_NOEUD *tabnoeud[2];
                tabnoeud[0]=lstnoeuddirect[mgseg->get_noeud1()->get_nouveau_numero()];
                tabnoeud[1]=lstnoeuddirect[mgseg->get_noeud2()->get_nouveau_numero()];
                FEM_SEGMENT2* seg=new FEM_SEGMENT2(mgseg,tabnoeud);
                ajouter_fem_segment(seg);
                }
        mgseg=maillage->get_suivant_segment(it_seg);
        }
}
if ((get_mg_geometrie()!=NULL) || ( (get_mg_geometrie()==NULL) && (dimsansgeo==2)))
{
LISTE_MG_TRIANGLE::iterator it_tri;
MG_TRIANGLE* mgtri=maillage->get_premier_triangle(it_tri);
while(mgtri)
        {
        if (mgtri->get_lien_topologie()!=NULL)  if (mgtri->get_lien_topologie()->get_dimension()!=2) {mgtri=maillage->get_suivant_triangle(it_tri);continue;}
        FEM_NOEUD *tabnoeud[3];
        tabnoeud[0]=lstnoeuddirect[mgtri->get_noeud1()->get_nouveau_numero()];
        tabnoeud[1]=lstnoeuddirect[mgtri->get_noeud2()->get_nouveau_numero()];
        tabnoeud[2]=lstnoeuddirect[mgtri->get_noeud3()->get_nouveau_numero()];
        FEM_TRIANGLE3* tri=new FEM_TRIANGLE3(mgtri,tabnoeud);
        ajouter_fem_triangle(tri);
        mgtri=maillage->get_suivant_triangle(it_tri);
        }
}
if ((get_mg_geometrie()!=NULL) || ( (get_mg_geometrie()==NULL) && (dimsansgeo==3)))
{
LISTE_MG_TETRA::iterator it_tetra;
MG_TETRA* mgtetra=maillage->get_premier_tetra(it_tetra);
while (mgtetra)
        {
        FEM_NOEUD *tabnoeud[4];
        tabnoeud[0]=lstnoeuddirect[mgtetra->get_noeud1()->get_nouveau_numero()];
        tabnoeud[1]=lstnoeuddirect[mgtetra->get_noeud2()->get_nouveau_numero()];
        tabnoeud[2]=lstnoeuddirect[mgtetra->get_noeud3()->get_nouveau_numero()];
        tabnoeud[3]=lstnoeuddirect[mgtetra->get_noeud4()->get_nouveau_numero()];
        FEM_TETRA4* tet=new FEM_TETRA4(mgtetra,tabnoeud);
        ajouter_fem_tetra(tet);
        mgtetra=maillage->get_suivant_tetra(it_tetra);
        }
}
}


void FEM_MAILLAGE::construire_quadratique(void)
{
LISTE_MG_NOEUD::iterator it_nd;
MG_NOEUD* mgnoeud=maillage->get_premier_noeud(it_nd);
vector<FEM_NOEUD*> lstnoeuddirect;
unsigned int i=0;
while(mgnoeud)
        {
        mgnoeud->change_nouveau_numero(i);
        FEM_NOEUD* femnoeud=new FEM_NOEUD(mgnoeud);
        femnoeud->change_numero(i+1);
        i++  ;
        ajouter_fem_noeud(femnoeud);
        mgnoeud=maillage->get_suivant_noeud(it_nd);
        lstnoeuddirect.insert(lstnoeuddirect.end(),femnoeud);
        }
int dimsansgeo=0;
if (get_mg_geometrie()==NULL)
        {
        if (maillage->get_nb_mg_tetra()!=0) dimsansgeo=3;
        else if (maillage->get_nb_mg_triangle()!=0) dimsansgeo=2;
        else if (maillage->get_nb_mg_segment()!=0) dimsansgeo=1;
        }
int nbmgnoeud = maillage->get_nb_mg_noeud();

LISTE_MG_SEGMENT::iterator it_seg;
MG_SEGMENT* mgseg=maillage->get_premier_segment(it_seg);
i=0;
while(mgseg)
        {
        mgseg->change_nouveau_numero(i++);
        MG_NOEUD *noeud1=mgseg->get_noeud1();
        MG_NOEUD *noeud2=mgseg->get_noeud2();
        double *xyz1=noeud1->get_coord();
        double *xyz2=noeud2->get_coord();
        FEM_NOEUD* femnoeud;
        if (mgseg->get_lien_topologie()!=NULL)
        if (mgseg->get_lien_topologie()->get_dimension()==1)
                {
                MG_ARETE* arete=(MG_ARETE*)(mgseg->get_lien_topologie());
                double t1,t2;
                arete->inverser(t1,xyz1);
                arete->inverser(t2,xyz2);
                if (arete->get_courbe()->est_periodique())
                        if (t1>t2) t2=t2+arete->get_courbe()->get_periode();
                double t=0.5*(t1+t2);
                double xyz[2];
                arete->evaluer(t,xyz);
                femnoeud= new FEM_NOEUD(mgseg,xyz[0],xyz[1],xyz[2]);
                }
        if (mgseg->get_lien_topologie()!=NULL)
        if (mgseg->get_lien_topologie()->get_dimension()==2)
                {
                MG_FACE* face=(MG_FACE*)(mgseg->get_lien_topologie());
                double uv1[2],uv2[2];
                OT_DECALAGE_PARAMETRE decalage(face->get_surface()->get_periode_u(),face->get_surface()->get_periode_v());
                face->inverser(uv1,xyz1);
                face->inverser(uv2,xyz2);
                double du=decalage.calcul_decalage_parametre_u(uv1[0]);
                double dv=decalage.calcul_decalage_parametre_v(uv1[1]);
                double u1=decalage.decalage_parametre_u(uv1[0],du);
                double v1=decalage.decalage_parametre_v(uv1[1],dv);
                double u2=decalage.decalage_parametre_u(uv2[0],du);
                double v2=decalage.decalage_parametre_v(uv2[1],dv);
                double ui=0.5*(u1+u2);
                double vi=0.5*(v1+v2);
                double uii=decalage.decalage_parametre_u(ui,-du);
                double vii=decalage.decalage_parametre_v(vi,-dv);
                double uv[2],xyz[3];
                uv[0]=uii;
                uv[1]=vii;
                face->evaluer(uv,xyz);
                femnoeud= new FEM_NOEUD(mgseg,xyz[0],xyz[1],xyz[2]);
                }
        if (mgseg->get_lien_topologie()!=NULL)
                {
                if (mgseg->get_lien_topologie()->get_dimension()==3)
                        {
                        double x=0.5*(xyz1[0]+xyz2[0]);
                        double y=0.5*(xyz1[1]+xyz2[1]);
                        double z=0.5*(xyz1[2]+xyz2[2]);
                        femnoeud= new FEM_NOEUD(mgseg,x,y,z);
                        }
                }
        else
                {
                double x=0.5*(xyz1[0]+xyz2[0]);
                double y=0.5*(xyz1[1]+xyz2[1]);
                double z=0.5*(xyz1[2]+xyz2[2]);
                femnoeud= new FEM_NOEUD(mgseg,x,y,z);
                }       
        femnoeud->change_numero(i+nbmgnoeud+1);
        ajouter_fem_noeud(femnoeud);
        mgseg=maillage->get_suivant_segment(it_seg);
        lstnoeuddirect.insert(lstnoeuddirect.end(),femnoeud);
        }

if ((get_mg_geometrie()!=NULL) || ( (get_mg_geometrie()==NULL) && (dimsansgeo==1)))
{

mgseg=maillage->get_premier_segment(it_seg);
while(mgseg)
        {
        if (mgseg->get_lien_topologie()!=NULL) if (mgseg->get_lien_topologie()->get_dimension()!=1) { mgseg=maillage->get_suivant_segment(it_seg);continue;}
                {
                FEM_NOEUD *tabnoeud[3];
                tabnoeud[0]=lstnoeuddirect[mgseg->get_noeud1()->get_nouveau_numero()];
                tabnoeud[1]=lstnoeuddirect[mgseg->get_nouveau_numero()+nbmgnoeud];
                tabnoeud[2]=lstnoeuddirect[mgseg->get_noeud2()->get_nouveau_numero()];
                FEM_SEGMENT3* seg=new FEM_SEGMENT3(mgseg,tabnoeud);
                ajouter_fem_segment(seg);
                }
        mgseg=maillage->get_suivant_segment(it_seg);
        }
}
if ((get_mg_geometrie()!=NULL) || ( (get_mg_geometrie()==NULL) && (dimsansgeo==2)))
{
LISTE_MG_TRIANGLE::iterator it_tri;
MG_TRIANGLE* mgtri=maillage->get_premier_triangle(it_tri);

while(mgtri)
        {
        if (mgtri->get_lien_topologie()!=NULL)  if (mgtri->get_lien_topologie()->get_dimension()!=2) {mgtri=maillage->get_suivant_triangle(it_tri);continue;}
                {
                FEM_NOEUD *tabnoeud[6];
                tabnoeud[0]=lstnoeuddirect[mgtri->get_noeud1()->get_nouveau_numero()];
                tabnoeud[1]=lstnoeuddirect[mgtri->get_segment1()->get_nouveau_numero()+nbmgnoeud];
                tabnoeud[2]=lstnoeuddirect[mgtri->get_noeud2()->get_nouveau_numero()];
                tabnoeud[3]=lstnoeuddirect[mgtri->get_segment2()->get_nouveau_numero()+nbmgnoeud];
                tabnoeud[4]=lstnoeuddirect[mgtri->get_noeud3()->get_nouveau_numero()];
                tabnoeud[5]=lstnoeuddirect[mgtri->get_segment3()->get_nouveau_numero()+nbmgnoeud];
                FEM_TRIANGLE6* tri=new FEM_TRIANGLE6(mgtri,tabnoeud);
                ajouter_fem_triangle(tri);
                }
        mgtri=maillage->get_suivant_triangle(it_tri);
        }
}
if ((get_mg_geometrie()!=NULL) || ( (get_mg_geometrie()==NULL) && (dimsansgeo==3)))
{

LISTE_MG_TETRA::iterator it_tetra;
MG_TETRA* mgtetra=maillage->get_premier_tetra(it_tetra);
while(mgtetra)
        {
        FEM_NOEUD *tabnoeud[10];
        tabnoeud[0]=lstnoeuddirect[mgtetra->get_noeud1()->get_nouveau_numero()];
        tabnoeud[1]=lstnoeuddirect[maillage->get_mg_segment(mgtetra->get_noeud1()->get_id(),mgtetra->get_noeud2()->get_id())->get_nouveau_numero()+nbmgnoeud];
        tabnoeud[2]=lstnoeuddirect[mgtetra->get_noeud2()->get_nouveau_numero()];
        tabnoeud[3]=lstnoeuddirect[maillage->get_mg_segment(mgtetra->get_noeud2()->get_id(),mgtetra->get_noeud3()->get_id())->get_nouveau_numero()+nbmgnoeud];
        tabnoeud[4]=lstnoeuddirect[mgtetra->get_noeud3()->get_nouveau_numero()];
        tabnoeud[5]=lstnoeuddirect[maillage->get_mg_segment(mgtetra->get_noeud3()->get_id(),mgtetra->get_noeud1()->get_id())->get_nouveau_numero()+nbmgnoeud];
        tabnoeud[6]=lstnoeuddirect[maillage->get_mg_segment(mgtetra->get_noeud1()->get_id(),mgtetra->get_noeud4()->get_id())->get_nouveau_numero()+nbmgnoeud];
        tabnoeud[7]=lstnoeuddirect[maillage->get_mg_segment(mgtetra->get_noeud2()->get_id(),mgtetra->get_noeud4()->get_id())->get_nouveau_numero()+nbmgnoeud];
        tabnoeud[8]=lstnoeuddirect[maillage->get_mg_segment(mgtetra->get_noeud3()->get_id(),mgtetra->get_noeud4()->get_id())->get_nouveau_numero()+nbmgnoeud];
        tabnoeud[9]=lstnoeuddirect[mgtetra->get_noeud4()->get_nouveau_numero()];
        FEM_TETRA10* tet=new FEM_TETRA10(mgtetra,tabnoeud);
        ajouter_fem_tetra(tet);
        mgtetra=maillage->get_suivant_tetra(it_tetra);
        }
}
}


void FEM_MAILLAGE::enregistrer(std::ostream& o)
{

  if (mg_geometrie==NULL) o << "%" << get_id() << "=FEM_MAILLAGE(" << degre << ",$" << maillage->get_id() << " );" << std::endl;
  else o << "%" << get_id() << "=FEM_MAILLAGE("<< degre << ",$" << maillage->get_id() << ",$" << mg_geometrie->get_id() <<");" << std::endl;
  LISTE_ENTITE::const_iterator i;
  for (i=lst_entite.begin();i!=lst_entite.end();i++)
     (*i)->enregistrer(o);
}


void FEM_MAILLAGE::exporter_cosmos(std::string& nomfic)
{
char chaine[500];
sprintf(chaine,"%s.ses",nomfic.c_str());
std::ofstream of1(chaine,std::ios::out|std::ios::trunc);
of1.precision(16);
of1.setf(std::ios::showpoint);
of1 << "TITLE, " << chaine << " : Fichier FEM de maillage MAGiC" << std::endl;
if (degre==1) of1 << "EGROUP,1,TETRA4,0,0,0,0,0,0,0,0;" << std::endl;
if (degre==2) of1 << "EGROUP,1,TETRA10,0,0,0,0,0,0,0,0;" << std::endl;
int nb_volume=mg_geometrie->get_nb_mg_volume();
int nbNoeuds = 0;


for (int i=0;i<nb_volume;i++)
        {
        MG_VOLUME *mgvol=mg_geometrie->get_mg_volume(i);
        if (mgvol->get_num_materiau()>=0)
                {
                MT_GESTIONNAIRE mtgest(mg_geometrie->get_gest_materiau());
                MT_MATERIAU* mat=mtgest.get_materiau(mgvol->get_num_materiau());
                for(int m=0; m<mat->get_nb_propriete();m++)
                    {
                    MT_PROPRIETE* prop=mat->get_propriete(m);
                    of1<< "MPROP, "<< mgvol->get_num_materiau()+1 << ", " << prop->get_nom()  <<", " << prop->get_valeur(m) << ";" << std::endl;
                    }
               /* MT_PROPRIETE* prop=mat->get_propriete("EX");
                of1<< "MPROP, "<< mgvol->get_num_materiau()+1 << ", EX, " << prop->get_valeur(2) << ";" << std::endl;
                prop=mat->get_propriete("NUXY");
                of1<< "MPROP, "<< mgvol->get_num_materiau()+1 << ", NUXY, " << prop->get_valeur(2) << ";" << std::endl;
                prop=mat->get_propriete("GXY");
                of1<< "MPROP, "<< mgvol->get_num_materiau()+1 << ", GXY, " << prop->get_valeur(2) << ";" << std::endl;
                prop=mat->get_propriete("DENS");
                of1<< "MPROP, "<< mgvol->get_num_materiau()+1 << ", DENS, " << prop->get_valeur(2) << ";" << std::endl;
                prop=mat->get_propriete("SIGXT");
                of1<< "MPROP, "<< mgvol->get_num_materiau()+1 << ", SIGXT, " << prop->get_valeur(2) << ";" << std::endl;
                prop=mat->get_propriete("SIGYLD");
                of1<< "MPROP, "<< mgvol->get_num_materiau()+1 << ", SIGYLD, " << prop->get_valeur(2) << ";" << std::endl;
                prop=mat->get_propriete("ALPX");
                of1<< "MPROP, "<< mgvol->get_num_materiau()+1 << ", ALPX, " << prop->get_valeur(2) << ";" << std::endl;
                prop=mat->get_propriete("KX");
                of1<< "MPROP, "<< mgvol->get_num_materiau()+1 << ", KX, " << prop->get_valeur(2) << ";" << std::endl;
                prop=mat->get_propriete("C");
                of1<< "MPROP, "<< mgvol->get_num_materiau()+1 << ", C , " << prop->get_valeur(2) << ";" << std::endl;*/
                }
        }
int nb_coques = mg_geometrie->get_nb_mg_coque();
for (int i=0;i<nb_coques;i++)
        {
        MG_COQUE *mgcoque=mg_geometrie->get_mg_coque(i);
        if (mgcoque->get_num_materiau()>=0)
                {
                MT_GESTIONNAIRE mtgest(mg_geometrie->get_gest_materiau());
                MT_MATERIAU* mat=mtgest.get_materiau(mgcoque->get_num_materiau());
                for(int m=0; m<mat->get_nb_propriete();m++)
                    {
                    MT_PROPRIETE* prop=mat->get_propriete(m);
                    of1<< "MPROP, "<< mgcoque->get_num_materiau()+1 << ", " << prop->get_nom()  <<", " << prop->get_valeur(m) << ";" << std::endl;
                    }
               }
        }
int nb_poutres = mg_geometrie->get_nb_mg_poutre();
for (int i=0;i<nb_poutres;i++)
        {
        MG_POUTRE *mgpoutre=mg_geometrie->get_mg_poutre(i);
        if (mgpoutre->get_num_materiau()>=0)
                {
                MT_GESTIONNAIRE mtgest(mg_geometrie->get_gest_materiau());
                MT_MATERIAU* mat=mtgest.get_materiau(mgpoutre->get_num_materiau());

                for(int m=0; m<mat->get_nb_propriete();m++)
                    {
                    MT_PROPRIETE* prop=mat->get_propriete(m);
                    of1<< "MPROP, "<< mgpoutre->get_num_materiau()+1 << ", " << prop->get_nom()  <<", " << prop->get_valeur(m) << ";" << std::endl;
                    }
               }
        }


of1 << "ACTSET,CS,0;" << std::endl;
of1 << "ACTSET,EG,1;" << std::endl;

// �riture des noeuds
int ii=0;
double valeurunite=mg_geometrie->get_valeur_unite();
for (LISTE_FEM_NOEUD::iterator i=lst_fem_noeud.begin();i!=lst_fem_noeud.end();i++)
         {
         FEM_NOEUD* noeud=((*i).second);
         of1 << "ND, " << noeud->get_numero() << ", " << noeud->get_x()*valeurunite << ", " << noeud->get_y()*valeurunite << ", " << noeud->get_z()*valeurunite << ";" << std::endl;
         char nom[3];
         int nb=noeud->get_lien_topologie()->get_nb_ccf();
         for (int k=0;k<nb;k++)
                {
                noeud->get_lien_topologie()->get_type_ccf(k,nom);
                if (strcmp(nom,"Da")==0) of1 << "DND , " << noeud->get_numero()  << ", AL, " << noeud->get_lien_topologie()->get_valeur_ccf(k) << ", " << noeud->get_numero() << ", 1;" << std::endl;
                if (strcmp(nom,"Dx")==0) of1 << "DND , " << noeud->get_numero()  << ", UX, " << noeud->get_lien_topologie()->get_valeur_ccf(k) << ", " << noeud->get_numero() << ", 1;" << std::endl;
                if (strcmp(nom,"Dy")==0) of1 << "DND , " << noeud->get_numero()  << ", UY, " << noeud->get_lien_topologie()->get_valeur_ccf(k) << ", " << noeud->get_numero() << ", 1;" << std::endl;
                if (strcmp(nom,"Dz")==0) of1 << "DND , " << noeud->get_numero()  << ", UZ, " << noeud->get_lien_topologie()->get_valeur_ccf(k) << ", " << noeud->get_numero() << ", 1;" << std::endl;
                if (strcmp(nom,"Fx")==0) of1 << "FND , " << noeud->get_numero()  << ", FX, " << noeud->get_lien_topologie()->get_valeur_ccf(k) << ", " << noeud->get_numero() << ", 1;" << std::endl;
                if (strcmp(nom,"Fy")==0) of1 << "FND , " << noeud->get_numero()  << ", FY, " << noeud->get_lien_topologie()->get_valeur_ccf(k) << ", " << noeud->get_numero() << ", 1;" << std::endl;
                if (strcmp(nom,"Fz")==0) of1 << "FND , " << noeud->get_numero()  << ", FZ, " << noeud->get_lien_topologie()->get_valeur_ccf(k) << ", " << noeud->get_numero() << ", 1;" << std::endl;
                }
         ii++;
         }
nbNoeuds = ii;
ii=0;

// �criture des t�tra�dres
int nummat=-10;   // num�o du mat�iau
for (LISTE_FEM_TETRA::iterator i=lst_fem_tetra.begin();i!=lst_fem_tetra.end();i++)
        {
        FEM_TETRA* tetra=((*i).second);
        MG_VOLUME* vol=(MG_VOLUME*)tetra->get_lien_topologie();
        int nummatac=vol->get_num_materiau();
        if ((nummatac!=nummat) && (nummatac!=-1))
                {
                nummat=nummatac;
                of1 << "ACTSET, MP," << nummat+1 << ";" << std::endl;
                }
        if (degre==1) of1 << "EL, " << ii+1 << ", VL, 1, 4, " << tetra->get_fem_noeud(0)->get_numero() << ", " << tetra->get_fem_noeud(1)->get_numero() << ", " << tetra->get_fem_noeud(2)->get_numero() << ", " << tetra->get_fem_noeud(3)->get_numero() << ", "  ;
        if (degre==2) of1 << "EL, " << ii+1 << ", VL, 1,10, " << tetra->get_fem_noeud(0)->get_numero() << ", " << tetra->get_fem_noeud(2)->get_numero() << ", " << tetra->get_fem_noeud(4)->get_numero() << ", " << tetra->get_fem_noeud(9)->get_numero() << ", " << tetra->get_fem_noeud(1)->get_numero() << ", " << tetra->get_fem_noeud(3)->get_numero() << ", " << tetra->get_fem_noeud(5)->get_numero() << ", " << tetra->get_fem_noeud(6)->get_numero() << ", " << tetra->get_fem_noeud(7)->get_numero() << ", " << tetra->get_fem_noeud(8)->get_numero() << ", "   ;
        ii++;
        MG_TETRA* tet=(MG_TETRA*)tetra->get_mg_element_maillage();
        if (tet->get_triangle1()->get_lien_topologie()->get_dimension()!=3) of1 << "1, "; else of1 << "0, ";
        if (tet->get_triangle2()->get_lien_topologie()->get_dimension()!=3) of1 << "1, "; else of1 << "0, ";
        if (tet->get_triangle3()->get_lien_topologie()->get_dimension()!=3) of1 << "1, "; else of1 << "0, ";
        if (tet->get_triangle4()->get_lien_topologie()->get_dimension()!=3) of1 << "1;" << std::endl ; else of1 << "0;" << std::endl;
        char nom[3];
        int nb=tet->get_triangle1()->get_lien_topologie()->get_nb_ccf();
        for (int k=0;k<nb;k++)
                {
                tet->get_triangle1()->get_lien_topologie()->get_type_ccf(k,nom);
                if (strcmp(nom,"Px")==0) of1 << "PEL , " << ii <<  ", " << tet->get_triangle1()->get_lien_topologie()->get_valeur_ccf(k) << ", 1, " << ii << ", 1, 1;" << std::endl;
                if (strcmp(nom,"Py")==0) of1 << "PEL , " << ii <<  ", " << tet->get_triangle1()->get_lien_topologie()->get_valeur_ccf(k) << ", 1, " << ii << ", 1, 2;" << std::endl;
                if (strcmp(nom,"Pz")==0) of1 << "PEL , " << ii <<  ", " << tet->get_triangle1()->get_lien_topologie()->get_valeur_ccf(k) << ", 1, " << ii << ", 1, 3;" << std::endl;
                if (strcmp(nom,"Pn")==0) of1 << "PEL , " << ii <<  ", " << tet->get_triangle1()->get_lien_topologie()->get_valeur_ccf(k) << ", 1, " << ii << ", 1, 4;" << std::endl;
                }
        nb=tet->get_triangle2()->get_lien_topologie()->get_nb_ccf();
        for (int k=0;k<nb;k++)
                {
                tet->get_triangle2()->get_lien_topologie()->get_type_ccf(k,nom);
                if (strcmp(nom,"Px")==0) of1 << "PEL , " << ii <<  ", " << tet->get_triangle2()->get_lien_topologie()->get_valeur_ccf(k) << ", 2, " << ii << ", 1, 1;" << std::endl;
                if (strcmp(nom,"Py")==0) of1 << "PEL , " << ii <<  ", " << tet->get_triangle2()->get_lien_topologie()->get_valeur_ccf(k) << ", 2, " << ii << ", 1, 2;" << std::endl;
                if (strcmp(nom,"Pz")==0) of1 << "PEL , " << ii <<  ", " << tet->get_triangle2()->get_lien_topologie()->get_valeur_ccf(k) << ", 2, " << ii << ", 1, 3;" << std::endl;
                if (strcmp(nom,"Pn")==0) of1 << "PEL , " << ii <<  ", " << tet->get_triangle2()->get_lien_topologie()->get_valeur_ccf(k) << ", 2, " << ii << ", 1, 4;" << std::endl;
                }
        nb=tet->get_triangle3()->get_lien_topologie()->get_nb_ccf();
        for (int k=0;k<nb;k++)
                {
                tet->get_triangle3()->get_lien_topologie()->get_type_ccf(k,nom);
                if (strcmp(nom,"Px")==0) of1 << "PEL , " << ii <<  ", " << tet->get_triangle3()->get_lien_topologie()->get_valeur_ccf(k) << ", 3, " << ii << ", 1, 1;" << std::endl;
                if (strcmp(nom,"Py")==0) of1 << "PEL , " << ii <<  ", " << tet->get_triangle3()->get_lien_topologie()->get_valeur_ccf(k) << ", 3, " << ii << ", 1, 2;" << std::endl;
                if (strcmp(nom,"Pz")==0) of1 << "PEL , " << ii <<  ", " << tet->get_triangle3()->get_lien_topologie()->get_valeur_ccf(k) << ", 3, " << ii << ", 1, 3;" << std::endl;
                if (strcmp(nom,"Pn")==0) of1 << "PEL , " << ii <<  ", " << tet->get_triangle3()->get_lien_topologie()->get_valeur_ccf(k) << ", 3, " << ii << ", 1, 4;" << std::endl;
                }
        nb=tet->get_triangle4()->get_lien_topologie()->get_nb_ccf();
        for (int k=0;k<nb;k++)
                {
                tet->get_triangle4()->get_lien_topologie()->get_type_ccf(k,nom);
                if (strcmp(nom,"Px")==0) of1 << "PEL , " << ii <<  ", " << tet->get_triangle4()->get_lien_topologie()->get_valeur_ccf(k) << ", 4, " << ii << ", 1, 1;" << std::endl;
                if (strcmp(nom,"Py")==0) of1 << "PEL , " << ii <<  ", " << tet->get_triangle4()->get_lien_topologie()->get_valeur_ccf(k) << ", 4, " << ii << ", 1, 2;" << std::endl;
                if (strcmp(nom,"Pz")==0) of1 << "PEL , " << ii <<  ", " << tet->get_triangle4()->get_lien_topologie()->get_valeur_ccf(k) << ", 4, " << ii << ", 1, 3;" << std::endl;
                if (strcmp(nom,"Pn")==0) of1 << "PEL , " << ii <<  ", " << tet->get_triangle4()->get_lien_topologie()->get_valeur_ccf(k) << ", 4, " << ii << ", 1, 4;" << std::endl;
                }
        }


        //*********************************************************
        // Les t�tras sont finis, maintenant le maillage des Coques
        //*********************************************************
//nb_coques=mg_geometrie->get_nb_mg_coque();  //variable d��d�lar� plus haut, elle devrait avoir conserv�sa valeur
if (nb_coques>0)
    {
    if (degre==1) of1 << "EGROUP,2,SHELL3,0,0,0,0,0,0,0,0;" << std::endl;  // Formulation coques minces... SHELL3T pout �aisses
    if (degre==2) of1 << "EGROUP,2,SHELL6,0,0,0,0,0,0,0,0;" << std::endl;
    }

if (nb_poutres>0 || nb_coques>0 )
{
    of1 << "EGROUP,3,BEAM3D,0,0,0,0,0,0,0,0;"<< std::endl;   //Egroup des mini-poutres
    of1 << "RCONST,3,3,1,3," << "1, 1 , 1; " << std::endl;
    of1 << "EGROUP,4,RBAR;"<< std::endl;  // EGroup des RBAR
    of1 << "EGROUP,5,BEAM3D,0,0,0,0,0,0,0,0;"<< std::endl;     //Egroup des mini-poutres infiniment rigide
    of1 << "RCONST,5,5,1,3," << "0.01, 0.01 , 0.01; " << std::endl;
    of1 << "EGROUP,6,BEAM3D,0,0,0,0,0,0,0,0;"<< std::endl;     //Egroup des mini-poutres avec inertie dans juste une direction (je sais, il n'est pas obligatoire, mais �a aide au d�buggage d'avoir un #5...)
    of1 << "RCONST,6,6,1,3," << "1e-16, 0.01 , 1e-16; " << std::endl;
}

int RCset =7;

LISTE_MG_COQUE::iterator ico;
MG_COQUE *mgcoque=mg_geometrie->get_premier_coque(ico);
while(mgcoque)
        {
        int nummatac=mgcoque->get_num_materiau();
        if ((nummatac!=nummat)&&(nummatac!=-1)) // si c'est = 0 alors on place le dernier mat�iau choisi.
                {
                nummat=nummatac;
                of1 << "ACTSET, MP," << nummat+1 << ";" << std::endl;
                }

        for (int numcoquille=0;numcoquille<mgcoque->get_nb_mg_coquille();numcoquille++)
            {
            MG_COQUILLE* mgcoquille = mgcoque->get_mg_coquille(numcoquille);
            for (int cf = 0; cf<mgcoquille->get_nb_mg_coface();cf++)
                {
                MG_COFACE* mgcoface = mgcoquille->get_mg_coface(cf);
                MG_FACE* mgface = mgcoface->get_face();

                double epaisseur=-1;
                for (int k=0;k<mgface->get_nb_ccf();k++)  // aller chercher l'�aisseur de la coque
                    {
                    char nomCCF[3];
                    mgface->get_type_ccf(k,nomCCF);
                    if (strcmp(nomCCF,"EP")==0)
                        epaisseur = mgface-> get_valeur_ccf(k);
                    }
                of1 << "ACTSET,EG,2;"<< std::endl;
                of1 << "RCONST,2,"<< ++RCset << ",1,1," << epaisseur << ";" << std::endl;       // le real constant, s'active automatiquement

                TPL_LISTE_ENTITE<FEM_ELEMENT_MAILLAGE*> * lien_maillage = mgface->get_lien_fem_maillage();
                int nbelem=lien_maillage->get_nb();
                for (int i=0;i<nbelem;i++)
                  {
                    FEM_TRIANGLE* femtriangle=(FEM_TRIANGLE*)lien_maillage->get(i);
                    if (!contient(femtriangle)) continue;
                    if (degre==1) of1 << "EL, " << ++ii <<  ", SF, 0, 3, " <<femtriangle->get_fem_noeud(0)->get_numero()<< ","<< femtriangle->get_fem_noeud(1)->get_numero()<<","<<femtriangle->get_fem_noeud(2)->get_numero()<< ";" <<  std::endl;
                    if (degre==2) of1 << "EL, " << ++ii <<  ", SF, 0, 6, " <<femtriangle->get_fem_noeud(0)->get_numero()<< ","<< femtriangle->get_fem_noeud(2)->get_numero()<<","<<femtriangle->get_fem_noeud(4)->get_numero()<< ","<<femtriangle->get_fem_noeud(1)->get_numero()<< ","<<femtriangle->get_fem_noeud(3)->get_numero()<<","<<femtriangle->get_fem_noeud(5)->get_numero()<<";"<< std::endl;
                    // dans les noeuds des sommets sont, 1, 2 et 3.  les noeuds sur les segments sont 4 5 et 6

                    long nbc;
                    MG_TRIANGLE* mgtri=(MG_TRIANGLE*)femtriangle->get_mg_element_maillage();
                    char nom[3];

                    nbc = mgtri->get_segment1()->get_lien_topologie()->get_nb_ccf();
                    for (int z =0 ; z< nbc; z++)    // le premier segment
                        { // il y a au moins une CCF "EP" qu'il faut ignorer
                        mgtri->get_segment1()->get_lien_topologie()->get_type_ccf(z,nom); //  ->get_triangle3()->get_lien_topologie()->get_type_ccf(k,nom);
                        if (strcmp(nom,"EP")==0) continue;
                        if (strcmp(nom,"Px")==0) of1 << "PEL, " << ii <<  ", " << mgtri->get_segment1()->get_lien_topologie()->get_valeur_ccf(z) << ", 1, " << ii << ", 1, 1;  C* CCF de Coque... " << std::endl;
                        if (strcmp(nom,"Py")==0) of1 << "PEL, " << ii <<  ", " << mgtri->get_segment1()->get_lien_topologie()->get_valeur_ccf(z)<< ", 1, " << ii << ", 1, 2;  C* CCF de Coque... " << std::endl;
                        if (strcmp(nom,"Pz")==0) of1 << "PEL, " << ii <<  ", " << mgtri->get_segment1()->get_lien_topologie()->get_valeur_ccf(z) << ", 1, " << ii << ", 1, 3;  C* CCF de Coque... "<< std::endl;
                        if (strcmp(nom,"Pn")==0) of1 << "PEL, " << ii <<  ", " << mgtri->get_segment1()->get_lien_topologie()->get_valeur_ccf(z) << ", 1, " << ii << ", 1, 4;  C* CCF de Coque... " << std::endl;
                        }
                    nbc = mgtri->get_segment2()->get_lien_topologie()->get_nb_ccf();
                    for (int z =0 ; z< nbc; z++)
                        {
                        mgtri->get_segment2()->get_lien_topologie()->get_type_ccf(z,nom);
                        if (strcmp(nom,"EP")==0) continue;
                        if (strcmp(nom,"Px")==0) of1 << "PEL, " << ii <<  ", " << mgtri->get_segment2()->get_lien_topologie()->get_valeur_ccf(z) << ", 2, " << ii << ", 1, 1;  C* CCF de Coque... " << std::endl;
                        if (strcmp(nom,"Py")==0) of1 << "PEL, " << ii <<  ", " << mgtri->get_segment2()->get_lien_topologie()->get_valeur_ccf(z)  << ", 2, " << ii << ", 1, 2;  C* CCF de Coque... " << std::endl;
                        if (strcmp(nom,"Pz")==0) of1 << "PEL, " << ii <<  ", " << mgtri->get_segment2()->get_lien_topologie()->get_valeur_ccf(z)  << ", 2, " << ii << ", 1, 3;  C* CCF de Coque... "<< std::endl;
                        if (strcmp(nom,"Pn")==0) of1 << "PEL, " << ii <<  ", " << mgtri->get_segment2()->get_lien_topologie()->get_valeur_ccf(z)  << ", 2, " << ii << ", 1, 4;  C* CCF de Coque... " << std::endl;
                        }
                    nbc = mgtri->get_segment3()->get_lien_topologie()->get_nb_ccf();
                    for (int z =0 ; z< nbc; z++)
                        {
                        mgtri->get_segment3()->get_lien_topologie()->get_type_ccf(z,nom);
                        if (strcmp(nom,"EP")==0) continue;
                        if (strcmp(nom,"Px")==0) of1 << "PEL, " << ii <<  ", " << mgtri->get_segment3()->get_lien_topologie()->get_valeur_ccf(z) << ", 3, " << ii << ", 1, 1;  C* CCF de Coque... " << std::endl;
                        if (strcmp(nom,"Py")==0) of1 << "PEL, " << ii <<  ", " << mgtri->get_segment3()->get_lien_topologie()->get_valeur_ccf(z) << ", 3, " << ii << ", 1, 2;  C* CCF de Coque... " << std::endl;
                        if (strcmp(nom,"Pz")==0) of1 << "PEL, " << ii <<  ", " << mgtri->get_segment3()->get_lien_topologie()->get_valeur_ccf(z) << ", 3, " << ii << ", 1, 3;  C* CCF de Coque... "<< std::endl;
                        if (strcmp(nom,"Pn")==0) of1 << "PEL, " << ii <<  ", " << mgtri->get_segment3()->get_lien_topologie()->get_valeur_ccf(z) << ", 3, " << ii << ", 1, 4;  C* CCF de Coque... " << std::endl;
                        }


                    nbc = mgtri->get_lien_topologie()->get_nb_ccf();   // et pour la face elle-m�e
                    for (int z =0 ; z< nbc; z++)
                        {
                        mgtri->get_lien_topologie()->get_type_ccf(z,nom);
                        if (strcmp(nom,"EP")==0) continue;
                        if (strcmp(nom,"Px")==0) of1 << "PEL, " << ii <<  ", " << mgtri->get_lien_topologie()->get_valeur_ccf(z) << ", 4, " << ii << ", 1, 1;  C* CCF de Coque... " << std::endl;
                        if (strcmp(nom,"Py")==0) of1 << "PEL, " << ii <<  ", " << mgtri->get_lien_topologie()->get_valeur_ccf(z) << ", 4, " << ii << ", 1, 2;  C* CCF de Coque... " << std::endl;
                        if (strcmp(nom,"Pz")==0) of1 << "PEL, " << ii <<  ", " << mgtri->get_lien_topologie()->get_valeur_ccf(z) << ", 4, " << ii << ", 1, 3;  C* CCF de Coque... "<< std::endl;
                        if (strcmp(nom,"Pn")==0) of1 << "PEL, " << ii <<  ", " << mgtri->get_lien_topologie()->get_valeur_ccf(z) << ", 4, " << ii << ", 1, 4;  C* CCF de Coque... " << std::endl;
                        }

                    }

                }  // next face
            } // next coquille
        mgcoque = mg_geometrie->get_suivant_coque(ico);
        }// next coque


// **************************
// maintenant les poutres
// **************************


//for (int i = 0 ; i<nb_poutres;i++)
LISTE_MG_POUTRE::iterator it_po;
MG_POUTRE* mgpoutre = mg_geometrie->get_premier_poutre(it_po);
while(mgpoutre)
    {
    //MG_POUTRE* mgpoutre = mg_geometrie->get_mg_poutre(i);
    int nummatac=mgpoutre->get_num_materiau();
    if ((nummatac!=nummat)&&(nummatac!=-1))
        {
        nummat=nummatac;
        of1 << "ACTSET, MP," << nummat+1 << ";" << std::endl;
        }

    for (int j = 0;j<mgpoutre->get_nb_mg_boucle();j++)
        {
        MG_BOUCLE* mgboucle = mgpoutre->get_mg_boucle(j);

        for (int k = 0 ; k<mgboucle->get_nb_mg_coarete();k++)
            {
            MG_COARETE* mgcoarete = mgboucle->get_mg_coarete(k);
            MG_ARETE* mgarete = mgcoarete->get_arete();

            double IXX = 0;
            double IYY = 0;
            double Aire = 0;
            double x3 = 0;
            double y3 = 0;
            double z3 = 0;

            for (int k=0;k<mgarete->get_nb_ccf();k++)
                {
                char nomCCF[3];
                mgarete->get_type_ccf(k,nomCCF);
                if (strcmp(nomCCF,"As")==0)
                    Aire = mgarete->get_valeur_ccf(k);
                if (strcmp(nomCCF,"IX")==0)
                    IXX = mgarete->get_valeur_ccf(k);
                if (strcmp(nomCCF,"IY")==0)
                    IYY = mgarete->get_valeur_ccf(k);
                if (strcmp(nomCCF,"NX")==0)
                    x3 = mgarete->get_valeur_ccf(k);
                if (strcmp(nomCCF, "NY")==0)
                    y3 = mgarete->get_valeur_ccf(k);
                if (strcmp(nomCCF, "NZ")==0)
                   z3 = mgarete->get_valeur_ccf(k);
                }

            if( Aire==-7) // -7 est un flag pour diff�rencier les mini-poutres...
                {     // mini-poutre
                of1 << "Actset,EG,5;"<<std::endl;
                of1 << "Actset,RC,5;"<<std::endl;
                of1 << "C* MINI - POUTRE" << std:: endl;

              for(int t = 0; t< mgarete->get_lien_fem_maillage()->get_nb();t++)
                    {
                    FEM_SEGMENT* femseg = (FEM_SEGMENT*)mgarete->get_lien_fem_maillage()->get(t);
                    of1 << "EL, " << ++ii <<  ", CR, 0, 2, " <<femseg->get_fem_noeud(0)->get_numero()<< ","<< femseg->get_fem_noeud(1)->get_numero()<<";" <<  std::endl;
                    }   // next segment
                }
            else
                {     // poutre normale
                of1 << "Actset,EG,3;"<<std::endl;
                of1 << "ND,"<< ++nbNoeuds/* NDMAX+1 */ << "," << x3*valeurunite << "," << y3*valeurunite << "," << z3*valeurunite << ";" << std::endl;          // cr�ation du trois�me point (qui est un noeud) ND,1,5,6,7
                of1 << "RCONST,3,"<< ++RCset << ",1,3," << Aire << "," << IXX << "," << IYY << ";" << std::endl;            //RCONST,3,66,1,3,99,98,97  ; EG associ�e RC , start location, nombre de constantes, aire, inertie XX, inertie YY

                for(int t = 0; t< mgarete->get_lien_fem_maillage()->get_nb();t++)
                    {
                    FEM_SEGMENT* femseg = (FEM_SEGMENT*)mgarete->get_lien_fem_maillage()->get(t);
                    if (!contient(femseg)) continue;
                    if (degre==1) of1 << "EL, " << ++ii <<  ", CR, 0, 3, " <<femseg->get_fem_noeud(0)->get_numero()<< ","<< femseg->get_fem_noeud(1)->get_numero()<<","<< nbNoeuds << ";" <<  std::endl;

                    for (int k=0;k<mgarete->get_nb_ccf();k++)
                        {
                        char nom[3];
                        mgarete->get_type_ccf(k,nom);
                        if (strcmp(nom,"Px")==0) of1 << "PEL, " << ii <<  ", " << femseg->get_lien_topologie()->get_valeur_ccf(k) << ", 4, " << ii << ", 1, 1;  C* CCF de Coque... " << std::endl;
                        if (strcmp(nom,"Py")==0) of1 << "PEL, " << ii <<  ", " << femseg->get_lien_topologie()->get_valeur_ccf(k) << ", 4, " << ii << ", 1, 2;  C* CCF de Coque... " << std::endl;
                        if (strcmp(nom,"Pz")==0) of1 << "PEL, " << ii <<  ", " << femseg->get_lien_topologie()->get_valeur_ccf(k) << ", 4, " << ii << ", 1, 3;  C* CCF de Coque... "<<  std::endl;
                        if (strcmp(nom,"Pn")==0) of1 << "PEL, " << ii <<  ", " << femseg->get_lien_topologie()->get_valeur_ccf(k) << ", 4, " << ii << ", 1, 4;  C* CCF de Coque... " << std::endl;  // l�� va en arracher!!!
                        }
                    }   // next segment
                }


            } // next arete
        } // next boucle
    mgpoutre=mg_geometrie->get_suivant_poutre(it_po);
    } // next poutre

// Cas sp�cial, une ar�te normale est flagg�e comme une mini-poutre (dans coupecot�)
LISTE_MG_ARETE::iterator ita;
MG_ARETE* mgarete = mg_geometrie->get_premier_arete(ita);

while(mgarete)
    {
    double Aire = 0;

    for (int k=0;k<mgarete->get_nb_ccf();k++)
        {
        char nomCCF[3];
        mgarete->get_type_ccf(k,nomCCF);
        if (strcmp(nomCCF,"As")==0)
            Aire = mgarete->get_valeur_ccf(k);
        }

    if( Aire==-7) // -7 est un flag pour diff�rencier les mini-poutres...  non orient�e
        {     // mini-poutre
        of1 << "Actset,EG,5;"<<std::endl;
        of1 << "Actset,RC,5;"<<std::endl;
        of1 << "C* MINI - POUTRE" << std:: endl;

        for(int t = 0; t< mgarete->get_lien_fem_maillage()->get_nb();t++)
            {
            FEM_SEGMENT* femseg = (FEM_SEGMENT*)mgarete->get_lien_fem_maillage()->get(t);
            if (!contient(femseg)) continue;
            of1 << "EL, " << ++ii <<  ", CR, 0, 2, " <<femseg->get_fem_noeud(0)->get_numero()<< ","<< femseg->get_fem_noeud(1)->get_numero()<<";" <<  std::endl;
            }
        }
        mgarete=mg_geometrie->get_suivant_arete(ita);
    }


    //************************************************************
    // Pour mettre des mini-poutres sur les faces concern�es *****
    //*************************************************************
// identifications des faces du qui sont des faces internes.

LISTE_MG_FACE::iterator itf;
MG_FACE* mgface = mg_geometrie->get_premier_face(itf);


while(mgface)
{
    for (int k=0;k<mgface->get_nb_ccf();k++)
    {
        char nomCCF[3];
        mgface->get_type_ccf(k,nomCCF);

        if (strcmp(nomCCF,"FI")==0)
        {
            TPL_LISTE_ENTITE<FEM_ELEMENT_MAILLAGE*> * lien_fem_maillage = mgface->get_lien_fem_maillage();
            FEM_ELEMENT_MAILLAGE* element;

            //si on a une faceinterne de poutre alors on a une CCF "Po" = 0
            // sinon la CCF Po est �gale � 1 et on a une coque
           int Poutre = -1;

           for (int z=0;z<mgface->get_nb_ccf();z++)
           {
                char nomCCF2[3];
                mgface->get_type_ccf(z,nomCCF);
                if(strcmp(nomCCF,"Po")==0)
                {
                    Poutre = mgface->get_valeur_ccf(z);
                }
           }


           if (Poutre ==0)
           { // mini-poutres pour des poutres
            of1 << "ACTSET,EG,5;  C* Face avec mini-poutres de POUTRE"<<std::endl;
            of1 << "ACTSET,RC,5;"<<std::endl;
                for (int i = 0 ; i< lien_fem_maillage->get_nb();i++)
                {
                    FEM_TRIANGLE* femtriangle=(FEM_TRIANGLE*)lien_fem_maillage->get(i) ;

                    FEM_NOEUD* nd0 =femtriangle->get_fem_noeud(0);
                    FEM_NOEUD* nd1 =femtriangle->get_fem_noeud(1);
                    FEM_NOEUD* nd2= femtriangle->get_fem_noeud(2);

                    // Pour des mini-poutres de rigidit� infinie    //EL,1,CR,0,3,1,2,3,0,0,0,0,0,0  num ��ent, CR, flag asso, nb noeuds, #1, #2, #3  ;
                    //si on veut orienter les propri�t�s,
                    of1 << "EL, " << ++ii <<  ", CR, 0, 2, " <<nd0->get_numero()<< "," << nd1->get_numero()<< ";" <<  std::endl;
                    of1 << "EL, " << ++ii <<  ", CR, 0, 2, " <<nd1->get_numero()<< "," << nd2->get_numero()<< ";" <<  std::endl;
                    of1 << "EL, " << ++ii <<  ", CR, 0, 2, " <<nd0->get_numero()<< "," << nd2->get_numero()<< ";" <<  std::endl;
                }
           }
           else if(Poutre==1)
           {  // mini-poutres pour des coques
            of1 << "ACTSET,EG,6;  C* Face avec mini-poutres de COQUE"<<std::endl;
            of1 << "ACTSET,RC,6;"<<std::endl;
                for (int i = 0 ; i< lien_fem_maillage->get_nb();i++)
                {
                    FEM_TRIANGLE* femtriangle=(FEM_TRIANGLE*)lien_fem_maillage->get(i) ;

                    FEM_NOEUD* nd0 =femtriangle->get_fem_noeud(0);
                    FEM_NOEUD* nd1 =femtriangle->get_fem_noeud(1);
                    FEM_NOEUD* nd2= femtriangle->get_fem_noeud(2);

                    of1 << "EL, " << ++ii <<  ", CR, 0, 3, " <<nd0->get_numero()<< "," << nd1->get_numero()<< "," << nd2->get_numero()<< ";" <<  std::endl;
                    of1 << "EL, " << ++ii <<  ", CR, 0, 3, " <<nd1->get_numero()<< "," << nd2->get_numero()<< "," << nd0->get_numero()<< ";" <<  std::endl;
                    of1 << "EL, " << ++ii <<  ", CR, 0, 3, " <<nd0->get_numero()<< "," << nd2->get_numero()<< "," << nd1->get_numero()<< ";" <<  std::endl;
                }
           }
        } // end if FI
    }  //   next CCF
        mgface= mg_geometrie->get_suivant_face(itf);
}   // end while face


//of1<< "NMERGE;" <<std::endl;
//of1<< "NCOMPRESS;" <<std::endl;

//of1<< "SETEPLOT,1,1,0,0;"<<std::endl;
//of1<< "ACTECLR,1,RC,1;"<<std::endl;

of1<< "EPLOT;"<<std::endl;
of1<< "PSCALE;"<<std::endl;
//lst_FaceInternes.vide();

}


double FEM_MAILLAGE::calcul_coef(void)
{
if (!deforme) return 0.;
LISTE_FEM_NOEUD::iterator  it;
double max=0.,dmax=0.;
for (FEM_NOEUD* noeud=get_premier_noeud(it);noeud;noeud=get_suivant_noeud(it)) 
        {
        double x=noeud->get_x();
        double y=noeud->get_y();
        double z=noeud->get_z();
        double dx=noeud->get_dx();
        double dy=noeud->get_dy();
        double dz=noeud->get_dz();
        if (fabs(x)>max) max=fabs(x);
        if (fabs(y)>max) max=fabs(y);
        if (fabs(z)>max) max=fabs(z);
        if (fabs(dx)>dmax) dmax=fabs(dx);
        if (fabs(dy)>dmax) dmax=fabs(dy);
        if (fabs(dz)>dmax) dmax=fabs(dz);
        }
if (dmax<1e-10) return 0;
double coef=0.1*max/dmax;
return coef;
}


void FEM_MAILLAGE::calcul_deforme(FEM_SOLUTION* solx,int numx,FEM_SOLUTION* soly,int numy,FEM_SOLUTION* solz,int numz)
{
LISTE_FEM_NOEUD::iterator  it;
solx->active_solution(numx);
for (FEM_NOEUD* noeud=get_premier_noeud(it);noeud;noeud=get_suivant_noeud(it))

        noeud->change_dx(noeud->get_solution()/get_mg_geometrie()->get_valeur_unite());

soly->active_solution(numy);
for (FEM_NOEUD* noeud=get_premier_noeud(it);noeud;noeud=get_suivant_noeud(it))

        noeud->change_dy(noeud->get_solution()/get_mg_geometrie()->get_valeur_unite());

solz->active_solution(numz);
for (FEM_NOEUD* noeud=get_premier_noeud(it);noeud;noeud=get_suivant_noeud(it))

        noeud->change_dz(noeud->get_solution()/get_mg_geometrie()->get_valeur_unite());

deforme=1;
}


int FEM_MAILLAGE::existe_deforme(void)
{
return deforme;
}

void FEM_MAILLAGE::active_deforme(void)
{
deforme=1;
}


void FEM_MAILLAGE::optimise_numerotation(void)
{
FEM_NOEUD* noeud=get_fem_noeud(0);
FEM_GRAPH_NOEUD *graph;
graph=new FEM_GRAPH_NOEUD(noeud,this);
int e=graph->get_excentricite();
TPL_MAP_ENTITE<FEM_NOEUD*> dernierniveau=graph->get_dernier_niveau();
int nbdernierniveau=dernierniveau.get_nb();
for (int i=0;i<nbdernierniveau;i++)
        {
        FEM_GRAPH_NOEUD* graphtmp=new FEM_GRAPH_NOEUD(dernierniveau.get(i),this);
        int etmp=graphtmp->get_excentricite();
        if (etmp>e)
                {
                delete graph;
                graph=graphtmp;
                dernierniveau=graph->get_dernier_niveau();
                nbdernierniveau=dernierniveau.get_nb();
                i=-1;
                e=etmp;
                }
        else delete graphtmp;
        }
int numero=get_nb_fem_noeud();
for (int i=0;i<e+1;i++)
        {
        TPL_MAP_ENTITE<FEM_NOEUD*> niveau=graph->get_niveau(i);
        int nbnoeud=niveau.get_nb();
        for (int j=0;j<nbnoeud;j++)
                {
                FEM_NOEUD* noeud=niveau.get(j);
                noeud->change_numero_opt(numero);
                numero--; 
                }
        }


delete graph;
}

void FEM_MAILLAGE::echange_numerotation(void)
{
int nbnoeud=get_nb_fem_noeud();
for (int i=0;i<nbnoeud;i++)
        {
        FEM_NOEUD* noeud=get_fem_noeud(i);
        int num=noeud->get_numero();
        noeud->change_numero(noeud->get_numero_opt());
        noeud->change_numero_opt(num);
        }

}


std::ostream& operator << (std::ostream& o,FEM_MAILLAGE& maillage)
{
maillage.enregistrer(o);
return o;
}


Revision:	253
Committed:	Tue Jul 13 19:40:46 2010 UTC (14 years, 10 months ago) by francois
File size:	59670 byte(s)
Log Message:	changement de hiearchie et utilisation de ccmake + mise a jour