Deformation and Failure of MXene Nanosheets

This work is aimed at the development of finite element models and prediction of the mechanical behavior of MXene nanosheets. Using LS-Dyna Explicit software, a finite element model was designed to simulate the nanoindentation process of a two-dimensional MXene Ti3C2Tz monolayer flake and to validate the material model. For the evaluation of the adhesive strength of the free-standing Ti3C2Tz-based film, the model comprised single-layered MXene nanosheets with a specific number of individual flakes, and the reverse engineering method with a curve fitting approach was used. The interlaminar shear strength, in-plane stiffness, and shear energy release rate of MXene film were predicted using this approach. The results of the sensitivity analysis showed that interlaminar shear strength and in-plane stiffness have the largest influence on the mechanical behavior of MXene film under tension, while the shear energy release rate mainly affects the interlaminar damage properties of nanosheets.


Introduction
A new class of two-dimensional (2D) nanomaterials, MXenes, was discovered in the last decade [1]. MXenes are transition metal carbides or nitrides produced by the etching of the A element from the MAX phases. Typically, nanomaterials can be divided into two groups: hydrophilic but not conductive, such as transition metal oxides or clays; and conductive but not hydrophilic, such as graphene. However, some MXenes (Ti 2 CT z , Ti 3 C 2 T z ) have the unique characteristics of both groups. Due to the combination of the electrical conductivity of transition metal carbides and the hydrophilicity of hydroxyl or oxygen-terminated surfaces, these MXenes behave as "conductive clays" [2].
MXenes have been widely investigated during the past few years. The main interest has been directed at the electric properties, their applications for sensors, energy storage (batteries, supercapacitors, hydrogen evolution reaction catalysts), harvesting, electromagnetic shielding, tribology, etc. [3][4][5][6]. Pioneer studies of the mechanical properties of MXenes showed promising results [1,2,7]. Elastic properties were obtained experimentally by nanoindentation with the tip of an atomic force microscope (AFM) and the elastic modulus of the most investigated MXene material, Ti 3 C 2 T z , was obtained at 0.33 ± 0.03 TPa [8]. According to classical molecular dynamics simulation [7], which does not take into account various material defects, the modulus was higher and equal to 0.502 TPa.
MXenes exhibit a high bending stiffness [9]. The critical deformations are much higher than the graphene ones, and this is an important feature of flexible electronics [4,10,11]. MXenes have good interactions with polymeric matrices for polymer composite applications [4,12,13]. For these reasons, MXene could be a good candidate to provide electrical conductivity for fiber-reinforced plastic composites without losing desirable mechanical properties and imparting additional self-sensing composites without losing desirable mechanical properties and imparting additional self-sensing functions. While some of the mechanical properties of MXene have been determined theoretically or experimentally, the mechanical behavior of MXene, and particularly its failure mechanisms, has not been studied well, and there is a huge lack of data that are needed for the development of polymer composites filled with MXene 2D nanosheets.
The aim of this study was to investigate the micromechanical behavior of MXene nanosheets by developing finite element (FE) computational models. The objectives were to (1) develop an FE model composed of single-layered MXenes nanosheets with a specific number of individual flakes for explicit analysis; and (2) identify the material parameters by FE simulation of interface shearing. The novelty of this study is that it shows first insights into the deformation and failure mechanisms of this new nanomaterial, as well as providing a basis for the future design of polymer composites reinforced with MXene nanosheets, and the development of MXene-polymer coatings with highdensity MXene-MXene interactions.

