/usr/include/InsightToolkit/Numerics/FEM/itkFEMElement2DStress.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 171 172 173 | /*=========================================================================
Program: Insight Segmentation & Registration Toolkit
Module: itkFEMElement2DStress.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 __itkFEMElement2DStress_txx
#define __itkFEMElement2DStress_txx
#include "itkFEMElement2DStress.h"
namespace itk {
namespace fem {
template<class TBaseClass>
Element2DStress<TBaseClass>
::Element2DStress() : Superclass(), m_mat(0) {}
//////////////////////////////////////////////////////////////////////////
/**
* Methods related to the physics of the problem.
*/
template<class TBaseClass>
void
Element2DStress<TBaseClass>
::GetStrainDisplacementMatrix(MatrixType& B, const MatrixType& shapeDgl) const
{
unsigned int p;
unsigned int Nn=this->GetNumberOfNodes();
B.set_size(3,2*Nn);
// Copy the shape function derivatives to the B matrix.
for (unsigned int i=0; i<Nn; i++) {
// Compute B index
p = i << 1;
// Compute B elements
B[0][p] = shapeDgl[0][i];
B[0][p+1] = 0;
B[1][p] = 0;
B[1][p+1] = shapeDgl[1][i];
B[2][p] = shapeDgl[1][i];
B[2][p+1] = shapeDgl[0][i];
}
}
template<class TBaseClass>
void
Element2DStress<TBaseClass>
::GetMassMatrix(MatrixType& Me) const
{
// Call the parent's get matrix function
Superclass::GetMassMatrix(Me);
// Since parent class doesn't have the material properties,
// we need to adjust Me matrix here for the density of the element.
Me=Me*m_mat->RhoC;
}
template<class TBaseClass>
void
Element2DStress<TBaseClass>
::GetMaterialMatrix(MatrixType& D) const
{
D.set_size(3,3);
/* Material properties matrix */
Float disot = (m_mat->h * m_mat->E)/(1.0 - (m_mat->nu*m_mat->nu));
D[0][0] = disot;
D[0][1] = disot * (m_mat->nu);
D[0][2] = 0.0;
D[1][0] = D[0][1];
D[1][1] = disot;
D[1][2] = 0.0;
D[2][0] = 0.0;
D[2][1] = 0.0;
D[2][2] = disot * (1.- m_mat->nu)/2.0;
}
template<class TBaseClass>
void
Element2DStress<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__,"Element2DStress::Read()",e.m_baseClassName,e.m_GN);
}
// Check if the material object was of correct class
if(!m_mat)
{
throw FEMExceptionWrongClass(__FILE__,__LINE__,"Element2DStress::Read()");
}
out:
if( !f )
{
throw FEMExceptionIO(__FILE__,__LINE__,"Element2DStress::Read()","Error reading FEM element!");
}
}
/*
* Write the element to the output stream.
*/
template<class TBaseClass>
void
Element2DStress<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__,"Element2DStress::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 __itkFEMElement2DStress_txx
|