node.c

Go to the documentation of this file.

/*
%M This modules contains the node functions and definitions
%a Joao Luiz Elias Campos.
%d September 2nd, 1998.
%r $Id: node.c,v 1.1 2004/06/22 05:29:59 joaoluiz Exp $
%w (C) COPYRIGHT 1995-1996, Eduardo Nobre Lages.
   (C) COPYRIGHT 1997-1998, Joao Luiz Elias Campos.
   All Rights Reserved
   Duplication of this program or any part thereof without the express
   written consent of the author is prohibited.
*/

/* Modified:  26-04-06  Juan Pablo Ibañez
 *   Modificada a funcão NodeApplyDisplacement, onde tipo DISPLACEMENT
 *   passa a ser setado no dof sem o controle de tolerancia. 
 *
 */

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

#include "load.h"
#include "elm.h"
#include "node.h"

/*
%V Element vector, node number and number of degrees of freedown per node
*/
sNode *NodeVector = 0L;
int    NumNodes   = 0;
int    NDof       = 2;


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

/*
%K Zero definition
*/
#ifndef ZERO
#define ZERO 0.000001
#endif


/*
** ------------------------------------------------------------------------
** Local functions:
*/

/*
** ------------------------------------------------------------------------
** Public functions:
*/

/*
%F This function initializes the prescribeds values
*/
void NodeInitPrescValues( void )
{
 int i;

/* Initialize prescribed values
 */
 for( i = 0; i < NumNodes; i++ )
 {
  NodeVector[i].dof.vpx = 0.0;
  NodeVector[i].dof.vpy = 0.0;
  NodeVector[i].dof.vpz = 0.0;
 }

} /* End of NodeInitPrescValues */


/*
%F This function release memory used by node vector
*/
void NodeFree( void )
{
/* Releases node vector
 */
 free( NodeVector );

} /* End of NodeFree */


/*
%F This function apply the nodal forces.
*/
void NodeApplyForce( int node, double fx, double fy, double fz )
{
/* Check to see if the node is not rezoned
 */
 if( NodeVector[node-1].rezone != NONE )
  return;

/* Check first degree of freedown
 */
 if( fabs(fx) > ZERO )
 {
  NodeVector[node-1].dof.x   = FORCE;
  NodeVector[node-1].dof.vpx = fx;
 }
/* Check second degree of freedown
 */
 if( fabs(fy) > ZERO )
 {
  NodeVector[node-1].dof.y   = FORCE;
  NodeVector[node-1].dof.vpy = fy;
 }
/* Check third degree of freedown
 */
 if( fabs(fz) > ZERO )
 {
  NodeVector[node-1].dof.z   = FORCE;
  NodeVector[node-1].dof.vpz = fz;
 }

} /* End of NodeApplyForce */
 

/*
%F This function apply the nodal displacements.
*/
void NodeApplyDisplacement( int node, double dx, double dy, double dz )
{
/* Check to see if the node is not rezoned
 */
 if( NodeVector[node-1].rezone != NONE )
  return;

/* Check first degree of freedown
 */
  NodeVector[node-1].dof.x   = DISPLACEMENT;
  NodeVector[node-1].dof.vpx = dx;

/* Check second degree of freedown
 */
  NodeVector[node-1].dof.y   = DISPLACEMENT;
  NodeVector[node-1].dof.vpy = dy;

/* Check third degree of freedown
 */
  NodeVector[node-1].dof.z   = DISPLACEMENT;
  NodeVector[node-1].dof.vpz = dz;

} /* End of NodeApplyDisplacement */
 

/*
%F This function apply the nodal velocities.
*/
void NodeApplyVelocity( int node, double vx, double vy, double vz )
{
/* Check to see if the node is not rezoned
 */
 if( NodeVector[node-1].rezone != NONE )
  return;

/* Check first degree of freedown
 */
 if( fabs(vx) > ZERO )
 {
  NodeVector[node-1].dof.x   = VELOCITY;
  NodeVector[node-1].dof.vpx = vx;
 }
/* Check second degree of freedown
 */
 if( fabs(vy) > ZERO )
 {
  NodeVector[node-1].dof.y   = VELOCITY;
  NodeVector[node-1].dof.vpy = vy;
 }
/* Check third degree of freedown
 */
 if( fabs(vz) > ZERO )
 {
  NodeVector[node-1].dof.z   = VELOCITY;
  NodeVector[node-1].dof.vpz = vz;
 }

} /* End of NodeApplyVelocity */
 

/*
%F This function initialize the force vector.
*/
void NodeInitForceVector( double *FVector )
{
 int i;

/* Initialize force vector
 */
 for( i = 0; i < NumNodes; i++ )
 {
  if( NodeVector[i].dof.x == FORCE )
   FVector[NDof*i]   = NodeVector[i].dof.vpx;
  else
   FVector[NDof*i+1] = 0.0;

  if( NodeVector[i].dof.y == FORCE )
   FVector[NDof*i]   = NodeVector[i].dof.vpy;
  else
   FVector[NDof*i+1] = 0.0;

  if( NDof == 3 )
  {
   if( NodeVector[i].dof.z == FORCE )
    FVector[NDof*i]   = NodeVector[i].dof.vpz;
   else
    FVector[NDof*i+2] = 0.0;
  }
 }

} /* End of NodeInitForceVector */


/* =======================================================  End of File  */

---