Mathematical Problems in Engineering

Mathematical Problems in Engineering / 2014 / Article

Research Article | Open Access

Volume 2014 |Article ID 482672 |

Hongbo Zhao, Zenghui Huang, Zhengsheng Zou, "Simulating the Stress-Strain Relationship of Geomaterials by Support Vector Machine", Mathematical Problems in Engineering, vol. 2014, Article ID 482672, 7 pages, 2014.

Simulating the Stress-Strain Relationship of Geomaterials by Support Vector Machine

Academic Editor: Delfim Soares Junior
Received19 Apr 2014
Revised13 Jul 2014
Accepted18 Jul 2014
Published14 Aug 2014


Stress-strain relationship of geomaterials is important to numerical analysis in geotechnical engineering. It is difficult to be represented by conventional constitutive model accurately. Artificial neural network (ANN) has been proposed as a more effective approach to represent this complex and nonlinear relationship, but ANN itself still has some limitations that restrict the applicability of the method. In this paper, an alternative method, support vector machine (SVM), is proposed to simulate this type of complex constitutive relationship. The SVM model can overcome the limitations of ANN model while still processing the advantages over the traditional model. The application examples show that it is an effective and accurate modeling approach for stress-strain relationship representation for geomaterials.

1. Introduction

Laboratory testing is the primary tool used by engineers and researchers for understanding geomaterials behavior. In the past few years, the use of artificial neural networks (ANN) has been introduced as an alternative approach to stress-strain relationship of geomaterials, and some of these ANN based constitutive relationships have successfully been applied to numerical analysis with improved accuracy [112]. However, ANN based model still has the following limitations [13].(i)ANN does not provide information about the relative importance of the various parameters.(ii)The knowledge acquired during the training of the ANN model is stored in an implicit manner and hence it is very difficult to have a reasonable interpretation of the overall structure of the network.(iii)In addition, ANN model has some drawbacks such as slow convergence speed, less generalizing performance, arriving at local minimum, and overfitting problems.Support vector machine seems to be a promising technique to circumvent these limitations. In recent years, support vector machine (SVM) methods have been rapidly developed for universal function approximations [14]. In geotechnical engineering, SVM has been applied to modeling nonlinear displacement time series of the high slope of the permanent Shiplock of the Three Gorges Project and a large landslide in China [15]. More applications of SVM can be seen in [1618].

Therefore, in this paper, SVM will be studied as an alternative tool for ANN to simulate the stress-strain relationship for geomaterials. In Section 2, SVM approach is introduced; in Section 3, application of SVM approach in simulating the stress-strain relationship of geomaterials is given; conclusions are made in Section 4.

2. Support Vector Machine (SVM)

The SVM was firstly proposed by Vapnik and is illustrated in Figure 1 [19]. It is used to train nonlinear relationships based on the structural risk minimization principle that seeks to minimize an upper bound of the generalization error rather than to minimize the empirical error implemented in neural networks. Merit of the SVM is that training it is a uniquely solvable quadratic optimization problem. The SVM uses nonlinear mapping based on an internal integral function to transform an input space to a high dimension space and then looks for a nonlinear relationship between inputs and outputs in that space. The SVM not only has theoretical support but also can find global optimum solutions for problems with small training samples, high dimensions, and nonlinear and local optima. A wide variety of applications such as pattern recognition and nonlinear regression have empirically shown the SVM’s ability of generalization.

Suppose that we are given a set of observation data (samples)  , ,  . For the regression problem based on the SVM, we can get the following regression function: where is the kernel function; that is, . , , and are obtained by solving the following quadratic programming problem:

maximize subject to The constant determines the tradeoff between the flatness of and the amount up to which deviations larger than are tolerated.

3. The Stress-Strain Model Based on SVM

In order to have a SVM-based stress-strain relationship for a geomaterial, firstly, laboratory or numerical test of the behavior of rocks and clay under different loading patterns is conducted, and then the obtained data are used to train a SVM model. If the training data contains enough relevant information, the trained SVM is supposed to be able to generalize the stress-strain relationship and predict the behavior of materials under new loading circumstances.

SVM modeling involves determination of input variables, output variables, and parameters of SVM (such as kernel function and constant mentioned above). An appropriate selection of the input and output is also important for a successful simulation using SVM. In strain-controlled triaxial compression tests, axial strains are given as input, and the corresponding stress responses are considered output. The simplest SVM-based stress-strain relationship of geomaterial is shown in Figure 2 as an example. Based on support vector machine, the stress-strain relationship can be expressed as The stress-strain relationship under stress-controlled circumstance can be written in the following form alternatively: where and denote values at different load time steps.

