KIII_2D_v2.m This code decompose the Stress intesity factors from strain maps directly without the need for integration Start with generating strain data using the calibration code and then use the output for the main function The code is self contained and does not need extra functions
There are two input desks for 2D/3D DIC data and HR-EBSD data
% for this function to work properly, spacing between potins in x and y % should be the same and the crack should be at the centre (this can be % done inside this fucntion also)
% this functions accept data from HR-EBSD and sttero-DIC % The crack needs to be on the centre exactly for the code to % work and give good results. This code already include the assumption of % sigma 33 free == 0 % stress in Pa, E in Pa, distance in m % use variable 'Stiffness' if you are using anistropic material or E and nu % for istroropic material % need some unit clbartion specially for mm
% Another option is to input the strain components as a vector matrix with % 9 columns the first three columns are the x, y and z coordinate in meters. z % cooridnate can be a zero column. the 4th to the 9th column are the strain % components arranged as % Maps = [X(:) Y(:) Z(:) E11(:) E12(:) E13(:) E21(:) E22(:) E23(:) E31(:) E32(:) E33(:)];
% if the map is a 2D strain map then zero all out of the plane components
% the material paramters as MatProp.E for Young's Modulus and MatProp.nu % for Possions ratio. or as a stifness matrix all in Pa % xEBSD will assume and solve for plane strain conditions with sigma33 == 0 + no % volumetric change.