material.c

Go to the documentation of this file.

/*
%M This modules contains the material super class methods and definitions
%a Joao Luiz Elias Campos.
%d September 2nd, 1998.
%r $Id: material.c,v 1.1 2004/06/22 05:29:59 joaoluiz Exp $
%w (C) COPYRIGHT 1995-1996, Eduardo Nobre Lages.
   (C) COPYRIGHT 1997-1999, 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.
*/

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

#include "load.h"
#include "elm.h"
#include "material.h"

/*
%V Material class, material list and number of material in the model
*/
sMatClass   MatClass[NumMatTypes];
sMaterial **MatList      = 0L;
int         NumMaterials = 0;


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

/*
** ------------------------------------------------------------------------
** Local functions:
*/
void MohrCoulombInit      ( void );
void MohrCoulombCutOffInit( void );
void MohrCoulombNACInit   ( void );
void MohrCoulombNAVInit   ( void );
void IsotropicInit        ( void );
void IsotropicFailInit    ( void );
void OrthotropicInit      ( void );
void HiperbolicInit       ( void );
void MisesInit            ( void );
void InterfaceMatInit     ( void );
void ItfMohrCoulombMatInit( void );
void KelvinInit           ( void );
void MaxwellInit          ( void );


/*
** ------------------------------------------------------------------------
** Public functions:
*/
void MaterialTimeStep( sMaterial *, double * );
void MaterialDensity ( sMaterial *, double * );


/*
%F This function returns the material density value.
%i Material identifier.
%o Density value.
*/
void MaterialDensity( sMaterial *mat, double *gamma )
{
/* Get density value
 */
 (*gamma) = mat->Gamma;

} /* End of MaterialDensity */


/*
%F
*/
void MaterialTimeStep( sMaterial *mat, double *dt )
{
 double e, nu, gamma;

/* material properties
 */
 MatEParameter( mat, &e );
 MatNuParameter( mat, &nu );
 MatDensity( mat, &gamma ); 

/* Compute time step
 */
 (*dt) = 0.5 * sqrt(gamma * (1.0 - (2.0 * nu)) * (1.0 + nu) / (e * (1.0 - nu)));

} /* End of MaterialTimeStep */


/*
%F This function initializes the globally shared "MatClass" variable. This
   variable holds the pointers to the necessary functions which implement 
   the material class methods.
*/
void MaterialInit( void )
{
 int i;

/* Define the subclasses initialization methods
 */
 MatClass[MOHR_COULOMB_CUTOFF].init = MohrCoulombCutOffInit;
 MatClass[MOHR_COULOMB].init        = MohrCoulombInit;
 MatClass[MOHR_COULOMB_NAC].init    = MohrCoulombNACInit;
 MatClass[MOHR_COULOMB_NAV].init    = MohrCoulombNAVInit;
 MatClass[ISOTROPIC].init           = IsotropicInit;
 MatClass[ISOTROPIC_FAIL].init      = IsotropicFailInit;
 MatClass[ORTHOTROPIC].init         = OrthotropicInit;
 MatClass[HIPERBOLIC].init          = HiperbolicInit;
 MatClass[MISES].init               = MisesInit;
 MatClass[INTERFACE_MAT].init       = InterfaceMatInit;
 MatClass[ITF_MOHR_COULOMB].init    = ItfMohrCoulombMatInit;
 MatClass[KELVIN].init              = KelvinInit;
 MatClass[MAXWELL].init             = MaxwellInit;

 for( i = 0; i < NumMatTypes; i++ )
  MatInit( i );

/* Allocate memory for element list
 */
 if( MatList != 0L )
  free( MatList );

 MatList = (sMaterial **)calloc( NumMaterials, sizeof(sMaterial *) );

} /* End of MaterialInit */


/*
%F This function release memory used by element super-class.
*/
void MaterialFree( void )
{
 int i;

/* Releases sub-classes
 */
 for( i = 0; i < NumMaterials; i++ )
  MatFree( MatList[i] );

/* Release element list and reset it
 */
 free( MatList );

 MatList = 0L;

} /* End of MaterialFree */

/*
%F
*/
int MaterialIsVisco( sMaterial *mat )
{
/* Check to see if material is visco
 */
 if( mat->type != MAXWELL )
  return 0;

 return 1;

} /* End of MaterialIsVisco */


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

---