A single camera three-dimensional digital image correlation system for the study of adiabatic shear bands
T. G. White J. R. W. Patten K. -H. Wan A. D. Pullen D. J. Chapman D. E. Eakins

A single camera 3D DIC system for the study of adiabatic shear bandsT. G. White et al.





Institute of Shock Physics, Imperial College London, London, UK 


We describe the capability of a high resolution three-dimensional digital image correlation (DIC) system specifically designed for high strain-rate experiments. Utilising open-source camera calibration and two-dimensional digital image correlation tools within the MATLAB framework, a single camera 3-D DIC system with sub-micron displacement resolution is demonstrated. The system has a displacement accuracy of up to 200 times the optical spatial resolution, matching that achievable with commercial systems. The surface strain calculations are benchmarked against commercially available software before being deployed on quasi-static tests showcasing the ability to detect both in- and out-of-plane motion. Finally, a high strain-rate (1.2×10^3 s^-1) test was performed on a top-hat sample compressed in a split-Hopkinson pressure bar in order to highlight the inherent camera synchronisation and ability to resolve the adiabatic shear band phenomenon. 




A single camera three-dimensional digital image correlation system for the study of adiabatic shear bands
    T. G. White, J. R. W. Patten, K. Wan, A. D. Pullen, D. J. Chapman and D. E. Eakins
    December 30, 2023
=========================================================================================================





§ INTRODUCTION



Full-field deformation measurements are an important component in the analysis of material behaviour when subjected to both quasi-static and high-rate loading.  Unlike traditional methods such as load cells, which provide global information, or strain gauges, which yield information averaged over a limited portion of the sample, optical techniques are a non-contact method able to provide point-wise information over a large area of the sample. This becomes particularly important when investigating inhomogeneous behaviour or systems in hostile environments that necessitate a non-contact measurement. 

Developed in the 1980s <cit.>, Digital Image Correlation (DIC) is a non-contact optical method capable of measuring the full-field strain distribution over a sample surface. The relatively lenient requirements necessary to perform DIC has lead to widespread use and popularity beyond competing optical techniques such as Moiré interferometry <cit.> or Holography <cit.>. Two-dimensional (2-D) DIC utilises a single camera to measure in-plane sample deformation through cross-correlation of sequential greyscale images and is capable of obtaining displacement vectors with up to 0.01 pixel precision <cit.>. Since its conception DIC has become an essential technique in experimental/solid mechanics as well as having uses in material characterization, architecture, aerospace and biology <cit.>. 



Two-dimensional DIC is limited to measuring the in-plane deformation of planar samples, hence requiring the camera to be located normal to the target surface. However, sample deformation is rarely limited to in-plane motion and more general target geometries require three-dimensional (3-D) position and displacement measurements. In this work this is achieved through 3-D DIC, that is stereoscopic imaging of the sample using two calibrated cameras. Combining knowledge of the relative camera angles with standard 2-D DIC reveals the out-of-plane deformation <cit.>. 

Many commercial systems for 3-D DIC now exist (VIC-3D <cit.>, ARAMIS <cit.>, DaVIS <cit.>), however, both the hardware and software can be typically prohibitively expensive for smaller research groups. In addition, these closed source codes prevent modification or improvement by the user and can lead to uncertainty over the precise details of the algorithms implemented. 

Performing 3-D DIC at high strain-rate introduces further complications and expense as this typically necessitates the use of two synchronized high speed cameras. In such systems the two greatest costs are the high frame-rate cameras and the bespoke software required for stereo camera calibration and digital image correlation, while the electronic synchronization of the two cameras provides the largest source of temporal error. To this end the development of a single camera 3D-DIC system, which eliminates issues surrounding high rate camera synchronization and offers improved stability, has received considerable attention. However, these attempts place stringent requirements on target geometry <cit.>, require precise camera and optical set-ups <cit.> or are analyzed with commercial software <cit.>. Here we aim to demonstrate a low-cost 3-D DIC system capable of achieving high rate, high spatial resolution, three-dimensional deformation measurements.

In the following study we develop a single camera 3-D DIC system which alleviates synchronization errors and is based on freely available open source codes. The camera calibration was carried out using the Caltech Calibration Toolbox for MATLAB <cit.> while the 2-D DIC is achieved using Ncorr developed at Georgia Institute of Technology <cit.>. The Ncorr system has previously demonstrated excellent agreement with commercial software <cit.>. 

