nfr.c

Vá para a documentação deste arquivo.

/*
** nfr.c - Reads a two-dimensional frame model from a neutral format file.
**
** Description:
** ---------------------------------------------------------------
**
** This file contains functions to read a two-dimensional frame
** model dada from a file and to store these data in the 
** adopted data structure.  
** This data structure is described in file "crs.c".
** The model  data is read from a stream file in the neutral file
** format adopted in Tecgraf and Civil Engineering Dept. of 
** PUC-Rio. Please, refer to specific documentation on this 
** format.
**
** The API of this module is described in file "nfr.h".
**
** ---------------------------------------------------------------
**
** Created:    07-Sep-2003    Luiz F. Martha
**
*/

/*
** ---------------------------------------------------------------
** Global variables and symbols:
*/
#include <stdlib.h>
#include <stdio.h>
#include <string.h>

#include "crs.h"
#include "nfr.h"

/*
** ---------------------------------------------------------------
** Local variables and symbols:
*/

typedef struct _matdata {
   double    E;      /* Young's modulus           */
   double    Nu;     /* Poisson's ratio           */
} MatData ;

static MatData *mat_data = NULL;

typedef struct _secprop {
   double Ax;
   double Ay;
   double Az;
   double Ix;
   double Iy;
   double Iz;
   double Cw;
   double Hy;
   double Hz;
} SecProp;

static SecProp *sec_prop = NULL;

typedef struct _endlib {
   int    idx;
   int    idy;
   int    idz;
   int    irx;
   int    iry;
   int    irz;
   int    jdx;
   int    jdy;
   int    jdz;
   int    jrx;
   int    jry;
   int    jrz;
} EndLib;

static EndLib *end_lib = NULL;


/*
** ---------------------------------------------------------------
** Private functions:
*/
static int  NextLabel( FILE *fd, char *label );
static void FreeData( void );
static int  ReadNodalSupports( FILE *fd );
static int  ReadMatData( FILE *fd );
static int  ReadSecProp( FILE *fd );
static int  ReadEndLib( FILE *fd );
static int  ReadUnifLoad( FILE *fd );
static void ReadBeamElems( FILE *fd );

/*=========================  NextLabel  =========================*/

static int  NextLabel( FILE *fd, char *label )
{
 int c;
 
 while( ( c = fgetc( fd ) ) != '%' )           /* finds next % */
 {
  if( c == EOF )
  return( 0 );
 }
 if( fscanf( fd, "%s", label ) != 1 )          /* scan label string */
  return( 0 );
 else
  return( 1 );
}

/*=========================  FreeData  =========================*/

static void FreeData( void )
{
/* Release memory used by auxiliary arrays.
 */
 if( mat_data )
  free( mat_data );
 mat_data = NULL;
 if( sec_prop )
  free( sec_prop );
 sec_prop = NULL;
 if( end_lib )
  free( end_lib );
 end_lib = NULL;
}

/*========================  ReadNodalSupports  =======================*/

static int ReadNodalSupports( FILE *fd )
{
 char      label[80];             /* sections labels             */
 int       found = 0;             /* flag for found nodal data   */
 int       nnodes;                /* read number of nodes        */
 int       node_id;               /* current node id             */
 int       i;                     /* counter                     */
 int       dx, dy, dz, rx, ry, rz;/* support codes (0=free,1=fix)*/

/* Find nodal support data section in neutral file.
 */
 rewind( fd );
 while( 1 )
 {
  if( NextLabel( fd, label ) == 0 )
  {
   printf( "\n Invalid neutral file or label END doesn't exist.\n\n" );
   return( 0 );
  }
  else if( strcmp( label, "NODE.SUPPORT" ) == 0 )
  {
   found = 1;
   break;
  }
  else if( strcmp( label, "END" ) == 0 )
   break;
  else             /* not a target label */
   continue;
 }
 if( !found )
 {
  return( 1 );
 }

/* Read number of nodes of this section.
 */
 fscanf( fd, "%d", &nnodes );

/* Read node support conditions.
 */
 for( i = 0; i < nnodes; i++ )
 {
  fscanf( fd, "%d  %d %d %d %d %d %d", &node_id, &dx, &dy, &dz, &rx, &ry, &rz );
  if( rz == 1 )
   crsFixNodeRot( node_id-1 );
 }

 return( 1 );
}

/*========================  ReadMatData  =======================*/

