Computers & Structures最新文献

{"title":"Nonlinear eigenvalue solver for spectral element of beam structures: An exponential matrix polynomial approximation with weighted residual method","authors":"Arindam Das ,&nbsp;Avisek Mukherjee ,&nbsp;Kamal Krishna Bera ,&nbsp;Arnab Banerjee","doi":"10.1016/j.compstruc.2025.107962","DOIUrl":"10.1016/j.compstruc.2025.107962","url":null,"abstract":"<div><div>The spectral element method (SEM) is a widely used frequency-domain technique for dynamic structural analysis. However, its frequency-dependent dynamic stiffness matrix leads to a nonlinear/transcendental eigenvalue problem (NLEP) for obtaining natural frequencies and mode shapes. Traditional numerical approaches, such as iterative root-finding and matrix polynomial linearization, are often used to solve NLEP. While linearization is robust, it becomes computationally expensive with increasing degrees of freedom and polynomial order. This study introduces a novel exponential matrix polynomial approximation of the dynamic stiffness matrix, combined with a weighted residual technique to compute polynomial coefficients without trigonometric or hyperbolic functions. The polynomial eigenvalue problem is then transformed into a generalized eigenvalue problem using a matrix pencil. The proposed method efficiently handles NLEP, even in cases with singularities and closely spaced modes. A lower-order matrix polynomial approximation improves computational efficiency while maintaining accuracy, outperforming Lagrange interpolating polynomials. The natural frequencies and mode shapes of thin-walled beams and frame structures with various boundary conditions are determined using the present NLEP solver, showing a close match with FEM results. However, in FEM the computational time increases exponentially with higher modes, whereas in SEM solved via NLEP, the increase is only marginal.</div></div>","PeriodicalId":50626,"journal":{"name":"Computers & Structures","volume":"318 ","pages":"Article 107962"},"PeriodicalIF":4.8,"publicationDate":"2025-09-15","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145060303","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":"Robust path-following and branch-switching in isogeometric nonlinear bifurcation analysis of variable angle tow panels with cutouts under compression","authors":"Xiaodong Chen","doi":"10.1016/j.compstruc.2025.107948","DOIUrl":"10.1016/j.compstruc.2025.107948","url":null,"abstract":"<div><div>In this paper, an isogeometric nonlinear analysis framework integrating robust path-following and branch-switching techniques is specifically developed to investigate the nonlinear bifurcation behaviors of variable angle tow composite panels containing cutouts under compressive loads. The framework integrates Reddy’s third-order shear deformation theory with von Kármán’s nonlinearity into an isogeometric analysis formulation. The inherent <span><math><msup><mi>C</mi><mrow><mn>1</mn></mrow></msup></math></span> continuity requirement of Reddy’s plate model is naturally satisfied via non-uniform rational B-splines basis functions. Geometric discontinuities arising from complex cutouts are efficiently addressed using the finite cell method with adaptive quadrature refinement. Nonlinear equilibrium equations under force- or displacement-controlled edge loads are first derived from the principle of virtual work and subsequently solved through robust path-following and branch-switching algorithms. The framework distinguishes itself through three core capabilities: (i) handling snap-through and snap-back instabilities using force- or displacement-controlled arc-length scheme; (ii) locating singular points with quadratic convergence via force- or displacement-controlled pinpointing scheme; and (iii) switching to secondary branches at bifurcation points without artificial perturbations. Its potential applicability can extend to more complex shell structures; however, the present study focuses exclusively on plate structures. The effectiveness and robustness of the proposed framework are validated against finite element solutions from ABAQUS. Effects of load condition, hole size and fiber angle on the nonlinear bifurcation behaviors of variable-stiffness composite panels with cutouts loaded in compression are also discussed in numerical examples. Results demonstrate that favourable variable-stiffness configurations retain their beneficial load redistribution mechanisms even in the presence of cutouts. These findings may provide valuable insights for the design of perforated thin-walled structures.</div></div>","PeriodicalId":50626,"journal":{"name":"Computers & Structures","volume":"318 ","pages":"Article 107948"},"PeriodicalIF":4.8,"publicationDate":"2025-09-13","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145049981","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 CNN-based crack propagation prediction in peridynamic concrete models under freeze-thaw cycles","authors":"Feng Nie ,&nbsp;Zhengzheng Wang ,&nbsp;Linfeng Liu ,&nbsp;Huili Wang ,&nbsp;Jiang Lin","doi":"10.1016/j.compstruc.2025.107959","DOIUrl":"10.1016/j.compstruc.2025.107959","url":null,"abstract":"<div><div>Freeze-thaw damage is a critical factor compromising the durability of concrete, rendering the prediction of crack evolution under freeze–thaw conditions essential for evaluating concrete service life. In this study, the relationship between pore frost heaving force and the frost heaving force state is derived, resulting in a peridynamic formulation for the concrete freeze–thaw problem. A three-dimensional convolutional neural network (3D CNN) prediction model, driven by peridynamic theory, is developed. The influence of porosity and pore frost heaving force on the freeze–thaw performance of concrete is systematically analyzed. Based on the proposed model, crack damage predictions over 20 freeze–thaw cycles are carried out. The results indicate that freeze–thaw-induced damage in concrete increases with higher porosity and greater pore frost heaving force. Furthermore, to enhance the accuracy of the 3D CNN in capturing the underlying physical mechanisms, it is recommended to appropriately reduce the number of pooling layers. The developed prediction model demonstrates excellent long-term prediction capability, achieving an accuracy exceeding 92.3%. Compared with traditional methods, the computational efficiency is improved. This study provides an approach for predicting freeze–thaw damage and the remaining service life of concrete in practical engineering applications.</div></div>","PeriodicalId":50626,"journal":{"name":"Computers & Structures","volume":"318 ","pages":"Article 107959"},"PeriodicalIF":4.8,"publicationDate":"2025-09-12","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145049980","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 2D multiresolution wavelet-based method for advanced simulations of highly transient responses in composite plates","authors":"Dimitris K. Dimitriou,&nbsp;Dimitris A. Saravanos","doi":"10.1016/j.compstruc.2025.107958","DOIUrl":"10.1016/j.compstruc.2025.107958","url":null,"abstract":"<div><div>The efficient and robust simulation of transient responses in laminated plates, which may encompass numerous wave modes, is a challenging computational problem with multiple applications. In this direction, an advanced computational method with additional localization capabilities is developed by employing the multiresolution approximation, resulting in the 2D multiresolution finite wavelet domain (MR-FWD) method. Daubechies scaling and wavelet functions are utilized, formulating a novel set of mass-decoupled multiresolution discretized equations of motion that involve four solution components: the coarse, horizontal fine, vertical fine and diagonal fine solution. The multiresolution discretization is combined with the first-order shear laminated plate theory for the effective modeling of composite plates. Numerical case studies focus on guided wave propagation in three lamination cases of increased complexity, showing the remarkable computational efficiency of the MR-FWD method compared to single-resolution approaches and time-domain spectral finite elements. Most importantly, the advanced localization properties of the proposed method are demonstrated, especially in the multi-wave response of an asymmetric composite plate, where each fine solution captures different wave responses depending on their directivity and wavenumbers. This unique feature can pave the way towards a new numerical analysis framework, in which multiple solution components contribute critical information about the structural response.</div></div>","PeriodicalId":50626,"journal":{"name":"Computers & Structures","volume":"318 ","pages":"Article 107958"},"PeriodicalIF":4.8,"publicationDate":"2025-09-12","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145049979","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 geometrically non-linear beam theory for active beams with arbitrary cross-section","authors":"Bram Seinhorst ,&nbsp;Wouter Hakvoort ,&nbsp;Marijn Nijenhuis","doi":"10.1016/j.compstruc.2025.107946","DOIUrl":"10.1016/j.compstruc.2025.107946","url":null,"abstract":"<div><div>Non-linear active beam theories can be used to model many types of prismatic structures that contain piezoelectric material. In this work, we show that the 3D governing continuum equations of active prismatic structures loaded only at their ends can, without any simplifying assumptions on the stress state or geometry of the cross-section, be decomposed into a large deflection active beam theory and deformation modes that exponentially decay in magnitude from both ends of the prismatic structure. This is a manifestation of Saint-Venant’s principle for active beams, where the contribution of the decaying solutions is only relevant near the ends of the beam and can thus be safely neglected if the structure is significantly longer than the dominant decay length. A finite element discretisation of the cross-section is used to handle arbitrary prismatic geometry. By directly discretising the cross-sectional Hamiltonian, the beam constitutive coefficients of the large deflection beam theory can be found efficiently by solving a sparse system of equations. Furthermore, the dominant decay length of the exponentially decaying solutions can be estimated with limited computational effort. The approach is validated against analytical solutions, other numerical cross-section analysis approaches and the 3D finite element software COMSOL for various validation cases.</div></div>","PeriodicalId":50626,"journal":{"name":"Computers & Structures","volume":"318 ","pages":"Article 107946"},"PeriodicalIF":4.8,"publicationDate":"2025-09-11","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145049975","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 efficient method for concurrent multiscale robust structural topology optimization under loading uncertainties","authors":"DingXin Du ,&nbsp;Dong Wang","doi":"10.1016/j.compstruc.2025.107960","DOIUrl":"10.1016/j.compstruc.2025.107960","url":null,"abstract":"<div><div>An efficient method is proposed for the concurrent robust topology optimization of a continuum structure in multiscale material distributions under the external load uncertainties of direction and magnitude such that the resultant design illustrates a prominent insensitivity to the loading variations. In this context, the loading direction and magnitude uncertainties are described in stochastic forms independently. By means of decomposition of the external load along the essential directions, the variation of the structural compliance is firstly represented by a multivariable quadratic Taylor series expansion. Then both the expectation and variance of the compliance can be evaluated efficiently through the statistical analyses‌. Furthermore, the concurrent design sensitivity analyses are performed readily upon the Taylor series formulations, and the multiscale robust topology optimization can be performed by a gradient-based strategy. Several benchmark examples are employed to demonstrate the feasibility of the proposed method, and the computational cost for the compliance statistical characteristics can be reduced remarkably in comparison with a stochastic simulation method. The obtained robust topology optimization designs show quite different from the conventional deterministic counterparts in both macro- and microscale structures. As the result, the variation of the structural compliance can be significantly reduced under the uncertain loading conditions.</div></div>","PeriodicalId":50626,"journal":{"name":"Computers & Structures","volume":"318 ","pages":"Article 107960"},"PeriodicalIF":4.8,"publicationDate":"2025-09-11","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145049977","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":"Eigenvalue buckling analysis through differential equations built on the beam cross section","authors":"Zaiwei Lin,&nbsp;Marco Morandini","doi":"10.1016/j.compstruc.2025.107945","DOIUrl":"10.1016/j.compstruc.2025.107945","url":null,"abstract":"<div><div>This study investigates the global and local buckling of prismatic beams using a semi-analytical approximation of the equilibrium equations, linearized with respect to a pre-stressed configuration. The cross section is discretized using finite elements, while the nodal displacement is modeled analytically along the beam axis. The resulting differential equation system enables tracking the dependence of general solutions on the pre-stress level. It is shown how above a critical stress threshold, some solutions become periodic. Buckling loads can be determined by assuming a periodic solution with a specific half wavelength, allowing for the solution of a simple eigenvalue problem without discretizing the beam along its length.</div></div>","PeriodicalId":50626,"journal":{"name":"Computers & Structures","volume":"318 ","pages":"Article 107945"},"PeriodicalIF":4.8,"publicationDate":"2025-09-11","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145049978","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":"Active learning Kriging method based on particle swarm optimization for reliability analysis with random and interval hybrid uncertainty","authors":"Cheng Yang ,&nbsp;Qingwei Liang ,&nbsp;Hancheng Huang ,&nbsp;Enrico Zio ,&nbsp;Yuxin Lin ,&nbsp;Shanshan Hu","doi":"10.1016/j.compstruc.2025.107947","DOIUrl":"10.1016/j.compstruc.2025.107947","url":null,"abstract":"<div><div>We propose an active learning Kriging reliability method, based on the particle swarm optimization algorithm, to solve structural reliability assessment problems in which both random variables and parameter interval uncertainty coexist. The method optimizes the selection of optimal training samples by using the U learning function as the optimization objective, combined with search space reduction and domain truncation techniques. An error-based stopping criterion is employed to ensure termination of the algorithm. The effectiveness and feasibility of the proposed method are validated through five specific examples, whose results demonstrate the capability of the method to improve accuracy and computational efficiency in the reliability assessment.</div></div>","PeriodicalId":50626,"journal":{"name":"Computers & Structures","volume":"318 ","pages":"Article 107947"},"PeriodicalIF":4.8,"publicationDate":"2025-09-11","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145049976","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 computational strategy for enhanced nonlinear structural stability analysis in Abaqus","authors":"Jiajia Shen ,&nbsp;Kai Wang ,&nbsp;Yibin Fu ,&nbsp;Chaofeng Lü","doi":"10.1016/j.compstruc.2025.107943","DOIUrl":"10.1016/j.compstruc.2025.107943","url":null,"abstract":"<div><div>The global pursuit of net-zero emissions intensifies the demand for sustainable and lightweight structures, heightening challenges associated with nonlinear buckling instabilities. Simultaneously, a paradigm shift is emerging: strategically leveraging these instabilities to unlock unprecedented functionalities in areas like soft robotics, 4D printing, and flexible electronics. This convergence creates a critical need for advanced, automated computational tools capable of systematic nonlinear bifurcation analysis. While commercial software like <span>Abaqus</span> provides robust equilibrium path tracing, it lacks integrated capabilities for explicit critical point pin-pointing and – crucially – automated branch switching at bifurcation points without resorting to symmetry-breaking imperfections. Addressing this significant gap, we introduce a novel computational strategy enabling rigorous bifurcation analysis directly within <span>Abaqus</span>. The core methodological innovation is a probing-based technique that exploits the critical eigenvector’s topology and multi-stability principles near bifurcations to traverse secondary equilibrium paths while rigorously preserving structural symmetry. We demonstrate the framework’s efficacy and broad applicability through three diverse case studies exhibiting complex instabilities: 1) buckling of 2D soft circular ring under uniform inward pressure; 2) interactive buckling of thin-walled rectangular hollow section strut; 3) a concentrically loaded shallow shell roof exhibiting complex buckling instabilities. Results are rigorously verified against analytical solutions or published benchmarks. This work provides a powerful, accessible, and symmetry-preserving approach within the ubiquitous <span>Abaqus</span> environment, enabling systematic explorations of complex post-buckling landscapes. By enabling precise critical point identification and automated branch switching, our framework fundamentally advances the understanding of post-buckling behaviour, leading to more robust lightweight load-bearing structures while unlocking transformative potential for next-generation designs where strategically controlled instabilities enable multifunctional performance.</div></div>","PeriodicalId":50626,"journal":{"name":"Computers & Structures","volume":"318 ","pages":"Article 107943"},"PeriodicalIF":4.8,"publicationDate":"2025-09-09","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145020438","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 thermodynamically-consistent chemo-mechanical-fracture coupled framework for stress corrosion cracking","authors":"Lang Min ,&nbsp;Xiaofei Hu ,&nbsp;Yichao Zhu ,&nbsp;Peng Zhang ,&nbsp;Tong Wang ,&nbsp;Shangtong Yang ,&nbsp;Weian Yao","doi":"10.1016/j.compstruc.2025.107938","DOIUrl":"10.1016/j.compstruc.2025.107938","url":null,"abstract":"<div><div>This article aims to provide a simulation tool for effective predictions over the full-life stress corrosion cracking (SCC) behavior of material and structures, for which experimental observation often proves prohibitively time-consuming. To this end, the SCC dynamics is modeled by means of a thermodynamic process. The proposed theory naturally captures the mechanics’ role in SCC development, that is, the high hydrostatic pressure gradient ahead a corrosion pit/crack enhances the moving tendency of the atoms in solid toward the corrosion environments, while the damage caused by such an atom loss in materials favors the crack advancement. The present theory is numerically realized with a phase-field description of the crack profile. To restore the mass transfer behavior near the smeared boundary, an equivalent sink term is adopted in this model. For its engineering predictability, a general strategy for parameter calibration is proposed and validated against experimental results of the crack growth rate (CGR), whose acquisition is far more feasible than that from the long transition period from pit to crack. Two cases bearing clear engineering origin are then studied with the calibrated model, and life-span predictions at a magnitude of years can be made.</div></div>","PeriodicalId":50626,"journal":{"name":"Computers & Structures","volume":"318 ","pages":"Article 107938"},"PeriodicalIF":4.8,"publicationDate":"2025-09-08","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145010727","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}