GetNodalFunctionsDerivativesParams

PURPOSE ^

GETNODALFUNCTIONSDERIVATIVESPARAMS - return the values of nodal functions derivatives on the reference element

SYNOPSIS ^

function dl1dl6=GetNodalFunctionsDerivativesParams(pentaelem,gauss_coord)

DESCRIPTION ^

GETNODALFUNCTIONSDERIVATIVESPARAMS - return the values of nodal functions derivatives on the reference element
    
   Usage:
      dl1dl6=GetNodalFunctionsDerivativesParams(pentaelem,gauss_coord)

   See also GETNODALFUNCTIONS, GETNODALFUNCTIONSDERIVATIVESBASIC, GETNODALFUNCTIONSDERIVATIVESPARAMSSTOKES

CROSS-REFERENCE INFORMATION ^

This function calls: This function is called by:

SOURCE CODE ^

0001 function dl1dl6=GetNodalFunctionsDerivativesParams(pentaelem,gauss_coord)
0002 %GETNODALFUNCTIONSDERIVATIVESPARAMS - return the values of nodal functions derivatives on the reference element
0003 %
0004 %   Usage:
0005 %      dl1dl6=GetNodalFunctionsDerivativesParams(pentaelem,gauss_coord)
0006 %
0007 %   See also GETNODALFUNCTIONS, GETNODALFUNCTIONSDERIVATIVESBASIC, GETNODALFUNCTIONSDERIVATIVESPARAMSSTOKES
0008 
0009 sqrt3=sqrt(3.0);
0010 
0011 %Gauss point coordinate in reference element
0012 r=gauss_coord(1,2) - gauss_coord(1,1);
0013 s= - 3.0/sqrt3 *(gauss_coord(1,1)+gauss_coord(1,2) - 2.0/3.0);
0014 zeta=gauss_coord(1,4);
0015 
0016 num_grids=6;
0017 dl1dl6=zeros(3,num_grids);
0018 
0019 %First nodal function:
0020 dl1dl6(1,1)=-1.0/2.0*(1-zeta)/2;
0021 dl1dl6(2,1)=-sqrt3/6*(1-zeta)/2;
0022 dl1dl6(3,1)=-1.0/2.0*(-1.0/2.0*r-sqrt3/6*s+1/3);
0023 
0024 %Second nodal function:
0025 dl1dl6(1,2)=1.0/2.0*(1-zeta)/2;
0026 dl1dl6(2,2)=-sqrt3/6*(1-zeta)/2;
0027 dl1dl6(3,2)=-1.0/2.0*(1.0/2.0*r-sqrt3/6*s+1/3);
0028 
0029 %Third nodal function:
0030 dl1dl6(1,3)=0;
0031 dl1dl6(2,3)=sqrt3/3.0*(1-zeta)/2;
0032 dl1dl6(3,3)=-1.0/2.0*(sqrt3/3.0*s+1/3);
0033 
0034 %Fourth nodal function:
0035 dl1dl6(1,4)=-1.0/2.0*(1+zeta)/2;
0036 dl1dl6(2,4)=-sqrt3/6*(1+zeta)/2;
0037 dl1dl6(3,4)=1.0/2.0*(-1.0/2.0*r-sqrt3/6*s+1/3);
0038 
0039 %Fifth nodal function:
0040 dl1dl6(1,5)=1.0/2.0*(1+zeta)/2;
0041 dl1dl6(2,5)=-sqrt3/6*(1+zeta)/2;
0042 dl1dl6(3,5)=1.0/2.0*(1.0/2.0*r-sqrt3/6*s+1/3);
0043 
0044 %Sixth nodal function:
0045 dl1dl6(1,6)=0;
0046 dl1dl6(2,6)=sqrt3/3.0*(1+zeta)/2;
0047 dl1dl6(3,6)=1.0/2.0*(sqrt3/3.0*s+1/3);
0048

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