Probabilistic analysis of tunnel face seismic stability in layered rock masses using Polynomial Chaos Kriging metamodel

Face stability is an essential issue in tunnel design and construction.Layered rock masses are typical and ubiquitous;uncertainties in rock properties always exist.In view of this,a comprehensive method,which combines the Upper bound Limit analysis of Tunnel face stability,the Polynomial Chaos Kriging,the Monte-Carlo Simulation and Analysis of Covariance method(ULT-PCK-MA),is proposed to investigate the seismic stability of tunnel faces.A two-dimensional analytical model of ULT is developed to evaluate the virtual support force based on the upper bound limit analysis.An efficient probabilistic analysis method PCK-MA based on the adaptive Polynomial Chaos Kriging metamodel is then implemented to investigate the parameter uncertainty effects.Ten input parameters,including geological strength indices,uniaxial compressive strengths and constants for three rock formations,and the horizontal seismic coefficients,are treated as random variables.The effects of these parameter uncertainties on the failure probability and sensitivity indices are discussed.In addition,the effects of weak layer position,the middle layer thickness and quality,the tunnel diameter,the parameters correlation,and the seismic loadings are investigated,respectively.The results show that the layer distributions significantly influence the tunnel face probabilistic stability,particularly when the weak rock is present in the bottom layer.The efficiency of the proposed ULT-PCK-MA is validated,which is expected to facilitate the engineering design and construction.

Load-bearing characteristics and energy evolution of fractured rock masses after granite and sandstone grouting

Experiments on grouting-reinforced rock mass specimens with different particle sizes and features were carried out in this study to examine the effects of grouting reinforcement on the load-bearing characteristics of fractured rock mass.The strength and deformation features of grouting-reinforced rock mass were analyzed under different loading manners;the energy evolution mechanism of grouting-reinforced rock mass specimens with different particle sizes and features was investigated;the energy dissipation ratio and post-peak stress decreasing rate were employed to evaluate the bearing stability of grouting-reinforced rock mass.The results show that the strength and ductility of granite-reinforced rock mass(GRM)under biaxial loading are higher than that of sandstone-reinforced rock mass(SRM)under uniaxial loading.Besides,the energy evolution characteristics of grouting-reinforced rock mass under uniaxial and biaxial loading mainly could be divided into early,middle,and late stages.In the early stage,total,elastic,and dissipation energies were quite small with flatter curves;in the middle stage,elastic energy increased rapidly,whereas dissipation energy increased slowly;in the late stage,dissipation energy increased sharply.The energy dissipation ratio was used to represent the pre-peak plastic deformation.Under uniaxial loading,this ratio increased as the particle size increased and the pre-peak plastic deformation of grouting-reinforced rock mass became larger;under biaxial loading,it dropped as the particle size increased,and the pre-peak plastic deformation of grouting-reinforced rock mass became smaller.The post-peak stress decline rate A_(v) was used to assess the post-peak bearing performance of grouting-reinforced rock mass.Under uniaxial loading,parameter A_(v) exhibited reduction as the particle size kept increasing,and the ability of post-peak of grouting-reinforced rock mass to allow deformation development was greater,and the bearing capacity was greater;under biaxial loading,A_(v) increased with the particle size,and the ability of post-peak of grouting-reinforced rock mass to allow deformation development was low and the bearing capacity was reduced.The findings are considered instrumental in improving the stability of the roadway-surrounding rock by granite and sandstone grouting.

Deformation characteristics and damage ontologies of soft and hard composite rock masses under impact loading

As one of the most common occurring geological landforms in deep rock formations, the dynamic mechanical properties of layered composite rock bodies under impact loading have been widely studied by scholars. To study the dynamic properties of soft and hard composite rocks with different thickness ratios, this paper utilizes cement, quartz sand and gypsum powder to construct soft and hard composite rock specimens and utilizes a combination of indoor tests, numerical calculations, and theoretical analyses to investigate the mechanical properties of soft and hard composite rock bodies. The test results reveal that:(1) When the proportion of hard rock increases from 20% to 50%, the strength of the combined rock body increases by 69.14 MPa and 87 MPa when the hard rock face and soft rock face are loaded, respectively;however, when the proportion of hard rock is the same, the compressive strength of the hard rock face impact is 9%-17% greater than that of the soft rock face impact;(2) When a specimen of soft and hard combined rock body is subjected to impact loading, the damage mode involves mixed tension and shear damage, and the cracks generally first appear at the ends of the specimen, then develop on the laminar surface from the impact surface, and finally end in the overall damage of the soft rock part. The development rate and the total number of cracks in the same specimen when the hard rock face is impacted are significantly greater than those when the soft rock face is impacted;(3) By introducing Weibull’s statistical strength theory to establish the damage variables of soft-hard combined rock bodies, combined with the DP strength criterion, the damage model and the Kelvin body are concatenated to obtain a statistical damage constitutive model, which can better fit the full stress-strain curve of soft-hard combined rock body specimens under a single impact load.

