Transformation of long-period stacking ordered structures in Mg-Gd-Y-Zn alloys upon synergistic characterization of first-principles calculation and experiment and its effects on mechanical properties

Based on experiments and first-principles calculations,the microstructures and mechanical properties of as-cast and solution treated Mg-10Gd-4Y-xZn-0.6Zr(x=0,1,2,wt.%)alloys are investigated.The transformation process of long-period stacking ordered(LPSO)structure during solidification and heat treatment and its effect on the mechanical properties of experimental alloys are discussed.Results reveal that the stacking faults and 18R LPSO phases appear in the as-cast Mg-10Gd-4Y-1Zn-0.6Zr and Mg-10Gd-4Y-2Zn-0.6Zr alloys,respectively.After solution treatment,the stacking faults and 18R LPSO phase transform into 14H LPSO phase.The Enthalpies of formation and reaction energy of 14H and 18R LPSO are calculated based on first-principles.Results show that the alloying ability of 18R is stronger than that of 14H.The reaction energies show that the 14H LPSO phase is more stable than the 18R LPSO.The elastic properties of the 14H and 18R LPSO phases are also evaluated by first-principles calculations,and the results are in good agreement with the experimental results.The precipitation of LPSO phase improves the tensile strength,yield strength and elongation of the alloy.After solution treatment,the Mg-10Gd-4Y-2Zn-0.6Zr alloy has the best mechanical properties,and its ultimate tensile strength and yield strength are 278.7 MPa and 196.4 MPa,respectively.The elongation of Mg-10Gd-4Y-2Zn-0.6Zr reaches 15.1,which is higher than that of Mg-10Gd-4Y0.6Zr alloy.The improving mechanism of elastic modulus by the LPSO phases and the influence on the alloy mechanical properties are also analyzed.

Predictive modeling for postoperative delirium in elderly patients with abdominal malignancies using synthetic minority oversampling technique

BACKGROUND Postoperative delirium,particularly prevalent in elderly patients after abdominal cancer surgery,presents significant challenges in clinical management.AIM To develop a synthetic minority oversampling technique(SMOTE)-based model for predicting postoperative delirium in elderly abdominal cancer patients.METHODS In this retrospective cohort study,we analyzed data from 611 elderly patients who underwent abdominal malignant tumor surgery at our hospital between September 2020 and October 2022.The incidence of postoperative delirium was recorded for 7 d post-surgery.Patients were divided into delirium and non-delirium groups based on the occurrence of postoperative delirium or not.A multivariate logistic regression model was used to identify risk factors and develop a predictive model for postoperative delirium.The SMOTE technique was applied to enhance the model by oversampling the delirium cases.The model’s predictive accuracy was then validated.RESULTS In our study involving 611 elderly patients with abdominal malignant tumors,multivariate logistic regression analysis identified significant risk factors for postoperative delirium.These included the Charlson comorbidity index,American Society of Anesthesiologists classification,history of cerebrovascular disease,surgical duration,perioperative blood transfusion,and postoperative pain score.The incidence rate of postoperative delirium in our study was 22.91%.The original predictive model(P1)exhibited an area under the receiver operating characteristic curve of 0.862.In comparison,the SMOTE-based logistic early warning model(P2),which utilized the SMOTE oversampling algorithm,showed a slightly lower but comparable area under the curve of 0.856,suggesting no significant difference in performance between the two predictive approaches.CONCLUSION This study confirms that the SMOTE-enhanced predictive model for postoperative delirium in elderly abdominal tumor patients shows performance equivalent to that of traditional methods,effectively addressing data imbalance.

Volume-averaged modeling of multiphase solidification with equiaxed crystal sedimentation in a steel ingot

Macrosegregation is a critical factor that limits the mechanical properties of materials.The impact of equiaxed crystal sedimentation on macrosegregation has been extensively studied,as it plays a significant role in determining the distribution of alloying elements and impurities within a material.To improve macrosegregation in steel connecting shafts,a multiphase solidification model that couples melt flow,heat transfer,microstructure evolution,and solute transport was established based on the volume-averaged Eulerian-Eulerian approach.In this model,the effects of liquid phase,equiaxed crystals,columnar dendrites,and columnar-to-equiaxed transition(CET)during solidification and evolution of microstructure can be considered simultaneously.The sedimentation of equiaxed crystals contributes to negative macrosegregation,where regions between columnar dendrites and equiaxed crystals undergo significant A-type positive macrosegregation due to the CET.Additionally,noticeable positive macrosegregation occurs in the area of final solidification in the ingot.The improvement in macrosegregation is beneficial for enhancing the mechanical properties of connecting shafts.To mitigate the thermal convection of molten steel resulting from excessive superheating,reducing the superheating during casting without employing external fields or altering the design of the ingot mold is indeed an effective approach to control macrosegregation.

