A high-cycle accumulation model for clay and its application to monopile foundations

A high-cycle accumulation (HCA) model predicting the accumulation of permanent strain or excess pore water pressure in clay under a large number of load cycles is presented. Data from an extensive laboratory testing program on kaolin under undrained cyclic loading has been analysed for that purpose. The influence of strain amplitude, void ratio, stress ratio, overconsolidation ratio and loading frequency on the accumulation rates is considered in the constitutive equations of the HCA model. The proposed model is validated first by the simulation of element tests. Subsequently, its application to offshore wind turbine foundations under long-term lateral cyclic loading is presented by the back-analysis of a centrifuge test on a monopile in soft clay. The results are in good accordance with the measurements in terms of pile displacement and bending moment versus number of applied cycles. It is concluded that the proposed model is feasible to describe the long-term behaviour of clay subjected to high-cyclic loading.

Parameter of the HCA model in function f N C N2 Parameter of the HCA model in function f N C N3 Parameter of the HCA model in function Parameter of the AVISA model a Parameter of the AVISA model M c Parameter of the AVISA model e i0 Parameter of the AVISA model f b0 Parameter of the AVISA model R Parameter of the AVISA model m R Parameter of the AVISA model d Parameter of the AVISA model b 0 Parameter of the AVISA model v 0 Parameter of the AVISA model v max Parameter of the AVISA model C a Parameter of the AVISA model H av Average horizontal load H ampl Amplitude of horizontal load Dp w Excess pore water pressure p 0 Initial mean effective stress De pl Increment of plastic strain D/ Increment of plastic potential # Scalar factor used for the calculation of the plastic strain rate p ref Reference mean effective stress for the AVISA model q s Grain density q w Water density

Introduction
Similar to sand, clay can show an accumulation of permanent strain due to cyclic loading with many repetitions. Such a high-cyclic loading is characterized by comparatively small strain amplitudes (larger, however, than for sand) but a considerable number of load cycles. Highcyclic loading is relevant for many structures, the most prominent and actual example being offshore wind turbine (OWT) foundations which are continuously loaded by wind and water waves during their lifetime. A reliable prediction of the permanent rotation of the OWT foundation caused by this loading in the design stage is important, as excessive tilting endangers its serviceability. In the offshore environment, the loading of the foundation will not be perfectly undrained as it is typically assumed in the element tests on clay. The accumulation of (volumetric) strain is thus strongly influenced by the consolidation process occurring simultaneously with the cyclic loading, making the problem at hand even more complex. To take into account these influencing factors and to allow for a reliable estimation of the permanent deformations of OWT foundations subjected to millions of load cycles, the application of numerical methods in combination with special constitutive models can represent a feasible approach. Up to now, only a small number of constitutive models for the prediction of the mechanical response of clay subjected to a large number of load cycles (N'10 3 À 10 4 ) exist (see e.g. the recently proposed model in [19]). In addition to the enormous and unacceptable time required to calculate this large number of load cycles, conventional constitutive models also fail due to a strong overestimation of the residual strains and a considerable numerical error at large numbers of load cycles, the first resulting from the constitutive equations and the latter owed to the large number of calculated increments [55]. In order to account for millions of load cycles, as faced in the offshore environment, special constitutive models such as the so-called high-cycle accumulation (HCA) models can be applied. These models are formulated in general stressstrain rate form, but do not predict the soil response during individual cycles. Only the average trend of the soil state during the cyclic loading is calculated. Opposite to existing methods used to predict the long-term behaviour of piles subjected to cyclic loading, such as p À y curves or strain wedge models, the application of the HCA models is not restricted to a certain boundary value problem. HCA models are suitable to study arbitrary problems with highcyclic loading. For the prediction of the high-cyclic response of sand, the HCA model by Niemunis et al. [55] is well established [39,44,57,63,64,68,79,83]. A corresponding model for clay is missing up to now. The development of a HCA model for unstructured clay is thus the aim of the study presented in this paper. The data of an extensive laboratory testing program on kaolin under undrained cyclic loading are analysed for this purpose. The development of the constitutive equations is presented and the model is validated by the back-analysis of the element tests. To show the applicability and robustness of the proposed model and to extend its validation to a larger scale, a monopile foundation for OWTs under cyclic lateral loading is studied in the framework of a back-analysis of a centrifuge model test.
2 Literature review on clay behaviour under cyclic loading The literature review in this section concentrates on three aspects being of relevance for the development and validation of the HCA model: 1) the behaviour of clay under cyclic loading in laboratory tests, which has to be reflected by the HCA model equations, 2) available conventional constitutive models for clay under cyclic loading which are needed for the low-cycle phases of the HCA simulations as described in Section 3, and 3) model tests with cyclic lateral loading on piles in clay, from which a specific centrifuge test is selected for the validation of the HCA model by means of a back analysis presented in Section 5.

Laboratory tests with cyclic loading on clay
Laboratory tests with cyclic loading on clay are usually performed under undrained conditions. The cyclic loading can lead to a build-up of excess pore water pressure (i.e. effective stress relaxation) and an accumulation of shear strain under those circumstances. Cumulative effects (i.e. build-up of excess pore water pressure and accumulation of shear strain) in clay start if the strain amplitude exceeds a certain threshold value of the strain amplitude e ampl . This threshold strain amplitude increases with increasing plasticity of the test material [32,33,36,76]. The intensity of effective stress relaxation or strain accumulation increases with the stress or strain amplitude [1, 3, 7, 12, 13, 17, 21, 28-30, 32-35, 51, 52, 71, 76, 87, 90]. In undrained cyclic triaxial tests with isotropic average stresses or a relatively small stress anisotropy, clay samples usually fail due to the development of large axial strain amplitudes, exceeding a pre-defined failure criterion (e.g. 10 %). In tests with larger stress ratios jg ¼ q=pj an excessive accumulation of permanent strain usually leads to failure [1,3,35,80]. Some experimental studies reported that a multidimensional cyclic loading (i.e. at least two stress components oscillate out-of-phase) causes larger cumulative effects than one-dimensional cycles (i.e. only one stress component oscillates or several stress components vary in-phase) [18,48,70,84]. At same stress amplitudes, the cumulative rates are lower in overconsolidated clays than in normally consolidated clays. If the stress amplitude is referred to the undrained shear strength s u from monotonic tests, however, the conclusions regarding the tendency are not unique in the literature [4,21,30,35,90]. Most laboratory studies performed in the framework of research projects were conducted on reconstituted samples, often prepared by reconsolidation out of a slurry. Based on a comparison of undisturbed and remoulded samples of several marine clays in undrained cyclic triaxial tests, Hyodo et al. [35] concluded that data of normally consolidated reconstituted samples define an intrinsic cyclic strength. This intrinsic cyclic strength is increased by aging effects, cementation and overconsolidation. There is a consensus in the literature that an increase in the loading frequency leads to a decrease in the accumulation rates [5,14,30,49,51,87,90]. The influence of frequency grows with the plasticity of the clay [1]. Generally, the rates of pore water pressure build-up and shear strain accumulation decrease with increasing plasticity [1,13,27,35,37].

