GetBprime_horiz

PURPOSE ^

GETBPRIME_HORIZ - compute the Bprime matrix for horizontal velocities in Pattyn's model

SYNOPSIS ^

function Bprime=GetBprime_horiz(pentaelem,xyz_list,gauss_coord)

DESCRIPTION ^

GETBPRIME_HORIZ - compute the Bprime matrix for horizontal velocities in Pattyn's model

   Compute B'  matrix. B'=[B1' B2' B3' B4' B5' B6'] where Bi' is of size 3*NDOF2. 
   For grid i, Bi' can be expressed in the basic coordinate system by:
   Bi_basic'=[ 2*dh/dx    dh/dy ]
             [   dh/dx   2*dh/dy]
             [   dh/dy   dh/dx  ]
             [   dh/dz      0   ]
             [   0       dh/dz  ]
   where h is the interpolation function for grid i.

   Usage:
      Bprime=GetBprime_horiz(pentaelem,xyz_list,gauss_coord)

   See also GETB_HORIZ, CREATEKMATRIX

CROSS-REFERENCE INFORMATION ^

This function calls: This function is called by:

SOURCE CODE ^

0001 function Bprime=GetBprime_horiz(pentaelem,xyz_list,gauss_coord)
0002 %GETBPRIME_HORIZ - compute the Bprime matrix for horizontal velocities in Pattyn's model
0003 %
0004 %   Compute B'  matrix. B'=[B1' B2' B3' B4' B5' B6'] where Bi' is of size 3*NDOF2.
0005 %   For grid i, Bi' can be expressed in the basic coordinate system by:
0006 %   Bi_basic'=[ 2*dh/dx    dh/dy ]
0007 %             [   dh/dx   2*dh/dy]
0008 %             [   dh/dy   dh/dx  ]
0009 %             [   dh/dz      0   ]
0010 %             [   0       dh/dz  ]
0011 %   where h is the interpolation function for grid i.
0012 %
0013 %   Usage:
0014 %      Bprime=GetBprime_horiz(pentaelem,xyz_list,gauss_coord)
0015 %
0016 %   See also GETB_HORIZ, CREATEKMATRIX
0017 
0018     num_grids=6;    
0019     NDOF2=2;
0020 
0021     dh1dh6_basic=GetNodalFunctionsDerivativesBasic(pentaelem,xyz_list,gauss_coord);
0022 
0023     %Build Bprime:
0024     Bprime=zeros(3,NDOF2*num_grids);
0025     for i=1:num_grids,
0026         Bprime(1,NDOF2*(i-1)+1)=2*dh1dh6_basic(1,i);
0027         Bprime(1,NDOF2*(i-1)+2)=  dh1dh6_basic(2,i);
0028         Bprime(2,NDOF2*(i-1)+1)=  dh1dh6_basic(1,i);
0029         Bprime(2,NDOF2*(i-1)+2)=2*dh1dh6_basic(2,i);
0030         Bprime(3,NDOF2*(i-1)+1)=dh1dh6_basic(2,i);
0031         Bprime(3,NDOF2*(i-1)+2)=dh1dh6_basic(1,i);        
0032         Bprime(4,NDOF2*(i-1)+1)=dh1dh6_basic(3,i);
0033         Bprime(4,NDOF2*(i-1)+2)=0;
0034         Bprime(5,NDOF2*(i-1)+1)=0;
0035         Bprime(5,NDOF2*(i-1)+2)=dh1dh6_basic(3,i);
0036 
0037     end

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