


GETB_PROG - compute Bprime 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=[h]
[h]
where h is the interpolation function for grid i.
Usage:
B=GetB_prog(triaelem,xyz_list,gauss_coord)
See also GETBPRIME_PROG, CREATEKMATRIX

0001 function B=GetB_prog(triaelem,xyz_list,gauss_coord) 0002 %GETB_PROG - compute Bprime 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=[h] 0007 % [h] 0008 % where h is the interpolation function for grid i. 0009 % 0010 % Usage: 0011 % B=GetB_prog(triaelem,xyz_list,gauss_coord) 0012 % 0013 % See also GETBPRIME_PROG, CREATEKMATRIX 0014 0015 num_grids=3; 0016 0017 %Get l1l2l3 in basic coordinate system: 0018 l1l2l3=GetNodalFunctions(triaelem,gauss_coord); 0019 0020 B=zeros(2,num_grids); 0021 for i=1:num_grids, 0022 B(1,i)=l1l2l3(i); 0023 B(2,i)=l1l2l3(i); 0024 end