Evaluation of slope stability through rock mass classification and kinematic analysis of some major slopes along NH-1A from Ramban to Banihal, North Western Himalayas

The network of Himalayan roadways and highways connects some remote regions of valleys or hill slopes,which is vital for India’s socio-economic growth.Due to natural and artificial factors,frequency of slope instabilities along the networks has been increasing over last few decades.Assessment of stability of natural and artificial slopes due to construction of these connecting road networks is significant in safely executing these roads throughout the year.Several rock mass classification methods are generally used to assess the strength and deformability of rock mass.This study assesses slope stability along the NH-1A of Ramban district of North Western Himalayas.Various structurally and non-structurally controlled rock mass classification systems have been applied to assess the stability conditions of 14 slopes.For evaluating the stability of these slopes,kinematic analysis was performed along with geological strength index(GSI),rock mass rating(RMR),continuous slope mass rating(CoSMR),slope mass rating(SMR),and Q-slope in the present study.The SMR gives three slopes as completely unstable while CoSMR suggests four slopes as completely unstable.The stability of all slopes was also analyzed using a design chart under dynamic and static conditions by slope stability rating(SSR)for the factor of safety(FoS)of 1.2 and 1 respectively.Q-slope with probability of failure(PoF)1%gives two slopes as stable slopes.Stable slope angle has been determined based on the Q-slope safe angle equation and SSR design chart based on the FoS.The value ranges given by different empirical classifications were RMR(37-74),GSI(27.3-58.5),SMR(11-59),and CoSMR(3.39-74.56).Good relationship was found among RMR&SSR and RMR&GSI with correlation coefficient(R 2)value of 0.815 and 0.6866,respectively.Lastly,a comparative stability of all these slopes based on the above classification has been performed to identify the most critical slope along this road.

Rock mass quality prediction on tunnel faces with incomplete multi-source dataset via tree-augmented naive Bayesian network

Rock mass quality serves as a vital index for predicting the stability and safety status of rock tunnel faces.In tunneling practice,the rock mass quality is often assessed via a combination of qualitative and quantitative parameters.However,due to the harsh on-site construction conditions,it is rather difficult to obtain some of the evaluation parameters which are essential for the rock mass quality prediction.In this study,a novel improved Swin Transformer is proposed to detect,segment,and quantify rock mass characteristic parameters such as water leakage,fractures,weak interlayers.The site experiment results demonstrate that the improved Swin Transformer achieves optimal segmentation results and achieving accuracies of 92%,81%,and 86%for water leakage,fractures,and weak interlayers,respectively.A multisource rock tunnel face characteristic(RTFC)dataset includes 11 parameters for predicting rock mass quality is established.Considering the limitations in predictive performance of incomplete evaluation parameters exist in this dataset,a novel tree-augmented naive Bayesian network(BN)is proposed to address the challenge of the incomplete dataset and achieved a prediction accuracy of 88%.In comparison with other commonly used Machine Learning models the proposed BN-based approach proved an improved performance on predicting the rock mass quality with the incomplete dataset.By utilizing the established BN,a further sensitivity analysis is conducted to quantitatively evaluate the importance of the various parameters,results indicate that the rock strength and fractures parameter exert the most significant influence on rock mass quality.

Effect of eccentric and inclined loading on the bearing capacity of strip footing placed on rock mass

