HCP施密特因子Schmid-factor-Mtalab计算.docVIP

clearclcratio = 1.624;% for Mgcrss_basal = 2.1;% CRSS for basal a in Mgcrss_pyr = 40;% CRSS for c+a in Mgcrss_ext = 11.5;% CRSS for extension twin in Mgcrss_comp = 125;% CRSS for compression twin in Mg%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%basal = [0,0,0,1];%basal planepyr1 = [-1,0,1,-1];%pyramidal I planepyr2 = [-1,-1,2,-2];%pyramidal II planeext = [1,0,-1,2];%extension twincomp = [-1,0,1,-1];%compression twinbasal_dir = [1,1,-2,0];%basal a Burgers vectorpyr1_dir = [1,1,-2,-3];%pyramidal c+a Burgers vectorpyr2_dir = [1,1,-2,-3];%pyramida