F Pierron

M Grédiac

Towards Material Testing 2.0

Keywords: Virtual Fields Method, Finite Element Model Updating, inverse identification, test design, optimization, full-field measurements, MT2.0

Full-field optical measurements like Digital Image Correlation or the Grid Method have brought a paradigm shift in the experimental mechanics community. While inverse identification techniques like Finite Element Model Updating (FEMU) or the Virtual Fields Method (VFM) have been the object of significant developments, current test methods, inherited from the age of strain gauges or LVDTs, are generally not well adapted to the rich information provided by these new measurement tools. This paper provides a review of the research dealing with the design and optimization of heterogeneous mechanical tests for identification of material parameters from full-field measurements, christened here Material Testing 2.0 (MT2.0). 1

Introduction

Full-field measurement techniques such as Digital Image Correlation (DIC, [START_REF] Sutton | Image Correlation for Shape, Motion and Deformation Measurements[END_REF]) or the Grid Method (GM, [START_REF] Grédiac | The grid method for in-plane displacement and strain measurement: a review and analysis[END_REF]) have become mature and are now widespread in the experimental mechanics community. This is mainly due to the fact that these techniques are able to provide, within certain limits, displacement and strain fields that occur at the surface of specimens of various shapes subjected to thermomechanical loadings. In a context of material characterization, such measurements make it possible to move from simple statically determinate tests to more complex ones, with a view to get more data from a single test, thus reducing the testing effort and increasing the quality of the material models. It should be noted that for mesoscopically-heterogeneous materials such as welds, large weave composites or graded foams, the traditional test methods only provide an averaged behaviour while full-field measurement have the potential to identify spatially-varying properties.

The problem with this approach however is that in general, no closed-form expression linking local kinematic measurements to constitutive parameters is available. Therefore, early research focused on developing robust and efficient techniques to identify material parameters from full-field measurements in that case. The need for optimizing test configurations adapted to this new paradigm emerged later, as the breadth of applications increased.

Using more complex tests, full-field measurements and inverse identification to identify mechanical constitutive parameters is a new paradigm in material testing that for the moment has been confined to research laboratories. There is no doubt that in the coming years, such approaches will develop further and will eventually lead to new test methods widely adopted by both academia and industry. For the sake of providing visibility to this technical area, the authors of this paper feel that a new denomination has to be introduced. It is proposed to refer to this as 'Material Testing 2.0' (MT2.0), capturing the fact that this will lead to a complete revamp of mechanical tests, eventually leading to new standards as cameras gradually replace strain gauges and extensometers in testing laboratories. This denomination, MT2.0, will be used in the rest of the paper for the sake of clarity and brevity.

The objective of this paper is to propose a review of the research dealing with the design and optimization of mechanical tests dedicated to MT2.0.

It is worth noting that optimizing test configurations has been a research topic for a long time. However, with the classical approach of mechanical testing, which solely uses a limited number of isolated sensors, the focus has generally been on optimizing the placement and orientation of the sensors to collect the most relevant information. Optimal placement of strain gauges [START_REF] Masroor | Designing an all-purpose force transducer[END_REF][START_REF] Gupta | Input load identification from optimally placed strain gages using D-optimal design and model reduction[END_REF], accelerometers [START_REF] Shah | A methodology for optimal sensor locations for identification of dynamic systems[END_REF][START_REF] Kirkegaard | Optimal Design of Experiments for Parametric Identification of Civil Engineering Structures[END_REF][START_REF] Kammer | Optimal placement of triaxial accelerometers for modal vibration tests[END_REF], or even GPS sensors on wide structures like bridges [START_REF] Meo | On the optimal sensor placement techniques for a bridge structure[END_REF] has been addressed in the literature. The proposed methodologies generally rely on statistical properties that have to be satisfied [START_REF] Ferodov | Theory of optimal experiments[END_REF]. They are of general purpose, and can be employed in other fields than experimental mechanics [START_REF] Oh | Determination of the minimal number and optimal sensor location in a nuclear system with fixed incore detectors[END_REF][START_REF] Naimimohasses | Sensor optimization using neural network sensitivity measures[END_REF][START_REF] Padula | Optimal sensor/actuator locations for active structural acoustic control[END_REF]. However, since full-field measurements provide spatially dense data, such methods devised for a small number of data points are generally irrelevant. The best placement of cameras for 3D measurements or the best location of displacement sensors still represents an open problem when full-field measurements are employed or when constitutive parameters are identified by inverse analysis, [START_REF] Vitse | Toward virtual design and optimization of a structural test monitored by a multi-view system[END_REF] and [START_REF] Huang | A study of the optimal measurement placement for parameter identification of orthotropic composites by the boundary element method[END_REF], respectively. But the present review will focus exclusively on designing specific specimen geometry and loading configuration for identification of constitutive parameters from either displacement or strain fields.

In this review, the papers have been classified according to four different categories following a progression of sophistication:

• Design by intuition • Design by strain state • Design by identification quality

• Design by full identification simulation

The first level, 'design by intuition', consists in devising test configurations arbitrarily based on previous experience or intuition. However, as the complexity of the models increases, together with the number of parameters to be identified, this approach generally leads to considerably sub-optimal configurations and more rational approaches to test design have to be proposed. A natural idea, referred to here as 'design by strain state', is to use the strain state as a design variable, as this bears witness of the material parameters. By trying to balance the different components, one can hope to improve identifiability. However, the dependence of the sought parameters to the strain components can be complex and a more relevant approach, called here 'design by identification quality', consists in using the quality of the identification as an optimization criterion, generally using kinematic data directly from numerical simulations. Nevertheless, this does not take into account important features coming from the measurements themselves, including hardware performance and image processing parameters. Therefore, the most sophisticated approach to date, 'design by full identification simulation', consists in simulating the whole identification chain, from image formation all the way to identification. This allows for full uncertainty propagation to be captured and for more realistic identification performance to be derived, which can then be used for test design and optimization.

The present paper is organized as follows. First, the two main identification techniques currently in use are briefly reviewed, namely, Finite Element Model Updating (FEMU) and the Virtual Fields Method (VFM). Then, representative articles in the four different categories of tests mentioned above are presented and discussed in separated sections. Some perspectives and recommendations are finally offered.

Main procedures used to extract constitutive parameters from full-field measurements

Before designing MT2.0 tests, the first step has been to elaborate procedures that could be used to extract constitutive parameters from heterogeneous strain fields, in cases for which no closed-form expression exists between these unknown parameters and the local deformation measurements. FEMU and the VFM have progressively imposed themselves as the two main techniques for this. The main features of both methods are briefly recalled hereafter.

FEMU

To the best of the authors' knowledge, the idea of updating a finite element model in order to identify material properties has been given for the first time in the early 70's in [START_REF] Kavanagh | Finite element applications in the characterization of elastic solids[END_REF], in the case of linear and non-linear elasticity. Updating is performed by minimizing a cost function with respect to the sought parameters, this cost function reflecting the "distance" between displacements measured at some points, and their counterparts estimated with FE computation at a given iteration. A cross-shaped specimen subjected to two perpendicular loads was used to this end in this reference. At that time, optimizing the test configuration was not the objective, but the underpinning idea was that the location of the applied load and the geometry of the specimen would be such that its response was "sufficiently" influenced by all the unknowns, so that they could be identified. A hole was added in order to increase heterogeneity. It is worth noting that FEMU does not necessarily require full-field measurements. Only numerical simulations were presented in [START_REF] Kavanagh | Finite element applications in the characterization of elastic solids[END_REF], but this paper served as a starting point to many subsequent studies where this idea was employed with actual measurements.

For instance, shortly after this first reference, its first author used input data provided by a limited number of strain gauges [START_REF] Kavanagh | Extension of classical experimental techniques for characterizing composite-material behavior[END_REF]. It seems that [START_REF] Courtade | Caracterisation élastique d'un stratifié résine-verre par la méthode des déplacements[END_REF] is the first example where FEMU input data are full-field measurements. Composite materials were tested in this reference. This was the starting point of many developments covering various materials and types of constitutive equations. Examples dealing with biological materials [START_REF] Oomens | A numerical-experimental method for a mechanical characterization of biological materials[END_REF][START_REF] Boyer | Assessment of the in-plane biomechanical properties of human skin using a finite element model updating approach combined with an optical full-field measurement on a new tensile device[END_REF], metals [START_REF] Meuwissen | Determination of the elasto-plastic properties of aluminium using a mixed numerical-experimental method[END_REF][START_REF] Kajberg | Characterisation of materials subjected to large strains by inverse modelling based on in-plane displacement fields[END_REF][START_REF] Pottier | Contribution of heterogeneous strain field measurements and boundary conditions modelling in inverse identification of material parameters[END_REF][START_REF] Réthoré | Robust identification of elasto-plastic constitutive law parameters from digital images using 3D kinematics[END_REF][START_REF] Mathieu | Estimation of elastoplastic parameters via weighted FEMU and integrated-DIC[END_REF][START_REF] Ruybalid | Comparison of the identification performance of conventional FEM updating and integrated DIC[END_REF],

