Soil-structure interaction with time-dependent behaviour of both concrete and soil

As the consideration of soil-structure interaction is increasingly being incorporated into design practice, there is a need for solutions that consider the most relevant aspects of the behaviour of both the concrete and the soil, some of which are time-dependent. The present work introduces a model that suits both the structure and the foundation – including the supporting soil – allowing the coupled analysis of hardening on cure, creep, shrinkage and cracking of concrete, and consolidation of soil. Both the structure and the foundation are modeled as one-dimensional finite elements. The time-dependent behaviour of the concrete and soil is modeled using Kelvin chains. For the structural elements, a mixed finite element formulation is used. Validation tests were conducted in the proposed modeling, comparing its results with available experimental and analytical data. The study of a reinforced concrete continuous beam supported by foundations on consolidating clay considering the time-dependent behaviour of the materials showed considerable changes in the effects of the soil-structure interaction.


INTRODUCTION
The interaction between the structure and its foundation is increasingly being considered and incorporated into design practice.Structural projects based on rigid supports are still made for small buildings or in cases where the foundations are certainly rigid (not settling).In structures founded on footings or rafts, or even on groups of piles that interact with each other (i.e.subject to group effect), the analysis including soil-structure interaction leads to different and more realistic results.In addition to changing the pattern of settlements, making them more uniform, it is possible to assess the redistribution of forces in the structure.
Foundation and structural standards around the world are beginning to include recommendations for interaction analyses both in working conditions and at failure, although the former is more usual.For the analysis, adequate models must be used, ideally in a coupled solution.Although being recognized as important factors in the behaviour of reinforced concrete structures, concrete shrinkage, hardening (aging of the elastic modulus), cracking and creep, as well as the presence of reinforcement are hardly included in soil-structure interaction analyses.
The main parameter to assess the implications of an interaction analysis is the relative stiffness between the structure and the foundation (Lee and Brown, 1972;Gusmão, 1990).Several studies were carried out to analyze the influence of different factors affecting this stiffness, such as three-dimensional effect (Majid and Cunnell, 1976), number of floors (Goschy, 1978;Lopes and Gusmão, 1991), construction sequence (Brown and Yu, 1986), soil nonlinearity (Jardine et al., 1986), soil constitutive model (Noorzaei et al., 1995), interface elements (Dalili et al., 2014), and even the presence of wall columns (Dalili et al., 2011) and brick panels (Hora, 2014).
Various studies considered the time-dependent soil behaviour (Viladkar et al., 1993), the viscoelastic behaviour of the structure (Gonçalves et al., 2016;Rosa et al., 2018) or both (Santa Maria et al., 1999;Fajman and Sejnoha, 2007).Nevertheless, no study could be found that also considered the effects of concrete hardening, cracking and tensionstiffening (increase in stiffness of a concrete section due to the transmission of stresses from rebars to the concrete between adjacent cracks) and the presence of reinforcement, all of which considerably influence the response of the structure (Espion and Halleux, 1990).
The present work presents a model that can be applied to both the structure and the foundation -including the soil -allowing a coupled analysis.The model allows the consideration, for the concrete, of hardening, shrinkage, creep, reinforcement, cracking and tension-stiffening.As for the soil, it enables the consideration of consolidation (optionally also secondary consolidation).The modeling uses, for both the structure and the foundation, one-dimensional finite elements (truss or beam elements).Kelvin chains are used to model the time-dependent behaviour of the materials, namely creep and consolidation.A mixed finite element formulation is used for the structural elements.
Validation tests were carried out after the models were programmed in a code referenced here by the aceonym SoilStruc, comparing its results with available experimental and analytical data.Then, a hypothetical case is studied.The case consists of a continuous beam supported by foundations on consolidating clay.Although simple, this case allows the examination the effects of the various phenomena on the behaviour of structure, alone and in combination.
This study shows that, although concrete hardening and reinforcement increase the stiffness of the structure, creep and cracking reduce it more significantly, causing both the uniformity of settlements and the redistribution of forcesexpected in a soil-structure interaction analysis -to be less significant.

