A Data Driven Modelling Approach for the Strain Rate Dependent 3D Shear Deformation and Failure of Thermoplastic Fibre Reinforced Composites: Experimental Characterisation and Deriving Modelling Parameters

: The 3D shear deformation and failure behaviour of a glass ﬁbre reinforced polypropylene in a shear strain rate range of ˙ γ = 2.2 × 10 − 4 to 3.4 1 / s is investigated. An I OSIPESCU testing setup on a servo-hydraulic high speed testing unit is used to experimentally characterise the in-plane and out-of-plane behaviour utilising three specimen conﬁgurations (12-, 13- and 31-direction). The experimental procedure as well as the testing results are presented and discussed. The measured shear stress–shear strain relations indicate a highly nonlinear behaviour and a distinct rate dependency. Two methods are investigated to derive according material characteristics: a classical engineering approach based on moduli and strengths and a data driven approach based on the curve progression. In all cases a J OHNSON –C OOK based formulation is used to describe rate dependency. The analysis methodologies as well as the derived model parameters are described and discussed in detail. It is shown that a phenomenologically enhanced regression can be used to obtain material characteristics for a generalising constitutive model based on the data driven approach.


Introduction
Fibre reinforced plastics (FRPs) have become a widely used material for a broad range of structural applications. Especially, the high capacity for energy dissipation is exploited in components exposed to crash and impact loading scenarios in the automotive [1,2] and aviation industries [3,4].
The key aspect for a reliable design process is the characterisation and accurate description of the material behaviour considering arbitrary three-dimensional loading cases as well as loading rates. Characterisation of material properties is particularly challenging for the through thickness (TT) shear behaviour of composites [5]. They typically exhibit pronounced nonlinear behaviour, which is strongly influenced and changed by the presence of other stress components [6,7]. The challenges of determining such behaviour arise mainly from two domains: (a) specimen design and experimental evaluation and (b) deduction of TT material characteristics and material modelling. Regarding test design and evaluation, many different approaches have been investigated. Jalai and Taheri developed the varying span method, where three point bending experiments with different span to thickness ratios are used in combination with extensive analytical models to determine the initial TT shear modulus [8]. Yoneyama et al. used a curved beam specimen, also subjected to three point bending, with permanent strain measurement throughout the experiment, enabled by digital image correlation (DIC) [9]. They determined the interlaminar shear modulus, as well as other elastic properties, at predefined load levels by utilising supporting finite element analyses (FEA). Several authors found a V-notched specimen in combination with DIC for strain measurement to be suitable for determining quasi-static shear characteristics [10][11][12].
Regarding the experimental determination of the shear properties of fibre reinforced plastic (FRP) at elevated velocities, only a few investigations have been published. Hsiao et al. used a combination of servo-hydraulic test machine, drop tower and split HOP-KINSON pressure bar (SHPB) to investigate a wide range of strain rates [13]. Zhao et al. employed a combination of a universal testing machine and an electromagnetic SHPB [14]. A modified IOSIPESCU-test setup on a servo-hydraulic high speed testing machine [15] was successfully used by Hufenbach et al. to determine rate dependent shear properties in composite materials [16]. The experimental investigations in the present study follow the latter approach.
Especially for FEA with respect to crash and impact applications, the in-depth knowledge of the characteristic 3D-stress-strain and -failure behaviour is of high importance. With respect to delamination failure in laminated textile reinforced composites, interlaminar shear properties are of special relevance for the determination of failure onset. Typical approaches to model the material behaviour are based on phenomenologically driven models, in combination with engineering constants. To capture the nonlinearity of the stress-strain relationship under TT shear, different approaches with varying experimental effort and potential for generalisation have been applied in literature. The approach requiring the least additional effort is to define parameter pairs of loading state and corresponding modulus [9]. It yields a good approximation close to the sampling points but is not well suited for extensive extrapolation. In regions of high nonlinearity, a high density of sampling points is required to capture the curve's progression. Alternatively, attempts to model the underlying effects leading to the deviation from linear behaviour have been made. Hassan et al. superposed inelastic behaviour, continuum damage and interface failure to analyse the energy dissipation during quasi-static failure of composite joints [17].
However, the impact of the strain rate on the TT shear behaviour of textile reinforced composites and the respective characteristics has not yet been investigated in detail. In order to add strain rate dependency to initially static models, often a JOHNSON-COOK [18] based approach is used. This enables the model's application to a wider range of problems. However, it doubles the number of model parameters necessary, significantly increasing the effort for parameter identification. This effect is further amplified by the additional interaction parameters often used to unify different loading scenarios [19][20][21]. Additionally, the simultaneous occurrence of different sources of nonlinearity during the loading of composites poses an enormous challenge to realise experimentally and to extract said parameters from the resulting curves. In the field of FRP, no standardised models or methods for parameter identification exist; it is rather up to the engineer to find a suitable model, familiarise oneself enough to understand the physical meaning of the material parameters, what kind of experimental procedure is necessary to trigger underlying effects and how to extract the necessary parameters from the experimental results.
Therefore, the aim of this paper is to develop the simplest model possible-with as few parameters as possible-still capable of accurately predicting the macroscopic material response. This type of model is chosen to provide a proof of concept for the proposed modelling approach as well as a stepping stone for the usage of more sophisticated data driven (DD) methods in the field of material modelling-in particular, considering that currently established strategies which focus on deriving 'stand alone' engineering parameter, based on e.g., ASTM or DIN-ISO standards, tend to disregard the majority of the acquired data. This is particularly relevant for highly nonlinear constitutive material behaviour, where conventional characterisation methodologies significantly reduce relevant information on the material behaviour initially observed experimentally, which may be important for subsequent modelling purposes. In addition, more comprehensive data and information are demanded by the rapidly increasing application of (semi-)automatised methods for data processing and analysis. The presented work aims to contribute to a paradigm shift: from material characteristics, which are evaluated at predefined points, towards a characteristic curve based one. To achieve this, here, a purely regression based model is derived from the development of the tangent modulus. The resulting DD model is compared with a classical modulus and strength based (MaS) approach, keeping track of the change in secant stiffness at predefined strain states.