MEMS [START_REF] Amiot | Identification of elastic property and loading fields from full-field displacement measurements[END_REF], composites [START_REF] He | Characterization of nonlinear shear properties for composite materials using digital image correlation and finite element analysis[END_REF][START_REF] Makeev | A method for measurement of multiple constitutive properties for composite materials[END_REF][START_REF] Passieux | Multiscale displacement field measurement using digital image correlation: Application to the identification of elastic properties[END_REF][START_REF] He | Characterization of stress-strain behavior of composites using digital image correlation and finite element analysis[END_REF], heterogeneous [START_REF] Kroon | A numerical framework for material characterisation of inhomogeneous hyperelastic membranes by inverse analysis[END_REF] and homogeneous [START_REF] Drapier | Identification strategy for orthotropic knitted elastomeric fabrics under large biaxial deformations[END_REF] hyperelastic materials, wood [START_REF] Réthoré | Multiscale digital image identification of heterogeneous elastic properties of softwoods[END_REF], and even crimped mineral wool [START_REF] Witz | Mechanical Properties of Crimped Mineral Wools: Identification From Digital Image Correlation[END_REF], can be found in the literature. Note finally that when DIC is the measuring technique, it is possible to merge the two steps involved in the procedure described above, namely measuring the displacement field on the one hand, and then applying FEMU on the other hand. This leads to the so-called Integrated DIC (IDIC) introduced in [START_REF] Mathieu | Estimation of elastoplastic parameters via weighted FEMU and integrated-DIC[END_REF]. In this reference, the authors find that merging or not these two steps gives equivalent results, while [START_REF] Ruybalid | Comparison of the identification performance of conventional FEM updating and integrated DIC[END_REF] claims the integrated version leads to more reliable results. It seems however that this approach has not yet been widely used, so we consider here that IDIC is a particular case of FEMU coupled with DIC measurements.

The main strength of FEMU is that it relies on the finite element method, which is a versatile tool widely used in the solid mechanics community. This is also its main weakness. Indeed, FE was initially developed for a different problem, namely finding the displacement/strain/stress fields by assuming that the constitutive equations and their governing parameters are known. Using this tool to solve the material identification problem automatically leads the identification procedure to be iterative, even in the simplest case of linear elasticity. As a consequence, FEMU is generally quite computationally intensive, particularly for geometrical or material non-linearities. Finally, boundary conditions are necessary to perform any FE calculation, and feeding FEMU with faithful experimental boundary conditions is often challenging, especially in terms of force distributions. DIC-measured displacements at the boundary of the field of view can be used but their noisy nature also creates specific problems.

The VFM

The drawbacks of FEMU have led researchers to elaborate alternative procedures which could straightforwardly (instead of iteratively) provide the unknown parameters. In this regard, a pioneering approach was proposed in 1986 in [START_REF] Foudjet | Méthodes de caractérisation de matériaux anisotropes[END_REF]. It consisted in using the Maxwell-Betti reciprocal work theorem in order to identify the six independent bending stiffnesses of wood plates. Indeed, assuming linear elasticity for the constitutive materials, considering a simple square plate subjected to three different loads and mixing the three corresponding deflection fields leads to six independent equations since the work done by the forces applied in test #i through the displacement of test #j is equal to the work done by the forces applied in test #j through the displacement of test #i. The system formed by these six equations being linear, the six unknowns are found after a mere inversion. The drawback is that even if only one specimen is sufficient to identify six independent parameters with this approach, three different tests must be performed on the same specimen to obtain a set of independent equations. Another point is that measurements are always noisy, and two measurement fields are involved in each equation, which makes this approach rather noise-sensitive.

These two limitations (i-three different tests, and ii-results highly sensitive to noise) have led a related approach to emerge in [START_REF] Grédiac | Principe des travaux virtuels et identification/principle of virtual work and identification[END_REF][START_REF] Grédiac | A new method for determination of bending rigidities of thin anisotropic plates[END_REF]. In these references, instead of using the Maxwell-Betti reciprocal work theorem, it was proposed to consider the principle of virtual work, which represents the weak form of equilibrium. Indeed, an interesting feature of the principle of virtual work is that it is rigorously satisfied for any continuous and piecewise differentiable virtual field, thus for an infinite number of such fields. The idea is therefore to combine a unique actual displacement field measured during a test with different but independent virtual displacement fields arbitrarily selected by the user. In the case of linear elasticity and if the actual displacement field is such that all the unknown parameters are involved in the response of the specimen, it can be shown that choosing at least as many independent virtual fields as unknowns leads to a system of linear equations, which directly provides the unknowns after inversion. Compared to the approach based on the Maxwell-Betti reciprocal work theorem, a single test is potentially sufficient with the VFM. In addition, only one actual measurement field affected by noise instead of two is involved, which makes this identification technique more robust. Last but not least, the principle of virtual work is valid whatever the type of constitutive equations, so even in the case of non-linearities such as large deformations or plasticity, which is not the case of the Maxwell-Betti reciprocal work theorem. These features have led VFM-related research to expand, with mainly three types of contributions. The first type deals with the improvement of the method itself, mainly concerning the choice of the virtual fields [START_REF] Grédiac | Special virtual fields for the direct determination of material parameters with the Virtual Fields Method. 1-Principle and definition[END_REF][START_REF] Avril | Sensitivity of the Virtual Fields Method to noisy data[END_REF][START_REF] Nguyen | A Fourier-series-based Virtual Fields Method for the identification of 2-d stiffness and traction distributions[END_REF][START_REF] Marek | Sensitivity-based virtual fields for the non-linear Virtual Fields Method[END_REF][START_REF] Marek | Extension of the sensitivity-based virtual fields to large deformation anisotropic plasticity[END_REF][START_REF] Tayeb | Sensitivity-based virtual fields for identifying hyperelastic constitutive parameters[END_REF][START_REF] Marek | Experimental validation of the sensitivity-based virtual fields for identification of anisotropic plasticity models[END_REF] since this choice directly influences the robustness of the method. The second type concerns applications involving more sophisticated constitutive equations. Anisotropic elasticity for composite materials [START_REF] Grédiac | A general method for data averaging of anisotropic elastic constants[END_REF][START_REF] Chalal | Experimental identification of a nonlinear model for composites using the grid technique coupled to the Virtual Fields Method[END_REF][START_REF] Moulart | Identification of the through-thickness rigidities of a thick laminated composite tube[END_REF] and wood [START_REF] Xavier | Measuring orthotropic bending stiffness components of pinus pinaster by the Virtual Fields Method[END_REF], elasticity or hyperlasticity of elastomers [START_REF] Promma | Application of the Virtual Fields Method to mechanical characterization of elastomeric materials[END_REF][START_REF] Guélon | A new characterization method for rubber[END_REF] and biological materials [START_REF] Zhang | Verification of a Virtual Fields Method to extract the mechanical properties of human optic nerve head tissues in vivo[END_REF][START_REF] Avril | Anisotropic and hyperelastic identification of in vitro human arteries from full-field optical measurements[END_REF][START_REF] Meng | 3d strain and elasticity measurement of layered biomaterials by optical coherence elastography based on digital volume correlation and virtual fields method[END_REF][START_REF] Miller | Estimation of transversely isotropic material properties from magnetic resonance elastography using the optimised virtual fields method[END_REF], elastoplasticity of plain metals [START_REF] Grédiac | Applying the Virtual Fields Method to the identification of plastic constitutive equations[END_REF][START_REF] Kim | Characterization of the post-necking strain hardening behavior using the Virtual Fields Method[END_REF][START_REF] Fu | Identification of nonlinear kinematic hardening constitutive model parameters using the Virtual Fields Method for advanced high strength steels[END_REF][START_REF] Rossi | Application of the Virtual Fields Method to large strain anisotropic plasticity[END_REF][START_REF] Martins | Calibration of anisotropic plasticity models using a biaxial test and the virtual fields method[END_REF][START_REF] Rossi | A general linear method to evaluate the hardening behaviour of metals at large strain with full-field measurements[END_REF][START_REF] Valeri | Determining the tensile response of materials at high temperature using DIC and the Virtual Fields Method[END_REF] or welds [START_REF] Le Louëdec | Identification of the local elasto-plastic behavior of FSW welds using the Virtual Fields Method[END_REF], are typical examples. The use of the VFM to characterize some materials which are less widely used in mechanical engineering, such as paperboard [START_REF] Considine | General anisotropy identification of paperboard with Virtual Fields Method[END_REF] or polymeric foams [START_REF] Wang | Optimised experimental characterisation of polymeric foam material using DIC and the Virtual Fields Method[END_REF], has also been reported in recent papers. The effect of time was also considered with vibrations [START_REF] Grédiac | Direct identification of elastic constants of anisotropic plates by modal analysis : theoretical and numerical aspects[END_REF][START_REF] Grédiac | Direct identification of elastic constants of anisotropic plates by modal analysis : experiments and results[END_REF]. This made it possible to characterize viscoelastic materials [START_REF] Giraudeau | Identification of stiffness and damping properties of thin isotropic vibrating plates using the Virtual Fields Method: theory and simulations[END_REF][START_REF] Giraudeau | Stiffness and damping identification from full field measurements on vibrating plates[END_REF]. Characterizing the response of various materials at high strain rates constitutes the most recent example of use of this identification method [START_REF] Seghir | A novel image-based ultrasonic test to map material mechanical properties at high strain-rates[END_REF][START_REF] Pierron | Beyond Hopkinson's bar[END_REF]. Full details on the VFM can be found in [START_REF] Pierron | The Virtual Fields Method[END_REF]. The third type of papers on the VFM deals with test design, which is the topic addressed herein. [START_REF] Gillard | The application of digital volume correlation (dvc) to study the microstructural behaviour of trabecular bone during compression[END_REF][START_REF] Connesson | Extension of the optimised virtual fields method to estimate viscoelastic material parameters from 3D dynamic displacement fields[END_REF]. This means that assumptions are required to deduce the through-thickness strain distribution from the measurements on the surface. Such assumptions are reasonable only in the case of thin specimens, while more complex 3D specimens can only be treated with FEMU. Another point is that representative boundary conditions have to be used with finite element calculations but they are not always available. On the other hand, no exact knowledge of boundary conditions is needed with the VFM if suitable virtual fields are selected. The quality of the results (in terms of robustness) obtained with both identification techniques is globally similar, but the main advantage of the VFM is that no resolution of the direct problem is necessary, which considerably reduces computation times. As an example, it was reported that for comparable identification performance, the VFM was 125 times faster than FEMU in the case of biomaterials exhibiting a hyperelatic response [START_REF] Zhang | Verification of a Virtual Fields Method to extract the mechanical properties of human optic nerve head tissues in vivo[END_REF] and 428 or 263 times faster for two cases of elasto-plasticity [START_REF] Martins | Comparison of inverse identification strategies for constitutive mechanical models using full-field measurements[END_REF].

