GetBprime

PURPOSE ^

GETBPRIME - compute the Bprime matrix for a diagnostic solution

SYNOPSIS ^

function Bprime=GetBprime(triaelem,xyz_list,gauss_l1l2l3)

DESCRIPTION ^

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

CROSS-REFERENCE INFORMATION ^

This function calls: This function is called by:

SOURCE CODE ^

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

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