PROPOSED MODELING FOR TIME-DEPENDET BEHAVIOUR OF REINFORCED CONCRETE BEAMS
In buildings, continuous (hyperstatic) beams are the main elements responsible for the redistribution of forces between columns when differential settlements occur.Therefore, it is crucial that its behaviour be adequately modeled for a correct force-displacement relationship of the structure.
Although many studies on soil-structure interaction consider the behaviour of structural elements as linear elastic (e.g. Brown and Yu, 1986;Masih, 1993;Viladkar et al., 1993;Arapakou and Papadopoulos, 2020), it is known that reinforced concrete beams tend to crack and develop a non-linear response (Cosenza and Greco, 1990).In addition, due to shrinkage on cure and creep of the concrete, the force-displacement relationship varies over time (Ghali, 1986), increasing the complexity of the analysis.
Building structures are usually modeled using one-dimensional elements.Due to the slenderness of these elements, the strains due to shear can be neglected and the Euler-Bernoulli hypothesis can be assumed.In this way, all strains are null, except the axial strain which is given by: where   is the strain at the cross-sectional reference point,  is the vertical distance from the fiber to the cross-sectional reference axis,  is the cross-sectional curvature,   and  are defined as: Latin American Journal of Solids and Structures, 2022, 19(5), e453 5/21 Besides the classical parameters, such as the elastic modulus of the material, geometric properties of the crosssection and element length, the phenomena that most affect beam strain are concrete shrinkage, hardening, creep and tension-stiffening (Espion and Halleux, 1990).Such phenomena can be modeled at different levels.While the first three can be modeled at the material level, the latter has to be modeled at the cross-section level.With the relationship between internal forces and cross-sectional strains, it is sufficient to integrate the strains to obtain displacements at the extremities of the element, disregarding rigid body movements.

Constitutive laws of the materials (a) Concrete
For a uniaxial condition, under constant stresses ( 0 ) applied at  0 , the total strain of the concrete () at time  can be given by the sum of four parts: where   ( 0 ) is the instantaneous strain due to loading,   () is the creep strain,  ℎ () is the shrinkage strain, and   () is the thermal strain.
For the in-service behaviour, the relation ( 0 ) x   () can be considered linear.Furthermore, assuming there is no temperature variation, (Eq.4) can be written as: where (,  0 ) is the creep function.
For a continuous loading function, the strain generated by the stress history is given by: The step-by-step method was selected to model creep as it is applicable for evaluating any stress history.In this method, creep can be represented either in integral or in differential form.The second representation was adopted because it is computationally more efficient as it is not necessary to store the entire stress history, but only a fixed number of variables that are updated at each step (Bazant and Jirásek, 2018).
The creep function can either be obtained in tests or taken from standards and recommendations.However, for the differential formulation, the creep function needs to be properly represented.For the case of hardening viscoelastic materials, Bazant and Jirásek (2018) suggest approximating the function by the following Dirichlet series: where   is the age-independent modulus, m is the number of Kelvin units in a chain,  is the index referring to the Kelvin unit in a chain,   is the parameter of μth hardening Kelvin unit,   is the μth retardation time and  is the Heaviside step function.
The least squares method can be applied to obtain the values of the   functions.To this end, the values of the retardation time   , the age at loading  ′ and the time increment  −  ′ must be chosen properly to ensure both that the system is not ill-conditioned, and that the solution has good accuracy (Bazant and Wu, 1973).The value of  is a function of the time interval for which the viscoelastic response is desired.The longer the range, the greater the number of Kelvin chains needed.
In the numerical solution, the problem is treated incrementally, and time is divided into intervals of size ∆  .Integration algorithms lose numerical stability or accuracy when such an interval is greater than the shortest retardation time of the Kelvin chain.Using the exponential algorithm for hardening materials proposed by Bazant (1971), this limitation is overcome, allowing greater flexibility in defining time intervals.
The incremental constitutive relation using the second order exponential algorithm in the time interval ∆  =  +1 −   can be written as: Latin American Journal of Solids and Structures, 2022, 19(5), e453 6/21 where the incremental modulus  �  , the incremental creep strain at constant stress Δ  ′′ and the incremental shrinkage strain Δ ℎ, are given by: (10) where  ,+1/2 is the mid-step stiffness modulus, and the factors  , ,  , and  , are given by: where  , is an internal variable representing the dimensionless strain rate, whose initial value is null and whose updated value at the end of each step is given by: The functions proposed by Comite Euro-International du Beton (CEB-FIP) in Model Code 1990(1993) were used to model the increase in the elastic modulus (hardening), creep, and shrinkage.