Material Configuration and Specimen Preparation
For all presented experiments, a multi-layered weft knitted fabric (MKF) reinforcement was used, which has previously been developed and manufactured at the TU Dresden [22]. A characteristic feature of the fibre architecture is the non-crimped fibres of the warp and weft threads without any undulations. Knitting loop threads secure the fibre interlock and prevent the delamination within each individual layer (detailed schematic and computed tomographic illustrations are given in [23]).
The MKF consists of commingled hybrid E-glass fibre polypropylene (GF/PP) rovings, with a fibre-fineness of 1400 tex in 1-(warp) and 2-direction (weft) and a loop yarn of 139 tex, representing a layer-wise 3D-reinforcement. The result is an equal mass share of 42 % (warp and weft) and 16.5 % (loop) and a total fibre volume fraction of 55 %.
The base plates have been manufactured using the hot pressing technology. A more detailed description of the manufacturing process is given in [15]. The thickness of a single consolidated textile layer equals 0.5 mm. Specimens have been cut from the plates by water jet cutting.

Experimental Setup and Testing Program
Three different specimen configurations have been used for the evaluation of the strain rate dependent shear behaviour of the GF/PP. Figure 1a illustrates the alignment of each shear specimen configuration: 12-(in-plane), 13-and 31-configuration (out-of-plane). The notation of the specimen configuration conforms to the orientation of the material axes related to the specimen dimensions whereupon the first index corresponds to the length dimension and the second to the width dimension. The specimen has a length of 78 mm, a width of 20 mm and is 45°V-notched in the middle with a notch tip distance of 13 mm. The thickness is 4 mm for 12-and 13-configuration and 10 mm for the 31-configuration. In case of the 31-configuration, a higher thickness is necessary since the material strength is very low, due to the layers being sheared off of one another. The determination of the shear properties has been performed with the IOSIPESCU shear testing device shown in Figure 2a. This test arrangement is in accordance with ASTM D 5379 [24], where the specimen, which is notched on both sides, is clamped and, when compressed, a zone of torque-free shear load between the notches is created [15,16]. The load is applied with a servo-hydraulic test system Instron VHS 160/20 with a load cell of 160 kN, enabling tests at deformation speeds of up to 20 m /s with a load cell accuracy of ±0.5 %. Deformation measurements and subsequent strain analysis were performed using high speed DIC, based on the stochastic grey scale pattern on the specimens' surface. For that, two high speed cameras with a maximum frame rate of 200,000 images per second in combination with the 3D DIC system Aramis by GOM mbH have been used.
The shear strain γ is calculated by averaging the values within a region of interest (ROI) and determined by the principal strains ε I and ε I I : Anticipatory to the experimental results and as illustrated in Figure 2b, the shear distribution is of inhomogeneous nature in the analysed area. Therefore, a constant ROI of 15 mm × 3 mm is chosen for all experiments, in which the shear deformation is analysed and averaged, resulting in one representative value. The chosen ROI covers the area between the V-notches of the specimen with a shear-dominated deformation ratio.
A current associated shear strain rateγ is determined during the the optical measurements byγ within the ROI, where the measurement time step ∆t corresponds to the camera acquisition rate. The associated shear stress τ is calculated with respect to the loading force at the smallest specimen cross sectional area between the notch tips. Each of the investigated configurations (12, 13 and 31) was tested at four nominal loading velocities: 0.00001, 0.0001, 0.01 and 0.1 m /s. For each of the three lower loading velocities, four repetitions were performed, three for the highest, leading to a total of 45 experiments. The nominal shear strain rateγ of each experiment is calculated between γ = 0.15 % and γ = 0.25 % and is seen as the reference value for the entire experiment, although the value varies slightly throughout the experiment. Subsequently, the shear strain rates determined for the respective experiments are averaged over all samples for each material configuration and loading velocity. An overview over the conducted experiments can be found in Table 1. The slight differences in strain rate at equal loading velocity and geometry are due to differences in the strain localisation behaviour for the different material configurations.
Tests per conf. 15 15 15 * In one experiment, the strain measurement was not triggered. It is excluded from the further analyses.

