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 "mc_gestionnaire.h"
#include "mc_propriete.h"
#include "mc_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_element3();
    supprimer_tout_fem_element2();
    supprimer_tout_fem_element1();
    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_ELEMENT1


int FEM_MAILLAGE::ajouter_fem_element1(FEM_ELEMENT1 *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_ELEMENT1*> tmp(femsegment->get_id(),femsegment);
    lst_fem_element1.insert(tmp);
    return OK;
}

FEM_ELEMENT1* FEM_MAILLAGE::get_fem_element1id(unsigned long num)
{
    LISTE_FEM_ELEMENT1::iterator i=lst_fem_element1.find(num);
    if (i==lst_fem_element1.end())
    {
//      afficheur << INEXISTE << enderr;
        return NULL;
    }
    return ((*i).second);
}

unsigned int FEM_MAILLAGE::get_nb_fem_element1(void)
{
    return lst_fem_element1.size();
}


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


int FEM_MAILLAGE::supprimer_fem_element1id(unsigned long num)
{
    FEM_ELEMENT1* femsegment=get_fem_element1id(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_ELEMENT1::iterator j=lst_fem_element1.find(num);
    lst_fem_element1.erase(j);
    delete femsegment;
    return OK;
}

int FEM_MAILLAGE::supprimer_fem_element1(unsigned int num)
{
    FEM_ELEMENT1* femsegment=get_fem_element1(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_ELEMENT1::iterator j=lst_fem_element1.find(num);
    lst_fem_element1.erase(j);
    delete femsegment;
    return OK;
}

void FEM_MAILLAGE::supprimer_tout_fem_element1(void)
{
    while (get_nb_fem_element1()!=0)
    {
        LISTE_FEM_ELEMENT1::iterator j=lst_fem_element1.begin();
        FEM_ELEMENT1* femsegment=(*j).second;
        MG_IDENTIFICATEUR* id=femsegment;
        LISTE_ENTITE::iterator i=lst_entite.find(id);
        lst_entite.erase(i);
        lst_fem_element1.erase(j);
        delete femsegment;
    }
}


FEM_ELEMENT1* FEM_MAILLAGE::get_premier_element1(LISTE_FEM_ELEMENT1::iterator &  it)
{
    it = lst_fem_element1.begin();
    if (it == lst_fem_element1.end())
        return NULL;
    return it->second ;
}

FEM_ELEMENT1* FEM_MAILLAGE::get_suivant_element1(LISTE_FEM_ELEMENT1::iterator &  it)
{
    it++;
    if (it==lst_fem_element1.end())
        return NULL;
    return it->second;
}


// ENTITE FEM_ELEMENT2


int FEM_MAILLAGE::ajouter_fem_element2(FEM_ELEMENT2 *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_ELEMENT2*> tmp(femtriangle->get_id(),femtriangle);
    lst_fem_element2.insert(tmp);
    return OK;
}

FEM_ELEMENT2* FEM_MAILLAGE::get_fem_element2id(unsigned long num)
{
    LISTE_FEM_ELEMENT2::iterator i=lst_fem_element2.find(num);
    if (i==lst_fem_element2.end())
    {
//      afficheur << INEXISTE << enderr;
        return NULL;
    }
    return ((*i).second);
}

unsigned int FEM_MAILLAGE::get_nb_fem_element2(void)
{
    return lst_fem_element2.size();
}


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


int FEM_MAILLAGE::supprimer_fem_element2id(unsigned long num)
{
    FEM_ELEMENT2* femtriangle=get_fem_element2id(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_ELEMENT2::iterator j=lst_fem_element2.find(num);
    lst_fem_element2.erase(j);
    delete femtriangle;
    return OK;
}

int FEM_MAILLAGE::supprimer_fem_element2(unsigned int num)
{
    FEM_ELEMENT2* femtriangle=get_fem_element2(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_ELEMENT2::iterator j=lst_fem_element2.find(num);
    lst_fem_element2.erase(j);
    delete femtriangle;
    return OK;
}

void FEM_MAILLAGE::supprimer_tout_fem_element2(void)
{
    while (get_nb_fem_element2()!=0)
    {
        LISTE_FEM_ELEMENT2::iterator j=lst_fem_element2.begin();
        FEM_ELEMENT2* femtriangle=(*j).second;
        MG_IDENTIFICATEUR* id=femtriangle;
        LISTE_ENTITE::iterator i=lst_entite.find(id);
        lst_entite.erase(i);
        lst_fem_element2.erase(j);
        delete femtriangle;
    }
}


FEM_ELEMENT2* FEM_MAILLAGE::get_premier_element2(LISTE_FEM_ELEMENT2::iterator &  it)
{
    it = lst_fem_element2.begin();
    if (it == lst_fem_element2.end())
        return NULL;
    return it->second ;
}

FEM_ELEMENT2* FEM_MAILLAGE::get_suivant_element2(LISTE_FEM_ELEMENT2::iterator &  it)
{
    it++;
    if (it==lst_fem_element2.end())
        return NULL;
    return it->second;
}


// ENTITE FEM_TETRA


int FEM_MAILLAGE::ajouter_fem_element3(FEM_ELEMENT3 *femele)
{
    gest->recherche_bonid(*femele);
    MG_IDENTIFICATEUR *id=femele;
    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_ELEMENT3*> tmp(femele->get_id(),femele);
    lst_fem_element3.insert(tmp);
    return OK;
}

FEM_ELEMENT3* FEM_MAILLAGE::get_fem_element3id(unsigned long num)
{
    LISTE_FEM_ELEMENT3::iterator i=lst_fem_element3.find(num);
    if (i==lst_fem_element3.end())
    {
//      afficheur << INEXISTE << enderr;
        return NULL;
    }
    return ((*i).second);
}

unsigned int FEM_MAILLAGE::get_nb_fem_element3(void)
{
    return lst_fem_element3.size();
}


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


int FEM_MAILLAGE::supprimer_fem_element3id(unsigned long num)
{
    FEM_ELEMENT3* femele=get_fem_element3id(num);
    if (femele==NULL)
    {
//      afficheur << INEXISTE2 << enderr;
        return FAIL;
    }
    MG_IDENTIFICATEUR* id=femele;
    LISTE_ENTITE::iterator i=lst_entite.find(id);
    lst_entite.erase(i);
    LISTE_FEM_ELEMENT3::iterator j=lst_fem_element3.find(num);
    lst_fem_element3.erase(j);
    delete femele;
    return OK;
}

int FEM_MAILLAGE::supprimer_fem_element3(unsigned int num)
{
    FEM_ELEMENT3* femele=get_fem_element3(num);
    if (femele==NULL)
    {
//      afficheur << INEXISTE2 << enderr;
        return FAIL;
    }
    MG_IDENTIFICATEUR* id=femele;
    LISTE_ENTITE::iterator i=lst_entite.find(id);
    lst_entite.erase(i);
    LISTE_FEM_ELEMENT3::iterator j=lst_fem_element3.find(num);
    lst_fem_element3.erase(j);
    delete femele;
    return OK;
}

void FEM_MAILLAGE::supprimer_tout_fem_element3(void)
{
    while (get_nb_fem_element3()!=0)
    {
        LISTE_FEM_ELEMENT3::iterator j=lst_fem_element3.begin();
        FEM_ELEMENT3* femele=(*j).second;
        MG_IDENTIFICATEUR* id=femele;
        LISTE_ENTITE::iterator i=lst_entite.find(id);
        lst_entite.erase(i);
        lst_fem_element3.erase(j);
        delete femele;
    }
}


FEM_ELEMENT3* FEM_MAILLAGE::get_premier_element3(LISTE_FEM_ELEMENT3::iterator &  it)
{
    it = lst_fem_element3.begin();
    if (it == lst_fem_element3.end())
        return NULL;
    return it->second ;
}

FEM_ELEMENT3* FEM_MAILLAGE::get_suivant_element3(LISTE_FEM_ELEMENT3::iterator &  it)
{
    it++;
    if (it==lst_fem_element3.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);
    std::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_element1(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_element2(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_element3(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);
    std::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_element1(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_element2(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_element3(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)
        {
            MC_GESTIONNAIRE mtgest(mg_geometrie->get_gest_materiau());
            MC_MATERIAU* mat=mtgest.get_materiau(mgvol->get_num_materiau());
            for (int m=0; m<mat->get_nb_propriete();m++)
            {
                MC_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)
        {
            MC_GESTIONNAIRE mtgest(mg_geometrie->get_gest_materiau());
            MC_MATERIAU* mat=mtgest.get_materiau(mgcoque->get_num_materiau());
            for (int m=0; m<mat->get_nb_propriete();m++)
            {
                MC_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)
        {
            MC_GESTIONNAIRE mtgest(mg_geometrie->get_gest_materiau());
            MC_MATERIAU* mat=mtgest.get_materiau(mgpoutre->get_num_materiau());

            for (int m=0; m<mat->get_nb_propriete();m++)
            {
                MC_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_ELEMENT3::iterator i=lst_fem_element3.begin();i!=lst_fem_element3.end();i++)
    {
        FEM_ELEMENT3* 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_ELEMENT2* femtriangle=(FEM_ELEMENT2*)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_ELEMENT1* femseg = (FEM_ELEMENT1*)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_ELEMENT1* femseg = (FEM_ELEMENT1*)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_ELEMENT1* femseg = (FEM_ELEMENT1*)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_ELEMENT2* femtriangle=(FEM_ELEMENT2*)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_ELEMENT2* femtriangle=(FEM_ELEMENT2*)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::desactive_deforme(void)
{
    deforme=0;
}


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:	309
Committed:	Tue Jan 24 15:25:16 2012 UTC (13 years, 3 months ago) by francois
File size:	60454 byte(s)
Log Message:	Preparation du maillage FEM pour inclure de nouveaux élements non necessairement triangle/tetra