


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

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