The aim of this work is to assess the applicability of using a single camera together with open-source software to achieve three-dimensional high-rate, high spatial and displacement resolution DIC, taking advantage of the improved temporal synchronisation and stability associated with single camera systems, the flexibility and transparency associated with using open-source codes and at a fraction of the cost of a dual-camera commercial system. 

This system has been tailored for a field of view of a few millimeters as this allows the study of adiabatic shear bands (ASBs), a widely encountered phenomenon in engineering alloys subjected to high strain-rate loading <cit.>. They form on the sub-millimeter scale and are of particular interest to the aerospace industry where shear bands often occur when a hard object impacts the rotating fan blades at high speed; the velocity of impact leads to low thermal conduction and an increase in temperature, thermal softening and ultimately in catastrophic failure from brittle-like fracture  <cit.>. The study of ASBs is one example of a situation that requires high-rate 3-D DIC, however the approach discussed here is broadly applicable, and can be tailored to different applications by simple exchange of camera lens and optical components. 

The need for 3-D DIC in the study of ASBs is not immediately obvious as they are typically a 2D phenomenon when viewed on the surface of the sample. However, with a 2-D system the camera must be placed orthogonal to the sample surface which is assumed to move only in-plane, this constraint is relaxed in the 3-D setup as the surface normal is explicitly measured. A 3-D system thus allows measurements to be taken despite sample motion (e.g. tilt, bulging, necking) and when non-planar sample geometries are used. Furthermore, ASBs are by definition a localization phenomenon and therefore marks the material deformation in a non-homogeneous manner with regions of large plastic deformation <cit.>. Measurement of the out-of-plane displacement allows the full finite-strain Lagrange strain tensor to be calculated <cit.>. Retaining all terms in the Lagrange strain tensor ensures the large finite strains present in ASBs are appropriately described.  





§ EXPERIMENTAL SET-UP





The general approach uses reflective optics to image a sample surface from two viewing-angles onto a single camera. The set-up consists of two planar 25 mm square mirrors along with a single knife-edge right-angled prism mirror; these were used to create two views of the target on a single CCD, see Figure <ref>. Both the prism mirror and side mirrors are separately mounted such that they can be moved independently to adjust the image on the camera. As with most commercial set-ups the system is relatively insensitive to the exact angles used due to the robustness of the camera calibration method, the angles do not need to be known accurately a-priori. However, the angle subtended by the two mirrors and the sample (marked θ in Figure <ref>) should remain between 20^∘-30^∘ to achieve a good balance between out-of-plane precision and image correlation. A larger angle allows the out-of-plane motion to be determined with greater accuracy while conversely making the image correlation between the two images more difficult due to increased image disparity <cit.>. The distances and angles used in this work are given in Figure <ref>, as are images showing a typical field of view as focus is shifted from the knife-edge to the target. Future work could involve extending the ideas presented here to more than two images, and hence viewing angles, helping to constrain the problem further and improve accuracy.

In order to demonstrate the ability to measure full-field displacement of both in- and out-of-plane motion from a single camera the system was set-up to observe a flat speckled surface. The speckle pattern was applied using an Iwata HP-B Plus airbrush with a 200 μm nozzle using Daler Rowney FW Acrylic Artist's Ink, Black 028. This ink was found to adhere best to the sample surface. The pressure was varied to alter the size of the droplets; a pressure of 1 bar was found to give an optimal speckle. In this case an optimal speckle is defined as being without large featureless areas which reduce the achievable spatial resolution and without fine marks below the optical resolution of the set-up. A nozzle to sample distance of 23 cm was used. A good review of optimal speckle patterns for DIC can be found in reference <cit.>. A typical speckle pattern can be seen in the insert of Figure <ref>. 