Conventional constitutive models for clay under low-cycle loading
The cyclic behaviour of fine-grained soils depends not only on the material state, such as the effective stress or the density, but also on the strain amplitude and the strain rate. The clay behaves elastically only under very small strain amplitudes of about e ampl \10 À4 . Under medium strain amplitudes 10 À4 \e ampl \10 À2 the stiffness increase at load reversals is accompanied by a reduction in the plastic strain rate [15,80]. With increasing distance to the reversal point, the stiffness decreases again, as observed on a typical stiffness degradation curve. Finally, when the soil is subjected to shearing with large strains jjDejj [ 10 À2 it tends asymptotically to the critical state corresponding to failure under monotonic loading and resulting in cyclic mobility under cyclic loading. These are some relevant aspects of clay behaviour which are expected to be considered by competent constitutive models. A single model capturing all these effects is rare to find if not impossible, and therefore, researchers recommend the usage of an ''appropriate'' model for each particular problem. For example, one could choose an elastoplastic model (e.g. Modified Cam Clay (MCC) model [11]), a hypoplastic model [46,74] or a barodesy model [50] for the simulation of a bearing capacity problem of a shallow foundation on clay. These models are only able to simulate monotonic loading for non-viscous clays, and as in detail discussed in [72], they would not capture the strain rate dependency typical for plastic clays (as is the case with soft marine clays). If the clay shows a high plasticity, a ratedependent model could be useful [10,86]. These models are, however, again restricted to monotonic loading conditions.
For cyclic loading with low to medium number of cycles, the constitutive model should certainly consider small-strain and memory effects. For that purpose, hypoplastic models can be extended by the intergranular strain concept [56] or the intergranular strain anisotropy (ISA-plasticity) [24]. For example, the visco-hypoplastic model of 2003 [53] includes the intergranular strain tensor, while the later version of 2009 [54] considers induced anisotropy, but no small-strain effects. Alternatively, the SANICLAY model [61] can be applied. A substantial drawback of most of these models is the incapability to adequately describe the soil behaviour under different loading amplitudes with the same set of material parameters. Intending to embrace a wide range of strain amplitudes with a single constitutive equation the ISA-plasticity for sand [24,25] was reformulated by Tafili & Triantafyllidis [73] to account for the specific behaviour of clays. This proposed anisotropic visco-ISA (AVISA) model incorporates the rate-and time-dependency of clays as well as inherent anisotropy. In 2020, the ISA-plasticity model was coupled with the hypoplastic model for finegrained soils developed by Mašín [26]. This model also accounts for the inherent anisotropy of fine-grained soils, but neglects rate-and time-dependent effects. In the same year, barodesy has been coupled with the intergranular strain concept in [9], however, without introducing the inherent anisotropy or the time-dependency into the model.

Piles in clay subjected to cyclic loading
Due to cyclic lateral loading, a reduction in stiffness of the fine-grained soil surrounding a pile with increasing number of cycles occurs, which is mainly a result of a build-up of excess pore water pressure. This has for instance been observed in field tests with up to one thousand lateral load cycles on piles (diameters D ¼ 0:273 m and D ¼ 1:22 m) in stiff overconsolidated clay [22]. Similar results were obtained in recent field tests on monopile foundations with a diameter of 2.2 m in soft clay reported in [91] as well as in the pioneer works of Matlock [47] (in soft clay) and Reese et al. [60] (in heavily overconsolidated clay). In [77] a strong degradation of the lateral stiffness of the pile (D ¼ 0:8 m, founded in soft clay) was observed during the cyclic loading as soon as the magnitude of the load was larger than a certain threshold level. The largest degradation occurred during the first five cycles followed by a stagnation of degradation. In case of cyclic loading with even higher load magnitude, the rate of reduction in pile stiffness increased with increasing number of applied cycles without showing the aforementioned stagnation. Similar observations were made for piles with a diameter of D ¼ 0:153 m in overconsolidated clay reported in [89]. Recently, Yang et al. [85] have performed centrifuge tests on monopile foundations (D ¼ 5:9 m in prototype scale) in soft kaolin clay subjected to lateral cyclic loading with hundreds of load cycles. These tests are considered for the back-analysis presented in Sect. 5 applying the HCA model for clay proposed hereinafter. In accordance with the observations documented in [58,59], the cyclic loading lead to an increase in the magnitude of the bending moment and a decrease in the pile stiffness. Small scale model tests on piles (D ¼ 50 mm) in marine clay reported in [41] with several thousand lateral load cycles showed an increasing permanent deflection of the pile for cyclic amplitudes in the range of 30 % to 60 % of the ultimate lateral resistance. Alternate cyclic loading with phases of reconsolidation between the cyclic events in centrifuge tests on monopiles (D ¼ 4 m and D ¼ 6 m in prototype scale) in soft clay revealed a partial recovery of the degraded lateral pile stiffness due to the resting period [40]. However, within each cyclic load package, the cumulative displacement of the pile increased with number of applied load cycles. Recently performed field tests on piles of different size (D ¼ 0:273 m, D ¼ 0:762 m and D ¼ 2 m) installed in an overconsolidated clay employing different rates of lateral monotonic loading showed that the piles exhibited a much higher resistance for higher load rates [16]. In addition, the formation of a gap between soil and pile was reported for all tests, which is believed to contribute to the reduction in the stiffness of the pile system when subjected to cyclic lateral loading.
A reduction in soil reaction force with ongoing cyclic lateral loading is incorporated in a simplified manner in the cyclic p À y curves of the American Petroleum Institute (API) [6] which are based on field tests [47]. Since the API cyclic p À y curves do not allow for a consideration of the number of applied load cycles or the load characteristics (e.g. the ratio of average and cyclic load), many modifications have been proposed over the span of the last decades. Such improved cyclic p À y curves have for instant recently been presented in [88], where the reduction in the soil reaction with ongoing cyclic loading has been incorporated using cyclic strain contour diagrams obtained from direct simple shear tests (see [2]). An extension of the p À y curves considering rate effects has been presented in [8].
Simulations of centrifuge tests on monopiles subjected to lateral cyclic loading in soft clay using a hypoplastic model with intergranular strain extension have been presented in [23,31]. However, due to the nature of conventional constitutive models, only a limited number of load cycles can be taken into account by such an approach. In order to consider millions of load cycles, special constitutive models such as the HCA model are necessary.
3 General concept of a HCA model and the HCA model for sand In the framework of a HCA model the strain resulting from cyclic loading is split in its oscillating part e ampl and its average part e av . The constitutive model works merely in terms of average values and calculates e av only. Thus, the time t can be replaced by the number of cycles N. For cyclic loading with a large number of cycles this is advantageous since it is not necessary to consider individual cycles during the simulation. This reduces not only the computational effort considerably but also avoids an accumulation of numerical errors within each increment, considering that millions of increments are necessary for the simulation of a large number of cycles. Two numerical strategies are combined in a simulation with the HCA model. They are termed the low-cycle (conventional, implicit) and the high-cycle (N-type, explicit) mode of calculation. The low-cycle mode uses a conventional constitutive model. Monotonic loading phases or a few complete cycles are calculated in this mode. The low cycle-mode is necessary to determine the field of the scalar strain amplitude e ampl , which is an important input parameter of the HCA model. The strain amplitude is determined in each integration point of each finite element from the tensorial strain path recorded during an individual load cycle calculated in the low-cycle mode. As depicted in Fig. 1, this can either be the second cycle or a so-called The strain amplitude e ampl is determined from the second loading cycle. It can be updated after a given number of cycles by interrupting the high-cycle phase by additional (update) cycles in the low-cycle mode update cycle. The first cycle should not be used for the calculation of the strain amplitude, because the first quarter of the first cycle represents an initial loading, while all further cycles represent un-and reloading. Based on the recorded tensorial strain path (six independent components of the strain tensor for 3D cases), the strain amplitude is defined as the Euclidean norm of six spans obtained by consecutive projection on (hyper-) planes (for details see the work presenting the HCA model for sand [55]). For one-dimensional strain paths this definition of a multidimensional strain amplitude simplifies to the classical definition e ampl ¼ ðe max À e min Þ=2. This amplitude definition is also adopted for the HCA clay model proposed in this work. During the high-cycle phases the strain amplitude is assumed constant. The field of the strain amplitude can be updated using an update cycle, which is simulated in the low-cycle mode again.
In the high-cycle mode the HCA model is used to predict the development of the average values of strain and effective stress. The basic equation of the HCA model for sand reads with the (objective) stress rate _ r of the effective Cauchy stress r (compression positive), the strain rate _ e (compression positive), the accumulation rate _ e acc , a plastic strain rate _ e pl (necessary only for stress paths touching the yield surface) and the pressure-dependent elastic stiffness E. In the framework of HCA models the dot over a symbol denotes a derivative with respect to the number of cycles N, i.e. _ t ¼ o t =oN. Note that when integrating Eq. (1) numerically in a finite-element (FE) analysis it has to be integrated over DN ¼ Dt=t cycle where t cycle is the duration of a single cycle (period) and Dt is the time increment applied by the FE program. In other words, although the constitutive equations of the HCA model are integrated with respect to N, the simulation is performed in time domain. Therefore, the modelling of transient processes is possible as with regular constitutive models. Depending on the boundary conditions, Eq. (1) predicts a change of stress ( _ r 6 ¼ 0) and / or an accumulation of strain ( _ e 6 ¼ 0). For the HCA model for sand, the rate of accumulated strain in Eq. (1) is calculated from: The factors f t of the intensity of accumulation _ e acc take into account the influence of the strain amplitude e ampl (function f ampl Þ, the cyclic preloading ( _ f N , using the preloading variable g A which weights the number N of applied cycles with the strain amplitude e ampl of these cycles), the void ratio e (f e ), the average mean effective stress p av (f p ) and the normalized average stress ratio Y av (f Y , Y av = 0 at isotropic stresses, Y av = 1 at critical stress ratio). Note that the factor f p describing the influence of changes of the cyclic loading direction used in previous publications on the model has been omitted here, because a recent experimental study proved it to be of minor importance [82]. m is the direction of accumulation (unit tensor), which determines the ratio of the rates of volumetric and deviatoric strain accumulation.

HCA model for clay
Equation (1) used in the HCA model for sand is adopted for the clay model as well. The stiffness E, the intensity of accumulation _ e acc , the plastic strain rate _ e pl and the direction of accumulation m have to be adjusted for the clay model, however. In the following, the development of suitable functions for _ e acc (Sect. 4.1) is presented. The investigation of m is presented in Sect. 4.2, while E and _ e pl are discussed in Sect. 4.3.
In order to study the influencing factors on the intensity of accumulation, an extensive laboratory testing program with undrained cyclic triaxial tests has been conducted on reconstituted kaolin clay (liquid limit w L ¼ 47:2%, plastic limit w P ¼ 12:2% and viscosity index I v ¼ 1:5%). All samples were consolidated out of a slurry (water content of w=w L ¼ 2:5). Dry kaolin powder was mixed with demineralized water under vacuum for at least 12 hours. The slurry was pre-consolidated in a consolidation apparatus. The axial loading during this pre-consolidation was applied pneumatically. The pre-consolidated samples had a diameter of d ¼ 130 mm and a height of h % 130 mm. Afterwards, triaxial samples (d ¼ h ¼ 50 mm) were cut out of the centre of the preconsolidated sample. In several test series the overconsolidation ratio OCR, the initial mean effective stress p 0 , the initial stress ratio g 0 ¼ q 0 =p 0 , the stress amplitude q ampl and the displacement rate _ s have been varied. If not varied, the cyclic loading was applied with a constant displacement rate _ s ¼ 0:1 mm/min. A more elaborated description of the triaxial tests as well as a detailed examination of the results can be found in [78,80]. In the following, only those parts of the results required for the development of the constitutive equations are presented.
The interval between the cycles 1 and 5 has been used for the evaluation of the rate of pore pressure accumulation _ u acc , i.e. _ u acc ¼ ½u acc ðN ¼ 5Þ À u acc ðN ¼ 1Þ=4. Since the average stress and the strain amplitudes continuously change with the number of load cycles N during the undrained cyclic loading, the analysis is restricted to the early stage of the tests. Compared to sand, where data from drained cyclic tests with almost constant values of r av and e ampl are usually available, the calibration of the HCA model for clay based on undrained cyclic test data is thus more challenging. It should be noted that the data of the first cycles is used to derive suitable equations and to undertake a first manual calibration of the material constants. This applies to all factors f t discussed in this and the following subsection. A further inspection of the equations and a fine-tuning of the material constants is done based on the element test simulations presented in Section 4.5.

Intensity of accumulation
Compared to sand, additional functions to take into account the overconsolidation ratio OCR (f OCR ) and the loading frequency (f f ) have to be added to the intensity of accumulation for the clay model. On the other hand, one of the functions f e and f p incorporated in the HCA model for sand can be omitted since in normally consolidated soils void ratio and pressure are directly interrelated to each other via the normal compression curve e(p). It has been decided to keep f e . The following multiplicative approach is used for the intensity of strain accumulation: In addition to the aforementioned modifications, the factor f Y of the model for sand has been replaced by f g in Eq. (3).
As has been shown in [78], the definition of this factor in terms of normalized average stress ratio g av is more convenient than using Y av employed for the HCA sand model.

Influence of the strain amplitude (f ampl )
The influence of the strain amplitude on the intensity of accumulation and the function f ampl describing this influence is inspected in Fig. 2, where the rates of pore pressure accumulation _ u acc measured during the initial phase of the triaxial tests are shown. Two different test series with different initial stress ratios g 0 = 0 and g 0 = 0.375 are considered. The initial mean effective stress was p 0 = 200 kPa in both test series, while the stress amplitudes lay in the ranges 30 kPa q ampl 70 kPa (g 0 = 0) or 15 kPa q ampl 40 kPa (g 0 = 0.375), respectively. The rate of pore pressure accumulation _ u acc and the intensity of strain accumulation _ e acc are proportional to each other ( _ u acc ¼ K _ e acc m v ). Since all tests within a series have been started at the same isotropic stress (p 0 = 200 kPa, g 0 = 0 or 0.375), the bulk modulus K and the volumetric component of the direction of accumulation m v are assumed constant within the first cycles. Therefore, _ u acc can be analysed instead of _ e acc . Figure 2 shows _ u acc as a function of the strain amplitude e ampl (mean value over the range N = 1 -4) for both the tests with isotropic (g 0 = 0, Fig. 2a) and anisotropic (g 0 = 0.375, Fig. 2b) initial stresses. Compared to sand, where values below 10 À3 are typically encountered, under similar initial and cyclic stress conditions, the strain amplitudes are much larger for the kaolin clay. Therefore, it is meaningful to use a larger reference strain amplitude e ampl ref ¼ 10 À3 in the HCA model for clay than in the sand model, where e ampl ref ¼ 10 À4 has been used. The increase in _ u acc with e ampl almost follows a square root function. If the data are fitted by the same exponent C ampl = 0.51 is obtained for both test series. The current study indicates that the amplitude dependence is thus much less pronounced than in case of sand for which values in the range of C ampl ¼ 1:3 to C ampl ¼ 2:4 are usually found [81]. Based on tests performed on sand, a suitable function to describe the increase in accumulated strain with increasing number of cycles N was found to be [55]: C N1 , C N2 and C N3 are material parameters. The function is adopted for the clay model as well. The required derivative of Eq. (5) with respect to N is given by: As has been demonstrated in [55], the definition of _ f N in Eq. (6) contradicts the Miner's rule. A certain number (say 10,000) of cycles with a strain amplitude tending to zero should lead to zero accumulated strain and have no effect on the accumulation of strain during a subsequent cyclic loading with larger strain amplitudes. However, this is not the case when using the formulation of Eq. (6), since the cyclic loading is solely quantified by the number of cycles, ignoring the corresponding strain amplitude. To resolve this, the intensity of cyclic loading has to be considered in the state variable quantifying cyclic preloading. Therefore, a so-called historiotropic state variable g A was introduced in the original model by Niemunis et al. [55], replacing f ampl Þ considers both the number of applied cycles and their corresponding intensity by f ampl . The same definition of the function _ f N as in the model for sand is used in the model for clay (for a detailed derivation of the equation the interested reader is referred to [55]). _ f N is split into an N-dependent portion _ f A N and a constant portion _ f B N .: with: Note that for a non-constant strain amplitude e ampl ðNÞ (and hence non-constant function f ampl ) the definition of _ f A N becomes [66]: N 0 is the number of the cycle at which the strain amplitude changed the last time. The rate of the preloading variable g A is defined by: The parameters C N1 , C N2 and C N3 are determined based on the back-analysis of element tests presented in Section 4.5.

Influence of the void ratio (f e )
Instead of f e and f p used in the HCA model for sand, the combined influence of void ratio and pressure is described solely by the void ratio function f e in the model for clay, considering that e and p are directly correlated in normally consolidated clays by the virgin compression curve. Tests with different initial void ratios e 0 obtained after consolidation at different mean effective initial stresses within the range 50 kPa p 0 300 kPa have been performed for this purpose, with g 0 ¼ 0 and a constant ratio of stress amplitude and initial pressure q ampl =p 0 = 0.2. Figure 3 displays the results of these tests. The influence of the strain amplitude has been eliminated by dividing _ u acc through f ampl . In addition, since the difference in mean effective stress also influences the bulk modulus K, the normalized rates _ u acc =f ampl have to be further divided by K, which corresponds to a conversion of the rates of pore pressure accumulation _ u acc into rates of volumetric strain Fig. 3 Calibration of function f e based on the undrained cyclic tests with different initial pressures p 0 Fig. 4 Calibration of function f g by means of the data from undrained cyclic triaxial tests with different initial stress ratios g 0 lying in either the compression or extension regime of the p-q plane accumulation _ e acc v ¼ _ u acc =K. The pressure-dependent bulk modulus has been calculated using: with atmospheric pressure p atm . Equation (12) has been derived from the unloading curve of oedometric compression tests reported in [78], considering that the stress relaxation caused by a cyclic loading represents a kind of unloading as well. The slight increase in _ u acc =ðKf ampl Þ with e visible in Figure 3 is captured by with C e = -0.97. The reference void ratio has been arbitrarily chosen as e ref = 1.25, lying in the middle of tested void ratios (in the model for sand the void ratio e max corresponding to the loosest possible state is applied). Obviously, an approximation of the data in Fig. 3 by a linear function would fit as well but Eq. (13) has been chosen in analogy to the HCA model for sand.

Influence of the stress ratio (f g )
The influence of stress ratio is investigated in Fig. 4, based on the data from tests with different initial stress ratios in the range À0:5 g 0 0:625, i.e. comprising the compression as well as the extension regime. All tests within this series were performed with p 0 = 200 kPa and q ampl = 30 kPa. The average stress ratio g av used on the horizontal axis in Fig. 4 is the average value of the g values during the first five cycles. M denotes the critical stress ratio in compression or extension, respectively, which has been calculated with a critical state friction angle of u c ¼ 19 (for it's determination see [78]). In general, for a certain average stress ratio g av , both an accumulation of pore water pressure and an accumulation of deviatoric strain take place. To analyse a single quantity, the rate of pore pressure accumulation _ u acc has been converted into a volumetric strain rate _ e acc v , before the total strain rate _ e acc is evaluated: with p av being the mean effective stress corresponding to g av . The intensity of accumulation calculated from Eq. (14) has been plotted versus the normalized stress ratio g av ¼ g av =M in Figure 4. The function f g describes the influence of the average stress ratio and is chosen as with the material constant C g . The data in Fig. 4 reveal that different parameters C g are necessary in order to fit the compression or extension test data. Therefore, an introduction of Lode's angle h into f g seems recommendable.
For such extension, however, the data from additional hollow cylinder tests with cycles applied at different h values would be advantageous (i.e. not only the special cases triaxial compression and extension should be tested).

Influence of the overconsolidation ratio (f OCR )
The necessity of a function describing the influence of the overconsolidation ratio OCR is demonstrated in Fig. 5a. Compared to the normally consolidated samples (OCR = 1, data overtaken from Fig. 3), the pore pressure accumulation rates in the overconsolidated clay (OCR values up to 2.5 have been tested) are much lower, despite having a similar void ratio. The higher OCR, the lower are the _ u acc =ðKf ampl Þ values. The decrease in the intensity of (b) (a) Fig. 5 Pore pressure accumulation rates _ u acc =ðKf ampl Þ or _ u acc =f ampl , respectively, from tests with different overconsolidation ratios OCR versus a) void ratio e and b) OCR. Note that all tests given in plot b) have been conducted with an initial mean effective stress p 0 ¼ 100 kPa for why no normalization with K is required accumulation by an overconsolidation is thus much stronger than expected from the change of the void ratio De caused by the preloading. It is thus impossible to describe the reduction in the cumulative rates with increasing overconsolidation by the void ratio function f e alone. Consequently, an additional function describing the OCR influence is indispensable. Figure 5b shows the normalized pore pressure accumulation rates from the initial phase of the tests with different overconsolidation ratios plotted versus OCR. An exponentially decreasing function has been chosen to describe the relationship: with a material constant C OCR = 1.1 derived from the data in Fig. 5b. It has to be considered, however, that OCR ¼ p preload =p increases during the undrained cyclic test, due to the relaxation of p. In Figure 5b the initial values OCR 0 ¼ p preload =p 0 have been plotted on the abscissa, which seems reasonable for the initial phase of the tests, where p av is not far from p 0 . For anisotropic average stresses a more general definition of OCR is necessary. For the time being the definition used in the isotropic visco-hypoplastic model of Niemunis [53], which originates from the MCC model, is adopted for the HCA model. The preloading surface is an ellipse described by with the critical stress ratio M. The equivalent pressure p þ e marks the intercept of the preloading surface with the paxis. The overconsolidation ratio is defined as where p e is calculated from the virgin compression line as observed in isotropic compression tests: e e0 and p e0 are the referential void ratio and pressure, respectively, and k is the inclination of the compression line in a diagram with log-log scale. The question, if Eq. (16) should be applied with the actual overconsolidation ratio OCR or its initial value OCR 0 is revisited based on element test simulations in Section 4.5.

