Computers and Geotechnics最新文献

{"title":"Influence of curtain parameter spatial variability on dam foundation seepage field considering regional heterogeneity","authors":"Lifeng Wen,&nbsp;Nannan He,&nbsp;Yanlong Li,&nbsp;Junrui Chai","doi":"10.1016/j.compgeo.2025.107435","DOIUrl":"10.1016/j.compgeo.2025.107435","url":null,"abstract":"<div><div>Seepage is a critical factor affecting safety and stability of dam foundations in hydraulic engineering. Anti-seepage curtain is an important measure to control seepage field in dam foundation. Due to the influence of construction quality and operating conditions, the parameters of the anti-seepage curtain inevitably exhibit spatial variability. The parameters spatial variability directly affects seepage field distribution and seepage control effectiveness in the dam foundation. This paper investigates influence of the parameters spatial variability of anti-seepage curtain on the seepage field of dam foundation. Firstly, a non-stationary random field discretization method considering regional heterogeneity in characteristics of the anti-seepage curtain is proposed by combining decomposed Karhunen-Loève (K-L) expansion and Sequential Gaussian Simulation (SGSIM) method. In order to consider regional heterogeneity in characteristics of the anti-seepage curtain that caused by the actual construction quality, the anti-seepage curtain is divided into three regions according to permeability coefficients difference obtained by borehole data. The decomposed K-L expansion method is employed to generate parameters random fields considering the regional heterogeneity by defining parameters mean in different regions. The SGSIM method is used to modify the parameters statistical characteristics of the borehole area to match the local parameters of borehole area with borehole data. Then, the influence of the parameters spatial variability in the anti-seepage curtain and stochastic field parameters on the seepage flow rate and hydraulic gradient mean of the dam foundation is discussed. The results indicate that the proposed parameters spatial variability simulation method could effectively characterize the non-stationary random field characteristics of the anti-seepage curtain permeability parameters. The spatial variability of the permeability parameters would not significantly change the distribution pattern of the seepage field in the dam foundation, but it has a significant impact on the statistical characteristics of the seepage field evaluation indexes, especially in the area near the anti-seepage curtain.</div></div>","PeriodicalId":55217,"journal":{"name":"Computers and Geotechnics","volume":"186 ","pages":"Article 107435"},"PeriodicalIF":5.3,"publicationDate":"2025-06-21","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144330615","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":1,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Viscoelastoplastic model for an integrated tunnel-track system in weak rock formations","authors":"Hafsa Farooq,&nbsp;Sanjay Nimbalkar","doi":"10.1016/j.compgeo.2025.107433","DOIUrl":"10.1016/j.compgeo.2025.107433","url":null,"abstract":"<div><div>Underground railway transportation systems offer a sustainable solution for space constraints, especially in urban areas, however, the construction of such systems embedded in weak rock presents significant challenges. This study introduces a novel computational approach to evaluate the dynamic response of integrated tunnel-track systems while considering the effect of moving trains. The effect of a moving train load is assessed for the tunnel inlay and the weak rock surrounding the tunnel. A viscoelastoplastic model, incorporating Hooke’s, Newton’s and St. Venant’s elements, is introduced to simulate complex interaction within the integrated tunnel-track system. The model utilises Hoek-Brown criteria for the rock mass surrounding the tunnel. Firstly, the dynamic effects of the moving train loads are assessed for the track and the tunnel inlay, and then the resulting stresses on the tunnel inlay are treated as internal support pressure to evaluate the response of the surrounding karstic rock. The model is validated against existing experimental and numerical investigations. This study provides a valuable computational tool for practising engineers to assess the long-term structural response of integrated tunnel-track systems constructed in challenging weak rock formations.</div></div>","PeriodicalId":55217,"journal":{"name":"Computers and Geotechnics","volume":"186 ","pages":"Article 107433"},"PeriodicalIF":5.3,"publicationDate":"2025-06-21","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144330701","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":1,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Micro-mechanism of shield-soil interaction during shield tunneling based on DEM-MBD coupling simulation","authors":"Hui Jin ,&nbsp;Enzhi Wang ,&nbsp;Xiaoli Liu ,&nbsp;Zhen Xu ,&nbsp;Mingyang Wang ,&nbsp;Dajun Yuan","doi":"10.1016/j.compgeo.2025.107439","DOIUrl":"10.1016/j.compgeo.2025.107439","url":null,"abstract":"<div><div>Shield interacts with soil based on its movement and dynamics, along with soil displacement and stress. Poor interaction management can harm construction quality and cause safety risks in strata and structures. This paper conducts coupling simulations of the discrete element method (DEM) and the multibody dynamics method (MBD) to facilitate the shield tunneling in soil. This method uncovers the interaction mechanism between the shield and soil at a mesoscopic scale. Results reveal that the particle stress and velocity fields are most prominent in front of the cutterhead during excavation. These fields expand forward from the cutterhead in a flared shape. Consequently, earth pressure and contact force on the cutterhead are significantly greater than those on the shell. Notably, high-stress zones form near the scrapers and center knives due to their penetration and cutting action on the soil. The cutterhead experiences the highest longitudinal load, followed by the front shield and then the tail shield, with a load ratio of approximately 7:2:1. The soil filling rate of the chamber and the shield attitude greatly influence the shield-soil interaction. Increasing the rotation speed of the screw conveyor decreases the soil filling rate of the chamber, alters particle flow direction, and changes the force transfer pattern within the chamber. This, in turn, reduces the total thrust and cutterhead torque. Conversely, increasing the propulsion speed of the bottom cylinders enhances the cutterhead penetration and the contact forces at the shield periphery, resulting in an increased total thrust and torque. Importantly, the heavy shield tends to move downward to gain ground resistance. Adjusting its posture alters the vertical bearing mode from lower-side to upper-side bearing, which can minimize total thrust and torque. Additionally, torsion of the shell observed during tunnel construction negatively impacts posture adjustment. Relevant results can guide high-quality shield construction.</div></div>","PeriodicalId":55217,"journal":{"name":"Computers and Geotechnics","volume":"186 ","pages":"Article 107439"},"PeriodicalIF":5.3,"publicationDate":"2025-06-21","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144330616","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":1,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Effect of bedding structures on hydraulic fracture propagation behavior investigated using a coupled thermo-hydraulic-mechanical numerical model based on the phase-field cohesive zone method","authors":"Xun Gong ,&nbsp;Zhijun Jin ,&nbsp;Xinhua Ma ,&nbsp;Yuyang Liu ,&nbsp;Chunyang Hong","doi":"10.1016/j.compgeo.2025.107427","DOIUrl":"10.1016/j.compgeo.2025.107427","url":null,"abstract":"<div><div>Shale contains numerous bedding structures, and the degree of communication with hydraulic fractures during the fracturing process affects the complexity of the generated fracture network. However, the propagation behavior of hydraulic fractures in bedding is complex and variable, and there are many influencing factors, limiting the effect of reservoir stimulation. Therefore, to improve the reservoir stimulation volume, this paper establishes a two-dimensional thermo-hydraulic-mechanical-phase field damage model based on the phase-field method and the cohesive zone method to simulate the propagation behavior of hydraulic fractures in bedding. The results show that three propagation patterns, namely, turning, turning and deviating and directly passing through, may occur when a hydraulic fracture encounters bedding. The specific conditions are affected by the mechanical properties of the rock, stress, bedding inclination angle, temperature, and other factors. During the fracturing process, increases in the ratio of the elastic modulus of the rock matrix to that of bedding, the bedding inclination angle, and the temperature are conducive to hydraulic fracture turning along bedding, and the bedding intersection angle is 0°, which is the optimal angle for hydraulic fractures to pass through bedding directly. In contrast, an increase in the principal stress difference is conducive to hydraulic fractures passing through bedding directly. Second, on the basis of the four influencing factors, namely, the shale matrix/bedding elastic modulus ratio, bedding inclination angle, principal stress difference, and temperature, a propagation criterion for hydraulic fractures in bedding under different influencing factors was established. The research results will guide the design and optimization of fracturing in shale with bedding.</div></div>","PeriodicalId":55217,"journal":{"name":"Computers and Geotechnics","volume":"186 ","pages":"Article 107427"},"PeriodicalIF":5.3,"publicationDate":"2025-06-21","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144330704","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":1,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Long-term deformation mechanism of reservoir bank slope using multi-objective inversion: A case study at Duonuo reservoir","authors":"Shu Jiang,&nbsp;Yingchao Gao,&nbsp;Wei Wei,&nbsp;Qinghui Jiang","doi":"10.1016/j.compgeo.2025.107423","DOIUrl":"10.1016/j.compgeo.2025.107423","url":null,"abstract":"<div><div>Modeling the long-term deformation of reservoir slope remains a challenging issue due to the difficulties in properly calibrating the creep properties for fractured rocks. Based on the deformation phenomenon shown by field measured data, this study adopts the modified Nishihara model to simulate the long-term deformation of the reservoir slope. Then, a combined procedure of orthogonal design (OD), forward numerical calculation, artificial neural network (ANN), and multi-objective genetic algorithm is adopted for the inverse analysis with high computational efficiency to estimate creep parameters of rock mass. The numerical results based on the inversion parameters well reproduce the long-term deformation process of the reservoir slope. The deformation of the reservoir slope is mainly triggered by rock mass in the water-level fluctuation zone (WLFZ) due to the deterioration effect of water-level fluctuation. During the filling stage, the rock masses in WLFZ are gradually saturated and softened, which leads to viscoelastic deformation of the slope. While, as the number of filling-drawdown cycles increases, the deterioration effect of water on the viscoelastic component of rock mass gradually decreases. During the drawdown stage, due to the delayed dissipation of pore water pressure within the slope, the plastic yields of rock masses in WLFZ occur and induce viscoplastic deformation. This study provides an efficient and effective approach to obtain an optimal solution in large-scale inverse modeling of long-term deformation of slope, further improving the safety assessment.</div></div>","PeriodicalId":55217,"journal":{"name":"Computers and Geotechnics","volume":"186 ","pages":"Article 107423"},"PeriodicalIF":5.3,"publicationDate":"2025-06-19","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144312891","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":1,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Three-dimensional discrete element simulation of a transversely isotropic rock considering initial cracks","authors":"Kaihui Li ,&nbsp;Jiezhen Chen ,&nbsp;Dongya Han ,&nbsp;Zhizhen Liu ,&nbsp;Kang Peng ,&nbsp;Jiangteng Li","doi":"10.1016/j.compgeo.2025.107434","DOIUrl":"10.1016/j.compgeo.2025.107434","url":null,"abstract":"<div><div>Shale, as a reservoir rock, exhibits obvious anisotropy in both its physical and mechanical properties. An in-depth understanding of its mechanical behaviors considering the pre-existing cracks is vital for the shale gas extraction. In this study, a three-dimensional discrete element model of transversely isotropic rock considering initial cracks was developed to investigate the nonlinear crack-closure behavior and failure characteristics of shale under uniaxial compression. The results show that the nonlinear crack-closure and failure behaviors of shale are greatly affected by the anisotropic angle (<em>β</em>), microcrack density and width in rock matrix. The crack-induced strain of shale specimen is mainly dependent on the change of microcrack width in soft rock matrix, which increases with <em>β</em>. Compared with the specimen without considering initial cracks, regardless of <em>β</em>, the specimen considering initial cracks is characterized by a lower crack-initiation stress and a higher brittleness index, which is consistent with the evolution of microcracks formed in specimen. The failure mechanism of specimen is further revealed by the stress transfer between stiff and soft rock matrices. With increasing <em>β</em>, the axial stress in soft rock matrix increases, whereas that in stiff rock matrix decreases first and then increases, resulting in the transition of failure mode of specimen. Fully considering the influences of weak planes and initial cracks, the proposed model can well capture the variation of rock brittleness and predict the critical breakout pressure of boreholes, thus providing a more reliable basis for analyzing engineering problems under complex stress conditions in shale gas extraction.</div></div>","PeriodicalId":55217,"journal":{"name":"Computers and Geotechnics","volume":"186 ","pages":"Article 107434"},"PeriodicalIF":5.3,"publicationDate":"2025-06-19","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144312888","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":1,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Regional slope stability and slope failure mechanics under strong seismic events in the upper reaches of Yellow River cascade reservoir dams via the SEM–FDEM–SPH coupling approach","authors":"Benbo Sun ,&nbsp;Shaowei Hu ,&nbsp;Juan Wang ,&nbsp;Jinjun Guo ,&nbsp;Shuyu Wu ,&nbsp;Pengfei Wang","doi":"10.1016/j.compgeo.2025.107334","DOIUrl":"10.1016/j.compgeo.2025.107334","url":null,"abstract":"<div><div>Early warning and hazardous point monitoring for large-scale slopes under strong ground motions are pivotal tools for safeguarding assets and stakeholders. In addition, the local-scale failure mechanisms and processes for earthquake-induced landslides can provide detailed information for decision making and planning. Providing regional seismic wave fields is critical for assessing regional slope stability and local-scale earthquake-induced landslides. Currently, most stable assessments and mechanism analyses of landslides use specified pseudo-accelerations or a small number of earthquake records because earthquake records are often very limited. Notably, these approaches cannot fully consider the spatial differences in seismology characteristics at various locations. In other words, there is always a generation gap between earthquake action and slope response and stability assessment. To reduce this gap, a hybrid simulation framework is created via the spectral element method (SEM), Bishop method, finite-discrete element method (FDEM), and smoothed-particle hydrodynamics (SPH) for large- and local-scale slopes in the upper reaches of Yellow River cascade reservoir dams. The SEM and Bishop methods enable straightforward estimation of seismic stability for regional slope systems subjected to earthquakes across 19 different scenarios by considering the sources, propagation paths, regional topographies, and geological conditions of the existing uncertainties. The weak coupling strategy of the time domain reduction method is implemented to achieve local-scale failure processes of earthquake-induced landslides via SEM–FDEM–SPH. The results indicate that the multiscale framework can realize regional earthquake landslide stability assessment and the chain disaster simulation of local unstable slopes. For projects involving regional-scale slope seismic stability assessment and local-scale mechanistic analysis, a physics-based mechanism framework for simulation studies and risk prediction is recommended.</div></div>","PeriodicalId":55217,"journal":{"name":"Computers and Geotechnics","volume":"186 ","pages":"Article 107334"},"PeriodicalIF":5.3,"publicationDate":"2025-06-19","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144312890","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":1,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Quantifying the contributions of multi-source uncertainties in dynamic soil-pile-nuclear power plant structure response via probability density evolution method","authors":"Zhenning Ba ,&nbsp;Chenyang Kuo ,&nbsp;Jianwen Liang","doi":"10.1016/j.compgeo.2025.107430","DOIUrl":"10.1016/j.compgeo.2025.107430","url":null,"abstract":"<div><div>Soil–pile–nuclear power plant (SPNPP) structures on non-bedrock sites exhibit complex seismic behaviors due to inherent uncertainties in ground motion, geotechnical properties, and structural materials. However, comprehensive uncertainty quantification remains limited. This work employs the Probability Density Evolution Method (PDEM), integrated with the Change of Probability Measure and Fréchet-derivative-based sensitivity analysis, to systematically assess the stochastic dynamic response of SPNPP structures. Taking the HPR1000 nuclear power plant as an engineering case study, a refined three-dimensional (3D) finite element model of the SPNPP system is established. The effects of 25 uncertain parameters, covering three uncertain sources (ground motion, geotechnical materials, and structural materials), are systematically analyzed. Based on 900 dynamic time-history simulations, the relative influence of each parameter is quantified. Results indicate that when all uncertainty sources are considered simultaneously, the coefficient of variation (COV) of the SPNPP structural response increases significantly. Ground motion uncertainty contributes most substantially to the overall structural response variation. For the top drift ratio response of the containment structure, uncertainties in soil shear wave velocity and structural elastic modulus both play key roles. Global sensitivity analysis (GSA) results for the 25 uncertain parameters are also presented. This study provides a practical framework for uncertainty quantification and performance-based optimization in seismic design of nuclear power structures on non-bedrock sites.</div></div>","PeriodicalId":55217,"journal":{"name":"Computers and Geotechnics","volume":"186 ","pages":"Article 107430"},"PeriodicalIF":5.3,"publicationDate":"2025-06-18","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144306972","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":1,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Nonlocal gradient viscoplastic analyses of interacting rate and scale effects in borehole failure","authors":"Dawei Xue ,&nbsp;Xilin Lü ,&nbsp;Maosong Huang ,&nbsp;Keng-Wit Lim","doi":"10.1016/j.compgeo.2025.107425","DOIUrl":"10.1016/j.compgeo.2025.107425","url":null,"abstract":"<div><div>While rate and scale effects influence the strength and inelastic behavior of geomaterials, their interplay in borehole failure problems remains largely unexplored, posing challenges for assessing the stability of small-diameter boreholes in time-sensitive rocks. This work proposes a nonlocal gradient viscoplastic model to numerically investigate this issue. Within the model, the viscoplastic shear strain is linked to a nonlocal variable through a non-standard balance equation, whose analogy to the heat equation enables a straightforward numerical implementation of the coupled governing equations. At the integration point level, the classical fully implicit integration scheme is extended to incorporate both rate dependence and nonlocal gradient regularization. The improved applicability of this framework over standard approaches is demonstrated through numerical examples. Simulations of loading-induced borehole failure, validated against experimental data, reveal that hole strength decreases as the system scale to material length scale ratio increases, indicating pronounced scale effects. Meanwhile, the failure pattern changes from a homogenous type to shear-bands-dominated one. Higher loading rates increase hole strength and attenuate scale effects, as evidenced by reduced discrepancies in strength and failure modes across different borehole sizes. Under creep conditions, smaller boreholes require higher far-field pressures to induce delayed failure. When subjected to identical, sufficiently high far-field pressures, larger boreholes fail in shorter creep times, while increasing the pressure mitigates these differences.</div></div>","PeriodicalId":55217,"journal":{"name":"Computers and Geotechnics","volume":"186 ","pages":"Article 107425"},"PeriodicalIF":5.3,"publicationDate":"2025-06-18","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144306971","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":1,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"A finite-strain procedure for calculating fictitious support pressure of a tunnel based on a re-developed GRC and LDP considering large deformation","authors":"Kai Guan ,&nbsp;Hongping Li ,&nbsp;Wancheng Zhu ,&nbsp;Xiongwei Zhang","doi":"10.1016/j.compgeo.2025.107426","DOIUrl":"10.1016/j.compgeo.2025.107426","url":null,"abstract":"<div><div>Fictitious support pressure on tunnel walls characterizes the face effect as the tunnel advances, necessitating quantification to fully mobilize self-supporting capacity, especially in rock masses undergoing large deformation. This study presents a method for calculating fictitious support pressure along a tunnel profile using a re-developed Ground Reaction Curve (GRC) and Longitudinal Deformation Profile (LDP). It employs a kinematic additive decomposition of the deformation rate tensor based on finite strain and hypoelastic–plastic theory, accounting for the motion of material points. The numerical implementation is provided, with examples for validation demonstrating the accuracy of the re-developed GRC and LDP in scenarios involving both slight deformations and significant large strains. Extensive efforts address practical issues affecting fictitious support pressure, including deformation formulation for large strain analysis, failure criteria, ground conditions, and rock behaviors governed by constitutive relationships and dilatancy. Results indicate that despite displacement deviations from small strain analysis and a simple linear failure criterion, these methods offer reliable approximations for predicting fictitious support pressure and tunnel convergence in suitable rock quality, thus eliminating the need for rigorous large strain elasto-plastic analyses and complex nonlinear failure criteria during the preliminary design stage of tunneling. An application case study is given to illustrate the practical use of the proposed finite-strain Confinement-Convergence Method framework for analyzing the support effect of steel sets during tunnel advancement.</div></div>","PeriodicalId":55217,"journal":{"name":"Computers and Geotechnics","volume":"186 ","pages":"Article 107426"},"PeriodicalIF":5.3,"publicationDate":"2025-06-18","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144306970","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":1,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}