stiffness.cpp

Go to the documentation of this file.

//********************************************************************************************************
// code:                          ###  ###  #####   ####  ##  ##
//                       ##  ##   ########  ##     ##  ## ##  ##
//                       ##  ##   ## ## ##  ###    ##     ######
//                       ##  ##   ##    ##  ##     ##  ## ##  ##
//                       #####    ##    ##  #####   ####  ##  ##
//                       ##
//
// name:           stiffness.cpp
// author(s):      Jan Novak
// language:       C, C++
// license:        This program is free software; you can redistribute it and/or modify
//                 it under the terms of the GNU Lesser General Public License as published by
//                 the Free Software Foundation; either version 2 of the License, or
//                 (at your option) any later version.
//
//                 This program is distributed in the hope that it will be useful,
//                 but WITHOUT ANY WARRANTY; without even the implied warranty of
//                 MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
//                 GNU Lesser General Public License for more details.
//
//                 You should have received a copy of the GNU Lesser General Public License
//                 along with this program; if not, write to the Free Software
//                 Foundation, Inc., 675 Mass Ave, Cambridge, MA 02139, USA.
//********************************************************************************************************

#include "stiffness.h"
#include "types.h"
#include <iostream>


namespace mumech {

void Stiffness :: giveReducedIsoStiffMatrix (double *C, double E, double nu, bool twodim)
{
  double mult = E / ( ( 1. + nu ) * ( 1. - 2. * nu ) );
  
  if (twodim) {
    C[0] = C[3] = mult * ( 1. - nu );
    C[1] = C[2] = mult * nu;
    C[4] = E / ( 1. + nu );
  }
  else {
    C[0] = C[4] = C[8] = mult * ( 1. - nu );
    C[1] = C[2] = C[3] = C[5] = C[6] = C[7] = mult * nu;
    C[9] = C[10] = C[11] = E / ( 1. + nu );
  }
}
void Stiffness :: giveReducedIsoComplMatrix (double M[12], double E, double nu)
{
  double mult = 1./E;
  
  M[0] = M[4] = M[8] = mult;
  M[1] = M[2] = M[3] = M[5] = M[6] = M[7] = -mult * nu;
  //mandel notation
  M[9] = M[10] = M[11] =  2. * mult * ( 1. + nu );
}//end of function: giveReducedIsoComplMatrix( )

} // end of namespace mumech

/*end of file*/