FE Model of Nanoindentation
The main purpose of this presimulation was to build an explicit FE model of a Ti3C2Tz MXene monolayer flake for analysis of the nanoindentation process in LS-Dyna software and, using the force vs. deflection curve, validate the FE model and material characteristics. The FE model ( Figure 1) was developed according to the experimental data presented by Lipatov et al. [8], where it was considered that the MXene flake has isotropic properties, and therefore, the membrane can be parametrized using Young's modulus, E and Poisson's ratio, v. The thickness of the Ti3C2Tz monolayer flake is an important parameter as it influences the results of the nanoindentation experiments. Using AFM for the determination of the thickness of monolayers of 2D materials has some limitations. The Ti3C2Tz MXene flake thicknesses obtained by AFM can differ significantly [8,14], and this directly affects the determined value of Young's modulus. MXene flakes with a thickness of 0.98 nm were modeled [1,2,11]. The mechanical properties of MXenes and the SiO2 support ring used in the model are based on the analysis of the literature data. The indenter was a diamond crystal with a modulus of 1 TPa. The mechanical properties of the materials used in the nanoindentation simulation are presented in Table 1. The thickness of the Ti 3 C 2 T z monolayer flake is an important parameter as it influences the results of the nanoindentation experiments. Using AFM for the determination of the thickness of monolayers of 2D materials has some limitations. The Ti 3 C 2 T z MXene flake thicknesses obtained by AFM can differ significantly [8,14], and this directly affects the determined value of Young's modulus. MXene flakes with a thickness of 0.98 nm were modeled [1,2,11]. The mechanical properties of MXenes and the SiO 2 support ring used in the model are based on the analysis of the literature data. The indenter was a diamond crystal with a modulus of 1 TPa. The mechanical properties of the materials used in the nanoindentation simulation are presented in Table 1. The 2D Ti 3 C 2 T z MXene flake was modeled with shell elements using a linear elastic material model. The size of the shell elements was 10 nm, while the center of the monolayer (contact zone with nanoindenter) was decreased to as little as 1.25 nm. The nanoindenter was defined as an elastic solid sphere (diameter 14 nm). The bottom of the MXene flakes was supported on the ring surface of SiO 2 , which was fully fixed. Two pinball-type contacts (AUTOMATIC_NODES_TO_SURFACE) were used: surface of SiO 2 -bottom flake, and surface of the nanoindenter-top flake. Two simulations were performed with initial pretension and without. As it was an explicit analysis, the initial pretension initiated the oscillation of the flake; therefore, the *CONTROL_DYNAMIC_RELAXATION must be activated.
The indentation was described by the dependence of the displacement on time. By linearly increasing the displacement, which gives a constant velocity, a small impact phenomenon was obtained in the model, and this initiated the vibrations of the MXene monolayer flake. Finally, a smooth increase in the displacement was chosen based on the assumption that the average speed of the nanoindenter is equal to~1 m/s, and the speed increases linearly from 0 up to the 2v_avg. The loading curve displacement vs. time was obtained by the function u(t) = v(t)dt.
A methodology for deriving the material parameters from experimental results, known as parameter identification, was applied here using the optimization procedure. The same nanoindentation FE model without initial pretension has been chosen to evaluate the sensitivity of the material and the geometric parameters to the mechanical behavior of the sample. The Young's modulus and flake thickness were chosen as variables. The Young's modulus varied between 200 and 400 GPa and the thickness from 0.4 to 1.5 nm. The mean square error was used as a curve fitting metric. The experimental force, F vs. deflection δ curve can be described by the following equation [8]: where σ 2D 0 = σ h represents the prestress in the membrane, E 2D = E h is the 2D elastic modulus (thickness h = 0.98 nm), and r is the radius of the well [8]. The dimensionless constant, q, is related to ν as q = 1/ 1.049 − 0.15ν − 0.16ν 2 = 0.9933. The first term in Equation (1) corresponds to the prestretched membrane regime and is valid for small loads. The second term for the nonlinear membrane behavior is characterized by a cubic F~δ 3 relationship with a coefficient of E 2D , which dominates at large loads. For comparison, only the second part of Equation (1) was used; therefore, the prestretching was not taken into account during the simulation.