This target was mounted on a Thorlabs PT3-Z8 translation stage providing XYZ motion while the motorised Z825B actuators provide up to 25 mm travel with a minimum resolution of 29 nm. A camera (either a digital SLR or high-speed camera) was attached to an Infinity K2 DistaMax long-distance microscope lens (standard configuration, with 4× magnification provided by two NTX tube amplifiers) to allow ultra-high resolution images to be taken. The entire experiment was performed on a vibrationally-damped optical table to isolate the set-up from any external vibrations. The high frame-rate Phantom is capable of taking 6600 frames a second utilising all 800x600 pixels and up to 50000 frames per second at a reduced pixel resolution of 240x184. Table <ref> lists the parameters of the three different set-ups evaluated in this work. 



§ RECONSTRUCTING THE 3-D WORLD



The majority of this work revolves around reconstruction of the 3-D world from a pair of stereoscopic images. The general method is outlined in Figure <ref>, with steps 1-3 representing calibration of the system, steps 4-5 digital image correlation between pairs of images, step 6 combines these results with the calibration to produce a 3-D surface and displacement map which are then both used to produce a surface strain map. This work utilises the pinhole camera model for the projective mapping from 3-D world coordinates, denoted 𝐏_𝐰=[X_w Y_w Z_w 1],  to the two-dimensional camera coordinates, denoted 𝐂_1=[u_1 v_1 1]. Mathematically this is expressed as the matrix multiplication,


    𝐂_1=
    [ u_1; v_1;   1 ] 
    =𝐊_1·[ 𝐑_1 | 𝐓_1 ]·𝐏_𝐰 ,
 where K_1 is the intrinsic matrix of the camera given by,


    𝐊_1= [ f_x   s x_0;   0 f_y y_0;   0   0   1 ].
 



The intrinsic matrix contains the coordinates of the principal point of the camera, (x_0, y_0), the horizontal and vertical focal lengths of the camera (f_x. f_y) given in units of pixel dimensions, and s=α_cf_x which allows for non-rectangular pixels through the angle α_c. The terms 𝐑_1 and 𝐓_1 together define the 3×4 extrinsic matrix of camera 1 and relates the position of the camera to a particular world coordinate frame, through a rotation followed by a rigid body translation. For a stereo system, i.e. a system where two cameras are used to resolve a single scene, an equivalent equation can be written for the second camera,


    𝐂_2= 𝐊_2·[ 𝐑_2 | 𝐓_2 ]·𝐏_𝐰.
 

Throughout this work the reference frame of camera 1 is used as the world reference frame, hence, 𝐑_1=𝐈, 𝐓_1=0, 𝐑_2=𝐑 and 𝐓_2=𝐓. Therefore equations (<ref>) and (<ref>) together represent a system of four simultaneous equations. Hence, it is possible with knowledge of the position of a single point on each camera, u_1,v_1,u_2,v_2 the equations can be solved to obtain the position of that point in the world reference frame, 𝐏_𝐰. The identification of equivalent points within each image is achieved through standard 2-D DIC.  

If the imaging system contains distortion then the reconstructed 3-D world, represented by the points 𝐏_𝐰, will also suffer from distortion. However, the camera calibration toolbox is naturally able to identify distortion in the imaging system through the chequerboard calibration method <cit.>. Radial distortion is parameterized through the two parameters κ_1 and κ_2 while tangential distortion is included through τ_1 and τ_2. The distorted world coordinates, 𝐏_𝐰, can then be related to the undistorted coordinates, 𝐏'_𝐰, through Equations (<ref>) and (<ref>), where x=X'_w/Z'_w and y=Y'_w/Z'_w are the reduced undistorted coordinates, r^2=x^2+y^2 and 𝐓_𝐝 the translational component of distortion.
 

    [ X_w/Z_w; Y_w/Z_w ] = [ x; y ] (1+κ_1r^2+κ_2r^4) + 𝐓_𝐝 ,
 
 

    𝐓_𝐝=[ 2τ_1xy + τ_2(r^2 2x^2); 2τ_2xy + τ_1(r^2 2y^2) ] .
 

