Carlos M. Patlán , Hugo Hernández-Barrios , Iván F. Huergo , Francisco Domínguez-Mota
{"title":"Time integration scheme for nonlinear structural dynamics, FAM, including structural vibration control","authors":"Carlos M. Patlán , Hugo Hernández-Barrios , Iván F. Huergo , Francisco Domínguez-Mota","doi":"10.1016/j.compstruc.2025.107645","DOIUrl":"10.1016/j.compstruc.2025.107645","url":null,"abstract":"<div><div>In this study, a method for the integration of the equation of motion for the inelastic analysis of structures utilizing the Force Analogy Method (FAM) and nonlinear control systems is proposed. The method is implicit, unconditionally stable, one-step scheme, multi-stage, with second-order precision, self-start capability, and high-frequency response filtering, exhibiting low overshooting. It enables consideration of sources of nonlinearity from the inelastic behavior of materials and the incorporation of control systems in structures. Four numerical examples are used to validate the proposed method, encompassing a diverse range of application scenarios, including varying numerical stiffness, dynamic load sources, and nonlinearity sources. The obtained results demonstrate excellent agreement with expected solutions, highlighting the method capacity to suppress high-frequency responses while maintaining solution accuracy. Our study suggests that the proposed method holds significant potential as a dynamic integration tool in analyzing complex systems. Its application in structures with nonlinear control systems and material nonlinearity represents a substantial contribution to the field, providing a robust and efficient solution for understanding structural response to dynamic actions.</div></div>","PeriodicalId":50626,"journal":{"name":"Computers & Structures","volume":"308 ","pages":"Article 107645"},"PeriodicalIF":4.4,"publicationDate":"2025-02-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143020309","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}
F. Kazemi , A. Ӧzyüksel Çiftçioğlu , T. Shafighfard , N. Asgarkhani , R. Jankowski
{"title":"RAGN-R: A multi-subject ensemble machine-learning method for estimating mechanical properties of advanced structural materials","authors":"F. Kazemi , A. Ӧzyüksel Çiftçioğlu , T. Shafighfard , N. Asgarkhani , R. Jankowski","doi":"10.1016/j.compstruc.2025.107657","DOIUrl":"10.1016/j.compstruc.2025.107657","url":null,"abstract":"<div><div>The utilization of advanced structural materials, such as preplaced aggregate concrete (PAC), fiber-reinforced concrete (FRC), and FRC beams has revolutionized the field of civil engineering. These materials exhibit enhanced mechanical properties compared to traditional construction materials, offering engineers unprecedented opportunities to optimize the design, construction, and performance of structures and infrastructures. This formal description elucidates the inherent mechanical properties of PAC, FRC, and FRC beams, explores their diverse applications in civil engineering projects. This research aims to propose a surrogate multi-subject ensemble machine-learning (ML) method (named RAGN-R) for estimating mechanical properties of aforementioned advanced materials. The proposed learning approach, RAGN-R, integrates Random forest, Adaptive boosting, and GradieNt boosting techniques, employing a Ridge regression framework for stacking the ensemble. For this purpose, three experimental dataset have been prepared to determine the capability of RAGN-R and the results of the study have been compared with six well-known ML models. It is noteworthy that the proposed RAGN-R has the ability of self-optimizing the hyperparameters, which facilitate the adoptability of the model with engineering problems. Moreover, three datasets have been investigated to show the ability of the RAGN-R for diverse problems. Different performance evaluation metrics have been conducted to present results and compare ML models, which confirms the highest performance of RAGN-R (i.e., 97.7% accuracy) in handling complex relationships and improving overall prediction accuracy.</div></div>","PeriodicalId":50626,"journal":{"name":"Computers & Structures","volume":"308 ","pages":"Article 107657"},"PeriodicalIF":4.4,"publicationDate":"2025-02-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143055292","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}
Álvaro López Varela , Daniel Dopico Dopico , Alberto Luaces Fernández
{"title":"An analytical approach to the sensitivity analysis of semi-recursive ODE formulations for multibody dynamics","authors":"Álvaro López Varela , Daniel Dopico Dopico , Alberto Luaces Fernández","doi":"10.1016/j.compstruc.2024.107642","DOIUrl":"10.1016/j.compstruc.2024.107642","url":null,"abstract":"<div><div>Sensitivity analysis is an extremely powerful tool in many applications such as in the optimization of the dynamics of multibody systems with gradient-based methods. Sensitivity calculations are computationally burdensome and, depending on the method chosen for differentiation and the set of dynamic equations, they could result highly inefficient. Semi-recursive dynamic methods are seldom studied analytically in terms of sensitivity analysis due to their complexity, even though their dynamic performance is usually among the most efficient.</div><div>This work explores the sensitivity analysis of a particular multibody-dynamics formulation, the semi-recursive Matrix R formulation, which is based on the nullspace of constraint equations and leads to a system of ordinary differential equations. As a result, two sets of sensitivity equations are proposed, one based on the direct differentiation method (DDM) and other on the Adjoint Variable Method (AVM), being these sensitivity formulations the main novelty of this work. The main derivatives required in the sensitivity equations are listed in this document, paying special attention to conciseness, correctness and completeness. The methods proposed have been implemented in the general purpose multibody library MBSLIM (<em>Multibody Systems in Laboratorio de Ingeniería Mecánica</em>), and their performance has been tested in two numerical experiments, a five-bar benchmark problem and a four-wheeled buggy vehicle.</div><div>A review and generalization of constrained and unconstrained kinematic problems in relative coordinates is provided as an introduction to the generation of the semi-recursive Matrix R equations of motion. Due to the importance of the selection of the set of independent coordinates, a more general description of the Matrix R method is presented as a novel contribution as well.</div></div>","PeriodicalId":50626,"journal":{"name":"Computers & Structures","volume":"308 ","pages":"Article 107642"},"PeriodicalIF":4.4,"publicationDate":"2025-02-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142986328","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":"An adaptive port technique for synthesising rotational components in component modal synthesis approaches","authors":"Xiang Zhao, My Ha Dao","doi":"10.1016/j.compstruc.2025.107655","DOIUrl":"10.1016/j.compstruc.2025.107655","url":null,"abstract":"<div><div>Component Modal Synthesis (CMS) is a reduced order modelling method widely used for large-scale complex systems. It can effectively approximate system-level models through component synthesis, in which the repetitive geometrical components are modelled once and synthesised together. However, the conventional CMS only applies to systems with stationary components connected by strictly compatible ports, limiting it from modelling systems with moving components. This paper presents an adaptive port (AP) technique to extend CMS approaches for modelling parametric systems with rotational parts. To demonstrate the capability of the AP technique, we apply it to the Static Condensation Reduced Basis Element (SCRBE), one widely used variant of CMS approaches. The AP-based SCRBE (AP-SCRBE) can enforce the synthesis of rotational-stationary components over a shared adaptive port when the connecting surfaces of two components are discretisation-wise incompatible, which happens when one component moves relative to the others. Numerical experiments on the NREL 5 MW wind turbine show that, in the context of rotational-stationary component synthesis, the AP-SCRBE can accurately and efficiently model the rotating rotor with pitch rotation of blades. It can produce almost identical results to a high-fidelity finite element model at two to three orders faster speeds.</div></div>","PeriodicalId":50626,"journal":{"name":"Computers & Structures","volume":"308 ","pages":"Article 107655"},"PeriodicalIF":4.4,"publicationDate":"2025-02-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143055226","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":"A coupled FD-SPH framework for the damage evaluation of ceramic-steel composite structures subjected to blast loading","authors":"Jian-Yu Chen , Xian-Zhao Song , Chong Peng","doi":"10.1016/j.compstruc.2025.107653","DOIUrl":"10.1016/j.compstruc.2025.107653","url":null,"abstract":"<div><div>Ceramic-composite structures are important for developing lightweight vehicles under the threat of high explosives or improvised explosive devices on battlefields. In this paper, the damage process of the ceramic-steel double-layered target subjected to blast loading is simulated by developing a coupled finite difference-smoothed particle hydrodynamics methodology. The shock wave propagation in air medium is simulated using the finite difference method, while the damage and fragmentation of ceramic medium is predicted using smoothed particle hydrodynamics. The information of different physical variables is transferred from finite difference to smoothed particle hydrodynamics by using the immersed boundary method. Firstly, the SPH solver was validated by simulating the high-velocity impact of a ceramic sphere on a ceramic-steel double-layered plate. Afterwards, the damage and fracture mechanisms of ceramic-steel multilayered plates with different types and thicknesses under blast loading were investigated. The numerical results obtained from the coupled finite difference-smoothed particle hydrodynamics approach were compared against the available experimental data, which demonstrates that the developed FD-SPH computational framework is capable of capturing crack propagation and damage patterns of ceramic composite structures well.</div></div>","PeriodicalId":50626,"journal":{"name":"Computers & Structures","volume":"309 ","pages":"Article 107653"},"PeriodicalIF":4.4,"publicationDate":"2025-01-31","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143125006","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}
A.A. Basmaji , A. Fau , R. Desmorat , U. Nackenhorst
{"title":"Cyclic nonlocal anisotropic damage modelling of concrete mesostructures with real-shaped aggregates","authors":"A.A. Basmaji , A. Fau , R. Desmorat , U. Nackenhorst","doi":"10.1016/j.compstruc.2025.107650","DOIUrl":"10.1016/j.compstruc.2025.107650","url":null,"abstract":"<div><div>Modelling and computing concrete mesostructures subjected to loads alternating between tension and compression are challenging. This paper presents a full computational model, from the random packing of real-shaped aggregates at the meso-scale to the FE computations with nonlocal anisotropic damage for alternate (cyclic) loading. <em>Concrete</em> is represented as a two-phase random heterogeneous material consisting of mortar and aggregates. Aggregates diversity and realism are reproduced by importing shapes from a laser-scanning database. Each elastic aggregate is described by several hundreds of facet vertices, packed using Oriented Bounding Boxes. The anisotropic nature of damage induced by loading in concrete is reproduced by assigning a nonlocal integral damage model to the mortar. A new cyclic dissymmetry (material) parameter enables the reduction in compressive strength after tensile loading to be modelled with modularity. Three-dimensional computations of the concrete mesostructures can then be performed for real-shaped aggregates without compromise, i.e., accounting for anisotropic damage of the mortar phase and key feature of the present work, dealing with alternate (cyclic) loading. The damage constitutive equations and their numerical implementation offer robustness up to high levels of damage and induced anisotropy in three-dimensional specimens. Concrete heterogeneities enhance damage evolution and induced anisotropic behaviour within the mortar matrix.</div></div>","PeriodicalId":50626,"journal":{"name":"Computers & Structures","volume":"309 ","pages":"Article 107650"},"PeriodicalIF":4.4,"publicationDate":"2025-01-30","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143125240","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":"On invariant and locking-free formulations for planar arbitrarily curved beams with Timoshenko-Ehrenfest beam model and peridynamic differential operator","authors":"Duc Van Nguyen , Pana Suttakul , Minh Ngoc Nguyen , Erdogan Madenci , Tinh Quoc Bui , Jaroon Rungamornrat , Duy Vo","doi":"10.1016/j.compstruc.2025.107658","DOIUrl":"10.1016/j.compstruc.2025.107658","url":null,"abstract":"<div><div>This study presents three formulations for the static, dynamic, and fracture simulations of planar arbitrarily curved beams using kinematic assumptions of Timoshenko-Ehrenfest beam model and peridynamic differential operator (PDDO). Displacements of the beam axis and rotation of the cross-section are considered as unknowns of the kinematic description. Two variants of the equations of motion are derived by means of the principle of virtual work, i.e., the first form is expressed in terms of cross-sectional stress resultants, whereas the second variant is written in teams of kinematic unknowns. PDDO is then incorporated into the two variants of the equations of motion and the principle of virtual work to convert them from differential to integral expressions. The driving force behind developing three formulations is to address the critical deficiency of literature, i.e., an invariant and locking-free PD beam formulation for the analysis of beams having complex geometry is not yet proposed. In this study, the invariant property of the proposed formulations is elucidated by theoretical means, and the locking effects are examined by numerical experiments. Several well-established examples are exhibited to assess the accuracy and robustness of the proposed formulations.</div></div>","PeriodicalId":50626,"journal":{"name":"Computers & Structures","volume":"309 ","pages":"Article 107658"},"PeriodicalIF":4.4,"publicationDate":"2025-01-30","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143125244","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":"Shape monitoring of morphing wing structures using the inverse Finite Element Method","authors":"Vincenzo Biscotti, Rinto Roy, Marco Gherlone","doi":"10.1016/j.compstruc.2025.107652","DOIUrl":"10.1016/j.compstruc.2025.107652","url":null,"abstract":"<div><div>This work presents a closed-loop control strategy for morphing wing structures where the feedback originates from monitoring the actual deformed shape of the morphed skin. The approach is based on the inverse Finite Element Method (iFEM), able to reconstruct the displacement field of a structure by minimizing, in a least squares sense, the error between the analytical strains and those experimentally measured in some discrete locations. Once the actual shape has been reconstructed, the actuation loads required to achieve the target shape are computed. The iFEM-based control strategy is assessed numerically on the example problem of a wing segment whose trailing-edge camber is modified via the morphing strategy. Actuation loads are represented by concentrated forces or by a distributed pressure, the effect of aerodynamic loads is taken into account, and strain data are measured on the top and bottom morphing skin. The results show accurate convergence to the target shape, thus demonstrating the potential of the proposed control-loop strategy.</div></div>","PeriodicalId":50626,"journal":{"name":"Computers & Structures","volume":"309 ","pages":"Article 107652"},"PeriodicalIF":4.4,"publicationDate":"2025-01-30","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143125235","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":"Buckling analysis of structures with local abnormality using non-uniform spline finite strip method","authors":"Hao Yu, Pizhong Qiao","doi":"10.1016/j.compstruc.2024.107597","DOIUrl":"10.1016/j.compstruc.2024.107597","url":null,"abstract":"<div><div>Significance of structural components with local abnormality in buckling analysis has drawn considerable interest from researchers. A versatile and effective non-uniform spline finite strip method (N-u SFSM) is developed to allow for mesh refinement in local zones, enabling a comprehensive analysis of buckling characteristics of structures with local abnormality. The inclusion of non-uniform spline functions facilitates the precise modeling and refinement of localized abnormal (e.g., damage or geometric change) regions, leading to more accurate forecasts of fluctuations in buckling modes compared to the original finite strip method. The convergence and validation studies confirm the capability of N-u SFSM to accurately predict buckling behaviors in structures with local abnormality. The computational efficiency is significantly enhanced through the spline knot reduction achieved by optimizing the spline interpolation points, resulting in time and resource savings. The versatility of the N-u SFSM is demonstrated through the successful applications in various scenarios involving the locally-damaged channels, plates, and cylinders with differing damage sizes, types, locations, degrees, numbers, and shapes. The results confirm that the developed N-u SFSM, as a highly efficient and practical numerical technique, can accurately predict the buckling behaviors of both intact and damaged structures exhibiting local abnormalities, thereby providing valuable insights for structural analysis and design.</div></div>","PeriodicalId":50626,"journal":{"name":"Computers & Structures","volume":"307 ","pages":"Article 107597"},"PeriodicalIF":4.4,"publicationDate":"2025-01-15","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142804465","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":"Design of compliant thermal actuators using topology optimization involving design-dependent thermal convection and pressure load","authors":"Shuya Onodera , Takayuki Yamada","doi":"10.1016/j.compstruc.2024.107600","DOIUrl":"10.1016/j.compstruc.2024.107600","url":null,"abstract":"<div><div>This study presents a topology optimization method for thermal actuators that accounts for boundary conditions influenced by variables such as thermal convection and pressure load. Thermal actuators with gripper-like designs are essential for handling hot and brittle materials. The objective of this study is to design actuator shapes that achieve an optimal balance between flexibility and stiffness in high-temperature environments. Unlike previous studies that consider load conditions imposed on fixed boundaries within the design domain, this research introduces a high-temperature fluid as the driving source and employs a novel approach to boundary condition setting by integrating fictitious physical problems. This approach allows for the precise specification of various boundary conditions across multiple domains. A weighted-sum method is applied to optimize three objective functions related to deformability, stiffness, and thermal diffusibility of the actuators. To address the issue of excessively thin structures compromising deformability and the poor convergence of the optimization process, stress constraints based on the optimization history are introduced. The proposed method is validated through numerical examples, demonstrating improvements in structural deformability while controlling deformation on the target surface plane. The numerical results confirm that the objective function decreases and stress is suppressed, verifying the effectiveness of the proposed approach.</div></div>","PeriodicalId":50626,"journal":{"name":"Computers & Structures","volume":"307 ","pages":"Article 107600"},"PeriodicalIF":4.4,"publicationDate":"2025-01-15","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142804466","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}