Engineering Computations最新文献

{"title":"Multivariate nonlinear regression analysis of hydraulic fracturing parameters based on hybrid FEM–DEM","authors":"Yang Li, Tianxiang Lan","doi":"10.1108/ec-06-2023-0270","DOIUrl":"https://doi.org/10.1108/ec-06-2023-0270","url":null,"abstract":"Purpose This paper aims to employ a multivariate nonlinear regression analysis to establish a predictive model for the final fracture area, while accounting for the impact of individual parameters. Design/methodology/approach This analysis is based on the numerical simulation data obtained, using the hybrid finite element–discrete element (FE–DE) method. The forecasting model was compared with the numerical results and the accuracy of the model was evaluated by the root mean square (RMS) and the RMS error, the mean absolute error and the mean absolute percentage error. Findings The multivariate nonlinear regression model can accurately predict the nonlinear relationships between injection rate, leakoff coefficient, elastic modulus, permeability, Poisson’s ratio, pore pressure and final fracture area. The regression equations obtained from the Newton iteration of the least squares method are strong in terms of the fit to the six sensitive parameters, and the model follow essentially the same trend with the numerical simulation data, with no systematic divergence detected. Least absolutely deviation has a significantly weaker performance than the least squares method. The percentage contribution of sensitive parameters to the final fracture area is available from the simulation results and forecast model. Injection rate, leakoff coefficient, permeability, elastic modulus, pore pressure and Poisson’s ratio contribute 43.4%, −19.4%, 24.8%, −19.2%, −21.3% and 10.1% to the final fracture area, respectively, as they increased gradually. In summary, (1) the fluid injection rate has the greatest influence on the final fracture area. (2)The multivariate nonlinear regression equation was optimally obtained after 59 iterations of the least squares-based Newton method and 27 derivative evaluations, with a decidability coefficient R2 = 0.711 representing the model reliability and the regression equations fit the four parameters of leakoff coefficient, permeability, elastic modulus and pore pressure very satisfactorily. The models follow essentially the identical trend with the numerical simulation data and there is no systematic divergence. The least absolute deviation has a significantly weaker fit than the least squares method. (3)The nonlinear forecasting model of physical parameters of hydraulic fracturing established in this paper can be applied as a standard for optimizing the fracturing strategy and predicting the fracturing efficiency in situ field and numerical simulation. Its effectiveness can be trained and optimized by experimental and simulation data, and taking into account more basic data and establishing regression equations, containing more fracturing parameters will be the further research interests. Originality/value The nonlinear forecasting model of physical parameters of hydraulic fracturing established in this paper can be applied as a standard for optimizing the fracturing strategy and predicting the fracturing efficiency in situ field ","PeriodicalId":50522,"journal":{"name":"Engineering Computations","volume":"53 48","pages":"0"},"PeriodicalIF":0.0,"publicationDate":"2023-11-13","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"134992736","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Biorthogonal flatlet multiwavelet collocation method for solving the singular nonlinear system with initial and boundary conditions","authors":"Maryam Mohseni, Davood Rostamy","doi":"10.1108/ec-12-2022-0730","DOIUrl":"https://doi.org/10.1108/ec-12-2022-0730","url":null,"abstract":"Purpose The numerical methods are of great importance for approximating the solutions of a system of nonlinear singular ordinary differential equations. In this paper, the authors present the biorthogonal flatlet multiwavelet collocation method (BFMCM) as a numerical scheme for a class of system of Lane–Emden equations with initial or boundary or four-point boundary conditions. Design/methodology/approach The approach is involved in combining the biorthogonal flatlet multiwavelet (BFM) with the collocation method. The authors investigate the properties and procedure of the BFMCM for first time on this class of equations. By using the BFM and the collocation points, the method is constructed and it transforms the nonlinear differential equations problem into a system of nonlinear algebraic equations. The unknown coefficients of the assuming solution are determined by solving the obtained system. Additionally, convergence analysis and numerical stability of the suggested method are provided. Findings According to the attained results, the proposed BFMCM has more accurate results in comparison with results of other methods. The maximum absolute errors are calculated by using the BFMCM for comparison purposes provided. Originality/value The key desirable properties of BFMCM are its efficiency, simple applicability and minimizes errors. Therefore, the proposed method can be used to solve nonlinear problems or problems with singular points.","PeriodicalId":50522,"journal":{"name":"Engineering Computations","volume":"56 51","pages":"0"},"PeriodicalIF":0.0,"publicationDate":"2023-11-13","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"134992716","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Multi-objective design optimization of reentrant auxetic model using Lichtenberg algorithm based on metamodel","authors":"Matheus Francisco, João Pereira, Lucas Oliveira, Sebastião Simões Cunha, G.F. Gomes","doi":"10.1108/ec-07-2022-0448","DOIUrl":"https://doi.org/10.1108/ec-07-2022-0448","url":null,"abstract":"Purpose The present paper aims at the multi-objective optimization of a reentrant hexagonal cell auxetic structure. In addition, a parametric analysis will be carried out to verify how each of the design factors impact each of the responses. Design/methodology/approach The multi-objective optimization of five different responses of an auxetic model was considered: mass, critical buckling load under compression effort, natural frequency, Poisson's ratio and failure load. The response surface methodology was applied, and a new meta-heuristic of optimization called the multi-objective Lichtenberg algorithm was applied to find the optimized configuration of the model. It was possible to increase the failure load by 26.75% in compression performance optimization. Furthermore, in the optimization of modal performance, it was possible to increase the natural frequency by 37.43%. Finally, all 5 responses analyzed simultaneously were optimized. In this case, it was possible to increase the critical buckling load by 42.55%, the failure load by 28.70% and reduce the mass and Poisson's ratio by 15.97 and 11%, respectively. This paper addresses something new in the scientific world to date when evaluating in a multi-objective optimization problem, the compression and modal performance of an auxetic reentrant model. Findings It was possible to find multi-objective optimized structures. It was possible to increase the critical buckling load by 42.82%, and the failure load in compression performance by 26.75%. Furthermore, in the optimization of modal performance, it was possible to increase the natural frequency by 37.43%, and decrease the mass by 15.97%. Finally, all 5 responses analyzed simultaneously were optimized. In this case, it was possible to increase the critical buckling load by 42.55%, increase the failure load by 28.70% and reduce the mass and Poisson's ratio by 15.97 and 11%, respectively. Originality/value There is no work in the literature to date that performed the optimization of 5 responses simultaneously of a reentrant hexagonal cell auxetic structure. This paper also presents an unprecedented statistical analysis in the literature that verifies how the design factors impact each of the responses.","PeriodicalId":50522,"journal":{"name":"Engineering Computations","volume":"5 6","pages":"0"},"PeriodicalIF":0.0,"publicationDate":"2023-11-07","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"135431904","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Efficient iterative schemes based on Newton's method and fixed-point iteration for solving nonlinear matrix equation <i>X^p</i> = <i>Q</i>±<i>A</i>(<i>X</i>⁻¹+<i>B</i>)⁻¹<i>A^T</i>","authors":"Raziyeh Erfanifar, Masoud Hajarian","doi":"10.1108/ec-07-2023-0322","DOIUrl":"https://doi.org/10.1108/ec-07-2023-0322","url":null,"abstract":"Purpose In this paper, the authors study the nonlinear matrix equation <m:math xmlns:m=\"http://www.w3.org/1998/Math/MathML\"><m:mrow><m:msup><m:mi>X</m:mi><m:mi>p</m:mi></m:msup><m:mo>=</m:mo><m:mo>Q</m:mo><m:mo>±</m:mo><m:mi>A</m:mi><m:mrow><m:mo>(</m:mo></m:mrow><m:msup><m:mi>X</m:mi><m:mrow><m:mi>-</m:mi><m:mn>1</m:mn></m:mrow></m:msup><m:mo>+</m:mo><m:mi>B</m:mi><m:mrow><m:mo>)</m:mo></m:mrow><m:msup><m:mrow /><m:mrow><m:mi>-</m:mi><m:mn>1</m:mn></m:mrow></m:msup><m:msup><m:mi>A</m:mi><m:mrow><m:mi>T</m:mi></m:mrow></m:msup><m:mo>,</m:mo></m:mrow></m:math> that occurs in many applications such as in filtering, network systems, optimal control and control theory. Design/methodology/approach The authors present some theoretical results for the existence of the solution of this nonlinear matrix equation. Then the authors propose two iterative schemes without inversion to find the solution to the nonlinear matrix equation based on Newton's method and fixed-point iteration. Also the authors show that the proposed iterative schemes converge to the solution of the nonlinear matrix equation, under situations. Findings The efficiency indices of the proposed schemes are presented, and since the initial guesses of the proposed iterative schemes have a high cost, the authors reduce their cost by changing them. Therefore, compared to the previous scheme, the proposed schemes have superior efficiency indices <m:math xmlns:m=\"http://www.w3.org/1998/Math/MathML\"><m:mo>.</m:mo></m:math> Originality/value Finally, the accuracy and effectiveness of the proposed schemes in comparison to an existing scheme are demonstrated by various numerical examples. Moreover, as an application, by using the proposed schemes, the authors can get the optimal controller state feedback of $x(t+1) = A x(t) + C v(t)$.","PeriodicalId":50522,"journal":{"name":"Engineering Computations","volume":"7 6","pages":"0"},"PeriodicalIF":0.0,"publicationDate":"2023-11-07","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"135431103","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Seismic reliability analysis of slope in spatially variable soils using multiple response surfaces","authors":"Chunli Li, Liang Li, Yungming Cheng, Liang Xu, Guangming Yu","doi":"10.1108/ec-04-2023-0165","DOIUrl":"https://doi.org/10.1108/ec-04-2023-0165","url":null,"abstract":"Purpose This paper aims to develop an efficient algorithm combining straightforward response surface functions with Monte Carlo simulation to conduct seismic reliability analysis in a systematical way. Design/methodology/approach The representative slip surfaces are identified and based on to calibrate multiple response surface functions with acceptable accuracy. The calibrated response surfaces are used to determine the yield acceleration in Newmark sliding displacement analysis. Then, the displacement-based limit state function is adopted to conduct seismic reliability analysis. Findings The calibrated response surface functions have fairly good accuracy in predicting the yield acceleration in Newmark sliding displacement analysis. The seismic reliability is influenced by such factors as PGA, spatial variability and threshold value. The proposed methodology serves as an effective tool for geotechnical practitioners. Originality/value The multiple sources of a seismic slope response can be effectively determined using the multiple response surface functions, which are easily implemented within geotechnical engineering.","PeriodicalId":50522,"journal":{"name":"Engineering Computations","volume":"4 5","pages":"0"},"PeriodicalIF":0.0,"publicationDate":"2023-11-07","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"135431734","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Iterated-local-search-based chaotic differential evolution algorithm for hybrid-load part feeding scheduling optimization in mixed-model assembly lines","authors":"Zhu Wang, Hongtao Hu, Tianyu Liu","doi":"10.1108/ec-07-2023-0369","DOIUrl":"https://doi.org/10.1108/ec-07-2023-0369","url":null,"abstract":"Purpose Driven by sustainable production, mobile robots are introduced as a new clean-energy material handling tool for mixed-model assembly lines (MMALs), which reduces energy consumption and lineside inventory of workstations (LSI). Nevertheless, the previous part feeding scheduling method was designed for conventional material handling tools without considering the flexible spatial layout of the robotic mobile fulfillment system (RMFS). To fill this gap, this paper focuses on a greening mobile robot part feeding scheduling problem with Just-In-Time (JIT) considerations, where the layout and number of pods can be adjusted. Design/methodology/approach A novel hybrid-load pod (HL-pod) and mobile robot are proposed to carry out part feeding tasks between material supermarkets and assembly lines. A bi-objective mixed-integer programming model is formulated to minimize both total energy consumption and LSI, aligning with environmental and sustainable JIT goals. Due to the NP-hard nature of the proposed problem, a chaotic differential evolution algorithm for multi-objective optimization based on iterated local search (CDEMIL) algorithm is presented. The effectiveness of the proposed algorithm is verified by dealing with the HL-pod-based greening part feeding scheduling problem in different problem scales and compared to two benchmark algorithms. Managerial insights analyses are conducted to implement the HL-pod strategy. Findings The CDEMIL algorithm's ability to produce Pareto fronts for different problem scales confirms its effectiveness and feasibility. Computational results show that the proposed algorithm outperforms the other two compared algorithms regarding solution quality and convergence speed. Additionally, the results indicate that the HL-pod performs better than adopting a single type of pod. Originality/value This study proposes an innovative solution to the scheduling problem for efficient JIT part feeding using RMFS and HL-pods in automobile MMALs. It considers both the layout and number of pods, ensuring a sustainable and environmental-friendly approach to production.","PeriodicalId":50522,"journal":{"name":"Engineering Computations","volume":"4 s9","pages":"0"},"PeriodicalIF":0.0,"publicationDate":"2023-11-07","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"135431735","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Robust adaptive analysis of dynamic responses of wave energy converters","authors":"Yingguang Wang","doi":"10.1108/ec-02-2023-0062","DOIUrl":"https://doi.org/10.1108/ec-02-2023-0062","url":null,"abstract":"Purpose The purpose of this paper is to exploit a new and robust method to forecast the long-term extreme dynamic responses for wave energy converters (WECs). Design/methodology/approach A new adaptive binned kernel density estimation (KDE) methodology is first proposed in this paper. Findings By examining the calculation results the authors has found that in the tail region the proposed new adaptive binned KDE distribution curve becomes very smooth and fits quite well with the histogram of the measured ocean wave dataset at the National Data Buoy Center (NDBC) station 46,059. Carefully studying the calculation results also reveals that the 50-year extreme power-take-off heaving force value forecasted based on the environmental contour derived using the new method is 3572600N, which is much larger than the value 2709100N forecasted via the Rosenblatt-inverse second-order reliability method (ISORM) contour method. Research limitations/implications The proposed method overcomes the disadvantages of all the existing nonparametric and parametric methods for predicting the tail region probability density values of the sea state parameters. Originality/value It is concluded that the proposed new adaptive binned KDE method is robust and can forecast well the 50-year extreme dynamic responses for WECs.","PeriodicalId":50522,"journal":{"name":"Engineering Computations","volume":"5 s11","pages":"0"},"PeriodicalIF":0.0,"publicationDate":"2023-11-07","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"135431905","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Numerical simulation of mechanical characteristics and tension-torsion coupling effect of tension-type anchor cable","authors":"Xin Meng, Qingyang Ren, Songqiang Xiao, Bin Chen, Hongfei Li","doi":"10.1108/ec-08-2023-0440","DOIUrl":"https://doi.org/10.1108/ec-08-2023-0440","url":null,"abstract":"Purpose The purpose of this paper is to simulate the tension process of tension-type anchor cable and to explore the mechanical characteristics and tension-torsion coupling effect of anchor cable subjected to tension. Design/methodology/approach ABAQUS numerical software is applied to construct the numerical models of tension-type anchor cables with different diameters. Through explicit contact, the characteristics of contact between grouting body-anchor cable and grouting body-rock mass are determined. Confining pressure is applied to the model through surface pressure, and drawing force is applied to the model by displacement loading so as to simulate the tension process of the anchor cable. Findings The results show that the stress is transmitted in both axial and radial directions in the anchorage section and distributed in a cone. The shear stress in the grouting body is unevenly distributed, and its peak value increases with the rise in confining pressure and anchor cable diameter. The stress characteristics of torque and axial force are basically consistent and evenly distributed in the free section; they gradually decrease in the anchorage section. Due to the tension-torsion coupling effect, the internal stress characteristics of the anchor cable structure vary. On average, the anchorage performance of each anchor cable model is improved by 6.19%. Originality/value The proposed method of numerical modelling is effective in addressing the interface contact between the anchor cable and the grouting body and in solving the problem with convergence of calculation. Compared with the indoor test, this method is more suited to collecting the internal mechanical data of the anchor body.","PeriodicalId":50522,"journal":{"name":"Engineering Computations","volume":"95 6","pages":"0"},"PeriodicalIF":0.0,"publicationDate":"2023-11-06","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"135584442","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Comparisons between direct interpolation and reciprocity techniques of the boundary element method for solving two-dimensional Helmholtz problems","authors":"Thiago Galdino Balista, Carlos Friedrich Loeffler, Luciano Lara, Webe João Mansur","doi":"10.1108/ec-06-2023-0290","DOIUrl":"https://doi.org/10.1108/ec-06-2023-0290","url":null,"abstract":"Purpose This work compares the performance of the three boundary element techniques for solving Helmholtz problems: dual reciprocity, multiple reciprocity and direct interpolation. All techniques transform domain integrals into boundary integrals, despite using different principles to reach this purpose. Design/methodology/approach Comparisons here performed include the solution of eigenvalue and response by frequency scanning, analyzing many features that are not comprehensively discussed in the literature, as follows: the type of boundary conditions, suitable number of degrees of freedom, modal content, number of primitives in the multiple reciprocity method (MRM) and the requirement of internal interpolation points in techniques that use radial basis functions as dual reciprocity and direct interpolation. Findings Among the other aspects, this work can conclude that the solution of the eigenvalue and response problems confirmed the reasonable accuracy of the dual reciprocity boundary element method (DRBEM) only for the calculation of the first natural frequencies. Concerning the direct interpolation boundary element method (DIBEM), its interpolation characteristic allows more accessibility for solving more elaborate problems. Despite requiring a greater number of interpolating internal points, the DIBEM has presented higher-quality results for the eigenvalue and response problems. The MRM results were satisfactory in terms of accuracy just for the low range of frequencies; however, the neglected higher-order primitives impact the accuracy of the dynamic response as a whole. Originality/value There are safe alternatives for solving engineering stationary dynamic problems using the boundary element method (BEM), but there are no suitable comparisons between these different techniques. This paper presents the particularities and detailed comparisons approaching the accuracy of the three important BEM techniques, aiming at response and frequency evaluation, which are not found in the specialized literature.","PeriodicalId":50522,"journal":{"name":"Engineering Computations","volume":"85 1","pages":"0"},"PeriodicalIF":0.0,"publicationDate":"2023-11-06","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"135585248","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Numerical simulation of the extrusion process of viscoplastic materials using a radial point interpolation method","authors":"Daniel E.S. Rodrigues, Jorge Belinha, Renato Natal Jorge","doi":"10.1108/ec-01-2023-0044","DOIUrl":"https://doi.org/10.1108/ec-01-2023-0044","url":null,"abstract":"Purpose Fused Filament Fabrication (FFF) is an extrusion-based manufacturing process using fused thermoplastics. Despite its low cost, the FFF is not extensively used in high-value industrial sectors mainly due to parts' anisotropy (related to the deposition strategy) and residual stresses (caused by successive heating cycles). Thus, this study aims to investigate the process improvement and the optimization of the printed parts. Design/methodology/approach In this work, a meshless technique – the Radial Point Interpolation Method (RPIM) – is used to numerically simulate the viscoplastic extrusion process – the initial phase of the FFF. Unlike the FEM, in meshless methods, there is no pre-established relationship between the nodes so the nodal mesh will not face mesh distortions and the discretization can easily be modified by adding or removing nodes from the initial nodal mesh. The accuracy of the obtained results highlights the importance of using meshless techniques in this field. Findings Meshless methods show particular relevance in this topic since the nodes can be distributed to match the layer-by-layer growing condition of the printing process. Originality/value Using the flow formulation combined with the heat transfer formulation presented here for the first time within an in-house RPIM code, an algorithm is proposed, implemented and validated for benchmark examples.","PeriodicalId":50522,"journal":{"name":"Engineering Computations","volume":"84 5","pages":"0"},"PeriodicalIF":0.0,"publicationDate":"2023-11-06","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"135585251","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}