(b) Steel
The behaviour of steel was considered as linear elastic for both compression and tension, as it is assumed that for working conditions, the yield stress is not reached.Therefore, the increase in stress in the steel Δ  is given by: where   is the elastic modulus of the steel and Δ  is the steel strain increment.

Constitutive law of the cross-section
For reinforced concrete elements in working conditions, the tension region can either be cracked or not.So, the section response is divided into non-cracked and cracked.In the first one, it is assumed that concrete withstands up to the tensile strength.In the second one, it is assumed that the tensile strength is null.Furthermore, it is well known that the concrete in tension between cracks has a significant contribution to the stiffness of the element (Beeby et al., 2005).This phenomenon is known as tension-stiffening.To consider such phenomenon, the proposal by Fédération Internationale du Béton (fib) in Fib Model Code for Concrete Structures 2010 (2013) was adopted.In the proposed formulation, the strain of the section at instant k is given by: where  1, is the non-cracked cross-sectional strain vector,  2, is the cracked cross-sectional strain vector, and  is an interpolation coefficient given by: Latin American Journal of Solids and Structures, 2022, 19(5), e453 7/21 where  is a coefficient that considers the influence of the loading time on the mean strain,   is the cracking moment and  represents the bending moment.
According to Beeby et al. (2005), the value of 0.5 for  can be used in the calculation of creep strain, since the loss of tension-stiffening occurs in a shorter period of time when compared to creep.To perform the analysis over time, Eq. ( 16) must be written in incremental form.Using the Taylor series expansion and removing the higher order terms, it follows: where Using the Euler-Bernoulli hypothesis, it is possible to define an expression relating the increase in internal forces with the increase in the strain of the section.In the case of viscoelastic materials, such expression is: where Δ �  is the increment of cross-sectional internal forces vector, Δ �  ′′ is the increment of constrained cross-sectional internal forces vector due to creep and shrinkage and  ̅ , is the stiffness matrix of the transformed cross-section.For the case of a section with n materials, the values of Δ �  ′′ and  ̅ , are given by: where Δ , ′′ is the increment of cross-sectional strain due to creep and shrinkage,  , is the cross-sectional area,  , is the section first moment of area and  , is the section second moment of area.
Considering the same stress increment Δ �  , the strain of the non-cracked section Δ 1, will be different from the strain of the cracked section Δ 2, , since the consideration or not of the concrete tensile strength affects the properties of area  , ,  , and  , of the cross-section.Thus, the response of the non-cracked and cracked section can be given, respectively, by: In the cracked section, it is necessary to define the height of the neutral line to calculate the properties of area.Such analysis can be done using an iterative process (Kawakami and Ghali, 1996).Furthermore, because of both creep and concrete shrinkage, the neutral line position tends to change over time.However, if this height is supposed constant over time, the analysis is simplified and errors resulting from this hypothesis are not significant (Ghali, 1986).In a comparison carried out by Torres (2001), the difference in section curvature was only 2%.