FE Model of Pure MXene Film
The stability of the stacked two-dimensional transition metal carbides and their interlayered friction in different configurations are comparatively studied by means of density functional theory. In recent years, the adhesive interactions of monolayers and few-layer 2D materials have been intensively investigated. At equilibrium, a larger interlayer distance corresponds to a smaller binding energy, suggesting an easier sliding between the layers [17]. Nanoindentation has been widely used to characterize the adhesion of thin films [18][19][20]. One of the important questions is to understand how the properties-in particular, adhesive strength-change when transitioning from bulk to 2D forms of the material. However, to our knowledge, no theoretical or experimental studies of the adhesive properties of MXenes have been reported to date. The results of direct AFM measurements of adhesion of two MXenes, Ti 3 C 2 T z and Ti 2 CT z , with a SiO 2 -coated Si spherical tip is one of the recent studies of adhesion properties [17].
One of the ways to analyze adhesive strength is to apply the reverse engineering method. Assuming that adhesion energy between surfaces of free-standing MXene nanosheets exists [21], then the strength of the interlayer surface has to depend on the overlapping area. For the study of adhesive interactions, the experimental tensile test data of the assembled free-standing Ti 3 C 2 T z -based films [4] were used. It was assumed that the single-layered nanosheet has a square form of 1 µm length [22][23][24][25], consisting of 18 Ti 3 C 2 T z individual flakes and an average thickness of 20 nm [25]. In total, 164 layers of single-layered nanosheets per film of 3.3 µm thickness were used in experimental testing [4]. As the overlapping of nanosheets has an essential influence on the strength of the interlayer surface, a 2D analysis of randomly placed rectangle nanosheets was performed using materials modeling software, Digimat. For the overlapping analysis, MXene nanosheets with dimensions of 1000 × 20 nm were chosen. The overall thickness of MXene film formed from these nanosheets was 3.3 µm, as it was in testing [4] (Figure 2a). The results show that the average overlapping of nanosheets is 20% (Figure 2b). This overlap value was used to create the FE model. of adhesion of two MXenes, Ti3C2Tz and Ti2CTz, with a SiO2-coated Si spherical tip is one of the recent studies of adhesion properties [17].
One of the ways to analyze adhesive strength is to apply the reverse engineering method. Assuming that adhesion energy between surfaces of free-standing MXene nanosheets exists [21], then the strength of the interlayer surface has to depend on the overlapping area. For the study of adhesive interactions, the experimental tensile test data of the assembled free-standing Ti3C2Tz-based films [4] were used. It was assumed that the single-layered nanosheet has a square form of 1 μm length [22][23][24][25], consisting of 18 Ti3C2Tz individual flakes and an average thickness of 20 nm [25]. In total, 164 layers of single-layered nanosheets per film of 3.3 μm thickness were used in experimental testing [4]. As the overlapping of nanosheets has an essential influence on the strength of the interlayer surface, a 2D analysis of randomly placed rectangle nanosheets was performed using materials modeling software, Digimat. For the overlapping analysis, MXene nanosheets with dimensions of 1000 × 20 nm were chosen. The overall thickness of MXene film formed from these nanosheets was 3.3 μm, as it was in testing [4] (Figure 2a). The results show that the average overlapping of nanosheets is 20% (Figure 2b). This overlap value was used to create the FE model. The FE model was designed taking into account the experimental setup data presented in the study [4]. The single-layered Ti3C2Tz MXene nanosheet was modeled with shell elements. The Young's modulus of the MXene nanosheet was set at 333 GPa; Poisson's ratio-0.227; the size of shell elements-5 nm. The nodes on the left-side edges had a fixed 4 degrees-of-freedom, allowing free contraction in the x-direction and rotation about the z-axis. For the right-side nodes, the displacement vs. time u(t) was applied with a speed of 1 mm/ms. The interface between the single-layered nanosheets was modeled using tiebreak contact, *CONTACT_AUTOMATIC_ONE_WAY_SURFACE_TO_SURFACE_TIEBREAK. The discrete crack model with power-law damage, which works with offset shell elements (option = 11), shown in Figure 4, was chosen. The parameters needed to describe tiebreak contact are presented in Table 2. The FE model was designed taking into account the experimental setup data presented in the study [4]. The single-layered Ti 3 C 2 T z MXene nanosheet was modeled with shell elements. The Young's modulus of the MXene nanosheet was set at 333 GPa; Poisson's ratio-0.227; the size of shell elements-5 nm. The nodes on the left-side edges had a fixed 4 degrees-of-freedom, allowing free contraction in the x-direction and rotation about the z-axis. For the right-side nodes, the displacement vs. time u(t) was applied with a speed of 1 mm/ms. The interface between the single-layered nanosheets was modeled using tiebreak contact, *CONTACT_AUTOMATIC_ONE_WAY_SURFACE_TO_SURFACE_TIEBREAK. The discrete crack model with power-law damage, which works with offset shell elements (option = 11), shown in Figure 4, was chosen. The parameters needed to describe tiebreak contact are presented in Table 2.    The internal force acting per nanosheet was calculated from the experimentally determined [4] tensile strength value as follows, where -experimentally determined [4] ultimate stress; ℎ-nanosheet thickness; -nanosheet width.
The strength of the assembled free-standing Ti3C2Tz-based film depends on the nanosheet interface shear strength, The minimum interlayer strength value achieved by an average overlapping length of 20% satisfies the experimental results [4]. For an FE model of three layers of MXene nanosheets = 2200 nm, the maximum resultant load acting on MXenes nanosheet cross-sections is the following: where -the number of interfaces between three nanosheets. This value was used as the criteria to validate the FE model for interface shear strength simulation. The graphical optimization tool LS-OPT was used for the identification of material constants. The material parameters were obtained using the curve fitting approach with the    The internal force acting per nanosheet was calculated from the experimentally determined [4] tensile strength value as follows, where -experimentally determined [4] ultimate stress; ℎ-nanosheet thickness; -nanosheet width. The minimum interlayer strength value achieved by an average overlapping length of 20% satisfies the experimental results [4]. For an FE model of three layers of MXene nanosheets = 2200 nm, the maximum resultant load acting on MXenes nanosheet cross-sections is the following: where -the number of interfaces between three nanosheets. This value was used as the criteria to validate the FE model for interface shear strength simulation. The graphical optimization tool LS-OPT was used for the identification of material constants. The material parameters were obtained using the curve fitting approach with the  The internal force acting per nanosheet was calculated from the experimentally determined [4] tensile strength value as follows, where σ u -experimentally determined [4] ultimate stress; h-nanosheet thickness; w-nanosheet width.
The strength of the assembled free-standing Ti 3 C 2 T z -based film depends on the nanosheet interface shear strength, where OVL -overlapping coefficient; L-nanosheet length; τ u = 22 (MPa)×20 (nm) 0.2×1000 (nm) = 2.2 MPa. The minimum interlayer strength value achieved by an average overlapping length of 20% satisfies the experimental results [4]. For an FE model of three layers of MXene nanosheets L ovl = 2200 nm, the maximum resultant load acting on MXenes nanosheet cross-sections is the following: where n l -the number of interfaces between three nanosheets. This value was used as the criteria to validate the FE model for interface shear strength simulation. The graphical optimization tool LS-OPT was used for the identification of material constants.
The material parameters were obtained using the curve fitting approach with the parameterized simulation of the physical tests and calibration to the test results. The objective was to minimize the mean squared error between the test results [4] and the FE simulation results.
The loading was described by a curve of linearly increased displacement vs. time. As was mentioned before, this gives a constant velocity, which initiates vibration of the structure and makes its behavior uncertain. Therefore, in the second step, a smooth increase of displacement was chosen for loading, from the assumption that an average speed of tensile loading is equal to~1 m/s, and speed is increasing linearly from 0 up to the 2v avg . The loading curve, displacement versus time, was obtained by integrating the linearly increasing velocity function: u(t) = v(t)dt (Figure 5a). The response F(t) was used for the curve fitting procedure and, for FE model validation, it was transformed into the F(u) curve (Figure 5b). Using Equation (2), the experimental tensile curve [4] was recalculated into a force versus displacement curve and used for model validation.
Materials 2020, 13, x FOR PEER REVIEW 6 of 10 parameterized simulation of the physical tests and calibration to the test results. The objective was to minimize the mean squared error between the test results [4] and the FE simulation results. The loading was described by a curve of linearly increased displacement vs. time. As was mentioned before, this gives a constant velocity, which initiates vibration of the structure and makes its behavior uncertain. Therefore, in the second step, a smooth increase of displacement was chosen for loading, from the assumption that an average speed of tensile loading is equal to ∼1 m/s, and speed is increasing linearly from 0 up to the 2 . The loading curve, displacement versus time, was obtained by integrating the linearly increasing velocity function: ( ) = ( ) (Figure 5a). The response ( ) was used for the curve fitting procedure and, for FE model validation, it was transformed into the ( ) curve (Figure 5b). Using Equation (2), the experimental tensile curve [4] was recalculated into a force versus displacement curve and used for model validation.