Background removal from global auroral images:Data-driven dayglow modeling

Global images of auroras obtained by cameras on spacecraft are a key tool for studying the near-Earth environment.However,the cameras are sensitive not only to auroral emissions produced by precipitating particles,but also to dayglow emissions produced by photoelectrons induced by sunlight.Nightglow emissions and scattered sunlight can contribute to the background signal.To fully utilize such images in space science,background contamination must be removed to isolate the auroral signal.Here we outline a data-driven approach to modeling the background intensity in multiple images by formulating linear inverse problems based on B-splines and spherical harmonics.The approach is robust,flexible,and iteratively deselects outliers,such as auroral emissions.The final model is smooth across the terminator and accounts for slow temporal variations and large-scale asymmetries in the dayglow.We demonstrate the model by using the three far ultraviolet cameras on the Imager for Magnetopause-to-Aurora Global Exploration(IMAGE)mission.The method can be applied to historical missions and is relevant for upcoming missions,such as the Solar wind Magnetosphere Ionosphere Link Explorer(SMILE)mission.

Surrogate modeling for unsaturated infiltration via the physics and equality-constrained artificial neural networks

Machine learning(ML)provides a new surrogate method for investigating groundwater flow dynamics in unsaturated soils.Traditional pure data-driven methods(e.g.deep neural network,DNN)can provide rapid predictions,but they do require sufficient on-site data for accurate training,and lack interpretability to the physical processes within the data.In this paper,we provide a physics and equalityconstrained artificial neural network(PECANN),to derive unsaturated infiltration solutions with a small amount of initial and boundary data.PECANN takes the physics-informed neural network(PINN)as a foundation,encodes the unsaturated infiltration physical laws(i.e.Richards equation,RE)into the loss function,and uses the augmented Lagrangian method to constrain the learning process of the solutions of RE by adding stronger penalty for the initial and boundary conditions.Four unsaturated infiltration cases are designed to test the training performance of PECANN,i.e.one-dimensional(1D)steady-state unsaturated infiltration,1D transient-state infiltration,two-dimensional(2D)transient-state infiltration,and 1D coupled unsaturated infiltration and deformation.The predicted results of PECANN are compared with the finite difference solutions or analytical solutions.The results indicate that PECANN can accurately capture the variations of pressure head during the unsaturated infiltration,and present higher precision and robustness than DNN and PINN.It is also revealed that PECANN can achieve the same accuracy as the finite difference method with fewer initial and boundary training data.Additionally,we investigate the effect of the hyperparameters of PECANN on solving RE problem.PECANN provides an effective tool for simulating unsaturated infiltration.

A methodology for damage evaluation of underground tunnels subjected to static loading using numerical modeling

We have proposed a methodology to assess the robustness of underground tunnels against potential failure.This involves developing vulnerability functions for various qualities of rock mass and static loading intensities.To account for these variations,we utilized a Monte Carlo Simulation(MCS)technique coupled with the finite difference code FLAC^(3D),to conduct two thousand seven hundred numerical simulations of a horseshoe tunnel located within a rock mass with different geological strength index system(GSIs)and subjected to different states of static loading.To quantify the severity of damage within the rock mass,we selected one stress-based(brittle shear ratio(BSR))and one strain-based failure criterion(plastic damage index(PDI)).Based on these criteria,we then developed fragility curves.Additionally,we used mathematical approximation techniques to produce vulnerability functions that relate the probabilities of various damage states to loading intensities for different quality classes of blocky rock mass.The results indicated that the fragility curves we obtained could accurately depict the evolution of the inner and outer shell damage around the tunnel.Therefore,we have provided engineers with a tool that can predict levels of damages associated with different failure mechanisms based on variations in rock mass quality and in situ stress state.Our method is a numerically developed,multi-variate approach that can aid engineers in making informed decisions about the robustness of underground tunnels.

Artificial intelligence-driven radiomics study in cancer:the role of feature engineering and modeling

