


GETNODALFUNCTIONSDERIVATIVESBASICSTOKES - return the values of nodal functions derivatives on the current element
Usage:
dh1dh7_basic=GetNodalFunctionsDerivativesBasicStokes(pentaelem,xyz_list,gauss_coord)
See also GETNODALFUNCTIONSSTOKES, GETNODALFUNCTIONSDERIVATIVESPARAMSSTOKES, GETNODALFUNCTIONSDERIVATIVESBASIC

0001 function dh1dh7_basic=GetNodalFunctionsDerivativesBasicStokes(pentaelem,xyz_list,gauss_coord) 0002 %GETNODALFUNCTIONSDERIVATIVESBASICSTOKES - return the values of nodal functions derivatives on the current element 0003 % 0004 % Usage: 0005 % dh1dh7_basic=GetNodalFunctionsDerivativesBasicStokes(pentaelem,xyz_list,gauss_coord) 0006 % 0007 % See also GETNODALFUNCTIONSSTOKES, GETNODALFUNCTIONSDERIVATIVESPARAMSSTOKES, GETNODALFUNCTIONSDERIVATIVESBASIC 0008 0009 num_grids=7; 0010 0011 %Get derivative values with respect to parametric coordinate system: 0012 dh1dh7_param=GetNodalFunctionsDerivativesParamsStokes(pentaelem,gauss_coord); 0013 0014 %Get Jacobian invert: 0015 Jinv=GetJacobianInvert(pentaelem,xyz_list, gauss_coord); 0016 0017 %Build dh1dh7_basic: 0018 % [dhi/dx]= Jinv'*[dhi/dr] 0019 % [dhi/dy] [dhi/ds] 0020 % [dhi/dz] [dhi/dzeta] 0021 0022 dh1dh7_basic=zeros(3,num_grids); 0023 0024 for i=1:num_grids, 0025 dh1dh7_basic(1,i)=Jinv(1,1)*dh1dh7_param(1,i)+Jinv(2,1)*dh1dh7_param(2,i)+Jinv(3,1)*dh1dh7_param(3,i); 0026 dh1dh7_basic(2,i)=Jinv(1,2)*dh1dh7_param(1,i)+Jinv(2,2)*dh1dh7_param(2,i)+Jinv(3,2)*dh1dh7_param(3,i); 0027 dh1dh7_basic(3,i)=Jinv(1,3)*dh1dh7_param(1,i)+Jinv(2,3)*dh1dh7_param(2,i)+Jinv(3,3)*dh1dh7_param(3,i); 0028 end