FE Simulation of Nanoindentation
The simulated deflection field of the indented MXene monolayer flake is presented in Figure 6. The highest deflection of 34 nm was reached at the center point of the MXene flake. The dependence of velocity of the nanoindenter and the center point of the MXene monolayer flake on time shows a linear behavior (Figure 7a). Analyzing the deflection of the nanoindenter and the center point of the MXene monolayer flake as a function of time, the dependence is observed to be nonlinear (Figure 7b).

FE Simulation of Nanoindentation
The simulated deflection field of the indented MXene monolayer flake is presented in Figure 6. The highest deflection of 34 nm was reached at the center point of the MXene flake. The dependence of velocity of the nanoindenter and the center point of the MXene monolayer flake on time shows a linear behavior (Figure 7a). Analyzing the deflection of the nanoindenter and the center point of the MXene monolayer flake as a function of time, the dependence is observed to be nonlinear (Figure 7b). Materials 2020, 13, x FOR PEER REVIEW 6 of 10 parameterized simulation of the physical tests and calibration to the test results. The objective was to minimize the mean squared error between the test results [4] and the FE simulation results. The loading was described by a curve of linearly increased displacement vs. time. As was mentioned before, this gives a constant velocity, which initiates vibration of the structure and makes its behavior uncertain. Therefore, in the second step, a smooth increase of displacement was chosen for loading, from the assumption that an average speed of tensile loading is equal to ∼1 m/s, and speed is increasing linearly from 0 up to the 2 . The loading curve, displacement versus time, was obtained by integrating the linearly increasing velocity function: ( ) = ( ) (Figure 5a). The response ( ) was used for the curve fitting procedure and, for FE model validation, it was transformed into the ( ) curve (Figure 5b). Using Equation (2), the experimental tensile curve [4] was recalculated into a force versus displacement curve and used for model validation.

