/usr/include/ITK-4.5/itkFEMElement2DC1Beam.h is in libinsighttoolkit4-dev 4.5.0-3.
This file is owned by root:root, with mode 0o644.
The actual contents of the file can be viewed below.
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 106 107 108 109 110 111 112 113 114 115 116 117 118 119 120 121 122 123 124 125 126 127 128 129 130 131 132 133 134 135 136 137 138 139 140 141 142 143 144 145 146 147 148 149 150 151 152 153 154 155 156 157 158 159 160 | /*=========================================================================
*
* Copyright Insight Software Consortium
*
* Licensed under the Apache License, Version 2.0 (the "License");
* you may not use this file except in compliance with the License.
* You may obtain a copy of the License at
*
* http://www.apache.org/licenses/LICENSE-2.0.txt
*
* Unless required by applicable law or agreed to in writing, software
* distributed under the License is distributed on an "AS IS" BASIS,
* WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
* See the License for the specific language governing permissions and
* limitations under the License.
*
*=========================================================================*/
#ifndef __itkFEMElement2DC1Beam_h
#define __itkFEMElement2DC1Beam_h
#include "itkFEMElementStd.h"
#include "itkFEMLoadElementBase.h"
#include "itkFEMMaterialLinearElasticity.h"
namespace itk
{
namespace fem
{
/**
* \class Element2DC1Beam
* \brief 1D Beam (spring that also bends) finite element in 2D space.
*
* The Displacements at each node are modelled with two translational,
* and one rotational, degree of freedom.
* \ingroup ITKFEM
*/
class Element2DC1Beam : public ElementStd<2, 2>
{
public:
/** Standard class typedefs. */
typedef Element2DC1Beam Self;
typedef ElementStd<2, 2> TemplatedParentClass;
typedef TemplatedParentClass Superclass;
typedef SmartPointer<Self> Pointer;
typedef SmartPointer<const Self> ConstPointer;
/** Method for creation through the object factory. */
itkSimpleNewMacro(Self);
/** Run-time type information (and related methods). */
itkTypeMacro(Element2DC1Beam, TemplatedParentClass);
/** CreateAnother method will clone the existing instance of this type,
* including its internal member variables. */
virtual::itk::LightObject::Pointer CreateAnother(void) const;
/**
* Default constructor only clears the internal storage
*/
Element2DC1Beam();
/**
* Construct an element by specifying two nodes and material
*/
Element2DC1Beam(Node::ConstPointer n1_, Node::ConstPointer n2_, Material::ConstPointer mat_);
// ////////////////////////////////////////////////////////////////////////
/*
* Methods related to the physics of the problem.
*/
/** Get the Stiffness matrix */
virtual void GetStiffnessMatrix(MatrixType & Ke) const;
/** Get the Mass matrix */
virtual void GetMassMatrix(MatrixType & Me) const;
/** Get the Strain Displacement matrix */
virtual void GetStrainDisplacementMatrix(MatrixType &, const MatrixType &) const
{
}
/** Get the Material matrix */
virtual void GetMaterialMatrix(MatrixType &) const
{
}
// ////////////////////////////////////////////////////////////////////////
/**
* Methods related to numeric integration
*/
enum { DefaultIntegrationOrder = 1 };
/** Get the Integration point and weight */
virtual void GetIntegrationPointAndWeight(unsigned int i, VectorType & pt, Float & w, unsigned int order = 0) const;
virtual unsigned int GetNumberOfIntegrationPoints(unsigned int order) const;
// ////////////////////////////////////////////////////////////////////////
/**
* Methods related to the geometry of an element
*/
/** Return the shape functions used to interpolate across the element */
virtual VectorType ShapeFunctions(const VectorType & pt) const;
/** Return the shape functions derivatives in the shapeD matrix */
virtual void ShapeFunctionDerivatives(const VectorType & pt, MatrixType & shapeD) const;
/** Convert from global to local coordinates */
virtual bool GetLocalFromGlobalCoordinates(const VectorType &, VectorType &) const
{
return false;
}
/** Return the determinate of the Jacobian */
virtual Float JacobianDeterminant(const VectorType & pt, const MatrixType *pJ) const;
/** Get the degrees of freedom for each node */
virtual unsigned int GetNumberOfDegreesOfFreedomPerNode(void) const
{
return 3;
}
/**
* Get/Set the material properties for the element
*/
virtual Material::ConstPointer GetMaterial(void) const
{
return dynamic_cast<const Material *>(m_mat);
}
virtual void SetMaterial(Material::ConstPointer mat_)
{
m_mat =
dynamic_cast<const MaterialLinearElasticity *>( mat_.GetPointer() );
}
/** No edges to populate in this class */
virtual void PopulateEdgeIds(void) { /* empty */ }
protected:
virtual void PrintSelf(std::ostream& os, Indent indent) const;
private:
/**
* Pointer to material properties of the element
*/
const MaterialLinearElasticity *m_mat;
};
}
} // end namespace itk::fem
#endif // #ifndef __itkFEMElement2DC1Beam_h
|