Stress-Strain Behaviour and Experimental Classification
The experimental stress-strain results for the investigated 12-, 13-and 31-configuration are illustrated in Figure 3.   The curves of the 12-and 13-configuration all exhibit an almost linear behaviour at very low strain. With increasing strain, the curves become strongly nonlinear until a relinearisation occurs. In case of the 12-configuration, a slight slope remains, whereas, in the 13-configuration, a plateau of constant stress is reached. The results of both configurations are clipped at a strain of 20 %. At higher strains, the curves begin to steepen quickly. This behaviour can be attributed to the fibres aligning with the load application. Thus, its cause can be found in the specimen's overall structure, not the material behaviour itself. Therefore, the steepening aspect is not considered in subsequent analyses.
A low pass filter was applied to the results from the 31-configuration, since high frequency oscillations dominated the raw measurements. The identifiable behaviour begins similar to the other two configurations, though the linear regime is limited to strains smaller than 0.5 %. Again, it is followed by strongly nonlinear behaviour. However, catastrophic failure occurs before re-linearisation can take place.
Due to the highly nonlinear nature of the stress-strain curves, the identification of single representative stiffness or failure parameters analogous to materials of brittle characteristic is not possible. Therefore, two approaches to characterise the behaviour are investigated in Section 4.

Failure and Fracture Behaviour
An exemplary post-experimental failure pattern for each investigated material configuration is shown in Figure 4. An influence of the strain rate on the observable failure patterns could not be identified. The displayed images are therefore considered to be representative across the investigated strain rate domains. For both the 12-and the 13-configuration, nonlinearities dominate the stress-deformation behaviour. In both the 13-and the 31-configuration, the predominant failure mode is of an interlaminar nature. This results in delaminations, which in turn lead to catastrophic failure and specimen separation in the 31-configuration. A more in-depth assessment of the observed phenomena and their origins is given in [15], where a 2/2-twill weave reinforcement is made of commingled hybrid GF/PP yarn and experimentally investigated in a similar manner, without addressing strain rate effects. Their findings are consistent with the ones obtained in this study for MKF reinforcement.

