


GETBPRIME - compute the Bprime matrix for a diagnostic solution
Compute B' matrix. B'=[B1' B2' B3'] 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 ]
where h is the interpolation function for grid i.
Usage:
Bprime=GetBprime(triaelem,xyz_list,gauss_l1l2l3)
See also GETB, CREATEKMATRIX

0001 function Bprime=GetBprime(triaelem,xyz_list,gauss_l1l2l3) 0002 %GETBPRIME - compute the Bprime matrix for a diagnostic solution 0003 % 0004 % Compute B' matrix. B'=[B1' B2' B3'] 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 % where h is the interpolation function for grid i. 0010 % 0011 % Usage: 0012 % Bprime=GetBprime(triaelem,xyz_list,gauss_l1l2l3) 0013 % 0014 % See also GETB, CREATEKMATRIX 0015 0016 num_grids=3; 0017 NDOF2=2; 0018 0019 dh1dh2dh3_basic=GetNodalFunctionsDerivativesBasic(triaelem,xyz_list,gauss_l1l2l3); 0020 0021 %Build Bprime: 0022 Bprime=zeros(3,NDOF2*num_grids); 0023 for i=1:num_grids, 0024 Bprime(1,NDOF2*(i-1)+1)=2*dh1dh2dh3_basic(1,i); 0025 Bprime(1,NDOF2*(i-1)+2)= dh1dh2dh3_basic(2,i); 0026 Bprime(2,NDOF2*(i-1)+1)= dh1dh2dh3_basic(1,i); 0027 Bprime(2,NDOF2*(i-1)+2)=2*dh1dh2dh3_basic(2,i); 0028 Bprime(3,NDOF2*(i-1)+1)=dh1dh2dh3_basic(2,i); 0029 Bprime(3,NDOF2*(i-1)+2)=dh1dh2dh3_basic(1,i); 0030 end