orth.c

Go to the documentation of this file.

/*
%M This modules contains the ORTHOTROPIC material sub-class methods and 
   definitions
%a Joao Luiz Elias Campos.
%d December 12nd, 1999.
%r $Id: orth.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.

   Modificacao: 28/04/2005    Alexandre A. Del Savio
     Foram substituídas todas as alocações dinâmicas feitas com malloc 
     por calloc.

*/

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

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


/* Parent methods
 */
void MaterialTimeStep( sMaterial *, double * );
void MaterialDensity ( sMaterial *, double * );


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

/*
%T ORTHOTROPIC material data definition
*/
typedef struct _orthdata
{
 double Ex;            /* Young modulos in the X direction       */
 double Ey;            /* Young modulos in the Y direction       */
 double Ez;            /* Young modulos in the Z direction       */
 double Nux;           /* Poisson coefficient in the X direction */
 double Nuy;           /* Poisson coefficient in the Y direction */
 double Nuz;           /* Poisson coefficient in the Z direction */
} sOrthData;


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


/*
** ------------------------------------------------------------------------
** Subclass methods:
*/
static void OrthotropicNew        ( int, sMaterial ** );
static void OrthotropicFree       ( sMaterial * );
static void OrthotropicRead       ( sMaterial * );
static void OrthotropicEParameter ( sMaterial *, double * );
static void OrthotropicNuParameter( sMaterial *, double * );
static void OrthotropicCMatrix    ( sMaterial *, double [6][6] );


/*
%F This method allocates memory for a ORTHOTROPIC material and fills its 
   data with the specific values.
%i Material label (number)
%o Material descriptor.
*/
static void OrthotropicNew( int label, sMaterial **mat )
{
 sOrthData *data = 0L;

/* Get memory for the material descriptor
 */
 (*mat) = (sMaterial *)calloc(1, sizeof(sMaterial));

/* Get memory for the ORTHOTROPIC material data
 */
 data = (sOrthData *)calloc(1, sizeof(sOrthData));

/* Fill ORTHOTROPIC material data
 */
 data->Ex  = 0.0;
 data->Ey  = 0.0;
 data->Ez  = 0.0;
 data->Nux = 0.0;
 data->Nuy = 0.0;
 data->Nuz = 0.0;

/* Fill material descriptor
 */
 (*mat)->type  = ORTHOTROPIC;
 (*mat)->label = label;
 (*mat)->Gamma = 0.0;
 (*mat)->data  = (void *)data;

/* Add to the material list
 */
 MatList[label-1] = (*mat);

} /* End of OrthotropicNew */


/*
%F This method frees the ORTHOTROPIC material data for a given material.
%i Material descriptor.
*/
static void OrthotropicFree( sMaterial *mat )
{
 sOrthData *data = 0L;

/* Get ORTHOTROPIC material data
 */
 data = (sOrthData *)mat->data;

/* Release allocated memory
 */
 free( data );

/* Reset material data
 */
 mat->data = 0L;

} /* End of OrthotropicFree */


/*
%F This method reads the ORTHOTROPIC material information.
%i Material descriptor.
*/
static void OrthotropicRead( sMaterial *mat )
{
 sOrthData *data = 0L;
 double     ex, ey, ez, nux, nuy, nuz;

/* Get the material data
 */
 data = (sOrthData *)mat->data;

/* Read the material parameters
 */
 fscanf( nf, "%lf %lf %lf %lf %lf %lf", &ex, &ey, &ez, &nux, &nuy, &nuz );

/* Fill material information
 */
 data->Ex  = ex;
 data->Ey  = ey;
 data->Ez  = ez;
 data->Nux = nux;
 data->Nuy = nuy;
 data->Nuz = nuz;

} /* End of OrthotropicRead */


/*
%F
*/
static void OrthotropicEParameter( sMaterial *mat, double *e )
{
 sOrthData *data = 0L;

/* Get material descriptor
 */ 
 data = (sOrthData *)mat->data;

/* Get Young modules parameters
 */
 (*e) = 0.5 * (data->Ex + data->Ey);

} /* End of OrthotropicEParameter */


/*
%F
*/
static void OrthotropicNuParameter( sMaterial *mat, double *nu )
{
 sOrthData *data = 0L;

/* Get material descriptor
 */ 
 data = (sOrthData *)mat->data;

/* Get Poisson coefficient parameters
 */
 (*nu) = 0.5 * (data->Nux + data->Nuy);

} /* End of OrthotropicNuParameter */


/*
%F This function computes the material constitutive matrix
*/
static void OrthotropicCMatrix( sMaterial *mat, double cm[6][6] )
{
 int        i, j;
 sOrthData *data = 0L;
 double     ex, ey, ez, nux, nuy, nuz;
 double     a, b, c, d, e;

/* Intialize matrix
 */
 for( i = 0; i < 6; i++ )
  for( j = 0; j < 6; j++ )
   cm[i][j] = 0.0;

/* Get material descriptor
 */
 data = (sOrthData *)mat->data;

/* Get element parameters
 */
 ex  = data->Ex;
 ey  = data->Ey;
 ez  = data->Ez;
 nux = data->Nux;
 nuy = data->Nuy;
 nuz = data->Nuz;

/* Compute matrix elements
 */
 if( NDof == 3 )
 {
  cm[0][0] = 
  cm[1][1] =
  cm[2][2] = 
  cm[1][0] =
  cm[2][0] = 
  cm[0][1] = 
  cm[2][1] =
  cm[0][2] =
  cm[1][2] = 
  cm[3][3] =
  cm[4][4] =
  cm[5][5] = 1.0;
 }
 else
 {
  a =  (1.0 - (nux * nuz))/ex;
  b = -(nuy * (1.0 + nuz))/ey;
  c = -(nuz * (1.0 + nuz))/ex;
  d =  (1.0 - (nuy * nuz))/ey;
  e = (a * d) - (b * c);
  cm[0][0] = d/e;
  cm[1][1] = a/e;
  cm[0][1] = -b/e;
  cm[1][0] = -c/e;
  cm[2][2] = (ex * ey)/(ey + (ex * (1.0 + (2.0 * nux))));
 }
 
} /* End of OrthotropicCMatrix */


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

/*
%F This function initializes the sub-class "ORTHOTROPIC".
*/
void OrthotropicInit( void );
void OrthotropicInit( void )
{
/* Initialize ORTHOTROPIC material sub-class
 */
 MatClass[ORTHOTROPIC].new       = OrthotropicNew;
 MatClass[ORTHOTROPIC].free      = OrthotropicFree;
 MatClass[ORTHOTROPIC].read      = OrthotropicRead;
 MatClass[ORTHOTROPIC].epar      = OrthotropicEParameter;
 MatClass[ORTHOTROPIC].nupar     = OrthotropicNuParameter;
 MatClass[ORTHOTROPIC].cmatrix   = OrthotropicCMatrix;
 MatClass[ORTHOTROPIC].updatestr = 0L;
 MatClass[ORTHOTROPIC].updatepar = 0L;
 MatClass[ORTHOTROPIC].timestep  = MaterialTimeStep;
 MatClass[ORTHOTROPIC].density   = MaterialDensity;
 MatClass[ORTHOTROPIC].vstrain   = 0L;

} /* End of OrthotropicInit */


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

---