Material Modelling and Property Identification
Two approaches to model the experimentally identified material behaviour are pursued: 1.
modulus and strength based (MaS): An analysis based on ASTM D5379 [24] and DIN EN ISO 14129 [25] with the determination of engineering constants at fixed strain values. 2.
data driven (DD): A closed formulation describing the entire experimental curves by coherent formulae with the determining parameters being the material's characteristics.
In both cases, the causes of the nonlinearities are not further investigated for the modelling, since all experiments were carried out with continuously increasing strain, making a differentiation of intrinsic effects impossible.

Determination of Engineering Constants
Four slopes, referred to as shear moduli, were determined for each experiment (see Table 1): three secant moduli G sec (between γ = 0.15 % and γ = 0.55, 1 and 2 %) and a tangent modulus G tan 5 % at γ = 5 % determined between γ = 4.8 and 5.2 %: The determined parameters are arithmetically averaged over the conducted experiments for each strain rate and material configuration. They are summarised in Table 2, and respective maximum stress values of the experiments are given additionally.

Strain Rate Dependency and Model Parameters
The effects of the shear strain rate on the material parameters G sec 0.55 % and stress levels τ max for the respective shear directions 12, 13 and 31 within the considered strain rate range are presented in Figure 5. The remaining values at 1 and 2 % exhibit similar behaviour.   In order to incorporate the increase of those parameters with the strain rate, a model originally proposed by JOHNSON and COOK [18], which is widely used in literature to model the strain rate dependency of model parameters, e.g., [26][27][28], is chosen: The strain rate dependency was found to be accurately described by this model. In (4), A G and A τ denote model parameters, which control the linear slope (in the log plot) of the model's predictions over the considered shear strain rates. These parameters are determined by a best fit approach. G re f and τ re f are the values at reference shear strain ratė γ re f . It is worth emphasising that this type of natural logarithm-formulation may result in unreasonable negative values when the model is evaluated far below the reference strain rate. Special attention has to be paid to this in implementations for FEA, since large local jumps in strain rate may occur, leading to negative moduli.
With the identified strain rate material constants A G 12,31/13 and A τ 12,31/13 for the corresponding 12, 31 and 13-configurations, the tendency of the rising values of shear modulus and shear stress respectively can be accurately estimated. An overview of the model constants for selected material parameters is given in Table 3. In this regard, it is emphasised that the shear moduli G 13 and G 31 are considered to be identical due to the typical assumption of a symmetric stiffness tensor for these textile reinforced composites. Reference values were taken at the lowest measured shear strain rates to avoid the aforementioned problems at strain rates lower than the reference one. The corresponding model predictions are presented for the initial shear modulus G sec 0.55 % (γ) and the maximum shear stress τ max (γ) in Figure 5. The resulting trajectories at different strains are presented in Figure 6 for the 12-configuration and strain rates of 2.2 × 10 −4 and 3 1 /s.

Determination of Material Characteristics
Within this subsection, the applied strategy to model the material's behaviour by continuous formulae is exemplary presented using the 12-configuration. For the other configurations, an identical procedure is used unless stated otherwise. In order to obtain a continuous description of the material's behaviour, the slopes along the curves-hereafter referred to as tangent modulus G tan (γ)-are investigated. Therefore, each experimental curve is considered as n pts subsequent stress-strain points and preprocessed by applying a low-pass filter. The tangent modulus at a given strain point γ k , with k < n pts − 1, is subsequently calculated by The results are presented in Figure 7. The curves show a finite value at zero strain and a rapid decay with an asymptotic approach to zero. Given these curve characteristics, an exponential law is chosen as model: with the initial tangent modulus G tan 0 and the decay parameter a. Integration of (6) with the initial condition τ(γ = 0) = 0 yields a closed form for the stress-strain relationship: From this, it becomes clear that the chosen model contains the asymptotic stress τ ∞ , which can be seen for the material configurations 12 and 13 in Figure 3 as Parameter identification for G tan 0 and a is carried out individually for each test by a standard least square fit. In all cases, both parameters exhibit a distinct dependency on the strain rate. The experimental curves in combination with the corresponding model predictions based on the previously identified parameters G tan 0 and a are presented in Figure 8 ("Model, standard").  The respective dependency of the identified model parameters in (7) on the strain rate can be found in Figure 9 ("Standard fit"). In the case of the decay parameter a (Figure 9b), it becomes apparent that there is no simple dependency of the material parameters on the strain rate. In particular, the simultaneous increase and decrease fromγ = 2.2 × 10 −4 1 /s to 3.6 × 10 −3 1 /s pose a major challenge to deterministic models. However, the usage of such an approach is highly desirable to obtain a more generalised and predictive model. Therefore, a second set of parameters is identified.
For this second methodology, additional bounds on the fitting parameters are applied: The fits at the lowest tested strain rate are carried out without any restrictions. For all subsequent fits, the lower bounds for the parameters are set to the means of the parameter at the previous strain rate. This ensures that the well-known tendency of stiffness parameters increasing with the strain rate [16,18,19,29,30] is preserved within the resulting parameter set. The model predictions for the stress-strain curves based on the parameters determined by this constrained method are illustrated in Figure 8 ("Model, constrained") and the respective model parameters in Figure 9 ("Constrained fit").