Beam element formulation
The set of equations that must be solved for a beam element are related to the strain-displacement relationship (Eq.1), constitutive cross-section relationship (Eqs. 16,23,24) and equilibrium equations given by: where   and   are the loadings per unit length of the beam in the x and y directions, respectively.
As can be seen in Eq. ( 16), due to tension-stiffening, the response of the cross-section is non-linear when the cracking moment is exceeded.To model such behaviour with finite elements based on the displacement field, it is necessary to subdivide the element so that the nonlinear behaviour of the material is satisfactorily approximated.When using elements based on the force field, such refinement is not necessary and, thus, elements with less degrees of freedom can be used, leading to greater computational efficiency.
In the present work, the mixed finite element formulation proposed by Taylor et al. (2003) was used.After discretizing the fields and calculating the integrals along the element, the equations obtained correspond to the balance of forces, to the displacement at the ends of the element, and to the balance of the internal forces along the element.Putting the three equations in matrix incremental form, follows: where index  indicates the value of the variable at a point of the quadrature,  is the vector of nodal displacements,  is the vector of nodal internal forces,  is the transformation matrix for eliminating rigid body movements,  is the interpolation function of the nodal internal forces,  is the nodal forces vector,   are the internal forces due to loading, and  � are the internal forces due to cross-sectional strain.These variables are given by: where  is the element length,   is the quadrature weight and length, and −1 ≤  ≤ 1. Substituting Equations ( 18), ( 23) and (24) in Equation ( 27), the incremental formulation for a viscoelastic material considering the contribution of the concrete between cracks is: Note that, as the constitutive relationship is linear and the force field along the element is known, for both noncracked and cracked sections it is possible to calculate the respective strain fields along the element and, therefore, the fields related to every element of the Kelvin chain.
The position of the cracked region changes depending on the load level and, thus, the positions of the integration points will also be modified, making it necessary to determine the values of the internal variables at the new points.Since the strain field is known, it is possible to calculate strains at any point.The solution algorithm used is given in appendix.

Solution algorithm
The solution of the structural system is obtained by iterations.The implemented algorithm follows the steps: 1.It is assumed that the process starts after the convergence of the previous step, which is represented by index k.Increment time according to  +1 =   + ∆  , nodal forces to  +1 =   + ∆  and loads along the element to  +1 =   + ∆  ; 2. For each material with viscoelastic behaviour, calculate:  ,+1/2 ,  ,+1 ,  ,+1 ,  � +1 , ∆ +1 ′′ ; 3. For each cross-section, calculate:  ̅ ,+1 , ∆ � +1 ′′ ,  ,+1 ; 4. In order to determine the state of the structure subject to the force vectors  +1 and  +1 , an iterative solution strategy is started aiming at reaching the global equilibrium with an iterative counter  = 0; the structure state in k is used as initial guess; 5.For each element, assuming that all nodal displacements increments are null, calculate:  , ,  , , ∆  ,  , ,  ,, ; 6. Create the stiffness matrix and the resisting nodal force vector and calculate the nodal displacements:

MODELING FOUNDATION SETTLEMENTS BY CONSOLIDATION
The magnitude of foundation settlements is fundamental for the analysis of the soil-structure interaction.In the proposed modeling, settlement development over time of a foundation derives from soil consolidation and a foundation element will be represented by a truss element (vertical) whose response over time is equivalent to that of the foundation.In this manner, the implementation of such analysis in frame-structure programs, which are used in most design offices, is simplified.

Long-term settlement calculation
The final settlement of a normally consolidated clay under 1D conditions can be obtained with Terzaghi's formulation: where   is the compression index,  0 is the initial void ratio,  0 ′ is the initial effective vertical stress, Δ 0 ′ is the stress increase (average in the layer) and ℎ is the layer thickness.
If the stress increment is small in relation to the initial stress, the following approximation can be made so as to linearize the function: In order to use a truss element to represent a consolidating clay, it is necessary to define its area , its length l, and the elastic modulus of the material.The length can be adopted as the clay layer thickness.The area can be defined as a function of the average stress increase in the layer, Δ 0 ′ , of the mean applied pressure , and of the dimensions  and  of the foundation with: The average stress increase in the consolidating layer can be calculated by the Theory of Elasticity.The average stress increase can be calculated from the profile of stress increases obtained with Boussinesq's equation.
Based on Eq. ( 37), the elastic modulus   relative to consolidation can be defined as:

