Finite Element Global Digital Volume Correlation (FE-DVC)

Finite-element-based global DVC method (guarantee global kinematic compatibility and decrease noise by adding regularization penalties).

This is a MATLAB 3D-volumetric finite-element-based DVC code I wrote to compare with our newly proposed ALDVC algorithm [1,2].

% =========================================

% Installation

FE_Global_DVC MATLAB code was tested on MATLAB versions later than R2018a. Please download and unzip the code to the MATLAB working path. Then, execute the mail file main_FE_Global_DVC.m.

% =========================================

% Code Manual (will be uploaded soon)

% =========================================

% Comments and references

[1] For full details, and to use this code, please cite our paper:

Yang, J., Hazlett, L., Landauer, A., Franck, C. Augmented Lagrangian Digital Volume Correlation. Experimental Mechanics, 2020 (https://link.springer.com/article/10.1007/s11340-020-00607-3).

Or request full text at:

https://www.researchgate.net/publication/343676441_Augmented_Lagrangian_Digital_Volume_Correlation_ALDVC

[2] Welcome to compare this code with our proposed ALDVC code, which combines both advantages of Local Subvolume DVC (fast, compute in parallel) and Global DVC (guarantee kinematic compatibility, see https://www.mathworks.com/matlabcentral/fileexchange/77019-augmented-lagrangian-digital-volume-correlation-aldvc.

[3] % ===================================

% ****** ATTENTION ******

% The "x,y,z" or "1-,2-,3-" coordinates in the ALDVC code always correspond to the 1st, 2nd and 3rd indices of Matlab workspace variable. For example, p_meas(:,1) and p_meas(:,2) are the x- & y-coordinates of scattered points.

%

% This is a little different from some MATLAB image processing functions. % For example, if a 3D image has size MxNxL, in this code, we always have the image size_x=M, size_y=N, size_z=L. If you use some Matlab computer vision/image post-processing function, for example, 'imagesc3D', or 'imshow3D', or 'surf', it will reads size_x=N, size_y=M, size_z=L.

%

% Please pay attention to this difference.

[4] I also attach my EASF webinar to introduce AL-DIC/DVC algorithm and review local and global DIC/DVC methods:

(Youtube) https://www.youtube.com/watch?v=-t61WrVagZ4

(Bilibili) https://www.bilibili.com/video/BV1ff4y1B71L/

% =========================================

Besides 3D-volumetric DVC, our new 2D-DIC code "AL-DIC" (augmented Lagrangian Digital Image Correlation) and FE-Global-DIC code to track 2D deformations are also available at

https://www.mathworks.com/matlabcentral/fileexchange/70499-augmented-lagrangian-digital-image-correlation-and-tracking

and

https://www.mathworks.com/matlabcentral/fileexchange/82873-2d-finite-element-global-digital-image-correlation-fe-dic

% =========================================

Contact & support:

I appreciate your comments and ratings to help me further improve this code. If you have other questions and suggestions, feel free to email me: [e-mail]

% ==================================================

References

[1] regularizeNd. https://www.mathworks.com/matlabcentral/fileexchange/61436-regularizend

[2] https://www.mathworks.com/matlabcentral/fileexchange/24049-streamcolor