Modelling of the Strain Rate Dependency
Given the excellent results of the JOHNSON-COOK-model for strain rate dependency in Section 4.1.2, the same approach is taken to model the strain rate dependency of the material parameters G tan 0 and a: with the values G tan 0,re f and a re f at the reference strain rateγ re f and the material constants A G and A a . The values A G and A a are obtained by a standard least square fit. Given the advantageous generalisation capability of the constrained methodology, the final model with the parametersγ re f , G tan 0,re f and a re f as well as the parameters A G and A a ( Table 4) are obtained by fitting the constrained model's parameters. In Figure 10, the resulting curves are plotted alongside the conducted tests. It can be seen that the model is in excellent agreement with the experimental results for the 12-and 13-configuration. In case of the 31-configurations, sound results are obtained as well, even though failure occurred already at shear strains of approximately 2%. The determined material parameters for the respective material configurations are presented in Table 4.

Determining the Range of Validity of the Material Modelling Approaches
A criterion to determine the range of validity for the MaS model-exemplary shown for the initial secant shear modulus G sec 0.55 % -is presented. It is intended to provide an indicator under which circumstances the additional effort required for the DD model is expected to yield a significant improvement in the quality of prediction. As a basis, the degree of nonlinearity and the subsequent deviation from a straight line is used for the assessment.
The transition from linear to nonlinear behaviour can be described by a defined percentage based deviation of the experimental curve from a given straight line with the initial modulus as slope. However, significant oscillations in the low strain area in some of the experiments at high strain rates are observed. Due to the low absolute values, small deviations with respect to the standard testing accuracies already lead to high percentage changes, resulting in an erroneous immediate prediction of the onset of nonlinearity. Therefore, the criterion has to be robust against such small absolute deviations. Thus, a criterion based on the coefficient of determination r 2 is employed, where a maximum value of 1 indicates perfect conformity.
For each experimental point k, the r 2 value is determined for the straight line with the slope G tan 0 approximating the actual experimental curve. In an extension of the definition for an entire data set, e.g., given in [31], this leads to as the coefficient of determination at every point along the curve. Significant nonlinearity is thereby defined as the strain at which the r 2 value reaches the last local maximum since the quality of approximation steadily decreases with proceeding strain. An example for this criterion is presented in Figure 11 up to a strain of 5 %.  This method yields good results for the material configurations 12 and 13 with one exception in the 13 case atγ = 0.34 1 /s. Therefore, the respective specimen is not taken into account for further investigations. Furthermore, it becomes clear that no exact boundaries for the range of validity can be specified. However, from a certain point on, the approximation quality of the linear model constantly decreases. Therefore, this point is chosen as a guideline, even though it must be noted that, in individual cases, the linear approximation might be valid for larger or smaller strains. The determined stresses at which initial nonlinearity becomes significant are presented in Figure 12. A distinct dependency on the shear rate may be observed. Similar to the onset of damage, the stress-strain curves show a distinct point at which the behaviour changes back to an almost linear one with a very small slope. This behaviour is attributed to the saturation of underlying effects. The point of nonlinearity saturation is determined analogously to the onset by linear regression, this time starting at the maximum considered strain of 20 % of the specimen and moving on to lower strains. In this case, the first local minimum of the mean squared error, calculated between the linear approximation and the experimental curve, has been found to be a well-suited criterion. The determined stress values R ∞ are shown in Figure 12. Furthermore, it becomes apparent that the stress at which the re-linearisation occurs increases with the strain rate. The area of the stress-strain curve following R ∞ yields great potential for use cases in which the energy absorption is crucial.