static int ReadMatData( FILE *fd )
{
 char      label[80];             /* sections labels             */
 int       found = 0;             /* flag for found material data*/
 int       nmats;                 /* read number of materials    */
 int       mat_id;                /* material id (not used)      */
 int       i;                     /* counter                     */

/* Find isotropic material data section in neutral file.
 */
 rewind( fd );
 while( 1 )
 {
  if( NextLabel( fd, label ) == 0 )
  {
   printf( "\n Invalid neutral file or label END doesn't exist.\n\n" );
   return( 0 );
  }
  else if( strcmp( label, "MATERIAL.ISOTROPIC" ) == 0 )
  {
   found = 1;
   break;
  }
  else if( strcmp( label, "END" ) == 0 )
   break;
  else             /* not a target label */
   continue;
 }
 if( !found )
 {
  printf( "\n Could not find material data in neutral file.\n\n" );
  return( 0 );
 }

/* Read number of isotropic materials.
 */
 fscanf( fd, "%d", &nmats );

/* Allocate memory for auxiliary array of material parameters.
 */
 mat_data = (MatData *)calloc( nmats, sizeof( MatData ) );

/* Read material data.
 */
 for( i = 0; i < nmats; i++ )
 {
  fscanf( fd, "%d %lg %lg", &mat_id, &mat_data[i].E, &mat_data[i].Nu );
 }

 return( 1 );
}

/*========================  ReadSecProp  =======================*/

static int ReadSecProp( FILE *fd )
{
 char      label[80];             /* sections labels                  */
 int       found = 0;             /* flag for found cross-section data*/
 int       nsecs;                 /* read number of cross-sections    */
 int       sec_id;                /* cross-section id (not used)      */
 int       i;                     /* counter                          */

/* Find cross-section data section in neutral file.
 */
 rewind( fd );
 while( 1 )
 {
  if( NextLabel( fd, label ) == 0 )
  {
   printf( "\n Invalid neutral file or label END doesn't exist.\n\n" );
   return( 0 );
  }
  else if( strcmp( label, "SECTION.PROPERTY" ) == 0 )
  {
   found = 1;
   break;
  }
  else if( strcmp( label, "END" ) == 0 )
   break;
  else             /* not a target label */
   continue;
 }
 if( !found )
 {
  printf( "\n Could not find cross-section data in neutral file.\n\n" );
  return( 0 );
 }

/* Read number of cross-sections
 */
 fscanf( fd, "%d", &nsecs );

/* Allocate memory for auxiliary array of cross-section parameters
 */
 sec_prop = (SecProp *)calloc( nsecs, sizeof( SecProp ) );

/* Read cross-section data.
 */
 for( i = 0; i < nsecs; i++ )
 {
  fscanf( fd, "%d %lg %lg %lg %lg %lg %lg %lg %lg %lg", &sec_id,
              &sec_prop[i].Ax, &sec_prop[i].Ay, &sec_prop[i].Az, 
              &sec_prop[i].Ix, &sec_prop[i].Iy, &sec_prop[i].Iz, 
              &sec_prop[i].Cw, &sec_prop[i].Hy, &sec_prop[i].Hz );
 }

 return( 1 );
}

/*========================  ReadEndLib  =======================*/

static int ReadEndLib( FILE *fd )
{
 char      label[80];             /* sections labels             */
 int       found = 0;             /* flag for found endlib data  */
 int       nlibs;                 /* read number of endlibs      */
 int       lib_id;                /* endlib id (not used)        */
 int       i;                     /* counter                     */

/* Find beam end-liberation data section in neutral file.
 */
 rewind( fd );
 while( 1 )
 {
  if( NextLabel( fd, label ) == 0 )
  {
   printf( "\n Invalid neutral file or label END doesn't exist.\n\n" );
   return( 0 );
  }
  else if( strcmp( label, "BEAM.END.LIBERATION" ) == 0 )
  {
   found = 1;
   break;
  }
  else if( strcmp( label, "END" ) == 0 )
   break;
  else             /* not a target label */
   continue;
 }
 if( !found )
 {
  return( 1 );
 }

/* Read number of beam end-liberations.
 */
 fscanf( fd, "%d", &nlibs );

/* Allocate memory for auxiliary array of beam end-liberations.
 */
 end_lib = (EndLib *)calloc( nlibs, sizeof( EndLib ) );

/* Read end-liberation data.
 */
 for( i = 0; i < nlibs; i++ )
 {
  fscanf( fd, "%d %d %d %d %d %d %d %d %d %d %d %d %d", &lib_id, 
              &end_lib[i].idx, &end_lib[i].idy, &end_lib[i].idz,
              &end_lib[i].irx, &end_lib[i].iry, &end_lib[i].irz,
              &end_lib[i].jdx, &end_lib[i].jdy, &end_lib[i].jdz,
              &end_lib[i].jrx, &end_lib[i].jry, &end_lib[i].jrz );
 }

 return( 1 );
}