Modern medicine is reliant on various medical imaging technologies for non-invasively observing patients’anatomy.However,the interpretation of medical images can be highly subjective and dependent on the expertise of clinicians.Moreover,some potentially useful quantitative information in medical images,especially that which is not visible to the naked eye,is often ignored during clinical practice.In contrast,radiomics performs high-throughput feature extraction from medical images,which enables quantitative analysis of medical images and prediction of various clinical endpoints.Studies have reported that radiomics exhibits promising performance in diagnosis and predicting treatment responses and prognosis,demonstrating its potential to be a non-invasive auxiliary tool for personalized medicine.However,radiomics remains in a developmental phase as numerous technical challenges have yet to be solved,especially in feature engineering and statistical modeling.In this review,we introduce the current utility of radiomics by summarizing research on its application in the diagnosis,prognosis,and prediction of treatment responses in patients with cancer.We focus on machine learning approaches,for feature extraction and selection during feature engineering and for imbalanced datasets and multi-modality fusion during statistical modeling.Furthermore,we introduce the stability,reproducibility,and interpretability of features,and the generalizability and interpretability of models.Finally,we offer possible solutions to current challenges in radiomics research.

First-principles study on stability and superconductivity of ternary hydride LaYH_(x)(x=2,3,6 and 8)

Recent studies have shown that the La-and Y-hydrides can exhibit significant superconducting properties under high pressures.In this paper,we investigate the stability,electronic and superconducting properties of LaYH_(x)(x=2,3,6 and 8)under 0-200 GPa.It is found that LaYH_(2) stabilizes in the C2/m phase at ambient pressure,and transforms to the Pmmn phase at 67 GPa.LaYH_(3) stabilizes in the C2/m phase at ambient pressure,and undergoes phase transitions of C2/m→P2_(1)/m→R3m at 12 GPa and 87 GPa,respectively.LaYH_(6) stabilizes in the P4_32_12 phase at ambient pressure,and undergoes phase transitions of P4_(3)2_(1)2→P4/mmm→Cmcm at 28 GPa and 79 GPa,respectively.LaYH_(8) stabilizes in the Imma phase at 60 GPa and transforms to the P4/mmm phase at 117 GPa.Calculations of the electronic band structures show that the P4/mmm-LaYH_(8) and all phases of LaYH_(2) and LaYH_(3) exhibit metallic character.For the metallic phases,we then study their superconducting properties.The calculated superconducting transition temperatures(T_c)are 0.47 K for C2/m-LaYH_(2) at 0 GPa,0 K for C2/m-LaYH_(3) at 0 GPa,and 55.51 K for P4/mmm-LaYH_(8) at 50 GPa.

The growth ofβphase in Mg-Gd-Y-Ni alloy by experimental and first-principles study

The paper reports on the atomic investigation aboutβphase in Mg_(96)Gd_(2)Y_(1)Ni_(1) alloy by using the first-principles study and the high-angle annular dark-field scanning transmission electron microscope(HAADF-STEM)corrected by atomic Cs.By using HAADF-STEM,the rectangularβphases were observed in the underage and peak aging stages in Mg_(96)Gd_(2)Y_(1)Ni_(1) alloy.Theβphase could be precipitated from the previously precipitatedβphase,and theβphase grew in steps when it was precipitated.A special transition structure of three atomic layer thicknesses was first observed at the edge of theβphase and the structure of this interface is probably as theβ/Mg_(1) interface for the analysis of thermodynamic characterization and electronic characterization.Theβ'phase and theβ_(H) structure were precipitated only at the edge of the length directions of theβphase.Theβ'phase continues to grow into aβphase directly without the formation ofβ_(1) phase,resulting in an increase in the length of theβphase,which is discovered for the first time.

Progressive fragmentation of granular assemblies within rockslides: Insights from discrete-continuous numerical modeling