This paper deals with the bearing capacity determination of strip footing on a rock mass in hilly area by considering the influence of inclined and eccentric loading. Applying the generalized HoekBrown failure criterion, the failure behavior of the rock mass is modeled with the help of the power cone programming in the lower bound finite element limit analysis framework. Using bearing capacity factor(Ns), the change in bearing capacity of the strip footing due to the occurrence of eccentrically inclined loading is presented. The variations of the magnitude of Ns are obtained by examining the effects of the Hoek-Brown rock mass strength parameters(uniaxial compressive strength(sci), disturbance factor(D), rock parameter(mi), and Geological Strength Index(GSI)) in the presence of different magnitudes of eccentricity(e) and inclination angle(λ) with respect to the vertical plane, and presented as design charts. Both the inclined loading modes, i.e., inclination towards the center of strip footing(+λ) and inclination away from the center of strip footing(-λ), are adopted to perform the investigation. In addition, the correlation between the input parameters and the corresponding output is developed by utilizing the artificial neural network(ANN). Additionally, from sensitivity analysis, it is observed that inclination angle(λ) is the most sensitive parameter. For practicing engineers, the obtained design equation and design charts can be beneficial to understand the bearing capacity variation in the existence of eccentrically inclined loading in mountain areas.

Numerical manifold method for thermo-mechanical coupling simulation of fractured rock mass

As a calculation method based on the Galerkin variation,the numerical manifold method(NMM)adopts a double covering system,which can easily deal with discontinuous deformation problems and has a high calculation accuracy.Aiming at the thermo-mechanical(TM)coupling problem of fractured rock masses,this study uses the NMM to simulate the processes of crack initiation and propagation in a rock mass under the influence of temperature field,deduces related system equations,and proposes a penalty function method to deal with boundary conditions.Numerical examples are employed to confirm the effectiveness and high accuracy of this method.By the thermal stress analysis of a thick-walled cylinder(TWC),the simulation of cracking in the TWC under heating and cooling conditions,and the simulation of thermal cracking of the SwedishÄspöPillar Stability Experiment(APSE)rock column,the thermal stress,and TM coupling are obtained.The numerical simulation results are in good agreement with the test data and other numerical results,thus verifying the effectiveness of the NMM in dealing with thermal stress and crack propagation problems of fractured rock masses.

Research on unloading nonlinear mechanical characteristics of jointed rock masses

Geological environments of rock mass projects are always very complicated, and further investigations on rock mechanical characteristics are needed. There are considerable distinctions in rock mechanical characteristics under unloading and loading conditions. A series of tests are conducted to study the stress-strain relationship of rock masses under loading and unloading conditions. Also, the anisotropy, the size effect, and the rheological property of unloading rock mass are investigated. The tests presented in the paper include model test and granite rheological test, which are conducted considering geological condition, rock mass structure, in-situ stress field of the permanent shiplock of the Three Gorges Project. The main differences between loading and unloading rock masses are stress paths, yield criteria, deformation and strength parameters, etc.. Different structural plane directions affect unloading rock mass evidently. With increasing size, the tensile strength, the compressive strength, the deformation modulus, the Poisson’s ratio and the anisotropy of rock mass all decrease. For sandstone samples with parallel bedding planes, the cohesion c increases but the internal friction angle ? decreases under unloading condition when compared with the values under loading condition. While for samples with vertical bedding planes, the trend is adverse. The rheological property of rocks has close relationship with the tensile stresses of rock masses. When the sandstone samples are tested under high stress condition, their rheological properties are very obvious with the unloading of confining pressure, and three typical rheological stages are shown. Rheological rate changes with the variations in axial stress and confining pressure.

Responses of jointed rock masses subjected to impact loading

Impact-induced damage to jointed rock masses has important consequences in various mining and civil engineering applications. This paper reports a numerical investigation to address the responses of jointed rock masses subjected to impact loading. It also focuses on the static and dynamic properties of an intact rock derived from a series of laboratory tests on meta-sandstone samples from a quarry in Nova Scotia, Canada. A distinct element code(PFC2D) was used to generate a bonded particle model(BPM) to simulate both the static and dynamic properties of the intact rock. The calibrated BPM was then used to construct large-scale jointed rock mass samples by incorporating discrete joint networks of multiple joint intensities into the intact rock matrix represented by the BPM. Finally, the impact-induced damage inflicted by a rigid projectile particle on the jointed rock mass samples was determined through the use of the numerical model. The simulation results show that joints play an important role in the impactinduced rock mass damage where higher joint intensity results in more damage to the rock mass. This is mainly attributed to variations of stress wave propagation in jointed rock masses as compared to intact rock devoid of joints.

