


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] and U=[u1;v1;..;u3;v3]
The matrix B is the same as in K=tBDBprime
Usage:
epsilon=GetStrainRate(pentaelem,vxvy_list,xyz_list,gauss_coord)
See also GETB

0001 function epsilon=GetStrainRate(triaelem,vxvy_list,xyz_list,gauss_l1l2l3) 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] and U=[u1;v1;..;u3;v3] 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 GETB 0013 0014 B=GetB(triaelem,xyz_list,gauss_l1l2l3); 0015 0016 vxvy_list_transp=vxvy_list'; 0017 U=[vxvy_list_transp(1:2,1);vxvy_list_transp(1:2,2);vxvy_list_transp(1:2,3)]; 0018 0019 epsilon=B*U;