VFM vs. FEMU

With a view to designing optimal testing conditions, it is clear that refining specimen shape or loading configuration requires many resolutions of the parameter identification problem. In that case, reduced computation times for the identification procedure becomes a significant asset.

Let us now examine in which way FEMU and the VFM have been used for designing tests in the literature.

Design by intuition

A simple route to obtain heterogeneous strain fields is to "recycle" pre-existing tests, by keeping them as such or by slightly adapting them. For instance, in one of the very first papers on this topic [START_REF] Meuwissen | Determination of the elasto-plastic properties of aluminium using a mixed numerical-experimental method[END_REF], it was proposed to perform a tensile test on a parallelepipedic specimen.

Heterogeneity was induced by two asymmetric notches machined on each side of the specimen.

The in-plane displacement was measured with a video tracking system, at 120 retro-reflective markers regularly deposited onto the front face of the specimen. From these measured displacements, the authors identified the yield stress of both the isotropic Von Mises and the orthotropic Hill yield criteria by using FEMU. A similar specimen geometry was discussed in [START_REF] Pottier | Contribution of heterogeneous strain field measurements and boundary conditions modelling in inverse identification of material parameters[END_REF] with the same objective, but DIC was used in this case to perform the measurements.

In [START_REF] Guner | Characterization of anisotropy of sheet metals employing inhomogeneous strain fields for Yld2000-2D yield function[END_REF], heterogeneity was induced by a varying cross-section and the planar anisotropy of the plastic response was identified. Another intuitive approach is to use an open-hole tensile test [START_REF] Lecompte | Identification of elastic orthotropic material parameters based on ESPI measurements[END_REF][START_REF] Silva | Identification of material properties using FEMU: Application to the open hole tensile test[END_REF][START_REF] Kowalczyk | Identification of mechanical parameters of composites in tensile tests using mixed numerical-experimental method[END_REF][START_REF] Gogu | Bayesian identification of elastic constants in multi-directional laminate from moiré interferometry displacement fields[END_REF] as the specimens are easy to manufacture and provide strain heterogeneity. The displacement field, measured with ESPI (Electronic Speckle Pattern Interferometry) in [START_REF] Lecompte | Identification of elastic orthotropic material parameters based on ESPI measurements[END_REF],

DIC in [START_REF] Silva | Identification of material properties using FEMU: Application to the open hole tensile test[END_REF][START_REF] Kowalczyk | Identification of mechanical parameters of composites in tensile tests using mixed numerical-experimental method[END_REF] and moiré interferometry in [START_REF] Gogu | Bayesian identification of elastic constants in multi-directional laminate from moiré interferometry displacement fields[END_REF], was used to feed FEMU with experimental data and identify the four in-plane stiffnesses of an orthotropic composite material. The authors of the latter reference note that various sources of uncertainty impair the quality of the results, and thus pay attention to quantifying the uncertainty with which the parameters are identified by using a Bayesian approach. A double notched specimen was also employed in [START_REF] Kowalczyk | Identification of mechanical parameters of composites in tensile tests using mixed numerical-experimental method[END_REF], and the authors observed that Poisson's ratio was more easily identified with this second geometry.

An extension of the classical tensile test is the biaxial tensile test, which gives rise to a heterogeneous strain field at the centre of cruciform specimens. This test has been used in various cases in the literature:

• composite materials have been characterized in [START_REF] Lecompte | Mixed numerical-experimental technique for orthotropic parameter identification using biaxial tensile tests on cruciform specimens[END_REF] by using DIC and FEMU. Results from specimens with and without a hole in the centre are compared. Interestingly, the authors report that results are better without the hole, presumably because of the bias due an insufficient spatial resolution arising from the concentrated strains around the hole. This remark justifies that a complete optimization procedure for test design should ideally account for the metrological performance of the measuring technique, as discussed in Section 6 below;

• the parameters of a hyperelastic law of an elastomeric material are identified in [START_REF] Promma | Application of the Virtual Fields Method to mechanical characterization of elastomeric materials[END_REF] with DIC and VFM. Compared to the preceding case, the strain level is much higher, which makes the signal-to-noise ratio more favourable;

• the parameters governing an elastoplastic law were identified in [START_REF] Cooreman | Identification of mechanical material behavior through inverse modeling and DIC[END_REF][START_REF] Bertin | Integrated digital image correlation applied to elastoplastic identification in a biaxial experiment[END_REF]. The geometry of the fillets between the branches of a cruciform specimens was adjusted in the latter reference to improve the quality of the identification procedure, and integrated DIC (IDIC, [START_REF] Mathieu | Estimation of elastoplastic parameters via weighted FEMU and integrated-DIC[END_REF]) was employed.

Characterizing plasticity parameters of a mild steel was proposed in [START_REF] Mahnken | A unified approach for parameter identification of inelastic material models in the frame of the finite element method[END_REF] by using a compact tension specimen (CT). This test is generally employed to study crack propagation in metals. In this reference, the authors took advantage of the strongly heterogeneous strain field which occurs in such specimens, in order to identify the parameters of a viscoplastic power law with isotropic hardening. The identification method is FEMU and displacement/strain fields were measured by tracking the intersection points between the lines of a grating deposited onto the specimen before the test. More recently, a 'bridge-like' specimen was used in conjunction with the VFM [START_REF] Fu | An identification method for anisotropic plastic constitutive parameters of sheet metals[END_REF] to identify a Hill48 anisotropic plasticity criterion. Unfortunately, the steps to reach this configuration are not detailed in the paper.

The three-point bending test has also been widely used to characterize materials in experimental mechanics. Short beam bending activates the whole set of parameters governing the in-plane response of orthotropic composites, leading to possible identification as illustrated in [START_REF] He | Characterization of nonlinear shear properties for composite materials using digital image correlation and finite element analysis[END_REF][START_REF] Makeev | A method for measurement of multiple constitutive properties for composite materials[END_REF][START_REF] Gras | Identification of the out-of-plane shear modulus of a 3D woven composites[END_REF][START_REF] He | Characterization of stress-strain behavior of composites using digital image correlation and finite element analysis[END_REF][START_REF] Liu | Feasibility and uncertainty analysis of constitutive parameters identification for composite materials using displacement field data[END_REF] with DIC and FEMU.

Another similar approach is to rely on a standardized shear test developed for composites, namely the Iosipescu, or double V-notch, shear test. The standard specimen geometry features a double symmetric notch in the gauge section to ensure a dominant and near-uniform strain field between the two notches [START_REF] Walrath | The losipescu shear test as applied to composite materials[END_REF]. In the same vein as for the tensile and bending tests discussed above, this shear test has been recycled by using the standard fixture, but removing the specimen notches in an attempt to create a combined state of bending and shear. This specimen thus behaves like a short double-clamped beam. As such and despite the dominance of the shear strain, the four stiffnesses of an orthotropic composite influence the strain field between the clamps. This property has been taken advantage of to characterize various orthotropic materials such as composites [START_REF] Pierron | Identification of the through-thickness moduli of thick composites from whole-field measurements using the Iosipescu fixture: theory and simulations[END_REF][START_REF] Chalal | Experimental identification of a nonlinear model for composites using the grid technique coupled to the virtual fields method[END_REF] or wood [START_REF] Xavier | Variation of transverse and shear stiffness properties of wood in a tree[END_REF], the VFM being the identification technique in these cases. It is worth noting that in the latter case, an off-axis configuration was used to enhance identifiability, as also favoured in [START_REF] Kretschmann | Optimized test design for identification of the variation of elastic stiffness properties of loblolly pine (pinus taeda) pith to bark[END_REF].