Influence of the loading frequency (f f )
The influence of the loading frequency on the accumulation rate and thus the necessity of a function f f is inspected in Fig. 6. The rates _ u acc from the initial phase of the tests with different displacement rates performed with p 0 = 200 kPa, g 0 = 0 and q ampl = 45 kPa are plotted versus the displacement rates. The _ u acc data have been divided by the amplitude function f ampl (known from Fig. 2) in order to purify them from the influence of different strain amplitudes. No clear dependence of _ u acc =f ampl on the displacement rate and thus on the loading frequency can be detected in Fig. 6. Therefore, the respective function can be set to f f ¼ 1 for kaolin. For materials with higher plasticity subjected to cyclic loading such as reported in [75], however, a function f f describing an increase in the cumulative rates with decreasing loading frequency will be indispensable.

Influence of the loading direction with respect
to the sedimentation direction  Comparison of normalized pore pressure accumulation rates _ u acc =f ampl from tests on triaxial samples cut out vertically or horizontally from the pre-consolidated kaolin kPa, g 0 = 0.5 and q ampl = 30 kPa. After the influence of strain amplitude has been eliminated by the division through f ampl , the data show almost no influence of the cutting direction. Therefore, the effect of anisotropy needs not be considered in Eq. (3). It should be captured, however, in the low-cycle model used to evaluate the field of the strain amplitude. In the present case the AVISA model is used for that purpose.

