ITK
4.1.0
Insight Segmentation and Registration Toolkit
|
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 __itkFEMElement2DC1Beam_h 00020 #define __itkFEMElement2DC1Beam_h 00021 00022 #include "itkFEMElementStd.h" 00023 #include "itkFEMLoadElementBase.h" 00024 #include "itkFEMMaterialLinearElasticity.h" 00025 00026 namespace itk 00027 { 00028 namespace fem 00029 { 00038 class Element2DC1Beam : public ElementStd<2, 2> 00039 { 00040 public: 00041 00043 typedef Element2DC1Beam Self; 00044 typedef ElementStd<2, 2> TemplatedParentClass; 00045 typedef TemplatedParentClass Superclass; 00046 typedef SmartPointer<Self> Pointer; 00047 typedef SmartPointer<const Self> ConstPointer; 00048 00050 itkSimpleNewMacro(Self); 00051 00053 itkTypeMacro(Element2DC1Beam, TemplatedParentClass); 00054 00057 virtual::itk::LightObject::Pointer CreateAnother(void) const; 00058 00059 00063 Element2DC1Beam(); 00064 00068 Element2DC1Beam(Node::ConstPointer n1_, Node::ConstPointer n2_, Material::ConstPointer mat_); 00069 00070 // //////////////////////////////////////////////////////////////////////// 00071 /* 00072 * Methods related to the physics of the problem. 00073 */ 00074 00076 virtual void GetStiffnessMatrix(MatrixType & Ke) const; 00077 00079 virtual void GetMassMatrix(MatrixType & Me) const; 00080 00082 virtual void GetStrainDisplacementMatrix(MatrixType &, const MatrixType &) const 00083 { 00084 } 00085 00087 virtual void GetMaterialMatrix(MatrixType &) const 00088 { 00089 } 00090 00091 // //////////////////////////////////////////////////////////////////////// 00096 enum { DefaultIntegrationOrder = 1 }; 00097 00099 virtual void GetIntegrationPointAndWeight(unsigned int i, VectorType & pt, Float & w, unsigned int order = 0) const; 00100 00101 virtual unsigned int GetNumberOfIntegrationPoints(unsigned int order) const; 00102 00103 // //////////////////////////////////////////////////////////////////////// 00109 virtual VectorType ShapeFunctions(const VectorType & pt) const; 00110 00112 virtual void ShapeFunctionDerivatives(const VectorType & pt, MatrixType & shapeD) const; 00113 00115 virtual bool GetLocalFromGlobalCoordinates(const VectorType &, VectorType &) const 00116 { 00117 return false; 00118 } 00119 00121 virtual Float JacobianDeterminant(const VectorType & pt, const MatrixType *pJ) const; 00122 00124 virtual unsigned int GetNumberOfDegreesOfFreedomPerNode(void) const 00125 { 00126 return 3; 00127 } 00128 00132 virtual Material::ConstPointer GetMaterial(void) const 00133 { 00134 return dynamic_cast<const Material *>(m_mat); 00135 } 00136 00137 virtual void SetMaterial(Material::ConstPointer mat_) 00138 { 00139 m_mat = 00140 dynamic_cast<const MaterialLinearElasticity *>( mat_.GetPointer() ); 00141 } 00142 00144 virtual void PopulateEdgeIds(void) { /* empty */ } 00145 protected: 00146 virtual void PrintSelf(std::ostream& os, Indent indent) const; 00147 00148 private: 00149 00153 const MaterialLinearElasticity *m_mat; 00154 00155 }; 00156 00157 } 00158 } // end namespace itk::fem 00159 00160 #endif // #ifndef __itkFEMElement2DC1Beam_h 00161