%************************************************
% Define elastic stiffness tensor (GPa)
%************************************************
%
% Reference Elastic constants
% Chedda T., Mookherjee, M., Moreira dos Santos, A.F., Molaison,J.,
% Chantel, J., Manthilake, G., and Mainprice, D. (2013) AGU abstract
%
% Reference Crystal Structure
% Hazen,R.M. and Burnham,C.W.(1973) The crystal structures of one-layer
% phlogopite and annite, Am. Mineral. 58,889-900.
% C2/m 1M KMg3(AiSi3)O10(OH)2 Z=2 d=2.872 g/cm3
%
% Define density (g/cm3)
rho= 2.872;
% N.B. Crystal structure is monoclinic
% z chosen normal to (001) planes (z||c*)
% Define Cartesian Tensor crystal symmetry & frame
cs_Tensor = crystalSymmetry('2/m',[  5.308  9.190 10.155],...
 [  90.00 100.08  90.00]*degree,'x||a','z||c*',...
 'mineral','Phlogopite 2014');
%
% elastic Cij stiffness tensor (GPa) as matrix M
 M =.... 
 [[  181.00   48.00   12.00    0.00  -16.00    0.00];...
 [    48.00  185.00   12.00    0.00   -5.00    0.00];...
 [    12.00   12.00   62.00    0.00   -1.00    0.00];...
 [     0.00    0.00    0.00   14.00    0.00   -6.00];...
 [   -16.00   -5.00   -1.00    0.00   20.00    0.00];...
 [     0.00    0.00    0.00   -6.00    0.00   68.00]];
%
% M as stiffness tensor C with MTEX tensor command
C = tensor(M,cs_Tensor);
%
