International Journal for Numerical Methods in Engineering最新文献

{"title":"A Unified Algorithm Framework for Phase-Field Based Multimaterial Topology Optimization With Various Graded Interfaces","authors":"Qian Yu,&nbsp;Qing Xia,&nbsp;Yibao Li,&nbsp;Chao Yang","doi":"10.1002/nme.70048","DOIUrl":"https://doi.org/10.1002/nme.70048","url":null,"abstract":"<div>\u0000 \u0000 <p>In this work, we present a unified algorithm framework for multimaterial topology optimization with various graded interfaces based on the phase field method. In the framework, the topology optimization problem is transformed into solving a set of linear partial differential equations (PDEs). We then introduce a scalar variable to grade the material property, resulting in a complex nonlinear interpolation operator in the effective elastic tensor. Following that, we define an auxiliary variable to substitute all the nonlinear transformations for a linear elliptic equation system. A second-order accurate Crank–Nicolson scheme is applied on the reformulated system, in which all high-order nonlinear terms are treated in a semi-explicit fashion. We prove that the proposed algorithm framework is unconditionally energy stable and demonstrate its stability as well as accuracy by numerical examples. A series of benchmark problems with different interfacial behaviors in topological design are investigated to verify the effectiveness of our method. The sensitivity of different parameters in the model is analyzed to evaluate their effects on the resulting structure.</p>\u0000 </div>","PeriodicalId":13699,"journal":{"name":"International Journal for Numerical Methods in Engineering","volume":"126 10","pages":""},"PeriodicalIF":2.7,"publicationDate":"2025-05-20","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144091845","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":"The Force Method Version of Goodman's Joint Element with Convergence Proof","authors":"Hong Zheng,&nbsp;Jia-han Jiang,&nbsp;Ming-kai Sun","doi":"10.1002/nme.70047","DOIUrl":"https://doi.org/10.1002/nme.70047","url":null,"abstract":"<div>\u0000 \u0000 <p>Interfaces widely present in nature and civil engineering are the most fundamental elements causing discontinuous deformation and failure. The Goodman joint element is incorporated into the vast majority of commercial and proprietary finite element programs due to its simplicity. However, the numerical properties of the Goodman element are not ideal. A lot of effort has gone into enhancing its numerical properties, only to make some empirical suggestions for adjusting mesh and computation parameters, yet there are still many examples that fail to reach convergence or incur drastic numerical oscillations of contact stress. Considering that the Goodman element is implemented within the framework of the displacement method, and poorlyposed problems in the displacement method are usually well-posed problems in the force method, this study proposes a force method version of the Goodman element, abbreviated as FMVGE. The primal unknowns in FMVGE are the contact stresses on the interface rather than nodal deformations. The core problem in FMVGE represents the contact conditions in the form of a quasi-variational inequality (QVI). By employing process iteration rather than state iteration used in the displacement method, at the same time, FMVGE precisely satisfies the contact conditions, with convergence being theoretically guaranteed. Analysis of classic examples and engineering cases indicates that the numerical properties of FMVGE are significantly superior to the widely adopted master–slave block method. Appendices A and B provide the proof of the convergence of the FMVGE solution and the core Matlab code for solving the QVI, respectively.</p>\u0000 </div>","PeriodicalId":13699,"journal":{"name":"International Journal for Numerical Methods in Engineering","volume":"126 10","pages":""},"PeriodicalIF":2.7,"publicationDate":"2025-05-15","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143950296","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":"Improved Conservative Phase Field Method for Contact Line Motion With Solidification","authors":"Mingguang Shen,&nbsp;Ben Q. Li","doi":"10.1002/nme.70036","DOIUrl":"https://doi.org/10.1002/nme.70036","url":null,"abstract":"<p>Phase field method is appealing due to its implicit capture of the interface. It mainly contains two types, one being the Allen–Cahn and the other being the Cahn-Hilliard. The latter conserves mass, but is fourth order, not as amenable to numerical methods as the former. Recently, a conservative Allen–Cahn phase field method has been put forward, improving the boundedness of the phase field. The new method conserves mass well and is only second order. This paper couples the conservative Allen–Cahn phase field model and the Navier–Stokes equation and the heat balance equation. Interfacial tension is expressed in a potential form. The model is discretized using a finite difference method on a half-staggered grid, and the velocity–pressure coupling is decoupled using an explicit projection method. The heat balance equation is solved using an iterative method. The model was tested against the classic static bubble case and was applied to impact in various conditions. Moreover, the model is parallelized using shared memory parallelism, OpenMP.</p>","PeriodicalId":13699,"journal":{"name":"International Journal for Numerical Methods in Engineering","volume":"126 9","pages":""},"PeriodicalIF":2.7,"publicationDate":"2025-05-14","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://onlinelibrary.wiley.com/doi/epdf/10.1002/nme.70036","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143944930","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 Novel Petrov-Galerkin 4-Node Quadrilateral Element With Radial Polynomial Interpolation for Linear Elastic Analysis","authors":"Lu-Zhen Dou,&nbsp;Ying-Qing Huang,&nbsp;Yuan-Fan Yang,&nbsp;Yan-Liang Ju,&nbsp;Hai-Bo Chen","doi":"10.1002/nme.70049","DOIUrl":"https://doi.org/10.1002/nme.70049","url":null,"abstract":"<div>\u0000 \u0000 <p>Based on the virtual work principle, a novel Petrov-Galerkin 4-node quadrilateral element is proposed in this article with different sets of test and trial functions. The virtual displacements are assumed by standard isoparametric interpolation which satisfies the interelement continuity requirement. It also ensures that the imposition of prescribed displacement boundary conditions and the calculation of equivalent nodal forces are the same as the conventional isoparametric elements. A support domain is formed for each element and all nodes within it are used to interpolate the actual displacements through the radial and polynomial basis functions. The nodal shape functions obtained by the radial polynomial interpolation possess the Kronecker delta property and sufficient completeness order for convergence. The resulted element stiffness matrix is unsymmetric, generally nonsquare, due to the Petrov-Galerkin formulation. However, the global unsymmetric stiffness matrix is square, sparse and structurally symmetric. Moreover, the determinant of the Jacobian matrix can be removed from the element stiffness matrix and this improves the immunity of the numerical accuracy to mesh distortion significantly. Numerical investigations demonstrate that the present element effectively combines the advantages of both finite element methods and radial basis meshless methods. Especially, the stress continuity and interelement smoothness are improved remarkedly.</p>\u0000 </div>","PeriodicalId":13699,"journal":{"name":"International Journal for Numerical Methods in Engineering","volume":"126 10","pages":""},"PeriodicalIF":2.7,"publicationDate":"2025-05-14","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143944707","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":"Cell Agglomeration Strategy for Cut Cells in eXtended Discontinuous Galerkin Methods","authors":"Teoman Toprak,&nbsp;Matthias Rieckmann,&nbsp;Florian Kummer","doi":"10.1002/nme.70013","DOIUrl":"https://doi.org/10.1002/nme.70013","url":null,"abstract":"<p>In this work, a cell agglomeration strategy for the cut cells arising in the eXtended discontinuous Galerkin (XDG) method is presented. Cut cells are a fundamental aspect of unfitted mesh approaches, where complex geometries or interfaces separating subdomains are embedded into structured background grids to facilitate the mesh generation process. In such methods, arbitrary small cells occur due to the intersections of background cells with embedded geometries and lead to discretization difficulties due to their diminutive sizes. Furthermore, temporal evolutions of these geometries may lead to topological changes across different time steps. Both of these issues, that is, small-cut cells and topological changes, can be addressed with a cell agglomeration technique, independent of discretization. However, cell agglomeration encounters significant difficulties in three dimensions due to the complexity of neighborship and issues like cycles and parallel agglomeration chains. The proposed strategy introduces a robust framework that mitigates these problems by incorporating methods for cycle prevention, chain agglomeration, and parallelization. Implemented in the open-source software package BoSSS, this strategy has been successfully tested on multiprocessor systems using dynamic multiphase test cases in both two and three dimensions, enabling simulations that were previously infeasible.</p>","PeriodicalId":13699,"journal":{"name":"International Journal for Numerical Methods in Engineering","volume":"126 9","pages":""},"PeriodicalIF":2.7,"publicationDate":"2025-05-08","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://onlinelibrary.wiley.com/doi/epdf/10.1002/nme.70013","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143926029","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 Fundamentally New Coupled Approach to Contact Mechanics via the Dirichlet-Neumann Schwarz Alternating Method","authors":"Alejandro Mota,&nbsp;Daria Koliesnikova,&nbsp;Irina Tezaur,&nbsp;Jonathan Hoy","doi":"10.1002/nme.70039","DOIUrl":"https://doi.org/10.1002/nme.70039","url":null,"abstract":"<div>\u0000 \u0000 <p>Contact phenomena are crucial for understanding the behavior of mechanical systems. However, existing computational approaches for simulating mechanical contact often face numerical challenges, such as inaccurate physical predictions, energy conservation errors, and unwanted oscillations. We introduce an alternative technique for simulating dynamic contact based on the non-overlapping Schwarz alternating method, originally developed for domain decomposition. In multibody contact scenarios, this method treats each body as a separate, non-overlapping domain and prevents interpenetration using an alternating Dirichlet–Neumann iterative process. This approach has a strong theoretical foundation, eliminates the need for contact constraints, and offers flexibility, making it ideal for multiscale and multiphysics applications. We conducted a numerical comparison between the Schwarz method and traditional methods, such as the Lagrange multiplier and penalty methods, focusing on a benchmark impact problem. Our results indicate that the Schwarz alternating method outperforms traditional methods in several key areas: it provides more accurate predictions for various measurable quantities and demonstrates exceptional energy conservation capabilities. To address unwanted oscillations in contact velocities and forces, we explored various algorithms and stabilization techniques, ultimately opting for the naïve-stabilized Newmark scheme for its simplicity and effectiveness. Additionally, we validated the efficiency of the Schwarz method in a three-dimensional impact problem, highlighting its inherent capacity to accommodate different mesh topologies, time-integration schemes, and time steps for each interacting body.</p>\u0000 </div>","PeriodicalId":13699,"journal":{"name":"International Journal for Numerical Methods in Engineering","volume":"126 9","pages":""},"PeriodicalIF":2.7,"publicationDate":"2025-05-08","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143919349","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":"Bridging Overlapping 2D Coarse and Fine Meshes Within the Phase Field Fracture Method","authors":"Zakaria Chafia,&nbsp;Julien Yvonnet,&nbsp;Jérémy Bleyer","doi":"10.1002/nme.70043","DOIUrl":"https://doi.org/10.1002/nme.70043","url":null,"abstract":"<p>A framework is proposed to bridge coarse and fine meshes in a single simulation within the phase field method for fracture. Fine meshes are used in the vicinity of localized defects to accurately capture crack initiation, while coarse meshes are used away from initial defects and include only crack propagation paths. This reduces the prohibitive computational times associated with uniformly fine meshes over the entire domain, or with the use of complex adaptive meshes in the phase field method. The coupling between the two overlapping meshes is achieved using a variational formulation in which the energies of the models associated with the fine and coarse meshes are weighted in the superposition zone. Two situations are considered. The first includes the resolution of the phase field problem only in the fine mesh, while the coarse mesh is limited to the undamaged elastic problem. In the second situation, cracks can propagate in both the fine and coarse mesh. Variational formulations and associated finite element implementations are detailed. Numerical examples are presented, showing the potential of this approach to significantly reduce computational costs in the phase field method for cracking without affecting the accuracy.</p>","PeriodicalId":13699,"journal":{"name":"International Journal for Numerical Methods in Engineering","volume":"126 9","pages":""},"PeriodicalIF":2.7,"publicationDate":"2025-05-07","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://onlinelibrary.wiley.com/doi/epdf/10.1002/nme.70043","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143914446","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":"Data-Driven Optimization for the Evolve-Filter-Relax Regularization of Convection-Dominated Flows","authors":"Anna Ivagnes,&nbsp;Maria Strazzullo,&nbsp;Michele Girfoglio,&nbsp;Traian Iliescu,&nbsp;Gianluigi Rozza","doi":"10.1002/nme.70042","DOIUrl":"https://doi.org/10.1002/nme.70042","url":null,"abstract":"&lt;p&gt;Numerical stabilization techniques are often employed in under-resolved simulations of convection-dominated flows to improve accuracy and mitigate spurious oscillations. Specifically, the evolve–filter–relax (EFR) algorithm is a framework that consists of evolving the solution, applying a filtering step to remove high-frequency noise, and relaxing through a convex combination of filtered and original solutions. The stability and accuracy of the EFR solution strongly depend on two parameters, the filter radius &lt;span&gt;&lt;/span&gt;&lt;math&gt;\u0000 &lt;semantics&gt;\u0000 &lt;mrow&gt;\u0000 &lt;mi&gt;δ&lt;/mi&gt;\u0000 &lt;/mrow&gt;\u0000 &lt;annotation&gt;$$ delta $$&lt;/annotation&gt;\u0000 &lt;/semantics&gt;&lt;/math&gt; and the relaxation parameter &lt;span&gt;&lt;/span&gt;&lt;math&gt;\u0000 &lt;semantics&gt;\u0000 &lt;mrow&gt;\u0000 &lt;mi&gt;χ&lt;/mi&gt;\u0000 &lt;/mrow&gt;\u0000 &lt;annotation&gt;$$ chi $$&lt;/annotation&gt;\u0000 &lt;/semantics&gt;&lt;/math&gt;. Standard choices for these parameters are usually fixed in time, and related to the full order model setting, that is, the grid size for &lt;span&gt;&lt;/span&gt;&lt;math&gt;\u0000 &lt;semantics&gt;\u0000 &lt;mrow&gt;\u0000 &lt;mi&gt;δ&lt;/mi&gt;\u0000 &lt;/mrow&gt;\u0000 &lt;annotation&gt;$$ delta $$&lt;/annotation&gt;\u0000 &lt;/semantics&gt;&lt;/math&gt; and the time step for &lt;span&gt;&lt;/span&gt;&lt;math&gt;\u0000 &lt;semantics&gt;\u0000 &lt;mrow&gt;\u0000 &lt;mi&gt;χ&lt;/mi&gt;\u0000 &lt;/mrow&gt;\u0000 &lt;annotation&gt;$$ chi $$&lt;/annotation&gt;\u0000 &lt;/semantics&gt;&lt;/math&gt;. The key novelties with respect to the standard EFR approach are: (i) time-dependent parameters &lt;span&gt;&lt;/span&gt;&lt;math&gt;\u0000 &lt;semantics&gt;\u0000 &lt;mrow&gt;\u0000 &lt;mi&gt;δ&lt;/mi&gt;\u0000 &lt;mo&gt;(&lt;/mo&gt;\u0000 &lt;mi&gt;t&lt;/mi&gt;\u0000 &lt;mo&gt;)&lt;/mo&gt;\u0000 &lt;/mrow&gt;\u0000 &lt;annotation&gt;$$ delta (t) $$&lt;/annotation&gt;\u0000 &lt;/semantics&gt;&lt;/math&gt; and &lt;span&gt;&lt;/span&gt;&lt;math&gt;\u0000 &lt;semantics&gt;\u0000 &lt;mrow&gt;\u0000 &lt;mi&gt;χ&lt;/mi&gt;\u0000 &lt;mo&gt;(&lt;/mo&gt;\u0000 &lt;mi&gt;t&lt;/mi&gt;\u0000 &lt;mo&gt;)&lt;/mo&gt;\u0000 &lt;/mrow&gt;\u0000 &lt;annotation&gt;$$ chi (t) $$&lt;/annotation&gt;\u0000 &lt;/semantics&gt;&lt;/math&gt;, and (ii) data-driven adaptive optimization of the parameters in time, considering a fully-resolved simulation as reference. In particular, we propose three different classes of optimized-EFR (&lt;b&gt;Opt-EFR&lt;/b&gt;) strategies, aiming to optimize one or both parameters. The new Opt-EFR strategies are tested in the under-resolved simulation of a turbulent flow past a cylinder at &lt;span&gt;&lt;/span&gt;&lt;math&gt;\u0000 &lt;semantics&gt;\u0000 &lt;mrow&gt;\u0000 &lt;mi&gt;R&lt;/mi&gt;\u0000 &lt;mi&gt;e&lt;/mi&gt;\u0000 &lt;mo&gt;=&lt;/mo&gt;\u0000 &lt;mn&gt;1&lt;/mn&gt;\u0000 &lt;mo&gt;,&lt;/mo&gt;\u0000 &lt;mn&gt;000&lt;/mn&gt;\u0000 &lt;/mrow&gt;\u0000 ","PeriodicalId":13699,"journal":{"name":"International Journal for Numerical Methods in Engineering","volume":"126 9","pages":""},"PeriodicalIF":2.7,"publicationDate":"2025-05-07","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://onlinelibrary.wiley.com/doi/epdf/10.1002/nme.70042","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143919452","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":"Component-Based Model-Order Reduction With Mortar Tied Contact for Nonlinear Quasi-Static Mechanical Problems","authors":"Stephan Ritzert,&nbsp;Jannick Kehls,&nbsp;Stefanie Reese,&nbsp;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,&nbsp;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}