Stephan Ritzert, Jannick Kehls, Stefanie Reese, Tim Brepols
{"title":"Component-Based Model-Order Reduction With Mortar Tied Contact for Nonlinear Quasi-Static Mechanical Problems","authors":"Stephan Ritzert, Jannick Kehls, Stefanie Reese, Tim Brepols","doi":"10.1002/nme.70041","DOIUrl":"https://doi.org/10.1002/nme.70041","url":null,"abstract":"<p>In this work, we present a model-order reduction technique for nonlinear structures assembled from components. The reduced- order model is constructed by reducing the substructures with proper orthogonal decomposition and connecting them by a mortar-tied contact formulation. The substructure projection matrices are computed by the proper orthogonal decomposition (POD) method from snapshots computed on the substructure level. The snapshots are computed using Latin hypercube sampling based on a parametrization of the boundary conditions. In numerical examples, we show the accuracy and efficiency of the method for nonlinear problems involving material and geometric nonlinearities as well as non-matching meshes. The method can predict solutions of new systems with varying boundary conditions and material behaviors.</p>","PeriodicalId":13699,"journal":{"name":"International Journal for Numerical Methods in Engineering","volume":"126 8","pages":""},"PeriodicalIF":2.7,"publicationDate":"2025-04-27","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://onlinelibrary.wiley.com/doi/epdf/10.1002/nme.70041","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143879987","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"KATO: Neural-Reparameterized Topology Optimization Using Convolutional Kolmogorov-Arnold Network for Stress Minimization","authors":"Shengyu Yan, Jasmin Jelovica","doi":"10.1002/nme.70034","DOIUrl":"https://doi.org/10.1002/nme.70034","url":null,"abstract":"<p>Topology optimization (TO) has been a cornerstone of advanced structural design for decades, yet it continues to face challenges in terms of convergence, optimality, and numerical stability, particularly for complex, non-convex problems like stress minimization. This paper introduces a novel approach to stress-based topology optimization through the development of neural-reparameterized topology optimization using the convolutional Kolmogorov-Arnold network (KATO). KATO uses the neural network to reparameterize the optimization problem, offering a unique solution to the challenges posed by stress minimization in TO. It also simplifies the penalization scheme by reducing sensitivity to certain parameters, which reduces the non-convexity of the stress minimization problem, enhancing convergence and stability. Our method demonstrates better performance in stress minimization compared to conventional approaches and a different neural network-based approach, achieving up to 10% lower maximum stress in common benchmark cases. KATO also shows remarkable efficiency, reducing computational time by up to 67% compared to conventional methods for stress minimization problems. We conduct a comprehensive analysis of KATO's performance, computational cost, scalability, and the impact of various neural network architectures. Our results indicate that KATO not only improves stress optimization but also offers insights into the relationship between neural network design and topology optimization performance, paving the way for more efficient and effective structural design processes.</p>","PeriodicalId":13699,"journal":{"name":"International Journal for Numerical Methods in Engineering","volume":"126 8","pages":""},"PeriodicalIF":2.7,"publicationDate":"2025-04-21","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://onlinelibrary.wiley.com/doi/epdf/10.1002/nme.70034","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143852806","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Fluid-Structure Interaction Simulation of Mitral Valve Structures in a Left Ventricle Model","authors":"Joel Kronborg, Johan Hoffman","doi":"10.1002/nme.70031","DOIUrl":"https://doi.org/10.1002/nme.70031","url":null,"abstract":"<p>Simulations of blood flow in patient-specific models of heart ventricles is a rapidly developing field of research, showing promise to improve future treatment of heart diseases. Fluid-structure interaction simulation of the mitral valve, with its complex structure including leaflets, chordae tendineae, and papillary muscles, provides additional prospects as well as challenges to such models. In this study, we combine a patient-specific model of the left ventricle with an idealized unified continuum fluid-structure interaction model of the mitral valve, to simulate the intraventricular diastolic blood flow. To the best of our knowledge, no monolithic fluid-structure interaction model, without the need for remeshing, has ever been used before to simulate the native mitral valve within the left ventricle. The chordae tendineae are simulated as a region of porous medium, to partially hinder the flow. Simulation results from this model are compared to those of a model with the same patient-specific left ventricle, but with the mitral valve simply modeled as a time-variant inflow boundary condition. The blood flow is analyzed with the E-wave propagation index, and by use of the triple decomposition of the velocity gradient tensor, which decomposes the flow into rigid body rotational flow, shearing flow, and irrotational straining flow. The triple decomposition enables analysis of the formation of initially large dominant flow features, such as the E-wave jet and the vortex ring around it, and their subsequent decay into smaller turbulent flow structures. This analysis of the development of flow structures over the duration of diastole appears to be in general agreement with the theory of the stability of rotation, shear, and strain structures. Elevated shear levels are investigated, but are found only in limited amounts that do not indicate significant risks of thrombus formation or other blood damage, which is to be expected in this healthy ventricle. The highest shear levels are localized at the leaflets in the fluid-structure interaction model, and at the ventricle wall in the planar model. The computed E-wave propagation indices are 1.21 for the fluid-structure interaction model and 1.90 for the planar valve model, which indicates proper washout in the apical region with no significant risk of blood stasis that could lead to left ventricular thrombus formation.</p>","PeriodicalId":13699,"journal":{"name":"International Journal for Numerical Methods in Engineering","volume":"126 8","pages":""},"PeriodicalIF":2.7,"publicationDate":"2025-04-21","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://onlinelibrary.wiley.com/doi/epdf/10.1002/nme.70031","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143852804","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Zhiyong Zhao, Hao Yuan, Jaan-Willem Simon, Lishuai Sun, Yujun Li
{"title":"Forward and Inverse Approaches for Uncertainty Quantification of the In-Plane Elastic Properties of Cellular Structures","authors":"Zhiyong Zhao, Hao Yuan, Jaan-Willem Simon, Lishuai Sun, Yujun Li","doi":"10.1002/nme.70035","DOIUrl":"https://doi.org/10.1002/nme.70035","url":null,"abstract":"<div>\u0000 \u0000 <p>Uncertainty quantification is essential to exploiting the complete potential of cellular structures. Forward methods allow quantifying the uncertainties in the mechanical properties of cellular structures by propagating the uncertainties of the input parameters, while inverse methods allow using experimental data to indirectly infer the uncertainty of the input parameters. In this paper, a closed-loop forward and inverse quantification of the in-plane elastic properties of cellular structures was proposed. A polynomial chaos expansion model was used in the forward model for uncertainty propagation and quantification of the in-plane elastic properties using the results from the Fast Fourier Transform simulations. The random input parameters field in the Fast Fourier Transform simulations, including geometry and material parameters of cellular structures, was involved by Karhunen–Loève expansion. Furthermore, the inverse uncertainty quantification was conducted in the framework of Markov Chain Monte Carlo sampling-based Bayesian inference using the constructed polynomial chaos surrogate model. The approach introduced was applied to analyze the uncertainty quantification in two types of cellular structures. The results showed that the thickness of the cell wall dramatically influences the effective in-plane elastic modulus of the cellular structures. The PCE could significantly reduce the iterations compared to the Monte Carlo simulation while ensuring the accuracy of uncertainty quantification of the in-plane elastic modulus. In addition, effective evaluation and calibration of the geometry and material parameters of the cellular structures based on the obtained posterior probability distribution have been achieved. This addresses the problem of uncertainty quantification of the in-plane elastic properties and the difficulty in measuring the geometry and material parameters of cellular structures.</p>\u0000 </div>","PeriodicalId":13699,"journal":{"name":"International Journal for Numerical Methods in Engineering","volume":"126 8","pages":""},"PeriodicalIF":2.7,"publicationDate":"2025-04-21","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143852863","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"A Novel Strain-Continuous Finite Element Formulation for Geometrically Exact Thin-Walled Beam on Lie Algebra","authors":"Ziheng Huang, Ju Chen, Shixing Liu, Yongxin Guo","doi":"10.1002/nme.70024","DOIUrl":"https://doi.org/10.1002/nme.70024","url":null,"abstract":"<div>\u0000 \u0000 <p>Most existing geometrically exact thin-walled beam formulations on Lie group SE(3)<span></span><math>\u0000 <semantics>\u0000 <mrow>\u0000 <mo>×</mo>\u0000 </mrow>\u0000 <annotation>$$ times $$</annotation>\u0000 </semantics></math><span></span><math>\u0000 <semantics>\u0000 <mrow>\u0000 <mi>ℝ</mi>\u0000 </mrow>\u0000 <annotation>$$ mathbb{R} $$</annotation>\u0000 </semantics></math> considering warping are <span></span><math>\u0000 <semantics>\u0000 <mrow>\u0000 <mi>C</mi>\u0000 </mrow>\u0000 <annotation>$$ C $$</annotation>\u0000 </semantics></math><span></span><math>\u0000 <semantics>\u0000 <mrow>\u0000 <msup>\u0000 <mo> </mo>\u0000 <mrow>\u0000 <mn>0</mn>\u0000 </mrow>\u0000 </msup>\u0000 </mrow>\u0000 <annotation>$$ {}^0 $$</annotation>\u0000 </semantics></math>-continuous. In this study, a novel strain-continuous element for geometrically exact thin-walled beam on Lie algebra <span></span><math>\u0000 <semantics>\u0000 <mrow>\u0000 <mi>s</mi>\u0000 </mrow>\u0000 <annotation>$$ mathfrak{s} $$</annotation>\u0000 </semantics></math><span></span><math>\u0000 <semantics>\u0000 <mrow>\u0000 <mi>e</mi>\u0000 </mrow>\u0000 <annotation>$$ mathfrak{e} $$</annotation>\u0000 </semantics></math>(3)<span></span><math>\u0000 <semantics>\u0000 <mrow>\u0000 <mo>×</mo>\u0000 </mrow>\u0000 <annotation>$$ times $$</annotation>\u0000 </semantics></math><span></span><math>\u0000 <semantics>\u0000 <mrow>\u0000 <mi>ℝ</mi>\u0000 </mrow>\u0000 <annotation>$$ mathbb{R} $$</annotation>\u0000 </semantics></math> is originally proposed, in which the torsion-related Wagner effect and warping are considered in the constitutive relations. The proposed beam element is not only locking-free intrinsically, but also is <span></span><math>\u0000 <semantics>\u0000 <mrow>\u0000 <mi>C</mi>\u0000 </mrow>\u0000 <annotation>$$ C $$</annotation>\u0000 </semantics></math><span></span><math>\u0000 <semantics>\u0000 <mrow>\u0000 <msup>\u0000 <mo> </mo>\u0000 <mrow>\u0000 <mn>1</mn>\u0000 </mrow>\u0000 </msup>\u0000 </mrow>\u0000 <annotation>$$ {}^1 $$</annotation>\u0000 </semantics></math>-continuous by using the proposed geometrical Hermite interpolation on Lie algebra. This interpolation based on Magnus expansion can simplify the description of interpolated element stress tensor and strain energy. Then, the sim","PeriodicalId":13699,"journal":{"name":"International Journal for Numerical Methods in Engineering","volume":"126 8","pages":""},"PeriodicalIF":2.7,"publicationDate":"2025-04-21","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143852803","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"A Deep Learning-Based Surrogate Model for Seismic Data Assimilation in Fault Activation Modeling","authors":"Caterina Millevoi, Claudia Zoccarato, Massimiliano Ferronato","doi":"10.1002/nme.70040","DOIUrl":"https://doi.org/10.1002/nme.70040","url":null,"abstract":"<p>Assessing the safety and environmental impacts of subsurface resource exploitation and management is critical and requires robust geomechanical modeling. However, uncertainties stemming from model assumptions, intrinsic variability of governing parameters, and data errors challenge the reliability of predictions. In the absence of direct measurements, inverse modeling and stochastic data assimilation methods can offer reliable solutions, but in complex and large-scale settings, the computational expense can become prohibitive. To address these challenges, this paper presents a deep learning-based surrogate model (SurMoDeL) designed for seismic data assimilation in fault activation modeling. The surrogate model leverages neural networks to provide simplified yet accurate representations of complex geophysical systems, enabling faster simulations and analyses essential for uncertainty quantification. The work proposes two different methods to integrate an understanding of fault behavior into the model, thereby enhancing the accuracy of its predictions. The application of the proxy model to integrate seismic data through effective data assimilation techniques efficiently constrains the uncertain parameters, thus bridging the gap between theoretical models and real-world observations.</p>","PeriodicalId":13699,"journal":{"name":"International Journal for Numerical Methods in Engineering","volume":"126 8","pages":""},"PeriodicalIF":2.7,"publicationDate":"2025-04-21","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://onlinelibrary.wiley.com/doi/epdf/10.1002/nme.70040","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143852802","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"A Chemo-Damage-Mechanical Coupled Phase-Field Model for Three-Dimensional Hydrogen-Assisted Dynamic Cracking","authors":"Hui Li, Shanyong Wang","doi":"10.1002/nme.70038","DOIUrl":"https://doi.org/10.1002/nme.70038","url":null,"abstract":"<p>This study develops a chemo-damage-mechanical coupled phase-field method for modeling two-dimensional and/or three-dimensional hydrogen-assisted transient dynamic cracking in metallic materials. In this method, hydrogen diffusion in solids is described by the evolution of bulk hydrogen concentration governed by the diffusion equation with an extended Fick's law. The hydrogen concentration and the inertial force of solids are incorporated into the governing equations of a non-standard quasi-brittle phase-field model (known as the phase-field regularized cohesive zone model, PF-CZM) capable of modeling complicated multiple crack nucleation, initiation, and propagation with insensitivity to mesh size and length scale. The resultant displacement-damage-hydrogen concentration coupled three-field equation is derived by the finite element Galerkin method and solved using a staggered Newton–Raphson iteration algorithm with an unconditionally stable implicit Newmark integration scheme. The new method was verified by four benchmark examples with hydrogen-free numerical/experimental results for comparison, including quasi-static/dynamic fracture of a notched plate under uniaxial tension, a dynamic crack branching experiment under constant traction, the Kalthoff-Winkler impact fracture experiment, and dynamic fragmentation of a cylinder under internal pressures, with the effects of loading velocity, initial bulk hydrogen concentration, and diffusion time on crack propagation investigated in detail. It is found that the present method is capable of modeling complex 2D and 3D hydrogen-assisted dynamic crack propagation and bifurcation under impact or internal pressure loadings, and thus holds the potential to be used for the structural design of hydrogen storage.</p>","PeriodicalId":13699,"journal":{"name":"International Journal for Numerical Methods in Engineering","volume":"126 8","pages":""},"PeriodicalIF":2.7,"publicationDate":"2025-04-21","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://onlinelibrary.wiley.com/doi/epdf/10.1002/nme.70038","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143852865","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Yiwei Feng, Lili Lv, Tiegang Liu, Kun Wang, Bangcheng Ai
{"title":"Investigation of Using High-Order Discontinuous Galerkin Methods in Adjoint Gradient-Based 3D Aerodynamic Shape Optimization","authors":"Yiwei Feng, Lili Lv, Tiegang Liu, Kun Wang, Bangcheng Ai","doi":"10.1002/nme.70033","DOIUrl":"https://doi.org/10.1002/nme.70033","url":null,"abstract":"<div>\u0000 \u0000 <p>Based on recent work by He et al. (<i>Computer Physics Communications</i> <b>286</b> (2023): 108660), this work develops a robust and efficient workflow for 3D aerodynamic shape optimization (ASO) based on the discontinuous Galerkin method (DGM) with solution remapping technique. It is theoretically found and numerically validated through 2D cases that the DGMs can provide more precise adjoint-enabled gradients even on a coarse mesh as compared with the FVMs under coequal computational costs. A number of 2D and 3D intricate cases are tested to illustrate the potential advantages of the DGM-based ASO workflow in terms of computational efficiency and final optimization results. The results demonstrate that the solution remapping technique can save around <span></span><math>\u0000 <semantics>\u0000 <mrow>\u0000 <mn>50</mn>\u0000 <mo>%</mo>\u0000 </mrow>\u0000 <annotation>$$ 50% $$</annotation>\u0000 </semantics></math> of the computational time. Moreover, the DGM-based workflow extends the capability to explore the design space, leading to aerodynamic shapes with superior performance.</p>\u0000 </div>","PeriodicalId":13699,"journal":{"name":"International Journal for Numerical Methods in Engineering","volume":"126 8","pages":""},"PeriodicalIF":2.7,"publicationDate":"2025-04-21","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143852805","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Adaptive Coupling of Peridynamic and Classical Continuum Mechanical Models Driven by Broken Bond/Strength Criteria for Structural Dynamic Failure","authors":"JiuYi Li, ShanKun Liu, Fei Han, Yong Mei, YunHou Sun, FengJun Zhou","doi":"10.1002/nme.70021","DOIUrl":"https://doi.org/10.1002/nme.70021","url":null,"abstract":"<div>\u0000 \u0000 <p>Peridynamics (PD) is widely used to simulate structural failure. However, PD models are time consuming. To improve the computational efficiency, we developed an adaptive coupling model between PD and classical continuum mechanics (PD-CCM) based on the Morphing method, driven by the broken bond or strength criteria. We derived the dynamic equation of the coupled models from the Lagrangian equation and then the discretized finite element formulation. An adaptive coupling strategy was introduced by determining the key position using the broken bond or strength criteria. The PD subdomain was expanded by altering the value of the Morphing function around the key position. Additionally, the PD subdomain was meshed by discrete elements (DEs) (i.e., nodes were not shared between elements), allowing the crack to propagate freely along the boundary of the DE. The remaining subdomains were meshed by continuous elements (CEs). Following the PD subdomain expansion, the CEs were converted into DEs, and new nodes were inserted. The displacement vector and mass matrix were reconfigured to ensure calculation consistency throughout the solving process. Furthermore, the relationship between the expansion radius of the PD subdomain and the speed of crack propagation was also discussed. Finally, the effectiveness, efficiency, and accuracy of the proposed model were verified via three two-dimensional numerical examples.</p>\u0000 </div>","PeriodicalId":13699,"journal":{"name":"International Journal for Numerical Methods in Engineering","volume":"126 7","pages":""},"PeriodicalIF":2.7,"publicationDate":"2025-04-08","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143793300","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Novel Study on Strain Modes-Based Interval Damage Identification Methodology Utilizing Orthogonal Polynomials and Collocation Theories","authors":"Lei Wang, Lihan Cheng, Qinghe Shi","doi":"10.1002/nme.70032","DOIUrl":"https://doi.org/10.1002/nme.70032","url":null,"abstract":"<div>\u0000 \u0000 <p>The sensitivity of strain modes to local stiffness changes within a structure underscores their potential as robust indicators of damage, enhancing the efficacy of damage identification processes. This study establishes sensitivity matrices of natural frequencies and strain modes to damage parameters, laying the groundwork for a novel fusion index that integrates both metrics to assess structural damage extent. In order to quantify the impact of uncertainty information on the results of damage identification processes, a non-probabilistic structural damage identification method rooted in the collocation methodology is proposed in this study. In consideration of computational efficiency, a two-step damage identification strategy encompassing localization and quantification is proposed. Initially, damage localization is achieved through the dynamic fingerprints, followed by the quantification of the uncertainty of damage extent. The proposed methodology is validated through a detailed numerical example, illustrating that the fusion index outperforms individual indices in terms of accuracy and computational efficiency. The non-probabilistic structural damage identification method based on collocation methodology can identify the damage extent and uncertainty interval even under the influence of uncertain factors.</p>\u0000 </div>","PeriodicalId":13699,"journal":{"name":"International Journal for Numerical Methods in Engineering","volume":"126 7","pages":""},"PeriodicalIF":2.7,"publicationDate":"2025-04-06","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143786701","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}