Steps 1-3 represent the process of camera resectioning, that is the technique by which both the intrinsic (𝐊 - focal length, pixel properties and distortion) as well as the extrinsic (𝐓 - translational and 𝐑 - rotational displacements) parameters of two or more cameras can be found. It is a crucial component of any stereoscopic imaging technique and this work utilises the Caltech Calibration Toolbox <cit.> as well as the inbuilt Stereo Calibration function in MATLAB <cit.>. Calibration in both methods is achieved using multiple images of a known, and well defined, two-dimensional pattern. Typically an asymmetrical chequerboard is used. Both methods used here utilise a non-linear minimisation on the re-projected chequerboard to find the unknown coefficients through a steepest gradient descent method <cit.>. The multiple chequerboard images enable the 3-D space to be mapped onto the 2-D images. Typically, the re-projection errors, i.e. the offset between a chequerboard corner and the re-projection of the corner using the pin-hole camera formalism are often found to be sub-pixel in all but the most distorted systems. The chequerboard calibration technique contains a single free parameter, the chequerboard width, which sets the scale of the image.



Steps 4-5 of the 3-D reconstruction is based on identifying the pixel position (𝐂 - u_1,v_1,u_2,v_2) of equivalent points in the set of two stereoscopic images, this is achieved through standard 2-D DIC techniques. Ncorr is an open-source MATLAB program for 2-D DIC which runs entirely inside the MATLAB environment; this enables easy compatibility with the camera calibration methods discussed above, but additionally utilises C++ and MEX to improve efficiency <cit.>. Ncorr uses the inverse compositional method which is fast, robust and accurate compared to more traditional Newton-Raphson techniques. During the 2-D DIC the image is divided into circular subsets which are used to correlate spatial regions between the two images, the larger this subset the greater precision with which displacement or position can be found but at the cost of spatial resolution, features smaller than the subset size are unable to be resolved. A circular subset of radius 21 pixels was found to work well for the current implementation, being the smallest subset possible without producing noisy artefacts in the data, and has been used in the displacement, low rate and high rate tests. The comparison with the DaVIS software was performed at a 31 subset pixel size. 



Step 6 involves reconstructing the 3-D world from the 2-D information. By using the known camera parameters (𝐑_1,𝐓_1,𝐊_1,𝐑_2,𝐓_2,𝐊_2) and equivalent pixel locations (𝐂_1,𝐂_2) to solve the resultant simultaneous equations given by (<ref>) and (<ref>) the world pixel locations, 𝐏_𝐖, are obtained. In order to reconstruct a 3-D scene at a single time the 2-D DIC needs to be run between the left and right stereoscopic images only. However, obtaining 3-D displacements is more involved as it is not possible to perform 2-D DIC between the stereoscopic images at time t_1 and separately at time t_0 and subtract the results as image correspondence is lost, i.e. pixels at two differing times may not be looking at the same point in the 3-D space. Instead 2-D DIC must be performed between the images in as shown by the lines in the dashed box in Figure <ref>. This enables the position of equivalent pixels in the each images at two differing times to be obtained and solves the sub-image correspondence problem, allowing the software to ascertain which parts of each image correspond to the same part in another <cit.>. This allows for subtraction of the world coordinates providing information on the 3-D displacement. A graphical user interface was developed in MATLAB which interfaced with the camera calibration toolbox and the Ncorr DIC system allowing the entire process to be carried out from within a single program. 



§ STRAIN CALCULATION



