Hi all,<br><br>In  Modules/Numerics/FEM/src/itkFEMElementBase.cxx   and in function GetStiffnessMatrix()  the element stiffness matrix (Ke) is calculated.<br>I believe that when the Linear tetrahedron element is used (class  Element3DC0LinearTetrahedronStrain)  then stiffness matrix equals to :<br>
<br>Ke = <b>(1/6)</b> * detJ * w * B.transpose() * D * B  <br>and not with  <br>Ke = detJ * w * B.transpose() * D * B<br><br>In case of the  hexahedron element the infinitesimal volume dV that we want to integrate the entries of the stiffness matrix equals to : <br>
dV = dx*dy*dz = det(J) * dr*ds*dt<br><br>However in case of thetrahedron is:<br>dV = dx*dy*dz = (1/6) * det(J) * dr*ds*dt<br><br>where det(J) is the determinant of the Jacobian matrix.<br><br>Regards<br>Fotis Drakopoulos<br>
<br><br>