Pub Date : 2021-11-26DOI: 10.1186/s40323-021-00211-7
A. Pasquale, A. Ammar, A. Falcó, S. Perotto, E. Cueto, J. Duval, F. Chinesta
{"title":"A separated representation involving multiple time scales within the Proper Generalized Decomposition framework","authors":"A. Pasquale, A. Ammar, A. Falcó, S. Perotto, E. Cueto, J. Duval, F. Chinesta","doi":"10.1186/s40323-021-00211-7","DOIUrl":"https://doi.org/10.1186/s40323-021-00211-7","url":null,"abstract":"","PeriodicalId":37424,"journal":{"name":"Advanced Modeling and Simulation in Engineering Sciences","volume":null,"pages":null},"PeriodicalIF":0.0,"publicationDate":"2021-11-26","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"65854082","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2021-11-03DOI: 10.1186/s40323-021-00210-8
de Gooijer, Boukje M., Havinga, Jos, Geijselaers, Hubert J. M., van den Boogaard, Anton H.
Surrogate modelling is a powerful tool to replace computationally expensive nonlinear numerical simulations, with fast representations thereof, for inverse analysis, model-based control or optimization. For some problems, it is required that the surrogate model describes a complete output field. To construct such surrogate models, proper orthogonal decomposition (POD) can be used to reduce the dimensionality of the output data. The accuracy of the surrogate models strongly depends on the (pre)processing actions that are used to prepare the data for the dimensionality reduction. In this work, POD-based surrogate models with Radial Basis Function interpolation are used to model high-dimensional FE data fields. The effect of (pre)processing methods on the accuracy of the result field is systematically investigated. Different existing methods for surrogate model construction are compared with a novel method. Special attention is given to data fields consisting of several physical meanings, e.g. displacement, strain and stress. A distinction is made between the errors due to truncation and due to interpolation of the data. It is found that scaling the data per physical part substantially increases the accuracy of the surrogate model.
{"title":"Evaluation of POD based surrogate models of fields resulting from nonlinear FEM simulations","authors":"de Gooijer, Boukje M., Havinga, Jos, Geijselaers, Hubert J. M., van den Boogaard, Anton H.","doi":"10.1186/s40323-021-00210-8","DOIUrl":"https://doi.org/10.1186/s40323-021-00210-8","url":null,"abstract":"Surrogate modelling is a powerful tool to replace computationally expensive nonlinear numerical simulations, with fast representations thereof, for inverse analysis, model-based control or optimization. For some problems, it is required that the surrogate model describes a complete output field. To construct such surrogate models, proper orthogonal decomposition (POD) can be used to reduce the dimensionality of the output data. The accuracy of the surrogate models strongly depends on the (pre)processing actions that are used to prepare the data for the dimensionality reduction. In this work, POD-based surrogate models with Radial Basis Function interpolation are used to model high-dimensional FE data fields. The effect of (pre)processing methods on the accuracy of the result field is systematically investigated. Different existing methods for surrogate model construction are compared with a novel method. Special attention is given to data fields consisting of several physical meanings, e.g. displacement, strain and stress. A distinction is made between the errors due to truncation and due to interpolation of the data. It is found that scaling the data per physical part substantially increases the accuracy of the surrogate model.","PeriodicalId":37424,"journal":{"name":"Advanced Modeling and Simulation in Engineering Sciences","volume":null,"pages":null},"PeriodicalIF":0.0,"publicationDate":"2021-11-03","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"138510033","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2021-10-30DOI: 10.21203/rs.3.rs-1050987/v1
Jhouben Cuesta-Ramirez, R. Riche, O. Roustant, G. Perrin, Cédric Durantin, A. Glière
Most real optimization problems are defined over a mixed search space where the variables are both discrete and continuous. In engineering applications, the objective function is typically calculated with a numerically costly black-box simulation. General mixed and costly optimization problems are therefore of a great practical interest, yet their resolution remains in a large part an open scientific question. In this article, costly mixed problems are approached through Gaussian processes where the discrete variables are relaxed into continuous latent variables. The continuous space is more easily harvested by classical Bayesian optimization techniques than a mixed space would. Discrete variables are recovered either subsequently to the continuous optimization, or simultaneously with an additional continuous-discrete compatibility constraint that is handled with augmented Lagrangians. Several possible implementations of such Bayesian mixed optimizers are compared. In particular, the reformulation of the problem with continuous latent variables is put in competition with searches working directly in the mixed space. Among the algorithms involving latent variables and an augmented Lagrangian, a particular attention is devoted to the Lagrange multipliers for which a local and a global estimation techniques are studied. The comparisons are based on the repeated optimization of three analytical functions and a beam design problem.
{"title":"A comparison of mixed-variables Bayesian optimization approaches","authors":"Jhouben Cuesta-Ramirez, R. Riche, O. Roustant, G. Perrin, Cédric Durantin, A. Glière","doi":"10.21203/rs.3.rs-1050987/v1","DOIUrl":"https://doi.org/10.21203/rs.3.rs-1050987/v1","url":null,"abstract":"Most real optimization problems are defined over a mixed search space where the variables are both discrete and continuous. In engineering applications, the objective function is typically calculated with a numerically costly black-box simulation. General mixed and costly optimization problems are therefore of a great practical interest, yet their resolution remains in a large part an open scientific question. In this article, costly mixed problems are approached through Gaussian processes where the discrete variables are relaxed into continuous latent variables. The continuous space is more easily harvested by classical Bayesian optimization techniques than a mixed space would. Discrete variables are recovered either subsequently to the continuous optimization, or simultaneously with an additional continuous-discrete compatibility constraint that is handled with augmented Lagrangians. Several possible implementations of such Bayesian mixed optimizers are compared. In particular, the reformulation of the problem with continuous latent variables is put in competition with searches working directly in the mixed space. Among the algorithms involving latent variables and an augmented Lagrangian, a particular attention is devoted to the Lagrange multipliers for which a local and a global estimation techniques are studied. The comparisons are based on the repeated optimization of three analytical functions and a beam design problem.","PeriodicalId":37424,"journal":{"name":"Advanced Modeling and Simulation in Engineering Sciences","volume":null,"pages":null},"PeriodicalIF":0.0,"publicationDate":"2021-10-30","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"44351456","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2021-10-04DOI: 10.1186/s40323-021-00207-3
Barabinot, Philippe, Scanff, Ronan, Ladevèze, Pierre, Néron, David, Cauville, Bruno
Digital Twins, which tend to intervene over the entire life cycle of products from early design phase to predictive maintenance through optimization processes, are increasingly emerging as an essential component in the future of industries. To reduce the computational time reduced-order modeling (ROM) methods can be useful. However, the spread of ROM methods at an industrial level is currently hampered by the difficulty of introducing them into commercial finite element software, due to the strong intrusiveness of the associated algorithms, preventing from getting robust and reliable tools all integrated in a certified product. This work tries to circumvent this issue by introducing a weakly-invasive reformulation of the LATIN-PGD method which is intended to be directly embedded into Simcenter Samcef $$^{hbox {TM}}$$ finite element software. The originality of this approach lies in the remarkably general way of doing, allowing PGD method to deal with not only a particular application but with all facilities already included in such softwares—any non-linearities, any element types, any boundary conditions...—and thus providing a new high-performance all-inclusive non-linear solver.
{"title":"Industrial Digital Twins based on the non-linear LATIN-PGD","authors":"Barabinot, Philippe, Scanff, Ronan, Ladevèze, Pierre, Néron, David, Cauville, Bruno","doi":"10.1186/s40323-021-00207-3","DOIUrl":"https://doi.org/10.1186/s40323-021-00207-3","url":null,"abstract":"Digital Twins, which tend to intervene over the entire life cycle of products from early design phase to predictive maintenance through optimization processes, are increasingly emerging as an essential component in the future of industries. To reduce the computational time reduced-order modeling (ROM) methods can be useful. However, the spread of ROM methods at an industrial level is currently hampered by the difficulty of introducing them into commercial finite element software, due to the strong intrusiveness of the associated algorithms, preventing from getting robust and reliable tools all integrated in a certified product. This work tries to circumvent this issue by introducing a weakly-invasive reformulation of the LATIN-PGD method which is intended to be directly embedded into Simcenter Samcef $$^{hbox {TM}}$$ finite element software. The originality of this approach lies in the remarkably general way of doing, allowing PGD method to deal with not only a particular application but with all facilities already included in such softwares—any non-linearities, any element types, any boundary conditions...—and thus providing a new high-performance all-inclusive non-linear solver.","PeriodicalId":37424,"journal":{"name":"Advanced Modeling and Simulation in Engineering Sciences","volume":null,"pages":null},"PeriodicalIF":0.0,"publicationDate":"2021-10-04","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"138510037","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2021-10-04DOI: 10.1186/s40323-021-00208-2
Kazemzadeh-Parsi, Mohammad Javad, Ammar, Amine, Duval, Jean Louis, Chinesta, Francisco
Space separation within the Proper Generalized Decomposition—PGD—rationale allows solving high dimensional problems as a sequence of lower dimensional ones. In our former works, different geometrical transformations were proposed for addressing complex shapes and spatially non-separable domains. Efficient implementation of separated representations needs expressing the domain as a product of characteristic functions involving the different space coordinates. In the case of complex shapes, more sophisticated geometrical transformations are needed to map the complex physical domain into a regular one where computations are performed. This paper aims at proposing a very efficient route for accomplishing such space separation. A NURBS-based geometry representation, usual in computer aided design—CAD—, is retained and combined with a fully separated representation for allying efficiency (ensured by the fully separated representations) and generality (by addressing complex geometries). Some numerical examples are considered to prove the potential of the proposed methodology.
{"title":"Enhanced parametric shape descriptions in PGD-based space separated representations","authors":"Kazemzadeh-Parsi, Mohammad Javad, Ammar, Amine, Duval, Jean Louis, Chinesta, Francisco","doi":"10.1186/s40323-021-00208-2","DOIUrl":"https://doi.org/10.1186/s40323-021-00208-2","url":null,"abstract":"Space separation within the Proper Generalized Decomposition—PGD—rationale allows solving high dimensional problems as a sequence of lower dimensional ones. In our former works, different geometrical transformations were proposed for addressing complex shapes and spatially non-separable domains. Efficient implementation of separated representations needs expressing the domain as a product of characteristic functions involving the different space coordinates. In the case of complex shapes, more sophisticated geometrical transformations are needed to map the complex physical domain into a regular one where computations are performed. This paper aims at proposing a very efficient route for accomplishing such space separation. A NURBS-based geometry representation, usual in computer aided design—CAD—, is retained and combined with a fully separated representation for allying efficiency (ensured by the fully separated representations) and generality (by addressing complex geometries). Some numerical examples are considered to prove the potential of the proposed methodology.","PeriodicalId":37424,"journal":{"name":"Advanced Modeling and Simulation in Engineering Sciences","volume":null,"pages":null},"PeriodicalIF":0.0,"publicationDate":"2021-10-04","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"138510048","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2021-09-28DOI: 10.1186/s40323-021-00206-4
Jessica Manganotti, F. Caforio, François Kimmig, P. Moireau, S. Imperiale
{"title":"Coupling reduced-order blood flow and cardiac models through energy-consistent strategies: modeling and discretization","authors":"Jessica Manganotti, F. Caforio, François Kimmig, P. Moireau, S. Imperiale","doi":"10.1186/s40323-021-00206-4","DOIUrl":"https://doi.org/10.1186/s40323-021-00206-4","url":null,"abstract":"","PeriodicalId":37424,"journal":{"name":"Advanced Modeling and Simulation in Engineering Sciences","volume":null,"pages":null},"PeriodicalIF":0.0,"publicationDate":"2021-09-28","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"65853475","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2021-09-21DOI: 10.1186/s40323-021-00205-5
Ghnatios, Chady, Barasinski, Anais
A nonparametric method assessing the error and variability margins in solutions depicted in a separated form using experimental results is illustrated in this work. The method assess the total variability of the solution including the modeling error and the truncation error when experimental results are available. The illustrated method is based on the use of the PGD separated form solutions, enriched by transforming a part of the PGD basis vectors into probabilistic one. The constructed probabilistic vectors are restricted to the physical solution’s Stiefel manifold. The result is a real-time parametric PGD solution enhanced with the solution variability and the confidence intervals.
{"title":"A nonparametric probabilistic method to enhance PGD solutions with data-driven approach, application to the automated tape placement process","authors":"Ghnatios, Chady, Barasinski, Anais","doi":"10.1186/s40323-021-00205-5","DOIUrl":"https://doi.org/10.1186/s40323-021-00205-5","url":null,"abstract":"A nonparametric method assessing the error and variability margins in solutions depicted in a separated form using experimental results is illustrated in this work. The method assess the total variability of the solution including the modeling error and the truncation error when experimental results are available. The illustrated method is based on the use of the PGD separated form solutions, enriched by transforming a part of the PGD basis vectors into probabilistic one. The constructed probabilistic vectors are restricted to the physical solution’s Stiefel manifold. The result is a real-time parametric PGD solution enhanced with the solution variability and the confidence intervals.","PeriodicalId":37424,"journal":{"name":"Advanced Modeling and Simulation in Engineering Sciences","volume":null,"pages":null},"PeriodicalIF":0.0,"publicationDate":"2021-09-21","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"138510034","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2021-09-04DOI: 10.1186/s40323-021-00203-7
Schein, Alexander, Gee, Michael W.
This work proposes a framework for projection-based model order reduction (MOR) of computational models aiming at a mechanical analysis of abdominal aortic aneurysms (AAAs). The underlying full-order model (FOM) is patient-specific, stationary and nonlinear. The quantities of interest are the von Mises stress and the von Mises strain field in the AAA wall, which result from loading the structure to the level of diastolic blood pressure at a fixed, imaged geometry (prestressing stage) and subsequent loading to the level of systolic blood pressure with associated deformation of the structure (deformation stage). Prestressing is performed with the modified updated Lagrangian formulation (MULF) approach. The proposed framework aims at a reduction of the computational cost in a many-query context resulting from model uncertainties in two material and one geometric parameter. We apply projection-based MOR to the MULF prestressing stage, which has not been presented to date. Additionally, we propose a reduced-order basis construction technique combining the concept of subspace angles and greedy maximin distance sampling. To further achieve computational speedup, the reduced-order model (ROM) is equipped with the energy-conserving mesh sampling and weighting hyper reduction method. Accuracy of the ROM is numerically tested in terms of the quantities of interest within given bounds of the parameter domain and performance of the proposed ROM in the many-query context is demonstrated by comparing ROM and FOM statistics built from Monte Carlo sampling for three different patient-specific AAAs.
{"title":"Greedy maximin distance sampling based model order reduction of prestressed and parametrized abdominal aortic aneurysms","authors":"Schein, Alexander, Gee, Michael W.","doi":"10.1186/s40323-021-00203-7","DOIUrl":"https://doi.org/10.1186/s40323-021-00203-7","url":null,"abstract":"This work proposes a framework for projection-based model order reduction (MOR) of computational models aiming at a mechanical analysis of abdominal aortic aneurysms (AAAs). The underlying full-order model (FOM) is patient-specific, stationary and nonlinear. The quantities of interest are the von Mises stress and the von Mises strain field in the AAA wall, which result from loading the structure to the level of diastolic blood pressure at a fixed, imaged geometry (prestressing stage) and subsequent loading to the level of systolic blood pressure with associated deformation of the structure (deformation stage). Prestressing is performed with the modified updated Lagrangian formulation (MULF) approach. The proposed framework aims at a reduction of the computational cost in a many-query context resulting from model uncertainties in two material and one geometric parameter. We apply projection-based MOR to the MULF prestressing stage, which has not been presented to date. Additionally, we propose a reduced-order basis construction technique combining the concept of subspace angles and greedy maximin distance sampling. To further achieve computational speedup, the reduced-order model (ROM) is equipped with the energy-conserving mesh sampling and weighting hyper reduction method. Accuracy of the ROM is numerically tested in terms of the quantities of interest within given bounds of the parameter domain and performance of the proposed ROM in the many-query context is demonstrated by comparing ROM and FOM statistics built from Monte Carlo sampling for three different patient-specific AAAs.","PeriodicalId":37424,"journal":{"name":"Advanced Modeling and Simulation in Engineering Sciences","volume":null,"pages":null},"PeriodicalIF":0.0,"publicationDate":"2021-09-04","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"138510049","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2021-09-04DOI: 10.1186/s40323-021-00204-6
M. Zahid, Khalid S. Syed
{"title":"Investigation of pollutants formation in a diesel engine using numerical simulation","authors":"M. Zahid, Khalid S. Syed","doi":"10.1186/s40323-021-00204-6","DOIUrl":"https://doi.org/10.1186/s40323-021-00204-6","url":null,"abstract":"","PeriodicalId":37424,"journal":{"name":"Advanced Modeling and Simulation in Engineering Sciences","volume":null,"pages":null},"PeriodicalIF":0.0,"publicationDate":"2021-09-04","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"42203236","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2021-08-04DOI: 10.1186/s40323-021-00202-8
Chevreuil, Mathilde, Slama, Myriam
The paper deals with approximations of periodic functions that play a significant role in harmonic analysis. The approach revisits the trigonometric polynomials, seen as combinations of functions, and proposes to extend the class of models of the combined functions to a wider class of functions. The key here is to use structured functions, that have low complexity, with suitable functional representation and adapted parametrizations for the approximation. Such representation enables to approximate multivariate functions with few eventually random samples. The new parametrization is determined automatically with a greedy procedure, and a low rank format is used for the approximation associated with each new parametrization. A supervised learning algorithm is used for the approximation of a function of multiple random variables in tree-based tensor format, here the particular Tensor Train format. Adaptive strategies using statistical error estimates are proposed for the selection of the underlying tensor bases and the ranks for the Tensor-Train format. The method is applied for the estimation of the wall pressure for a flow over a cylinder for a range of low to medium Reynolds numbers for which we observe two flow regimes: a laminar flow with periodic vortex shedding and a laminar boundary layer with a turbulent wake (sub-critic regime). The automatic re-parametrization enables here to take into account the specific periodic feature of the pressure.
{"title":"A generalized Fourier transform by means of change of variables within multilinear approximation","authors":"Chevreuil, Mathilde, Slama, Myriam","doi":"10.1186/s40323-021-00202-8","DOIUrl":"https://doi.org/10.1186/s40323-021-00202-8","url":null,"abstract":"The paper deals with approximations of periodic functions that play a significant role in harmonic analysis. The approach revisits the trigonometric polynomials, seen as combinations of functions, and proposes to extend the class of models of the combined functions to a wider class of functions. The key here is to use structured functions, that have low complexity, with suitable functional representation and adapted parametrizations for the approximation. Such representation enables to approximate multivariate functions with few eventually random samples. The new parametrization is determined automatically with a greedy procedure, and a low rank format is used for the approximation associated with each new parametrization. A supervised learning algorithm is used for the approximation of a function of multiple random variables in tree-based tensor format, here the particular Tensor Train format. Adaptive strategies using statistical error estimates are proposed for the selection of the underlying tensor bases and the ranks for the Tensor-Train format. The method is applied for the estimation of the wall pressure for a flow over a cylinder for a range of low to medium Reynolds numbers for which we observe two flow regimes: a laminar flow with periodic vortex shedding and a laminar boundary layer with a turbulent wake (sub-critic regime). The automatic re-parametrization enables here to take into account the specific periodic feature of the pressure.","PeriodicalId":37424,"journal":{"name":"Advanced Modeling and Simulation in Engineering Sciences","volume":null,"pages":null},"PeriodicalIF":0.0,"publicationDate":"2021-08-04","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"138510050","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
京公网安备 11010802042870号
京ICP备2023020795号-1
扫码关注我们
求助内容：
应助结果提醒方式：