There is still a growing increase in interest in the use of mathematical models in geotechnical design. Most commonly, these models are developed based on finite element method (FEM).
A successful simulation requires a careful model construction. A very important part of constructing the model is the calibration of model parameters. The aim of this paper is to present a methodology of automatic soil model calibration employing inverse analysis, namely genetic algorithms (GA). Inverse problem requires that mathematical model is closely related to the search algorithm. Building blocks of the calibration system are briefly described in this paper.
The objective methods of parameter calibration are usually based on inverse analysis. The general assumption here is that the estimation should lead to the maximum possible imitation of the experiment by the calibrated model. With the use of inverse methods, the obtained parameters values would, in general, have rather ‘effective’ than ‘physical’ meaning. This is because of the maximal correspondence between the experiment and the simulation result would cause, as the not modelled phenomena would also manifest in the obtained parameters values. This could be a drawback for theoretical analysis, but for practical (engineering) modelling this is actually an advantage, as it allows us to include in the results aspects that are not formally included in the model.
Inverse calibration techniques are based on some search algorithms. Typical implementation would search through the calibrated parameters space to find a set of parameters leading to maximal concurrency with experimental data. The correspondence between both must be expressed as a search criteria and the most common is the least square error between experimental data and simulation results. The typical estimation procedure is presented in Figure 1.
Our experiments are based on a model of multiphase medium, assuming that the solid has hydraulically connected pores or microcracks, which allow for the flow of liquid and gas filtration. A mathematical model of the creep of a porous medium defined as a multicomponent diphase body was first introduced by Biot [13,14]. In a very general way the model has been analyzed on the basis of asymptotic homogenization theory for the case when the diphase medium’s pores are filled with a poorly compressible liquid. For this purpose [3,5,12,21] used the methods for periodic structures while [20,23] used statistical methods. For the case when the medium’s pores are filled with a gas, a mathematical model based on the asymptotic homogenization method was presented by [4]. One of the fundamental assumptions of porous medium mechanics is the porous medium continuity postulate. It is a considerable simplification when this theory is applied to soils. Despite this assumption, models describing the creep of the Biot body represent well the actual processes of deformation of cohesive soils (argillaceous and clayey soils and mudclays). The processes are described also taking into account temperature changes by [26]. In [9] PhD thesis and in [8,10] one can find consolidation equations based on the Biot model with the Kelvin-Voight rheological skeleton, showing that besides volumetric compressibility and shape elasticity, the viscosity of the soil skeleton also plays a significant role.
The model calibration was performed employing inverse analysis using genetic algorithms (GA). The objective of the analysis is the estimation of filtration coefficient
The research was performed on clays obtained from the mine dump soil KWB Turow. An alternative way was presented to determine the coefficient of filtration than the previously commonly used method based on the Terzahgi model. That method is not complete because it does not take into account the soil skeleton deformation. The work appears to confirm the thesis that genetic algorithms are a highly effective tool for automatic analysis, based on simple rules and simultaneous calibration of all the parameters. It shows how to determine the basic parameters of the models using the example of the Biot model. The presented paper is a continuation and development of research concerning the determination of the parameters of complex rheological models [9,11].
The classical constitutive relations for a poroelastic body in isothermal problems were given by [14]. The effective parameters of the Biot model (Biot’s constants) were described for the first time in 1957 [15] as
where
Consolidation equations were derived by [14] and have the following form:
where:
In his postdoctoral thesis [27] put forward an idea of using the genetic computation technique to solve the selected geotechnical engineering problems (the stability of scarps and slopes, the interpretation of the results of triaxial compression tests, etc.). He demonstrated that the state-of-the-art numerical techniques can restore the full functionality of some solutions that have existed in geotechnics for decades, but which have not been applied in practice because of, for example, too high computing costs. Biot’s model calibration based on inverse analysis, carried out by us using the genetic algorithm, has proved the belief that genetic algorithms constitute a highly effective calibration tool for the automatic (based on simple rules and simultaneous for all the parameters) determination of effective parameters, also for the more complex problems of mechanics. The results presented here can be used as the basis for further research on the application of genetic algorithms to the analysis of the effective parameters of the Biot model with the viscosity of the soil skeleton taken into account.
Genetic algorithms (GAs) can be defined as a numerical technique inspired by natural processes of information processing [17]. They are classified as resistant optimization methods, or methods showing similar efficacy for a wide range of issues [19]. This flexibility derives from the fact that the optimization is performed based only on the value of the objective function. The algorithm does not have any more information of the model or of a modeled process [27]. This is a very important property of GA, as it means that for different models geometries, processes, parameters sets and experiments, they would have similar effectiveness and can be used without modification.
The methodology is, in fact, the inverse estimation. During the estimation process, the algorithm generates sets of parameters, which are used for process simulation. Based on results obtained, the estimator least-squared error is applied for valuation of the obtained solutions (Figure 1). Based on error values, more sets are generated and the process is continued until satisfying results are obtained.
GAs have proved to be very successful in many scientific and engineering areas [2,25]. A great theoretical introduction for the calibration with GAs was given by [18]. The complicated model calibration was used by [1,22,28] and many others.
The terminology connected with genetic algorithms includes, beside typical optimization terms, also terms derived from genetics. This is linked with the way genetic algorithms work, which, as mentioned above, imitates natural evolutionary processes.
GAs operate on variables represented by code combinations called
The data structures (code combinations) on which GA operates are converted into values of the each variable. Then the variables are inserted into the calibrated model and used for simulation and then evaluation by means of the objective function calculation (least square error estimator). The set of parameter values and the evaluation constitute a
Figure 2 shows a block diagram illustrating how genetic algorithms work. Seven principal steps can be distinguished in the diagram [24].
The first step in the operation of GAs is the selection of an initial population. In our case, the initialization was done through the random selection of the values of consecutive genes.
The adaptation of the chromosomes in the population is evaluated on the basis of the adaptation function values of all the individuals in the current generation. This process takes place in several stages. First, using an appropriate decoding method, the chromosomes are converted into values of variables. Then the variables are used to calculate the objective function for all the individuals in a given population.
The stop condition is defined on the basis of the optimization requirements. In the simplest case, the performance of a set number of iterations, the testing of a proper number of points of space, or time out can be the condition. In our case, the algorithm was stopped when there was no significant improvement in the results.
The selection of chromosomes is conducted in order to isolate chromosomes for a parental pool, from which a new population (a set of points) is created using genetic operators. We used the roulette method with the scaling function. Regardless of the method used, selection is always made solely on the basis of the adaptation function values for the individuals.
The chromosomes of the individuals selected for a parental pool are subjected to the operation of genetic operators, as a result of which a population of descendants is obtained (Figure 3). We used two classic genetic operators, i.e. the crossover operator and the mutation operator. The crossover operator operates simultaneously on two chromosomes selected from a parental pool, causing the exchange of fragments of their chromosomes. Mutation is performed on one gene, with an assigned probability, causing a change in its state.
As a result of the operation of the genetic operators, a new population is formed from the individuals. This population becomes a current population and objective function and adaptation function values are calculated for it. Then the stop condition is checked and if the computations need to be continued, another cycle is started. When the stop condition is reached, the result is given as the output. The result is a set of values of the decision variables, calculated for the chromosome with the highest (when a maximum is sought) value of the objective function.
Turoszow clays, obtained from the mine dump soil KWB Turow, due to its structural complexity, were used in the presented research. Because of its structure, this material is extremely interesting for research, and problematic from the perspective of engineering practice. The presence of water associated with the boundary surface of clay particles causes that the mechanical processes occurring in these soils are much more complex and should take into account the influence of clay particles’ contact through bound water. Therefore, defining individual parameters with the greatest accuracy is a key task in the context of the modeling of engineering processes of this soil.
Basic soil parameters are:
liquid limit: plasticity limit: density: density of solid particles initial sample soil moisture: porosity
The consolidation chamber used in the studies allowed the application of vertical load
The presented consolidation chamber was used for the oedometer test. The pressures applied were: 12.5, 25, 50 and 100 kPa. Initial sample height
Based on the results, the bulk modulus was calculated:
Oedometer test results were used as calibration data in the investigations. The calibration itself was carried out using a program implementing the classical genetic algorithm, designed by the authors, and the FlexPDE software for running FEM simulation needed for an inverse analysis. Thanks to this software combination, calibration can be performed exactly on the model, which is later used to solve actual engineering problems. In addition, owing to the FlexPDE software capabilities a complex experiment design involving the gradual loading of an oedometer sample can be used. A mathematical model was implemented in the FlexPDE program.
The algorithm was implemented in the form of an independent program written in the C# language in the .NET environment. The program was written to make it possible to calibrate any mathematical model with the help of an external program supplying simulation results.
The scheme of C# software operation in connection with the FlexPDE program is presented in Fig. 8. At the beginning of its operation, the GA reads laboratory measurement results from a file. During each cycle, GA generates parameter values and inserts it to basic pde file. Then it starts FlexPDE and instructs it to run simulation from this file and save its results. Next, GA reads simulation results and compares it with the data read earlier. The comparison is based on least square error estimator (objective function). The operation is stopped when completion criteria is met.
The calibration was performed with the use of a specially developed software implementing genetic algorithm and FlexPDE. The experimental data from four-stage oedometer test (as described in chapter 6) were used for the calibration. The mathematical model was a third Biot model with a rheological Kelvin–Voight skeleton. The subject of calibration was only one parameter . To check the reproducibility of the method, 20 independent estimations were performed.
The results of the calibration are presented in Table 1. The values of hydraulic conductivities obtained varies, but the relatively small standard deviation
Parameter estimation results.Exp. no Exp. no 1 70.52 11 62.43 2 63.06 12 75.83 3 61.44 13 61.23 4 66.63 14 61.12 5 62.43 15 65.01 6 71.05 16 60.39 7 60.70 17 61.85 8 65.01 18 62.43 9 66.58 19 62.43 10 61.85 20 72.94
In the present paper, model calibration was carried out employing an inverse analysis using the still novel method of genetic algorithms for determining the hydraulic conductivity
It has been demonstrated that genetic algorithms are a highly effective tool enabling automatic calibration based on simple rules. Experiments with genetic algorithms prove its utility for model parameter estimation. When least-squared estimator of simulation quality is used, values of hydraulic conductivity are of at least one order smaller than values obtained with the Terzaghi model. This is in the agreement with common observations that this method often leads to overestimations of
The obtained results will form the basis for further research as part of which still other rheological models and sets of experimental data can be analyzed.