Computers & Structures最新文献

{"title":"Dynamic characterization of cross-physics coupling strengths, a methodology to coupling and reordering partitioned solvers for multiphysics applications","authors":"Christopher Nahed ,&nbsp;Jacques de Lamare","doi":"10.1016/j.compstruc.2024.107615","DOIUrl":"10.1016/j.compstruc.2024.107615","url":null,"abstract":"<div><div>The role of dimensionless ratios in engineering and physics is ubiquitous; but their utility in the multiphysics community is sometimes overlooked. Notably, in the multiphysics modelling community, coupling methods are often discussed and developed without an explicit monitoring of the various dimensionless ratios of the various inter-physics coupling terms. However, it is evident that the varying strengths of the coupling terms in a multiphysics model of <em>k</em> physics solvers/modules will influence either the convergence rate, the stability of the coupling scheme and the program execution speed. In fact, it is well known that the “ordering” of the predictor physics modules is primordial to the performance characteristics of a multiphysics coupling scheme. However, the question of “how to order” (who came first, the chicken or the egg?) the <em>k</em> physics modules remains vaguely discussed. In fact, physics ordering is generally based on the scientist's experience or on problem specific stability analyses performed on academic computational configurations. In the case of generic multiphysics coupling, where volume, interface and/or surface coupling terms can manifest, the optimal ordering of the physics modules may strongly vary along simulation time (for the same application) and/or across applications. Motivated to find an approximate measure that does not resort to cumbersome and problem specific stability analyses, we borrow the concept of dimensionless numbers from physics and apply it to the algebraic systems that manifest in multiphysics computational models. The “chicken-egg” algorithm is based on a dimensionless methodology that serves to “reorder” the Jacobian matrix of an exact Newton-Raphson implicit scheme. The method poses a dimensionless preconditioner that estimates the different strengths of the various coupling terms found in the multiphysics application. The chicken-egg algorithm estimates at every given time step the order of magnitude of coupling terms and correspondingly orders the <em>k</em> partitioned physics solvers automatically. This algorithm is tested for the first time on a thermo-hygro-corrosive multiphysics model and shows promising results. Benchmarking against monolithic and diagonalised calculation strategies, the first numerical tests show a significant reduction in iterations before convergence and thus over a 1.7-fold improvement in program execution time.</div></div>","PeriodicalId":50626,"journal":{"name":"Computers & Structures","volume":"307 ","pages":"Article 107615"},"PeriodicalIF":4.4,"publicationDate":"2025-01-15","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142841244","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Modular-topology optimization for additive manufacturing of reusable mechanisms","authors":"Marek Tyburec ,&nbsp;Martin Doškář ,&nbsp;Michael Somr ,&nbsp;Martin Kružík ,&nbsp;Jan Zeman","doi":"10.1016/j.compstruc.2024.107630","DOIUrl":"10.1016/j.compstruc.2024.107630","url":null,"abstract":"<div><div>Modular designs have gained popularity because they can generally address manufacturing efficiency, reusability, and sustainability concerns. Here, we contribute to the growing field by proposing a fully automatic design method for modules utilized in several products. Our manufacturing-aware procedure is composed of three consecutive steps: (i) free-material optimization for obtaining the optimal spatially distributed elasticity tensors, (ii) hierarchical clustering of the stiffness tensors directly into a given number of modules while allowing for their flipping, and (iii) single-scale topology optimization with manufacturing aspects to design the final topology of mechanically compatible modules. These aspects include connectivity constraints to ensure the integrity of all modules and the three-field projection to account for manufacturing inaccuracies. We illustrate the entire procedure with the design and fabrication of eight different reusable modules assembled into well-functioning modular inverter and gripper mechanisms. These mechanisms were 3D printed and subjected to mechanical testing using an in-house testing machine and digital image correlation. The experimental results show excellent agreement between the predicted and observed behavior and highlight the potential of the method for scalable additive manufacturing.</div></div>","PeriodicalId":50626,"journal":{"name":"Computers & Structures","volume":"307 ","pages":"Article 107630"},"PeriodicalIF":4.4,"publicationDate":"2025-01-15","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142884294","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}

{"title":"Bayesian updating using accelerated Hamiltonian Monte Carlo with gradient-enhanced Kriging model","authors":"Qiang Li,&nbsp;Pinghe Ni,&nbsp;Xiuli Du,&nbsp;Qiang Han,&nbsp;Kun Xu,&nbsp;Yulei Bai","doi":"10.1016/j.compstruc.2024.107598","DOIUrl":"10.1016/j.compstruc.2024.107598","url":null,"abstract":"<div><div>Bayesian methods have been widely used to improve the accuracy of finite element model in civil engineering. However, Bayesian methods generally suffer from the computational complexity involved in accurately identifying the posterior distribution. To address this issue, this paper proposes a novel method by combining the Hamiltonian Monte Carlo (HMC) algorithm with the gradient-enhanced Kriging (GEK) model, termed HMC-GEK, for more efficient model updating. The proposed method uses the potential function and gradient information generated during the burn-in phase of the HMC to train the GEK model. By replacing high-cost potential function with the GEK model, the original HMC sampling process is accelerated. An eight-story frame structure and a Y-shaped arch bridge are used to validate the accuracy and efficiency of the proposed method. Furthermore, the HMC-GEK method has been employed to identify damage of a real eight-story steel frame structure. Compared with the HMC method with the traditional Kriging model, HMC-GEK makes more full use of the gradient information of the potential function and significantly improves the sample acceptance rate and computational efficiency. In addition, the successful application of the method in damage identification of the real structure demonstrates its value for engineering applications.</div></div>","PeriodicalId":50626,"journal":{"name":"Computers & Structures","volume":"307 ","pages":"Article 107598"},"PeriodicalIF":4.4,"publicationDate":"2025-01-15","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142804473","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}

{"title":"Non-parametric ground motion model for displacement response spectra and Fling for Himalayan region using machine learning","authors":"Jyothi Yedulla,&nbsp;Ravi Kanth Sriwastav,&nbsp;S.T.G. Raghukanth","doi":"10.1016/j.compstruc.2024.107626","DOIUrl":"10.1016/j.compstruc.2024.107626","url":null,"abstract":"<div><div>Displacement response spectra (DRS) are crucial for seismic design as earthquake damage correlates more with displacements than forces. Previous efforts to develop attenuation relations for DRS have been largely approximate. Permanent displacement or Fling poses significant design, repair and rehabilitation challenges. Consideration of DRS and Fling in seismic design and performance assessment necessitates its accurate estimation. This paper presents the two Artificial Neural Network (ANN)-based non-parametric Ground Motion Models (GMMs). The first model predicts DRS for horizontal and vertical spectral ordinates. The second model focuses on predicting the Fling step in fault-parallel, fault-normal and vertical components. Both the models are developed for the Himalayan region. Given the limited availability of recorded data, ground motion recorded in tectonically similar regions is also utilized to develop DRS GMM. The sparsely recorded Fling data in the Himalayan region is supplemented by additional Fling values simulated using a physics-based approach, alongside data recorded from tectonically similar regions. The simulated Fling values are validated against recorded Fling data. The performance of developed GMMs is compared with existing GMMs and seismic codes which demonstrated its satisfactory performance. The correlation coefficient for ordinates of DRS and Fling are reported to be greater than 0.86 and 0.80, respectively.</div></div>","PeriodicalId":50626,"journal":{"name":"Computers & Structures","volume":"307 ","pages":"Article 107626"},"PeriodicalIF":4.4,"publicationDate":"2025-01-15","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142841190","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}

{"title":"Improved hexahedral mesh generation from quadrilateral surface meshes","authors":"Jingchen Gao,&nbsp;Zhoufang Xiao,&nbsp;Shuwei Shen,&nbsp;Chenhao Xu,&nbsp;Jingjing Cai,&nbsp;Gang Xu","doi":"10.1016/j.compstruc.2024.107620","DOIUrl":"10.1016/j.compstruc.2024.107620","url":null,"abstract":"<div><div>The quadrilateral surface mesh modification method based on dual cycle operations shows promising advantages in hexahedral mesh generation. However, as only simple cycle eliminations are considered, the existing methods can not handle complex surface meshes. In this study, an improved method based on cycle elimination is proposed for high-quality hexahedral mesh generation from a given quadrilateral mesh. For complex cycles, the dual cycles of a quadrilateral mesh are classified into different types, i.e., concave cycles, flat cycles and convex cycles, and processed accordingly. Then, a novel reversible chord removal operation is proposed to transform the concave cycles into convex cycles. Besides, the interlocked convex cycles are also split to new cycles suitable for elimination with the proposed chord removal operation. After that, the convex cycles can be eliminated. To find the best cycle for elimination, the convex cycles are ranked with geometric and topology information and the cycle with the highest rank is chosen to be eliminated. The above steps repeat until only a hexahedron exists and the final hexahedral mesh can be obtained by reversing the order of the surface modification operations. The effectiveness of the proposed method is validated with dozens of surface meshes.</div></div>","PeriodicalId":50626,"journal":{"name":"Computers & Structures","volume":"307 ","pages":"Article 107620"},"PeriodicalIF":4.4,"publicationDate":"2025-01-15","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142841193","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}

{"title":"An enrichment of Q4γs plate finite element using incomplete quadratic functions, an assumed energy orthogonality of Bergan’s free formulation, and mixed transverse shear strains","authors":"Andi Makarim Katili ,&nbsp;Kai-Uwe Bletzinger ,&nbsp;Irwan Katili","doi":"10.1016/j.compstruc.2024.107619","DOIUrl":"10.1016/j.compstruc.2024.107619","url":null,"abstract":"<div><div>This paper introduces a new quadrilateral plate element named DSPM4, which improves upon the previous Q4γ<em>_s</em> element. The DSPM4 element has twelve DOFs and four temporary DOFs at the mid-sides of the element. The rotation functions βs are modified by adding an incomplete quadratic function to improve the bending performance. An orthogonality condition between the lower and higher-order bending energy is assumed to satisfy the constant bending patch test. The element development involves kinematic and mixed transverse shear strains. In the local normal-tangential (<em>n-s</em>) coordinate system of each side, the mixed transverse shear strains are obtained by combining the constitutive law for shear forces, equilibrium equations, the constitutive law for bending, and curvatures. The relationship between the temporary DOFs and the final DOFs at the corner nodes is established by equating the constant kinematic and mixed transverse shear strains in the direction cosines of each side of the element. As in the Q4γ<em>_s</em> element, the discrete shear projection method is applied to avoid the shear-locking effect. Numerical analyses are conducted for thin and thick isotropic plates, and the results are compared with those of the DKMQ and MITC4 elements. Numerical results prove that the DSPM4 element is free of shear locking in thin and thick plates, has a proper rank, and provides better convergence behaviour than the Q4γ_s element. Its performance in static and free vibration analyses of isotropic plates is highly comparable to the DKMQ element.</div></div>","PeriodicalId":50626,"journal":{"name":"Computers & Structures","volume":"307 ","pages":"Article 107619"},"PeriodicalIF":4.4,"publicationDate":"2025-01-15","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142841246","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}

{"title":"On the use of the inverse finite element method to enhance knowledge sharing in population-based structural health monitoring","authors":"Giulia Delo ,&nbsp;Rinto Roy ,&nbsp;Keith Worden ,&nbsp;Cecilia Surace","doi":"10.1016/j.compstruc.2024.107635","DOIUrl":"10.1016/j.compstruc.2024.107635","url":null,"abstract":"<div><div>Efficient Structural Health Monitoring (SHM) is critical for ensuring safety and improving the operation and maintenance of aerospace structures. This study focusses on advanced shape-sensing methods, such as the inverse Finite Element Method (iFEM), which can estimate the complete displacement field of a structure based on a restricted number of strain measurements, fostering continuous and real-time monitoring. This approach additionally provides valuable insights into the dynamic behaviour of a structure by extracting its Frequency Response Functions (FRFs) and modal properties to perform vibration-based SHM. However, effectively extending SHM to a fleet or population of structures would require a significant amount of data for each one, which may be unavailable or incomplete. A population-based Structural Health Monitoring (PBSHM) strategy can solve data scarcity by sharing knowledge between similar structures via transfer-learning algorithms. In PBSHM, handling data from diverse sources is paramount for achieving accurate results. Therefore, this study integrates iFEM into the PBSHM framework, enhancing knowledge transfer by harmonising fibre-optic strain measurements to vibration-based features and providing reliable source data to inform diagnostics on similar structures. The proposed approach is validated on a population of laboratory-scale steel aircraft subjected to specific operating and damage conditions tested using three different sensor setups.</div></div>","PeriodicalId":50626,"journal":{"name":"Computers & Structures","volume":"307 ","pages":"Article 107635"},"PeriodicalIF":4.4,"publicationDate":"2025-01-15","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143095111","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Open-source implementations and comparison of explicit and implicit crystal-plasticity finite element methods","authors":"Hassan M. Asadkandi ,&nbsp;Tomáš Mánik ,&nbsp;Bjørn Holmedal ,&nbsp;Odd Sture Hopperstad","doi":"10.1016/j.compstruc.2024.107621","DOIUrl":"10.1016/j.compstruc.2024.107621","url":null,"abstract":"<div><div>In this study, two state-of-the-art implementations of the rate-dependent Crystal Plasticity Finite Element Method (CPFEM) as user material subroutines in the finite element solvers Abaqus/Explicit and Abaqus/Standard (Implicit) are presented. Adaptive substepping in the explicit solver and line-search stabilized implementation in the implicit solver enable fast and stable calculations also for small strain-rate sensitivities. A comparative analysis of simulation results is made and computation times of the two implementations are compared and discussed, including the use of reduced-integration elements and mass scaling for the explicit solver. It is identified under which conditions the explicit or the implicit solver is preferred.</div></div>","PeriodicalId":50626,"journal":{"name":"Computers & Structures","volume":"307 ","pages":"Article 107621"},"PeriodicalIF":4.4,"publicationDate":"2025-01-15","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142841392","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Simultaneous optimization of topology and bi-material distribution of three-dimensional structures for addressing local heat accumulation in layer-upon-layer additive manufacturing process","authors":"Nima Yaghoobi,&nbsp;Mohammad Hossein Abolbashari","doi":"10.1016/j.compstruc.2024.107632","DOIUrl":"10.1016/j.compstruc.2024.107632","url":null,"abstract":"<div><div>This paper introduces a novel approach based on a topology optimization (TO) model to efficiently distribute material phases for minimizing structural compliance and enhance local heat evacuation in additive manufacturing (AM). The approach simultaneously optimizes the structure for its intended function and behavior during layer-by-layer production. While AM allows intricate, topologically optimal multi-material designs, it often induces high temperatures and heat fluxes, risking part failure and compromising mechanical properties. To address this, a 3D gradient-based bi-material and functionally graded TO method is presented, considering total volume percentage and a thermal constraint based on temperatures of local sub-domains. The methodology involves density-based multi-material TO, interpolating elastic moduli, thermal conductivity, and heat flux based on proposed extensions of solid isotropic material with penalization method. Subsequently, a steady state analysis is performed in each sub-domain whose top element layer is exposed to a heat flux, simulating AM process. Both 2D and 3D numerical results demonstrate the contribution of the presented approach in preventing localized heating-induced geometrical patterns in AM. Additionally, the proposed method proves effective in producing superior designs without the need for sacrificial support structures in bi-material and functionally graded material designs, offering self-supported structures.</div></div>","PeriodicalId":50626,"journal":{"name":"Computers & Structures","volume":"307 ","pages":"Article 107632"},"PeriodicalIF":4.4,"publicationDate":"2025-01-15","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142884291","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}

{"title":"Corrigendum to “Topology optimization of lattice structures for target band gaps with optimum volume fraction via Bloch-Floquet theory” [Comput. Struct. 307 (2025) 107601]","authors":"F. Gómez-Silva ,&nbsp;R. Zaera ,&nbsp;R. Ortigosa ,&nbsp;J. Martínez-Frutos","doi":"10.1016/j.compstruc.2025.107649","DOIUrl":"10.1016/j.compstruc.2025.107649","url":null,"abstract":"","PeriodicalId":50626,"journal":{"name":"Computers & Structures","volume":"307 ","pages":"Article 107649"},"PeriodicalIF":4.4,"publicationDate":"2025-01-15","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143095110","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}