GetStrainRateStokes

PURPOSE ^

GETSTRAINRATESTOKES - return the strain rate for pentaelem with stokes' disgnostic

SYNOPSIS ^

function epsilon=GetStrainRateStokes(pentaelem,vxvyvz_list,xyz_list,gauss_coord)

DESCRIPTION ^

GETSTRAINRATESTOKES - return the strain rate for pentaelem with stokes' disgnostic
 
   This function return the strain rate needed for Glen's flow law for the element pentaelem
   The strain rate epsilon is calculated thanks to the formula epsilon=B*U
   Where epsilon=[exx eyy ezz exy exz eyz] and U=[u1;v1;..;u6;v6]
   The matrix B is the same as in K=tBDBprime

   Usage:
      epsilon=GetStrainRateStokes(pentaelem,vxvyvz_list,xyz_list,gauss_coord)
 
   See also GETSTRAINRATE, GETBSTOKES

CROSS-REFERENCE INFORMATION ^

This function calls: This function is called by:

SOURCE CODE ^

0001 function epsilon=GetStrainRateStokes(pentaelem,vxvyvz_list,xyz_list,gauss_coord)
0002 %GETSTRAINRATESTOKES - return the strain rate for pentaelem with stokes' disgnostic
0003 %
0004 %   This function return the strain rate needed for Glen's flow law for the element pentaelem
0005 %   The strain rate epsilon is calculated thanks to the formula epsilon=B*U
0006 %   Where epsilon=[exx eyy ezz exy exz eyz] and U=[u1;v1;..;u6;v6]
0007 %   The matrix B is the same as in K=tBDBprime
0008 %
0009 %   Usage:
0010 %      epsilon=GetStrainRateStokes(pentaelem,vxvyvz_list,xyz_list,gauss_coord)
0011 %
0012 %   See also GETSTRAINRATE, GETBSTOKES
0013 
0014 B=GetBStokes(pentaelem,xyz_list,gauss_coord);
0015 B_reduced=B(1:6,[1:3 5:7 9:11 13:15 17:19 21:23]);
0016 
0017 vxvyvz_list_transp=vxvyvz_list';
0018 U=[vxvyvz_list_transp(:,1);vxvyvz_list_transp(:,2);vxvyvz_list_transp(:,3);vxvyvz_list_transp(:,4);vxvyvz_list_transp(:,5);vxvyvz_list_transp(:,6)];
0019 
0020 epsilon=B_reduced*U; %for the six principal grids

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