Finite Elements in Analysis and Design最新文献

{"title":"Stress-based topology optimization of continuum structures incorporating a piecewise P-norm stabilization strategy","authors":"Shi-An Zhou,&nbsp;Song Yao","doi":"10.1016/j.finel.2025.104401","DOIUrl":"10.1016/j.finel.2025.104401","url":null,"abstract":"<div><div>In this paper, a combined Floating Projection Topology Optimization (FPTO) method is applied to stress-based topology optimization of continuum structures. To balance the need for a stable optimization process under a low P-norm and the goal of achieving a uniform stress distribution at a higher P-norm, the FPTO method here is combined with piecewise P-norm strategy. A new parameter, solid rate, is introduced, serving not only as an additional convergence criterion but also as a decision criterion for adaptively increasing the P-norm value. To validate the effectiveness of the proposed FPTO method incorporating both piecewise P-norm stabilization and the solid rate control strategy, this study conducts stiffness maximization optimization with stress constraints on various typical structures. The results show that the method effectively reduces computational oscillations in stress-based topology optimization for structures with stress concentrations, which helps to the discovery of superior designs.</div></div>","PeriodicalId":56133,"journal":{"name":"Finite Elements in Analysis and Design","volume":"250 ","pages":"Article 104401"},"PeriodicalIF":3.5,"publicationDate":"2025-06-26","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144480638","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"A finite element implementation of the SRTD algorithm for an Oldroyd 3-parameter viscoelastic fluid model","authors":"Christian Austin ,&nbsp;Sara Pollock ,&nbsp;L. Ridgway Scott","doi":"10.1016/j.finel.2025.104399","DOIUrl":"10.1016/j.finel.2025.104399","url":null,"abstract":"<div><div>In this paper, we discuss a finite element implementation of the SRTD algorithm described by Girault and Scott for the steady-state case of a certain 3-parameter subset of the Oldroyd models. We compare it to the well-known EVSS method, which, though originally described for the upper-convected Maxwell model, can easily accommodate the Oldroyd 3-parameter model. We obtain numerical results for both methods on two benchmark problems: the lid-driven cavity problem and the journal-bearing, or eccentric rotating cylinders, problem. We find that the resulting finite element implementation of SRTD is stable with respect to mesh refinement and is generally faster than EVSS, though is not capable of reaching as high a Weissenberg number as EVSS.</div></div>","PeriodicalId":56133,"journal":{"name":"Finite Elements in Analysis and Design","volume":"250 ","pages":"Article 104399"},"PeriodicalIF":3.5,"publicationDate":"2025-06-25","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144471163","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Investigation on an optimal aggregation level for a parallel meshless multigrid method based on domain decomposition method","authors":"Sang Truong Ha ,&nbsp;Hyeong Cheol Park ,&nbsp;Han Young Yoon ,&nbsp;Hyoung Gwon Choi","doi":"10.1016/j.finel.2025.104402","DOIUrl":"10.1016/j.finel.2025.104402","url":null,"abstract":"<div><div>We developed three parallel algorithms for a meshless geometric multigrid (GMG) method proposed for the finite element discretization of elliptic partial differential equation. These methods for parallel multigrid (PMG) are based on the message passing interface (MPI) for domain decomposition and coarse matrix aggregation (CMA) algorithm for coarser levels. Using coarse matrices obtained by a parallel Galerkin condition for the present meshless GMG, we proposed a parameter by which an optimal aggregation level is determined. This parameter is defined as the ratio of total number of external interface nodes from all the subdomains before aggregation to the number of non-zero entries of gathered matrix after aggregation. Three methods <strong>—M1, M2, and M3—</strong> are classified depending on how the coarsest matrix is solved and the number of coarser levels for which CMA is applied. M1 (M2) solves the coarsest matrix via an iterative (direct) solver applying CMA only for the coarsest level, whereas M3 determines the multigrid levels with CMA based on the parameter and employs a direct solver for the coarsest matrix. We found that M3 is more efficient than the others and much more efficient in the case of complicated geometry because communication overhead is reduced compared to the other methods. Furthermore, the present PMG could achieve super-linear scalability owing to the cache effect for a large problem.</div></div>","PeriodicalId":56133,"journal":{"name":"Finite Elements in Analysis and Design","volume":"250 ","pages":"Article 104402"},"PeriodicalIF":3.5,"publicationDate":"2025-06-18","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144307470","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Multi-material structure topology optimization considering length scale control based on isogeometric analysis approach","authors":"Jianping Zhang,&nbsp;Ou Guo,&nbsp;Yaping Zhao,&nbsp;Zhijian Zuo,&nbsp;Mintao Chen,&nbsp;Haishan Lu,&nbsp;Shuguang Gong","doi":"10.1016/j.finel.2025.104391","DOIUrl":"10.1016/j.finel.2025.104391","url":null,"abstract":"<div><div>The multi-material structure topology optimization method considering length scale control (LSC) using the isogeometric analysis approach is put forward. The density distribution function (DDF) is applied for improving structural smoothness, and the resulting structure has smoother and clearer boundaries compared to the conventional finite element method. The alternating active-phase algorithm and gradient algorithm are utilized for building the multi-material interpolation model and updating the design variable, respectively. The effects of different LSC scheme, the maximum length scale control (MaxLSC) domain radius <span><math><mrow><msub><mi>R</mi><mi>max</mi></msub></mrow></math></span>, the radius of DDF influence domain <span><math><mrow><msub><mi>r</mi><mtext>fil</mtext></msub></mrow></math></span> and the aggregation factor <span><math><mrow><msub><mi>p</mi><mi>n</mi></msub></mrow></math></span> on the structure performance are investigated. The results show the structure with both MaxLSC and MinLSC applied exhibits a more uniform material distribution, the structural manufacturability is effectively guaranteed. When <span><math><mrow><msub><mi>r</mi><mtext>fil</mtext></msub><mo>≤</mo><mn>3</mn><mi>h</mi></mrow></math></span>, the jagged boundaries appear in the structure, and when <span><math><mrow><msub><mi>r</mi><mtext>fil</mtext></msub><mo>≥</mo><mn>5</mn><mi>h</mi></mrow></math></span>, the branching structures decrease. The topological structure obtained when <span><math><mrow><msub><mi>r</mi><mtext>fil</mtext></msub><mo>=</mo><mn>3.5</mn><mi>h</mi><mo>−</mo><mn>4</mn><mi>h</mi></mrow></math></span> has the relatively uniform material distribution. When <span><math><mrow><msub><mi>R</mi><mi>max</mi></msub><mo>=</mo><mn>5</mn><mi>h</mi></mrow></math></span>, the island phenomenon appears in the structure. When <span><math><mrow><msub><mi>R</mi><mi>max</mi></msub><mo>=</mo><mn>10</mn><mi>h</mi></mrow></math></span>, the branching structure is reduced and thickened simultaneously, the recommended range for <span><math><mrow><msub><mi>R</mi><mi>max</mi></msub></mrow></math></span> is <span><math><mrow><mn>6</mn><mi>h</mi><mo>−</mo><mn>8</mn><mi>h</mi></mrow></math></span>. When <span><math><mrow><msub><mi>p</mi><mi>n</mi></msub><mo>=</mo><mn>80</mn><mo>−</mo><mn>110</mn></mrow></math></span>, the topological structure exhibits more branching structure in both materials. It is proved that the effectiveness of the LSC method can still be guaranteed in the three-dimensional problem and curved edge structure.</div></div>","PeriodicalId":56133,"journal":{"name":"Finite Elements in Analysis and Design","volume":"250 ","pages":"Article 104391"},"PeriodicalIF":3.5,"publicationDate":"2025-06-14","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144280979","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Profilometry-based indentation plastometry (PIP) integrated finite element (FE) model validation using neutron diffraction","authors":"Gaurav Vithalani ,&nbsp;Stuart Bell ,&nbsp;Rezwanul Haque ,&nbsp;Edmund Pickering ,&nbsp;Mark Reid ,&nbsp;Geoffrey Will ,&nbsp;Richard Clegg ,&nbsp;Theodore A. Steinberg","doi":"10.1016/j.finel.2025.104400","DOIUrl":"10.1016/j.finel.2025.104400","url":null,"abstract":"<div><div>Stress assisted corrosion can lead to premature failure of structural alloys such as stainless steel 316. The mechanisms governing the various failure modes are poorly understood due to the complexities involved in testing the combined effects of stress and corrosion at high temperatures. Modified Turski (MT) specimens are an invaluable tool for studying these mechanisms as they can be pre-stressed to induce residual stresses and exposed to a corrosive environment without needing specialised equipment. However, it is challenging to measure the stresses experimentally without altering the stress state. To overcome this, finite element (FE) modelling can be used, but validated models are scarce, particularly for MT specimen geometry. In this work, an FE model was developed based on the non-destructive characterisation of material using profilometry-based indentation Plastometry (PIP) and validated for residual stresses using enhanced neutron diffraction (ND) measurement. The model results showed significant correlations with experimentally observed deflection and residual stress profile. The deviation observed between the model and experiment was analysed and discussed to understand the evolution of strain field in materials. To conclude, this validated model, and approach is reliable and enables rapid material compatibility assessment for stress assisted corrosion, which is particularly important in emerging industries like concentrated solar thermal (CST) and hydrogen energy, amongst many others.</div></div>","PeriodicalId":56133,"journal":{"name":"Finite Elements in Analysis and Design","volume":"250 ","pages":"Article 104400"},"PeriodicalIF":3.5,"publicationDate":"2025-06-11","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144263687","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Parallel computing on GPU with CuPy and vectorized SpMV for large-scale topology optimization","authors":"Jiangnan Hou ,&nbsp;Jiajie Li ,&nbsp;Shengfeng Zhu ,&nbsp;Xindi Hu ,&nbsp;Zeyang Yu","doi":"10.1016/j.finel.2025.104388","DOIUrl":"10.1016/j.finel.2025.104388","url":null,"abstract":"<div><div>This paper presents a vectorized programming approach to Sparse Matrix-Vector multiplication (SpMV) based on a novel decomposition of the finite element matrix for GPU. The approach is employed for large-scale topology optimization. Unlike the existing matrix-free SpMV strategies (e.g., EbE, NbN, and DbD), our framework transforms SpMV into vector operations. This improvement enhances accessibility for non-expert users in Python environment (e.g., CuPy) while maintaining computational efficiency. Numerical examples of topology optimization on GPU, such as heat transfer and structural design, in both 2D and 3D with up to 63 million elements illustrate the effectiveness and efficiency of the proposed approach. Codes could be seen at <span><span>https://github.com/HouJiangnan81/CuPy-for-TO</span><svg><path></path></svg></span>.</div></div>","PeriodicalId":56133,"journal":{"name":"Finite Elements in Analysis and Design","volume":"250 ","pages":"Article 104388"},"PeriodicalIF":3.5,"publicationDate":"2025-06-11","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144255005","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"A physically based Proper Orthogonal Decomposition for structural mechanics involving dissipative processes","authors":"L. Cimatti Lucarelli ,&nbsp;N. Blal ,&nbsp;A. Platzer ,&nbsp;A. Gravouil ,&nbsp;D. Iampietro","doi":"10.1016/j.finel.2025.104386","DOIUrl":"10.1016/j.finel.2025.104386","url":null,"abstract":"<div><div>This paper proposes a framework for a physically based POD (Proper Orthogonal Decomposition), denoted <span><math><mi>φ</mi></math></span>POD. The proposed Reduced Order Model (ROM) is set to respect one of the evolution equations of an elasto-viscoplastic material. The <span><math><mi>φ</mi></math></span>POD is obtained through the solution of a constrained minimization problem over the reduced coordinates of the POD. Furthermore, a physically based Reduced Order Surrogate Model (ROSM) is constructed based on <em>Interpolation on a Tangent Space to the Grassmann manifold</em> (ITSGM). This paper uses the 2D Cook’s membrane under plane strain hypothesis as an application case. The performance metrics to compare the <span><math><mi>φ</mi></math></span>POD to the regular POD are the residual of the constrained equation, projection energies, and interpolation error for mechanical fields.</div></div>","PeriodicalId":56133,"journal":{"name":"Finite Elements in Analysis and Design","volume":"250 ","pages":"Article 104386"},"PeriodicalIF":3.5,"publicationDate":"2025-06-09","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144239800","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"A parallel FreeFEM framework for topology optimization of structures into three spatial dimensions","authors":"J.M.M. Luz Filho,&nbsp;A.T.A. Gomes,&nbsp;A.A. Novotny","doi":"10.1016/j.finel.2025.104384","DOIUrl":"10.1016/j.finel.2025.104384","url":null,"abstract":"<div><div>This article presents a simple and compact parallel FreeFEM implementation for topology optimization of structures into three spatial dimensions. The topology optimization algorithm relies on the topological derivative concept together with a level-set domain representation method, in which the topological derivative is used as a steepest descent direction evolving the level-set function within the optimization procedure. In addition, adaptive mesh refinement is performed in the optimization scheme for enhancing the boundary representation of the resulting structure and reducing the computational cost of structured mesh refinements. In addition, a rich and thorough discussion of the main aspects of the parallel FreeFEM implementation is given in full detail. In particular, a version of the resulting code is provided for the reader convenience, including two benchmark examples for structural stiffness maximization and for the design of compliant mechanisms, showing the simplicity and effectiveness of the proposed implementation. Besides, we also show how the resulting code can be easily converted into a SIMP-based topology design algorithm. Since all the background knowledge on the topological derivative method is presented and the complete FreeFEM implementation is provided as supplementary material, this paper may also be used as educational/pedagogical material for graduate students and/or newcomers to the field of topology optimization.</div></div>","PeriodicalId":56133,"journal":{"name":"Finite Elements in Analysis and Design","volume":"250 ","pages":"Article 104384"},"PeriodicalIF":3.5,"publicationDate":"2025-06-05","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144221768","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Temperature fields calculation in heat exchangers using the finite element method","authors":"Jose M. Chaquet ,&nbsp;Pedro Galán del Sastre","doi":"10.1016/j.finel.2025.104385","DOIUrl":"10.1016/j.finel.2025.104385","url":null,"abstract":"<div><div>Heat exchanger (HEX) design is an optimization process that seeks to maximize heat transfer between two fluids while minimizing pressure drops. There are several conceptual design methods based on integral equations that only work with specific temperature values at the inlet and outlet of the HEX. However, it is very interesting to obtain approximate temperature distributions in these early stages of analysis to verify that the design criteria are met. To do this, it is necessary to solve systems of differential equations depending on the HEX configuration. Under certain assumptions, these equations have an analytical solution. However, in most cases it is only possible to obtain a numerical approximation. This work presents the solution of these equations for 6 HEX arrangements based on the finite element method. After validating the results in the cases whose analytical solution is known, the proposed method is applied to two realistic cases. Firstly, the effectiveness of a double-pass crossflow heat exchanger is studied. Since the inlet distributions in the second HEX module are not constant, no analytical solution is available. The numerical solution allows to analyze under what conditions the second pass is not effective due to thermal inversion. Secondly, a simplified 2D geometry of a compact intercooler HEX for hydrogen-fueled aero engine is solved. Specifically, an analysis is carried out to locate possible malfunctions due to obstructions of the fluid flows making use only of the metal temperature distributions at the outlets.</div></div>","PeriodicalId":56133,"journal":{"name":"Finite Elements in Analysis and Design","volume":"249 ","pages":"Article 104385"},"PeriodicalIF":3.5,"publicationDate":"2025-05-30","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144170140","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Optimization of point-melting strategies for the Electron Beam Melting process","authors":"Geovane Augusto Haveroth ,&nbsp;Carl-Johan Thore ,&nbsp;Roberto Federico Ausas ,&nbsp;Stefan Jakobsson ,&nbsp;José Alberto Cuminato ,&nbsp;Maicon Ribeiro Correa","doi":"10.1016/j.finel.2025.104356","DOIUrl":"10.1016/j.finel.2025.104356","url":null,"abstract":"<div><div>This study proposes an optimization methodology to find optimal heat source paths for point-melting in Electron Beam Melting (EBM) Powder Bed Fusion (PBF) processes, aiming to reduce the need for support structures and improve print quality. The building process is simulated using a time-dependent, one-way coupled, non-linear thermo-mechanical model, assuming negligible molten flow, with elastoplastic behavior and temperature-dependent material parameters. The goal of the optimization problem is to find heat source paths that minimize a global temperature measure with a penalty on excessive local temperatures. The numerical methodology is based on solving the non-linear partial differential equations via the Finite Element Method (FEM) and is applied in numerical examples for printing with titanium alloy Ti6Al4V. Metrics related to heat, residual displacement, and residual stresses are considered to assess the performance of different point-melting strategies and to compare optimized and conventional paths. The feasibility of the proposed optimization methodology for practical applications and alternatives towards future methodological advancements are discussed. The study provides a Python-based, MPI-parallelized implementation using open-source libraries and is made available for further research and applications.</div></div>","PeriodicalId":56133,"journal":{"name":"Finite Elements in Analysis and Design","volume":"249 ","pages":"Article 104356"},"PeriodicalIF":3.5,"publicationDate":"2025-05-30","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144170141","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}