/usr/include/InsightToolkit/Numerics/FEM/itkFEMElement1DStress.txx is in libinsighttoolkit3-dev 3.20.1-1.
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 161 162 163 164 165 166 167 168 169 170 | /*=========================================================================
Program: Insight Segmentation & Registration Toolkit
Module: itkFEMElement1DStress.txx
Language: C++
Date: $Date$
Version: $Revision$
Copyright (c) Insight Software Consortium. All rights reserved.
See ITKCopyright.txt or http://www.itk.org/HTML/Copyright.htm for details.
This software is distributed WITHOUT ANY WARRANTY; without even
the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR
PURPOSE. See the above copyright notices for more information.
=========================================================================*/
#ifndef __itkFEMElement1DStress_txx
#define __itkFEMElement1DStress_txx
#include "itkFEMElement1DStress.h"
namespace itk {
namespace fem {
template<class TBaseClass>
Element1DStress<TBaseClass>
::Element1DStress() : Superclass(), m_mat(0) {}
//////////////////////////////////////////////////////////////////////////
/**
* Methods related to the physics of the problem.
*/
template<class TBaseClass>
void
Element1DStress<TBaseClass>
::GetStrainDisplacementMatrix(MatrixType& B, const MatrixType& shapeDgl) const
{
B.set_size(1,2);
// Copy the shape function derivatives to the B matrix.
B[0][0] = shapeDgl[0][0];
B[0][1] = shapeDgl[0][1];
}
template<class TBaseClass>
void
Element1DStress<TBaseClass>
::GetMaterialMatrix(MatrixType& D) const
{
D.set_size(1,1);
D.fill(0.0);
// Material properties matrix is a scalar
D[0][0] = (m_mat->E*m_mat->A);
}
template<class TBaseClass>
void
Element1DStress<TBaseClass>
::GetStiffnessMatrix(MatrixType& Ke) const
{
const unsigned int Ndims=this->GetNumberOfSpatialDimensions();
const unsigned int Nn=this->GetNumberOfNodes();
// First we obtain the Ke by calling the parent's
// GetStiffnessMatrix function. This is the stiffness matrix
// with only one DOF per node.
Superclass::GetStiffnessMatrix(Ke);
// Calculate the nodal displacement transformation matrix according
// to the number of dimensions in global coordinate system
MatrixType T(2,2*Ndims,0.0);
VectorType d=this->GetNodeCoordinates(1)-this->GetNodeCoordinates(0);
d=d/d.magnitude();
for(unsigned int i=0;i<Ndims;i++)
{
for(unsigned int n=0;n<Nn;n++)
{
T[n][n*Ndims+i]=d[i];
}
}
// Apply the nodal displacement transformation matrix to original
// element stiffness matrix to obtain the element stiffness
// matrix in global coordinates.
Ke=T.transpose()*Ke*T;
}
template<class TBaseClass>
void
Element1DStress<TBaseClass>
::Read( std::istream& f, void* info )
{
int n;
/*
* Convert the info pointer to a usable objects
*/
ReadInfoType::MaterialArrayPointer mats=static_cast<ReadInfoType*>(info)->m_mat;
/* first call the parent's read function */
Superclass::Read(f,info);
try
{
/*
* Read and set the material pointer
*/
this->SkipWhiteSpace(f); f>>n; if(!f) goto out;
m_mat=dynamic_cast<const MaterialLinearElasticity*>( &*mats->Find(n));
}
catch ( FEMExceptionObjectNotFound e )
{
throw FEMExceptionObjectNotFound(__FILE__,__LINE__,"Element1DStress::Read()",e.m_baseClassName,e.m_GN);
}
// Check if the material object was of correct class
if(!m_mat)
{
throw FEMExceptionWrongClass(__FILE__,__LINE__,"Element1DStress::Read()");
}
out:
if( !f )
{
throw FEMExceptionIO(__FILE__,__LINE__,"Element1DStress::Read()","Error reading FEM element!");
}
}
/**
* Write the element to the output stream.
*/
template<class TBaseClass>
void
Element1DStress<TBaseClass>
::Write( std::ostream& f ) const
{
// First call the parent's write function
Superclass::Write(f);
/*
* then write the actual data (material number)
* We also add some comments in the output file
*/
f<<"\t"<<m_mat->GN<<"\t% MaterialLinearElasticity ID\n";
// check for errors
if (!f)
{
throw FEMExceptionIO(__FILE__,__LINE__,"Element1DStress::Write()","Error writing FEM element!");
}
}
#ifdef _MSC_VER
// Declare a static dummy function to prevent a MSVC 6.0 SP5 from crashing.
// I have no idea why things don't work when this is not declared, but it
// looks like this declaration makes compiler forget about some of the
// troubles it has with templates.
static void Dummy( void );
#endif // #ifdef _MSC_VER
}} // end namespace itk::fem
#endif // #ifndef __itkFEMElement1DStress_txx
|