Rock fragmentation plays a critical role in rock avalanches,yet conventional approaches such as classical granular flow models or the bonded particle model have limitations in accurately characterizing the progressive disintegration and kinematics of multi-deformable rock blocks during rockslides.The present study proposes a discrete-continuous numerical model,based on a cohesive zone model,to explicitly incorporate the progressive fragmentation and intricate interparticle interactions inherent in rockslides.Breakable rock granular assemblies are released along an inclined plane and flow onto a horizontal plane.The numerical scenarios are established to incorporate variations in slope angle,initial height,friction coefficient,and particle number.The evolutions of fragmentation,kinematic,runout and depositional characteristics are quantitatively analyzed and compared with experimental and field data.A positive linear relationship between the equivalent friction coefficient and the apparent friction coefficient is identified.In general,the granular mass predominantly exhibits characteristics of a dense granular flow,with the Savage number exhibiting a decreasing trend as the volume of mass increases.The process of particle breakage gradually occurs in a bottom-up manner,leading to a significant increase in the angular velocities of the rock blocks with increasing depth.The simulation results reproduce the field observations of inverse grading and source stratigraphy preservation in the deposit.We propose a disintegration index that incorporates factors such as drop height,rock mass volume,and rock strength.Our findings demonstrate a consistent linear relationship between this index and the fragmentation degree in all tested scenarios.

Discontinuity development patterns and the challenges for 3D discrete fracture network modeling on complicated exposed rock surfaces

Natural slopes usually display complicated exposed rock surfaces that are characterized by complex and substantial terrain undulation and ubiquitous undesirable phenomena such as vegetation cover and rockfalls.This study presents a systematic outcrop research of fracture pattern variations in a complicated rock slope,and the qualitative and quantitative study of the complex phenomena impact on threedimensional(3D)discrete fracture network(DFN)modeling.As the studies of the outcrop fracture pattern have been so far focused on local variations,thus,we put forward a statistical analysis of global variations.The entire outcrop is partitioned into several subzones,and the subzone-scale variability of fracture geometric properties is analyzed(including the orientation,the density,and the trace length).The results reveal significant variations in fracture characteristics(such as the concentrative degree,the average orientation,the density,and the trace length)among different subzones.Moreover,the density of fracture sets,which is approximately parallel to the slope surface,exhibits a notably higher value compared to other fracture sets across all subzones.To improve the accuracy of the DFN modeling,the effects of three common phenomena resulting from vegetation and rockfalls are qualitatively analyzed and the corresponding quantitative data processing solutions are proposed.Subsequently,the 3D fracture geometric parameters are determined for different areas of the high-steep rock slope in terms of the subzone dimensions.The results show significant variations in the same set of 3D fracture parameters across different regions with density differing by up to tenfold and mean trace length exhibiting differences of 3e4 times.The study results present precise geological structural information,improve modeling accuracy,and provide practical solutions for addressing complex outcrop issues.

The relationship between compartment models and their stochastic counterparts:A comparative study with examples of the COVID-19 epidemic modeling

Deterministic compartment models(CMs)and stochastic models,including stochastic CMs and agent-based models,are widely utilized in epidemic modeling.However,the relationship between CMs and their corresponding stochastic models is not well understood.The present study aimed to address this gap by conducting a comparative study using the susceptible,exposed,infectious,and recovered(SEIR)model and its extended CMs from the coronavirus disease 2019 modeling literature.We demonstrated the equivalence of the numerical solution of CMs using the Euler scheme and their stochastic counterparts through theoretical analysis and simulations.Based on this equivalence,we proposed an efficient model calibration method that could replicate the exact solution of CMs in the corresponding stochastic models through parameter adjustment.The advancement in calibration techniques enhanced the accuracy of stochastic modeling in capturing the dynamics of epidemics.However,it should be noted that discrete-time stochastic models cannot perfectly reproduce the exact solution of continuous-time CMs.Additionally,we proposed a new stochastic compartment and agent mixed model as an alternative to agent-based models for large-scale population simulations with a limited number of agents.This model offered a balance between computational efficiency and accuracy.The results of this research contributed to the comparison and unification of deterministic CMs and stochastic models in epidemic modeling.Furthermore,the results had implications for the development of hybrid models that integrated the strengths of both frameworks.Overall,the present study has provided valuable epidemic modeling techniques and their practical applications for understanding and controlling the spread of infectious diseases.

First-principles study on the diffusion behavior of Cs and I in Cr coating

Cs and I can migrate through fuel-cladding interfaces and accelerate the cladding corrosion process induced by the fuel-cladding chemical interaction.Cr coating has emerged as an important candidate for mitigating this chemical interaction.In this study,first-principles calculations were employed to investigate the diffusion behavior of Cs and I in the Cr bulk and grain boundaries to reveal the microscopic interaction mitigation mechanisms at the fuel-cladding interface.The interaction between these two fission products and the Cr coating were studied systematically,and the Cs and I temperature-dependent diffusion coefficients in Cr were obtained using Bocquet’s oversized solute-atom model and Le Claire’s nine-frequency model,respectively.The results showed that the Cs and I migration barriers were significantly lower than that of Cr,and the Cs and I diffusion coefficients were more than three orders of magnitude larger than the Cr self-diffusion coefficient within the temperature range of Generation-IV fast reactors(below 1000 K),demonstrating the strong penetration ability of Cs and I.Furthermore,Cs and I are more likely to diffuse along the grain boundary because of the generally low migration barriers,indicating that the grain boundary serves as a fast diffusion channel for Cs and I.

