Pub Date : 2026-01-01Epub Date: 2025-11-23DOI: 10.1016/j.advengsoft.2025.104072
Yidong Zhao , Xuan Li , Chenfanfu Jiang , Jinhyun Choo
The material point method (MPM), a hybrid Lagrangian–Eulerian particle method, is increasingly used to simulate large-deformation and history-dependent behavior of geomaterials. While explicit time integration dominates current MPM implementations due to its algorithmic simplicity, such schemes are unsuitable for quasi-static and long-term processes typical in geomechanics. Implicit MPM formulations are free of these limitations but remain less adopted, largely due to the difficulty of computing the Jacobian matrix required for Newton-type solvers, especially when consistent tangent operators should be derived for complex constitutive models. In this paper, we introduce GeoWarp—an implicit MPM framework for geomechanics built on NVIDIA Warp—that exploits GPU parallelism and reverse-mode automatic differentiation to compute Jacobians without manual derivation. To enhance efficiency, we develop a sparse Jacobian construction algorithm that leverages the localized particle–grid interactions intrinsic to MPM. The framework is verified through forward and inverse examples in large-deformation elastoplasticity and coupled poromechanics. Results demonstrate that GeoWarp provides a robust, scalable, and extensible platform for differentiable implicit MPM simulation in computational geomechanics.
{"title":"GeoWarp: An automatically differentiable and GPU-accelerated implicit MPM framework for geomechanics based on NVIDIA Warp","authors":"Yidong Zhao , Xuan Li , Chenfanfu Jiang , Jinhyun Choo","doi":"10.1016/j.advengsoft.2025.104072","DOIUrl":"10.1016/j.advengsoft.2025.104072","url":null,"abstract":"<div><div>The material point method (MPM), a hybrid Lagrangian–Eulerian particle method, is increasingly used to simulate large-deformation and history-dependent behavior of geomaterials. While explicit time integration dominates current MPM implementations due to its algorithmic simplicity, such schemes are unsuitable for quasi-static and long-term processes typical in geomechanics. Implicit MPM formulations are free of these limitations but remain less adopted, largely due to the difficulty of computing the Jacobian matrix required for Newton-type solvers, especially when consistent tangent operators should be derived for complex constitutive models. In this paper, we introduce GeoWarp—an implicit MPM framework for geomechanics built on NVIDIA Warp—that exploits GPU parallelism and reverse-mode automatic differentiation to compute Jacobians without manual derivation. To enhance efficiency, we develop a sparse Jacobian construction algorithm that leverages the localized particle–grid interactions intrinsic to MPM. The framework is verified through forward and inverse examples in large-deformation elastoplasticity and coupled poromechanics. Results demonstrate that GeoWarp provides a robust, scalable, and extensible platform for differentiable implicit MPM simulation in computational geomechanics.</div></div>","PeriodicalId":50866,"journal":{"name":"Advances in Engineering Software","volume":"212 ","pages":"Article 104072"},"PeriodicalIF":5.7,"publicationDate":"2026-01-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145618142","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2026-01-01Epub Date: 2025-10-30DOI: 10.1016/j.advengsoft.2025.104051
Taegeun Kim , Hyoungsoo Ko , Jaewon Jang , Sejin Kim , Sunyoung Park , Jae Myung Choe , Young-gu Kim , Dae Sin Kim , Sangseung Lee , Donghyun You
An efficient multiscale coupling method is proposed for simulations of reactor-scale chemical vapor deposition (CVD) with microstructural features. Reactor-scale and microstructure-resolved feature-scale models are coupled through an effective reaction rate formalism, enabling high-resolution deposition simulations while significantly reducing computational cost. A parameterized microstructural model is introduced, in which the relationship between the effective reaction rate and local species consumption rates in the reactor-scale model is directly mapped using precomputed Monte Carlo simulation data. This eliminates the need for iterative calculations or direct numerical simulations of the surface reaction across all the discretized grid points on the wafer, ensuring predictive accuracy while enhancing computational efficiency. Furthermore, an adaptive time-stepping method is developed, dynamically adjusting the time-step size for the feature-scale model based on variations in the effective reaction rate. Through this approach, simulation time is reduced by more than one-third compared to conventional fixed time-step methods, while preserving the accuracy of the effective reaction rate model. The proposed method enables practical and scalable multiscale CVD simulations applicable to industrial reactor design and process optimization, establishing a computationally efficient strategy for integrating reactor-scale and microstructure-resolved feature-scale models.
{"title":"An efficient multiscale coupling method for simulations of reactor-scale chemical vapor deposition with microstructural features","authors":"Taegeun Kim , Hyoungsoo Ko , Jaewon Jang , Sejin Kim , Sunyoung Park , Jae Myung Choe , Young-gu Kim , Dae Sin Kim , Sangseung Lee , Donghyun You","doi":"10.1016/j.advengsoft.2025.104051","DOIUrl":"10.1016/j.advengsoft.2025.104051","url":null,"abstract":"<div><div>An efficient multiscale coupling method is proposed for simulations of reactor-scale chemical vapor deposition (CVD) with microstructural features. Reactor-scale and microstructure-resolved feature-scale models are coupled through an effective reaction rate formalism, enabling high-resolution deposition simulations while significantly reducing computational cost. A parameterized microstructural model is introduced, in which the relationship between the effective reaction rate and local species consumption rates in the reactor-scale model is directly mapped using precomputed Monte Carlo simulation data. This eliminates the need for iterative calculations or direct numerical simulations of the surface reaction across all the discretized grid points on the wafer, ensuring predictive accuracy while enhancing computational efficiency. Furthermore, an adaptive time-stepping method is developed, dynamically adjusting the time-step size for the feature-scale model based on variations in the effective reaction rate. Through this approach, simulation time is reduced by more than one-third compared to conventional fixed time-step methods, while preserving the accuracy of the effective reaction rate model. The proposed method enables practical and scalable multiscale CVD simulations applicable to industrial reactor design and process optimization, establishing a computationally efficient strategy for integrating reactor-scale and microstructure-resolved feature-scale models.</div></div>","PeriodicalId":50866,"journal":{"name":"Advances in Engineering Software","volume":"212 ","pages":"Article 104051"},"PeriodicalIF":5.7,"publicationDate":"2026-01-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145419535","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2026-01-01Epub Date: 2025-10-23DOI: 10.1016/j.advengsoft.2025.104042
Nastaran Movahedi, Jongyeop Kim, Kyoungsoo Park
This study proposes a straightforward method for computation of stress intensity factors (SIFs) for three-dimensional cracks featuring curved crack fronts and curved crack surfaces, without concerning the mesh topology around a crack tip. The key idea is to utilize a virtual grid-based stress recovery (VGSR) method on arbitrary nonplanar crack geometries along with the interaction energy integral to extract SIFs correspond to different modes. Since the VGSR technique reduces errors related to numerical differentiation and integral domain, it provides accurate computation of SIFs even with arbitrary unstructured meshes. Validation against benchmark problems shows excellent agreement with analytical solutions, highlighting the efficacy of this conjugated approach for precise SIF evaluation. Computational results demonstrate the convergence to analytical solutions while effectively reducing the pointwise oscillations of SIFs under mesh refinement. Additionally, parametric studies are comprehensively performed concerning the characteristics associated with virtual grid domain size, virtual grid element size, finite element sizes, and number of numerical integration points.
{"title":"Stress intensity factor evaluation for non-planar cracks using virtual grid stress recovery (VGSR) and interaction integral methods","authors":"Nastaran Movahedi, Jongyeop Kim, Kyoungsoo Park","doi":"10.1016/j.advengsoft.2025.104042","DOIUrl":"10.1016/j.advengsoft.2025.104042","url":null,"abstract":"<div><div>This study proposes a straightforward method for computation of stress intensity factors (SIFs) for three-dimensional cracks featuring curved crack fronts and curved crack surfaces, without concerning the mesh topology around a crack tip. The key idea is to utilize a virtual grid-based stress recovery (VGSR) method on arbitrary nonplanar crack geometries along with the interaction energy integral to extract SIFs correspond to different modes. Since the VGSR technique reduces errors related to numerical differentiation and integral domain, it provides accurate computation of SIFs even with arbitrary unstructured meshes. Validation against benchmark problems shows excellent agreement with analytical solutions, highlighting the efficacy of this conjugated approach for precise SIF evaluation. Computational results demonstrate the convergence to analytical solutions while effectively reducing the pointwise oscillations of SIFs under mesh refinement. Additionally, parametric studies are comprehensively performed concerning the characteristics associated with virtual grid domain size, virtual grid element size, finite element sizes, and number of numerical integration points.</div></div>","PeriodicalId":50866,"journal":{"name":"Advances in Engineering Software","volume":"212 ","pages":"Article 104042"},"PeriodicalIF":5.7,"publicationDate":"2026-01-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145340732","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2026-01-01Epub Date: 2025-11-20DOI: 10.1016/j.advengsoft.2025.104069
Dongyang Li , Zhen Chen , Chao Sun
As the fundamental component of ship structures, the hull plate made of steel may be subjected to multi-directional loading and potential cracking damage. The joint action of these factors is extremely harmful to the ultimate strength. In the current paper, a novel strategy is proposed to decouple the implicit interaction relationship of ultimate stresses of hull plates with a through-thickness crack under biaxial compression. The evolution mechanisms of load-carrying capacity influenced by crack faces contact are clarified using a shell-solid mixed finite element model. Then, a simplified approach to establish the FE model of cracked plates is employed in the benchmark parametric analysis, which is validated via experimental results. Subsequently, extensive collapse analysis is conducted to investigate the complicated characteristics of residual ultimate strength considering material and geometric nonlinearities. The coupling effect of plate aspect ratio and slenderness ratio, crack length, angle and location together with in-plane compressive loads is dealt with synthetically. Based on 3360 (1680 × 2) sample points derived from numerical calculation, a set of empirical formulae are reported to predict the reduction factor of axial ultimate strength. Combining these formulae with a projection approach for angular cracks, a generalized closed-form approach is proposed to accurately model the interaction relationships of residual ultimate strength. The generalization of these formulae is verified through independent databases with more than 480 sample points.
{"title":"Closed-form neural network solutions for biaxial compressive strength prediction of cracked steel plates","authors":"Dongyang Li , Zhen Chen , Chao Sun","doi":"10.1016/j.advengsoft.2025.104069","DOIUrl":"10.1016/j.advengsoft.2025.104069","url":null,"abstract":"<div><div>As the fundamental component of ship structures, the hull plate made of steel may be subjected to multi-directional loading and potential cracking damage. The joint action of these factors is extremely harmful to the ultimate strength. In the current paper, a novel strategy is proposed to decouple the implicit interaction relationship of ultimate stresses of hull plates with a through-thickness crack under biaxial compression. The evolution mechanisms of load-carrying capacity influenced by crack faces contact are clarified using a shell-solid mixed finite element model. Then, a simplified approach to establish the FE model of cracked plates is employed in the benchmark parametric analysis, which is validated via experimental results. Subsequently, extensive collapse analysis is conducted to investigate the complicated characteristics of residual ultimate strength considering material and geometric nonlinearities. The coupling effect of plate aspect ratio and slenderness ratio, crack length, angle and location together with in-plane compressive loads is dealt with synthetically. Based on 3360 (1680 × 2) sample points derived from numerical calculation, a set of empirical formulae are reported to predict the reduction factor of axial ultimate strength. Combining these formulae with a projection approach for angular cracks, a generalized closed-form approach is proposed to accurately model the interaction relationships of residual ultimate strength. The generalization of these formulae is verified through independent databases with more than 480 sample points.</div></div>","PeriodicalId":50866,"journal":{"name":"Advances in Engineering Software","volume":"212 ","pages":"Article 104069"},"PeriodicalIF":5.7,"publicationDate":"2026-01-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145571699","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2026-01-01Epub Date: 2025-11-25DOI: 10.1016/j.advengsoft.2025.104067
Tianyong Jiang , Jun Tang , Chunjun Hu , Ke Huang , Xiang Tian , Lei Wang
Fiber-reinforced polymer (FRP) can substitute for steel bars to improve the durability problem of reinforced concrete (RC) beams attributed to corrosion. But A high-precision and interpretable prediction method for the flexural strength of FRP-RC beams has not yet been constructed. This study proposed a genetic algorithm optimized artificial neural network (GA-ANN) model to predict the flexural strength of FRP-RC beams. A database of 166 samples was established to train and validate the model. The input parameters include the FRP reinforcement area, FRP ultimate tensile strength, FRP type, elastic modulus of FRP, concrete compressive strength, beam width, and beam depth. The prediction accuracy and practicability of the GA-ANN model were assessed by comparison with other machine learning (ML) models and design guidelines. A parametric sensitivity analysis was performed based on the proposed model. Finally, the SHapley Additive exPlanation (SHAP) was introduced to investigate the intrinsic mechanisms and the parameter contribution of the ML prediction. The results revealed the GA-ANN model achieves superior prediction performance, with a coefficient of determination (R2) on the validation set of 0.992, which is 1.74% to 6.43% higher than that of other models. Moreover, the trends of flexural strength with the input parameters can be well captured, which is highly consistent with the design guidelines. Interpretability analysis shows that the beam depth and the FRP reinforcement area are the dominant factors affecting flexural strength. This study provides reliable support for the accurate prediction of flexural strength and effective reference for engineering applications.
{"title":"Data-driven prediction for flexural strength of FRP bars reinforced concrete beams based on optimized machine learning and SHAP method","authors":"Tianyong Jiang , Jun Tang , Chunjun Hu , Ke Huang , Xiang Tian , Lei Wang","doi":"10.1016/j.advengsoft.2025.104067","DOIUrl":"10.1016/j.advengsoft.2025.104067","url":null,"abstract":"<div><div>Fiber-reinforced polymer (FRP) can substitute for steel bars to improve the durability problem of reinforced concrete (RC) beams attributed to corrosion. But A high-precision and interpretable prediction method for the flexural strength of FRP-RC beams has not yet been constructed. This study proposed a genetic algorithm optimized artificial neural network (GA-ANN) model to predict the flexural strength of FRP-RC beams. A database of 166 samples was established to train and validate the model. The input parameters include the FRP reinforcement area, FRP ultimate tensile strength, FRP type, elastic modulus of FRP, concrete compressive strength, beam width, and beam depth. The prediction accuracy and practicability of the GA-ANN model were assessed by comparison with other machine learning (ML) models and design guidelines. A parametric sensitivity analysis was performed based on the proposed model. Finally, the SHapley Additive exPlanation (SHAP) was introduced to investigate the intrinsic mechanisms and the parameter contribution of the ML prediction. The results revealed the GA-ANN model achieves superior prediction performance, with a coefficient of determination (<em>R</em><sup>2</sup>) on the validation set of 0.992, which is 1.74% to 6.43% higher than that of other models. Moreover, the trends of flexural strength with the input parameters can be well captured, which is highly consistent with the design guidelines. Interpretability analysis shows that the beam depth and the FRP reinforcement area are the dominant factors affecting flexural strength. This study provides reliable support for the accurate prediction of flexural strength and effective reference for engineering applications.</div></div>","PeriodicalId":50866,"journal":{"name":"Advances in Engineering Software","volume":"212 ","pages":"Article 104067"},"PeriodicalIF":5.7,"publicationDate":"2026-01-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145618141","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2026-01-01Epub Date: 2025-11-25DOI: 10.1016/j.advengsoft.2025.104073
Pengfei Sun , Jiantao Bai , Ran Zhang , Fei Cheng , Xiaojiang Zhang , Wenjie Zuo
Fiber-reinforced composite structures (FRCs) are widely used in engineering. The design of the structural topology, together with fiber paths and cross-sectional size is important for enhancing their structural performance. Therefore, a topology optimization method is proposed that considers fiber orientation and cross-sectional size for FRCs. A bar-embedded model is employed to model the FRCs. The Solid Isotropic Material with Penalization method is applied to optimize the structural topology, whereas the Normal Distribution Fiber Optimization method is used to optimize the fiber orientation and cross-sectional size. The objective is to minimize compliance subject to prescribed matrix and fiber volume fractions. Numerical examples are provided to validate the effectiveness of the proposed method.
{"title":"Topology optimization of fiber-reinforced composite structures considering fiber orientation and cross-sectional size","authors":"Pengfei Sun , Jiantao Bai , Ran Zhang , Fei Cheng , Xiaojiang Zhang , Wenjie Zuo","doi":"10.1016/j.advengsoft.2025.104073","DOIUrl":"10.1016/j.advengsoft.2025.104073","url":null,"abstract":"<div><div>Fiber-reinforced composite structures (FRCs) are widely used in engineering. The design of the structural topology, together with fiber paths and cross-sectional size is important for enhancing their structural performance. Therefore, a topology optimization method is proposed that considers fiber orientation and cross-sectional size for FRCs. A bar-embedded model is employed to model the FRCs. The Solid Isotropic Material with Penalization method is applied to optimize the structural topology, whereas the Normal Distribution Fiber Optimization method is used to optimize the fiber orientation and cross-sectional size. The objective is to minimize compliance subject to prescribed matrix and fiber volume fractions. Numerical examples are provided to validate the effectiveness of the proposed method.</div></div>","PeriodicalId":50866,"journal":{"name":"Advances in Engineering Software","volume":"212 ","pages":"Article 104073"},"PeriodicalIF":5.7,"publicationDate":"2026-01-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145618139","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
This study presents vibration modeling and analysis of functionally graded porous (FGP) corrugated plates using isogeometric analysis (IGA) and the first-order shear deformation theory (FSDT). FGP sinusoidally corrugated plates (SCPs) and arc corrugated plates (ACPs) with porosity distributions in the thickness and width directions are studied for the first time. The corrugated plates are modeled by using the non-uniform rational B-splines (NURBS) and multi-patch technique. After building the discretized model for every patch, global mass and stiffness matrices are derived by employing a coordinate transformation between the global and local coordinate systems. The convergence and accuracy of the presented method are validated through comparison with other available data. Then, free vibration behaviors of FGP SCPs and ACPs are analyzed, with a focus on the effects of boundary conditions, porosity distributions, and geometric parameters, such as half-period, half-amplitude, and plate thickness. The results provide benchmark data for future research and offer valuable insights into the advanced structural design and optimization for FGP corrugated plates.
{"title":"Free vibration analysis of functionally graded porous corrugated plates with porosity distributions in the thickness and width directions","authors":"Yaqiang Xue , Chunyu Zhang , Kangkang Shi , Yuan Gao , Zhenyang Gao","doi":"10.1016/j.advengsoft.2025.104066","DOIUrl":"10.1016/j.advengsoft.2025.104066","url":null,"abstract":"<div><div>This study presents vibration modeling and analysis of functionally graded porous (FGP) corrugated plates using isogeometric analysis (IGA) and the first-order shear deformation theory (FSDT). FGP sinusoidally corrugated plates (SCPs) and arc corrugated plates (ACPs) with porosity distributions in the thickness and width directions are studied for the first time. The corrugated plates are modeled by using the non-uniform rational B-splines (NURBS) and multi-patch technique. After building the discretized model for every patch, global mass and stiffness matrices are derived by employing a coordinate transformation between the global and local coordinate systems. The convergence and accuracy of the presented method are validated through comparison with other available data. Then, free vibration behaviors of FGP SCPs and ACPs are analyzed, with a focus on the effects of boundary conditions, porosity distributions, and geometric parameters, such as half-period, half-amplitude, and plate thickness. The results provide benchmark data for future research and offer valuable insights into the advanced structural design and optimization for FGP corrugated plates.</div></div>","PeriodicalId":50866,"journal":{"name":"Advances in Engineering Software","volume":"212 ","pages":"Article 104066"},"PeriodicalIF":5.7,"publicationDate":"2026-01-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145520410","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2026-01-01Epub Date: 2025-11-12DOI: 10.1016/j.advengsoft.2025.104065
Zhen Liu, Liang Xia
Lattice structures have attracted extensive research interest due to their hierarchical architecture and multi-functional properties, enabling unprecedented design flexibility across diverse engineering fields. In general, lattice structure modeling employs two primary methods: boundary representation (B-rep) and implicit representation. The latter is distinguished by its ability to generate lattice structures with more intricate geometries and more diverse functions compared to the former. However, the generated surface mesh the implicitly represented lattice structures is accomplished by feature distortion, non-manifold meshes, and self-intersecting meshes. This not only results in the failure of the generation of body-fitted meshes for finite element analysis (FEA) but also render the performance of additive manufacturing (AM) using the STL model built from the surface mesh impossible. To address these challenges, this work proposes a novel framework of feature-preserving meshing strategies by extending the dual contouring algorithm. The enhanced algorithm outperforms the dual contouring algorithm by ensuring generated surface meshes strictly adhere to topological validity requirements (manifold, closed, oriented), completely eliminating self-intersections, and faithfully preserving sharp geometric features. Subsequently, the remeshing of the surface mesh is performed to optimize the shape and reduce the count of triangles with preserved sharp geometric features, followed by the generation of body-fitted tetrahedral meshes, as depicted in Fig. 1. Finally, the proposed closed-loop mesh generation workflow generates a finite element (FE) model in the standard .inp file format, ensuring compatibility with commercial computational mechanics software (e.g., ABAQUS, ANSYS). Numerical examples show that the proposed meshing workflow is feasible and effective.
{"title":"Feature-preserving mesh generation and simulation for implicitly represented lattice structures","authors":"Zhen Liu, Liang Xia","doi":"10.1016/j.advengsoft.2025.104065","DOIUrl":"10.1016/j.advengsoft.2025.104065","url":null,"abstract":"<div><div>Lattice structures have attracted extensive research interest due to their hierarchical architecture and multi-functional properties, enabling unprecedented design flexibility across diverse engineering fields. In general, lattice structure modeling employs two primary methods: boundary representation (B-rep) and implicit representation. The latter is distinguished by its ability to generate lattice structures with more intricate geometries and more diverse functions compared to the former. However, the generated surface mesh the implicitly represented lattice structures is accomplished by feature distortion, non-manifold meshes, and self-intersecting meshes. This not only results in the failure of the generation of body-fitted meshes for finite element analysis (FEA) but also render the performance of additive manufacturing (AM) using the STL model built from the surface mesh impossible. To address these challenges, this work proposes a novel framework of feature-preserving meshing strategies by extending the dual contouring algorithm. The enhanced algorithm outperforms the dual contouring algorithm by ensuring generated surface meshes strictly adhere to topological validity requirements (manifold, closed, oriented), completely eliminating self-intersections, and faithfully preserving sharp geometric features. Subsequently, the remeshing of the surface mesh is performed to optimize the shape and reduce the count of triangles with preserved sharp geometric features, followed by the generation of body-fitted tetrahedral meshes, as depicted in Fig. 1. Finally, the proposed closed-loop mesh generation workflow generates a finite element (FE) model in the standard .inp file format, ensuring compatibility with commercial computational mechanics software (e.g., ABAQUS, ANSYS). Numerical examples show that the proposed meshing workflow is feasible and effective.</div></div>","PeriodicalId":50866,"journal":{"name":"Advances in Engineering Software","volume":"212 ","pages":"Article 104065"},"PeriodicalIF":5.7,"publicationDate":"2026-01-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145520355","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2025-12-01Epub Date: 2025-09-17DOI: 10.1016/j.advengsoft.2025.104031
Lingkuan Xuan , Gonghao Zhao , Jingfeng Gong , Shengli Su , Yin Yan
This study presents an efficient frequency domain cell vertex finite volume method (FD-CVFVM) to predict muffler transmission loss (TL). The heterogeneous Helmholtz equation is discretized based on FD-CVFVM. A control volume is constructed around each vertex. Acoustic pressure is stored at each mesh vertex. Shape functions are used to describe the acoustic pressure distribution. A calculation program implementing the FD-CVFVM is developed using the C++ language. The TL of a simple expansion chamber muffler, a resistive muffler, and a perforated resistive muffler are computed using different mesh models. The results are compared and analyzed against those obtained from commercial FEM software. The numerical results demonstrate that the FD-CVFVM predictions are in good agreement with those of the FEM results. It is found that the computational efficiency of the FD-CVFVM is significantly superior to that of commercial FEM software. The maximum computation time is reduced by approximately 78.2 %. An analysis of the sparsity pattern of the coefficient matrix is accomplished to reveal the reason of the superior computational speed over the commercial FEM software. This method is anticipated to offer a novel numerical approach for predicting muffler TL.
{"title":"Frequency domain cell-vertex finite volume method for muffler transmission loss prediction","authors":"Lingkuan Xuan , Gonghao Zhao , Jingfeng Gong , Shengli Su , Yin Yan","doi":"10.1016/j.advengsoft.2025.104031","DOIUrl":"10.1016/j.advengsoft.2025.104031","url":null,"abstract":"<div><div>This study presents an efficient frequency domain cell vertex finite volume method (FD-CVFVM) to predict muffler transmission loss (TL). The heterogeneous Helmholtz equation is discretized based on FD-CVFVM. A control volume is constructed around each vertex. Acoustic pressure is stored at each mesh vertex. Shape functions are used to describe the acoustic pressure distribution. A calculation program implementing the FD-CVFVM is developed using the <em>C</em>++ language. The TL of a simple expansion chamber muffler, a resistive muffler, and a perforated resistive muffler are computed using different mesh models. The results are compared and analyzed against those obtained from commercial FEM software. The numerical results demonstrate that the FD-CVFVM predictions are in good agreement with those of the FEM results. It is found that the computational efficiency of the FD-CVFVM is significantly superior to that of commercial FEM software. The maximum computation time is reduced by approximately 78.2 %. An analysis of the sparsity pattern of the coefficient matrix is accomplished to reveal the reason of the superior computational speed over the commercial FEM software. This method is anticipated to offer a novel numerical approach for predicting muffler TL.</div></div>","PeriodicalId":50866,"journal":{"name":"Advances in Engineering Software","volume":"211 ","pages":"Article 104031"},"PeriodicalIF":5.7,"publicationDate":"2025-12-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145096626","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2025-12-01Epub Date: 2025-10-14DOI: 10.1016/j.advengsoft.2025.104043
Gerhard Venter, Melody Neaves
Digital Image Correlation (DIC) is a powerful non-contact optical technique used to measure displacement and strain fields. While commercial tools exist, their high cost and lack of flexibility can hinder innovation and accessibility. This paper introduces SUN-DIC, a Python-based open-source DIC tool designed for extensibility, clarity, and ease of use. We outline the theoretical foundations and describe the software architecture. SUN-DIC’s utility is demonstrated through selected examples and the code is successfully benchmarked against the DIC Challenge 2.0 for 2D analyses.
{"title":"SUN-DIC: A Python-based open-source software tool for Digital Image Correlation","authors":"Gerhard Venter, Melody Neaves","doi":"10.1016/j.advengsoft.2025.104043","DOIUrl":"10.1016/j.advengsoft.2025.104043","url":null,"abstract":"<div><div>Digital Image Correlation (DIC) is a powerful non-contact optical technique used to measure displacement and strain fields. While commercial tools exist, their high cost and lack of flexibility can hinder innovation and accessibility. This paper introduces SUN-DIC, a Python-based open-source DIC tool designed for extensibility, clarity, and ease of use. We outline the theoretical foundations and describe the software architecture. SUN-DIC’s utility is demonstrated through selected examples and the code is successfully benchmarked against the DIC Challenge 2.0 for 2D analyses.</div></div>","PeriodicalId":50866,"journal":{"name":"Advances in Engineering Software","volume":"211 ","pages":"Article 104043"},"PeriodicalIF":5.7,"publicationDate":"2025-12-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145321138","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}