Strain Rate and Configuration Dependent Experimental Characterisation
The investigated GF/PP-MKF shows similar behaviour under the three shear loading conditions 12, 13 and 31 (see Figure 13). In the beginning of the linear regime, the stiffness under 12, 13 and 31 loading is almost identical. This supports the often made assumption of all these kinds of loading being identical for balanced MKF. However, this equality of the curves only holds up to γ ≈ 0.5 %, at which point nonlinearity appears in the 31 case. At γ ≈ 1 %, the 12 and 13 cases show nonlinear behaviour as well. All specimens of the 31-configuration fail in this nonlinear phase. In the other configurations, a re-linearisation occurs. This phenomenon is much more pronounced and progresses much quicker in the 13 case than in the 12. In the 13 case, an almost constant value is reached at γ ≈ 5 %, whereas the nonlinear regime in the 12 case reaches up to γ ≈ 10 %. In both cases, the behaviour stays linear up to strains of approximately 20 %.
The experimental design does not allow for a final assessment of the underlying effects leading to the change in principal material behaviour. This would require unloading of the specimen at intermediate loading states, which currently cannot be implemented with the presented test set-up. However, possible causes are discussed in the context of the developed models and existing literature in more detail in Section 5.3.
Furthermore, the GF/PP's behaviour shows a distinct dependency on the strain rate. In particular, the initial material stiffness increases strongly with the strain rate, and the overall shape of the curve is only influenced slightly. This unequal impact of strain rate on the different aspects of the material behaviour is also captured by the DD-model: For the 12-and 13-configurations, the ratio of the rate dependency parameters A G /A a is around 3, indicating a three times stronger rate dependency of the stiffness.

Presented Modelling Approaches
Two approaches to capture the highly nonlinear material behaviour exhibited by the GF/PP under different shear-loading conditions have been investigated.

Modulus and Strength Based (MaS) Approach
The MaS method is very close to well-known linear elasticity models and can therefore be used almost directly in all FEA-programs. It makes use of established analysis techniques and provides a simple framework. Furthermore, the identification of "material parameters" for constant strain rates is as easy as reading values of experimental data. These aspects greatly facilitate the usage of this approach.
It was shown that the results obtained by the MaS-model can yield good results, as long as the analysis is limited to small strains. Within this range, the model captures the dependency of the material's behaviour on the strain rate with high accuracy. However, it leads to non-smooth curves, with the deviation between experimental and numerical curve depending highly on the number of evaluation points. This necessitates a multitude of "material parameters" for high fidelity analyses, due to the large nonlinear regime within the stress-strain curve. The determination of how many parameters are necessary is aggravated by the fact that the model contains no information about its own applicability, i.e., up to which strain the error is acceptable. In addition, the generalisation to arbitrary strain rates becomes increasingly difficult with every "material parameter", since each requires additional parameters to model this dependency.

Data Driven (DD) Approach
On the other hand, the presented DD approach poses a comparatively high initial hurdle. The underlying model is not necessarily implemented in every FEA-program, requiring a potential user to implement the constitutive model via some subroutine. Furthermore, the material parameters cannot be read directly from the stress-strain curve. It is rather necessary to perform a nonlinear regression analysis of the curve.
In this work, two regression techniques were applied and the resulting parameter sets compared. In Figure 8, it can be seen that the model predictions for both parameter identification methodologies are in good agreement with the experimental data and among themselves. Solely at very high strain rates of 3.0 1 /s and high strains of more than 15 %, the constrained model predictions differ notably from the approximation without additional assumptions. Additionally, the trend of rising material parameters with increased strain rate can be captured more accurately by the constrained fit (cf. Figure 9). This allows for the model to be more generalising and predictive. Furthermore, the two physically based curve parameters G tan 0 and τ ∞ barely differ between both fitting approaches. Therefore, the parameters identified by the constrained methodology are suggested for future works.
The DD approach tremendously outperforms the MaS method in terms of accuracyas long as the number of data points for the MaS method remains somewhat cappedextensibility to arbitrary strain rates and therefore overall generalisability. This is mainly due to the fact that, within the presented framework, the strain rates dependency of every single material parameter is modelled individually with an additional parameter, doubling the necessary parameters if rate effects are to be considered. Furthermore, it provides a closed and differentiable form of the stress-strain relationship. This is highly beneficial for the numeric time integration schemes implemented in commercial FEA-programs.

