Next Article in Journal
Scientometric Analysis for Mechanical Performance of Broken-Line Long-Span Steel Structure in Construction Considering Geometric Nonlinearity
Next Article in Special Issue
Certain Applications of Generalized Kummer’s Summation Formulas for 2F1
Previous Article in Journal
Connectivity of Semiring Valued Graphs
Previous Article in Special Issue
Equivalent Properties of Two Kinds of Hardy-Type Integral Inequalities
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

The Rheological Analytical Solution and Parameter Inversion of Soft Soil Foundation

College of Water Conservancy and Hydropower Engineering, Hohai University, Nanjing 210098, China
*
Author to whom correspondence should be addressed.
Submission received: 4 June 2021 / Revised: 5 July 2021 / Accepted: 7 July 2021 / Published: 8 July 2021
(This article belongs to the Special Issue Various Approaches for Generalized Integral Transforms)

Abstract

:
In soft soil engineering projects, the building loads are always required to be symmetrically distributed on the surface of the foundation to prevent uneven settlement. Even if the buildings and soft clay are controlled by engineers, it can still lead to the rheology of the foundation. The analytical solution based on the Laplace integral transformation method has positive significance for providing a simple and highly efficient way to solve engineering problems, especially in the long-term uneven settlement deformation prediction of buildings on soft soil foundations. This paper proposes an analytical solution to analyze the deformation of soft soil foundations. The methodology is based on calculus theory, Laplace integral transformation, and viscoelastic theory. It combines an analytical solution with finite theory to solve the construction sequences and loading processes. In addition, an improved quantum genetic algorithm is put forward to inverse the parameters of soft soil foundations. The analytical solution based on Laplace integral transformation is validated through an engineering case. The results clearly illustrate the accuracy of the method.

1. Introduction

“Symmetry” has a precise definition in mathematics. It is usually used to refer to an object that is invariant under some transformations, including translation, reflection, rotation, and scaling. In architecture, whether a “symmetrical layout of blocks, masses, and structures” or “wings and balance of masses” is more artistic has been a popular topic of discussion. For structural engineers, a symmetrical structure is not only artistic, but it also guarantees safety and durability, especially in bad geological conditions. In soft soil engineering projects, high water content, large void ratio, low shear strength, high compressibility, and poor permeability are the typical properties of soft clay [1]. Even if the building is of the required symmetry and the soft clay is reinforced, it can still lead to the rheology of the foundation. A soft soil foundation is mainly composed of soft soil with fine particles of clay and silt, organic soil with large pores, peat, and loose sand. The ground on which the filling and structures stand is unstable and subsided due to high groundwater levels. Long-term monitoring data show that most buildings on soft soil foundations in Chinese coastal areas still maintain a settlement rate of about 1 cm per year; some buildings even reach a settlement rate of 3–4 cm per year. The rheology property results in significant alterations of the stability of a building [2]. The long-term uneven settlement deformation prediction of buildings on soft soil foundations is the primary work that can assure the safety and durability of engineering [3]. In addition, the soft soil parameters are accurately obtained based on engineering, which also faces great challenges.
The roots of research on rheology go well back into the 1930s [4]. During that period, the time dependencies of the stress–strain behavior of materials were calculated following the Bingham theory [5]. The applicability of Bingham’s law to soil was verified by Genie and Zongji in 1948 [6]. After the third ICSMFE conference, the rheology study of soil became popular. The nonlinear rheology deformation characteristics of soil were clear after a lot of experimental research [7]. At present, the rheological models for soft soils mainly include component models, yield surface models, damage rheology model, and empirical models. The component models are comprised of varying series-parallel of Hooke and Kelvin units, such as the Kelvin model, the Maxwell model, and Burger’s model [8]. Su et al. [9] analyzed the rheological behavior of soft soil using the Maxwell model. Wang and Wang [10] presented a semianalytical method to analyze the creep and thermal consolidation behaviors of layered saturated clays. The yield surface is changed in yield surface models when soil creeps. A large number of scholars have studied the change of yield surfaces. These models have features for choosing a correlation or noncorrelation criterion [11], isotropy or anisotropy [12], a single yield surface or a double yield surface [13], and over-consolidated soils or normally consolidated soils [14]. Different damage evolution equations are applied in damage rheology models. Singh adopted a creep rate and a creep function to construct a famous empirical creep model [15]. These rheology studies of soft soil have made great contributions to engineering. Although many studies have been published on the rheological models of soft soil, there is still a lack of analytical solutions.
Deformation monitoring of soft soil foundations has made a great contribution to forward analysis and inversion analysis [16]. Forward analysis is aimed at predicting the future status of engineering by building a regression model between environmental loads and displacement [17,18]. Inversion analysis is devoted to checking stability using the mechanical parameters of the foundation [19]. The latter is more widely accepted because it allows further research on the constitutive models. To achieve parameter inversion, heuristic algorithms are the main methods of optimizing parameters in the feasible region. The quantum genetic algorithm (QGA) was proposed by Narayanan in 1996 [20]. It uses quantum theory in a genetic algorithm (GA) and has excellent searchability for a global optimal solution [21]. The QGA has attracted wide attention and has been developed by scholars. Draa et al. [22] presented a quantum-inspired differential evolution algorithm for solving the N-queens problem. Ma and Jin [23] presented a parallel quantum genetic algorithm. Laboudi and Chikhi [24] made a comparison of GA and QGA. These results show that QGA is a promising tool for exploring search spaces. However, QGA has not been used in the parameter inversion of soft soil foundations.
Laplace integral transformation is an important integral transformation method. It has a good performance in finding analytical solutions of partial differential equations [25]. The research on this method is of concern to mathematicians and engineers. Mathematicians are devoted to modifying the Laplace transform and applying this method to get analytical solutions to classical or modern mathematical problems [26,27]. These modified Laplace transformations always have more useful properties and formulas in theorems [28]. Engineers concentrate on applying the method to engineering projects, such as the rheological property of soil and rock [29], thermal problems [30], thermal viscoelastic problems [31], and so on.
In this paper, we will first show the connection of elastic solutions under normal and tangential distributed force to that under concentrated force. Then, the viscoelastic solutions based on Laplace integral transformation are obtained. The influence of construction sequences and loading is considered through finite element theory. The computational efficiency of the QGA is improved through the features of soft engineering. An engineering case is verified by this method.
Compared with previous research on analytical solutions for soft soil foundations, our contributions in this paper are as follows:
(1)
A three-dimensional analytical viscoelastic deformation solution of soft soil foundations is obtained by us. Some previous studies only considered one-dimensional or two-dimensional conditions. The analytical solution in this paper is simple and has few parameters. The construction sequences and loading processes can also be considered. It shows good adaptability to engineering applications.
(2)
On the basis of our model, an improved quantum genetic algorithm is put forward to inverse the parameters of soft soil foundations. It provides a simple and highly efficient way to predict the long-term deformation of soft soil foundations.

