<div class="ajy"><img class="ajz" id=":8a" tabindex="0" src="https://mail.google.com/mail/u/0/images/cleardot.gif" alt=""></div>Hi all,<br><br>In  Modules/Numerics/FEM/src/<div id=":80">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<div class="yj6qo ajU">
<div id=":7n" class="ajR" tabindex="0">Fotis Drakopoulos<br></div></div><span class="HOEnZb adL"><font color="#888888"></font></span></div>