Settlement development over time
The development of settlement by consolidation can be modelled as a viscoelastic behaviour (Fajman and Sejnoha, 2007).For the 1D case, the consolidation equation is: where  is the excess pore-pressure and   is the coefficient of consolidation.Assuming that both upper and lower boundaries are draining and 2ℎ  apart, and that the load application is instantaneous, the settlement   can be calculated by: Comparing Eq. ( 41) with the Dirichlet series (Eq.7), it can be observed that the term   tends to infinity, that is, the material does not present instantaneous responses to loading.The terms   and   can be written as: Latin American Journal of Solids and Structures, 2022, 19(5), e453 11/21 The analogy described above allows the rational determination of the number of elements in the Kelvin chain and their respective parameters to accurately represent soil consolidation.It is interesting to note that the delay time value   decreases with , and that it is the smallest value of   that indicates the capacity of the model to accurately represent the shortest time interval under analysis.Thus, considering that a finite number  of elements in the Kelvin chain will be used in the simulation, it is necessary to introduce the term   , which is given by: As the viscous model used for the soil is the same as that used for the concrete, it is possible to compare the values of the delay time of each model to define how many Kelvin elements should be used and which minimum time step should be initially adopted in the analysis.

Experimental validation of concrete creep
For the experimental validation of creep, Benchmark no. 1 was used, with tests of one span reinforced concrete beams conducted by Subcommittee 3 of the Technical Committee TC 114 of RILEM (Espion, 1993).
The characteristics of the beams and loading can be seen in Figure 1 and Table 1.A total of 7 beams were tested with 5 load levels, 2 of them repeated in order to assess the consistency of the results.The main objective of the tests was to determine the development of the displacement at the middle of the span over time.The beams were cured under constant moisture and temperature, thus, allowing to disregard thermal variation.The loading was applied instantly at 28 days and kept constant until the test was concluded.
Some parameters required for modeling had to be estimated, based on standards and recommendations, as creep tests of the concrete were not available.Thus, a perfect match between the experimental results and the numerical predictions should not expect, but rather the latter should reproduce the response of the structure with acceptable accuracy.
The experimental and numerical results can be compared in Figure 2.For the non-cracked case (C11) or the case with higher loading (C15), there is a better agreement with test results than for the other cases, in which the size of the crack formation region is larger.Overall, it can be concluded that the model provides satisfactory results.

Validation of soil consolidation modeling
In order to verify the accuracy of the method presented in item 3.2, the analytical solution was compared with results using 1, 2, 5 and 7 elements in the Kelvin chain.Figure 3 presents the percentage of settlement as a function of the time factor, T, for the analytical solution (U) and for different numbers of elements in the Kelvin chain (U1, U2, U5 and U7).The parameters were adjusted using Eq. ( 44) as to obtain 100% of settlement at the end of the analysis.Note that for  > 0.1, there is no difference between the models used, whereas for smaller values of T, a significant difference Latin American Journal of Solids and Structures, 2022, 19(5), e453 12/21 between U1 and U can be observed.The difference between U5 and U7 is very small and both models were able to produce the expected response with good accuracy.

STUDY OF A CONTINUOUS BEAM ON SETTLING FOOTINGS DUE TO SOIL CONSOLIDATION
To assess time-dependent effects on soil-structure interaction, a continuous beam with 3 supports on footings was analyzed.This case was chosen as the beams of a structure are the main elements responsible for load redistribution, and for keeping the problem simple in order to evaluate how each of the intervening phenomena influences the problem.