Compression of a circular disk (referred to as the Brazilian test) is a well-known test for brittle materials such as rock or concrete. In [START_REF] Cao | Application of moiréinterferometry to the characterization of orthotropic materials in the antisymmetric configuration using the Virtual Fields Method[END_REF], the authors showed that the heterogeneous strain field observed with moiré interferometry during this test could be used to characterize an orthotropic carbon fibre composite. In [START_REF] Huchzermeyer | The application of full-field techniques to estimate both tensile and fracture properties: an investigation into modifications to standard sample geometries[END_REF], it was proposed to use slotted disks to increase heterogeneity in the strain field, and to deduce the Young's modulus and Poisson's ratio of an isotropic material.

In conclusion, while sidetracking standard tests to adapt them to the new paradigm of MT2.0 is a first natural approach, it generally falls short when more complex behaviours have to be identified as it is less intuitive to ensure that all parameters are involved in the response.

Also, even when identification is successful, the configurations are generally sub-optimal and a more rational approach to test design has to be devised.

Design by strain state

The aim here is to present some examples of test design found in the literature, which are governed by the fact that certain types of strain states should take place within a specimen to activate all the unknown parameters for successful identification. The design variables are mainly some parameters which govern the shape of the specimen after suitable parametrization, the fixture being then adjusted to enable attachment of and load transfer to the specimen.

Compared to the preceding examples, non-conventional specimens and/or loading systems are used, and finite element computations are employed to optimize design variables and check that the desired strain states are obtained. This category has been divided into two subcategories: test configurations designed intuitively on strain distribution; and test configurations obtained through actual shape optimization with strain-distribution-based cost function.

Intuition-led design based on strain state

A few examples of test designs based on strain states can be found in the literature, where unusual specimen shapes have been intuitively defined based on considerations on strain states. One of the examples deals with a punch test designed to identify an anisotropic Hill48 surface coupled to a Ludwig strain hardening model [START_REF] Pottier | Out-of-plane testing procedure for inverse identification purpose: application in sheet metal plasticity[END_REF]. The specimen was designed in such a way that various relevant strain states coexisted, each of them occurring in a different part of the specimen. These strain states are biaxial tension, uniaxial tension and shear. Figure 1 shows the geometry of the specimen as well as the strain field obtained with a finite element model. The outer periphery of the specimen is clamped while the load is applied at the centre, perpendicular to the plate, as in a punch test. The specimen is circular, so measurements along both the rolling and transverse directions are possible with the same specimen. DIC was used to measure the 3D displacement and strain at the surface of the specimen, and six constitutive parameters were identified with displacement-and force-based FEMU, two for the Ludwig hardening (the initial yield stress was obtained separately), and four for the Hill48 criterion. This test configuration was used recently by other authors to calibrate a Yld2000

anisotropic criterion [START_REF] Wang | Identification of 7b04 aluminum alloy anisotropy yield criteria with conventional test and Pottier test at elevated temperature[END_REF].

Figure 1: Specimen designed for punch test, after [START_REF] Pottier | Out-of-plane testing procedure for inverse identification purpose: application in sheet metal plasticity[END_REF]. E xx , E yy and E xy are the in-plane components of the Green-Lagrange strain tensor.

In the same spirit, hyperelastic materials were considered in [START_REF] Guélon | A new characterization method for rubber[END_REF]. Classically, three homogeneous tests are performed to characterize hyperelastic materials, namely uniaxial tensile (UT), pure shear (PS) and equibiaxial tensile (ET) tests [START_REF] Holzapfel | Nonlinear Solid Mechanics[END_REF]. Note that the three corre-sponding strain states coexist in the biaxial test discussed in the preceding section but in [START_REF] Guélon | A new characterization method for rubber[END_REF],

it was proposed to combine these three types of strain states in a single test which derives from the classic tensile test, merely by adding to the specimen a third branch perpendicular to the two of the classic geometry of any tensile specimen, see Figure 2-a. Compared to the biaxial test on cruciform specimens, the benefit is to obtain these three strain states with a uniaxial tensile machine, which is more readily available in testing labs than a biaxial one.

This third branch was attached to a grip fixed to the frame of the test machine, see Fig- An appealing aspect of MT2.0 is that specimens featuring non-standard geometries can be tested, thus enabling the measurement of actual properties of components instead of those measured on specimens elaborated specifically for testing purposes. This is especially relevant in the case of materials which properties are sensitive to the manufacturing process.

In [START_REF] Moulart | Identification of the through-thickness rigidities of a thick laminated composite tube[END_REF] for instance, it was proposed to measure the elastic through-thickness parameters of thick composites tubes. It is clear that in this case, even cutting parallelepipedic specimens through the thickness of such tubes and testing them with classical homogeneous tests is irrelevant because of the curvature of the fibres. In [START_REF] Moulart | Identification of the through-thickness rigidities of a thick laminated composite tube[END_REF], it was proposed to cut a thin ring in such tubes, to test it under compression, and to consider only a portion of this ring to perform identification, see the grey zone in Figure 3 below. The location of this zone on the ring has been chosen in such a way that the amplitude of the three in-plane strain components are balanced on average. The aim was to ensure better use of the camera pixels as imaging a larger part of the ring, or even its entirety, would result in the loss of many pixels thus challenging spatial resolution. The same procedure was employed later to obtain the equivalent orthotropic elastic properties of superconducting winding rings (or 'pancakes'), constitutive of a large MRI magnet [START_REF] Kim | Characterizing elastic properties of superconducting windings by simulations and experiments[END_REF].

Figure 3: Compression test on a ring cut in a composite tube, after [START_REF] Moulart | Identification of the through-thickness rigidities of a thick laminated composite tube[END_REF].

Finally, a D-shaped specimen was designed in [START_REF] Jones | Parameter covariance and non-uniqueness in material model calibration using the virtual fields method[END_REF] to identify the parameters governing a visco-plastic law. The specimen is plane and its D-shaped geometry is unconventional, as can be seen in Figure 4. It was loaded vertically in tension using a standard uniaxial testing machine. This geometry was defined manually but iteratively, starting with a shape chosen by intuition, then by performing FE calculations with a set of constitutive parameters equal to the reference ones, and finally by analyzing the results in such a way that various criteria were satisfied in the best possible way. These criteria are that stress and strain rate heterogeneities are maximized, gradients in stress or strain near sample edges are minimized, and buckling is prevented. The authors processed FE simulated data. They found a significant difference between the identified and reference parameters, because of uniqueness issues. This problem was fixed by either reducing the number of unknown parameters or, more satisfactorily, by considering a second test in addition to the first one. This second test was performed by applying a displacement to the moving grip with a higher velocity, which increased the richness of the simulated data input into the identification process. No experimental validation has been published so far.

Figure 4: Specimen geometry proposed in [START_REF] Jones | Parameter covariance and non-uniqueness in material model calibration using the virtual fields method[END_REF] to identify visco-plastic parameters.

Optimization-led design based on strain state

To the best of the authors' knowledge, the first published application of this approach concerns the optimization of the dimensions of a T-shaped specimen for orthotropic composite materials [START_REF] Grédiac | A T-shaped specimen for the direct characterization of orthotropic materials[END_REF]. The idea here was to combine a bending test on short beam and a tensile test. The former activates the Young's modulus in the horizontal direction as well as the in-plane shear modulus, the latter the Young's modulus along the vertical direction as well as one of the Poisson's ratios, thus making it possible to identify the four parameters of an orthotropic material on a single test specimen. In [START_REF] Grédiac | A T-shaped specimen for the direct characterization of orthotropic materials[END_REF], the dimensions of the specimen a, b, c, d

(see Figure 5) were adjusted in such a way that the strain levels over the different zones were balanced. This was achieved by minimizing the following cost function with respect to the design variables: after [START_REF] Grédiac | A T-shaped specimen for the direct characterization of orthotropic materials[END_REF].

f ( b a , c a , d a ) = (ε yy -ε xx ) 2 (ε yy + ε xx ) + (ε yy -2ε xy ) 2 (ε yy + 2ε xy ) + (ε yy -2ε xy ) 2 (ε yy + 2ε xy ) (1)
with a acting as a scaling parameter. ε xx is the average of ε xx over area 1, ε yy the average of ε yy over area 2 and ε yy the average of ε yy over area 3, see Figure 5.

The optimization was performed directly in ANSYS 5.1 using their native APDL language.

Unsurprisingly, it was shown that optimal dimensions differed depending on the material orthotropy ratios. A further study proposed an experimental implementation of this test [START_REF] Grédiac | Novel procedure for complete in-plane composite characterization using a T-shaped specimen[END_REF] where elastic constants obtained from actual T-shaped tests were in good agreement with reference parameters obtained with standard homogeneous tests. The grid method was used for the strain maps, and the identification procedure was the VFM with manually-defined virtual fields.

