ITK  4.1.0
Insight Segmentation and Registration Toolkit
itkFEMElement3DMembrane1DOF.h
Go to the documentation of this file.
00001 /*=========================================================================
00002 *
00003 * Copyright Insight Software Consortium
00004 *
00005 * Licensed under the Apache License, Version 2.0 (the "License");
00006 * you may not use this file except in compliance with the License.
00007 * You may obtain a copy of the License at
00008 *
00009 * http://www.apache.org/licenses/LICENSE-2.0.txt
00010 *
00011 * Unless required by applicable law or agreed to in writing, software
00012 * distributed under the License is distributed on an "AS IS" BASIS,
00013 * WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
00014 * See the License for the specific language governing permissions and
00015 * limitations under the License.
00016 *
00017 *=========================================================================*/
00018 
00019 #ifndef __itkFEMElement3DMembrane1DOF_h
00020 #define __itkFEMElement3DMembrane1DOF_h
00021 
00022 #include "itkFEMElementBase.h"
00023 #include "itkFEMMaterialLinearElasticity.h"
00024 
00025 namespace itk
00026 {
00027 namespace fem
00028 {
00044 template <class TBaseClass = Element>
00045 class Element3DMembrane1DOF : public TBaseClass
00046 {
00047 public:
00048 
00050   typedef Element3DMembrane1DOF    Self;
00051   typedef TBaseClass               Superclass;
00052   typedef SmartPointer<Self>       Pointer;
00053   typedef SmartPointer<const Self> ConstPointer;
00054 
00056   itkTypeMacro(Element3DMembrane1DOF, TBaseClass);
00057 
00058   // Repeat the required typedefs and enums from parent class
00059   typedef typename Superclass::Float      Float;
00060   typedef typename Superclass::MatrixType MatrixType;
00061   typedef typename Superclass::VectorType VectorType;
00062 
00066   Element3DMembrane1DOF();
00067 
00068   // ////////////////////////////////////////////////////////////////////////
00069   /*
00070    * Methods related to the physics of the problem.
00071    */
00072 
00076   virtual void GetStrainDisplacementMatrix(MatrixType & B, const MatrixType & shapeDgl) const;
00077 
00081   virtual void GetMaterialMatrix(MatrixType & D) const;
00082 
00086   void GetMassMatrix(MatrixType & Me) const;
00087 
00091   virtual unsigned int GetNumberOfDegreesOfFreedomPerNode(void) const
00092   {
00093     return 3;
00094   }
00095 
00097   virtual void GetStiffnessMatrix(MatrixType & Ke) const;
00098 
00102   virtual Material::ConstPointer GetMaterial(void) const
00103   {
00104     return dynamic_cast<const Material *>(m_Mat.GetPointer());
00105   }
00106 
00107   virtual void SetMaterial(Material::ConstPointer mat_)
00108   {
00109     m_Mat =
00110       dynamic_cast<const MaterialLinearElasticity *>( mat_.GetPointer() );
00111   }
00112 
00113 protected:
00114 
00115   virtual void PrintSelf(std::ostream& os, Indent indent) const;
00116 
00120   MaterialLinearElasticity::ConstPointer m_Mat;
00121 
00122 };  // class Element3DMembrane1DOF
00123 
00124 }
00125 }  // end namespace itk::fem
00126 
00127 #ifndef ITK_MANUAL_INSTANTIATION
00128 #include "itkFEMElement3DMembrane1DOF.hxx"
00129 #endif
00130 
00131 #endif  // #ifndef __itkFEMElement3DMembrane1DOF_h
00132