Step 7 involves calculation of the strain map from the surface and displacement maps. Calculation of the full field strains is performed over a local volume element, centred on each reconstructed point and containing those points which lie within a predefined radius. Judicial choice of this radius is required to optimise the balance between accuracy and smoothness. The displacements are rotated into a local coordinate system with the z axis locally normal to the surface. The displacement gradient tensor,


    𝐇=[ ∂ u_x/∂ x ∂ u_x/∂ y ∂ u_x/∂ z; ∂ u_y/∂ x ∂ u_y/∂ y ∂ u_y/∂ z; ∂ u_z/∂ x ∂ u_z/∂ y ∂ u_z/∂ z ]
 can then be calculated through linear regression. The deformation gradient tensor (𝐅=𝐇+𝐈), Cauchy strain tensor (𝐂=(𝐅^T𝐅) and Lagrange strain tensor (𝐄=1/2(𝐂-𝐈)) are defined in the usual way. It should be noted that the derivatives with respect to the surface normal direction, here defined as z, cannot be obtained from measurements of the surface displacement without further assumptions which leads to only the E_xx, E_yy, E_xy and E_yx components of the Lagrange strain tensor being defined, with the in-plane local directions x and y chosen in advance. 



§ DISPLACEMENT ACCURACY



An initial demonstration of the accuracy of any DIC system is to replicate known linear motions both in- and out-of-plane. For this purpose we drove a flat target with the speckle pattern shown in the inset of Fig. <ref> distances of up to 160 μm in the x-, y- and z-directions, with the z-direction corresponding to out-of-plane displacement. Using Thorlabs Z825BV actuators the speckle pattern was driven to a specified distance, always moving in the direction of positive travel to avoid backlash. At each distance the sample was held stationary and an image taken. The same speckle pattern was used for each of the tests. Figure <ref> shows the displacements measured with the high resolution Canon set-up and the read-out from the linear actuators; the results show exceptional agreement between the calculated and actual displacements with an RMS error of 0.5 μm. In order to test the performance of the reconstruction algorithm during high rate experiments the Canon EOS 600D was replaced with a Phantom V7.3 high speed camera using just 240x184 pixels of the CCD. At this resolution the camera is capable of running at 50000 frames per second, the results are shown in Figure <ref> and demonstrate a RMS error of 1.6 μm. The RMS error was found to be similar, i.e. within 20%, in each of the orthogonal directions.

However, it is not clear if the reduction in precision can be attributed to the reduced optical resolution, the applicability of the speckle pattern at this resolution or instability introduced through the cooling fan in the Phantom V7.3. To test for sources of vibration the displacement of subsequent images of a stationary target were measured and found to be 1.6 μm when taken 1 s apart and 1.5 μm when taken 0.2 ms apart, suggesting some component of the resolution is due to vibrational motion. This becomes more pronounced when running the Phantom with the full CCD area (800x600 pixels) at 6600 frames per second. In this case the displacement of subsequent images of a stationary target were found to be 1.5 μm when taken 1 s apart and 0.5 μm when taken 0.15 ms apart. This discrepancy suggests a source of vibration in the measurements with the Phantom 7.3. The lack of vibrational motion in the high speed data does not preclude the possibility that the source of error is from fluctuations in the target environment, the lack of such motion with the Canon is highly suggestive that the cause is internal to the Phantom. Indeed, the Phantom V7.3 uses a continuous cooling fan while newer models (e.g. v2511) contain a `Quiet Fans' mode which turns off the fans for a short period specifically for vibration-sensitive applications. 

It should be noted that for optimal performance the idealised speckle pattern should be modified to match the field of view and resolution <cit.>, however this work utilised the same pattern in each case, therefore it may be possible to improve these results further. However, the achieved accuracy is consistently around 200 times greater than the optical resolution, comparable to other DIC systems <cit.>, and as such we would expect this to give marginal improvements.





§ COMPARISON TO DAVIS SOFTWARE


In order to demonstrate the applicability of the software we have performed a comparison of our system to a commercially available code, in this case the DaVIS software by LaVision<cit.>. A comparison test was performed on a short length of cold-formed galvanised steel C-section with the dimensions given in Figure <ref>. A two camera system compatible with the DaVIS software was used and were calibrated using their respective techniques, in the case of the DaVIS software this was a 3D calibration board, while our system was calibrated using the chequerboard technique. The same images were fed into both codes to remove ambiguity around illumination, speckle quality, camera resolution or lens quality. Speckle was applied using standard black spray paint over a white base coating to give maximum contrast, see Figure <ref>. Compression was applied to the sample at a constant velocity using an Instron SATEC with a 600 kN loadcell, at 10 mm/min. 

The results of the compression test are given in Figure <ref>. Figure <ref> shows the displacement of the sample along the dashed-dotted line in Figure <ref>, results are shown for times corresponding to 1 minute snapshots from the beginning of the compression. Due to the lack of a fiducial the results from the two codes have been shifted horizontally to find the closest agreement and found to match to within 200 μm. 



Comparison of the local strain calculations between the two systems was carried out by calculating the maximum normal surface strain across the same line out. The maximum normal surface strain was calculated from the eigenvalues of the Lagrange strain tensor. All points within a 1 mm radius sphere are used to calculate each strain value. Results are shown in Fig. <ref>. The two techniques agree to within 0.2% strain. A projection of the maximum normal surface strain onto the image of the deformed sampled is shown in Figure <ref> for clarity. 



§ LOW STRAIN RATE TESTS



In order to demonstrate the performance of the system a quasi-static compression test and a quasi-static tension test were performed. The first of these involved investigating the shear strain induced in millimeter sized top-hat specimens when subjected to a compressive load. This test utilised the Canon EOS with full resolution and is principally designed to demonstrate the high resolution aspect of the system. 



A flat faced top-hat shaped sample with dimensions given in Figure <ref> was used. The samples were made from near-alpha titanium alloy IMI834, where high rate material behaviour has important relevance with regards to aero-engine discs, and in particular shear-band formation <cit.>. The characteristic shape and size of the sample leads to a small region of high shear strain localised in the region where the top-hat leg intersects the body <cit.>. A compressive force was applied to the top of the sample while the legs of the sample were held still against an anvil. Figure <ref>a-d shows an overlay of the calculated y-displacement (that is displacement in the vertical direction) plotted against global strain, where global strain is defined as ΔL/L with ΔL sample compression and L the initial length. Figure <ref>e shows a spatial average of the shear strain across the high shear region. Here, the shear strain is defined as E_xy, where 𝐄 is the Lagrange strain tensor defined above and x and y correspond to local, in-plane directions shown in Figure <ref>a. As the global compressive strain is increased a shear band can be seen to develop, the width of this shear band, approximately 400 μm, necessitates a high resolution DIC system which in this case is achieved through the use of the K2 lens. Towards the end of the test the shear strain increases further to a maximum observed shear strain of 0.12 and ultimately leads to sample failure. 

In order to test the 3-D aspect of the DIC system a quasti-static tension test was used. A dog-bone shaped sample of low-carbon steel was chosen due to its tendency to neck under tension, a schematic of the sample is given Figure <ref>. The sample was 33.81 mm long with a depth of 5.92 mm and width of 2.96 mm. The DIC system was set-up in order to image the thinner 2.96 mm face, this was chosen to maximise the out-of-plane motion during the test. In order to have a field-of-view that encompasses the full length of the sample the additional 4X magnification used throughout this work was removed. In this case it was found the surface roughness of the sample was such that a speckle pattern did not need to be applied. Tension was applied to the sample at a constant velocity using an Instron 5584, with a 100 kN loadcell, at 1 mm/min. Images of the resulting necking are shown in Figure <ref>, the out-of-plane motion averaged over the highlighted area is given in the insert. Each measurement contains less than 50 μm displacement variation across the sample width. 





§ HIGH STRAIN RATE TESTS





The experimental set-up discussed above has been designed to offer a field of view of a few millimeters with a spatial resolution on the order of a few microns. The system was specifically designed to study the spatial and temporal characteristics of adiabatic shear bands, a phenomena typically associated with high strain-rate deformation and attributed to poor heat dissipation, thermal softening and ultimately material failure. Understanding this mechanism for failure is of the utmost importance in many industrial applications <cit.>.   

To further demonstrate the ease-of-use of the system we used a Phantom 2511 high-rate camera capable of recording 120 000 frames per second (8.3 μs inter-frame time) at a resolution of 240×896 pixels to record a top-hat sample under high strain-rate compression. The high-rate loading was performed on a split-Hopkinson pressure bar <cit.> utilising 6 mm diameter maraging steel bars held in position with frictionless air bearings. A 40 cm long striker capable of supplying a 140 μs loading pulse is used to compress millimeter-sized metallic samples at strain rates ranging from 10^2 to 10^4 s^-1. A schematic of the set-up is shown in Figure <ref>. 

As before flat faced top-hat shaped sample of IMI834 with dimensions given in Figure <ref> was used. A striker velocity of 4.8±0.1 m/s, measured optically, was selected to produce an average global strain rate of 1.2×10^3 s^-1. The 3-D DIC system was used to measure the surface deformation, and in particular the shear strain. Figure <ref>a-c shows the sample with an overlay of the calculated x-displacement (that is displacement in the compressive direction) while Figure <ref>d shows a spatial average of the shear strain across the high shear region. As the global compressive strain is increased a region of high shear can be seen to develop. In Figure <ref> the deformation of the target causes the speckle pattern to become distorted to such a degree that the image correlation technique is no longer able to identify the correct regions causing spatial anomalies in the displacement map. The width of the shear band, approximately 250 μm, is less than observed in the low strain-rate tests. This is expected as the shear band characteristic width is a combination of thermal diffusivity and loading time <cit.>, in the high rate test heat is localised close to the region of maximum shear. As before, towards the end of the test the shear strain increases further to a maximum observed value of 0.16 and ultimately leads to sample failure. The width of the shear bands found in this work are larger than those found previously in Ti alloy samples <cit.> and those predicted by crystal plasticity modelling <cit.> which are typically in the range 10-50 μm, however, the top-hat setup is sensitive to geometric effects, in particular, the overlap of the leg with the body of the top-hat <cit.>. Future work will investigate how the shear-band width varies with parameters such as strain-rate, geometry and microstructure.  



§ CONCLUSIONS



Experimental mechanics is constantly investigating the performance of new materials under load as attempts are made to improve cost, efficiency and performance. A crucial metric is obtaining the stress-strain data from such tests, and the far field measurements typically obtained are supplemented with surface strain measurements. This data is then fed into constitutive material models to aid future development, and can often involve performing tests on non-uniform samples or unique sample geometries. A full 3-D DIC displacement map allows a complete and direct comparison with finite element analysis allowing greater confidence in obtaining elastoplastic parameters such as strength or strain hardening exponents <cit.>. 



A method for obtaining the full-field displacement of both in and out-of-plane motion from a single high frame-rate camera is demonstrated. Using a standard digital SLR (Canon EOS) and a high frame rate camera (Phantom V7.3) together with Infinity K2 DistaMax long-distance microscope lens enabled high resolution images to be taken. 

To achieve this we utilised two open-source codes. The Ncorr two dimensional DIC code is used to perform the correlations between a pair of 2-D images while the MATLAB camera calibration or CalTech camera calibration toolbox is used to calibrate the system through the chequerboard method. The single camera system provides improved temporal synchronisation and spatial stability while the open-source software provides flexibility and transparency without resorting to expensive commercial software or necessitating multiple high speed cameras. 

Tests were performed using a high resolution Canon EOS with 5184x3456 pixels as well as with the Phantom high speed camera at full (800x600) and reduced (240x184) pixel number. In each case the obtainable maximum accuracy in the DIC displacement tests was found to be around 200-300 times better than the optical resolution of the set-up, in line with previous work <cit.>. The software was then compared to the DaVis commercially available code using the same image sets. This removed any ambiguity around illumination, speckle quality, camera resolution or lens quality. The code reproduced the displacement to within 200 um and the maximum normal surface strain to within 0.2. 

The performance of the system was then demonstrated on two quasi-static strain tests. In the first test a millimeter  sized top-hat sample was compressed to induce large shear strains in the region between the legs and body of the top-hat sample. The high resolution of the DIC system enabled a 400 μm band of high shear strain to be identified and quantified before failure. The second quasi-static test utilised a dog-bone sample in a tension test to observe necking in the Z direction. Large displacements were observed and the necking was measured to within 50 μm uncertainty.  Finally, a high-strain test on the same top-hat samples as compressed quasi-statically was carried out on a split-Hopkinson pressure bar at a global strain rate of 1.2×10^3 s^-1. Images captured at 120 000 Hz showed a thinner region of high-shear approximately 250 μm wide, suggesting adiabatic conditions were reach in the high rate tests.

Future work will utilise the single camera system and codes to further investigate the physics surrounding dynamic events on the sub-millimeter scale such as grains, shear bands and localised surface strains. 



§ ACKNOWLEDGEMENTS


This research was supported by EPSRC grant Heterogeneous Mechanics in Hexagonal Alloys across Length and Time Scales (EP/K034332/1). In addition, the authors are very grateful to David Rugg and Rolls-Royce plc for provision of Ti material.



9

chu
T. Chu, W. Ranson, M. Sutton, Exp. Mech. 25(3), (1985) 232244

sutton
M. A. Sutton, M. Cheng, W. H. Peters, Y. J. Chao, S. R. McNeill, Image Vision Comput. 4(3), (1986) 143150 

interferometry
D. Post, Exp. Mech. 23(2), (1983) 203-210

holography
T. D. Dudderar, H. J. Gorman, Exp. Mech. 13(4), (1973) 145-149

suttonbook
H. Schreier, J. J. Orteu, M. A. Sutton, Digital Image Correlation for Shape and Deformation Measurements (Springer, Verlag 2009) 

luo
P. Luo, Y. Chao, M. Sutton, W. Peters, Exp. Mech. 33(2), (1993) 123132

vic3d
 VIC-3D system (Correlated Solutions, Columbia, South Carolina)

aramis
ARAMIS system (GOM Company, Braunschweig, Germany) 

davis
DaVis system (LaVision, Goettingen, Germany)

singlecamera2
C. J. Tay, C. Quan, Y. H. Huang, Y. Fu, Opt. Commun. 251, (2005) 2336

singlecamera3
C. Quan, C. J. Tay, W. Sun, X. He, Appl. Opt. 47, (2008) 583593

singlecamera4
W. Sun, E. Dong, X. He, Proc. SPIE 6723, (2007) 37230

Prentice
H. J. Prentice, W. G. Proud, AIP Conf. Proc. 845 (2006) 1275-1278

singlecamera1
M. Pankow, B. Justusson, A. M. Waas, Applied Optics 49(17), (2010) 3418-3427 

CameraCal
J. Y. Bouguet, Camera Calibration Toolbox for Matlab (2013) Available at: http://www.vision.caltech.edu/bouguetj/calib_doc/index.html

Ncorr1
J. Blaber, A. Antoniou, Ncorr: Open-source 2D Digital Image Correlation Matlab Software (2014)  Available at: http://www.ncorr.com

Ncorr3
J. Blaber, B. Adair, A. Antoniou, Exp. Mech. 55 (2015) 1105

Ncorr2
R. Harilal, M. Ramji, Adaptation of Open Source 2D DIC Software Ncorr for Solid Mechanics Applications 9th International Symposium on Advanced Science and Technology in Experimental Mechanics (2014)

bdodd
B. Dodd and Y. Bai Adiabatic Shear Localization (Elsevier (2012))

mechconmed
L. E. Malvern Introduction to the Mechanics of a Continuous Medium (New Jersey: Prentice-Hall, (1969))

aerospace
I. Balasundara, T. Raghua, B. P. Kashyapb, Prog. Nat. Sci.23(6) (2013) 598607

speckle
T. A. Berfield, J. K. Patel, R. G. Shimmin, P. V. Braun, J. Lambros, N. R. Sottos, Exp. Mech. 47, (2007) 51-62

matlab
MATLAB version 8.3.0.532 (R2014a). Natick, Massachusetts: The MathWorks Inc., 2014.

PrenticeThesis
H. J. Prentice,  Development of Stereoscopic Speckle Photography Techniques for Studies of Dynamic Plate Deformation (Magdalene College, University of Cambridge, 2006)

shpb1
H. Kolsky, Proc. Phys. Soc. London, B62, (1949) 676.

shpb2
B. A. Gama, S. L. Lopatnikov and J. W. Gillespie, Appl. Mech. Rev 57(4), (2004) 223-250 

shpb3
A. Mohammadhosseini, S. H. Masood, D. Fraser and M. Jahedi, Adv. Manuf. (3), (2015) 232

tophat1
J. Peirs, P. Verleysen, J. Degrieck, and F. Coghe, International Journal of Impact Engineering 37(6), (2010) 703-714

tophat2
R. Clos, U. Schreppel, U., and P. Veit, Journal de Physique, 110(4), (2003) 111-116

tophat3
X. Teng, T. Wierzbicki, and H. Couque, Mechanics of Materials, 39(2), (2007) 107-125

stevesmith
S. W. Smith The Scientist and Engineer's Guide to Digital Signal Processing (California Technical Pub; 1st edition (1997))

3ddicfea
G. Li, F. Xu, G. Sun, Q. Li, Int. J. Adv. Manuf. Technol. 74 (2014) 893905

zhen
Z. Zhang, D. E. Eakins, F. P. E. Dunne, International Journal of Plasticity 79 (2016) 196-216

bandwidth1
S. Liao, J. Duffy, Journal of the Mechanics and Physics of Solids 46(11) (1998) 2201-2231

bandwidth2
F. Coghe, L. Rabet, L. Kestens, Journal De Physique. IV 134 (2006) 845-850