/*========================  ReadUnifLoad  =======================*/

static int ReadUnifLoad( FILE *fd )
{
 char      label[80];        /* sections labels                    */
 int       found = 0;        /* flag for found unif.load data      */
 int       nunfs;            /* read number of uniform loads       */
 int       member_id;        /* current member id                  */
 int       dir;              /* load direction (local or global)   */
 double    a1, a2;           /* ratio of starting and end load dist*/
 double    q[3];             /* load force components              */
 double    m[3];             /* load moment components             */
 int       i;                /* counter                            */

/* Find beam uniform load data section in neutral file.
 */
 rewind( fd );
 while( 1 )
 {
  if( NextLabel( fd, label ) == 0 )
  {
   printf( "\n Invalid neutral file or label END doesn't exist.\n\n" );
   return( 0 );
  }
  else if( strcmp( label, "LOAD.CASE.BEAM.UNIFORM" ) == 0 )
  {
   found = 1;
   break;
  }
  else if( strcmp( label, "END" ) == 0 )
   break;
  else             /* not a target label */
   continue;
 }
 if( !found )
 {
  printf( "\n Could not find beam uniform load data in neutral file.\n\n" );
  return( 0 );
 }

/* Read number of beam uniform loads
 */
 fscanf( fd, "%d", &nunfs );

/* Read beam uniform load data.
 */
 for( i = 0; i < nunfs; i++ )
 {
  fscanf( fd, "%d", &member_id );
  fscanf( fd, "%d", &dir );
  fscanf( fd, "%lg %lg", &a1, &a2 );
  fscanf( fd, "%lg %lg %lg", &q[0], &q[1], &q[2] );
  fscanf( fd, "%lg %lg %lg", &m[0], &m[1], &m[2] );
  crsMemberUnifLoad( member_id-1, q[1] );
 }

 return( 1 );
}

/*=====================  ReadBeamElems  =======================*/

static void ReadBeamElems( FILE *fd )
{
 int       nelems;                /* read number of elements     */
 int       dummy_id;              /* dummy element id            */
 int       elem_id;               /* current element id          */
 int       node0_id;              /* element initial node id     */
 int       node1_id;              /* element final node id       */
 int       mat_id;                /* element material id         */
 int       sec_id;                /* element cross-section id    */
 int       ori_id;                /* non used elem. orientat. id */
 int       lib_id;                /* element end-liberation id   */
 int       i;                     /* counters                    */

 fscanf( fd, "%d", &nelems );

 for( i = 0; i < nelems; i++ )
 {
/* Read element number, material index, section index, orientation 
 * index (not considered), and end-liberation index.
 */
  fscanf( fd, "%d %d %d %d %d", &dummy_id, &mat_id, &sec_id, &ori_id, &lib_id );

/* Read element node incidence.
 * It is assumed that the read incidences contains numbers 
 * (index is the number decreased by one).
 */
  fscanf( fd, "%d %d", &node0_id, &node1_id );

/* Add this element to project data structure.
 */
  elem_id = crsAddMember( node0_id-1, node1_id-1 );

/* Set the elasticity modulus of this element.
 */
  crsMemberElasticity( elem_id, mat_data[mat_id-1].E );

/* Set the moment of inertia of this element.
 */
  crsMemberInertia( elem_id, sec_prop[sec_id-1].Iz );

/* Set the end liberations of this member.
 */
  if( end_lib[lib_id-1].irz == 0 )
   crsFreeMemberEndRot( elem_id, 0 );
  if( end_lib[lib_id-1].jrz == 0 )
   crsFreeMemberEndRot( elem_id, 1 );
 }
}


/*
** ---------------------------------------------------------------
** Entry points begin here:
*/

/*=======================  nfrCrsModel  =======================*/