ZONAL DISINTEGRATION MECHANISM OF DEEP CRACK-WEAKENED ROCK MASSES UNDER DYNAMIC UNLOADING

Size and quantity of fractured zone and non-fractured zone are controlled by cracks contained in deep rock masses. Zonal disintegration mechanism is strongly dependent on the interaction among cracks. The strong interaction among cracks is investigated using stress superposition principle and the Chebyshev polynomials expansion of the pseudo-traction. It is found from numerical results that crack nucleation, growth and coalescence lead to failure of deep crack- weakened rock masses. The stress redistribution around the surrounding rock mass induced by unloading excavation is studied. The effect of the excavation time on nucleation, growth, interaction and coalescence of cracks was analyzed. Moreover, the influence of the excavation time on the size and quantity of fractured zone and non-fractured zone was given. When the excavation time is short, zonal disintegration phenomenon may occur in deep rock masses. It is shown from numerical results that the size and quantity of fractured zone increase with decreasing excavation time, and the size and quantity of fractured zone increase with the increasing value of in-situ geostress.

Strength of massive to moderately jointed hard rock masses

The Hoek-Brown (HB) failure criterion and the geological strength index (GSI) were developed for the estimation of rock mass strength in jointed and blocky ground where rock mass failure is dominated by sliding along open joints and rotation of rock blocks. In massive, veined and moderately jointed rock in which rock blocks cannot form without failure of intact rock, the approach to obtain HB parameters must be modified. Typical situations when these modifications are required include the design of pillars, excavation and cavern stability, strainburst potential assessment, and tunnel support in deep underground conditions (around σ1/σci > 0.15, where σ1 is the major principal compressive stress and σci is the unconfined compressive strength of the homogeneous rock) in hard brittle rocks with GSI ≥ 65. In this article, the strength of massive to moderately jointed hard rock masses is investigated, and an approach is presented to estimate the rock mass strength envelope using laboratory data from uniaxial and triaxial compressive strength tests without reliance on the HB-GSI equations. The data from tests on specimens obtained from massive to moderately jointed heterogeneous (veined) rock masses are used to obtain the rock and rock mass strengths at confining stress ranges that are relevant for deep tunnelling and mining;and a methodology is presented for this purpose from laboratory data alone. By directly obtaining the equivalent HB rock mass strength envelope for massive to moderately jointed rock from laboratory tests, the HB-GSI rock mass strength estimation approach is complemented for conditions where the GSIequations are not applicable. Guidance is also provided on how to apply the proposed approach when laboratory test data are not or not yet available.

Advances in statistical mechanics of rock masses and its engineering applications

To efficiently link the continuum mechanics for rocks with the structural statistics of rock masses,a theoretical and methodological system called the statistical mechanics of rock masses(SMRM)was developed in the past three decades.In SMRM,equivalent continuum models of stressestrain relationship,strength and failure probability for jointed rock masses were established,which were based on the geometric probability models characterising the rock mass structure.This follows the statistical physics,the continuum mechanics,the fracture mechanics and the weakest link hypothesis.A general constitutive model and complete stressestrain models under compressive and shear conditions were also developed as the derivatives of the SMRM theory.An SMRM calculation system was then developed to provide fast and precise solutions for parameter estimations of rock masses,such as full-direction rock quality designation(RQD),elastic modulus,Coulomb compressive strength,rock mass quality rating,and Poisson’s ratio and shear strength.The constitutive equations involved in SMRM were integrated into a FLAC3D based numerical module to apply for engineering rock masses.It is also capable of analysing the complete deformation of rock masses and active reinforcement of engineering rock masses.Examples of engineering applications of SMRM were presented,including a rock mass at QBT hydropower station in northwestern China,a dam slope of Zongo II hydropower station in D.R.Congo,an open-pit mine in Dexing,China,an underground powerhouse of Jinping I hydropower station in southwestern China,and a typical circular tunnel in Lanzhou-Chongqing railway,China.These applications verified the reliability of the SMRM and demonstrated its applicability to broad engineering issues associated with jointed rock masses.

Influences of strain softening and seepage on elastic and plastic solutions of circular openings in nonlinear rock masses

Considering the influence of strain softening, the solutions of stress, displacement, plastic softening region radius and plastic residual region radius were derived for circular openings in nonlinear rock masses subjected to seepage. The radial stress distribution curve, ground reaction curve, and relation curve between plastic softening region radius and supporting force in three different conditions were drawn respectively. From the comparisons among these results for different conditions, it is found that when the supporting force is the same, the displacement of tunnel wall considering both seepage and strain softening is 85.71% greater than that only considering seepage. The increase values of radial displacement at 0.95 m and plastic softening region radius at 6.6 m show that the seepage and strain softening have the most unfavorable effects on circular opening stability in strain softening rock masses.

Earthquake-induced collapse mechanism of two types of dangerous rock masses

As the economy of China develops, an increasing number of key traffic projects have been undertaken in the west of China, where there are high, steep rock slopes. The collapse of dangerous rock masses, especially following a strong earthquake, is one of common geological disasters known in rock slope engineering. Therefore, it is important to study the collapse mechanism of dangerous rock masses induced by an earthquake and the analysis approach of its stability. This study provides a simple and convenient method to determine the collapse mechanisms of two types of dangerous rock masses （i.e. cantilever and upright） associated with the definition and calculation of the safety factor, which is based on the flexure theory of a constant-section beam by combining with the maximum tensile-stress criterion to depict the process of crack propagation caused by seismic waves. The calculation results show that there are critical crack depths in each form of the dangerous rock masses. Once the accumulated depth of the crack growth during an earthquake exceeds the critical depth, the collapse will occur. It is also demonstrated that the crack extension amount of each step is not a constant value, and is closely associated with the current accumulated crack depth. The greater the cumulative crack depth, the more easily the crack propagates. Finally, the validity and applicability of the proposed method are verified through two actual engineering examples.

Determination of hydraulic conductivity of fractured rock masses:A case study for a rock cavern project in Singapore

In order to reduce the risk associated with water seepage in an underground rock cavern project inSingapore, a reliable hydro-geological model should be established based on the in situ investigationdata. The key challenging issue in the hydro-geological model building is how to integrate limitedgeological and hydro-geological data to determine the hydraulic conductivity of the fractured rockmasses. Based on the data obtained from different stages （feasibility investigation stage, constructionstage, and post-construction stage）, suitable models and methods are proposed to determine the hydraulicconductivities at different locations and depths, which will be used at other locations in thefuture. 2015 Institute of Rock and Soil Mechanics, Chinese Academy of Sciences. Production and hosting byElsevier B.V. All rights reserved.

Shear strength criteria for rock,rock joints,rockfill and rock masses:Problems and some solutions

Although many intact rock types can be very strong,a critical confining pressure can eventually be reached in triaxial testing,such that the Mohr shear strength envelope becomes horizontal.This critical state has recently been better defined,and correct curvature or correct deviation from linear Mohr-Coulomb（MC） has finally been found.Standard shear testing procedures for rock joints,using multiple testing of the same sample,in case of insufficient samples,can be shown to exaggerate apparent cohesion.Even rough joints do not have any cohesion,but instead have very high friction angles at low stress,due to strong dilation.Rock masses,implying problems of large-scale interaction with engineering structures,may have both cohesive and frictional strength components.However,it is not correct to add these,following linear M-C or nonlinear Hoek-Brown（H-B） standard routines.Cohesion is broken at small strain,while friction is mobilized at larger strain and remains to the end of the shear deformation.The criterion ＇c then σn tan φ＇ should replace ＇c plus σn tan φ＇ for improved fit to reality.Transformation of principal stresses to a shear plane seems to ignore mobilized dilation,and caused great experimental difficulties until understood.There seems to be plenty of room for continued research,so that errors of judgement of the last 50 years can be corrected.

Face stability analysis of circular tunnels in layered rock masses using the upper bound theorem

An analysis of tunnel face stability generally assumes a single homogeneous rock mass.However,most rock tunnel projects are excavated in stratified rock masses.This paper presents a two-dimensional(2D)analytical model for estimating the face stability of a rock tunnel in the presence of rock mass stratification.The model uses the kinematical limit analysis approach combined with the block calculation technique.A virtual support force is applied to the tunnel face,and then solved using an optimization method based on the upper limit theorem of limit analysis and the nonlinear Hoek-Brown yield criterion.Several design charts are provided to analyze the effects of rock layer thickness on tunnel face stability,tunnel diameter,the arrangement sequence of weak and strong rock layers,and the variation in rock layer parameters at different positions.The results indicate that the thickness of the rock layer,tunnel diameter,and arrangement sequence of weak and strong rock layers significantly affect the tunnel face stability.Variations in the parameters of the lower layer of the tunnel face have a greater effect on tunnel stability than those of the upper layer.

Computing in-situ strength of rock masses based upon RQD and modified joint factor: Using pressure and damage sensitive constitutive relationship

In th is study, a n e w m odel w as p re se n te d for com p u tin g stre n g th o f rock m asses based u p o n in -situo bservations o f RQD p o pularly kno w n as rock quality d esignation. This m odel links u p th e rock m assp aram eters from in -situ investigations w ith th e stre n g th p a ram eters o f jo in ted rocks obtain ed fromlaboratory scale ex p erim en tal observations. Using th e co n stitu tiv e relation, th e a u th o r derived a p ressu reand d am age sensitive plastic p a ra m e te r to d ete rm in e stre n g th o f rock m asses for varied ex te n ts ofd isco n tin u ity an d p ressu re induced dam age. The te s t results show th a t plasticity characterized byhard en in g an d softening inclusive o f dam ag e invariably d e p en d s u p o n m ean p ressu re an d e x te n t ofdefo rm atio n s alread y experien ced by rock m asses. The p re se n t w ork explores th e te s t d a ta th a t revealth e d ep en d en c e o f in -situ stren g th on increm ental jo in t p ara m e te rs o b tain ed from th e jo in t num ber,jo in t orien tatio n , jo in t roughness, gouge p a ram eters an d w a te r pressure. S ubstituting th e relationshipb e tw e e n th e RQD and m odified jo in t factor w ith th a t b e tw e e n m odulus ratio an d stren g th ratio, th em odel show s successfully th a t using d am age inclusive plastic p a ra m e te r an d RQD provides a relationshipfor estim atin g th e stre n g th o f rock m asses. One o f th e m ain objectives o f this w ork is to illustrate th a t th ep re se n t m odel is sensitive to p la s tic ity a n d dam ag e to g e th e r in estim atin g in -situ stre n g th o f rock m assesin foundations, u n d e rg ro u n d excavation an d tunnels.

THE CONSTITUTIVE RELATION OF CRACK-WEAKENED ROCK MASSES UNDER AXIAL-DIMENSIONAL UNLOADING

An accurate and efficient numerical method for solving the crack-crack interaction problem is presented. The method is mainly by means of the dislocation model, stress superposition principle and Chebyshev polynomial expansion of the pseudo-traction. This method can be applied to compute the stress intensity factors of multiple kinked cracks and multiple rows of periodic cracks as well as the overall strains of rock masses containing multiple kinked cracks under complex loads. Many complex computational examples are given. The dependence of the crack-crack interaction on the crack configuration, the geometrical and physical parameters, and loads pattern, is investigated. By comparison with numerical results under confining pressure unloading, it is shown that the crack-crack interaction under axial-dimensional unloading is weaker than those under confining pressure unloading. Numerical results for single faults and crossed faults show that the single faults are more unstable than the crossed faults. It is found from numerical results for different crack lengths and different crack spacing that the interaction among kinked cracks decreases with an increase in length of the kinked cracks and the crack spacing under axial-dimensional unloading.

New design equations for estimation of ultimate bearing capacity of shallow foundations resting on rock masses

Rock masses are commonly used as the underlying layer of important structures such as bridges, dams and transportation constructions. The success of a foundation design for such structures mainly depends on the accuracy of estimating the bearing capacity of rock beneath them. Several traditional numerical approaches are proposed for the estimation of the bearing capacity of foundations resting on rock masses to avoid performing elaborate and expensive experimental studies. Despite this fact, there still exists a serious need to develop more robust predictive models. This paper proposes new nonlinear prediction models for the ultimate bearing capacity of shallow foundations resting on non-fractured rock masses using a novel evolutionary computational approach, called linear genetic programming. A comprehensive set of rock socket, centrifuge rock socket, plate load and large-scaled footing load test results is used to develop the models. In order to verify the validity of the models, the sensitivity analysis is conducted and discussed. The results indicate that the proposed models accurately characterize the bearing capacity of shallow foundations. The correlation coefficients between the experimental and predicted bearing capacity values are equal to 0.95 and 0.96 for the best LGP models. Moreover, the derived models reach a notably better prediction performance than the traditional equations.