3.1. Stress-Strain Relationship of Clay Based on SVM

In order to examine the ability of SVM model, a numerical test of consolidated-drained (CD) triaxial compress of normally consolidated clays was calculated using modified Cam clay model. The test was used to calculate the consolidated-drained triaxial behavior of a normally consolidated clay specimen subjected to a confining pressure 137.8 kPa [20]. The soil parameters are bulk modulus , Lame constant , and initial void ratio . It is well known that stress-strain relationship of geomaterials is greatly influenced by such important factors as stress path. Therefore past history of stress and strain is part of the input data.

Finally, 19 samples were generated based on modified Cam model (Table 1). Figure 3 shows the excellent generalization performance of SVM.

(kPa) (kPa)


3.2. Stress-Strain Relationship of Rock Based on SVM

SVM representation of the stress-strain relationship for rock is investigated in this section.

To provide the input data for training and validating a SVM model, experimental tests were conducted in laboratory. In laboratory, rock specimens (taken from a Liangbei coal mine site in Henan, China) were subjected to triaxial compression tests.

Triaxial compression tests under various confined pressures were carried out (see Table 2). For the purpose of demonstration, stress-strain relationships are shown in Figure 4. This indicates high nonlinearity between axial stress and strain.

Test IDDiameter (mm)Height (mm)Confined pressure (Mpa)


The test data sets B-2-5 were used as the training data to build the support vector machine model. There are 558 test records in specimen B-2-5 and compose 556 sample for support vector machine. The sample was randomly separated into two groups, that is, training samples and testing samples. There are 277 training samples and 279 testing samples. When the stress-strain relationship was built based on test data form specimen B-2-5, the model can predict the stress-strain relationship of other specimens.

SVM model is obtained by being trained with the stress and strain data generated from triaxial compression test, and this SVM model is supposed to represent the relationship between stresses and strains. Here, the SVM model was tested with both trained and untrained data to examine its performance in generalization and prediction. Figures 5 and 6 show the stress-strain relationship measured in the experiments and generated by the SVM model for trained and untrained samples. It can be observed from the figure that SVM model is capable of describing the complex, nonlinear stress-strain relationship of geomaterials such as rock. Figures 7, 8, 9, and 10 demonstrate further the prediction capability of the SVM model under various confining pressures by comparisons with experimental results. It can be seen that the obtained SVM model can represent the stress-strain relationship of geomaterials. To the special phenomenon such as soften feature, it needs more data (information) to build the relationship.

From the preceding procedure, it is evident that with the SVM modeling one does not need to find a series of material parameters describing mathematical equations associated with a constitutive model. This is one of the important advantages of the SVM model over a traditional constitutive model, since the computation of material parameters is usually a very tedious and difficult process and therefore prone to error.

4. Conclusions and Discussions

Support vector machine provides an effective alternative for modeling mechanical behavior of geomaterials by overcoming some drawbacks of ANN model. The excellent performance of the SVM model is demonstrated by successful simulation and prediction of stress-strain relationship of geomaterials under various confining pressures. Meanwhile, similar to ANN model, the SVM model has some advantages over traditional approaches. Firstly, the SVM model is essentially based on experimental data only. No assumptions are made, which allows the model to become more objective. In other words, the SVM model is not to be influenced by the shape of stress-strain curves. This feature is of particular significance in dealing with geomaterials constitutive behavior. Secondly, the SVM model is set up without any calculation of parameters required by a mathematical constitutive model. Therefore, the SVM model is simple and effective for stress-strain relationship modeling, if appropriate experimental data are available for the geomaterial. In this paper, it only proposes a general frame of stress-strain using SVM which can be used like the traditional stress-strain relationship. To the application of model in numerical analysis, it will be considered in later research.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.


