Pub Date : 2026-01-15DOI: 10.1016/j.compstruc.2026.108109
F. Zanelli, A. Galimberti, N. Debattisti, S. Negri, G. Tomasini
Condition monitoring is becoming an essential tool in the railway industry, increasing the efficiency of vehicle maintenance. This is particularly critical in the case of freight trains since most wagons are still not equipped with sensors or wired. The proposed research work aims at developing a monitoring system with two main purposes: to apply to different wagon typologies and to be an integrated solution for the monitoring of the different mechanical subsystems of the vehicle. The focus is put on the energy harvesting to suitably power supply the wireless sensor nodes and on the performance enhancement for the identification of possible faults in the braking plant. The design of the monitoring system has been driven by an empirical model of the braking plant realized also to support the analysis of experimental data collected by the monitoring devices in the diagnostic stage. Moreover, a CFD analysis was performed to optimize energy harvesters (i.e. micro wind turbines) positioning on the wagon to enhance their efficiency. In the paper, it is shown how the numerical tools allowed to suitably design the wireless monitoring system, which is adopted in an experimental campaign aiming at collecting a database for the validation of condition monitoring algorithms.
{"title":"A modular and integrated on-board system for freight train condition monitoring: design approach and testing","authors":"F. Zanelli, A. Galimberti, N. Debattisti, S. Negri, G. Tomasini","doi":"10.1016/j.compstruc.2026.108109","DOIUrl":"10.1016/j.compstruc.2026.108109","url":null,"abstract":"<div><div>Condition monitoring is becoming an essential tool in the railway industry, increasing the efficiency of vehicle maintenance. This is particularly critical in the case of freight trains since most wagons are still not equipped with sensors or wired. The proposed research work aims at developing a monitoring system with two main purposes: to apply to different wagon typologies and to be an integrated solution for the monitoring of the different mechanical subsystems of the vehicle. The focus is put on the energy harvesting to suitably power supply the wireless sensor nodes and on the performance enhancement for the identification of possible faults in the braking plant. The design of the monitoring system has been driven by an empirical model of the braking plant realized also to support the analysis of experimental data collected by the monitoring devices in the diagnostic stage. Moreover, a CFD analysis was performed to optimize energy harvesters (i.e. micro wind turbines) positioning on the wagon to enhance their efficiency. In the paper, it is shown how the numerical tools allowed to suitably design the wireless monitoring system, which is adopted in an experimental campaign aiming at collecting a database for the validation of condition monitoring algorithms.</div></div>","PeriodicalId":50626,"journal":{"name":"Computers & Structures","volume":"321 ","pages":"Article 108109"},"PeriodicalIF":4.8,"publicationDate":"2026-01-15","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145961764","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-15DOI: 10.1016/j.compstruc.2026.108104
Jesús Flores Escribano, Eduardo Salete Casino, Juan José Benito Muñoz, Eduardo Conde López
Local surface sedimentary structures are made up of stacking strata. Layers of rock or sediment characterized by certain lithological properties that distinguish them from adjacent layers, separated by bedding planes here called interfaces. These strata are typically parallel, but their geometry can be complex, with inclined layers of variable thicknesses produced by tectonic movements or erosion processes. In this work, a formulation is proposed for the treatment of interfaces in seismic wave propagation problems with 3D domains using generalized finite difference method. This capacity is of great relevance, since it is the complexity of the soil structures what makes this meshless method interesting in comparison to the computationally efficient finite difference method. To validate the proposal, a set of examples are solved and their results are compared to those derived from analytical expressions or through finite element method models.
{"title":"New treatment of interfaces for 3D seismic wave propagation problem using generalized finite difference method","authors":"Jesús Flores Escribano, Eduardo Salete Casino, Juan José Benito Muñoz, Eduardo Conde López","doi":"10.1016/j.compstruc.2026.108104","DOIUrl":"10.1016/j.compstruc.2026.108104","url":null,"abstract":"<div><div>Local surface sedimentary structures are made up of stacking strata. Layers of rock or sediment characterized by certain lithological properties that distinguish them from adjacent layers, separated by bedding planes here called interfaces. These strata are typically parallel, but their geometry can be complex, with inclined layers of variable thicknesses produced by tectonic movements or erosion processes. In this work, a formulation is proposed for the treatment of interfaces in seismic wave propagation problems with 3D domains using generalized finite difference method. This capacity is of great relevance, since it is the complexity of the soil structures what makes this meshless method interesting in comparison to the computationally efficient finite difference method. To validate the proposal, a set of examples are solved and their results are compared to those derived from analytical expressions or through finite element method models.</div></div>","PeriodicalId":50626,"journal":{"name":"Computers & Structures","volume":"321 ","pages":"Article 108104"},"PeriodicalIF":4.8,"publicationDate":"2026-01-15","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145976685","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 paper details a new implicit gradient-enhanced damage model designed for simulating fracture in rubber-like materials under large or finite deformations. This model extends previous work by incorporating a rate-dependent crack propagation feature and reformulating the theory for finite strain contexts using a neo-Hookean hyperelastic model and an energy limiter concept. The methodology aims to create a stable and mesh-insensitive numerical tool for fracture analysis, which is implemented and solved using the Finite Element Method (FEM) with a staggered algorithm. Numerical examples validate the model’s accuracy, confirming its ability to predict complex crack growth and demonstrating the importance of the rate-dependent term for stabilizing simulations at large strains.
{"title":"A finite deformation energy limiter-based rate-dependent gradient damage model for fracture analysis","authors":"Tinh Quoc Bui , Hung Thanh Tran , Jaroon Rungamornrat","doi":"10.1016/j.compstruc.2025.108096","DOIUrl":"10.1016/j.compstruc.2025.108096","url":null,"abstract":"<div><div>This paper details a new implicit gradient-enhanced damage model designed for simulating fracture in rubber-like materials under large or finite deformations. This model extends previous work by incorporating a rate-dependent crack propagation feature and reformulating the theory for finite strain contexts using a neo-Hookean hyperelastic model and an energy limiter concept. The methodology aims to create a stable and mesh-insensitive numerical tool for fracture analysis, which is implemented and solved using the Finite Element Method (FEM) with a staggered algorithm. Numerical examples validate the model’s accuracy, confirming its ability to predict complex crack growth and demonstrating the importance of the rate-dependent term for stabilizing simulations at large strains.</div></div>","PeriodicalId":50626,"journal":{"name":"Computers & Structures","volume":"321 ","pages":"Article 108096"},"PeriodicalIF":4.8,"publicationDate":"2026-01-09","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145924954","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-08DOI: 10.1016/j.compstruc.2026.108098
Ran Zheng , Bing Yi , Gil Ho Yoon , Wenlong Liu , Long Liu , Xiang Peng
Although additive manufacturing has advantages in the fabrication of complicated structures and has been applied in many fields, large-scale structures often need to be partitioned into smaller components to comply with the size limitations of the printer, which compromises their overall structural performance. This paper presents a two-scale concurrent topology optimization method for multiple assembled structures, which can be fabricated with an additive manufacturing machine under maximum size limitations and further assembled via conventional joining processes. At the macroscale, the topology design of the macroscale structure and the partitioning of the overall structure into multiple components are realized by incorporating component size constraints into the Solid Isotropic Material with Penalization (SIMP) topology optimization framework. At the microscale, the topology of the self-connected microstructure unit located in the macroscale components and the bolted joint positions of the assembled microstructure unit located in the macroscale joints between different components are optimized based on the homogenization method. The smooth connection between the self-connected microstructure and the assembled microstructure is ensured by a geometric constraint. Finally, several numerical examples and a printing example are provided to illustrate the effectiveness of the proposed method, and the effects of some design parameters on the optimization results are analyzed.
{"title":"Multiscale concurrent topology optimization for large-scale assembled structures","authors":"Ran Zheng , Bing Yi , Gil Ho Yoon , Wenlong Liu , Long Liu , Xiang Peng","doi":"10.1016/j.compstruc.2026.108098","DOIUrl":"10.1016/j.compstruc.2026.108098","url":null,"abstract":"<div><div>Although additive manufacturing has advantages in the fabrication of complicated structures and has been applied in many fields, large-scale structures often need to be partitioned into smaller components to comply with the size limitations of the printer, which compromises their overall structural performance. This paper presents a two-scale concurrent topology optimization method for multiple assembled structures, which can be fabricated with an additive manufacturing machine under maximum size limitations and further assembled via conventional joining processes. At the macroscale, the topology design of the macroscale structure and the partitioning of the overall structure into multiple components are realized by incorporating component size constraints into the Solid Isotropic Material with Penalization (SIMP) topology optimization framework. At the microscale, the topology of the self-connected microstructure unit located in the macroscale components and the bolted joint positions of the assembled microstructure unit located in the macroscale joints between different components are optimized based on the homogenization method. The smooth connection between the self-connected microstructure and the assembled microstructure is ensured by a geometric constraint. Finally, several numerical examples and a printing example are provided to illustrate the effectiveness of the proposed method, and the effects of some design parameters on the optimization results are analyzed.</div></div>","PeriodicalId":50626,"journal":{"name":"Computers & Structures","volume":"321 ","pages":"Article 108098"},"PeriodicalIF":4.8,"publicationDate":"2026-01-08","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145924952","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-08DOI: 10.1016/j.compstruc.2026.108105
Heng Zhang , Gaoyao Feng , Dan Huang , Ding Chen , Mengxin Liu
In this paper, a new fully coupled peridynamic thermo-mechanical interface model is proposed for the bimaterial fracture analysis with the perfect or imperfect material interface. First, a fully coupled peridynamic thermo-mechanical model is presented, in which a novel formation of the peridynamic free energy density is defined, and the peridynamic mechanical deformation model and heat conduction model with their coupling terms are systematically deduced. Then, the peridynamic thermo-mechanical coupled interface models are proposed for both perfect and imperfect interfaces, and corresponding bond failure criteria are given for the thermally induced interface fracture analysis. Four examples, i.e. a square plate under the shock tension loading, a bimaterial plate subjected to the heat flux, double layers beams (DLB) under the thermal loading, and three-point bending (TPB) tests of a bimaterial specimen, are analyzed with the proposed model and compared to the FEM and experimental solutions. The results show that the proposed model can well capture interaction characteristics of the heat conduction, mechanical deformation and crack growth of the bimaterial with the perfect or imperfect interface.
{"title":"Thermo-mechanical fracture analysis of the bimaterials with a fully coupled peridynamic interface model","authors":"Heng Zhang , Gaoyao Feng , Dan Huang , Ding Chen , Mengxin Liu","doi":"10.1016/j.compstruc.2026.108105","DOIUrl":"10.1016/j.compstruc.2026.108105","url":null,"abstract":"<div><div>In this paper, a new fully coupled peridynamic thermo-mechanical interface model is proposed for the bimaterial fracture analysis with the perfect or imperfect material interface. First, a fully coupled peridynamic thermo-mechanical model is presented, in which a novel formation of the peridynamic free energy density is defined, and the peridynamic mechanical deformation model and heat conduction model with their coupling terms are systematically deduced. Then, the peridynamic thermo-mechanical coupled interface models are proposed for both perfect and imperfect interfaces, and corresponding bond failure criteria are given for the thermally induced interface fracture analysis. Four examples, i.e. a square plate under the shock tension loading, a bimaterial plate subjected to the heat flux, double layers beams (DLB) under the thermal loading, and three-point bending (TPB) tests of a bimaterial specimen, are analyzed with the proposed model and compared to the FEM and experimental solutions. The results show that the proposed model can well capture interaction characteristics of the heat conduction, mechanical deformation and crack growth of the bimaterial with the perfect or imperfect interface.</div></div>","PeriodicalId":50626,"journal":{"name":"Computers & Structures","volume":"321 ","pages":"Article 108105"},"PeriodicalIF":4.8,"publicationDate":"2026-01-08","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145924953","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-05DOI: 10.1016/j.compstruc.2025.108095
Zhiyuan Tong, Mauricio Ponga
We present a novel meshfree method for accurate and efficient solutions to solid mechanics problems involving large deviatoric deformation, fracture, and fragmentation. Central to the approach is a newly developed meshfree shape function that satisfies the Kronecker delta property, exhibits first-order consistency, is non-negative, and achieves smoothness. The formulation employs nodal integration which is desirable for problems with large topological changes. The combination of nodal integration and the Kronecker delta property leads to a naturally diagonal explicit equilibrium equation even in the presence of complex boundary and contact conditions. A key innovation is the introduction of shadow nodes, which, in conjunction with a local triangle removal strategy, enables the seamless handling of complex geometries and evolving discontinuities without explicit boundary representations. The method demonstrates excellent convergence and high accuracy across a range of linear and nonlinear benchmark problems. Its robustness and versatility are further illustrated through challenging simulations involving extreme fracture and fragmentation.
{"title":"A new meshfree method for accurate and efficient solutions to solid mechanics problems involving large deviatoric deformation, fracture, and fragmentation","authors":"Zhiyuan Tong, Mauricio Ponga","doi":"10.1016/j.compstruc.2025.108095","DOIUrl":"10.1016/j.compstruc.2025.108095","url":null,"abstract":"<div><div>We present a novel meshfree method for accurate and efficient solutions to solid mechanics problems involving large deviatoric deformation, fracture, and fragmentation. Central to the approach is a newly developed meshfree shape function that satisfies the Kronecker delta property, exhibits first-order consistency, is non-negative, and achieves <span><math><msup><mrow><mi>C</mi></mrow><mn>1</mn></msup></math></span> smoothness. The formulation employs nodal integration which is desirable for problems with large topological changes. The combination of nodal integration and the Kronecker delta property leads to a naturally diagonal explicit equilibrium equation even in the presence of complex boundary and contact conditions. A key innovation is the introduction of shadow nodes, which, in conjunction with a local triangle removal strategy, enables the seamless handling of complex geometries and evolving discontinuities without explicit boundary representations. The method demonstrates excellent convergence and high accuracy across a range of linear and nonlinear benchmark problems. Its robustness and versatility are further illustrated through challenging simulations involving extreme fracture and fragmentation.</div></div>","PeriodicalId":50626,"journal":{"name":"Computers & Structures","volume":"321 ","pages":"Article 108095"},"PeriodicalIF":4.8,"publicationDate":"2026-01-05","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145902509","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-02DOI: 10.1016/j.compstruc.2025.108086
Nasrin Talebi, Knut Andreas Meyer, Magnus Ekh
Simulation of many loading cycles with traditional time-domain material models, requiring discretization of each cycle with several time steps, can result in high computational cost. One effective approach to speed up cyclic simulations is employing cycle-domain material models. Finite element simulations of rails subjected to many wheel passages are a relevant application of such models. Proposing a per-cycle evolution equation for plastic strains in cycle-domain models is, however, a challenge. To address this, we investigate the feasibility and accuracy of using machine learning models as tools for formulating such an equation. Specifically, we enforce our knowledge from constitutive modeling for elasticity and formulate the evolution law by employing feed-forward neural networks with different inputs, as well as symbolic regression to discover an interpretable expression. Training, validation, and test data have been generated using a cyclic time-domain plasticity model considering pulsating uniaxial stress loadings with constant and variable strain ranges. The obtained results demonstrate the potential of cycle-domain plasticity modeling using both uninterpretable and interpretable data-driven machine learning as an alternative to time-domain material modeling. Furthermore, both approaches have revealed reasonably good extrapolation performance beyond the training regime.
{"title":"Cycle-domain plasticity modeling using neural networks and symbolic regression","authors":"Nasrin Talebi, Knut Andreas Meyer, Magnus Ekh","doi":"10.1016/j.compstruc.2025.108086","DOIUrl":"10.1016/j.compstruc.2025.108086","url":null,"abstract":"<div><div>Simulation of many loading cycles with traditional time-domain material models, requiring discretization of each cycle with several time steps, can result in high computational cost. One effective approach to speed up cyclic simulations is employing cycle-domain material models. Finite element simulations of rails subjected to many wheel passages are a relevant application of such models. Proposing a per-cycle evolution equation for plastic strains in cycle-domain models is, however, a challenge. To address this, we investigate the feasibility and accuracy of using machine learning models as tools for formulating such an equation. Specifically, we enforce our knowledge from constitutive modeling for elasticity and formulate the evolution law by employing feed-forward neural networks with different inputs, as well as symbolic regression to discover an interpretable expression. Training, validation, and test data have been generated using a cyclic time-domain plasticity model considering pulsating uniaxial stress loadings with constant and variable strain ranges. The obtained results demonstrate the potential of cycle-domain plasticity modeling using both uninterpretable and interpretable data-driven machine learning as an alternative to time-domain material modeling. Furthermore, both approaches have revealed reasonably good extrapolation performance beyond the training regime.</div></div>","PeriodicalId":50626,"journal":{"name":"Computers & Structures","volume":"321 ","pages":"Article 108086"},"PeriodicalIF":4.8,"publicationDate":"2026-01-02","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145884327","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-02DOI: 10.1016/j.compstruc.2025.108094
Yopi P. Oktiovan , Francesco Messali , Bora Pulatsu , Satyadhrik Sharma , José V. Lemos , Jan G. Rots
This paper presents a cyclic joint constitutive model within a Distinct Element Method framework to simulate the in-plane response of unreinforced masonry structures. The model combines multi-surface failure criteria, including tensile cut-off, Coulomb friction, and an elliptical compression cap. It incorporates exponential softening, a unified damage scalar for stiffness degradation, and a hardening–softening law for compression. Shear-induced dilatancy is captured via an uplift-correction mechanism with an exponential dilatancy-decay law, while stiffness degradation governs energy dissipation. The model is validated at both material and structural scales. Material-level simulations of cyclic compression and shear tests show close agreement with experimental data. Structural-scale validation on full-height calcium-silicate walls under combined compression and cyclic lateral loading demonstrates the ability to reproduce rocking-dominated, shear-dominated, and hybrid failure mechanisms. The model successfully replicated global hysteretic force–drift loops, capturing stiffness decay and energy dissipation, as well as local failures like cracking, sliding, and toe crushing. The model also reproduced the drift-dependent transition from rocking to friction-controlled sliding, a key mechanism for earthquake assessment. By integrating these features into a single, efficient framework, the proposed constitutive model provides a robust tool for evaluating seismic performance and conserving heritage.
{"title":"Cyclic constitutive model for masonry joint damage and energy dissipation using the distinct element method","authors":"Yopi P. Oktiovan , Francesco Messali , Bora Pulatsu , Satyadhrik Sharma , José V. Lemos , Jan G. Rots","doi":"10.1016/j.compstruc.2025.108094","DOIUrl":"10.1016/j.compstruc.2025.108094","url":null,"abstract":"<div><div>This paper presents a cyclic joint constitutive model within a Distinct Element Method framework to simulate the in-plane response of unreinforced masonry structures. The model combines multi-surface failure criteria, including tensile cut-off, Coulomb friction, and an elliptical compression cap. It incorporates exponential softening, a unified damage scalar for stiffness degradation, and a hardening–softening law for compression. Shear-induced dilatancy is captured via an uplift-correction mechanism with an exponential dilatancy-decay law, while stiffness degradation governs energy dissipation. The model is validated at both material and structural scales. Material-level simulations of cyclic compression and shear tests show close agreement with experimental data. Structural-scale validation on full-height calcium-silicate walls under combined compression and cyclic lateral loading demonstrates the ability to reproduce rocking-dominated, shear-dominated, and hybrid failure mechanisms. The model successfully replicated global hysteretic force–drift loops, capturing stiffness decay and energy dissipation, as well as local failures like cracking, sliding, and toe crushing. The model also reproduced the drift-dependent transition from rocking to friction-controlled sliding, a key mechanism for earthquake assessment. By integrating these features into a single, efficient framework, the proposed constitutive model provides a robust tool for evaluating seismic performance and conserving heritage.</div></div>","PeriodicalId":50626,"journal":{"name":"Computers & Structures","volume":"321 ","pages":"Article 108094"},"PeriodicalIF":4.8,"publicationDate":"2026-01-02","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145884457","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-02DOI: 10.1016/j.compstruc.2025.108093
Ziling Song , Smriti , Sundararajan Natarajan , Tiantang Yu
Structural shape optimization plays a crucial role in finding aesthetically pleasing shape designs with reasonable mechanical performance. Isogeometric analysis offers a promising approach for such optimization due to its advantage of unifying the design and analysis models. This paper presents a comprehensive shape optimization methodology for free-form surfaces within isogeometric analysis framework, addressing both compliance and buckling problems. The analytical solution of Kirchhoff-Love shell is derived to enable efficient gradient-based optimization. For complex free-form surfaces modeled with multiple non-uniform rational B-spline (NURBS) patches, a gradient-free optimization strategy is employed to ensure robustness. Continuity constraints across multi-patch interfaces are enforced through Nitsche’s method. The proposed method is validated through several numerical examples, demonstrating simultaneous improvement in structural compliance and buckling stability.
{"title":"Multi-patch isogeometric shape optimization of complex free-form surfaces with buckling constraints","authors":"Ziling Song , Smriti , Sundararajan Natarajan , Tiantang Yu","doi":"10.1016/j.compstruc.2025.108093","DOIUrl":"10.1016/j.compstruc.2025.108093","url":null,"abstract":"<div><div>Structural shape optimization plays a crucial role in finding aesthetically pleasing shape designs with reasonable mechanical performance. Isogeometric analysis offers a promising approach for such optimization due to its advantage of unifying the design and analysis models. This paper presents a comprehensive shape optimization methodology for free-form surfaces within isogeometric analysis framework, addressing both compliance and buckling problems. The analytical solution of Kirchhoff-Love shell is derived to enable efficient gradient-based optimization. For complex free-form surfaces modeled with multiple non-uniform rational B-spline (NURBS) patches, a gradient-free optimization strategy is employed to ensure robustness. Continuity constraints across multi-patch interfaces are enforced through Nitsche’s method. The proposed method is validated through several numerical examples, demonstrating simultaneous improvement in structural compliance and buckling stability.</div></div>","PeriodicalId":50626,"journal":{"name":"Computers & Structures","volume":"321 ","pages":"Article 108093"},"PeriodicalIF":4.8,"publicationDate":"2026-01-02","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145884329","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-30DOI: 10.1016/j.compstruc.2025.108092
Chunfa Wang , Libang Hu , Yudong Li , Hengxiao Lu , Yan Li , Han Hu , Zhiqiang Feng
Shell-like elastic structures—such as hemispherical shells and semi-cylinders—exhibit heightened susceptibility to buckling instability under axial compressive loading, a phenomenon particularly critical in confined assemblies, stacked components, or adhesive-bonded structures. In such practical engineering scenarios, frictional contact and adhesion interactions with adjacent objects emerge as pivotal factors influencing buckling behavior. These interfacial forces, stemming from localized contact pressures, confinement effects, or adhesive bonding, induce a complex coupling between contact mechanics and buckling phenomena, thereby fundamentally altering the structural response under compressive stress. To analyze these effects, we present a computational framework that integrates a point-to-segment (PTS) contact formulation and an exponential cohesive zone model for adhesion. This unified framework enables the simultaneous simulation of friction, adhesion, and buckling, including large deformations and sliding. Implemented on an in-house isogeometric analysis platform, the framework is rigorously validated against theoretical, experimental, and numerical benchmarks. Numerical experiments demonstrate its robustness under challenging conditions, revealing key bidirectional couplings: (1) friction and adhesion suppress buckling by resisting compressive stresses within the shell-like structures, thereby increasing the critical buckling load; (2) buckling-induced geometric nonlinearities dynamically alter contact areas and pressure distributions, which in turn modulate interfacial friction and adhesion strengths.
{"title":"Analysis of friction-adhesion-buckling interactions in shell-like elastic structures via an isogeometric point-to-segment contact formulation","authors":"Chunfa Wang , Libang Hu , Yudong Li , Hengxiao Lu , Yan Li , Han Hu , Zhiqiang Feng","doi":"10.1016/j.compstruc.2025.108092","DOIUrl":"10.1016/j.compstruc.2025.108092","url":null,"abstract":"<div><div>Shell-like elastic structures—such as hemispherical shells and semi-cylinders—exhibit heightened susceptibility to buckling instability under axial compressive loading, a phenomenon particularly critical in confined assemblies, stacked components, or adhesive-bonded structures. In such practical engineering scenarios, frictional contact and adhesion interactions with adjacent objects emerge as pivotal factors influencing buckling behavior. These interfacial forces, stemming from localized contact pressures, confinement effects, or adhesive bonding, induce a complex coupling between contact mechanics and buckling phenomena, thereby fundamentally altering the structural response under compressive stress. To analyze these effects, we present a computational framework that integrates a point-to-segment (PTS) contact formulation and an exponential cohesive zone model for adhesion. This unified framework enables the simultaneous simulation of friction, adhesion, and buckling, including large deformations and sliding. Implemented on an in-house isogeometric analysis platform, the framework is rigorously validated against theoretical, experimental, and numerical benchmarks. Numerical experiments demonstrate its robustness under challenging conditions, revealing key bidirectional couplings: (1) friction and adhesion suppress buckling by resisting compressive stresses within the shell-like structures, thereby increasing the critical buckling load; (2) buckling-induced geometric nonlinearities dynamically alter contact areas and pressure distributions, which in turn modulate interfacial friction and adhesion strengths.</div></div>","PeriodicalId":50626,"journal":{"name":"Computers & Structures","volume":"321 ","pages":"Article 108092"},"PeriodicalIF":4.8,"publicationDate":"2025-12-30","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145884326","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}