FE Simulation of Nanoindentation
The simulated deflection field of the indented MXene monolayer flake is presented in Figure 6. The highest deflection of 34 nm was reached at the center point of the MXene flake. The dependence of velocity of the nanoindenter and the center point of the MXene monolayer flake on time shows a linear behavior (Figure 7a). Analyzing the deflection of the nanoindenter and the center point of the MXene monolayer flake as a function of time, the dependence is observed to be nonlinear (Figure 7b).  The influence of MXene monolayer flake thickness ( Figure 8a) and Young's modulus (Figure 8b) on the force vs. deflection curve was analyzed. The results of thickness influence showed that if the Young's modulus of the MXene monolayer flake is set as E = 333 GPa, the best fit to the experimental result is obtained when the thickness is equal to 1.1 nm, which is higher than results presented in the literature [1,2,11]. On the other hand, if we use a thickness of 0.98 nm, the best fit appears at a Young's modulus of 380 MPa. To achieve the same results as the experimental ones, for finite element modeling, the thickness of the MXene monolayer flake should be increased from 0.98 to 1.1 nm, or the Young's modulus should be increased from 333 to 380 GPa.

FE Simulation of Pure MXene Film
To simulate the behavior of pure MXene film under tension conditions, three constants were selected for calibration and sensitivity analysis: interlaminar shear strength, stiffness, and the energy release rate used for damage calculation. The results obtained by the curve fitting approach are shown in Figure 9. The best fitted FE curve was obtained using an interlaminar shear strength interl = 2.2 MPa, in-plane stiffness interl = 0.26  The influence of MXene monolayer flake thickness ( Figure 8a) and Young's modulus (Figure 8b) on the force vs. deflection curve was analyzed. The results of thickness influence showed that if the Young's modulus of the MXene monolayer flake is set as E = 333 GPa, the best fit to the experimental result is obtained when the thickness is equal to 1.1 nm, which is higher than results presented in the literature [1,2,11]. On the other hand, if we use a thickness of 0.98 nm, the best fit appears at a Young's modulus of 380 MPa. To achieve the same results as the experimental ones, for finite element modeling, the thickness of the MXene monolayer flake should be increased from 0.98 to 1.1 nm, or the Young's modulus should be increased from 333 to 380 GPa. The influence of MXene monolayer flake thickness ( Figure 8a) and Young's modulus (Figure 8b) on the force vs. deflection curve was analyzed. The results of thickness influence showed that if the Young's modulus of the MXene monolayer flake is set as E = 333 GPa, the best fit to the experimental result is obtained when the thickness is equal to 1.1 nm, which is higher than results presented in the literature [1,2,11]. On the other hand, if we use a thickness of 0.98 nm, the best fit appears at a Young's modulus of 380 MPa. To achieve the same results as the experimental ones, for finite element modeling, the thickness of the MXene monolayer flake should be increased from 0.98 to 1.1 nm, or the Young's modulus should be increased from 333 to 380 GPa.