Hybrid modeling for carbon monoxide gas-phase catalytic coupling to synthesize dimethyl oxalate process

Ethylene glycol(EG)plays a pivotal role as a primary raw material in the polyester industry,and the syngas-to-EG route has become a significant technical route in production.The carbon monoxide(CO)gas-phase catalytic coupling to synthesize dimethyl oxalate(DMO)is a crucial process in the syngas-to-EG route,whereby the composition of the reactor outlet exerts influence on the ultimate quality of the EG product and the energy consumption during the subsequent separation process.However,measuring product quality in real time or establishing accurate dynamic mechanism models is challenging.To effectively model the DMO synthesis process,this study proposes a hybrid modeling strategy that integrates process mechanisms and data-driven approaches.The CO gas-phase catalytic coupling mechanism model is developed based on intrinsic kinetics and material balance,while a long short-term memory(LSTM)neural network is employed to predict the macroscopic reaction rate by leveraging temporal relationships derived from archived measurements.The proposed model is trained semi-supervised to accommodate limited-label data scenarios,leveraging historical data.By integrating these predictions with the mechanism model,the hybrid modeling approach provides reliable and interpretable forecasts of mass fractions.Empirical investigations unequivocally validate the superiority of the proposed hybrid modeling approach over conventional data-driven models(DDMs)and other hybrid modeling techniques.

Mechanical Properties and Electronic Structures of M(M=Ti,V,Cr,Mn and Fe)Dopedβ-Si_(3)N_(4) from First-Principle

The structures,mechanical properties and electronic structures of M metals(M=Ti,V,Cr,Mn and Fe)dopedβ-Si_(3)N_(4) were investigated by First-principles calculations within CASTEP.The calculated lattice parameters ofβ-Si_(3)N_(4) were consistent with previous date.The cohesive energy and formation enthalpy show that initialβ-Si_(3)N_(4) has the highest structural stability.The calculated elastic constant and the Voigt-Reuss-Hill approximation indicate that elastic moduli ofβ-Si_(3)N_(4) are slightly reduced by M doping.Based on Poisson’s and Pugh’s ratio,β-Si_(3)N_(4) is a ductile material and the toughness ofβ-Si_(3)N_(4) increases with M doping,and Fe doping exhibited the best toughness.The results of density of states,charge distributions and overlapping populations indicate thatβ-Si_(3)N_(4) has the strong covalent and ionic bond strength between N and Si.

Exploring Capillary Fringe Flow:Quasilinear Modeling with Kirchhoff Transforms and Gardner Model

Recent studies have underscored the significance of the capillary fringe in hydrological and biochemical processes.Moreover,its role in shallow waters is expected to be considerable.Traditionally,the study of groundwater flow has centered on unsaturated-saturated zones,often overlooking the impact of the capillary fringe.In this study,we introduce a steady-state two-dimensional model that integrates the capillary fringe into a 2-D numerical solution.Our novel approach employs the potential form of the Richards equation,facilitating the determination of boundaries,pressures,and velocities across different ground surface zones.We utilized a two-dimensional Freefem++finite element model to compute the stationary solution.The validation of the model was conducted using experimental data.We employed the OFAT(One_Factor-At-Time)method to identify the most sensitive soil parameters and understand how changes in these parameters may affect the behavior and water dynamics of the capillary fringe.The results emphasize the role of hydraulic conductivity as a key parameter influencing capillary fringe shape and dynamics.Velocity values within the capillary fringe suggest the prevalence of horizontal flow.By variation of the water table level and the incoming flow q0,we have shown the correlation between water table elevation and the upper limit of the capillary fringe.

Total ionizing dose effect modeling method for CMOS digital-integrated circuit