The problem in study
The structure in study is shown in Figure 4a.Circular footings had dimensions defined by for an allowable stress of 50 kN/m 2 , with the loads obtained with rigid supports: 112.5 kN and 375 kN.The ground profile is shown in Figure 4b.The water table is at the base of the sand, and the unit weights of the soils are   = 19 kN/m³ and   = 14.7 kN/m³.The initial vertical effective stress in the middle of the clay layer is   ′ = 113 kN/ m².The stress increase profile in the soft clay layer -obtained with Boussinesq's Equation -is shown in Figure 4b.
To predict settlements due to the consolidation of the clay, the vertical stress increases along its thickness were considered.Performing numerical integration, the average stress increases are 1.91 kN/m 2 and 5.76 kN/m 2 under the side and center footings, respectively.Considering a normally consolidated clay and with   = 1.0 and  0 = 2.5, the final settlements are 1.24 cm and 3.70 cm, respectively.The contribution of the dense sand to the settlements of the footings can be neglected and the clay consolidation settlements will be the (long term) footing settlements.The structural model is shown in Figure 5, where the horizontal element represents the beam and the vertical elements represent the footings (more specifically, the soil behaviour under the footing).The characteristics of the material and cross-section of the beam are shown in Table 2.The cross-sectional area of the side vertical element is 65.1 m 2 and of the center element is 58.9 m 2 , and correspond to the average stress increase in the soft clay layer, i.e., the area of the elements was adopted so that the stresses remain equal to those indicated above (1.91 kN/m 2 and 5.76 kN/m 2 ).

Influence of the viscoelastic behaviour of materials
To assess the effect of the viscoelastic behaviour of materials − namely concrete creep and soil consolidation − 6 cases were analyzed, as indicated in Table 3. E 0 represents the elastic modulus of the material at the initial moment, E ∞ represents the modulus at infinite time, and the viscoelastic material presents a transition from E 0 to E ∞ over time.In relation to concrete, the instantaneous elastic behaviour is represented by cases 1, 2 and 4. Creep is considered in cases 3 and 5; and the elastic behaviour after material creep is considered in case 6.For the soil, the elastic behaviour after consolidation is considered in cases 2, 3 and 6; and consolidation is considered in cases 4 and 5.At the beginning of the analysis, case 3 is equal to 2, and cases 4 and 5 are both equal to 1.It is also noted that there is a smooth transition over time.Hence, at the end of the analysis, case 4 is equal to case 2, and cases 3 and 5 are both Latin American Journal of Solids and Structures, 2022, 19(5), e453 15/21 equal to case 6. Cases 4 and 2 present the same result at the end of the analysis because in both cases the material of the beam is elastic, and that the only difference is the development of the supports' settlement over time.The equality of cases 3, 5 and 6 indicates that the viscoelastic material presents a smooth transition from instantaneous behaviour to the behaviour at infinity time.Therefore, both the behaviour at the beginning and at the end of the analysis can be obtained using the respective elastic modulus, rendering unnecessary an analysis over time.
Knowing that the differential settlement calculated with rigid supports (case 1) is 24.6mm, it can be noted that the consideration of flexible supports (case 2) causes three different reductions: in differential settlement (36%), in the reaction of support (16%) and in the negative (tension in upper fiber) bending moment (81%).Those results were expected since the soil-structure interaction tends to cause more uniform settlements and to transfer loads from central to the external columns.However, when considering the elastic modulus of the concrete in infinite time together with the flexible support (case 6), it can be seen that the reductions in differential settlement (11%), in support reaction (6%), and in bending moment (30%) are much less when compared to case 2. Thus, it can be stated that creep tends to reduce the effects of the soil-structure interaction, since it reduces the elastic modulus of the structure, making it more flexible.
It is interesting to note that the redistribution of internal forces is quite significant and that these forces are used for the design of the beams.Therefore, not considering the interaction can lead to results on the unsafe side.

Influence of the coefficient of consolidation
To evaluate the influence of settlement rate on internal forces, analyzes were performed by varying the value of the soil coefficient of consolidation,   , and comparing them to the cases of elastic and rigid supports.Evaluating the results for the central support (Figure 7), it can be noted that at  = 0, all analyzes indicate the same response as that of rigid supports, as expected.Over time, there is a smooth transition so that by the end of the analysis all results coincide with that of the elastic supports.It can be concluded that the   of the soil does not influence neither the initial nor final response, but rather the way in which the transition between these two states occurs.Furthermore, it seems that the cases of rigid and elastic supports act as an envelope for all other cases.