FE Simulation of Pure MXene Film
To simulate the behavior of pure MXene film under tension conditions, three constants were selected for calibration and sensitivity analysis: interlaminar shear strength, stiffness, and the energy release rate used for damage calculation. The results obtained by the curve fitting approach are shown in Figure 9. The best fitted FE curve was obtained using an interlaminar shear strength interl = 2.

FE Simulation of Pure MXene Film
To simulate the behavior of pure MXene film under tension conditions, three constants were selected for calibration and sensitivity analysis: interlaminar shear strength, stiffness, and the energy release rate used for damage calculation. The results obtained by the curve fitting approach are shown in Figure 9. The best fitted FE curve was obtained using an interlaminar shear strength τ interl = 2.2 MPa, in-plane stiffness E interl = 0.26 GPa µm , and a shear energy release rate G IIc = 3.8 × 10 −2 J m 2 . The results show that the behavior of tensioned film of free-standing MXene nanosheets is most sensitive to their interlaminar shear strength and in-plane stiffness, while the shear energy release rate mainly influences the interlaminar damage properties of nanosheets. The material parameters identified by FE modeling were based on 20% of the random average overlapping of nanosheets. In the validated FE model, we changed the overlapping of nanosheets by up to 50% and evaluated the behavior of the ideal overlapping case (Figure 10a). By the FE simulations obtained, the force versus displacement relations were transformed into stress-strain curves and compared with the experimental tensile curve [4] (Figure 10b). The simulation results, based on previously identified material parameters and maximum overlapping (50%) compared with the experimental tensile curve, show the increase of the free-standing Ti3C2Tz film strength, stiffness, and failure strain when the overlapping is increased.
It is clearly seen (Figure 10b, red curve) that the random overlapping of nanosheets decreases the strength and stiffness of films assembled from free-standing MXene nanosheets. The FE models with curve fitting approaches can be used for material constants and adhesion energy identification, and only a proper statistical interpretation of the geometrical parameters of the nanosheets is needed.