void *nfrCrsModel( FILE *fd )
{
 int       num_nodes = 0;
 int       num_elems = 0;
 char      label[80];             /* sections labels             */
 void      *model;                /* created model */

/* Read total number of nodes and elements.
 */
 rewind( fd );
 while( 1 )
 {
  if( NextLabel( fd, label ) == 0 )
  {
   printf( "\n Invalid neutral file or label END doesn't exist.\n\n" );
   return( NULL );
  }
  else if( strcmp( label, "NODE" ) == 0 )
  {
   fscanf( fd,"%d", &num_nodes );
   if( (num_nodes > 0) && (num_elems > 0) )
    break;
  }
  else if( strcmp( label, "ELEMENT" ) == 0 )
  {
   fscanf( fd,"%d", &num_elems );
   if( (num_nodes > 0) && (num_elems > 0) )
    break;
  }
  else if( strcmp( label, "END" ) == 0 )
   break;
  else             /* not a target label */
   continue;
 }

 if( (num_nodes < 1) && (num_elems < 1) )
 {
  printf( "\n Could not find mesh dimensions in neutral file.\n\n" );
  return( NULL );
 }

/* Initialiaze project database and return the created model pointer.
 */
 model = crsInitModel( num_nodes, num_elems );

 return( model );
}

/*========================  nfrCrsNodes  =======================*/

int nfrCrsNodes( FILE *fd )
{
 char      label[80];             /* sections labels             */
 int       found = 0;             /* flag for found nodal data   */
 int       nnodes;                /* read number of nodes        */
 int       dummy_id;              /* dummy node id               */
 int       node_id;               /* current node id             */
 int       i;                     /* counter                     */
 double    x, y, z;               /* coordinates of a node       */

/* Find nodal coordinate data section in neutral file.
 */
 rewind( fd );
 while( 1 )
 {
  if( NextLabel( fd, label ) == 0 )
  {
   printf( "\n Invalid neutral file or label END doesn't exist.\n\n" );
   return( 0 );
  }
  else if( strcmp( label, "NODE.COORD" ) == 0 )
  {
   found = 1;
   break;
  }
  else if( strcmp( label, "END" ) == 0 )
   break;
  else             /* not a target label */
   continue;
 }
 if( !found )
 {
  printf( "\n Could not find nodal coordinate data in neutral file.\n\n" );
  return( 0 );
 }

/* Read number of nodes of this section.
 */
 fscanf( fd, "%d", &nnodes );

/* Read node number (not considered) and nodal coordinates.
 * It is assumed that the nodes are numbered in consecutive
 * order.
 */
 for( i = 0; i < nnodes; i++ )
 {
  fscanf( fd, "%d %lf %lf %lf", &dummy_id, &x, &y, &z );
  node_id = crsAddNode( x, y );
 }

/* Read nodal support data.
 */
 if( ! ReadNodalSupports( fd ) )
  return( 0 );

 return( 1 );
}

/*=======================  nfrCrsMembers  =====================*/

int nfrCrsMembers( FILE *fd )
{
 char      label[80];             /* sections labels             */
 int       found = 0;             /* flag for found element data */

/* Read arrays of material, cross-section, and end-liberation data.
 */
 if( ! ReadMatData( fd ) )
 {
  FreeData( );
  return( 0 );
 }
 if( ! ReadSecProp( fd ) )
 {
  FreeData( );
  return( 0 );
 }
 if( ! ReadEndLib( fd ) )
 {
  FreeData( );
  return( 0 );
 }

/* Find element data sections in neutral file and read data.
 */
 rewind( fd );
 while( 1 )
 {
  if( NextLabel( fd, label ) == 0 )
  {
   printf( "\n Invalid neutral file or label END doesn't exist.\n\n" );
   return( 0 );
  }
  else if( strcmp( label, "ELEMENT.BEAM" ) == 0 )
  {
   ReadBeamElems( fd );
   found = 1;
  }
  else if( strcmp( label, "END" ) == 0 )
   break;
  else             /* not a target label */
   continue;
 }
 if( !found )
 {
  printf( "\n Could not find element data in neutral file.\n\n" );
  return( 0 );
 }

/* Read beam uniform load data.
 */
 if( ! ReadUnifLoad( fd ) )
 {
  FreeData( );
  return( 0 );
 }

/* Release memory used by auxiliary arrays.
 */
 FreeData( );

 return( 1 );
}

---