Direction of accumulation
The definition of the direction of accumulation m of the MCC model is also used for the HCA model for clay: where t ! ¼ t=k t k denotes the normalization of a tensorial quantity and t Ã ¼ t À trðtÞ=31 denotes the deviatoric part of a tensor. 1 is the second-order unit tensor. The suitability of this definition for m is discussed in [78,80]. Therein an alternative anisotropic definition of m is investigated as well. It was concluded that the isotropic definition of m is sufficient.

Incorporation of plastic strain and hypoelastic stiffness
Analogous to the HCA model for sand, incorporation of a plastic strain rate is required since the stress can take unphysical values otherwise (e.g. far outside the failure locus). The isotropic yield surface of the MCC model introduced in Eq. (17) is used for this purpose. Using an associative flow rule, the plastic strain can be calculated by where _ / [ 0 is the rate of plastic potential and F is the scalar value of the yield surface in Eq. (17). oF or is equivalent to m introduced in Eq. (20). The plastic strain is calculated iteratively using Newton's method. The mathematical procedures are described in Appendix A.
A simple isotropic stiffness is used for E in Eq. (1) where the scalar factors are defined by: In Eq. (23) j is the swelling index and m is the Poisson's ratio. I in Eq. (22) is the symmetric fourth-order unit tensor. As has been outlined in [55], a hyper-elastic definition of stiffness is inessential for the HCA model since such models do not describe the course of stress and strain during individual cycles. Thus, no accumulation of artificial strain during closed stress cycles or vice versa occurs despite the hypo-elastic definition of stiffness.

Final set of equations and calibration of the parameters
The definition of the functions of the factors f t used for the calculation of the intensity of strain accumulation (Eq. (3)) and their parameters are summarized in Table 1. The model needs a total of 7 material constants and 2 reference quantities (apart from the parameters for the direction of accumulation and elastic stiffness). For the calibration, undrained cyclic triaxial tests with variations of  • stress amplitude (parameter C ampl ) • initial void ratio (parameter C e ) • average stress ratio (parameter C g ) • initial overconsolidation ratio (parameter C OCR ) are necessary while keeping the remaining parameters constant within each series. The parameters C N1 , C N2 and C N3 can be finally determined based on the data from all tests.

Element test simulations
All cyclic tests performed on kaolin have been simulated with an element test program. It calculates the curves of accumulated pore water pressure u acc ðNÞ and accumulated deviatoric strain e acc q ðNÞ, using the initial stress (p 0 , g 0 ), the initial void ratio e 0 , the preloading pressure p preload and the measured strain amplitudes e ampl ðNÞ (changing from cycle to cycle) as input. The overconsolidation ratio entering f OCR was either evaluated as OCR ¼ p e =p þ e with p þ e continuously updated for the actual stress state or simply set to the initial value OCR 0 . In the first case, OCR increases during the undrained cyclic tests also for normally consolidated samples, due to the reduction in p þ e . The latter approach leads to OCR ¼ 1 for normally consolidated samples throughout the test.
Starting from the values estimated from Figs. 2, 3, 4 and 5, the parameters of the intensity functions in Table 1 have been iteratively optimized in the element test simulations. The optimized parameters are summarized in Table 2.
Despite Fig. 4, the same value C g has been applied for the triaxial compression and extension tests. The curves u acc ðNÞ and e acc q ðNÞ obtained from the simulations with the variable OCR are compared to the experimental data in Fig. 8. With the exception of some discrepancies for single tests, a satisfying reproduction of the data from the cyclic laboratory tests can be concluded. The results look similar for the approach assuming a constant OCR (OCR ¼ OCR 0 ). However, different material constants are needed for the different approaches. If the actual OCR is applied in f OCR , the increasing OCR values lead to a decrease in f OCR with N and thus the predicted cumulative rates _ u acc and _ e acc q get lower. This can be compensated by larger parameters C ampl and C N1 , while C N2 has to be reduced ( Table 2).

Back-analysis of centrifuge tests on monopiles in soft clay
In order to show the feasibility of the proposed constitutive model for the study of OWT foundations subjected to lateral cyclic loading, a back-analysis of centrifuge tests performed by Yang et al. [85] on soft kaolin clay is presented in the following. Note that a back-analysis of these tests has been presented in [20] as well, using an elastoplastic bounding surface model and assuming locally undrained conditions.

Centrifuge tests by Yang et al. (2019)
The considered centrifuge tests have been performed at the Zhejiang University in 2019. All quantities are given in model scale in the following. A strongbox with internal dimensions of 1.2 m x 0.9 m x 1.0 m and Malaysia kaolin clay (plastic limit w P ¼ 35% and liquid limit w L ¼ 80%), which is a standard material of the National University of Singapore, have been used. The soil sample was prepared by mixing dry kaolin powder with water using a soil-mixing machine and subjecting the slurry to 8 hours of consolidation in the centrifuge at 100 g. For these tests, water and not a fluid with a higher viscosity has been used as pore fluid. Following the consolidation of the clay, the centrifuge was stopped and the pile was pushed into the soil. The aluminium pile had an outer diameter of 0.059 m, a wall thickness of 0.0022 m and an embedded length of 0.55 m. Once the pile had been installed, the soil including the pile was subjected to 100 g for one hour in order to allow for any settlement around the pile.
Following the resting period, the pile was laterally loaded at 0.23 m above the ground surface. Three cyclic load packages with varying amplitude and average loading with 100 cycles each and a frequency of 0.2 Hz have then been applied. 30 s of consolidation between the packages has been considered. Only the first two cyclic loading packages are considered for the back-analysis since the third package with the largest amplitudes led to large pile head displacements (approximately 1 m in prototype scale) for which a large-deformation numerical method should be applied.

Numerical model of the centrifuge tests
The finite-element analysis was performed using the program numgeo written by the first two authors (Machaček & Staubach, see [43,45,63,69] and numgeo.de). numgeo is a finite-element code specialized on geotechnical boundary value problems. numgeo has been applied for the simulation of vibratory pile driving in saturated sand reported in [45,65] and for the investigation of the influence of the installation on the behaviour of piles subjected to lateral cyclic loading using the HCA model for sand in [63]. The finite element model of the centrifuge test is displayed in Fig. 9. The simulations were performed in model scale. In order to incorporate the change in pore water pressure as well as consolidation effects, u-p elements were used. u-p elements discretize the displacement of the solid phase u and the pore water pressure p w . As has been shown in [62], u-U elements (discretizing the displacement of solid phase u and the displacement of water u w ) could be advantageous in case of low hydraulic conductivities and small time increments since the matrix conditioning of the resulting left-hand side is superior to the one of u-p elements. The well-known oscillations in pore water pressure occurring when using u-p elements are absent using u-U elements. However, since the investigation of the constitutive model is the main target in the current study, the 'conventional' u-p elements are judged to be sufficient for the present purpose. u27p8 brick elements, with 27 nodes discretizing the displacement using bi-quadratic Lagrangian interpolation functions and 8 nodes discretizing the pore water pressure using linear interpolation functions, were utilized.
A surface-to-surface (STS) method was used for the contact discretization. Using the STS technique, the contact stress is determined in each contact surface integration point by calculating the shortest distance to the paired surface by evaluation of the convective coordinate. Due to the bi-quadratic shape functions of the u27p8 element the derivatives are continuous and the normal vector is definable at every point of the surface. Details on the employed STS technique can be found in [67]. A stiffness-dependent penalty factor (30 times the trace of the stiffness tensor of the adjacent continuum) was used to enforce the normal contact constraints. A Coulomb friction model with a friction coefficient of 0.2 was applied.
A bulk modulus of the pore water of 2.2 GPa was assumed. The permeability of the Malaysian kaolin clay is K Perm ¼ g w k w =c w ¼ 2 Á 10 À15 m 2 (assuming the dynamic viscosity of water to be g w ¼ 1 Á 10 À6 kPas). With the increased gravity in the centrifuge test, this resulted in a hydraulic conductivity of k w % 2 Á 10 À6 m/s during the testing.
Installation effects were not taken into account in the simulations since in the experiment the pile has been jacked into the soil at 1 g. Due to the increase in gravity to 100 g, the influence of the installation on the initial soil state is assumed to decrease significantly. The initial overconsolidation ratio prior to the spin-up of the centrifuge was assumed to be one. The lateral stress coefficient was set to 0.5.
For the low-cycle part of the simulations, the anisotropic visco-ISA (AVISA) model [73] was used. This model is a rate-dependent model for fine-grained soils able to reproduce small-strain effects, inherent anisotropy as well as the influence of the overconsolidation ratio. The strain rate dependency is reproduced by incorporating a third strain rate mechanism (in addition to the elastic and hypoplastic strain rate), which can be switched off in case of lowplasticity fine-grained soils. To cover a wide range of strain/stress amplitudes as well as to account for the specific behaviour of clays the ISA plasticity approach of [24] has been revised and slightly reformulated in [73]. Furthermore, to account for initially anisotropically consolidated states (induced anisotropy) the kinematic hardening mechanism of [24] has been modified according to the experimental findings gained in [78]. In addition, a loading surface has been incorporated to define a threedimensional overconsolidation ratio and to account for its effects on the accumulation behaviour. To capture the inherent anisotropy exhibited by some clays, a transversal hypoelastic stiffness has been further introduced into the mathematical formulation of AVISA. The model is also able to reproduce the cyclic tests presented in Figure 8 well, as shown in [73], and is thus technically applicable for the complete simulation of the centrifuge tests (without the HCA model). However, the computational effort necessary to simulate all cycles conventionally exceeded the available computational resources. Therefore, the application of the HCA model seems compulsory for the present simulations.
The parameters of the AVISA model have been partially determined for the Malaysian kaolin clay ðw P ¼ 35%; w L ¼ 80%Þ used in the centrifuge tests (parameters of the MCC model: k, j, m h and M c obtained from [42]) or have been assumed to lie between the parameters calibrated for the kaolin clay considered in the previous sections ðw P ¼ 12:2%; w L ¼ 47:2%Þ and those for the more plastic Lower Rhine Clay ðw P ¼ 34%; w L ¼ 56:1%Þ reported in [71,73]. The utilized set of parameters is given in Table 3, whereby further details about the correlations used for the calibration as well as the numerical calculation of an oedometric compression test are provided in Appendix B. It is worth mentioning that using the HCA model for clay, the influence of the conventional model is much lower than in case of the HCA model for sand since the influence of the strain amplitude, calculated by the intrinsic model, is much lower (see lower C ampl values discussed in Section 4.1.1).
The soil inside the pile has been modelled elastically since it is only involved in a rigid deformation when the pile rotates and does not influence the pile response. Preliminary simulations showed that an elastic stiffness of 10,000 kPa and a Poisson's ratio of 0.3 are appropriate values.
The parameters of the HCA model for clay determined in the previous Section given in Table 2 (MCC flow rule with OCR 0 ) were applied for the simulations. The same reference values were used as well. For the calculation of f g and the direction of accumulation the critical friction angle used in Sect. 4.1.4 was employed.
During the high-cycle phase of the HCA model, instead of introducing update cycles as shown in Figure 1, the change of the strain amplitude e ampl due to the change in the soil stiffness caused by the cyclic loading has been taken into account using a so-called adaptive strain amplitude definition, which has been recently proposed in [66]. The strain amplitude is periodically updated by linearly scaling it with respect to the change in soil stiffness compared to the stiffness at the end of the second load cycle. A non-local smoothing algorithm is used to avoid localization of large values of strain amplitude in some integration points due to self-reinforcement effects similar to the approaches used in non-local plasticity. For a detailed description of this procedure, the interested reader is referred to [66].
The simulations were performed in the following steps:  deformations due to the cyclic loading were predicted by the HCA model. Taking into account the frequency of the cyclic loading of 0.2 Hz, 500 s were simulated in the high-cycle phase. 6. 30 s of consolidation were allowed before application of the next cyclic package. Only H av has been applied in this phase. 7. Repetition of steps 3-6 using the load magnitudes of the second package of cyclic loads (H av ¼ 110 N and H ampl ¼ 65 N).

Results of the back-analysis
As has been done in the study by Yang et al. [85], the results are evaluated in prototype scale. The pile head displacement measured at the point of load application is given in Figure 10 for the measurements as well as for the results of the simulation. The results of the experiment are displayed as mean pile head displacement beginning after the end of the second load cycle. The measured displacement after the first two cycles is slightly underestimated by The subsequent HCA phase starts at t ¼ 1; 300 s (N ¼ 2) and first predicts a very fast accumulation which slows down after approximately % 1; 500 s. Compared to the experiment the accumulation rate of permanent pile head displacement is slightly too small in this latter phase. The pile head displacement predicted at the end of the first cyclic load package fits well to the measurements, however. For the second load package two individual cycles are simulated again using AVISA. Due to the larger loading magnitudes, larger displacement amplitudes are observed compared to the cyclic loading of the first package (see plot on right-hand side of Fig. 10). Consequently, the strain amplitude is larger leading to a higher accumulation rate of the HCA model in the subsequent high-cycle phase. The accumulation of pile head displacement predicted by the HCA model during this phase fits well to the values recorded in the experiment.
The curves of bending moment versus the height of the pile after N ¼ 1 and N ¼ 90 cycles are given in Figure 11. The mudline is located at a vertical coordinate of 0 m. The distribution after one cycle shows that the simulation underestimates the magnitude of the bending moment slightly but captures the general shape of the measured curve very well. After the application of 90 cycles both the bending moment measured in the experiment and the one obtained from the simulation have increased in magnitude. This increase is more pronounced in case of the simulation. The accordance at N ¼ 90 between the measurements and the results of the simulation is almost perfect. Figure 12 displays the spatial distribution of the strain amplitude e ampl at N ¼ 2 and N ¼ 90 for the first cyclic load package (top part of figure). The largest strain amplitudes with values above 0.1 % occur at the soil surface in the vicinity of the pile. Thus, the function f ampl is largest at this location. As has been explained earlier, the change in the strain amplitude is taken into account using an adaptive strain amplitude definition. The spatial distribution of the strain amplitude following the last update performed during the first cyclic load package (at N ¼ 90) shows that the strain amplitude has changed only slightly due to the application of 90 cycles. Moderately larger values are obtained at the left-hand side of the pile (the average loading of the pile is applied to the right) compared to the distribution at N ¼ 2. The strain amplitude is much larger for the second cyclic load package as is visible from the lower part of Fig. 12. During the second cyclic load package, the change of the strain amplitude is larger compared to the first cyclic load package as is visible from the distribution at N ¼ 90. Therefore, the consideration of the change in strain amplitude in the simulation with the HCA model is judged mandatory for the cyclic loading with larger amplitudes but is of less importance for the cyclic load package with lower load magnitude. This is reasonable since the cyclic loading with higher magnitude causes larger increase in excess pore water pressure and   figure) for the first cyclic load package consequently a more pronounced reduction in soil stiffness leading to larger strain amplitudes.
The function f g shows large values alongside the lefthand side of the pile at the start of the HCA phase as is visible from Fig. 13. Since the pile is loaded to the right, the horizontal effective stress is reduced on its left, which leads to large average stress ratios g av . From the factors influencing the rate of strain accumulation, f g varies the most in magnitude over the calculated domain. While the factor f e takes values in the range of 0.8 to 1.3, f g varies from 0 to values above 60. f ampl has a comparably modest influence with values in the range of 0 to 2 (for the first cyclic load package). f OCR takes values reaching from 1 to 2 due to the change in the overconsolidation ratio prior to the HCA phase caused by the first two cycles of loading.
The fields of excess pore water pressure ratio (Dp w =p 0 ) at the beginning of the HCA phase (at N ¼ 2) and at the end of the simulation (at N ¼ 100) are given in Figure 14 for the first cyclic load package. At N ¼ 2 an increased pore water pressure in the upper half of the soil on the right-hand side of the pile (where the soil is pushed by the pile) is evident, while a decrease is observed in the soil on the left-hand side (where the soil is unloaded). After the application of 100 cycles, a significant increase in the soil area at the upper part of the pile affected by large excess pore water pressures is observed. The increase in excess pore water pressure is due to the tendency of the soil to compact around the pile when subjected to cyclic loading. The consolidation process taking place simultaneously is, however, not fast enough to allow free drainage of the pore water.
6 Summary, conclusions and outlook A HCA model for clay which predicts the accumulation of permanent strain under cyclic loading with a large number of cycles has been presented. The data of an extensive laboratory testing program on kaolin clay under undrained cyclic loading has been analysed for this purpose. Compared to sand under cyclic loading, a lower influence of the strain amplitude on the strain accumulation was found. The rates of strain accumulation were found to increase with increasing void ratio and increasing stress ratio. As expected, samples with larger initial overconsolidation ratios showed less accumulation. No significant influence of the frequency of loading on the rate of accumulation was observed for the studied kaolin, due to its rather low plasticity.
Building upon these results, constitutive equations of the HCA model for clay have been proposed. The concept of the constitutive model is based on the HCA model for sand by Niemunis et al. [55]. The model has then been validated by the simulation of element tests and the back-analysis of a centrifuge test on a monopile in soft clay under cyclic lateral loading. The results of the simulation of the centrifuge test are in good accordance with the measured data regarding the development of pile displacement with the number of applied cycles. The measured bending moment of the pile after N ¼ 1 and N ¼ 90 cycles could also be captured accurately by the simulation.
Future work will focus on the further development and validation of the constitutive equations for clays with higher plasticity. An influence of the loading frequency on the accumulation rate is expected in this case. In addition, other relations for the overconsolidation ratio as for example the equivalent time line model proposed by [38] could be adopted. In future work it will be investigated whether the incorporation of this concept can cover both the overconsolidation effects as well as the frequency dependence of plastic fine-grained soils. The concept of [38] can even be generalized for three-dimensional stress states considering the formulation developed in e.g. [73] and cover the influence of the stress ratio as well. Thus, the factor f g could maybe be omitted. Since the proposed model was developed on the basis of tests on reconstituted samples only, it is intended to study the accumulation behaviour of undisturbed natural samples in future work, where effects of destructuring have to be additionally considered.
The proposed model can also be used for a parametric study on monopile foundations in clay similar to the study considering sandy soil [68] and up to 1 million load cycles. In addition, the influence of the installation of the pile prior to the cyclic lateral loading should be studied. Investigations considering sand showed that the installation induced changes in the soil state influence the pile response to cyclic lateral loading considerably [64]. So far, it is not clear how large the installation (method) influences the pile response to subsequent loading in cohesive soils.

