International Journal for Numerical and Analytical Methods in Geomechanics最新文献

{"title":"Effect of Mineral Heterogeneity on the Triaxial Compressive Behavior of Granular Rock: A PFC3D‐GBM Numerical Study","authors":"Zi‐ming Song, Qing‐bin Meng, Hai Pu, Jiang‐yu Wu, Jiang‐feng Liu","doi":"10.1002/nag.70077","DOIUrl":"https://doi.org/10.1002/nag.70077","url":null,"abstract":"The mechanical properties of granular rocks are intrinsically linked to the heterogeneity of mineral grains. In this study, a 3D grain‐based model (PFC3D‐GBM) incorporating a hybrid contact model was established. Triaxial compression numerical simulations were conducted to investigate how confining pressure and mineral heterogeneity—characterized by spatial distribution, volume fraction, geometric size, and boundary strength of mineral grains—affect the fracture behavior and mechanical properties of coarse‐grained limestone. The results show that force chains serve as load‐bearing paths that distribute the external load. The orientation of these chains remains uniform, irrespective of the load magnitude. The ratio of high‐strength force chains to weak‐strength force chains (<jats:italic>R</jats:italic><jats:sub>HF/WF</jats:sub>) and the ratio of high‐strength force chains to microcracks (<jats:italic>R</jats:italic><jats:sub>HF/c</jats:sub>) serve as quantitative indicators of bearing capacity and fracture resistance. At lower confining pressure (<jats:italic>σ</jats:italic><jats:sub>3</jats:sub> ≤ 10 MPa), the specimens exhibit reduced peak strength and increased brittleness, with failure typically occurring along a single macroscopic fracture plane. As <jats:italic>σ</jats:italic><jats:sub>3</jats:sub> increases, the formation of “X”‐shaped failure becomes more pronounced. At the grain scale, a higher dolomite volume facilitates the formation of additional high‐strength force chains within the specimen, along with an increase in Tran‐D‐cs at peak stress. As Tran‐D‐cs require a higher concentration of stress, <jats:italic>σ</jats:italic><jats:sub>p</jats:sub> increases by 13%. Moreover, stronger grain boundary strength enhances the stress concentration needed for fracture by modifying the proportion of Trans‐c, leading to a 70 MPa increase in <jats:italic>σ</jats:italic><jats:sub>p</jats:sub>. Larger grain sizes, in turn, accommodate more force chains, necessitating a higher <jats:italic>σ</jats:italic><jats:sub>p</jats:sub> for fracture.","PeriodicalId":13786,"journal":{"name":"International Journal for Numerical and Analytical Methods in Geomechanics","volume":"52 1","pages":""},"PeriodicalIF":4.0,"publicationDate":"2025-09-16","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145072035","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Control Mechanism of Large Deformation in Soft Rock Tunnel Based on Polyurethane Foam Buffer Layer","authors":"Jimeng Feng, Junfu Fu, Jiadai Song, Junru Zhang, Longyan Duan, Bo Wang, Zhiqiang Sun","doi":"10.1002/nag.70072","DOIUrl":"https://doi.org/10.1002/nag.70072","url":null,"abstract":"Large deformations of the surrounding rock often occur in tunnels excavated in soft rocks with high geostress, which can lead to the destruction of the lining structure. Therefore, the adoption of an appropriate support structure to control soft rock deformation is particularly important to ensure the safety of the tunnel in the preconstruction and post‐operation. In this paper, a yielding support system containing polyurethane (PU) foam energy‐absorbing buffer layer is proposed, and based on the yielding support mechanism, combined with numerical simulation and model experiments, the thickness of the energy‐absorbing buffer layer and the location of the arrangement are specifically studied. The findings indicate that the yielding layer, with a PU foam energy‐absorbing buffer layer, exhibits a substantial yielding phase compared to a rigid support structure. The strain‐dependent adjustment of the elastic modulus method can better simulate the stress characteristics of the yielding layer in numerical calculations than the case of a constant modulus of elasticity. In this study, the optimum thickness of the PU foam energy‐absorbing buffer layer was 20 cm, and the optimum location was between the initial support structure and the secondary lining structure. Scale model testing further corroborates these results, demonstrating that the rigid support system's secondary lining structure develops cracks at the arched roof and elevated arch locations under a 2.0‐MPa load. In contrast, the yielding support system absorbs energy from surrounding rock deformation through compression, maintaining stability even under an 8.0‐MPa load, thereby enhancing the overall performance of the support structure system. The scale model test analysis serves as additional validation for the numerical simulation results.","PeriodicalId":13786,"journal":{"name":"International Journal for Numerical and Analytical Methods in Geomechanics","volume":"36 1","pages":""},"PeriodicalIF":4.0,"publicationDate":"2025-09-15","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145072037","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Analytical Solutions for Evaluating the Efficacy of the Clay Shock Method in Shield Tunnels With Space Curves","authors":"Jiannan Xie, Pengfei Li, Zhaoguo Ge, Shuang Chen, Shaohua Li, Fei Jia","doi":"10.1002/nag.70076","DOIUrl":"https://doi.org/10.1002/nag.70076","url":null,"abstract":"The Clay Shock method has been widely applied in shield tunneling engineering, attributed to its strong adaptability and prominent characteristic of rapid solidification for stratum reinforcement. As a supplementary measure to conventional grouting and filling techniques, it can further compensate for stratum losses. However, the impact of employing Clay Shock on the stratum to which the tunnel belongs and its efficacy in stratum control remain elusive. This paper presents a theoretical analysis method for evaluating the efficacy of the Clay Shock method in shield tunnels with space curves. The analytical method includes two shield cross‐section models. By applying the elasticity general solution, a comprehensive analytical solution is derived in this study to quantify the variations in the ground displacement field and additional stress field induced by shield tunnel boring after implementing the Clay Shock method. After conducting an analysis of the relevant parameters in accordance with the specific project, the method is validated through numerical methods and field data. It was found that the relationship between the grouting method on the outside of the pipe sheet and the mechanical disturbance of the formation is more complex. It is inaccurate to use only the simple formation loss rate to represent the formation deformation in the unit case. We propose a corresponding grouting design for the above problems. This study not only elucidates the specific impacts of the Clay Shock method on soil behavior, but also establishes a comprehensive theoretical foundation for its application principle in shield tunnelling.","PeriodicalId":13786,"journal":{"name":"International Journal for Numerical and Analytical Methods in Geomechanics","volume":"171 1","pages":""},"PeriodicalIF":4.0,"publicationDate":"2025-09-15","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145072036","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Assessment of Heat Transfer‐Induced Contaminant Transport for a GMB/GCL/CCL Liner System","authors":"Jin‐Wei Qiu, Jia‐Jun Pan, Ding‐Bao Song, Bo Hu, Jun Tong, Xian Zhou","doi":"10.1002/nag.70071","DOIUrl":"https://doi.org/10.1002/nag.70071","url":null,"abstract":"This paper presents a numerical investigation for full transient coupled heat transfer and contaminant transport in a landfill bottom liner system consisting of a geomembrane (GMB), a geosynthetic clay liner (GCL), and a compacted clay liner (CCL), considering the effect of three‐phases in saturated soil (i.e., solid particles, mobile pore fluid, and immobile pore fluid). Effects of effective porosity, GMB defects and wrinkles, and temperature‐induced changes in effective molecular diffusion coefficient, permeability coefficient, and distribution coefficient are considered. The proposed numerical solution is verified by comparing against the available analytical solutions, experimental results, and numerical model. Using the verified solution, simulations are performed to evaluate the effects of the temperature at the top boundary, GMB defects, GCL permeability, effective porosity of the GCL and CCL, and CCL thickness on contaminant transport. All the simulation results indicate that, regardless of the range of parameter values, heat transfer has a significant impact on contaminant transport through the GMB/GCL/CCL liner system, not only during the period of heat transfer but also decades after heat transfer is complete. For the condition considered in this study, neglecting the effect of heat transfer can underestimate the cumulative contaminant mass outflow by at least 1.487 times and, therefore, may lead to unconservative liner design and adverse environmental impacts.","PeriodicalId":13786,"journal":{"name":"International Journal for Numerical and Analytical Methods in Geomechanics","volume":"104 1","pages":""},"PeriodicalIF":4.0,"publicationDate":"2025-09-15","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145072168","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Transient State Analysis of Rainfall Infiltration Into Layered Vegetated Soils","authors":"Cheng Yuan, Changbing Qin, Liang Li, Xiaoqin Lei","doi":"10.1002/nag.70068","DOIUrl":"https://doi.org/10.1002/nag.70068","url":null,"abstract":"Vegetation plays a pivotal role in altering the hydraulic properties of soils, and its effect is stemmed from plant transpiration and also undermined by rainfall. As for the intricate transient behavior of vegetated soils under rainfall, many extant studies either presumed homogeneous soil properties or concentrated on a steady‐state analysis, neglecting the moisture distribution with time in layered vegetated soils. This study proposes a transient analytical model for rainfall infiltration in layered vegetated soils based on Richard's equation. This model is developed by incorporating a depth‐dependent root water uptake sink term and utilizing variable substitution, separation of variables, and the transfer matrix method to derive the pore water pressure (PWP) distribution across different soil layers. The results calculated from the proposed model coincide with those from numerical simulations, with a maximum error of 2.5% in negative PWP in the case study, demonstrating the robustness of this model. Moreover, a parametric study is performed for a better understanding of the impacts of rainfall intensity, thickness of vegetated soil, and transpiration rate on the hydraulic properties of layered vegetated soils. This study well addresses the gap in transient analytical solutions of rainfall infiltration into layered vegetated soils and provides a benchmark for related numerical models.","PeriodicalId":13786,"journal":{"name":"International Journal for Numerical and Analytical Methods in Geomechanics","volume":"29 1","pages":""},"PeriodicalIF":4.0,"publicationDate":"2025-09-12","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145072039","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Data‐Driven Detection of Internal Erosion Initiation in Gap‐Graded Soils: Combining Particle‐Scale CFD‐DEM Simulation With 3D Convolutional Autoencoder","authors":"Jie Qi, Negin Yousefpour, Guillermo A. Narsilio, Mehdi Pouragha","doi":"10.1002/nag.70063","DOIUrl":"https://doi.org/10.1002/nag.70063","url":null,"abstract":"Internal erosion in gap‐graded soils poses significant risks to water‐retaining structures such as earth dams. However, its underlying mechanisms at the particle scale remain poorly understood. This study couples the discrete element method (DEM) with computational fluid dynamics (CFD) to simulate internal erosion in gap‐graded soil assemblies and employs data‐driven techniques to detect early‐stage erosion. Particle‐scale parameters, such as contact forces, particle velocity and fluid velocity, are extracted from the transient CFD‐DEM simulations. These features are transformed into multi‐dimensional voxel‐based tensors representing the particle–fluid interactions, which are used to train deep learning models. Autoencoder models with 3D convolutional neural network (CNN) layers as encoder and decoder are developed to investigate the micro‐scale patterns within the particle‐fluid assembly. Through sequential training techniques, the temporal evolution of anomalies is captured, enabling identification of the initiation point of internal erosion. The results reveal how microscale behaviours, such as particle motion, contact forces and contact number, contribute to macroscale erosion processes. The outcome of this research can inspire further research into AI‐based early detection techniques for internal erosion in earth dams.","PeriodicalId":13786,"journal":{"name":"International Journal for Numerical and Analytical Methods in Geomechanics","volume":"6 1","pages":""},"PeriodicalIF":4.0,"publicationDate":"2025-09-11","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145035524","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Evolution Characteristics of Principal Stress Transmission Path of Surrounding Rock in Adjacent Parallel Roadway Excavation","authors":"Yu Xiong, Dezhong Kong, Yujun Zuo, Zhijie Wen, Zhenhao Wang, Baizhi Mao, Qipeng Zhong","doi":"10.1002/nag.70069","DOIUrl":"https://doi.org/10.1002/nag.70069","url":null,"abstract":"In order to solve the problems that the three‐dimensional principal stress state and stress transmission path of surrounding rock are unclear in the process of parallel roadway excavation (PRE), based on field stress monitoring, elastic mechanics and numerical simulation, the principal stress data processing program is developed to study the magnitude gradual change and direction rotation law of maximum principal stress <jats:italic>σ</jats:italic><jats:sub>1</jats:sub>, intermediate principal stress <jats:italic>σ</jats:italic><jats:sub>2</jats:sub> and minimum principal stress <jats:italic>σ</jats:italic><jats:sub>3</jats:sub> in the process of adjacent PRE. The evolution characteristics of the principal stress transmission path of the roadway under superimposed disturbance are quantitatively described. Furthermore, the evolution difference of stress path and failure degree of surrounding rock between excavation direction parallel to <jats:italic>σ</jats:italic><jats:sub>1</jats:sub> and perpendicular to <jats:italic>σ</jats:italic><jats:sub>1</jats:sub> is discussed. The research results show that, during the PRE from left to right along <jats:italic>σ</jats:italic><jats:sub>1</jats:sub> direction, the stress rotation trajectories of the left and right roadsides have the characteristics of asymmetric distribution, and the stress peak value and rotation angle of the left roadside are larger than those of the right roadside. The principal stress at the shallow part of the roof and the roadside rotates approximately about <jats:italic>σ</jats:italic><jats:sub>2</jats:sub> and <jats:italic>σ</jats:italic><jats:sub>1</jats:sub> directions, respectively. If the excavation direction of the roadway is perpendicular to <jats:italic>σ</jats:italic><jats:sub>1</jats:sub>, the magnitude change amplitude and the direction rotation angle of principal stress are larger, which leads to more severe failure and instability of the surrounding rock. The research results emphasize the significance of reasonable excavation direction to control surrounding rock failure from the perspective of principal stress path.","PeriodicalId":13786,"journal":{"name":"International Journal for Numerical and Analytical Methods in Geomechanics","volume":"310 1","pages":""},"PeriodicalIF":4.0,"publicationDate":"2025-09-10","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145035525","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"A Unified Bond‐Based Peridynamic Model With Insights Into High‐Frequency Elastodynamic Problems in Geotechnical and Structural Engineering","authors":"Luyu Wang, Zhen‐Yu Yin","doi":"10.1002/nag.70062","DOIUrl":"https://doi.org/10.1002/nag.70062","url":null,"abstract":"The nonlocal effects in high‐frequency dynamic problems cannot be captured by classical continuum mechanics (CM), thereby introducing several compelling topics that warrant further exploration. Peridynamics (PD) offers a novel perspective for investigating these issues. This study proposes a unified bond‐based peridynamic (UBB‐PD) model, with an emphasis on high‐frequency dynamics that account for nonlocal properties. The UBB‐PD model incorporates a general criterion for constructing the micromodulus function. Then, the eigenfunction method is introduced to solve the UBB‐PD governing equations. The proposed model can naturally reduce into three different versions: CM, local PD model, and nonlocal PD model. The equivalence between PD and CM can be achieved by selecting an appropriate length‐scale parameter , with the wave frequency serving as a bridge connecting the two theories. Simulation results reveal that the local PD model perfectly reproduces the elastodynamic behavior found in CM across all frequencies. However, the nonlocal PD model inherently exhibits wave dispersion, arising from its nonlocal nature, which cannot be eliminated at high frequencies. PD stresses are affected by wave dispersion at high frequencies, with dissipative forces arising from inappropriate potentially inducing non‐conservation of mechanical energy. These results reveal findings not previously reported in relevant literature.","PeriodicalId":13786,"journal":{"name":"International Journal for Numerical and Analytical Methods in Geomechanics","volume":"66 1","pages":""},"PeriodicalIF":4.0,"publicationDate":"2025-09-09","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145035106","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"A Rate‐Dependent Stillinger‐Weber Potential‐Based Discretized Virtual Internal Bond Model for Rock Fracture Simulation Under Different Loading Rates","authors":"Yuezong Yang, Zihan Liu, Yujie Wang","doi":"10.1002/nag.70067","DOIUrl":"https://doi.org/10.1002/nag.70067","url":null,"abstract":"Strain rate has a notable influence on the mechanical properties of rocks. Understanding these properties and how they influence rock deformation and failure under various loading rates is critical for the engineering considerations in tunnel excavation, energy resource extraction, and mining operations. A rate‐dependent Stillinger‐Weber (SW) potential is constructed on the foundation of the trilinear damage potential, considering the strain rate effect on material properties. The developed potential is implanted into the discretized virtual internal bond (DVIB) model for rock fracture simulation at varying loading rates. A number of rock tests were simulated to see how accurately the developed model could simulate rock fracture over a range of strain rates. These rock tests included direct tension, uniaxial compression, three‐point bending, and semi‐circular bending tests. The findings reveal that both tensile and compressive strengths exhibit a linear increase with the logarithmic increase in strain rate. In the uniaxial compression tests, the failure modes transition from single‐plane shear to X‐shaped shear, then to extension, and ultimately to split failure as the strain rate escalates. In the semi‐circular bending tests, the mode I fracture toughness was observed to decrease with an increasing crack inclination angle, whereas the mode II fracture toughness exhibited the trend of increasing first and then decreasing. The rate‐dependent SW‐DVIB model has been demonstrated to accurately simulate the mechanical properties and failure behaviors of rocks across varying strain rates, as corroborated by the congruence between simulated and experimental results.","PeriodicalId":13786,"journal":{"name":"International Journal for Numerical and Analytical Methods in Geomechanics","volume":"20 1","pages":""},"PeriodicalIF":4.0,"publicationDate":"2025-09-08","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145017465","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Numerical Investigation on Consolidation Characteristics of Geotextile‐Encased Recycled Aggregate Columns in Soft Soil","authors":"Ling Zhang, Xiaocong Cai, Jingpeng Tan, Yaqian Shi, Binbing Mao","doi":"10.1002/nag.70064","DOIUrl":"https://doi.org/10.1002/nag.70064","url":null,"abstract":"The drainage and consolidation are important for the bearing capacity and deformation of geotextile‐encased recycled aggregate columns (GERACs) in soft soil. To investigate the consolidation characteristics of GERAC in soft soil, this paper develops a coupled numerical model incorporating the permeability of soft soil, the discreteness of RAs, and column drainage. The drainage process, vertical (<jats:italic>v</jats:italic><jats:sub>z</jats:sub>) and radial (<jats:italic>v</jats:italic><jats:sub>r</jats:sub>) permeable velocity, excess pore water pressure ratio (<jats:italic>R</jats:italic><jats:sub>EPWP</jats:sub>), vertical effective stress (<jats:italic>σ</jats:italic><jats:sub>zes</jats:sub>), and bulging deformation are methodically analyzed. Furthermore, the effect of eight essential factors on the consolidation characteristics is examined. Simulation results show that an influence boundary of bulging on the PWP and permeable velocity vector is 0.5 times the column length (<jats:italic>L</jats:italic>) in the depth direction and three times the column diameter (<jats:italic>D</jats:italic>) in the radial direction. The <jats:italic>v</jats:italic><jats:sub>r</jats:sub> relies on the bulging deformation of the column. The position of maximum bulging deformation changes from a depth of 0 to 0.67<jats:italic>D</jats:italic> and finally maintains at 0.67<jats:italic>D</jats:italic>. A slight bulging is observed at deeper positions of 1.50–3.00<jats:italic>D</jats:italic>. The gradation slightly affects the bulging deformation, but relatively significantly impacts the EPWP and <jats:italic>σ</jats:italic><jats:sub>zes</jats:sub>. A <jats:italic>L</jats:italic>/<jats:italic>D</jats:italic> of 5–9 and a relative density of not less than 55% is recommended.","PeriodicalId":13786,"journal":{"name":"International Journal for Numerical and Analytical Methods in Geomechanics","volume":"28 1","pages":""},"PeriodicalIF":4.0,"publicationDate":"2025-09-05","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145003093","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}