{"title":"Proper generalized decomposition surrogate modeling with application to the identification of Rayleigh damping parameters","authors":"Clément Vella , Serge Prudhomme","doi":"10.1016/j.compstruc.2025.107826","DOIUrl":"10.1016/j.compstruc.2025.107826","url":null,"abstract":"<div><div>This paper extends the Proper Generalized Decomposition framework to develop a reduced-order model parameterized by Rayleigh damping coefficients. The developed method incorporates damping modes to construct a damped surrogate model effectively. A novel method is introduced for treating the problem in space: during the offline phase, the spatial problem is initially projected onto the subspace spanned by the Ritz vectors of the system to provide an efficient prediction of the spatial modes. The prediction is then refined using a MinRes iterative solver. This two-step, prediction–correction process reduces the computational cost of a full-order solution while improving the accuracy of the reduced model. The resulting Proper Generalized Decomposition surrogate is subsequently employed within a Particle Swarm Optimization algorithm to determine optimal damping coefficients based on a given snapshot. Numerical experiments demonstrate the effectiveness of the developed method.</div></div>","PeriodicalId":50626,"journal":{"name":"Computers & Structures","volume":"315 ","pages":"Article 107826"},"PeriodicalIF":4.4,"publicationDate":"2025-05-29","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144169039","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":"Efficient parametric model order reduction in contact mechanics","authors":"Ganesh S. Pawar, Salil S. Kulkarni","doi":"10.1016/j.compstruc.2025.107838","DOIUrl":"10.1016/j.compstruc.2025.107838","url":null,"abstract":"<div><div>Contact problems are inherently non-linear and present significant computational challenges in simulations. Traditional proper orthogonal decomposition-based non-linear system reduction often proves inefficient due to the complexity of handling full-scale models. This article presents a generalized parametric model order reduction framework tailored for dynamic contact problems involving arbitrarily-shaped inclusions, incorporating various specialized approaches. The proposed framework features a two-tier reduction process: the first tier applies a proper orthogonal decomposition-based model order reduction approach, while the second tier employs either the discrete empirical interpolation method or energy conserving sampling and weighting to address the reduction of non-linear terms. Discrete empirical interpolation method identifies nodes in contact as indices for the reduction of non-linear terms, and energy conserving sampling and weighting highlight contact elements as the most contributing elements to virtual work associated with non-linear forces. These results align with physics of the problem, as in this article, everything is linear except for the contact non-linearity. Two construction strategies are evaluated: one utilizing a single global basis and another using interpolated local bases. The framework leveraging local bases interpolation on the tangent space to the Grassmannian manifold demonstrates better accuracy compared to the global basis approach. The effectiveness of the different approaches within the generalized parametric model order reduction framework is evaluated through standard 2D dynamic contact problems involving perfectly bonded inclusions.</div></div>","PeriodicalId":50626,"journal":{"name":"Computers & Structures","volume":"315 ","pages":"Article 107838"},"PeriodicalIF":4.4,"publicationDate":"2025-05-27","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144137285","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}
Rodolfo Pina-Torres, Dong Zhao, Johannes Storm, Michael Kaliske
{"title":"Time homogenization: An acceleration scheme for phase-field modeling of fatigue","authors":"Rodolfo Pina-Torres, Dong Zhao, Johannes Storm, Michael Kaliske","doi":"10.1016/j.compstruc.2025.107824","DOIUrl":"10.1016/j.compstruc.2025.107824","url":null,"abstract":"<div><div>This study proposes a novel time homogenization scheme designed for phase-field formulations in fatigue fracture analysis. Inspired by the methodologies for evaluating the long-term behavior of asphalt pavements, this study builds upon a phase-field formulation that accounts for material degradation due to fatigue and the Representative Crack Element formulation as an energy split. The novelty of this approach lies in its application of a methodology based on computational homogenization to fracture mechanics. The time homogenization accelerates the simulation by upscaling micro-scale behaviors over extended macro-time periods. Through comparisons between high-fidelity simulations and those employing the homogenization scheme, the study evaluates the potential and accuracy of the approach. Numerical examples demonstrate the model’s effectiveness in capturing crack growth under diverse loading conditions. While the homogenized model approximates damage progression observed in detailed simulations, it exhibits a delay in response during advanced stages of crack propagation, suggesting areas for further refinement. Overall, this research validates time homogenization within phase-field formulations as a tool for fracture analysis. Further model potential is discussed and necessary development for the time homogenization method is concluded, setting the stage for future development and optimization efforts.</div></div>","PeriodicalId":50626,"journal":{"name":"Computers & Structures","volume":"315 ","pages":"Article 107824"},"PeriodicalIF":4.4,"publicationDate":"2025-05-24","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144130813","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":"Vershinin–Bai–Wierzbicki plastic model for mild steel and accurate prediction of structural plastic response and failure behavior","authors":"Ya-Chao Hu , Feng Xi , Feng Liu , Ying-Hua Tan","doi":"10.1016/j.compstruc.2025.107834","DOIUrl":"10.1016/j.compstruc.2025.107834","url":null,"abstract":"<div><div>The isotropic plastic hardening behavior of ductile metals is commonly characterized by the accumulated equivalent plastic strain. Advanced plasticity models further incorporate dependencies on hydrostatic pressure and Lode angle to more accurately represent material behavior under complex stress states. However, such models have rarely been applied to mild steels or structural failure analyses. This study evaluates the applicability of the Vershinin–Bai–Wierzbicki (VBW) model to mild steel, such as Q235 in China. Constitutive parameters are calibrated using a parameter identification approach that integrates inverse modeling with multi-objective optimization. Computational results of the VBW and Mises models are compared, with particular attention to differences in the large-deformation plastic stage. In addition, the VBW model incorporating the Lode angle-modified void growth model (LMVGM) is evaluated for its capability to predict fracture behavior and structural failure responses. Simulation results indicate that plastic deformation and fracture behavior of mild steel under large strains are highly sensitive to stress triaxiality and Lode angle. The conventional Mises model fails to accurately capture plastic evolution influenced by these stress-state parameters, whereas the VBW model, coupled with LMVGM, provides accurate predictions of both large-strain plasticity and fracture responses in mild steel materials and structures.</div></div>","PeriodicalId":50626,"journal":{"name":"Computers & Structures","volume":"315 ","pages":"Article 107834"},"PeriodicalIF":4.4,"publicationDate":"2025-05-22","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144107664","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":"Time-domain dynamic analysis of structures equipped with fractional viscoelastic solid and fluid dampers via improved pseudo-force approach","authors":"Federica Genovese , Giuseppe Muscolino","doi":"10.1016/j.compstruc.2025.107816","DOIUrl":"10.1016/j.compstruc.2025.107816","url":null,"abstract":"<div><div>A numerical method for the time-domain dynamic analysis of structures with viscoelastic energy dissipation dampers, modeled using fractional derivatives, is presented. Two fractional viscoelastic models are considered: the <em>fractional Kelvin-Voigt model</em> and another one referred to here as the <em>fractional simplified Maxwell model</em>, to distinguish it from the widely used fractional Maxwell model, where two fractional derivatives appear. In the context of the co<em>nvolution integral method</em> in state variables, the proposed approach, called the <em>improved pseudo-force method</em> involves: <em>i</em>) discretization of<!--> <!-->fractional derivatives using the <em>Grünwald-Letnikov approximation</em>; <em>ii</em>) piecewise linear interpolation of both the excitation function and the pseudo-force; <em>iii</em>) evaluation of the response in the modal subspace through recursive relations, using operators computed only at the beginning of the procedure; <em>iv</em>) evaluation of the nodal response by the modal superposition method. The accuracy and stability of the method are demonstrated through applications to a Single-Degree-of-Freedom (SDoF) oscillator and a five-story building equipped with viscoelastic dampers.</div></div>","PeriodicalId":50626,"journal":{"name":"Computers & Structures","volume":"315 ","pages":"Article 107816"},"PeriodicalIF":4.4,"publicationDate":"2025-05-22","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144107663","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":"Multi-vehicle responses for high-resolution bridge mode shape identification integrating Kalman filter and compressive sensing","authors":"Yi He , Judy P. Yang","doi":"10.1016/j.compstruc.2025.107837","DOIUrl":"10.1016/j.compstruc.2025.107837","url":null,"abstract":"<div><div>This study introduces a three-step procedure for identifying high-resolution bridge mode shapes using responses from a limited number of test vehicles. First, contact-point displacements are retrieved from the vehicle responses using the generalized Kalman filter with unknown input algorithm. Second, the sparse bridge response matrix, populated with contact-point displacements, is completed using the spatial compressive sensing theory. Third, high-resolution mode shapes are extracted by applying singular value decomposition in the completed response matrix. An illustrative example shows that the first two mode shapes of a 60-m bridge can be well identified using the responses of eight test vehicles, achieving a spatial resolution of 0.5 m. The performance of the procedure is further evaluated by considering practical factors, including bridge boundary conditions, environmental noises, and the number of test vehicles. Additionally, the subtraction strategy has successfully removed the effect of pavement irregularity for mode shape construction. The capability of the procedure for accurately and effectively identifying the high-resolution bridge mode shapes is therefore demonstrated.</div></div>","PeriodicalId":50626,"journal":{"name":"Computers & Structures","volume":"315 ","pages":"Article 107837"},"PeriodicalIF":4.4,"publicationDate":"2025-05-21","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144107662","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":"3D-printed concrete fracture: Effects of cohesive laws, mixes, and print parameters in 3D eXtended FEM","authors":"Faisal Mukhtar","doi":"10.1016/j.compstruc.2025.107822","DOIUrl":"10.1016/j.compstruc.2025.107822","url":null,"abstract":"<div><div>Unlike conventional concrete fractures, few models of 3D-printed concrete (3DPC) fractures have been reported; moreover, systematic validation across diverse tests, materials, and laboratories is lacking. This paper first reviews existing 3DPC fracture simulations against experiments, noting mixed performance in most cases. Additionally, current models often require excessive material parameters that can be difficult to measure and interpret, along with a scarcity of 3D modeling. To address these issues, the paper develops a robust 3D validation framework using the generalized/eXtended finite element method to simulate both interlayer bonds and intralayer fractures in 3DPC/cementitious materials. Four softening models were assessed by capturing linear and nonlinear fracture responses under tensile, bending, and shear tests. Different 3DPC mixtures (plain, fiber-reinforced, and limestone-calcined clay) from existing experimental studies across different laboratories were used for validation, showing good agreement with experimental load–displacement behavior. The influences of layer print interval and nozzle standoff distance from two separate tests were investigated, revealing that the bilinear cohesive model provided the most consistent predictions for plain and fiber-reinforced 3DPC under tension and bending. In contrast, the Park–Paulino–Roesler and Xu–Needleman cohesive models better captured fracture behavior for the limestone-calcined clay mix, while the linear cohesive model was most suitable for shear mode. A fracture in a topologically optimized 3DPC girder was also simulated, illustrating the method’s applicability to large-scale structural elements, the importance of the material mix, and highlighting a scenario where 2D approximations are inapplicable.</div></div>","PeriodicalId":50626,"journal":{"name":"Computers & Structures","volume":"315 ","pages":"Article 107822"},"PeriodicalIF":4.4,"publicationDate":"2025-05-19","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144083776","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}
Yu Sun , Sen Zhang , Mengting Xing , Zhiyong Shi , Pengfei Liu
{"title":"Three-dimensional train-periodic slab track-subgrade dynamics model based on the iterative solution and Green’s function method","authors":"Yu Sun , Sen Zhang , Mengting Xing , Zhiyong Shi , Pengfei Liu","doi":"10.1016/j.compstruc.2025.107784","DOIUrl":"10.1016/j.compstruc.2025.107784","url":null,"abstract":"<div><div>This paper develops an iterative solution model for the efficient and accurate simulation of the dynamics of a three-dimensional (3D) train-periodic slab track-subgrade (TPSTS) system. The entire system is divided into the train-rail subsystem and the periodic slab-subgrade subsystem. An ordinary differential equation (ODE) model of the train-rail system is established, and a step-moving strategy is employed to enhance the computational efficiency. A frequency-domain Green’s function method (GFM) model is employed for modelling the periodic slab-subgrade. Considering the periodicity, symmetry, and attenuation characteristics of the track structure, the Green’s function of the slab-subgrade structure is obtained by conducting harmonic response analysis on a finite element model (FEM) of a quarter slab-subgrade structure. The proposed model is utilized to investigate the vibration characteristics of the high-speed train-track-subgrade system. The results indicate that the model exhibits good convergence, accuracy, and efficiency. The periodicity of the track structure has a significant impact on the statistical values of subgrade displacement and stress amplitudes along the track’s longitudinal direction. At the subgrade surface, the amplitudes of displacement, stress, and acceleration exhibit fluctuations in the lateral direction of the track. However, with increasing depth, these amplitudes generally attenuate, leading to a more uniform distribution.</div></div>","PeriodicalId":50626,"journal":{"name":"Computers & Structures","volume":"315 ","pages":"Article 107784"},"PeriodicalIF":4.4,"publicationDate":"2025-05-19","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144083775","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}
Qiong Liu , Lars Vabbersgaard Andersen , Min Wu , Mingzhong Zhang , Didier Snoeck
{"title":"Numerical investigations of the abrasion behavior of concrete based on a coupled Eulerian–Lagrangian approach","authors":"Qiong Liu , Lars Vabbersgaard Andersen , Min Wu , Mingzhong Zhang , Didier Snoeck","doi":"10.1016/j.compstruc.2025.107808","DOIUrl":"10.1016/j.compstruc.2025.107808","url":null,"abstract":"<div><div>This paper presents numerical investigations of the abrasion behavior of concrete for hydraulic structures considering concrete structural characteristics as well as various hydraulic conditions. Three-dimensional mesoscale models of concrete composed of aggregates, mortar, and interfacial transition zones are developed using in-house Python 2 codes and the commercial finite-element software Abaqus 2021. The coupled Eulerian–Lagrangian approach is employed to simulate the complex interaction effects between the hydrodynamic field and the concrete structure with a focus on abrasion material loss and energy response as functions of the governing hydraulic parameters (flow velocity and sediment concentration). The results indicate that the concrete abrasion behavior is greatly influenced by the flow velocity and sediment concentration, which are highly associated with the kinetic energy in the flow exerted on the concrete structure.</div></div>","PeriodicalId":50626,"journal":{"name":"Computers & Structures","volume":"315 ","pages":"Article 107808"},"PeriodicalIF":4.4,"publicationDate":"2025-05-17","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144071217","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}
Ductho Le , Haidang Phan , Hoai Nguyen , Mauricio D. Sacchi , Lawrence H. Le
{"title":"Guided waves in anisotropic and viscoelastic stratified plates: Application to bone quantitative ultrasound","authors":"Ductho Le , Haidang Phan , Hoai Nguyen , Mauricio D. Sacchi , Lawrence H. Le","doi":"10.1016/j.compstruc.2025.107820","DOIUrl":"10.1016/j.compstruc.2025.107820","url":null,"abstract":"<div><div>This paper presents a numerical investigation into the behavior of ultrasonic guided waves in a multilayered bone plate characterized by anisotropic and viscoelastic properties. A spectral collocation formula for bone structures is developed, enabling accurate computation of three-dimensional dispersion curves, wave structures, and attenuation of guided waves. Additionally, a mode classification scheme is introduced to identify individual modes, facilitating mode selection during bone inversion. Our results are benchmarked against alternative wave computation methods and experimental data to validate their robustness and reliability. The effects of incorporating soft tissue and marrow into the bone models are explored in relation to mode number, trajectory, and attenuation. With its fast computational speed and mode-classification features, the proposed methodology serves as an effective forward modeling tool for enhanced ultrasonic bone assessment.</div></div>","PeriodicalId":50626,"journal":{"name":"Computers & Structures","volume":"315 ","pages":"Article 107820"},"PeriodicalIF":4.4,"publicationDate":"2025-05-12","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143935005","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}