


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

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