2. The Viscoelastic Solutions to Soft Soil Foundations Based on Laplace Integral Transformation

2.1. Assumptions

The elastic solution is based on elastic mechanics. It has these assumptions:
(1)
It is assumed that the object is continuous, that is, the entire volume of the object is filled by the medium that composes the object, leaving no gaps and maintaining its continuity throughout the deformation process.
(2)
It is assumed that the object is completely elastic. After removing the external force, the object can completely restore its original shape and size. The deformation of the object corresponds to the external force it receives.
(3)
It is assumed that the object is uniform. All parts of the entire object have the same elastic properties.
(4)
It is assumed that the object is isotropic. The elastic properties are the same in all directions, regardless of the direction of investigation.
(5)
It is assumed that the object only has infinitesimal motion. The displacement of the object is much smaller than the size of the object.
The viscoelastic solution is put forward based on the elastic solution. In addition to the above elastic assumptions, the viscoelasticity assumptions are as follows:
(1)
Boltzmann’s superposition principle. The viscoelastic displacement of an object has linear viscoelastic behavior. The creep of the object is a function of the entire loading history. The contribution of the load applied at each stage to the final deformation is independent.
(2)
Elastic–viscoelastic correspondence principle. The linear viscoelastic problem is the transformed linear elastic problem in the Laplace transformed state. The viscoelastic solution can be obtained through the elastic solution and Laplace transformation.

2.2. Elastic Solution under Normal and Tangential Distributed Force

