GetStrainRate

PURPOSE ^

GETSTRAINRATE - return the strain rate of a pentaelem

SYNOPSIS ^

function epsilon=GetStrainRate(pentaelem,vxvy_list,xyz_list,gauss_coord)

DESCRIPTION ^

GETSTRAINRATE - return the strain rate of a pentaelem

   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 exy exz eyz] and U=[u1;v1;..;u6;v6]
   The matrix B is the same as in K=tBDBprime

   Usage:
      epsilon=GetStrainRate(pentaelem,vxvy_list,xyz_list,gauss_coord)

   See also GETSTRAINRATESTOKES, GETB_HORIZ

CROSS-REFERENCE INFORMATION ^

This function calls: This function is called by:

SOURCE CODE ^

0001 function epsilon=GetStrainRate(pentaelem,vxvy_list,xyz_list,gauss_coord)
0002 %GETSTRAINRATE - return the strain rate of a pentaelem
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 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=GetStrainRate(pentaelem,vxvy_list,xyz_list,gauss_coord)
0011 %
0012 %   See also GETSTRAINRATESTOKES, GETB_HORIZ
0013 
0014 B=GetB_horiz(pentaelem,xyz_list,gauss_coord);
0015 vxvy_list_transp=vxvy_list';
0016 U=[vxvy_list_transp(1:2,1);vxvy_list_transp(1:2,2);vxvy_list_transp(1:2,3);vxvy_list_transp(1:2,4);vxvy_list_transp(1:2,5);vxvy_list_transp(1:2,6)];
0017 
0018 epsilon=B*U;

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