Influence of concrete hardening
Analyzing the cases considering concrete hardening on cure (Figure 8), it can be observed that when   is high (  = 1,7 • 10 7 and   = 1,7 • 10 6  2 /), consolidation occurs in such a short time interval that hardening does not make a difference, the result remaining equal to those of the elastic supports.For the other cases, it can be observed that the smaller the coefficient of consolidation   , the longer the differential settlement take to occur.Hence, the structure will be more rigid, causing both a reduction in the final differential settlement, and an increase in the internal forces redistribution.Therefore, it can be concluded that the rate of consolidation affects the differential settlements and the internal forces.Furthermore, it is interesting to note that the cases of elastic and rigid supports cease to function as an envelope for possible results.

Influence of cracking with tension-stiffening
The reinforcement adopted for the analysis of the cracked structure was obtained by designing the cross-section with the usual safety factors (e.g., from Brazilian standard NBR 6118:2014).As the beam is subjected to positive and negative bending moments, it was divided into three sections with two different cross-sections as indicated in Table 4.
Table 4 Reinforcement of the beam (Section 1 along the span and Section 2 at the center).

Bending Moment
Reinforcement (cm 2 ) The 6 cases shown in Figure 6 are also presented in Figure 9, this time considering cracking and tension-stiffening.Comparing both figures, it is notable that there is a reduction of approximately 50% in the difference in support reaction and bending moment for cases 1 and 2, leaving the two analyzes with close results.This reduction in the redistribution of forces is due to the increased flexibility of the elements when considering concrete cracking.Analyzing cases 2 and 6, it is observed that both the difference in support reaction and in bending moment were reduced by approximately 70%.As such cases represent the initial and the final behaviour of the structure, and the creep of the concrete defines the transition between them, it can be concluded that in the cracked structure creep has a smaller effect on the redistribution of internal forces.This stems from the fact that, due to cracking, the area of concrete that is mechanically contributing is smaller.
The results of cases 2 and 5 are compared in Figure 10 for the central support of the non-cracked (nc) and cracked (c) beams.Note that, without cracking, the difference in the support reaction varies from 15.8% at the initial time to 10.8% at the final time.Considering cracking, the variation goes from 7.6% to 2.8% over time.The same behaviour is observed for the bending moment.Such reductions indicate that there is less force redistribution in the structure.Furthermore, it is observed that both creep and cracking have a tendency to reduce the redistribution of forces.
Observing the influence of the   value on the results for the central support (Figure 11), it can be noticed that the behaviour of the cracked beam is similar to that of the non-cracked one, as the values tend to remain the same as the results of the elastic support when time goes to infinity, with the exception of rigid support.In this case, such a difference is due to the redistribution of forces between the concrete and the reinforcement caused by creep, which ends up inducing a change in the internal forces.However, the magnitude of the redistribution of forces in the cracked beam is much lower than in the non-cracked one.
Including the concrete hardening effect (Figure 12), it is possible to observe that this phenomenon influences the response of the structure for cases where the   is smaller, similar to the observation made from Figure 8.However, it is noted that the contribution of such an effect is smaller in the cracked structure, where the variation of the support reaction and of the bending moment were 0.2% and 5.9%, respectively.As for the non-cracked structure, the values were 4.1% and 25.4%, in the same order.