Simulating the total ionizing dose(TID)of an electrical system using transistor-level models can be difficult and expensive,particularly for digital-integrated circuits(ICs).In this study,a method for modeling TID effects in complementary metaloxide semiconductor(CMOS)digital ICs based on the input/output buffer information specification(IBIS)was proposed.The digital IC was first divided into three parts based on its internal structure:the input buffer,output buffer,and functional area.Each of these three parts was separately modeled.Using the IBIS model,the transistor V-I characteristic curves of the buffers were processed,and the physical parameters were extracted and modeled using VHDL-AMS.In the functional area,logic functions were modeled in VHDL according to the data sheet.A golden digital IC model was developed by combining the input buffer,output buffer,and functional area models.Furthermore,the golden ratio was reconstructed based on TID experimental data,enabling the assessment of TID effects on the threshold voltage,carrier mobility,and time series of the digital IC.TID experiments were conducted using a CMOS non-inverting multiplexer,NC7SZ157,and the results were compared with the simulation results,which showed that the relative errors were less than 2%at each dose point.This confirms the practicality and accuracy of the proposed modeling method.The TID effect model for digital ICs developed using this modeling technique includes both the logical function of the IC and changes in electrical properties and functional degradation impacted by TID,which has potential applications in the design of radiation-hardening tolerance in digital ICs.

Data-Driven Earthquake Multi-impact Modeling:A Comparison of Models

In this study,a broad range of supervised machine learning and parametric statistical,geospatial,and non-geospatial models were applied to model both aggregated observed impact estimate data and satellite image-derived geolocated building damage data for earthquakes,via regression-and classification-based models,respectively.For the aggregated observational data,models were ranked via predictive performance of mortality,population displacement,building damage,and building destruction for 375 observations across 161 earthquakes in 61 countries.For the satellite image-derived data,models were ranked via classification performance(damaged/unaff ected)of 369,813 geolocated buildings for 26 earthquakes in 15 countries.Grouped k-fold,3-repeat cross validation was used to ensure out-of-sample predictive performance.Feature importance of several variables used as proxies for vulnerability to disasters indicates covariate utility.The 2023 Türkiye-Syria earthquake event was used to explore model limitations for extreme events.However,applying the AdaBoost model on the 27,032 held-out buildings of the 2023 Türkiye-Syria earthquake event,predictions had an AUC of 0.93.Therefore,without any geospatial,building-specific,or direct satellite image information,this model accurately classified building damage,with significantly improved performance over satellite image trained models found in the literature.

Modeling and analysis of hysteresis using the Maxwell-slip model for variable stiffness actuators

Hysteresis non-linearity in variable stiffness actuators(VSAs)causes significant torque errors and reduces the stability of the actuators,leading to poor human–computer interaction performance.At present,fewer hysteresis compensation models have been developed for compliant drives,so it is necessary to establish a suitable hysteresis model for compliant actuators.In this work,a new model with a combination of the Maxwell-slip model and virtual deformation is proposed and applied to an elbow compliant actuator.The method divides the periodic variation of the actuator into three parts:an ascending phase,a descending phase,and a transition phase.Based on the concept of virtual deformation,the nonlinear hysteresis curve is transformed into a polyline,and the output torque is estimated using the revised Maxwell-slip model.The simulation results are compared with the experimental data.Its torque error is controlled within 0.2Nm,which validates the model.An inverse model is finally established to calculate the deformation deflection angle for hysteresis compensation.The results show that the inverse model has high accuracy,and the deformation deflection is less than 0.15 rad.

Modularized and Parametric Modeling Technology for Finite Element Simulations of Underground Engineering under Complicated Geological Conditions

The surrounding geological conditions and supporting structures of underground engineering are often updated during construction,and these updates require repeated numerical modeling.To improve the numerical modeling efficiency of underground engineering,a modularized and parametric modeling cloud server is developed by using Python codes.The basic framework of the cloud server is as follows:input the modeling parameters into the web platform,implement Rhino software and FLAC3D software to model and run simulations in the cloud server,and return the simulation results to the web platform.The modeling program can automatically generate instructions that can run the modeling process in Rhino based on the input modeling parameters.The main modules of the modeling program include modeling the 3D geological structures,the underground engineering structures,and the supporting structures as well as meshing the geometric models.In particular,various cross-sections of underground caverns are crafted as parametricmodules in themodeling program.Themodularized and parametric modeling program is used for a finite element simulation of the underground powerhouse of the Shuangjiangkou Hydropower Station.This complicatedmodel is rapidly generated for the simulation,and the simulation results are reasonable.Thus,this modularized and parametric modeling program is applicable for three-dimensional finite element simulations and analyses.