Interpretability of the Nonlinearities
The respective transitions between the aforementioned material states cannot be distinctly attributed to a material phenomenon, based on the experiments conducted. However, other works on similar materials allow for a preliminary analysis of the underlying effects.
In the case of the first transition, the determined stress value R 0 marks the point from which, going forward, the prediction accuracy of linear models significantly decreases. Therefore, such models should only be used in cases where it is previously known that stresses will not significantly exceed R 0 . However, given the plateau-like nature of the local maximum, minor overstepping of the boundary could only result in slight discrepancies between the models. On the material side, the fibres' support in the form of the matrix begins to degrade and behave nonlinearly. It has been shown that GF/PP exhibits strongly nonlinear viscoelastic behaviour [32,33], which in itself can lead to significant deviations from a linear stress-strain behaviour at strains of less than 1 % [34]. Another common cause for the nonlinearity in the behaviour of FRP can be found in the occurrence of micro-or macroscopic cracks within the matrix, which lead to a decrease in the material's stiffness [35][36][37]. Additionally, PP is known to accumulate plastic strain even at very low levels of stress [38]. This in turn weakens the support of the fibres and therefore limits the material's load bearing capabilities [39]. Recent studies have shown that, in the case of GF/PP, the viscoelastic behaviour is nonlinear from the beginning and that damage and plasticity are triggered at the same time, for the in plane behaviour at shear stresses as little as 1.76 MPa [19].
The re-linearisation occurring at higher strains suggests that the underlying effects begin saturating. In the case of inter-fibre fracture, this behaviour has been thoroughly investigated in recent years [35,37,40,41]. In that case, the material's stiffness has been significantly reduced by the formation of a multitude of cracks; however, the initiation of new cracks is energetically not feasible. Therefore, some loads can still be transmitted even though they lead to drastically increasing defamation, with a return to the linear stress-strain relation.

Conclusions and Outlook
The TT behaviour of GF/PP with MKF reinforcement at various strain rates has been determined experimentally by using a lightweight IOSIPESCU test setup in combination with high speed DIC. It was subsequently modelled using a MaS and a DD method.
The MaS model has been shown to be mostly suitable for preliminary analyses, reaching only small strain values. Otherwise, the amount of model parameters necessary quickly becomes unfeasible.
On the other hand, the DD model represents the entire stress-strain curve by only two parameters. Furthermore, it requires only two additional parameters to generalise to arbitrary strain rates, if the phenomenologically enhanced regression method is used.
Therefore, the resulting model is applicable to high fidelity analyses with monotonically increasing loading, a high number of model and gradient evaluations (i.e., in the case of a large number of elements or time steps) and varying or previously unknown strain rate and state.
Future research of the following two aspects is expected to offer great potential: firstly, the investigation of the constitutive material behaviour describing the orthotropic shear stress-strain relations. In that respect, models incorporating damage and failure as well as multi scale approaches to directly address the influence of the textile architecture should be studied in detail. In particular, the influence of viscoelasticity and plasticity of the composite constituent materials glass and PP [19,42] are identified to be promising to describe the observable strain rate sensitivity; secondly, the pursuit of DD methods in the field of complex material behaviour. Treating experimental results as pure data with known input and output and thus applying methods from the field of supervised machine learning to them is expected to further facilitate model generation. Additionally, such methods could be employed to completely automatise the parameter identification process and thereby significantly lower the initial hurdle for using new material models.  Acknowledgments: The algorithms used for fitting and optimisation are implemented in the open source Python library scipy [43]. Most of the figures presented here were created using the open source Python library matplotlib [44].

Institutional Review Board Statement: Not applicable.
Informed Consent Statement: Not applicable.

Data Availability Statement:
The data presented in this study are available on request from the corresponding author.

Conflicts of Interest:
The authors declare no conflict of interest.

Abbreviations
The following abbreviations are used in this manuscript: