GetJacobian

PURPOSE ^

GETJACOBIAN- computes the jacobian for a tetraelem

SYNOPSIS ^

function J=GetJacobian(tetraelem,xyz_list,gauss_coord)

DESCRIPTION ^

GETJACOBIAN- computes the jacobian for a tetraelem

   Usage:
      J=GetJacobian(tetraelem,z_list,gauss_coord)

   See also GETJACOBIANDETERMINANT

CROSS-REFERENCE INFORMATION ^

This function calls: This function is called by:

SOURCE CODE ^

0001 function J=GetJacobian(tetraelem,xyz_list,gauss_coord)
0002 %GETJACOBIAN- computes the jacobian for a tetraelem
0003 %
0004 %   Usage:
0005 %      J=GetJacobian(tetraelem,z_list,gauss_coord)
0006 %
0007 %   See also GETJACOBIANDETERMINANT
0008 
0009 %The Jacobian is constant over the element, discard the gaussian points.
0010 
0011 J=zeros(3,3);
0012 sqrt3=sqrt(3.0);
0013 sqrt6=sqrt(6);
0014 
0015 %grids coordinates
0016 x1=xyz_list(1,1);
0017 y1=xyz_list(1,2);
0018 z1=xyz_list(1,3);
0019 
0020 x2=xyz_list(2,1);
0021 y2=xyz_list(2,2);
0022 z2=xyz_list(2,3);
0023 
0024 x3=xyz_list(3,1);
0025 y3=xyz_list(3,2);
0026 z3=xyz_list(3,3);
0027 
0028 x4=xyz_list(4,1);
0029 y4=xyz_list(4,2);
0030 z4=xyz_list(4,3);
0031 
0032 %Computation of the jacobian on this gauss point
0033 J(1,1)=1.0/2.0* (-x1+x2);
0034 J(2,1)=1.0/2.0* (-y1+y2);
0035 J(3,1)=1.0/2.0* (-z1+z2);
0036 
0037 J(1,2)=1.0/(2*sqrt3)* (-x1+x2+2*x3);
0038 J(2,2)=1.0/(2*sqrt3)* (-y1+y2+2*y3);
0039 J(3,2)=1.0/(2*sqrt3)* (-z1+z2+2*z3);
0040 
0041 J(1,3)=1.0/(2*sqrt6)* (-x1-x2-x3+x4);
0042 J(2,3)=1.0/(2*sqrt6)* (-y1-y2-y3+y4);
0043 J(3,3)=1.0/(2*sqrt6)* (-z1-z2-z3+z4);

Generated on Sun 29-Mar-2009 20:22:55 by m2html © 2003