itkFEMElement3DStrain.h

Go to the documentation of this file.
00001 /*=========================================================================
00002 
00003   Program:   Insight Segmentation & Registration Toolkit
00004   Module:    $RCSfile: itkFEMElement3DStrain.h,v $
00005   Language:  C++
00006   Date:      $Date: 2009-01-29 20:09:13 $
00007   Version:   $Revision: 1.8 $
00008 
00009   Copyright (c) Insight Software Consortium. All rights reserved.
00010   See ITKCopyright.txt or http://www.itk.org/HTML/Copyright.htm for details.
00011 
00012      This software is distributed WITHOUT ANY WARRANTY; without even 
00013      the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR 
00014      PURPOSE.  See the above copyright notices for more information.
00015 
00016 =========================================================================*/
00017 
00018 #ifndef __itkFEMElement3DStrain_h
00019 #define __itkFEMElement3DStrain_h
00020 
00021 #include "itkFEMElementBase.h"
00022 #include "itkFEMMaterialLinearElasticity.h"
00023 
00024 namespace itk {
00025 namespace fem {
00026 
00041 template<class TBaseClass=Element>
00042 class Element3DStrain : public TBaseClass
00043 {
00044 FEM_ABSTRACT_CLASS(Element3DStrain,TBaseClass)
00045 public:
00046 
00047   // Repeat the required typedefs and enums from parent class
00048   typedef typename Superclass::Float      Float;
00049   typedef typename Superclass::MatrixType MatrixType;
00050   typedef typename Superclass::VectorType VectorType;
00051 
00055   virtual void Read( std::istream&, void* info );
00056 
00060   virtual void Write( std::ostream& f ) const;
00061 
00065   Element3DStrain ();
00066 
00068 
00075   virtual void GetStrainDisplacementMatrix(MatrixType& B, const MatrixType& shapeDgl) const;
00076 
00080   virtual void GetMaterialMatrix(MatrixType& D) const;
00081 
00085   virtual unsigned int GetNumberOfDegreesOfFreedomPerNode( void ) const
00086   { return 3; }
00087 
00088 public:
00089 
00093   MaterialLinearElasticity::ConstPointer m_mat;
00094   virtual Material::ConstPointer GetMaterial(void) const { return m_mat; }
00095   virtual void SetMaterial(Material::ConstPointer mat_ ) { m_mat=dynamic_cast<const MaterialLinearElasticity*>(&*mat_); }
00097 
00098 }; // class Element3DStrain
00099 
00100 #ifdef _MSC_VER
00101 // Declare a static dummy function to prevent a MSVC 6.0 SP5 from crashing.
00102 // I have no idea why things don't work when this is not declared, but it
00103 // looks like this declaration makes compiler forget about some of the
00104 // troubles it has with templates.
00105 static void Dummy( void );
00106 #endif // #ifdef _MSC_VER
00107 
00108 }} // end namespace itk::fem
00109 
00110 #ifndef ITK_MANUAL_INSTANTIATION
00111 #include "itkFEMElement3DStrain.txx"
00112 #endif
00113 
00114 #endif  // #ifndef __itkFEMElement3DStrain_h
00115