


GETBPRIME_PROG - compute the B matrix for a prognostic
Compute B matrix. B=[B1 B2 B3 ] where Bi is of size numgrids*NDOF
For grid i, Bi can be expressed in the basic coordinate system by:
Bi_basic=[dh/dx]
[dh/dy]
where h is the interpolation function for grid i.
Usage:
B=GetBprime_prog(triaelem,xyz_list,gauss_coord)
See also GETB_PROG, CREATEKMATRIX

0001 function B=GetBprime_prog(triaelem,xyz_list,gauss_coord) 0002 %GETBPRIME_PROG - compute the B matrix for a prognostic 0003 % 0004 % Compute B matrix. B=[B1 B2 B3 ] where Bi is of size numgrids*NDOF 0005 % For grid i, Bi can be expressed in the basic coordinate system by: 0006 % Bi_basic=[dh/dx] 0007 % [dh/dy] 0008 % where h is the interpolation function for grid i. 0009 % 0010 % Usage: 0011 % B=GetBprime_prog(triaelem,xyz_list,gauss_coord) 0012 % 0013 % See also GETB_PROG, CREATEKMATRIX 0014 0015 num_grids=3; 0016 0017 dh1dh6_basic=GetNodalFunctionsDerivativesBasic(triaelem,xyz_list,gauss_coord); 0018 %Build B: 0019 B=zeros(2,num_grids); 0020 for i=1:num_grids, 0021 B(1,i)=dh1dh6_basic(1,i); 0022 B(2,i)=dh1dh6_basic(2,i); 0023 end