Computer Methods in Applied Mechanics and Engineering最新文献

{"title":"A machine-learning enabled digital-twin framework for next generation precision agriculture and forestry","authors":"","doi":"10.1016/j.cma.2024.117250","DOIUrl":"10.1016/j.cma.2024.117250","url":null,"abstract":"<div><p>This work utilizes the modern synergy between flexible, rapid, simulations and quick assimilation of data in order to develop next-generation tools for precise biomass management of large-scale agricultural and forestry systems. Additionally, when integrated with satellite and drone-based digital elevation technologies, the results lead to digital replicas of physical systems, or so-called digital-twins, which offer a powerful framework by which to optimally manage agricultural and forestry assets. Specifically, this enables the investigation of inverse problems seeking to ascertain ideal parameter combinations, such as the number of plants/trees, plant/tree spacing, light intensity, water availability, soil resources, available planting surface area, initial seedling size, genetic variation, etc. to obtain optimal system performance. Towards this goal, a digital-twin framework is developed, consisting of a rapid computational physics engine to simulate an agricultural installation, containing thousands of growing, interacting, plants/trees. This model is then driven by a machine-learning algorithm to obtain optimal parameter sets that match observed statistical representations of a time series of growing agricultural canopy surfaces, measured by digital elevation models. Model simulations are provided to illustrate the approach and to show how such a tool can be used for large-scale biomass management.</p></div>","PeriodicalId":55222,"journal":{"name":"Computer Methods in Applied Mechanics and Engineering","volume":null,"pages":null},"PeriodicalIF":6.9,"publicationDate":"2024-08-26","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142076277","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":1,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Heteroscedastic Gaussian Process Regression for material structure–property relationship modeling","authors":"","doi":"10.1016/j.cma.2024.117326","DOIUrl":"10.1016/j.cma.2024.117326","url":null,"abstract":"<div><p>Uncertainty quantification is a critical aspect of machine learning models for material property predictions. Gaussian Process Regression (GPR) is a popular technique for capturing uncertainties, but most existing models assume homoscedastic aleatoric uncertainty (noise), which may not adequately represent the heteroscedastic behavior observed in real-world datasets. Heteroscedasticity arises from various factors, such as measurement errors and inherent variability in material properties. Ignoring heteroscedasticity can lead to lower model performance, biased uncertainty estimates, and inaccurate predictions. Existing Heteroscedastic Gaussian Process Regression (HGPR) models often employ complicated structures to capture input-dependent noise but may lack interpretability. In this paper, we propose an HGPR approach that combines GPR with polynomial regression-based noise modeling to capture and quantify uncertainties in material property predictions while providing interpretable noise models. We demonstrate the effectiveness of our approach on both synthetic and physics-based simulation datasets, including mechanical properties (effective stress) of porous materials. We also introduce an approximated expected log predictive density method for model selection, which eliminates the need for retraining the model during leave-one-out cross-validation, allowing for efficient hyperparameter tuning and model evaluation. By capturing heteroscedastic behavior, enhancing interpretability, and improving model selection, our approach contributes to the development of more robust and reliable machine learning models for material property predictions, enabling informed decision-making in material design and optimization.</p></div>","PeriodicalId":55222,"journal":{"name":"Computer Methods in Applied Mechanics and Engineering","volume":null,"pages":null},"PeriodicalIF":6.9,"publicationDate":"2024-08-26","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142076273","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":1,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"An optimally convergent Fictitious Domain method for interface problems","authors":"","doi":"10.1016/j.cma.2024.117327","DOIUrl":"10.1016/j.cma.2024.117327","url":null,"abstract":"<div><p>We introduce a novel Fictitious Domain (FD) unfitted method for interface problems associated with a second-order elliptic linear differential operator, that achieves optimal convergence without the need for adaptive mesh refinements nor enrichments of the Finite Element spaces. The key aspect of the proposed method is that it extends the solution into the fictitious domain in a way that ensures high global regularity. Continuity of the solution across the interface is enforced through a boundary Lagrange multiplier. The subdomains coupling, however, is not achieved by means of the duality pairing with the Lagrange multiplier, but through an <span><math><msup><mrow><mi>L</mi></mrow><mrow><mn>2</mn></mrow></msup></math></span> product with the <span><math><msup><mrow><mi>H</mi></mrow><mrow><mn>1</mn></mrow></msup></math></span> Riesz representative of the latter, thus avoiding gradient jumps across the interface. Thanks to the enhanced regularity, the proposed method attains an increase, with respect to standard FD methods, of up to one order of convergence in energy norm. The Finite Element formulation of the method is presented, followed by its analysis. Numerical tests on a model problem demonstrate its effectiveness and its superior accuracy compared to standard unfitted methods.</p></div>","PeriodicalId":55222,"journal":{"name":"Computer Methods in Applied Mechanics and Engineering","volume":null,"pages":null},"PeriodicalIF":6.9,"publicationDate":"2024-08-26","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://www.sciencedirect.com/science/article/pii/S0045782524005826/pdfft?md5=76bbb3321bd111da71ce05588f536f5b&pid=1-s2.0-S0045782524005826-main.pdf","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142076271","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":1,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Shape optimization of non-matching isogeometric shells with moving intersections","authors":"","doi":"10.1016/j.cma.2024.117322","DOIUrl":"10.1016/j.cma.2024.117322","url":null,"abstract":"<div><p>While shape optimization using isogeometric shells exhibits appealing features by integrating design geometries and analysis models, challenges arise when addressing computer-aided design (CAD) geometries comprised of multiple non-uniform rational B-splines (NURBS) patches, which are common in practice. The intractability stems from surface intersections within these CAD models. In this paper, we develop an approach for shape optimization of non-matching isogeometric shells incorporating intersection movement. Separately parametrized NURBS surfaces are modeled using Kirchhoff–Love shell theory and coupled using a penalty-based formulation. The optimization scheme allows shell patches to move without preserving relative location with other members during the shape optimization. This flexibility is achieved through an implicit state function, and analytical sensitivities are derived for the relative movement of shell patches. The introduction of differentiable intersections expands the design space and overcomes challenges associated with large mesh distortion, particularly when optimal shapes involve significant movement of patch intersections in physical space. Throughout optimization iterations, all members within the shell structures maintain the NURBS geometry representation, enabling efficient integration of analysis and design models. The optimization approach leverages the multilevel design concept by selecting a refined model for accurate analysis from a coarse design model while maintaining the same geometry. We adopt several example problems to verify the effectiveness of the proposed scheme and demonstrate its applicability to the optimization of the internal stiffeners of an aircraft wing.</p></div>","PeriodicalId":55222,"journal":{"name":"Computer Methods in Applied Mechanics and Engineering","volume":null,"pages":null},"PeriodicalIF":6.9,"publicationDate":"2024-08-26","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://www.sciencedirect.com/science/article/pii/S0045782524005772/pdfft?md5=995c0a51bc4050740a320a45c51adef2&pid=1-s2.0-S0045782524005772-main.pdf","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142076276","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":1,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Incorporating interface effects into multi-material topology optimization by improving interface configuration: An energy-based approach","authors":"","doi":"10.1016/j.cma.2024.117325","DOIUrl":"10.1016/j.cma.2024.117325","url":null,"abstract":"<div><p>Interfaces between structural multi-materials generally exhibit asymmetric resistance to tension and compression. Given this interface behavior, this work suggests an energy-based approach to improve the interface configuration for multi-material topology optimization. Based on the strain spectral decomposition, we decompose the structural elastic strain energy into tensile and compressive portions. In the density-based topology optimization framework, we use the gradient-based method to track the interface between multiple materials. Then, we construct an interface-associated scalar field to penalize the tensile portion of the strain energy, causing a pseudo-degradation of the strain energy at the interface region. Finally, within limited material usages and by minimizing the linear weighted structural strain energy and its pseudo-degradation, multi-material topology optimization with improved interface configuration is achieved. Several 2D and 3D numerical examples are investigated, by which the effectiveness and robustness of the suggested approach are fairly validated.</p></div>","PeriodicalId":55222,"journal":{"name":"Computer Methods in Applied Mechanics and Engineering","volume":null,"pages":null},"PeriodicalIF":6.9,"publicationDate":"2024-08-24","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142049702","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":1,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Discontinuous Galerkin schemes for hyperbolic systems in non-conservative variables: Quasi-conservative formulation with subcell finite volume corrections","authors":"","doi":"10.1016/j.cma.2024.117311","DOIUrl":"10.1016/j.cma.2024.117311","url":null,"abstract":"<div><p>We present a novel quasi-conservative arbitrary high order accurate ADER (Arbitrary-Derivative) discontinuous Galerkin method allowing to efficiently use a non-conservative form of the considered partial differential system, so that the governing equations can be solved directly in the most physically relevant set of variables. This is particularly interesting for multi-material flows with moving interfaces and steep, large magnitude contact discontinuities, as well as in presence of highly non-linear thermodynamics. However, the non-conservative formulation of course introduces a conservation error which would normally lead to a wrong approximation of shock waves. Hence, from the theoretical point of view, we give a formal definition of the conservation defect of non-conservative schemes and we analyze this defect providing a local quasi-conservation condition, which allows us to prove a <em>modified Lax–Wendroff theorem</em>. Within this formalism, we also reformulate classical results concerning smooth solutions, contact discontinuities and moving interfaces. Then, to deal with shock waves in practice, we exploit the framework of the so-called <em>a posteriori</em> subcell finite volume (FV) limiter, so that, in troubled cells appropriately detected, we can incorporate a local conservation correction. Our corrected FV update entirely removes the local conservation defect, allowing, at least formally, to fit in the hypotheses of the proposed modified Lax–Wendroff theorem. Here, the shock-triggered troubled cells are detected by combining physical admissibility criteria, a discrete maximum principle and a shock sensor inspired by Lagrangian hydrodynamics.</p><p>To prove the capabilities of our novel approach, first, we show that we are able to recover the same results given by conservative schemes on classical benchmarks for the single-fluid Euler equations. We then conclude the presentation by showing the improved reliability of our scheme on the <em>multi-fluid Euler system</em> on examples like the interaction of a shock with a helium bubble for which we are able to avoid the development of any spurious oscillations.</p></div>","PeriodicalId":55222,"journal":{"name":"Computer Methods in Applied Mechanics and Engineering","volume":null,"pages":null},"PeriodicalIF":6.9,"publicationDate":"2024-08-24","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://www.sciencedirect.com/science/article/pii/S004578252400567X/pdfft?md5=d0858bec5c60dcd6967539ab1f864eb2&pid=1-s2.0-S004578252400567X-main.pdf","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142058395","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":1,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"A conforming frictional beam contact model","authors":"","doi":"10.1016/j.cma.2024.117310","DOIUrl":"10.1016/j.cma.2024.117310","url":null,"abstract":"<div><p>We develop a model for predicting the mechanical behavior of a system of slender one-dimensional bodies (fibers or beams) interacting via frictional contact. Relying on an integral penalty-based formulation, it can robustly capture the behavior in the case of conforming contact occurring over regions of finite size. Two formulations of the model are presented and validated against fully resolved continuum finite element simulations. Overall, the proposed framework is an effective tool in exploring the mechanical behavior of fabrics, textiles as well as three-dimensional frictional architected solids, as demonstrated by the simulation of the effective response of a periodic intertwined metamaterial.</p></div>","PeriodicalId":55222,"journal":{"name":"Computer Methods in Applied Mechanics and Engineering","volume":null,"pages":null},"PeriodicalIF":6.9,"publicationDate":"2024-08-23","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://www.sciencedirect.com/science/article/pii/S0045782524005668/pdfft?md5=68c10bcc21ee2e4de9723ebe33ca95d0&pid=1-s2.0-S0045782524005668-main.pdf","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142044826","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":1,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Stress-related discrete variable topology optimization with handling non-physical stress concentrations","authors":"","doi":"10.1016/j.cma.2024.117293","DOIUrl":"10.1016/j.cma.2024.117293","url":null,"abstract":"<div><p>The accuracy of stress calculation with a fixed mesh significantly affects the stress-based topology optimization, due to potential non-physical stress concentrations in voxel-based topology descriptions. This paper proposes a novel problem-independent machine learning enhanced high-precision stress calculation method (PIML-HPSCM) to address this challenge. As an immersed analysis method, PIML-HPSCM combines the high efficiency of fixed mesh with the accuracy of body-fitted mesh, without complex integration schemes of other immersed methods. PIML-HPSCM utilizes the extended multiscale finite element method to depict the material heterogeneity within high-resolution boundary elements. The accurate stress field can then be calculated conveniently by establishing stress evaluation matrices of high-resolution boundary elements. Moreover, the PIML is independent of problem settings and is applicable for various problems with the same governing equation type. Invoking the offline-trained neural network online can enhance stress calculation efficiency by 10–20 times. The stress-based discrete variable topology optimization, which naturally avoids singular stress phenomenon, is efficiently addressed by the sequential approximate integer programming method with PIML-HPSCM. Results from 2D and 3D examples demonstrate that the stresses calculated by PIML-HPSCM are consistent with those by body-fitted mesh, and optimized designs effectively eliminate stress concentrations of initial designs and have uniform stress distributions.</p></div>","PeriodicalId":55222,"journal":{"name":"Computer Methods in Applied Mechanics and Engineering","volume":null,"pages":null},"PeriodicalIF":6.9,"publicationDate":"2024-08-23","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142044825","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":1,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Model order reduction of time-domain acoustic finite element simulations with perfectly matched layers","authors":"","doi":"10.1016/j.cma.2024.117298","DOIUrl":"10.1016/j.cma.2024.117298","url":null,"abstract":"<div><p>This paper presents a stability-preserving model reduction method for an acoustic finite element model with perfectly matched layers (PMLs). PMLs are often introduced into an unbounded domain to simulate the Sommerfeld radiation condition. These layers act as anisotropic damping materials to absorb the scattered field, of which the material properties are frequency- and coordinate-dependent. The corresponding time-domain model size is often very large due to this frequency-dependent property and the number of elements needed per wavelength. Therefore, to enable efficient transient simulations, this paper proposes a two-step method to generate stable reduced order models (ROMs) of such systems. Firstly, the modified and stable version of PMLs is projected by a one-sided split basis, which gives a stable intermediate ROM. Secondly, the intermediate ROM is modified to satisfy the stability-preserving condition by applying the modal transformation. Applying any one-sided model order reduction method on this modified model leads to a stable and small ROM. This two-step method is further extended to account for the locally-conformal PML model by reformulating it in curvilinear coordinates, which works for arbitrary convex truncated domains. The proposed method is successfully verified by several numerical simulations.</p></div>","PeriodicalId":55222,"journal":{"name":"Computer Methods in Applied Mechanics and Engineering","volume":null,"pages":null},"PeriodicalIF":6.9,"publicationDate":"2024-08-23","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142045835","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":1,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"A formulation for fluid–structure interaction problems with immersed flexible solids: Application to splitters subjected to flow past cylinders with different cross-sections","authors":"","doi":"10.1016/j.cma.2024.117306","DOIUrl":"10.1016/j.cma.2024.117306","url":null,"abstract":"<div><p>In the finite element method framework, a fluid–structure formulation is developed by coupling an Eulerian fixed-mesh fluid approach with a Lagrangian deforming-mesh description for a flexible solid. The coupled formulation is solved using a staggered scheme during time. For the fluid solution stage, the solid walls are considered as a time-variable internal boundary. The velocity and pressure fields are obtained by solving the weak form of the fluid dynamic equations in which the solid velocity is imposed on the internal boundary via a penalization term. For the solid solution stage, the displacement field is obtained by solving the discrete solid dynamic equations which consider traction forces computed by integrating pressures and viscous stresses on the nodes belonging to the solid walls. This novel technique is firstly applied to analyze a flexible splitter under the shedding of a flow past square cylinder due to this problem is considered as a benchmark in the literature. The present solutions agree with those computed using body-fitted techniques, thus validating the proposal. Secondly, flexible splitter motions under the shedding of flow past cylinders with different cross-sections and splitter lengths are comprehensively studied. Overall, the computed results confirmed that the hydrodynamic coefficients on the cylinders were reduced because of the presence of the splitter.</p></div>","PeriodicalId":55222,"journal":{"name":"Computer Methods in Applied Mechanics and Engineering","volume":null,"pages":null},"PeriodicalIF":6.9,"publicationDate":"2024-08-23","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142044827","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":1,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}