The most advanced implementation of a strain-based test design arises from a series of publications by a Portuguese-French team [108, [START_REF] Souto | Mechanical design of a heterogeneous test for material parameters identification[END_REF][START_REF] Aquino | Design of heterogeneous mechanical tests: Numerical methodology and experimental validation[END_REF][START_REF] Martins | Calibration of anisotropic plasticity models with an optimized heterogeneous test and the virtual fields method[END_REF]. This concerns the identification of the parameters governing the anisotropic Yld2004-18p [START_REF] Barlat | Linear transfomation-based anisotropic yield functions[END_REF] plasticity model. The test was designed by maximizing five quantities defined in [108], which reflect some properties that the specimen should ideally exhibit. The first one is that the strain state range within the specimen should be the widest possible. Each strain state was described by the ratio ε 2 ε 1 between the minor (ε 2 ) and major (ε 1 ) strains of an in-plane stress state. Indeed this ratio represents different strain states (see Figure 6-a), and the widest possible difference between the maximum and minimum values of this ratio represents a wanted feature. This criterion is reflected by C 1 in Equation 2 below.

                                         C 1 = ε 2 ε 1 R = ε 2 ε 1 max - ε 2 ε 1 min C 2 = Std( ε 2 ε 1 ) = n i ε 2 ε 1 i -µ( ε 2 ε 1 ) 2 n -1 C 3 = Std(ε p ) = n i (ε p i -µ(ε p )) 2 n -1 C 4 = ε p M ax = ε p test + ε p tens + ε p shear + ε p biaxial + ε p comp + ε p plane 6 C 5 = Av(ε p ) = n i ε p i v i v t . (2)
A second feature is that the ε 2 ε 1 ratio should be spread out over a large range of values between the extrema for this ratio shown in Figure 6-a, and not be concentrated around the mean value. This property is reflected by the standard deviation of the ratio ε 2 ε 1 throughout the specimen, which gives C 2 in Equation 2, where n is the number of pixels.

The heterogeneity of the plastic strain distribution should also be as wide as possible, which means that the standard deviation of the equivalent plastic strain ε p defined by C 3 in Equation 2 should be maximized.

Since the specimen was designed for sheet metal forming processing, a similar deformation level as the one recorded in actual sheet metal forming should be reached during the test.

Hence the authors proposed to consider a global maximum plastic strain defined by the mean value of the maximum plastic strain achieved during the test on the one hand, and the maximum values of the most relevant strain states expected to occur during the test on the other hand, namely, equibiaxial tension, pure shear, uniaxial tension, plane strain tension and uniaxial compression. This gives quantity C 4 in Equation 2.

Finally, the average plastic strain reached within the specimen was also considered because it reflects the global plastic strain level. The points characterized by the highest value are expected to contain the richest information. This quantity is defined by C 5 in Equation 2.

The global criterion was a scaled and weighted sum of the five quantities C i , i = 1 • • • 5, thus

I t = 5 1 w ri C i w ai (3)
where w ri are scaling parameters and w ai are weighting parameters.

This indicator was then used to rank several test configurations from the literature.

In [START_REF] Souto | Mechanical design of a heterogeneous test for material parameters identification[END_REF], the process was taken one step further by designing a specimen through the minimization of a cost function built up with this indicator. The specimen edges were defined by control points of a spline function parameterizing the shape, and minimization was performed with respect to the coordinates of these points, using a Nelder-Mead direct search algorithm.

In the same spirit, the authors proposed to complete this approach by also designing the shape of the grips through which the load was to be applied to the specimen [START_REF] Aquino | Design of heterogeneous mechanical tests: Numerical methodology and experimental validation[END_REF]. This procedure led to the so-called butterfly specimen (see Figure 6-b), which has been validated

experimentally in [START_REF] Aquino | Design of heterogeneous mechanical tests: Numerical methodology and experimental validation[END_REF] with DIC and FEMU (force and strain cost function) on the Yld2004-18p plasticity model with a Voce hardening function, leading to 18 parameters to identify. In [START_REF] Martins | Calibration of anisotropic plasticity models with an optimized heterogeneous test and the virtual fields method[END_REF], the VFM was used with the same specimen to identify the Hill48 model with a Swift hardening law.

The same group [START_REF] Almeida | On the development of a heterogeneous mechanical test specimen using topology optimization[END_REF][START_REF] Barroqueiro | Design of mechanical heterogeneous specimens using topology optimization[END_REF] recently extended this approach with a topology optimization tool known as Solid Isotropic Material Penalization (SIMP). As opposed to the approach in [START_REF] Souto | Mechanical design of a heterogeneous test for material parameters identification[END_REF] where the external geometry was controlled by splines, this starts with a symmetrical rectangular shape, and uses a density parameter to give more or less weight to different parts of the sample, with the aim of achieving a binary distribution of density, either 0 or 1, defining the shape of the specimen. This requires a scheme to avoid unrealistic solutions in the form of high frequency variations of the density. Some interesting solutions were found.

Design by identification quality

The procedures described above only focus on strain distributions. While they clearly influence the quality of identification, it is difficult to establish a quantitative link between strain distributions and identified quantities without performing the identification. This is why the researchers have been trying to use metrics on the quality of identification as a target for test design. In this case, the cost functions driving the test optimization involve identification errors or uncertainties for VFM-based identifications, or FE parameter sensitivities for FEMU identification. The price to pay is that the identification procedure is called in each loop of the iterative design procedure. This type of optimization is therefore much more computationally intensive compared to the use of just strains. Therefore, the computation time of the identification technique becomes a key parameter in this case.

Sensitivities to constitutive parameters in FEMU

A natural way to ensure that constitutive parameters are well activated in a particular test is to compute maps of sensitivities of displacement or strain fields to each constitutive parameter. This is done in practice by varying each parameter individually and deriving the sensitivity field by finite difference. This requires at least a number of FE computations equal to twice the number of parameters for a single sensitivity evaluation. It should be noted that in IDIC, the sensitivities are an integral part of the identification process.

To the best knowledge of the authors, the first attempt at test optimization using such sensitivities can be found in [START_REF] Magorou | Identification of constitutive laws for wood-based panels by means of an inverse method[END_REF]. The authors considered bending tests on square wood panels. The location of the supports and applied loads was optimized by maximizing the sensitivity of the deflection of the specimen at the nodes of a regular mesh with respect each unknown bending rigidity, deflection being estimated with finite element computations. The sensitivity S(D ij ) of each bending stiffness component D ij is defined by:

S(D ij ) = || dU dD ij ||/ || U D ij ||, ij = xx, yy, xy, ss. (4)
where ||X|| denotes the L 2 norm of vector X. D xx , D yy , D xy and D ss are the four orthotropic bending stiffness components. Maximization was performed by using a genetic algorithm. The location of these points being found, an experiment was performed on an Okoume plywood panel. The authors reported that the stiffnesses identified with this procedure were in good agreement with theoretical expectations. The measuring technique for the deflection field was fringe projection, and the identified method FEMU with a cost function based on deflections.

More recently, Bertin et al. [START_REF] Bertin | Optimization of a cruciform specimen geometry for the identification of constitutive parameters based upon full-field measurements[END_REF] used sensitivity maps to optimize the fillet radii of a biaxial cruciform specimen for elastic and elasto-plastic identification. It should be noted that for elasticity, a scaling procedure was introduced, which here was the maximum admissible von

Mises stress to ensure that the specimen remained fully elastic. Indeed, different geometries exhibit different global stiffnesses and the sensitivities may be influenced by that. This need for scaling is a recurring issue to optimize tests for elasticity, and will be encountered in several studies reported later in this article.

The work presented in [START_REF] Chamoin | Coupling between topology optimization and digital image correlation for the design of specimen dedicated to selected material parameters identification[END_REF] considers a more advanced shape optimization methodology, Solid Isotropic Material Penalization (SIMP) but still based on FE sensitivities. It was proposed to use IDIC [START_REF] Mathieu | Estimation of elastoplastic parameters via weighted FEMU and integrated-DIC[END_REF] to identify the shear modulus of linear elastic orthotropic materials.

The initial specimen shape was rectangular and the load was a uniaxial compression. The specimen was divided into small elements, each of them being affected by a density ρ ranging between 1 and a scalar α, which is positive and close to 0 (and not rigorously 0 to keep the problem well-posed). 0 reflects the total absence of material and 1 its presence. The authors defined the geometry of the specimen by considering the sensitivity of the displacement field to the spatial distribution of ρ. The cost function to be minimized involved both the sensitivity to the shear modulus and the sensitivity to the spatial distribution of ρ. At convergence, a value of ρ close to zero means that the corresponding element should be removed, a value close to one that it should be kept. Numerical simulations carried with three completely different initial spatial distributions of ρ showed that similar final patterns were obtained at convergence in the three cases. A relatively marked contrast in the distribution of ρ then facilitated the definition of a threshold value beyond which the small ρ-value elements were removed. Only the shear modulus was considered as unknown, the three other constitutive parameters being given a priori. Experiments carried out on an optimized orthotropic wood specimen showed that the shear modulus could be identified with this approach, the sharpness of the cost function around the identified value illustrating the robustness of the procedure.

Noise sensitivities in the VFM

In the VFM, FE sensitivities are less relevant as the procedure does not involve any FE computation. In elasticity, the relevant quality indicators are the noise sensitivity parameters η ij , which relate the standard deviation of the strain noise to the standard deviation of each Q ij stiffness parameter as in [START_REF] Avril | Sensitivity of the virtual fields method to noisy data[END_REF].

Std(Q ij ) = η ij γ, ij = xx, yy, xy, ss. ( 5
)
where γ is the standard deviation of the strain noise, assumed here to be Gaussian, white and uncorrelated between the three strain components. The latter assumption is clearly not rigorous but it is necessary to solve the problem of virtual field optimization; even if the virtual fields are not optimal stricto sensu, they are still much better than any manually-defined ones.

The first study using these noise sensitivity parameters is reported in [START_REF] Pierron | Identification of the orthotropic elastic stiffnesses of composites with the Virtual Fields Method: sensitivity study and experimental validation[END_REF], where it was proposed to optimize the unnotched Iosipescu specimen discussed in Section 3 with respect to two design variables, the length of the gauge section L and the orientation of the fibres in the specimen, θ, see Figure 7-a. The cost function to be minimized is such that the sensitivities η ij are balanced. This cost function is defined by:

f (L, θ) = (η xx -η yy ) 2 + (η xx -η ss ) 2 -(η yy -η ss ) 2 (η xx + η yy + η ss ) 2 (6)
η xy was not considered here because its value is always much larger than the others and would therefore have overwhelmed the cost function. The idea of using a normalized cost function is to avoid the problem of scaling. Indeed, since the problem is fundamentally a signalto-noise ratio one, the procedure would naturally converge towards the stiffest configuration when displacements are imposed, thus leading to the largest possible strain values, regardless of whether they are feasible or not (onset of non-linearity, fracture, poor aspect ratio). In [START_REF] Rossi | Effect of DIC spatial resolution, noise and interpolation error on identification results with the VFM[END_REF], this problem was solved by using a maximum stress failure criterion, playing the same role as the von Mises stress in [START_REF] Bertin | Optimization of a cruciform specimen geometry for the identification of constitutive parameters based upon full-field measurements[END_REF]. But here, a scaled cost function was used to mitigate this issue for the sake of simplicity.

Results show that the sensitivity is much higher for θ than for L, as illustrated in Figure 5b where the cost function is nearly flat along the horizontal direction. L was therefore fixed at a reasonable value and θ selected according to the cost function. Experimental results are presented in [START_REF] Pierron | Identification of the orthotropic elastic stiffnesses of composites with the Virtual Fields Method: sensitivity study and experimental validation[END_REF]. The grid method was used as strain measurement technique.

(a) Schematic view of specimen and fixture, after [START_REF] Pierron | Identification of the orthotropic elastic stiffnesses of composites with the Virtual Fields Method: sensitivity study and experimental validation[END_REF] (b) Cost function and two design points, after [START_REF] Pierron | Identification of the orthotropic elastic stiffnesses of composites with the Virtual Fields Method: sensitivity study and experimental validation[END_REF] Figure 7: Design of an unnotched Iosipescu specimen, after [START_REF] Pierron | Identification of the orthotropic elastic stiffnesses of composites with the Virtual Fields Method: sensitivity study and experimental validation[END_REF].

One of the main limitations of this approach is that it does not take into account the resolution of the camera. In fact, preliminary optimizations tended to converge towards very elongated specimens, a trend that can be seen in Figure 7-b. This forced the authors to put an upper limit on the specimen length in the optimization process. This is one of the motivations to move to image deformation as presented in the last section.

Other cases of optimal test design based on this approach have been proposed in the literature. For instance, [START_REF] Wang | Identification of material parameters of PVC foams using digital image correlation and the Virtual Fields Method[END_REF] addressed the design of a test for orthotropic foam specimens.

The specimen was placed in an Arcan fixture, and both the angle of the fixture with respect to the vertical load and the orientation of the specimen in this fixture were optimized. Data are often missing near the free edges of a test specimen, so the sensitivity of the identified parameters to these missing data was also considered in the optimization procedure. The cost function was similar to that in Eq. 6 [START_REF] Pierron | Identification of the orthotropic elastic stiffnesses of composites with the Virtual Fields Method: sensitivity study and experimental validation[END_REF], which meant that no force scaling was necessary as it was embedded in the cost function. It was found that the four stiffnesses (and not only the shear modulus) could be successfully identified with this test when optimized parameters were considered.

In the same spirit, the VFM sensitivity to noise was employed in [START_REF] Muhammad | Characterization of composite plates using the virtual fields method using optimized loading conditions[END_REF] to design a bending test where the six bending stiffnesses of a composite plate were the unknowns. The shape of the plate was parameterized with a spline function. The design variables were the coordinates of the control points of this function, the coordinates of the three supports of the plate and those of the point where the force was to be applied. This led to some non intuitive shape for the specimen, as well as non-intuitive locations of the supports and applied force, as illustrated in Figure 8. No experimental validation was carried out. This example also shows the limitations of this approach which does not take into account the measurement process.

The shape in Figure 8 will lead to the loss of many pixels of the camera, thus challenging the spatial resolution of the measurements.

Figure 8: Optimization of a bending test for composite plates, after [START_REF] Syed-Muhammad | Optimization of a mechanical test on composite plates with the virtual fields method[END_REF].

Finally, in [START_REF] Zhang | Experimental identification of statistical correlation between mechanical properties of FRP composite[END_REF], it was proposed to optimize an open-hole off-axis specimen subjected to a tensile test in order to identify the four elastic parameters of a unidirectional orthotropic carbon-epoxy composite. The statistical correlation existing between these parameters was also studied. The cost function to be minimized was the sum of the relative difference between reference and identified parameters, and the design variables were the fibre angle and the diameter of the hole. DIC was used to measure the strain field, and the VFM with special optimized virtual fields [START_REF] Avril | Sensitivity of the virtual fields method to noisy data[END_REF]. The optimal configuration they identified was with the fibres aligned with the specimen length and a hole of 7 mm for a width of 36 mm and a length of

In the examples of the preceding section, the data used in the optimization procedure were directly provided by finite element models. However, these data are not exactly equal to the displacement or strain measurements that would be returned by a white-light full-field measurement system such as DIC. Indeed the metrological performance of full-field measurement techniques is now better known, and various types of random and systematic errors impairing the measurements have been characterized in many recent papers. In particular, a single pixel cannot provide displacement information; it is necessary to rely on the information at the pixels located in a zone surrounding this point. In classical DIC, this zone is a subset, which shape is generally quadrilateral [START_REF] Sutton | Image Correlation for Shape, Motion and Deformation Measurements[END_REF]. When processing regular patterns, various strategies can be employed [START_REF] Grédiac | Comparing several spectral methods used to extract displacement and strain fields from checkerboard images[END_REF] but the most popular one is based on the windowed Fourier transform restricted to a single frequency equal to the nominal frequency of the regular pattern, and the window may have various shapes [START_REF] Sur | Influence of the analysis window on the metrological performance of the grid method[END_REF]. With both techniques, the consequence is that the displacement returned at a given point also depends on the displacement returned at points located in its vicinity, and the wider the subset or the window, the more marked this dependence. This introduces a bias in the displacement and strain maps, which manifests itself by a blur, thus by the fact that the amplitude of the details and strain gradients are attenuated, as explained for instance in [START_REF] Schreier | Systematic errors in digital image correlation due to undermatched subset shape functions[END_REF][START_REF] Sur | On biases in displacement estimation for image registration, with a focus on photomechanics[END_REF] for DIC, and [START_REF] Sur | Towards deconvolution to enhance the grid method for in-plane strain measurement[END_REF] for the GM. It also means that considering directly the results of finite element simulations as synthetic experimental data may potentially lead to misleading conclusions when designing optimal testing configurations.

For instance, considering these simulations alone may lead to an optimal solution featuring high strain gradients containing rich information, but such gradients may not be correctly rendered by an actual full-field measurement system. Reducing as much as possible the size of the subset for DIC or the window for the GM limits this bias, but the random error due to the fact that images are unavoidably noisy increases in proportion. A trade-off has therefore to be found, and various parameters can be adjusted for that in addition to the size discussed above, such as the degree of the matching function used to model the actual displacement field within the subsets and the step (defined by the distance between two remote subsets)

for DIC. The best trade-off should ideally be obtained by adding these settings to the design variables.

Before optimizing all these variables at the same time, intermediate steps have been proposed in the literature. The first one was to simulate the whole chain going from the images to the identified parameters in order to study the influence of various design variables. The first attempt was published in 2012 [START_REF] Rossi | On the use of simulated experiments in designing tests for material characterization from full-field measurements[END_REF] where synthetic grid images were considered. Displacement and strain fields were extracted from these images, and then considered as input data for the VFM. This study was the direct follow-up of that presented in the previous section, [START_REF] Pierron | Identification of the orthotropic elastic stiffnesses of composites with the Virtual Fields Method: sensitivity study and experimental validation[END_REF], where image deformation was not taken into account. The results here showed that a natural minimum appeared in the cost function without having to restrict the range of specimen length. As the procedure naturally takes into account the resolution of the camera, high aspect ratio specimens were penalized and a length of 30 mm was identified as optimal for the unnotched Iosipescu specimen. Since measured displacements are noisy, it is generally necessary to apply some spatial or temporal smoothing before proceeding to differentiation to obtain strains. Such smoothing acts as an additional low pass filter in space and it is necessary to consider how this operation can be performed in an optimal manner to minimize the bias arising from it. Such optimal displacement smoothing parameters can be selected thanks to the simulator introduced in [START_REF] Rossi | On the use of simulated experiments in designing tests for material characterization from full-field measurements[END_REF], and uncertainty bounds for the identified stiffness components can be derived from the process. This simulator was then extended to DIC [START_REF] Rossi | Effect of DIC spatial resolution, noise and interpolation error on identification results with the VFM[END_REF],

unsurprisingly showing the superiority of the performances of the grid method, consistent with recent results obtained with a rigorous approach [START_REF] Grédiac | A critical comparison of some metrological parameters characterizing local digital image correlation and grid method[END_REF]. This simulator with the Grid Method has also been used in transient dynamics to optimize post-processing parameters for the so-called Image-Based Inertial Impact (IBII) test [START_REF] Fletcher | An image-based inertial impact (IBII) test for tungsten carbide cermets[END_REF]. It was used in [START_REF] Bouda | A computational approach to design new tests for viscoplasticity characterization at high strain-rates[END_REF] to compare three IBII test configurations (a plain one, and two others with a hole or notches) for viscoplastic identification of a titanium alloy. The effect of impact speed was also explored, together with spatial and smoothing parameters.

All the links of the chain spanning between synthetic images and final constitutive parameters being available through suitable models, the next stage is to optimize the settings and design variables governing each of these links. This is however quite a complex problem because of i-the number of parameters to be considered, which is greater than in Section 5 above, and ii-the fact that some variables such as the subset size for DIC and the degree of the matching function are discrete, which makes it necessary to use more involved optimization procedures like genetic algorithms. This obviously induces significant computation times, the identification procedure being run within each iteration of the optimization procedure.

A first attempt was proposed in [START_REF] Rossi | A procedure for specimen optimization applied to material testing in plasticity with the virtual fields method[END_REF] for elasto-plastic identification, where the geometry of a double notch machined in a tensile specimen was optimized. Two subset sizes and 25 different combinations of variables defining the notches were considered. The cost function used here was the squared difference between the reference and identified stress-strain curves.

The reason for this choice is that the stress-strain curve can be more or less sensitive to some parameters, so only the most influencing ones were considered if choosing this cost function.

The best design was the one minimizing the distance between these two stress-strain curves.

Interestingly, the best option was shown to be the one leading the a diffuse plastic strain distribution while non-optimal ones exhibited high strain concentration close to the notches, which was probably due to the fact that high strain gradients were poorly rendered by the DIC system. As in [START_REF] Bouda | A computational approach to design new tests for viscoplasticity characterization at high strain-rates[END_REF], this approach was however only a comparison between potential candidates defined a priori, so the procedure probably missed some better configurations.

A more sophisticated approach was presented in [START_REF] Gu | Towards the design of a new standard for composite stiffness identification[END_REF] to tackle this problem by using a two-step approach. Indeed the geometry and the load were first optimized by using fixed (but reasonable) values for the DIC settings. The first set of design variables defining the geometry and load were then fixed after this first step, and DIC parameters were optimized in a second step. Different tests were considered to identify the parameters governing the elastic behaviour of orthotropic materials like composites, namely an off-axis tensile test on short specimen (SOAT), a Brazilian test performed on off-axis disks (OABD) and a Iosipescu test on an unnotched off-axis specimen (OAUIT). A hole was also added at the centre of the first type of specimen to examine the benefit of an increased heterogeneity due to this strain concentrator (SOAHT). These different configurations are shown in Figure 9-a. Compared to the open-hole tensile specimen proposed in [START_REF] Lecompte | Identification of elastic orthotropic material parameters based on ESPI measurements[END_REF][START_REF] Silva | Identification of material properties using FEMU: Application to the open hole tensile test[END_REF][START_REF] Kowalczyk | Identification of mechanical parameters of composites in tensile tests using mixed numerical-experimental method[END_REF][START_REF] Gogu | Bayesian identification of elastic constants in multi-directional laminate from moiré interferometry displacement fields[END_REF] and discussed in Section 3 above, the aspect ratio of the specimen considered here is lower in order to increase heterogeneity, even in the absence of hole. The orientation of the specimen and off-axis angle were first optimized in addition to the diameter of the hole. The C 3 cost function below was used for this optimization, combining systematic (C 1 ) and random (C 2 ) errors in C 3 (maximal error within a 95% confidence interval):

                   C 1 = 1 4 ij |Q id ij -Q ref ij | Q ref ij C 2 = 1 4 ij σ id ij Q ref ij
, ij = 11, 22, 12, 66

C 3 = C 1 + 2C 2 (7)
where id denotes an identified stiffness and ref the reference one. σ id ij is the standard deviation of Q id ij estimated from images with 30 different copies of a white Gaussian grey level noise. The same maximum stress criterion scaling strategy as in [START_REF] Rossi | Effect of DIC spatial resolution, noise and interpolation error on identification results with the VFM[END_REF] was used. Figures 9(b)-(e) show C 3 as a function of the design variables chosen for each configuration. The main conclusion is that a unique minimum exists in each case, and the variables corresponding to this minimum are those retained to define the optimal test design.

OAUIT was found to be the best candidate after the first optimization step, closely followed by SOAHT. A second step was then designed to optimize DIC processing parameters (subset size and strain window, with step size fixed at 50% of the subset size) to try and reduce the values of the C 3 cost function. Interestingly, the results from the second step change this ranking with SOAHT stepping ahead of OAUIT. This result shows that constraints coming from the measuring system should be taken into account during the design process of optimal tests. Another interesting conclusion of this second step is that quadratic matching functions (instead of bilinear ones) and steps of one pixel between successive subsets should systematically be used, leading to lower values of C 3 . The superiority the quadratic shape functions has been documented in previous work [START_REF] Schreier | Systematic errors in digital image correlation due to undermatched subset shape functions[END_REF][START_REF] Blaysat | Towards criteria characterizing the metrological performance of full-field measurement techniques. application to the comparison between local and global versions of DIC[END_REF], while step size has, regrettably, not been the focus of much attention in the past. It is worth noting that the orientation of the fibres in OAUIT has also been optimized with this type of approach in the particular case of wood, see [START_REF] Kretschmann | Optimized test design for identification of the variation of elastic stiffness properties of loblolly pine (pinus taeda) pith to bark[END_REF].

In the same spirit as in the preceding study, a refinement of the work presented in [START_REF] Wang | Identification of material parameters of PVC foams using digital image correlation and the Virtual Fields Method[END_REF] on orthotropic polymeric foams (and discussed in the preceding section) tested with an Arcan fixture was proposed in [START_REF] Wang | Optimised experimental characterisation of polymeric foam material using DIC and the Virtual Fields Method[END_REF]. The experimental device is shown in Figure 10. The specimen was bonded onto the surfaces of the Arcan rig which was in turn mounted onto a uniaxial test machine. Optimization was conducted again in two steps. DIC was the measurement method, so reasonable settings were fixed for a first optimization of both the angle of the fixture with respect to the vertical load (α) and the orientation of the specimen in the fixture (θ). Synthetic reference and deformed speckle images were used to obtain displacement and strain maps. These maps were then fed into the optimisation procedure. Force scaling was performed using a maximum strain criterion. The second step consisted in considering the optimized configuration defined during the first step, and then in optimizing the DIC settings.

This led to an optimal configuration which was different from the one given in [START_REF] Wang | Identification of material parameters of PVC foams using digital image correlation and the Virtual Fields Method[END_REF]. The stiffnesses obtained with this approach were similar but the uncertainty with which they were provided could be predicted only with this more refined approach. Experimental results showed that the standard deviation reflecting this uncertainty were predicted within less than a factor two, which can be considered as satisfactory since other sources of errors than the camera sensor exist when actual experiments are performed, such as microvibrations, light fluctuation, geometrical defects, defects in load introduction and alignment etc.

Figure 10: Arcan fixture used to test an orthotropic foam specimen, after [START_REF] Wang | Identification of material parameters of PVC foams using digital image correlation and the Virtual Fields Method[END_REF]. Both the orientation of the fixture (angle α) and the off-axis angle of the specimen, θ, were optimized, as well as the DIC and smoothing parameters used to extract strain fields from the speckle images by DIC.

In the examples discussed above, the two-step approach was used only in order to simplify the problem, but it is clear that all the variables should be considered at the same time while optimizing them so that any possible coupling is not missed. A preliminary study has been carried out in [START_REF] Gu | Full optimization of the T-shaped tensile test using genetic algorithm[END_REF] to reach this goal. The design of the T-shaped specimen addressed in Section 4 was revisited. The total height and width of the specimen were fixed to correspond to the camera sensor aspect ratio, so that the lowest number of pixels are lost. This was a first approach that would need to be refined by including a buckling criterion to check that the design is not limited by the wide horizontal bar. Experimental validation of this design is still needed.

Figure 11: T-shaped specimens, after [START_REF] Gu | Full optimization of the T-shaped tensile test using genetic algorithm[END_REF], with orthotropy direction at 27 • from the horizontal axis for the optimal T, and 70 • for the sub-optimal one.

Discussion and conclusion

The following remarks can be drawn from the examples shown in the four preceding sections. • While early efforts towards MT2.0 focused on the development of full-field measurements and efficient inverse analysis tools (like FEMU and the VFM) to identify constitutive parameters from measured heterogeneous kinematic fields, gradually, researchers are turning towards the design of optimal tests fully adapted to the philosophy of MT2.0.

The distribution of the contributions over the years shows that this activity is indeed quite recent, see Figure 12. There is only a limited number of attempts before 2000, but this research field is ramping up. It can be seen that designing by intuition is by far the subject of the largest number of studies. This is logical since this is the first step towards the more sophisticated approaches which have progressively appeared over the years. Interestingly, even recent contributions belong to this first category, certainly because of the arrival of newcomers to this field fuelled by the increased interest from the scientific community in MT2.0. • Most of the examples deal with anisotropic elasticity, and a few with hyperelasticity and elasto-plasticity, which makes room for future improvement in the last two cases. This is especially true for elasto-plasticity, where advanced constitutive laws (like the Yld2004-18p anisotropic plasticity model in [START_REF] Barlat | Linear transfomation-based anisotropic yield functions[END_REF]) involve an increasingly larger number of parameters.

• Most of the cases discussed above concern the VFM, and this is particularly true when considering the most sophisticated procedures presented in the last cases, where FEMU was used in one case only. This is certainly due to the difference in computational efficiency between FEMU and the VFM. It could be argued that IDIC optimization based on sensitivity maps may be equivalent to the VFM identification simulator as IDIC may automatically incorporate the effects of the systematic errors caused by the lack of camera spatial resolution; though the errors arising from the interplay between the FE and DIC meshes is not entirely clear when these meshes are different.

• The nature of the cost function used to optimize the design parameters (and DIC or GM parameters in Section 6) varies from one example to another. A study solely devoted to the influence of the nature of the cost function on the quality of the results and on the convergence speed should be undertaken, along similar lines to those discussed for the minimization of the optical residual in DIC [START_REF] Tong | An evaluation of digital image correlation criteria for strain mapping applications[END_REF].

• The last three types of design procedures rely on finite element simulations. Such simulations are performed for a given constitutive law and for a given set of constitutive parameters. It means that the optimal solutions they provide are theoretically optimal for these parameters only, and that changing either the type of constitutive law and/or its constitutive parameters should lead to other solutions, thus limiting the scope of any optimal test obtained for certain parameters and not others. For instance, one particular test can be optimal for glass/epoxy composites used during the simulations, but keeping the same geometry and load with carbon/epoxy composites is probably no longer optimal. This advocates for a systematic exploration of the design space for mechanical tests before new materials are tested if their properties are too different from previously considered ones. This could be performed through a commercialized software platform dedicated to test design and uncertainty quantification. Some efforts towards this goal are currently underway in a few teams worldwide but there is no doubt that such tools will be readily available in the near future. This is the adaptation of the notion of digital virtual twin to the MT2.0 paradigm, where each test has its virtual numerical twin used to assess performances and provide uncertainty.

• It is worth noting that only kinematic (displacement, acceleration, slope or strain) measurements are involved in the examples discussed in the four sections above. Thermal fields are also widely used in experimental mechanics, either alone (for instance when performing Thermoelastic Stress Analysis or estimating heat sources at the surface of loaded specimens) or in combination with displacement and strain measurements. Designing tests to exhibit simultaneous kinematic and thermal heterogeneities in order to enrich the experimental database is certainly a promising route.

• Finally, the problem of the assessment of quality of the constitutive model is the large gaping hole in the current state of the art, as also nicely pointed out in [START_REF] Roux | Optimal procedure for the identification of constitutive parameters from experimentally measured displacement fields[END_REF]. Novel nonparametric approaches are emerging [START_REF] Dalémat | Measuring stress field without constitutive equation[END_REF][START_REF] Leygue | Data-based derivation of material response[END_REF][START_REF] Leygue | Non-parametric material state field extraction from full field measurements[END_REF][START_REF] Réthoré | Computational measurements of stress fields from digital images[END_REF], piecewise-defined models have been explored [START_REF] Rossi | A general linear method to evaluate the hardening behaviour of metals at large strain with full-field measurements[END_REF] but in spite of these promising efforts, there is a need for ways of either discriminating between models or formulating models based on MT2.0 methodologies.

This is still an open problem.

In conclusion, it can be said that designing mechanical tests for MT2.0 is an emerging research topic, which lies at the frontier between image processing, experimental and computational mechanics. This is confirmed by the increasing number of recent studies on this topic published in the literature, which are briefly reviewed in the present paper. The examples discussed here show that a certain graduation can be identified, with studies ranging between a mere adaptation of pre-existing classical tests to much more sophisticated ones, where even constraints imposed by the measurement technique are accounted for. The price to pay is the increasing computational resources needed when going from the former to the latter. Recent advances on inverse techniques devoted to identification from full-field measurements on the one hand, and on the metrological performance of full-field measurement techniques on the other hand, now forge a solid footing for MT2.0, which means that optimizing mechanical tests dedicated to this type of approach should attract more researchers in the near future and will hopefully penetrate industry, and lead to new standards in the future.

  ure 2-b. Simulations confirmed the coexistence of these three strain states at the centre of the specimen, see Figure 2-c. According to [101], processing the test images with DIC and identifying the parameters with the VFM provided constitutive parameters in agreement with those obtained with classical homogeneous tests.(a) Model of the specimen (b) Testing device (c) Three strain states in the specimen

Figure 2 :

 2 Figure 2: Test design for elastomeric materials where uniaxial tensile (UT), pure shear (PS) and equibiaxial tensile (ET) strain coexist, after [101].

Figure 5 :

 5 Figure 5: Schematic view of the T-shaped specimen (left), and areas for strain balancing (right) after [106].

( a )

 a Desired strain states to occur in the specimen, after [108] (b) Butterfly specimen, after[START_REF] Aquino | Design of heterogeneous mechanical tests: Numerical methodology and experimental validation[END_REF]

Figure 6 :

 6 Figure 6: Designing a specimen for elasto-plasticity, after [108, 110].

Figure 9 :

 9 Figure 9: Four tests for composite characterization and cost function C 3 as a function of various of design variables, after [135].

  Five design variables were considered to define the geometry of the specimen and the orientation of the material: three geometrical parameters b and c as defined in Figure5, the radius of the fillet joining the horizontal and vertical bars of the 'T', the orthotropy angle, the DIC settings, namely the subset size, the step size, and the degree of the matching functions (linear or quadratic), and the size of the Gaussian filter used to smooth the displacement field prior to differentiation to obtain the strain components. A genetic algorithm was used to minimize cost function C 3 defined above because variables 5 to 7 are discrete while the others are continuous. The maximum stress was again used for scaling purposes.The optimal shape is shown in Figure11. It is very interesting to note how the initial idea of the T-shaped specimen (horizontal and shear moduli active in the horizontal bar, and transverse modulus and Poisson's ratio in the vertical bar) has been totally lost in the converged design. Indeed, the identifiability of all orthotropic stiffness components is now ensured by the 'flow' of a tensile stress applied by the short vertical bar into what can be seen as a semi-infinite plane, with the off-axis angle ensuring the richness of the test. Compared to the bad configuration on the right hand-side of the figure, one can first see how the optimization tends to ensure maximal use of the camera pixels to maximize spatial resolution, and avoid high strain concentrations at the corners of the T. It is also worth mentioning that the algorithm systematically converged towards a step size of 1 and quadratic shape functions.

Figure 12 :

 12 Figure 12: Distribution over time of the studies dealing with optimization of tests for Material Testing 2.0, which are cited in the present paper.

•

  With the maturation of camera-based full-field deformation measurements (Digital Image Correlation, the Grid Method), a paradigm shift in the mechanical testing of ma-

	terials is emerging. Traditional test methods developed for point sensors like strain
	gauges or extensometers have largely relied on statically determinate configurations, for
	which a closed-form expression linking the unknown constitutive parameters and local
	measurements exists. Unfortunately, these tests exhibit strong limitations like the re-
	liance on well-controlled boundary conditions, restricted stress-strain information and
	difficulties to deal with heterogeneous materials and complex constitutive models that
	require many tests to be calibrated. In contrast, the wealth of information contained
	in full-field (or more rigorously, spatially-dense) measurements opens up opportunities
	to design new statically indeterminate tests for more efficient and cost-effective model
	calibration. It was chosen here to give visibility to this new paradigm by christening it
	'Materials Testing 2.0' (MT2.0).