The elastic solution of a semi-infinite space, subject to a normal and tangential concentrated force on the boundary, is as follows [32]:
{ u z z = ( 1 + μ ) P 2 π E R [ 2 ( 1 μ ) + z 2 R 2 ] u z r = ( 1 + μ ) P 2 π E R [ r z R 2 ( 1 2 μ ) r R + z ]
{ u x x = ( 1 + μ ) P 2 π E R { 1 + x 2 R 2 + ( 1 2 μ ) [ R R + z x 2 ( R + z ) 2 ] } u x y = ( 1 + μ ) P 2 π E R [ x y R 2 ( 1 2 μ ) x y ( R + z ) 2 ] u x z = ( 1 + μ ) P 2 π E R [ x y R 2 + ( 1 2 μ ) x R + z ]
where u z z , u z r are the normal and radial deformation; P is the concentrated force; E is the elasticity modulus; μ is Poisson’s ratio; R is the distance from a point to the origin of coordinates, R 2 = r 2 + z 2 = x 2 + y 2 + z 2 . u x x , u x y , u x z are the deformation values in the x, y, z directions under tangential force.
The elastic solutions, subject to concentrated force, can be generalized to the distributed force, which is shown as:
d P = d η d ξ a b
where a , b are the length and width of the distributed force; ξ , η are the length and width of the differential distributed force.
On the surface of the foundation, z is zero. From Equations (1)–(3), the elastic solutions under normal and tangential distributed force are shown as:
{ δ x x = 1 + μ a b π E [ ( 1 μ ) η = x a 2 η = x + a 2 ξ = y b 2 ξ = y + b 2 d ξ d η ξ 2 + η 2 + μ x 2 η = x a 2 η = x + a 2 ξ = y b 2 ξ = y + b 2 d ξ d η ( ξ 2 + η 2 ) 3 2 ] δ x y = ( 1 + μ ) μ a b π E x y η = x a 2 η = x + a 2 ξ = y b 2 ξ = y + b 2 d ξ d η ( ξ 2 + η 2 ) 3 2 δ x z = ( 1 + μ ) ( 1 2 μ ) 2 a b π E x η = x a 2 η = x + a 2 ξ = y b 2 ξ = y + b 2 d ξ d η ( ξ 2 + η 2 ) δ z x = ( 1 + μ ) ( 2 μ 1 ) a b π E η = x a 2 η = x + a 2 ξ = y b 2 ξ = y + b 2 d ξ d η ξ 2 + η 2 x x 2 + y 2 δ z y = ( 1 + μ ) ( 2 μ 1 ) a b π E η = x a 2 η = x + a 2 ξ = y b 2 ξ = y + b 2 d ξ d η ξ 2 + η 2 y x 2 + y 2 δ z z = 1 μ 2 a b π E η = x a 2 η = x + a 2 ξ = y b 2 ξ = y + b 2 d ξ d η ξ 2 + η 2
where δ z x , δ z y , δ z z are the deformation in the x, y and z directions under normal distribution force; δ x x , δ x y , δ x z are the deformation in the x, y and z directions under tangential distribution force.

2.3. Viscoelastic Solutions Based on Laplace Integral Transformation

Viscoelastic models usually consist of one or several Hooke elastomers and Newton viscous bodies. The equations are shown as [33]:
{ k = 0 l p k d k d t k s i j = k = 0 r q k d k d t k e i j k = 0 l 1 p k d k d t k σ i i = k = 0 r 1 q k d k d t k ε i i
where σ = 3 K e , S i j = 2 G e i j , S i j is the deviatoric tensor of stress, e i j is the deviatoric tensor of strain, K is the bulk modulus, G is the shear modulus; p k and q k are soil property parameters. Applying the Laplace transform to Equation (5), one obtains:
{ k = 0 l p k s k S ¯ i j ( s ) = k = 0 r q k s k e ¯ i j ( s ) k = 0 l 1 p k s k σ ¯ i i ( s ) = k = 0 r 1 q k s k ε ¯ i i ( s )
where k = 0 l p k s k = P ¯ ( s ) , k = 0 r q k s k = Q ¯ ( s ) , k = 0 l 1 p k s k = P ¯ ( s ) , k = 0 r 1 q k s k = Q ¯ ( s ) . The time dependence of the elasticity modulus and Poisson’s ratio in the Maxwell model are shown as:
{ E ¯ ( s ) = 9 K · x 6 K + x μ ¯ ( s ) = 3 K x 6 K + x , x = 2 ( G 3 + q 1 s ) ( 1 + p 1 s ) , p 1 = η 1 G 2 , q 1 = ( G 3 G 2 + 1 ) η 1
where G 1   G 2   G 3 are the parameters of shear stiffness; K is the bulk modulus; η 1 is the coefficient of viscosity. From Equations (3) and (7), the viscoelastic deformation based on the Maxwell model and Laplace transforms can be expressed by us, as follows:
{ δ x x = 4 G 3 F i + 3 K F i 2 G 3 F g x 2 + 3 K F g x 2 4 G 3 a b π ( G 3 + 3 K ) 3 η 1 e t ( G 2 G 3 + 3 G 2 K ) G 2 η 1 + G 3 η 1 + 3 K η 1 ( G 2 F i G 2 F g x 2 ) 4 a b π ( G 3 + 3 K ) ( G 2 η 1 + G 3 η 1 + 3 K η 1 ) η 1 e G 2 G 3 t G 2 η 1 + G 3 η 1 ( G 2 F g x 2 + G 2 F i ) 4 G 3 a b π ( G 2 η 1 + G 3 η 1 ) δ x y = 3 G 2 η 1 F g x y e t ( G 2 G 3 + 3 G 2 K ) G 2 η 1 + G 3 η 1 + 3 K η 1 4 a b π ( G 3 + 3 K ) ( G 2 η 1 + G 3 η 1 + 3 K η 1 ) F g x y ( 2 G 3 3 K ) 4 G 3 a b π ( G 3 + 3 K ) G 2 η 1 F g xye G 2 G 3 t G 2 η 1 + G 3 η 1 4 G 3 a b π ( G 2 η 1 + G 3 η 1 ) δ x z = 3 F k x 4 a b π ( G 3 + 3 K ) 3 G 2 η 1 F k x e t ( G 2 G 3 + 3 G 2 K ) G 2 η 1 + G 3 η 1 + 3 K η 1 4 a b π ( G 3 + 3 K ) ( G 2 η 1 + G 3 η 1 + 3 K η 1 ) δ z x = 3 G 2 η 1 F i x e t ( G 2 G 3 + 3 G 2 K ) G 2 η 1 + G 3 η 1 + 3 K η 1 4 a b π x 2 + y 2 ( G 3 + 3 K ) ( G 2 η 1 + G 3 η 1 + 3 K η 1 ) 3 F i x 4 a b π x 2 + y 2 ( G 3 + 3 K ) δ z y = 3 G 2 η 1 F i y e t ( G 2 G 3 + 3 G 2 K ) G 2 η 1 + G 3 η 1 + 3 K η 1 4 a b π x 2 + y 2 ( G 3 + 3 K ) ( G 2 η 1 + G 3 η 1 + 3 K η 1 ) 3 F i y 4 a b π x 2 + y 2 ( G 3 + 3 K ) δ z z = F i ( 4 G 3 + 3 K ) 4 G 3 a b π ( G 3 + 3 K ) 3 G 2 η 1 F i e t ( G 2 G 3 + 3 G 2 K ) G 2 η 1 + G 3 η 1 + 3 K η 1 4 a b π ( G 3 + 3 K ) ( G 2 η 1 + G 3 η 1 + 3 K η 1 ) G 2 η 1 F k e G 2 G 3 t G 2 η 1 + G 3 η 1 4 G 3 a b π ( G 2 η 1 + G 3 η 1 )
where δ x x   δ x y   δ x z are the deformation in the x, y, and z directions under unit distribution force in the x direction; δ z x   δ z y   δ z z are the deformation in the x, y, and z directions under unit distribution force in the z direction; a, b are the length and width of the unit distribution force; x, y are the surface positions of the soft soil foundation; F i , F k ,   F g are the integral formulas. It expresses the influence of distributed force:
{ F i = η = x a 2 η = x + a 2 ξ = y b 2 ξ = y + b 2 d ξ d η ξ 2 + η 2 F k = η = x a 2 η = x + a 2 ξ = y b 2 ξ = y + b 2 d ξ d η ξ 2 + η 2 F g = η = x a 2 η = x + a 2 ξ = y b 2 ξ = y + b 2 d ξ d η ( ξ 2 + η 2 ) 3 / 2

3. The Influence of Construction Sequences and Loading

Equation (8) is the analytical solution under unit distribution force. According to the superposition principle of structural mechanics, the rheology deformation of any point on the surface of a soft soil foundation can be obtained. Figure 1 is the schematic diagram of the deformation of any point on the surface under distribution force. The soft soil foundation acts on the irregularly shaped distribution force. The rheology deformation of any time is calculated using the following steps:
(1)
The irregularly shaped region is dispersed into n little regions. The total time is dispersed into an l time step.
(2)
The irregularly shaped distribution force is dispersed into n distribution forces. The distribution forces in each region are shown as q x 1 ¯ q x n ¯ , q y 1 ¯ q y n ¯ , q z 1 ¯ q z n ¯ , where q x 1 ¯ q x n ¯ are the distribution forces in each region in the x direction; q y 1 ¯ q y n ¯ are the distribution forces in each region in the y direction; q z 1 ¯ q z n ¯ are the distribution forces in each region in the z direction.
(3)
The coordinate systems are set up. It takes each distributed load center point as the origin. The coordinates of the calculated position in each coordinate system are recorded as (x1,y1), (x2,y2),…, (xn,yn).
(4)
According to Equation (8), the deformation of each unit distribution load under each dispersed region is calculated.
(5)
The deformation of the calculated position on the surface under distribution force is shown as:
U x =   q x 1 ¯ δ x x 1 + + q x n ¯ δ x x n   +   q y 1 ¯ δ y x + q y n ¯ δ y x n   +   q z 1 ¯ δ z x + q z n ¯ δ z x n U y =   q x 1 ¯ δ x y 1 + + q x n ¯ δ x y n   +   q y 1 ¯ δ y y + q y n ¯ δ y y n   +   q z 1 ¯ δ z y + q z n ¯ δ z y n U z =   q x 1 ¯ δ x z 1 + + q x n ¯ δ x z n   +   q y 1 ¯ δ y z + q y n ¯ δ y z n   +   q z 1 ¯ δ z z + q z n ¯ δ z z n
where U x , U y , U z are the deformation of the position in the three directions; δ x x 1 , δ x y 1 , δ x z 1 are the deformation of the calculated position in the x, y, and z directions under the first unit distribution force in the x direction; δ y x 1 , δ y y 1 , δ y z 1 are the deformation of the calculated position in the x, y, and z directions under the first unit distribution force in the y direction; δ z x 1 , δ z y 1 , δ z z 1 are the deformation of the calculated position in the x, y, and z directions under the first unit distribution force in the z direction; δ x x n , δ x y n , δ x z n are the deformation of the calculated position in the x, y, and z directions under the nth unit distribution force in the x direction; δ y x n , δ y y n , δ y z n are the deformation of the calculated position in the x, y, and z directions under the nth unit distribution force in the y direction; δ z x n , δ z y n , δ z z n are the deformation of the calculated position in the x, y, and z directions under the nth unit distribution force in the z direction.
Figure 2 shows the flow chart of soft soil foundation surface deformation. The process mainly includes three parts: pre-processing, computational process, and post-processing. In pre-processing, the parameters of the soft soil foundation, the discretization of time, the discretization of forces, and the deformation of each basic coordinate system are calculated. In the computational process, the deformation of each time step is calculated. In post-processing, the results are processed into a graph.

4. The Parameter Inversion of Soft Soil Foundations

4.1. Quantum Genetic Algorithm

Inspired by the theory of quantum mechanics, Narayanan proposed a quantum genetic algorithm (QGA) in 1996. The QGA uses qubit coding instead of traditional chromosome coding. Therefore, it has a parallelism property and advantages in finding the global optimal solution.
The chromosome is shown by quantum states in the QGA method. A quantum state can be represented as the superposition of two states (0, 1). It is shown as:
| Ψ = α | 0 + β | 1
where | Ψ is a quantum state vector; α ,   β are the probability of | 0 and | 1 , α 2 + β 2 = 1 . Then, the individual qubit Γ m t is shown as:
Γ m t = [ τ 1 t τ 2 t   τ m t ] = [ α 1 t α 2 t   α m t β 1 t β 2 t β m t   ]
where τ 1 t ,   τ 2 t ,     τ m t are each quantum states of the individual qubit. The probability of an individual qubit | 0000 001 is ( α 1 t ) 2 + ( α 2 t ) 2 + + ( α m 1 t ) 2 + ( β m t ) 2 . A random number in the interval of [0, 1] can be used to assign binary 0s and 1s by comparing it with the probability of the quantum-bit when the chromosome is converted to a binary format.
The evolution of the population is completed by a quantum revolving gate. It is shown as:
[ α m t β m t ] = [ c o s θ t s i n θ t s i n θ t c o s θ t ] [ α m t β m t ]
where θ t is the parameter of the rotation angle, θ t = S ( α , β ) δ θ ; S ( α , β ) is the direction of rotation; δ θ is the magnitude of the rotation angle.

4.2. Soft Soil Foundation Parameter Inversion Model

The rheology of soft soil foundations is a typical physical phenomenon. The deformation can be obtained by in-situ monitoring. The range of soft soil foundation parameters is shown as:
{ E 1 m i n E 1 E 1 m a x E 2 m i n E 2 E 2 m a x η 1 m i n η 1 η 1 m a x μ m i n μ μ m a x
The initial parameters are M i n i = ( M m i n + M m a x ) / 2 , M = E 1 , E 2 , μ , η 1 . From Equations (8)–(10), the deformation S ( x i , y i , t i ) of the monitoring positions are calculated. The optimal fitness for each generation is shown as:
f = i = 1 n S ( x i , y i , t i ) S ( x i , y i , t i ) 2 2
where S ( x i , y i , t i ) are the in-situ monitoring data. The termination conditions are shown as:
f f k   or   G G m a x
where f k is the termination value of the fitness value; G m a x is the maximum number of generations. If Equation (16) is not true, the quantum revolving gate is used for population evolution, and the fitness is calculated again until the termination conditions are met.

4.3. Improved Quantum Genetic Algorithm

In engineering projects, parameter inversion always needs to be calculated with a large amount of data. The QGA method has low computational efficiency in actual engineering calculations. It needs to be improved. The quantum genetic algorithm improves this in two ways. The algorithm processes of the QGA and IQGA are shown in Figure 3.
(1)
In each generation calculation, some individuals, similar to the optimal individuals of the previous generation, are formed so that individuals with certain guiding effects are produced in the next generation. The number of similar individuals, N, can be decided as shown:
N = [ ( 1 G G m a x ) 10 ] , [ ]   i s   i n t e g r a l   f u n c t i o n
where G is the number of generations.
(2)
The rotation angle is not a constant. The rotation angle is a big value when the number of generations is small. The rotation angle is a small value when the number of generations is big. It can be shown as:
{ δ θ = 0.1 π ,   G G m a x < 0.5 ; δ θ = 0.01 π ,   G G m a x 0.5 ;
Overall, computational efficiency can be accelerated in these theories by generating optimal individuals with high similarity and controlling the mutation process.

5. Case Study

In hydraulic engineering, ship locks are always built on lakes and rivers. These foundations are mainly composed of clay and silt. Affected by the gravity of the structure, the uneven settlement of these soft soil foundations is serious. The monitoring data and construction processes of ship locks are shown in [34]. The ship lock diagram and construction processes are shown in Figure 4. The vertical view and side elevation of a ship lock and the measurement position are shown in Figure 5. The construction processes mainly contain: ① a bottom plate, ② a corridor and an empty box, ③ the first gravel, ④ the second gravel, and ⑤ store water. The length and width of gravel soil in the west are 25 and 30 m. The length and width of gravel soil in the east are 20 and 30 m. The length and width of the ship lock are 40 and 30 m. The bottom plate is finished on the 25th day. The first gravel soil is finished on the 125th day. The corridor and the empty box are finished on the 225th day. The second gravel soil is finished on the 250th day. The store water is finished on the 450th day. The water level is 6 m. The construction processes and simplified loading distribution are shown in Table 1. The density of plate, corridor, and empty box are 2400 kg/m3. The density of gravel soil is 1700 kg/m3. The density of store water is 1000 kg/m3. According to these materials’ densities, the gravity of each part is changed to the area loads.
The range of parameters is shown in Table 2.
The inversion parameters, using the QGA and IQGA methods, are shown as:
QGA : { E 1 = 1.87013 × 10 7   Pa E 2 = 1.59643 × 10 7   Pa η 1 = 5.20775 × 10 13   Pa · s μ = 0.3520   ; I Q G A : { E 1 = 1.82161 × 10 7   Pa E 2 = 1.56227 × 10 7   Pa η 1 = 4.82339 × 10 13   Pa · s μ = 0.3718
The inversion parameters of the two methods are similar and can be used as soft soil foundation parameters. Figure 6 is the convergence of the best fitness value for each generation. The figure shows that the optimal fitness value increases gradually with the increase in generations. The IQGA method converges to the optimal solution in the 20th generation. The QGA method converges to the optimal solution in the 40th generation. IQGA has faster convergence speed than QGA. The calculation efficiency is greatly improved.
Table 3 shows the comparison between the measured values and calculated values of each position. It can be seen from the table that the calculated values at each location are in good agreement with the measured values. The parameters of the Maxwell viscoelastic soft soil foundation obtained by inversion in this paper are reliable. Figure 7 is the prediction displacement of the soft soil foundation in the z direction at each time. The gravity of gravel soil significantly changes the deformation distribution. The maximum deformation position is under the corridor. The creep of soil causes a 0.05–0.08 m settlement of the structure at 51,840,000 s. The displacement of each position in the z direction is shown in Figure 8. As the construction continues, the displacement of the foundation will increase rapidly after each process. At the end of calculation time, the displacement also has a tendency to increase. The measurement displacement can reach −0.075 m.

6. Conclusions

This paper reports a methodology to analyze the rheology of soft soil. Viscoelastic solutions under normal and tangential distributed force were obtained. The influence of construction sequences and loading was considered. An improved quantum genetic algorithm method was used to inverse the parameters of soft soil foundations. An engineering case was used to verify the accuracy of the method. The following important conclusions are summarized:
(1)
The viscoelastic solution of soft soil is based on calculus theory, Laplace integral transformation, and finite theory. The model is simple and has few parameters. The error between the predicted value and the measured value is acceptable.
(2)
The improved quantum genetic algorithm has a faster convergence speed than the quantum genetic algorithm. It is well applied to engineering.
(3)
The case shows that this method can be well appropriated for inversing the parameters of soft soil foundations and predicting the long-term uneven settlement deformation of buildings on soft soil foundations. It can be widely used in soft soil engineering.

Author Contributions

Conceptualization, H.Z. and C.S.; software, H.Z. and J.B.; writing—original draft preparation H.Z. and R.Y.; writing—review and editing, Y.M. and W.W.; All authors have read and agreed to the published version of the manuscript.

Funding

This work was funded by the National Natural Science Foundation of China (no. 51579089).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Arora, K.R. Soil Mechanics and Foundation Engineering, 6th ed.; Lomus Offset Press: Delhi, India, 2006; pp. 587–635. [Google Scholar]
  2. Keedwell, M. Rheology and Soil Mechanics; Elsevier Applied Science: London, UK, 1984; pp. 67–69. [Google Scholar]
  3. Wang, R.; Xu, W.; Meng, Y.; Chen, H.; Zhou, Z. Numerical analysis of long-term stability of left bank abutment high slope at Jinping I hydropower station. Chin. J. Rock Mech. Eng. 2014, 33, 105–3113. [Google Scholar] [CrossRef]
  4. Barden, L. Consolidation of clay with non-linear viscosity. Geotechnique 1965, 15, 345–362. [Google Scholar] [CrossRef]
  5. Bjerrum, L. Engineering geology of Norwegian normally-consolidated marine clays as related to settlements of buildings. Geotechnique 1967, 17, 83–118. [Google Scholar] [CrossRef] [Green Version]
  6. Tjong-Kie, T. Determination of the rheological parameters and the hardening coefficients of Clays. In Rheology and Soil Mechanics/Rhéologie et Mécanique des Sols; Springer: Berlin/Heidelberg, Germany, 1966; pp. 256–272. [Google Scholar]
  7. Berre, T.; Iversen, K. Oedometer test with different specimen heights on a clay exhibiting large secondary compression. Geotechnique 1972, 22, 53–70. [Google Scholar] [CrossRef]
  8. Perzyna, P. Fundamental problems in viscoplasticity. Adv. Appl. Mech. 1966, 9, 243–377. [Google Scholar] [CrossRef]
  9. Su, C.; Zhang, H.; Hu, S.; Bai, J.; Dai, J. Combining Finite Element and Analytical methods to Contact Problems of 3D Structure on Soft Foundation. Math. Probl. Eng. 2020, 2020, 8827681. [Google Scholar] [CrossRef]
  10. Wang, L.; Wang, L. Semianalytical analysis of creep and thermal consolidation behaviors in layered saturated clays. Int. J. Geomech. 2020, 20, 06020001. [Google Scholar] [CrossRef]
  11. Borja, R.I. Generalized creep and stress relaxation model for clays. J. Geotech. Eng. 1992, 118, 1765–1786. [Google Scholar] [CrossRef]
  12. Yin, Z.Y.; Chang, C.S.; Karstunen, M.; Hicher, P.-Y. An anisotropic elastic–viscoplastic model for soft clays. Int. J. Solids Struct. 2010, 47, 665–677. [Google Scholar] [CrossRef] [Green Version]
  13. Alshamrani, M.A.; Sture, S. A time-dependent bounding surface model for anisotropic cohesive soils. Soils Found. 1998, 38, 61–76. [Google Scholar] [CrossRef] [Green Version]
  14. Adachi, T.; Oka, F.; Mimura, M. Mathematical structure of an overstress elasto-viscoplastic model for clay. Soils Found. 1987, 27, 31–42. [Google Scholar] [CrossRef] [Green Version]
  15. Singh, A.; Mitchell, J.K. General stress-strain-time function for soils. J. Soil Mech. Found. Div. 1968, 94, 21–46. [Google Scholar] [CrossRef]
  16. Ji, W.; Liu, X.; Qi, H.; Liu, X.; Lin, C.; Li, T. Mechanical Parameter Identification of Hydraulic Engineering with the Improved Deep Q-Network Algorithm. Math. Probl. Eng. 2020, 2020, 6404819. [Google Scholar] [CrossRef]
  17. Lin, C.; Li, T.; Chen, S.; Lin, C.; Liu, X.; Gao, L.; Sheng, T. Structural identification in long-term deformation characteristic of dam foundation using meta-heuristic optimization techniques. Adv. Eng. Softw. 2020, 148, 102870. [Google Scholar] [CrossRef]
  18. Chen, S.; Gu, C.; Lin, C.; Zhao, E.; Song, J. Safety monitoring model of a super-high concrete dam by using RBF neural network coupled with kernel principal component analysis. Math. Probl. Eng. 2018, 2018, 1712653. [Google Scholar] [CrossRef] [Green Version]
  19. Lin, C.; Li, T.; Chen, S.; Liu, X.; Lin, C.; Liang, S. Gaussian process regression-based forecasting model of dam deformation. Neural Comput. Appl. 2019, 31, 8503–8518. [Google Scholar] [CrossRef]
  20. Narayanan, A.; Moore, M. Quantum-inspired genetic algorithms. In Proceedings of the IEEE International Conference on Evolutionary Computation, Nagoya, Japan, 20–22 May 1996; pp. 61–66. [Google Scholar] [CrossRef]
  21. Grover, L.K. A fast quantum mechanical algorithm for database search. In Proceedings of the Twenty-Eighth Annual ACM Symposium on Theory of Computing, Philadelphia, PA, USA, 22–24 May 1996; Association for Computing Machinery: New York, NY, USA, 1996. [Google Scholar]
  22. Draa, A.; Meshoul, S.; Talbi, H.; Batouche, M. A quantum-inspired differential evolution algorithm for solving the N-queens problem. Neural Netw. 2011, 1, 21–27. [Google Scholar] [CrossRef] [Green Version]
  23. Ma, S.; Jin, W. A new parallel quantum genetic algorithm with probability-gate and its probability analysis. In Proceedings of the 2007 International Conference on Intelligent Systems and Knowledge Engineering, Chengdu, China, 15–16 October 2007; Atlantis Press: Paris, France, 2007. [Google Scholar] [CrossRef] [Green Version]
  24. Laboudi, Z.; Chikhi, S. Comparison of genetic algorithm and quantum genetic algorithm. Int. Arab J. Inf. Technol. 2012, 9, 243–249. [Google Scholar]
  25. Zhiber, A.V.; Startsev, S.Y. Integrals, solutions, and existence problems for Laplace transformations of linear hyperbolic systems. Math. Notes 2003, 74, 803–811. [Google Scholar] [CrossRef]
  26. Yüce, A.; Tan, N. On the approximate inverse Laplace transform of the transfer function with a single fractional order. Trans. Inst. Meas. Control 2021, 43, 1376–1384. [Google Scholar] [CrossRef]
  27. Ender, A.Y.; Énder, I.A.; Bakaleinikov, L.A.; Flegontova, E.Y. Construction of kernels of the nonlinear collision integral in the Boltzmann equation using laplace transformation. Tech. Phys. 2012, 57, 735–742. [Google Scholar] [CrossRef]
  28. Kim, Y.; Kim, B.M.; Jang, L.-C.; Kwon, J. A Note on Modified Degenerate Gamma and Laplace Transformation. Symmetry 2018, 10, 471. [Google Scholar] [CrossRef] [Green Version]
  29. Zhou, F.; Wang, L.; Liu, H. A fractional elasto-viscoplastic model for describing creep behavior of soft soil. Acta Geotech. 2021, 16, 67–76. [Google Scholar] [CrossRef]
  30. Varghese, V.; Bhoyar, S.; Khalsa, L. Thermoelastic response of a nonhomogeneous elliptic plate in the framework of fractional order theory. Arch. Appl. Mech. 2021, 91, 3223–3246. [Google Scholar] [CrossRef]
  31. Magdy, A.E. Analytical study of two-dimensional thermo-mechanical responses of viscoelastic skin tissue with temperature-dependent thermal conductivity and rheological properties. Mech. Based Des. Struct. Mach. 2021, 1–19. [Google Scholar] [CrossRef]
  32. Davis, R.; Selvadurai, A.; Pacheo, M. Elasticity and Geomechanics; Cambridge University Press: Cambridge, UK, 1998; p. 157. [Google Scholar]
  33. Chen, E.Y.; Pan, E.; Green, R. Surface loading of a multilayered viscoelastic pavement: Semianalytical solution. J. Eng. Mech. 2009, 135, 517–528. [Google Scholar] [CrossRef] [Green Version]
  34. Su, C.; Jiang, H.D.; Tan, E.H. An inverse analysis method for foundation parameters and its application based on viscoelastic foundation beam computation. Chin. J. Geotech. Eng. 2000, 22, 186–189. [Google Scholar]
Figure 1. The schematic diagram of the deformation of any point on the surface under distribution force.
Figure 1. The schematic diagram of the deformation of any point on the surface under distribution force.
Symmetry 13 01228 g001
Figure 2. Flow chart of soft soil foundation surface deformation.
Figure 2. Flow chart of soft soil foundation surface deformation.
Symmetry 13 01228 g002
Figure 3. The algorithm processes of QGA and IQGA.
Figure 3. The algorithm processes of QGA and IQGA.
Symmetry 13 01228 g003
Figure 4. The ship lock diagram and construction processes.
Figure 4. The ship lock diagram and construction processes.
Symmetry 13 01228 g004
Figure 5. The vertical view and side elevation of the ship lock and the measurement position (descriptions of ①, ②, ③, ④, and ⑤ are shown in Figure 4).
Figure 5. The vertical view and side elevation of the ship lock and the measurement position (descriptions of ①, ②, ③, ④, and ⑤ are shown in Figure 4).
Symmetry 13 01228 g005
Figure 6. Convergence of the optimal fitness for each generation.
Figure 6. Convergence of the optimal fitness for each generation.
Symmetry 13 01228 g006
Figure 7. The prediction displacement of the soft soil foundation in the z direction at each time.
Figure 7. The prediction displacement of the soft soil foundation in the z direction at each time.
Symmetry 13 01228 g007
Figure 8. The displacement of each position in the z direction (descriptions of ①, ②, ③, ④ and ⑤ are shown in Figure 4).
Figure 8. The displacement of each position in the z direction (descriptions of ①, ②, ③, ④ and ⑤ are shown in Figure 4).
Symmetry 13 01228 g008
Table 1. The construction processes and simplified loading distribution.
Table 1. The construction processes and simplified loading distribution.
Construction ProcessesTimeLoads
Backfill in West (N/m2)Bottom Plate in West (N/m2)Middle Plate (N/m2)Bottom Plate in East (N/m2)Backfill in East (N/m2)
① Bottom plate2,160,000 s072,00072,00072,0000
② Corridor and empty-box10,800,000 s0236,00072,000200,0000
③ First gravel19,440,000 s68,000236,00072,000200,00068,000
④ Second gravel21,600,000 s102,000236,00072,000200,000102,000
⑤ Store water38,880,000 s102,000236,000132,000200,000102,000
Table 2. The range of parameters.
Table 2. The range of parameters.
ParametersRange
E 1 1.2 × 10 7 P a 2 × 10 7   Pa
E 2 1.2 × 10 7 P a 2 × 10 7   Pa
η 1 5.18 × 10 13 P a · s 7.78 × 10 13   Pa · s
μ 0.15 ~ 0.40
Population size40
Maximum generation50
Table 3. The comparison of calculated values and measured values (unit: m).
Table 3. The comparison of calculated values and measured values (unit: m).
PositionDataTime
2,160,000 s10,800,000 s19,440,000 s21,600,000 s38,880,000 s51,840,000 s
1Measured−0.0192−0.0399−0.0542−0.0548−0.0609−0.0721
QGA−0.0171−0.0381−0.0527−0.0541−0.0607−0.0744
IQGA−0.0169−0.0378−0.0526−0.0541−0.0608−0.0744
2Measured−0.0191−0.0378−0.0619−0.0631−0.0793−0.0850
QGA−0.0176−0.0333−0.0461−0.0473−0.0529−0.0694
IQGA−0.0173−0.0331−0.0461−0.0473−0.0530−0.0694
3Measured−0.0173−0.0302−0.0404−0.0400−0.0467−0.0599
QGA−0.0161−0.0345−0.0479−0.0492−0.0550−0.0698
IQGA−0.0160−0.0343−0.0479−0.0492−0.0550−0.0699
4Measured−0.0177−0.0314−0.0430−0.0435−0.0589−0.0638
QGA−0.0145−0.0387−0.0546−0.0562−0.0633−0.0751
IQGA−0.0145−0.0386−0.0547−0.0563−0.0635−0.0754
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Zhang, H.; Su, C.; Bai, J.; Yuan, R.; Ma, Y.; Wang, W. The Rheological Analytical Solution and Parameter Inversion of Soft Soil Foundation. Symmetry 2021, 13, 1228. https://0-doi-org.brum.beds.ac.uk/10.3390/sym13071228

AMA Style

Zhang H, Su C, Bai J, Yuan R, Ma Y, Wang W. The Rheological Analytical Solution and Parameter Inversion of Soft Soil Foundation. Symmetry. 2021; 13(7):1228. https://0-doi-org.brum.beds.ac.uk/10.3390/sym13071228

Chicago/Turabian Style

Zhang, Heng, Chao Su, Jiawei Bai, Rongyao Yuan, Yujun Ma, and Wenjun Wang. 2021. "The Rheological Analytical Solution and Parameter Inversion of Soft Soil Foundation" Symmetry 13, no. 7: 1228. https://0-doi-org.brum.beds.ac.uk/10.3390/sym13071228

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop