Theoretical analysis of hydrogen solubility in direct coal liquefaction solvents

The cyclic hydrogenation technology in a direct coal liquefaction process relies on the dissolved hydrogen of the solvent or oil participating in the hydrogenation reaction.Thus,a theoretical basis for process optimization and reactor design can be established by analyzing the solubility of hydrogen in liquefaction solvents.Experimental studies of hydrogen solubility in liquefaction solvents are challenging due to harsh reaction conditions and complex solvent compositions.In this study,the composition and content of liquefied solvents were analyzed.As model compounds,hexadecane,toluene,naphthalene,tetrahydronaphthalene,and phenanthrene were chosen to represent the liquefied solvents in chain alkanes and monocyclic,bicyclic,and tricyclic aromatic hydrocarbons.The solubility of hydrogen X(mol/mol)in pure solvent components and mixed solvents(alkanes and aromatics mixed in proportion to the chain alkanes+bicyclic aromatic hydrocarbons,bicyclic saturated aromatic hydrocarbons+bicyclic aromatic hydrocarbons,and bicyclic aromatic hydrocarbons+compounds containing het-eroatoms composed of mixed components)are determined using Aspen simulation at temperature and pressure conditions of 373–523 K and 2–10 MPa.The results demonstrated that at high temperatures and pressures,the solubility of hydrogen in the solvent increases with the increase in temperature and pressure,with the pressure having a greater impact.Further-more,the results revealed that hydrogen is more soluble in straight-chain alkanes than in other solvents,and the solubility of eicosanoids reaches a maximum of 0.296.The hydrogen solubility in aromatic ring compounds decreased gradually with an increase in the aromatic ring number.The influence of chain alkanes on the solubility of hydrogen predominates in a mixture of solvents with different mixing ratios of chain alkanes and aromatic hydrocarbons.The solubility of hydrogen in mixed aromatic solvents is less than that in the corresponding single solvents.Hydrogen is less soluble in solvent compounds containing heteroatoms than in compounds without heteroatoms.

Numerical and experimental study on the falling film flow characteristics with the effect of co-current gas flow in hydrogen liquefaction process

Liquid hydrogen storage and transportation is an effective method for large-scale transportation and utilization of hydrogen energy. Revealing the flow mechanism of cryogenic working fluid is the key to optimize heat exchanger structure and hydrogen liquefaction process(LH2). The methods of cryogenic visualization experiment, theoretical analysis and numerical simulation are conducted to study the falling film flow characteristics with the effect of co-current gas flow in LH2spiral wound heat exchanger.The results show that the flow rate of mixed refrigerant has a great influence on liquid film spreading process, falling film flow pattern and heat transfer performance. The liquid film of LH2mixed refrigerant with column flow pattern can not uniformly and completely cover the tube wall surface. As liquid flow rate increases, the falling film flow pattern evolves into sheet-column flow and sheet flow, and liquid film completely covers the surface of tube wall. With the increase of shear effect of gas-phase mixed refrigerant in the same direction, the liquid film gradually becomes unstable, and the flow pattern eventually evolves into a mist flow.

Liquefaction susceptibility and deformation characteristics of saturated coral sandy soils subjected to cyclic loadings-a critical review

Coral sandy soils widely exist in coral island reefs and seashores in tropical and subtropical regions.Due to the unique marine depositional environment of coral sandy soils,the engineering characteristics and responses of these soils subjected to monotonic and cyclic loadings have been a subject of intense interest among the geotechnical and earthquake engineering communities.This paper critically reviews the progress of experimental investigations on the undrained behavior of coral sandy soils under monotonic and cyclic loadings over the last three decades.The focus of coverage includes the contractive-dilative behavior,the pattern of excess pore-water pressure(EPWP)generation and the liquefaction mechanism and liquefaction resistance,the small-strain shear modulus and strain-dependent shear modulus and damping,the cyclic softening feature,and the anisotropic characteristics of undrained responses of saturated coral sandy soils.In particular,the advances made in the past decades are reviewed from the following aspects:(1)the characterization of factors that impact the mechanism and patterns of EPWP build-up;(2)the identification of liquefaction triggering in terms of the apparent viscosity and the average flow coefficient;(3)the establishment of the invariable form of strain-based,stress-based,or energy-based EPWP ratio formulas and the unique relationship between the new proxy of liquefaction resistance and the number of cycles required to reach liquefaction;(4)the establishment of the invariable form of the predictive formulas of small strain modulus and strain-dependent shear modulus;and(5)the investigation on the effects of stress-induced anisotropy on liquefaction susceptibility and dynamic deformation characteristics.Insights gained through the critical review of these advances in the past decades offer a perspective for future research to further resolve the fundamental issues concerning the liquefaction mechanism and responses of coral sandy sites subjected to cyclic loadings associated with seismic events in marine environments.

Effect of Wave Nonlinearity on the Instantaneous Seabed Liquefaction

The nonlinear variation of wave is commonly seen in nearshore area,and the resulting seabed response and liquefaction are of high concern to coastal engineers.In this study,an analytical formula considering the nonlinear wave skewness and asymmetry is adopted to provide wave pressure on the seabed surface.The liquefaction depth attenuation coefficient and width growth coefficient are defined to quantitatively characterize the nonlinear effect of wave on seabed liquefaction.Based on the 2D full dynamic model of wave-induced seabed response,a detailed parametric study is carried out in order to evaluate the influence of the nonlinear variation of wave loadings on seabed liquefaction.Further,new empirical prediction formulas are proposed to fast predict the maximum liquefaction under nonlinear wave.Results indicate that(1)Due to the influence of wave nonlinearity,the vertical transmission of negative pore water pressure in the seabed is hindered,and therefore,the amplitude decreases significantly.(2)In general,with the increase of wave nonlinearity,the liquefaction depth of seabed decreases gradually.Especially under asymmetric and skewed wave loading,the attenuation of maximum seabed liquefaction depth is the most significant among all the nonlinear wave conditions.However,highly skewed wave can cause the liquefaction depth of seabed greater than that under linear wave.(3)The asymmetry of wave pressure leads to the increase of liquefaction width,whereas the influence of skewedness is not significant.(4)Compared with the nonlinear waveform,seabed liquefaction is more sensitive to the variation of nonlinear degree of wave loading.

The Analysis of the Correlation between SPT and CPT Based on CNN-GA and Liquefaction Discrimination Research

The objective of this study is to investigate themethods for soil liquefaction discrimination. Typically, predicting soilliquefaction potential involves conducting the standard penetration test (SPT), which requires field testing and canbe time-consuming and labor-intensive. In contrast, the cone penetration test (CPT) provides a more convenientmethod and offers detailed and continuous information about soil layers. In this study, the feature matrix based onCPT data is proposed to predict the standard penetration test blow count N. The featurematrix comprises the CPTcharacteristic parameters at specific depths, such as tip resistance qc, sleeve resistance f s, and depth H. To fuse thefeatures on the matrix, the convolutional neural network (CNN) is employed for feature extraction. Additionally,Genetic Algorithm (GA) is utilized to obtain the best combination of convolutional kernels and the number ofneurons. The study evaluated the robustness of the proposed model using multiple engineering field data sets.Results demonstrated that the proposed model outperformed conventional methods in predicting N values forvarious soil categories, including sandy silt, silty sand, and clayey silt. Finally, the proposed model was employedfor liquefaction discrimination. The liquefaction discrimination based on the predicted N values was comparedwith the measured N values, and the results showed that the discrimination results were in 75% agreement. Thestudy has important practical application value for foundation liquefaction engineering. Also, the novel methodadopted in this research provides new ideas and methods for research in related fields, which is of great academicsignificance.

Mechanisms to explain soil liquefaction triggering,development,and persistence during an earthquake

Mechanisms have been proposed to explain the triggering,development,and persistence of soil liquefaction.The mechanism explaining the horizontal failure plane(triggering)and its depth below the phreatic surface is governed by the flux properties and effective stress at that plane.At the failure plane,the pore water pressure was higher than the effective stress,and the volume change was the highest.The pore water pressure is a function of the soil profile features(particularly the phreatic zone width)and bedrock motion(horizontal acceleration).The volume change at the failure plane is a function of the intrinsic permeability of the soil and bedrock displacement.The failure plane was predicted to occur during the oscillation with the highest amplitude,disregarding further bedrock motion,which was consistent with low seismic energy densities.Two mechanisms were proposed to explain the persistence of soil liquefaction.The first is the existence of low-permeability layers in the depth range in which the failure planes are predicted to occur.The other allows for the persistence and development of soil liquefaction;it is consistent with homogeneous soils and requires water inflow from bedrock water springs.The latter explains many of the features of soil liquefaction observed during earthquakes,namely,surficial effects,“instant”liquefaction,and the occurrence of short-and long-term changes in the level of the phreatic surfaces.This model(hypothesis),the relationship between the flux characteristics and loss of soil shear strength,provides self-consistent constraints on the depth below the phreatic surfaces where the failure planes are observed(expected to occur).It requires further experimental and observational evidence.Similar reasoning can be used to explain other saturated soil phenomena.

Liquefaction and post-liquefaction behaviors of sands affected by immersion-induced degradation of crushed mudstone

A series of undrained triaxial tests was conducted to investigate the effect of crushed mudstone with the immersion-induced degradation on the liquefaction and post-liquefaction properties,and the undrained shearing behavior without precedent cyclic-loading histories of sands containing crushed mudstone.The tested materials with a main particle diameter of 2-0.85 mm were prepared by mixing sands and crushed mudstone to reach the prescribed mudstone content defined by dry mass ranging from 0% to 50%.The mixtures were subjected to immersion under a certain stress level and were subsequently tested.In addition,one-dimensional compression tests were also supplementally performed to visually observe the immersion-induced degradation of crushed mudstone.The test results mainly showed that: (1) the liquefaction resistance,the post-liquefaction undrained strength,and the undrained strength without a precedent cyclic-loading history decreased significantly with increasing mudstone content,M c ,up to 20%;(2) even a small amount of crushed mudstone affected these strengths;(3) the above-mentioned large reductions in the strengths were attributed to the immersion-induced degradation of crushed mudstone;(4) at M_(c) >20%,the liquefaction resistance increased while the significant increase in the undrained static strengths with and without precedent cyclic-loading histories was not observed;and (5) the increase in the liquefaction resistance at M_(c) >20% may have been attributed to both the gradual increase in the plasticity and the formation of the soil aggregates among deteriorated crushed mudstone,while the increase in the specimen density did not play an important role in such behavior.

Effect of particle composition and consolidation degree on the wave-induced liquefaction of soil beds

The wave-induced liquefaction of seabed is responsible for causing damage to marine structures.Particle composition and consolidation degree are the key factors affecting the pore water pressure response and liquefaction behavior of the seabed under wave action.The present study conducted wave flume experiments on silt and silty fine sand beds with varying particle compositions.Furthermore,a comprehensive analysis of the differences and underlying reasons for liquefaction behavior in two different types of soil was conducted from both macroscopic and microscopic perspectives.The experimental results indicate that the silt bed necessitates a lower wave load intensity to attain the liquefaction state in comparison to the silty fine sand bed.Additionally,the duration and development depth of liquefaction are greater in the silt bed.The dissimilarity in liquefaction behavior between the two types of soil can be attributed to the variation in their permeability and plastic deformation capacity.The permeability coefficient and compression modulus of silt are lower than those of silty fine sand.Consequently,silt is more prone to the accumulation of pore pressure and subsequent liquefaction under external loading.Prior research has demonstrated that silt beds with varying consolidation degrees exhibit distinct initial failure modes.Specifically,a dense bed undergoes shear failure,whereas a loose bed experiences initial liquefaction failure.This study utilized discrete element simulation to examine the microscopic mechanisms that underlie this phenomenon.

Storm Liquefaction Deposits:A Possibility of Time Reversal in Sedimentary Strata of an Estuarial Coastal Area

In the sedimentary strata dating of estuarial coastal areas,it is often found that there is phenomenon of time-reversal in strata.The seabed sediments could be liquefied under storm waves.A laboratory wave flume experiment demonstrated that storm-induced liquefaction deposits are formed by the oscillations of liquefied sediments.In this paper,the particle size distribution and ^(210)Pb_(ex) specific activity of the sediment samples from the liquefaction disturbed zone and adjacent stable zone of the Yellow River Delta were tested.The stratigraphic divisions based on storm liquefaction deposit sequence can effectively explain the vertical changes in particle size and ^(210)Pb_(ex) specific activity.Due to the differentiation of particles during the storm induced liquefaction,coarse and fine particles regrouped,which could explain the phenomenon of time-reversal in dating data.

Assessment of liquefaction potential based on shear wave velocity:Strain energy approach

Liquefaction assessment based on strain energy is significantly superior to conventional stress-based methods.The main purpose of the present study is to investigate the correlation between shear wave velocity and strain energy capacity of silty sands.The dissipated energy until liquefaction occurs was calculated by analyzing the results of three series of comprehensive cyclic direct simple shear and triaxial tests on Ottawa F65,Nevada,and Firoozkuh sands with varying silt content by weight and relative densities.Additionally,the shear wave velocity of each series was obtained using bender element or resonant column tests.Consequently,for the first time,a liquefaction triggering criterion,relating to effective overburden normalized liquefaction capacity energy(WL=s’c)to effective overburden stresscorrected shear wave velocity(eVs1)has been introduced.The accuracy of the proposed criteria was evaluated using in situ data.The results confirm the ability of shear wave velocity as a distinguishing parameter for separating liquefied and non-liquefied soils when it is calculated against liquefaction capacity energy(WL=s’c).However,the proposed WL=s’c-Vs1 curve,similar to previously proposed cyclic resistance ratio(CRR)-Vs1 relationships,should be used conservatively for fields vulnerable to liquefaction-induced lateral spreading.

The Effect of Preloading on the Cyclic Liquefaction Strength Measured in the Laboratory

The effect of preloading on the liquefaction cyclic strength was investigated by cyclic shear tests where horizontal shear stress oscillated about a zero mean value on sands with varying fines content and at varying prestress ratios, densities and verticalstresses. Test results showed a marked increase of the cyclic soil strength with the prestress ratio. The effect is more pronounced for the looser specimens. An empirical expression predicting this effect is proposed. This expression is validated from results of a field test.

Constitutive modelling of fabric effect on sand liquefaction

Sand liquefaction under static and dynamic loading can cause failure of embankments,slopes,bridges and other important infrastructure.Sand liquefaction in the seabed can also cause submarine landslides and tsunamis.Fabric anisotropy related to the internal soil structure such as particle orientation,force network and void space is found to have profound influence on sand liquefaction.A constitutive model accounting for the effect of anisotropy on sand liquefaction is proposed.Evolution of fabric anisotropy during loading is considered according to the anisotropic critical state theory for sand.The model has been validated by extensive test results on Toyoura sand with different initial densities and stress states.The effect of sample preparation method on sand liquefaction is qualitatively analysed.The model has been used to investigate the response of a sand ground under earthquake loading.It is shown that sand with horizontal bedding plane has the highest resistance to liquefaction when the sand deposit is anisotropic,which is consistent with the centrifuge test results.The initial degree of fabric anisotropy has a more significant influence on the liquefaction resistance.Sand with more anisotropic fabric that can be caused by previous loading history or compaction methods has lower liquefaction resistance.

Liquefaction proneness of stratified sand-silt layers based on cyclic triaxial tests

Most studies on liquefaction have addressed homogeneous soil strata using sand or sand with fine content without considering soil stratification.In this study,cyclic triaxial tests were conducted on the stratified sand specimens embedded with the silt layers to investigate the liquefaction failures and void-redistribution at confining stress of 100 kPa under stress-controlled mode.The loosening of underlying sand mass and hindrance to pore-water flow caused localized bulging at the sand-silt interface.It is observed that at a silt thickness of 0.2H(H is the height of the specimen),nearly 187 load cycles were required to attain liquefaction,which was the highest among all the silt thicknesses with a single silt layer.Therefore,0.2H is assumed as the optimum silt thickness(t_(opt)).The silt was placed at the top,middle and bottom of the specimen to understand the effect of silt layer location.Due to the increase in depth of the silt layer from the top position(capped soil state)to the bottom,the cycles to reach liquefaction(N_(cyc,L))increased 2.18 times.Also,when the number of silt layers increased from single to triple,there was an increase of about 880%in N_(cyc,L).The micro-characterization analysis of the soil specimens indicated silty materials transported in upper sections of the specimen due to the dissipated pore pressure.The main parameters,including thickness(t),location(z),cyclic stress ratio(CSR),number of silt layers(n)and modified relative density(D_(r,m)),performed significantly in governing the lique-faction resistance.For this,a multilinear regression model is developed based on critical parameters for prediction of N_(cyc,L).Furthermore,the developed constitutive model has been validated using the data from the present study and earlier findings.

Novel Soft ComputingModel for Predicting Blast-Induced Ground Vibration in Open-Pit Mines Based on the Bagging and Sibling of Extra Trees Models

This study considered and predicted blast-induced ground vibration(PPV)in open-pit mines using bagging and sibling techniques under the rigorous combination of machine learning algorithms.Accordingly,four machine learning algorithms,including support vector regression(SVR),extra trees(ExTree),K-nearest neighbors(KNN),and decision tree regression(DTR),were used as the base models for the purposes of combination and PPV initial prediction.The bagging regressor(BA)was then applied to combine these base models with the efforts of variance reduction,overfitting elimination,and generating more robust predictive models,abbreviated as BA-ExTree,BAKNN,BA-SVR,and BA-DTR.It is emphasized that the ExTree model has not been considered for predicting blastinduced ground vibration before,and the bagging of ExTree is an innovation aiming to improve the accuracy of the inherently ExTree model,as well.In addition,two empirical models(i.e.,USBM and Ambraseys)were also treated and compared with the bagging models to gain a comprehensive assessment.With this aim,we collected 300 blasting events with different parameters at the Sin Quyen copper mine(Vietnam),and the produced PPV values were also measured.They were then compiled as the dataset to develop the PPV predictive models.The results revealed that the bagging models provided better performance than the empirical models,except for the BA-DTR model.Of those,the BA-ExTree is the best model with the highest accuracy(i.e.,88.8%).Whereas,the empirical models only provided the accuracy from 73.6%–76%.The details of comparisons and assessments were also presented in this study.

Seismic Liquefaction Resistance Based on Strain Energy Concept Considering Fine Content Value Effect and Performance Parametric Sensitivity Analysis

Liquefaction is one of the most destructive phenomena caused by earthquakes,which has been studied in the issues of potential,triggering and hazard analysis.The strain energy approach is a common method to investigate liquefaction potential.In this study,two Artificial Neural Network(ANN)models were developed to estimate the liquefaction resistance of sandy soil based on the capacity strain energy concept(W)by using laboratory test data.A large database was collected from the literature.One group of the dataset was utilized for validating the process in order to prevent overtraining the presented model.To investigate the complex influence of fine content(FC)on liquefaction resistance,according to previous studies,the second database was arranged by samples with FC of less than 28%and was used to train the second ANN model.Then,two presented ANN models in this study,in addition to four extra available models,were applied to an additional 20 new samples for comparing their results to show the capability and accuracy of the presented models herein.Furthermore,a parametric sensitivity analysis was performed through Monte Carlo Simulation(MCS)to evaluate the effects of parameters and their uncertainties on the liquefaction resistance of soils.According to the results,the developed models provide a higher accuracy prediction performance than the previously publishedmodels.The sensitivity analysis illustrated that the uncertainties of grading parameters significantly affect the liquefaction resistance of soils.

Application of machine learning to the Vs-based soil liquefaction potential assessment

Earthquakes can cause violent liquefaction of the soil, resulting in unstable foundations that can cause serious damage to facilities such as buildings, roads, and dikes. This is a primary cause of major earthquake disasters. Therefore, the discrimination and prediction of earthquake-induced soil liquefaction has been a hot issue in geohazard research. The soil liquefaction assessment is an integral part of engineering practice. This paper evaluated a dataset of 435 seismic sand liquefaction events using machine learning algorithms. The dataset was analyzed using seven potential assessment parameters. Ten machine learning algorithms are evaluated for their ability to assess seismic sand liquefaction potential, including Linear Discriminant Analysis(LDA), Quadratic Discriminant Analysis(QDA), Naive Bayes(NB), KNearest Neighbor(KNN), Artificial Neural Network(ANN), Classification Tree(CT), Support Vector Machine(SVM), Random Forest(RF), e Xtreme Gradient Boosting(XGBoost), Light Gradient Boosting Machine(Light GBM). A 10-fold cross-validation(CV) method was used in the modeling process to verify the predictive performance of the machine learning models. The final percentages of significant parameters that influenced the prediction results were obtained as Cyclic Stress Ratio(CSR) and Shear-Wave Velocity( VS1) with 56% and 38%, respectively. The final machine learning algorithms identified as suitable for seismic sand liquefaction assessment were the CT, RF, XGBoost algorithms, with the RF algorithm performing best.

Kinetics and Process Studies of the Potential for Transformation of Biogas to Biomethane and Liquefaction using Cryogenic Liquid for Domestic Applications

The present work dealt with the generation, purifying and liquefaction of biomethane to improve energy density using local materials for domestic applications. Cow dung was sourced at JKUAT dairy farm and experiments were conducted at JKUAT Bioenergy laboratory using biogas generated in laboratory scale 1 m3 bioreactors. Experiments were done in triplicates and repeated under different conditions to get the optimal conditions. The results showed that enhanced cow dung substrate displayed an improved fermentation process with increased biogas yields. Purified biogas optimized methane content from 56% ± 0.18% for raw biogas to 95% ± 0.98% for biomethane which was ideal for liquefaction.

Numerical Study on Mechanism of Blast-Induced Damage Considering Guiding Effect of Water Jet Slot

Damage is one of the most important characteristics of rock failure.Studying the damage mechanism of rock blasting under the guiding effect of the water jet slot and revealing the mechanism of controlled blasting with water jet assistance are crucial.In this study,a rock-like material was chosen as the research object for the calibration experiment of the numerical model.The numerical simulation models were then established by ANSYS/LS-DYNA,and the blastinduced damage mechanism under the guiding effect of the water jet slot was analyzed according to the blasting theory.The results indicated that explosive energy accumulates toward the direction of the slot as the guiding effect of the water jet slot,which allows the rock mass in the direction of the slot bear more damage.Meanwhile,the rock mass in the middle of the connection line between two blast-holes bears more damage under the combination of the effect of the explosion stress wave and guiding effect of water the jet slot on the detonation gas during double-slotted borehole blasting,which results in the formation of a gourd-shaped blast-induced damage area.In addition,the influence of the water jet slot on blast-induced damage varies depending on the blasting-process stage.

Effects of non-liquefiable crust layer and superstructure mass on the response of 2×2 pile groups to liquefaction-induced lateral spreading

In this research,two shake table experiments were conducted to study the effects of non-liquefiable crust layer and superstructure mass on the responses of two sets of 22 pile groups to liquefactioninduced lateral spreading.In this regard,an inclined base layer overlain by a very loose liquefiable layer was constructed in both models;while only in one model,a non-liquefiable crust layer was built.A lumped mass,being representative of a superstructure,was attached to the cap of one pile group in both models.The models were fully instrumented with various sensors,including acceleration,displacement,and pore water pressure transducers.Also,the piles were instrumented with pair strain gauges to measure pure bending moments induced by cyclic and monotonic loadings associated with ground shaking and lateral spreading,respectively.The results showed that the existence of the non-liquefiable crust layer increases both the maximum and residual soil displacements at the free field and also the maximum bending moments in the piles.The results of the experiments indicated that the crust layer induces a high kinematic lateral soil pressure and force on the piles which are not present in the crustless case.The crust layer increases the pile cap displacement before liquefaction,albeit decreases it after liquefaction,due to the elastic rebound of the piles in the liquefiable layer.The crust layer postpones both liquefaction triggering and dissipation of excess pore water pressure.The existence of the superstructure mass on the pile caps decreases the acceleration amplitude of the pile caps,while increases their maximum displacement.

Calibration of an elastoplastic model of sand liquefaction using the swarm intelligence with a multi-objective function

According to post-seismic observations,spectacular examples of engineering failures can be ascribed to the occurrence of sand liquefaction,where a sandy soil stratum could undergo a transient loss of shear strength and even behave as a“liquid”.Therefore,correct simulation of liquefaction response has become a challenging issue in geotechnical engineering field.In advanced elastoplastic models of sand liquefaction,certain fitting parameters have a remarkable effect on the computed results.However,the identification of these parameters,based on the experimental data,is usually intractable and sometimes follows a subjective trial-and-error procedure.For this,this paper presented a novel calibration methodology based on an optimization algorithm(particle swarm optimization(PSO))for an advanced elastoplastic constitutive model.A multi-objective function was designed to adjust the global quality for both monotonic and cyclic triaxial simulations.To overcome computational problem probably appearing in simulation of the cyclic triaxial test,two interrupt mechanisms were designed to prevent the particles from wasting time in searching the unreasonable space of candidate solutions.The Dafalias model has been used as an example to demonstrate the main programme.With the calibrated parameters for the HN31 sand,the computed results were highly consistent with the laboratory experiments(including monotonic triaxial tests under different confining pressures and cyclic triaxial tests in two loading modes).Finally,an extension example is given for Ottawa sand F65,suggesting that the proposed platform is versatile and can be easily customized to meet different practical needs.