itkFEMElement3DStrain.h Source File

00001 /*=========================================================================
00002 
00003   Program:   Insight Segmentation & Registration Toolkit
00004   Module:    $RCSfile: itkFEMElement3DStrain.h,v $
00005   Language:  C++
00006   Date:      $Date: 2002/09/13 21:02:45 $
00007   Version:   $Revision: 1.6 $
00008 
00009   Copyright (c) 2002 Insight 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 
00027 
00028 
00029 
00044 template<class TBaseClass=Element>
00045 class Element3DStrain : public TBaseClass
00046 {
00047 FEM_ABSTRACT_CLASS(Element3DStrain,TBaseClass)
00048 public:
00049 
00050   // Repeat the required typedefs and enums from parent class
00051   typedef typename Superclass::Float Float;
00052   typedef typename Superclass::MatrixType MatrixType;
00053   typedef typename Superclass::VectorType VectorType;
00054 
00058   virtual void Read( std::istream&, void* info );
00059 
00063   virtual void Write( std::ostream& f ) const;
00064 
00068   Element3DStrain ();
00069 
00070 
00071 
00073   /*
00074    * Methods related to the physics of the problem.
00075    */
00076 
00080   virtual void GetStrainDisplacementMatrix(MatrixType& B, const MatrixType& shapeDgl) const;
00081 
00085   virtual void GetMaterialMatrix(MatrixType& D) const;
00086 
00090   virtual unsigned int GetNumberOfDegreesOfFreedomPerNode( void ) const
00091   { return 3; }
00092 
00093 public:
00094 
00098   MaterialLinearElasticity::ConstPointer m_mat;
00099   virtual Material::ConstPointer GetMaterial(void) const { return m_mat; }
00100   virtual void SetMaterial(Material::ConstPointer mat_ ) { m_mat=dynamic_cast<const MaterialLinearElasticity*>(&*mat_); }
00101 
00102 }; // class Element3DStrain
00103 
00104 
00105 
00106 
00107 #ifdef _MSC_VER
00108 // Declare a static dummy function to prevent a MSVC 6.0 SP5 from crashing.
00109 // I have no idea why things don't work when this is not declared, but it
00110 // looks like this declaration makes compiler forget about some of the
00111 // troubles it has with templates.
00112 static void Dummy( void );
00113 #endif // #ifdef _MSC_VER
00114 
00115 }} // end namespace itk::fem
00116 
00117 #ifndef ITK_MANUAL_INSTANTIATION
00118 #include "itkFEMElement3DStrain.txx"
00119 #endif
00120 
00121 #endif  // #ifndef __itkFEMElement3DStrain_h