CONCLUSIONS
The incorporation of hardening on cure, shrinkage, creep, cracking and tension-stiffening of the concrete, and the consolidation of the soil in the soil-structure interaction allowed an assessment of these phenomena not found in previous publications.The modelling proposed to consider all those phenomena in a coupled analysis also contains some novelties.
The results obtained in the present study enabled the following conclusions: (1) When considering the creep of concrete, the structure tends to become more flexible, reducing the uniformity of settlements and the redistribution of internal forces.The results observed for bending moment and support reactions are found between the results obtained in the analyses with elastic support and with rigid support.
(2) The value of the coefficient of consolidation influences the rate of transition between the behaviour with rigid supports and the behaviour with elastic supports.For high values of the coefficient of consolidation, settlement occurs so quickly that the effects of concrete hardening are not noticeable.However, for smaller values, there is a gain in the structure's stiffness and, consequently, a change in differential settlements and internal forces.
(3) Cracking has a relevant contribution in reducing the stiffness of the structure, significantly reducing the redistribution of forces.Creep effects tend to be attenuated by cracking, as the concrete area that is mechanically stressed is reduced.
(4) Considering both cracking and creep, the reduction in the stiffness of the structure can be considerable, so that settlements and internal forces can be considerably different from those obtained in the more usual elastic soilstructure analysis.
From the above results, it can be said that the consideration of both creep and cracking of the concrete significantly reduces the stiffness of the structure, attenuating the effects of soil-structure interaction.Thus, it is recommended that such phenomena be considered in analyses, for more realistic results.

Figure 1
Figure 1 Loading and beam cross-section (dimensions in mm).

Figure 2
Figure 2 Comparison of experimental results with numerical predictions.

Figure 3
Figure 3 Percentage of settlement as a function of time factor for exact solution (U) and for solutions using Kelvin chains.

Figure 4
Figure 4 (a) Analyzed structure and (b) ground profile with the vertical stress increases.

Figure 6
Figure 6 shows (a) the differential settlement, (b) the reaction force at the central support (normal force in the element representing the foundation), and (c) the bending moment in the beam at the central support.

Figure 6
Figure 6Result for the central support, over time, for the 6 cases analyzed.

Figure 7
Figure 7 Result for the central support, over time, for different values of   .

Figure 8
Figure 8 Result for the central support, over time, for different values of cv, considering concrete hardening.

Figure 9
Figure 9 Results for the central support, over time, for the 5 cases considering tension-stiffening

Figure 10
Figure 10 Comparison of results for the central support of cracked and non-cracked structure.

Figure 11 21 Figure 12
Figure 11 Result for the central support of the cracked structure, over time, for different values of   .
, �; 7. Using the displacement increment in each element Δ ,+1 , determine the displacements in local coordinates without rigid body movement:  +1 =   +   Δ ,+1 ; 8. Using the displacement increment in each element ∆  , calculate the force vector increment: ∆  =  , �∆  −   � ̅ ,+1 � Considering the displacement increment ∆ +1 , calculate the stiffness matrix of the structure  +1 and the resisting nodal force vector  ,+1 ; 14. Calculate the residual of forces  +1 =  +1 −  ,+1 .if the residual norm is less than the established tolerance or the number of iterations is greater than the maximum limit, update the state variables in  + 1 with the values of the variables in  + 1 and proceed to step 15; otherwise, increment in  and return to step 6; 15.Calculate the strain increment Δ  and the total strain  +1 =   + Δ  ; 16.For each viscoelastic material of the cross-section of each element, calculate: ∆  , ∆  ,  ,,+1 ,  ,+1 ,  ,+1 e  +1 ; 17. Increment the value of k and go back to step 1.
9. Start iteration counter  = 0 and use system state at  as initial attempt; 10.Update the values of the nodal forces applied to the elements:   =   + ∆  ; 11.For each force field   , calculate:  ,+1 ,  +1 ,  ,+1 ; 12. Calculate the nodal force increment using line search to determine the value of α that minimizes the norm � +1 −  +1 �: ∆ +1 =  ,+1 � +1 −  +1 � e ∆ ,+1 = ∆ +1 .If the residual strain norm is less than the specified tolerance or the number of iterations is greater than the established maximum limit, use the states in  + 1 for the states in  + 1 and proceed to the next step; otherwise, increment  and return to step 10; 13.

Table 1
Applied loads and average compressive strength of concrete in each test.

Table 2
Material and cross-section characteristics of the beam.

Table 3
Behaviour of the material for the different cases analyzed.