Calculation of plastic strain
The increment in plastic strain is calculated by using the definition of the yield criterion F given in Eq. (17). D/ is determined by computing the zero of Eq. (17) using Newton's method: Equation (25) is solved in every iteration i until F is sufficiently close to zero. Using the chain rule, the derivative of F with respect to D/ is obtained: The partial derivatives of p, q and p þ e are (see [11]): In Eq. (29) # is defined by: Calibration of AVISA parameters The AVISA model requires the calibration of 15 parameters: • Four parameters are needed for the description of the transversal (hypo)elastic stiffness: the compression index k, the swelling index j, the Poisson's ratio m h and the anisotropic coefficient a. k, j and m h have been adopted from the MCC model parameters reported for Malaysia kaolin clay in [42]. The accuracy of these values for the AVISA model has been proven by simulating the oedometric compression test from [42] with one unloading-reloading cycle, as shown in Fig. 15. a has been chosen to lie between the values reported for the kaolin considered in the present study and the highly plastic Lower Rhine Clay reported in [71]. • For the determination of the critical as well as the loading surface, three parameters are required: the triaxial compression critical state slope M c ¼ 6 sinðu c Þ=½3 À sinðu c Þ adopted from the MCC value reported in [42]; the loading surface factor f b0 determined in the same manner as a and the maximum void ratio e i0 at p ref ¼ 1 kPa calibrated by means of the empirically derived correlation function e i0 ¼ 0:9 q s =q w [71] with q s ¼ 2:6 g/cm 3 from [85]. • One parameter is responsible for the rate dependency of the model, the so-called viscosity index I v defining explicitly the creep intensity. It has been assumed identical to the value calibrated for kaolin clay. However, a more accurate calibration of this parameter could improve significantly the predictive capabilities of the model. • The seven parameters which define the small-strain stiffness have been assumed in the same manner as a.
Since only two cycles are simulated with AVISA, these parameters are of secondary importance and could even be switched off due to the hierarchical composition of the model.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons. org/licenses/by/4.0/.