/Users/jack/Code/basso_dev/inc/NastranInput.h

Go to the documentation of this file.

#ifndef _NASTRAN_INPUT_BASSO_H_
#define _NASTRAN_INPUT_BASSO_H_

// I/O includes
#include <iostream>
#include <iomanip>
#include <fstream>

// STL includes
#include <vector>
#include <list>
#include <set>
#include <map>

// Basso includes
#include "basso.h"
#include "FEInput.h"

namespace basso {


class NastranCard 
{

public:
        NastranCard() : mFields(25) { mNumFelds=0; }
        
        bool ReadNextGrid( ifstream &inFile )    { }
        bool ReadNextElement( ifstream &inFile ) { }
        
        int NumFields() const { return mNumFields; }
        
        const int &operator ( int f ) const { }
        const Numeric &operator ( int f ) const { }
        const String &operator ( int f ) const { }
        
private:
        
        int mNumFields;
        Array<String> mFields;
        
};

class NastranInput : public FEInput {
        
public:
        

                NastranInput( ) : Input( ) { fProcessed=false; }
                NastranInput( const String &file ) : Input(file) { fProcessed=false; }
        virtual ~NastranInput() {}
        
        // MEMEBER FUNCTIONS
        
        virtual void ReadNodes( Array<Node> &nodes );
        
        virtual void AddElements( list<Element> &elem, int pid=-999 ) const;
        
        virtual void AddNodeSet( set<int> &nid, int pid ) const;

protected:
        
        int ReadNextCard( list<String> &fields, ifstream &inFile ) const;
        
        int ConvertField( int fid, const list<String> &fields ) const;
        
        Numeric ConvertField( int fid, const list<String> &fields ) const;
        
        BasisType GetBasisType( const String &etype, int nn ) const;

        int NumNodes() const;
        
protected:


                
};

int NastranInput::NumNodes() const
{
        // open Nastran file
        std::ifstream inFile;
        if ( OpenFile(inFile) ) 
                WarningMessage("NastranInput::NumNodes","could not open file");
        
        int nn=0;       
        
        while ( !GotoNextSection("GRID") )
                ++nn;
                
        CloseFile(inFile);
        
        return nn;
}

void NastranInput::ReadNodes( Array<Node> &node )  
{

        // open Nastran file
        std::ifstream inFile;
        if ( OpenFile(inFile) ) 
                WarningMessage("NastranInput::ReadNodes","could not open file");
 
        //  get the number of nodes
        int nn=NumNodes();
        if ( node.size()!=nn )
                resize( node, nn );
                
        // read the nodes
        int i=0;
        while ( !GotoNextSection("GRID") )
        {
                int id, cid;
                char card(16);
                Numeric x, y, z;
                inFile >> card >> id >> cid >> x >> y >> z;
                node[i]=Node(i,x,y,z);
                mNodeMap[id]=i++;
        }

        CloseFile(inFile);
        
}

void NastranInput::AddElements( list<Element> &elem, int pid ) const
{
        // open Nastran file
        std::ifstream inFile;
        if ( OpenFile(inFile) ) 
                WarningMessage("NastranInput::ReadElements","could not open file");
        
        WarningMessage("NastranInput::AddElements","not yet implemented");

        while ( GotoNextSection("C") )
        {
                if ( fields)
        }
        
        CloseFile(inFile);
         
}

void NastranInput::AddNodeSet( set<int> &nset, int pid ) const
{
        // open Nastran file
        std::ifstream inFile;
        if ( OpenFile(inFile) ) 
                WarningMessage("NastranInput::ReadNodeSet","could not open file");
         
        WarningMessage("NastranInput::AddNodeSet","not yet implemented");

        CloseFile(inFile);
}

        

BasisType NastranInput::GetBasisType( const String &etype, int nn ) const
{
        switch (etype) {
                
                case "CBAR":
                return kLINE2;
                
                case "CTRIA3":
                return kTRIA3;
                
                case "CQUAD4":
                return kQUAD4;
                
                case "CTETRA":
                return kTETRA4;
                
                case "CHEXA ":
                return kHEXA8;
                
                case "CPENTA":
                return kPRISM6;
                
                case "CTRIA6":
                return kTRIA6;
                
                case "CQUAD8":
                return kQUAD8;
                
                case "CTETRA":
                return kTETRA10;
                
                case "CHEXA ":
                return kHEXA20;
                
                case "CPENTA":
                return kPRISM18;
                
                case "CMASS1":
                return kPOINT1;
                
                case "CPENTA":
                return kPRISM15;
                
                default:
                WarningMessage("BasisType NastranInput::GetBasisType( int etype )","Unknown etype");
                return kNONE;
        }
}





} // end Basso namepsace 
#endif  

---