material.h

Go to the documentation of this file.

/*
%M This module contains the material class and data structure definitions
%a Joao Luiz Elias Campos.
%d September 2nd, 1998.
%r $Id: material.h,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:  Agosto-06  André Luis Muller
*    Criada a variável para material, Lambda
*/


#ifndef _MATERIAL_H
#define _MATERIAL_H

/*
%T Material type definition
*/
typedef enum _mattype
{
 MOHR_COULOMB_CUTOFF,      /* Mohr coulomb material with traction cutoff   */
 MOHR_COULOMB,             /* Mohr coulomb material                        */
 MOHR_COULOMB_NAC,         /* Mohr coulomb nonassociated constant material */
 MOHR_COULOMB_NAV,         /* Mohr coulomb nonassociated variable material */
 ISOTROPIC,                /* Isotropic material                           */
 ISOTROPIC_FAIL,           /* Isotropic fail material                      */
 ORTHOTROPIC,              /* Orthotropic material                         */
 HIPERBOLIC,               /* Hiperbolic material                          */
 MISES,                    /* Von Mises material                           */
 INTERFACE_MAT,            /* Interface material                           */
 ITF_MOHR_COULOMB,         /* Mohr coulomb iterface material               */
 KELVIN,                   /* Kelvin vico-elastic material                 */
 MAXWELL,                  /* Maxwell vico-elastic material                */
 NumMatTypes
} eMatType;

/*
%T Material data structure definition
*/
typedef struct _material
{
 eMatType  type;              /* The material type       */
 int       label;             /* Material label (number) */
 double    Gamma;             /* Material density        */
 double    Rhof;              /* Fluid density           */
 void     *data;              /* Material data           */
 double    Lambda;            /* Poro pressure parameter        */
} sMaterial;

/*
%T Material class definition
*/
typedef struct _matclass
{
 void   (*init)     ( void );
 void   (*new)      ( int, sMaterial ** );
 void   (*free)     ( sMaterial * );
 void   (*read)     ( sMaterial * );
 void   (*epar)     ( sMaterial *, double * );
 void   (*nupar)    ( sMaterial *, double * );
 void   (*cmatrix)  ( sMaterial *, double [6][6] );
 void   (*cpsmatrix)( sMaterial *, double [6][6] ); 
 void   (*updatestr)( sMaterial *, double, double *, double *, double *,
                                   double * );
 void   (*updatepar)( sMaterial *, double );
 void   (*timestep) ( sMaterial *, double * );
 void   (*density)  ( sMaterial *, double * );
 void   (*vstrain)  ( sMaterial *, double, sTensor *, sTensor * );
 void   (*rf)( sMaterial *, double , double , double * );
 
} sMatClass;

/* Material class
 */
extern sMatClass MatClass[NumMatTypes];

/* Material list
 */
extern sMaterial **MatList;

/* Number of material in the model
 */
extern int NumMaterials;

/*
** ------------------------------------------------------------------------
** Macro definitions:
*/

#define MatInit(type)                           \
        if(MatClass[type].init != 0L)           \
         (*MatClass[type].init)( )

#define MatNew(type,label,mat)                  \
        if(MatClass[type].new != 0L)            \
         (*MatClass[type].new)( label, mat )

#define MatFree(mat)                            \
        if(MatClass[mat->type].free != 0L)      \
         (*MatClass[mat->type].free)( mat )

#define MatRead(mat)                            \
        if(MatClass[mat->type].read != 0L )     \
         (*MatClass[mat->type].read)( mat )

#define MatEParameter(mat,e)                    \
        if(MatClass[mat->type].epar != 0L )     \
         (*MatClass[mat->type].epar)( mat, e )

#define MatNuParameter(mat,nu)                  \
        if(MatClass[mat->type].nupar != 0L )    \
         (*MatClass[mat->type].nupar)( mat, nu )

#define MatConstitutiveMatrix(mat,cm)           \
        if(MatClass[mat->type].cmatrix != 0L )  \
         (*MatClass[mat->type].cmatrix)( mat, cm )
         
#define MatPSMatrix(mat,cm)                             \
        if(MatClass[(mat)->type].cpsmatrix != NULL )    \
         (*MatClass[(mat)->type].cpsmatrix)( (mat), (cm) )       

#define MatUpdateStress(mat,dt,y,effd,str,def)  \
        if(MatClass[mat->type].updatestr != 0L )\
         (*MatClass[mat->type].updatestr)( mat, dt, y, effd, str, def )

#define MatUpdateParameter(mat,str)             \
        if(MatClass[mat->type].updatepar != 0L )\
         (*MatClass[mat->type].updatepar)( mat, str )

#define MatTimeStep(mat,dt)                     \
        if(MatClass[mat->type].timestep != 0L ) \
         (*MatClass[mat->type].timestep)( mat, dt )

#define MatDensity(mat,gamma)                   \
        if(MatClass[mat->type].density != 0L )  \
         (*MatClass[mat->type].density)( mat, gamma )

#define MatViscoStrain(mat,ts,str,strv)         \
        if(MatClass[mat->type].vstrain != 0L )  \
         (*MatClass[mat->type].vstrain)( mat, ts, str, strv )
         
#define MatRF(mat,sig1,sig3, rfv)       \
        if(MatClass[mat->type].rf != 0L )\
         (*MatClass[mat->type].rf)( mat, sig1, sig3, rfv )       
          

/* Public functions
*/
void MaterialInit   ( void );
void MaterialFree   ( void );
int  MaterialIsVisco( sMaterial * );

#endif

---