0001 function J=GetJacobian(tetraelem,xyz_list,gauss_coord)
0002
0003
0004
0005
0006
0007
0008
0009
0010
0011 J=zeros(3,3);
0012 sqrt3=sqrt(3.0);
0013 sqrt6=sqrt(6);
0014
0015
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
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);