The authors wish to acknowledge the anonymous reviewer for his/her helpful comments on the paper. The financial support from the Program for New Century Excellent Talents in University (China) (no. NCET-08-0662), National Science Foundation of China (41072224), and Innovation Group of Henan Polytechnic University (no. T2014-2) is greatly acknowledged.


  1. J. Ghaboussi, J. H. Garrett, and X. Wu, “Knowledge-based modeling of material behavior with neural networks,” Journal of Engineering Mechanics, vol. 117, no. 1, pp. 132–153, 1991. View at: Publisher Site | Google Scholar
  2. D. E. Sidarta and J. Ghaboussi, “Coonstitutive modeling of geomaterials from non-uniform material tests,” Computers and Geotechnics, vol. 22, no. 1, pp. 53–71, 1998. View at: Publisher Site | Google Scholar
  3. J. Ghaboussi and D. E. Sidarta, “New nested adaptive neural networks (NANN) for constitutive modeling,” Computers and Geotechnics, vol. 22, no. 1, pp. 29–52, 1998. View at: Publisher Site | Google Scholar
  4. S. Jung and J. Ghaboussi, “Neural network constitutive model for rate-dependent materials,” Computers and Structures, vol. 84, no. 15-16, pp. 955–963, 2006. View at: Publisher Site | Google Scholar
  5. Y. M. A. Hashash, S. Jung, and J. Ghaboussi, “Numerical implementation of a neural network based material model in finite element analysis,” International Journal for Numerical Methods in Engineering, vol. 59, no. 7, pp. 989–1005, 2004. View at: Publisher Site | Google Scholar | Zentralblatt MATH
  6. Y. M. Najjar and C. Huang, “Simulating the stress-strain behavior of Georgia kaolin via recurrent neuronet approach,” Computers and Geotechnics, vol. 34, no. 5, pp. 346–361, 2007. View at: Publisher Site | Google Scholar
  7. A. Garaga and G. M. Latha, “Intelligent prediction of the stress-strain response of intact and jointed rocks,” Computers and Geotechnics, vol. 37, no. 5, pp. 629–637, 2010. View at: Publisher Site | Google Scholar
  8. D. Penumadu and R. Zhao, “Triaxial compression behavior of sand and gravel using artificial neural networks (ANN),” Computers and Geotechnics, vol. 24, no. 3, pp. 207–230, 1999. View at: Publisher Site | Google Scholar
  9. M. Banimahd, S. S. Yasrobi, and P. K. Woodward, “Artificial neural network for stress-strain behavior of sandy soils: knowledge based verification,” Computers and Geotechnics, vol. 32, no. 5, pp. 377–386, 2005. View at: Publisher Site | Google Scholar
  10. S. Pernot and C. H. Lamarque, “Application of neural networks to the modelling of some constitutive laws,” Neural Networks, vol. 12, no. 2, pp. 371–392, 1999. View at: Publisher Site | Google Scholar
  11. J. H. Zhu, M. M. Zaman, and S. A. Anderson, “Modeling of shearing behavior of residual soil with recurrent neural network,” International Journal for Numerical and Analysis Method in Geomechnics, vol. 22, pp. 671–687, 1998. View at: Google Scholar
  12. J. H. Zhu, M. M. Zaman, and S. A. Anderson, “Modeling of soil behavior with a recurrent neural network,” Canadian Geotechnical Journal, vol. 35, no. 5, pp. 858–872, 1998. View at: Publisher Site | Google Scholar
  13. D. Park and L. R. Rilett, “Forecasting freeway link travel times with a multilayer feedforward neural network,” Computer-Aided Civil and Infrastructure Engineering, vol. 14, no. 5, pp. 357–367, 1999. View at: Publisher Site | Google Scholar
  14. V. N. Vapnik, S. E. Golowich, and A. J. Smola, “Support vector method for function approximation, regression estimation, and signal processing,” Advanced Neural Information Processing Systems, vol. 9, pp. 281–287, 1996. View at: Google Scholar
  15. X. Feng, H. Zhao, and S. Li, “Modeling non-linear displacement time series of geo-materials using evolutionary support vector machines,” International Journal of Rock Mechanics and Mining Sciences, vol. 41, no. 7, pp. 1087–1107, 2004. View at: Publisher Site | Google Scholar
  16. H. B. Zhao, “Slope reliability analysis using a support vector machine,” Computers and Geotechnics, vol. 35, no. 3, pp. 459–467, 2008. View at: Publisher Site | Google Scholar
  17. H. Zhao and S. Yin, “Geomechanical parameters identification by particle swarm optimization and support vector machine,” Applied Mathematical Modelling, vol. 33, no. 10, pp. 3997–4012, 2009. View at: Publisher Site | Google Scholar | Zentralblatt MATH
  18. A. T. C. Goh and S. H. Goh, “Support vector machines: their use in geotechnical engineering as illustrated using seismic liquefaction data,” Computers and Geotechnics, vol. 34, no. 5, pp. 410–421, 2007. View at: Publisher Site | Google Scholar
  19. J. Shen, Y. Syau, and E. S. Lee, “Support vector fuzzy adaptive network in regression analysis,” Computers & Mathematics with Applications, vol. 54, no. 11-12, pp. 1353–1366, 2007. View at: Publisher Site | Google Scholar | MathSciNet
  20. S. Helwany, Applied Soil Mechanics with Abaqus Applications, John Wiley & Sons, Hoboken, NJ, USA, 2007.

Copyright © 2014 Hongbo Zhao et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

More related articles

 PDF Download Citation Citation
 Download other formatsMore
 Order printed copiesOrder

Related articles