00001 /*========================================================================= 00002 00003 Program: Insight Segmentation & Registration Toolkit 00004 Module: $RCSfile: itkFEMElement1DStress.h,v $ 00005 Language: C++ 00006 Date: $Date: 2003-09-10 14:29:39 $ 00007 Version: $Revision: 1.6 $ 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 __itkFEMElement1DStress_h 00019 #define __itkFEMElement1DStress_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 Element1DStress : public TBaseClass 00046 { 00047 FEM_ABSTRACT_CLASS(Element1DStress,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 Element1DStress(); 00069 00070 00071 00072 00074 /* 00075 * Methods related to the physics of the problem. 00076 */ 00077 00081 virtual void GetStrainDisplacementMatrix(MatrixType& B, const MatrixType& shapeDgl) const; 00082 00086 virtual void GetMaterialMatrix(MatrixType& D) const; 00087 00093 virtual void GetStiffnessMatrix(MatrixType& Ke) const; 00094 00102 virtual unsigned int GetNumberOfDegreesOfFreedomPerNode( void ) const 00103 { return 2; } 00104 00105 00106 00107 public: 00108 00112 MaterialLinearElasticity::ConstPointer m_mat; 00113 virtual Material::ConstPointer GetMaterial(void) const { return m_mat; } 00114 virtual void SetMaterial(Material::ConstPointer mat_ ) { m_mat=dynamic_cast<const MaterialLinearElasticity*>(&*mat_); } 00116 00117 00118 }; // class Element1DStress 00119 00120 00121 00122 00123 #ifdef _MSC_VER 00124 // Declare a static dummy function to prevent a MSVC 6.0 SP5 from crashing. 00125 // I have no idea why things don't work when this is not declared, but it 00126 // looks like this declaration makes compiler forget about some of the 00127 // troubles it has with templates. 00128 static void Dummy( void ); 00129 #endif // #ifdef _MSC_VER 00130 00131 }} // end namespace itk::fem 00132 00133 #ifndef ITK_MANUAL_INSTANTIATION 00134 #include "itkFEMElement1DStress.txx" 00135 #endif 00136 00137 #endif // #ifndef __itkFEMElement1DStress_h 00138