Conclusions
The explicit FE model was developed in LS-Dyna to simulate the nanoindentation process of the 2D Ti3C2Tz MXene monolayer flake and validate the material model. Sensitivity analysis of Young's modulus and the flake thickness was performed. The obtained force versus deflection curves showed The material parameters identified by FE modeling were based on 20% of the random average overlapping of nanosheets. In the validated FE model, we changed the overlapping of nanosheets by up to 50% and evaluated the behavior of the ideal overlapping case (Figure 10a). By the FE simulations obtained, the force versus displacement relations were transformed into stress-strain curves and compared with the experimental tensile curve [4] (Figure 10b). The simulation results, based on previously identified material parameters and maximum overlapping (50%) compared with the experimental tensile curve, show the increase of the free-standing Ti 3 C 2 T z film strength, stiffness, and failure strain when the overlapping is increased. The material parameters identified by FE modeling were based on 20% of the random average overlapping of nanosheets. In the validated FE model, we changed the overlapping of nanosheets by up to 50% and evaluated the behavior of the ideal overlapping case (Figure 10a). By the FE simulations obtained, the force versus displacement relations were transformed into stress-strain curves and compared with the experimental tensile curve [4] (Figure 10b). The simulation results, based on previously identified material parameters and maximum overlapping (50%) compared with the experimental tensile curve, show the increase of the free-standing Ti3C2Tz film strength, stiffness, and failure strain when the overlapping is increased.
It is clearly seen (Figure 10b, red curve) that the random overlapping of nanosheets decreases the strength and stiffness of films assembled from free-standing MXene nanosheets. The FE models with curve fitting approaches can be used for material constants and adhesion energy identification, and only a proper statistical interpretation of the geometrical parameters of the nanosheets is needed.

Conclusions
The explicit FE model was developed in LS-Dyna to simulate the nanoindentation process of the 2D Ti3C2Tz MXene monolayer flake and validate the material model. Sensitivity analysis of Young's modulus and the flake thickness was performed. The obtained force versus deflection curves showed It is clearly seen (Figure 10b, red curve) that the random overlapping of nanosheets decreases the strength and stiffness of films assembled from free-standing MXene nanosheets. The FE models with curve fitting approaches can be used for material constants and adhesion energy identification, and only a proper statistical interpretation of the geometrical parameters of the nanosheets is needed.

Conclusions
The explicit FE model was developed in LS-Dyna to simulate the nanoindentation process of the 2D Ti 3 C 2 T z MXene monolayer flake and validate the material model. Sensitivity analysis of Young's modulus and the flake thickness was performed. The obtained force versus deflection curves showed lower results than the experimental ones. To achieve similar results to the experimental ones, the thickness of the MXene flake should be increased from 0.98 up to 1.1 nm, or the Young's modulus should be increased from 333 up to 380 GPa.
The reverse engineering method with the curve fitting approach was applied to evaluate the adhesive strength of the assembled free-standing Ti 3 C 2 T z -based films. The interlaminar shear strength, stiffness, and energy release rate were selected as the main variables; these variables affect the tensile strength of the assembled free-standing Ti 3 C 2 T z -based film. The experimental tensile test data has been used for the curve fitting approach. A simulation was performed using three and nineteen layers of single-layered nanosheets with 20% overlap. The best fitted FE curve was obtained using an interlaminar shear strength τ interl = 2.2 MPa, in-plane stiffness E interl = 0.26 GPa µm , and shear energy release rate G IIc = 3.8 × 10 −2 J m 2 . The results of the sensitivity analysis showed that the largest influences on the behavior of tensioned film of free-standing MXene nanosheets are interlaminar shear strength and in-plane stiffness, while the shear energy release rate mainly affects the interlaminar damage properties of nanosheets. The simulation results based on identified material parameters showed an increase in the free-standing Ti 3 C 2 T z film strength, stiffness, and failure strain when the overlapping was increased by 50%.
The developed FE models with curve fitting approaches are of general purpose and can be used to determine material strength and stiffness properties and adhesion energy identification in different cases; however, a proper statistical interpretation of the single-layered nanosheet geometrical parameters is needed.
The obtained values of interlaminar shear strength, stiffness, and energy release rate in MXenes are very important parameters for further investigation of MXene-polymer composites by finite element modeling. Areas for further research include analysis of the influence of multilayered or agglomerated nanosheets on the composite properties as well as development of MXene-polymer coatings with high-density MXene-MXene interactions.