Computers & Structures最新文献_第7页

An extended Generalized beam theory formulation for large deformation analyses of perforated Thin-Walled members 多孔薄壁构件大变形分析的广义梁理论扩展公式

IF 4.8 2区工程技术 Q1 COMPUTER SCIENCE, INTERDISCIPLINARY APPLICATIONS

Computers & Structures

Pub Date : 2026-01-15 Epub Date: 2025-12-02 DOI: 10.1016/j.compstruc.2025.108052

Liping Duan , Ji Miao , Jincheng Zhao

This paper presents a geometrically exact beam finite element based on the classic Generalized Beam Theory (GBT) for large-deformation analyses of perforated thin-walled sections. The proposed formulation enables the beam to undergo arbitrary cross-sectional deformations, such as in-plane distortion, out-of-plane warping, transverse bending, and extension of the cross-sectional walls. The key novelty of the proposed approach lies in the extension/improvement of a previous GBT formulation for handling linear buckling analyses of perforated thin-section beams, proposed by the authors Duan et al. (2022), to the realm of large-deformation analysis. The geometrically exact kinematic description is built by enriching the cross-section kinematics of a spatially rotated Euler-Bernoulli/Vlasov thin-section beam with some hierarchic deformation modes and using the enrichment functions constructed in terms of signed distance functions to approximate cross-section discontinuities caused by holes. The Neo-Hookean hyperelastic material model is imposed on the proposed beam FE formulation, where an iterative procedure is performed to ensure the plane stress assumption for the stress fields in cross-section walls. The discrete equations are obtained by using the displacement Galerkin method. Some illustrative examples are presented to demonstrate the potential of the proposed beam FE, and the results demonstrate its capability to perform the GBT modal decomposition of buckled configurations of perforated thin-section members.

本文基于经典的广义梁理论（GBT），提出了一种用于多孔薄壁截面大变形分析的几何精确梁有限元。所提出的公式使梁能够经受任意截面变形，如面内变形、面外翘曲、横向弯曲和截面壁的延伸。所提出方法的关键新颖之处在于，将Duan等人（2022）提出的用于处理穿孔薄截面梁线性屈曲分析的先前GBT公式扩展/改进到大变形分析领域。将空间旋转的Euler-Bernoulli/Vlasov薄截面梁的截面运动学丰富为若干层次变形模式，并利用符号距离函数构造的丰富函数近似孔洞引起的截面不连续，从而建立几何精确的运动学描述。将Neo-Hookean超弹性材料模型应用于所提出的梁有限元公式中，并进行迭代过程以确保截面壁面应力场的平面应力假设。采用位移伽辽金法得到离散方程。本文给出了一些实例来证明所提出的梁有限元分析的潜力，结果表明它能够执行带孔薄壁构件屈曲结构的GBT模态分解。

{"title":"An extended Generalized beam theory formulation for large deformation analyses of perforated Thin-Walled members","authors":"Liping Duan , Ji Miao , Jincheng Zhao","doi":"10.1016/j.compstruc.2025.108052","DOIUrl":"10.1016/j.compstruc.2025.108052","url":null,"abstract":"<div><div>This paper presents a geometrically exact beam finite element based on the classic Generalized Beam Theory (GBT) for large-deformation analyses of perforated thin-walled sections. The proposed formulation enables the beam to undergo arbitrary cross-sectional deformations, such as in-plane distortion, out-of-plane warping, transverse bending, and extension of the cross-sectional walls. The key novelty of the proposed approach lies in the extension/improvement of a previous GBT formulation for handling linear buckling analyses of perforated thin-section beams, proposed by the authors Duan et al. (2022), to the realm of large-deformation analysis. The geometrically exact kinematic description is built by enriching the cross-section kinematics of a spatially rotated Euler-Bernoulli/Vlasov thin-section beam with some hierarchic deformation modes and using the enrichment functions constructed in terms of signed distance functions to approximate cross-section discontinuities caused by holes. The Neo-Hookean hyperelastic material model is imposed on the proposed beam FE formulation, where an iterative procedure is performed to ensure the plane stress assumption for the stress fields in cross-section walls. The discrete equations are obtained by using the displacement Galerkin method. Some illustrative examples are presented to demonstrate the potential of the proposed beam FE, and the results demonstrate its capability to perform the GBT modal decomposition of buckled configurations of perforated thin-section members.</div></div>","PeriodicalId":50626,"journal":{"name":"Computers & Structures","volume":"321 ","pages":"Article 108052"},"PeriodicalIF":4.8,"publicationDate":"2026-01-15","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145657745","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}

引用次数: 0

Experimental and numerical analysis of angle-ply laminated plates based on a refined higher-order theory 基于改进高阶理论的角层合板试验与数值分析

IF 4.8 2区工程技术 Q1 COMPUTER SCIENCE, INTERDISCIPLINARY APPLICATIONS

Computers & Structures

Pub Date : 2026-01-15 Epub Date: 2025-12-23 DOI: 10.1016/j.compstruc.2025.108078

Jie Mei , Shengbo Ling , Yunhao Huang

Accurately predicting the deformation of composite laminates is challenging due to the zig-zag effect in in-plane displacements, which arises from significant stiffness differences between layers. To address these limitations, this study proposes a novel Angle-ply Higher-order Deformation Theory coupled with a computationally efficient refined finite element, to investigate the bending response of angle-ply hybrid composite laminates. The AHDT is established via a global–local concept, utilizing Legendre polynomials to capture the zig-zag effect and interlaminar stress continuity phenomenon. The corresponding triangular element implements this model using C⁰ and C¹ continuities, ensuring high accuracy and computational efficiency. The accuracy and computational efficiency of the model are rigorously assessed through comparisons with three-dimensional elastic solutions and a three-dimensional finite element model. Furthermore, validation is provided via three-point bending tests on carbon/glass fiber hybrid composites, with full-field deformation measured using Digital Image Correlation. Both theoretical and experimental results confirm that the stacking sequence is the predominant factor affecting flexural stiffness, outweighing both material type and ply orientation. A change in the stacking sequence alone can reduce deflection by up to 70.77%. Consequently, this study conclusively shows that a rational design of the stacking sequence, material combination, and ply orientation is critical for optimizing the flexural performance of composite laminates.

由于层间显著的刚度差异导致的平面内位移呈锯齿状，因此准确预测复合材料层合板的变形具有挑战性。为了解决这些限制，本研究提出了一种新的角铺层高阶变形理论，结合计算效率高的精细有限元，来研究角铺层混杂复合材料层合板的弯曲响应。AHDT是通过全局-局部概念建立的，利用勒让德多项式来捕捉锯齿效应和层间应力连续性现象。相应的三角形单元使用C⁰和C1连续性实现该模型，确保高精度和计算效率。通过与三维弹性解和三维有限元模型的比较，严格评估了模型的精度和计算效率。此外，通过碳/玻璃纤维混合复合材料的三点弯曲测试提供了验证，并使用数字图像相关测量了全场变形。理论和实验结果均证实，叠层顺序是影响抗弯刚度的主要因素，其影响大于材料类型和铺层方向。仅改变堆叠顺序就可以减少高达70.77%的挠度。因此，本研究最终表明，合理设计叠层顺序、材料组合和层向是优化复合材料层合板抗弯性能的关键。

{"title":"Experimental and numerical analysis of angle-ply laminated plates based on a refined higher-order theory","authors":"Jie Mei , Shengbo Ling , Yunhao Huang","doi":"10.1016/j.compstruc.2025.108078","DOIUrl":"10.1016/j.compstruc.2025.108078","url":null,"abstract":"<div><div>Accurately predicting the deformation of composite laminates is challenging due to the zig-zag effect in in-plane displacements, which arises from significant stiffness differences between layers. To address these limitations, this study proposes a novel Angle-ply Higher-order Deformation Theory coupled with a computationally efficient refined finite element, to investigate the bending response of angle-ply hybrid composite laminates. The AHDT is established via a global–local concept, utilizing Legendre polynomials to capture the zig-zag effect and interlaminar stress continuity phenomenon. The corresponding triangular element implements this model using C⁰ and C<sup>1</sup> continuities, ensuring high accuracy and computational efficiency. The accuracy and computational efficiency of the model are rigorously assessed through comparisons with three-dimensional elastic solutions and a three-dimensional finite element model. Furthermore, validation is provided via three-point bending tests on carbon/glass fiber hybrid composites, with full-field deformation measured using Digital Image Correlation. Both theoretical and experimental results confirm that the stacking sequence is the predominant factor affecting flexural stiffness, outweighing both material type and ply orientation. A change in the stacking sequence alone can reduce deflection by up to 70.77%. Consequently, this study conclusively shows that a rational design of the stacking sequence, material combination, and ply orientation is critical for optimizing the flexural performance of composite laminates.</div></div>","PeriodicalId":50626,"journal":{"name":"Computers & Structures","volume":"321 ","pages":"Article 108078"},"PeriodicalIF":4.8,"publicationDate":"2026-01-15","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145823166","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}

引用次数: 0

Enhanced load-shape integrated measurement of plate structures based on matched section analysis 基于匹配截面分析的板式结构增强荷载-形状综合测量

IF 4.8 2区工程技术 Q1 COMPUTER SCIENCE, INTERDISCIPLINARY APPLICATIONS

Computers & Structures

Pub Date : 2026-01-15 Epub Date: 2025-12-24 DOI: 10.1016/j.compstruc.2025.108083

Yanhao Guo , Yuhui Hu , Yan Zeng , Feifei Zhao , Hong Bao

Load and deformation measurement of plate structures are fundamental for structural design and health monitoring. However, existing methods are limited to independently reconstructing either deformation or load, without being capable of integrating sensing both the deformation and load in a unified sensing formulation. To address this gap, a matched-section discretization method is proposed, transforming the plate into two orthogonal sets of equivalent beam elements. The structural displacement field is formulated in terms of cross-sectional deformation parameters, internal force/moment parameters, and uniformly distributed loads. Mapping equations between discrete strain measurements and these parameters are then established. By inputting real-time discrete strain data, the least-squares inverse method is employed to simultaneously reconstruct the deformation and load distributions. The proposed method is validated through simulations and experiments on plate structures with varying thicknesses, boundary conditions, and load cases. In the simulation results, the maximum deformation sensing error is limited to 3.1%, and the load sensing error remains below 8.1% with one sensing element. These findings demonstrate the high accuracy and practical potential of the proposed approach in structural health monitoring and intelligent structural design.

板结构的荷载和变形测量是结构设计和健康监测的基础。然而，现有的方法仅限于单独重建变形或载荷，无法将变形和载荷集成到统一的传感公式中。为了解决这一缺陷，提出了一种匹配截面离散化方法，将板转化为两组正交的等效梁单元。结构位移场由截面变形参数、内力/弯矩参数和均布荷载组成。然后建立离散应变测量值与这些参数之间的映射方程。通过输入实时离散应变数据，采用最小二乘反算法同时重构变形和载荷分布。通过对不同厚度、边界条件和载荷情况下的板结构进行仿真和实验，验证了该方法的有效性。在仿真结果中，变形传感的最大误差控制在3.1%以内，负载传感的最大误差控制在8.1%以内。这些结果证明了该方法在结构健康监测和智能结构设计中的高精度和实用潜力。

{"title":"Enhanced load-shape integrated measurement of plate structures based on matched section analysis","authors":"Yanhao Guo , Yuhui Hu , Yan Zeng , Feifei Zhao , Hong Bao","doi":"10.1016/j.compstruc.2025.108083","DOIUrl":"10.1016/j.compstruc.2025.108083","url":null,"abstract":"<div><div>Load and deformation measurement of plate structures are fundamental for structural design and health monitoring. However, existing methods are limited to independently reconstructing either deformation or load, without being capable of integrating sensing both the deformation and load in a unified sensing formulation. To address this gap, a matched-section discretization method is proposed, transforming the plate into two orthogonal sets of equivalent beam elements. The structural displacement field is formulated in terms of cross-sectional deformation parameters, internal force/moment parameters, and uniformly distributed loads. Mapping equations between discrete strain measurements and these parameters are then established. By inputting real-time discrete strain data, the least-squares inverse method is employed to simultaneously reconstruct the deformation and load distributions. The proposed method is validated through simulations and experiments on plate structures with varying thicknesses, boundary conditions, and load cases. In the simulation results, the maximum deformation sensing error is limited to 3.1%, and the load sensing error remains below 8.1% with one sensing element. These findings demonstrate the high accuracy and practical potential of the proposed approach in structural health monitoring and intelligent structural design.</div></div>","PeriodicalId":50626,"journal":{"name":"Computers & Structures","volume":"321 ","pages":"Article 108083"},"PeriodicalIF":4.8,"publicationDate":"2026-01-15","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145823164","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}

引用次数: 0

Development of isogeometric-finite element unified numerical methodology for nonlinear planar instability analysis of dual-phase SiO2-graphene nanoplatetlets reinforced shallow curved microbeams 双相二氧化硅-石墨烯纳米片增强浅弯曲微梁非线性平面失稳分析等几何-有限元统一数值方法的发展

IF 4.8 2区工程技术 Q1 COMPUTER SCIENCE, INTERDISCIPLINARY APPLICATIONS

Computers & Structures

Pub Date : 2026-01-15 Epub Date: 2025-12-24 DOI: 10.1016/j.compstruc.2025.108082

Saeid Sahmani , Reza Ansari , Mohammad Kazem Hassanzadeh-Aghdam , Mona Zareichian

The main intention of the proposed multiscale framework is to employ an isogeometric analysis formulation for the size-dependent nonlinear planar instability analysis of dual-phase nanocomposite shallow curved beams at the microscale. A finite element–based micromechanics model is developed at the representative volume element level to capture the material properties. The homogenized properties obtained from the representative volume element-level finite element analysis are directly incorporated into the isogeometric model. This coupling enables accurate surveying of the small scale-dependent nonlinear in-plane stability characteristics of uniformly laterally loaded dual-phase inhomogeneous shallow curved microbeams reinforced with SiO₂ nanoparticles and graphene nanoplatelets, while embracing distinct strain gradient tensors. In this regard, cuboid-shaped representative volume elements are employed. This enables consideration of the interphase between the dual-phase nanofillers and the polymer, as well as the critical role of nanofiller agglomeration, in order to create an accurate multiscale correlation. Additionally, non-uniform rational B-splines are utilized in the relevant discretization process. This process involves distinct microstructural-dependent strain gradient tensors. The numerical results reveal that increasing the SiO₂ nanoparticle volume fraction significantly enhances both the upper and lower limit loads by nearly 69.5%. This increase does not markedly affect the axial resultant load or the lateral deflection. Conversely, increasing the SiO₂ nanoparticle diameter at a fixed volume fraction notably decreases the load-bearing capacity by about 49.5%. Similarly, a rise in graphene nanoplatelet thickness leads to an approximately 61.1% reduction in the stability limits. The inclusion of the interphase region between the nanofillers and the matrix improves the upper and lower limit loads by around 17.4%, demonstrating its reinforcing influence. Furthermore, aligning nanofillers along the beam’s longitudinal direction increases the limit loads by roughly 48.1% compared to the random dispersion case. In contrast, agglomeration has the opposite effect, reducing the load-carrying capacity by about 12.4%.

所提出的多尺度框架的主要目的是采用一种等几何分析公式来分析微尺度下双相纳米复合材料浅弯曲梁的尺寸相关非线性平面不稳定性。在具有代表性的体元水平上建立了基于有限元的细观力学模型来捕捉材料的特性。从具有代表性的体积单元级有限元分析中获得的均质性直接纳入等几何模型。这种耦合可以精确测量均匀横向加载的双相非均匀浅弯曲微梁的小尺度非线性面内稳定性特性，这些微梁由SiO2纳米颗粒和石墨烯纳米片增强，同时包含不同的应变梯度张量。在这方面，采用长方体形状的代表性体积元素。这可以考虑双相纳米填料和聚合物之间的界面，以及纳米填料团聚的关键作用，以便创建精确的多尺度关联。另外，在相关的离散化过程中采用了非均匀有理b样条。这一过程涉及不同的微观结构相关的应变梯度张量。数值结果表明，增加SiO2纳米颗粒体积分数，上、下限载荷均显著提高近69.5%。这种增加并不显著影响轴向合成载荷或横向挠度。相反，在一定体积分数下，增加SiO2纳米颗粒直径，其承载能力显著降低约49.5%。同样，石墨烯纳米板厚度的增加导致稳定性极限降低约61.1%。纳米填料与基体间相区域的加入使材料的上下限载荷提高了17.4%左右，表明了纳米填料的强化作用。此外，与随机分散的情况相比，沿梁的纵向排列纳米填料增加了大约48.1%的极限载荷。而结块则相反，其承载能力降低约12.4%。

{"title":"Development of isogeometric-finite element unified numerical methodology for nonlinear planar instability analysis of dual-phase SiO2-graphene nanoplatetlets reinforced shallow curved microbeams","authors":"Saeid Sahmani , Reza Ansari , Mohammad Kazem Hassanzadeh-Aghdam , Mona Zareichian","doi":"10.1016/j.compstruc.2025.108082","DOIUrl":"10.1016/j.compstruc.2025.108082","url":null,"abstract":"<div><div>The main intention of the proposed multiscale framework is to employ an isogeometric analysis formulation for the size-dependent nonlinear planar instability analysis of dual-phase nanocomposite shallow curved beams at the microscale. A finite element–based micromechanics model is developed at the representative volume element level to capture the material properties. The homogenized properties obtained from the representative volume element-level finite element analysis are directly incorporated into the isogeometric model. This coupling enables accurate surveying of the small scale-dependent nonlinear in-plane stability characteristics of uniformly laterally loaded dual-phase inhomogeneous shallow curved microbeams reinforced with SiO<sub>2</sub> nanoparticles and graphene nanoplatelets, while embracing distinct strain gradient tensors. In this regard, cuboid-shaped representative volume elements are employed. This enables consideration of the interphase between the dual-phase nanofillers and the polymer, as well as the critical role of nanofiller agglomeration, in order to create an accurate multiscale correlation. Additionally, non-uniform rational B-splines are utilized in the relevant discretization process. This process involves distinct microstructural-dependent strain gradient tensors. The numerical results reveal that increasing the SiO<sub>2</sub> nanoparticle volume fraction significantly enhances both the upper and lower limit loads by nearly 69.5%. This increase does not markedly affect the axial resultant load or the lateral deflection. Conversely, increasing the SiO<sub>2</sub> nanoparticle diameter at a fixed volume fraction notably decreases the load-bearing capacity by about 49.5%. Similarly, a rise in graphene nanoplatelet thickness leads to an approximately 61.1% reduction in the stability limits. The inclusion of the interphase region between the nanofillers and the matrix improves the upper and lower limit loads by around 17.4%, demonstrating its reinforcing influence. Furthermore, aligning nanofillers along the beam’s longitudinal direction increases the limit loads by roughly 48.1% compared to the random dispersion case. In contrast, agglomeration has the opposite effect, reducing the load-carrying capacity by about 12.4%.</div></div>","PeriodicalId":50626,"journal":{"name":"Computers & Structures","volume":"321 ","pages":"Article 108082"},"PeriodicalIF":4.8,"publicationDate":"2026-01-15","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145823163","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}

引用次数: 0

A finite strain elastic-corrector rate formulation for void growth models 孔洞生长模型的有限应变弹性修正率公式

IF 4.8 2区工程技术 Q1 COMPUTER SCIENCE, INTERDISCIPLINARY APPLICATIONS

Computers & Structures

Pub Date : 2026-01-15 Epub Date: 2025-12-25 DOI: 10.1016/j.compstruc.2025.108072

Meijuan Zhang , Guadalupe Vadillo , Francisco J. Montáns

Recently, a novel framework for large strain elastoplasticity has been introduced. This approach preserves the multiplicative decomposition of the deformation gradient, avoiding Green-like elastoplastic decompositions of total strains, but resulting in additive predictor-corrector algorithmic implementations as in small strains, and using any form of hyperelastic stored energy. Using this approach, large strain complexities reduce to non-iterative kinematical pre- and post-processors. In this work, we extend the elastic-corrector rates formulation for dealing with non-isochoric void-growth plasticity. In contrast to other large strain formulations, in the present work the relation of the void growth with the plastic flow of the matrix is not assumed, but obtained from the multiplicative decomposition, preserving exactly the kinematics of the volumetric parts. We use the Gurson-Tvergaard-Needleman (GTN) function to demonstrate the use of the framework as well as the implementation of a fully implicit Newton algorithm. The novel GTN formulation is hyperelastic, does not need the exponential mapping, and keeps an additive algorithmic update as in the infinitesimal formulation. Other non-isochoric functions can also be employed with the same kinematics.

近年来，提出了一种新的大应变弹塑性框架。这种方法保留了变形梯度的乘法分解，避免了总应变的格林样弹塑性分解，但导致了小应变的加性预测校正算法的实现，并使用了任何形式的超弹性存储能量。使用这种方法，大的应变复杂性减少到非迭代的运动学预处理和后处理。在这项工作中，我们扩展了弹性校正率公式来处理非等时空隙生长塑性。与其他大应变公式不同的是，在本工作中，孔洞生长与基体塑性流动的关系不是假设的，而是从乘法分解中得到的，准确地保留了体积部分的运动学。我们使用Gurson-Tvergaard-Needleman （GTN）函数来演示该框架的使用以及完全隐式牛顿算法的实现。新的GTN公式是超弹性的，不需要指数映射，并且像无穷小公式一样保持加性算法更新。其他的非等时函数也可以采用相同的运动学。

{"title":"A finite strain elastic-corrector rate formulation for void growth models","authors":"Meijuan Zhang , Guadalupe Vadillo , Francisco J. Montáns","doi":"10.1016/j.compstruc.2025.108072","DOIUrl":"10.1016/j.compstruc.2025.108072","url":null,"abstract":"<div><div>Recently, a novel framework for large strain elastoplasticity has been introduced. This approach preserves the multiplicative decomposition of the deformation gradient, avoiding Green-like elastoplastic decompositions of total strains, but resulting in additive predictor-corrector algorithmic implementations as in small strains, and using any form of hyperelastic stored energy. Using this approach, large strain complexities reduce to non-iterative kinematical pre- and post-processors. In this work, we extend the elastic-corrector rates formulation for dealing with non-isochoric void-growth plasticity. In contrast to other large strain formulations, in the present work the relation of the void growth with the plastic flow of the matrix is not assumed, but obtained from the multiplicative decomposition, preserving exactly the kinematics of the volumetric parts. We use the Gurson-Tvergaard-Needleman (GTN) function to demonstrate the use of the framework as well as the implementation of a fully implicit Newton algorithm. The novel GTN formulation is hyperelastic, does not need the exponential mapping, and keeps an additive algorithmic update as in the infinitesimal formulation. Other non-isochoric functions can also be employed with the same kinematics.</div></div>","PeriodicalId":50626,"journal":{"name":"Computers & Structures","volume":"321 ","pages":"Article 108072"},"PeriodicalIF":4.8,"publicationDate":"2026-01-15","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145823155","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}

引用次数: 0

An energy-limited gradient damage approach for 3D thermal fracture analysis 三维热断裂分析的能量限制梯度损伤方法

IF 4.8 2区工程技术 Q1 COMPUTER SCIENCE, INTERDISCIPLINARY APPLICATIONS

Computers & Structures

Pub Date : 2026-01-15 Epub Date: 2025-12-28 DOI: 10.1016/j.compstruc.2025.108084

Cancan Su , Dechun Lu , Xiaoying Zhuang , Timon Rabczuk , Xin Zhou , Qin He , Xiuli Du

This study presents a thermo-mechanically coupled nonlocal damage model for brittle fracture that is practically implementable and extensible to multi-physics applications. The proposed damage formulation preserves thermodynamic consistency within a variational framework and employs an implicit gradient scheme for efficient regularization. A local damage variable is defined to represent different fracture modes, while a Helmholtz-type PDE is introduced to mitigate mesh dependence. The model is first verified through mechanical benchmark problems, which confirm reduced mesh sensitivity. Its predictive capability is then demonstrated in thermo-mechanical simulations involving diverse geometries, initial conditions, material properties, and spatial dimensions (2D and 3D). The results consistently indicate that the proposed approach delivers robust and accurate predictions of thermal fracture processes while maintaining high computational efficiency.

本研究提出了一种可实际实现并可扩展到多物理场应用的脆性断裂热-机械耦合非局部损伤模型。提出的损伤公式在变分框架内保持热力学一致性，并采用隐式梯度格式进行有效正则化。定义了一个局部损伤变量来表示不同的断裂模式，同时引入了Helmholtz-type PDE来减轻网格依赖性。首先通过力学基准问题对模型进行验证，证实了网格灵敏度的降低。然后在热机械模拟中展示了其预测能力，包括不同的几何形状、初始条件、材料特性和空间尺寸（2D和3D）。结果一致表明，该方法在保持较高的计算效率的同时，提供了对热破裂过程的鲁棒性和准确性预测。

引用次数: 0

Development of a hybrid dynamic stiffness method adapted to Framed structures 适用于框架结构的混合动力刚度方法的发展

IF 4.8 2区工程技术 Q1 COMPUTER SCIENCE, INTERDISCIPLINARY APPLICATIONS

Computers & Structures

Pub Date : 2026-01-15 Epub Date: 2025-12-30 DOI: 10.1016/j.compstruc.2025.108088

A. Brisson , E. Perrey-Debain , J.D. Chazot , C. Guigou-Carter , C. Fillol , P. Jean

Mitigating noise and vibration nuisances in engineering design requires efficient and rapid simulation tools during the design and pre-project phases. In civil engineering, post-and-beam structures which consist of skeletal frameworks composed of slabs, beams, and posts, semi-analytical methods such as the popular Dynamic Stiffness Method can be employed. The method avoids heavy computational resources but it relies on certain restrictive modeling simplifications. In contrast, the 3D Finite Element Method allows for taking into account all geometrical details of the structure, resulting in precise and reliable outcomes. In this paper, we propose a hybrid approach that combines both methods to enhance the accuracy of the Dynamic Stiffness Method while maintaining fast computation times. Specifically, the connections between posts and beams are modeled precisely with the FEM while the rest of the building structure is modeled with the Dynamic Stiffness Method. An additional Craig-Bampton reduction step is also applied to the finite element domain to further decrease the problem size. This article details the development of the hybrid dynamic stiffness method, demonstrating its effectiveness in terms of accuracy and model size reduction through comparisons with a full finite element model. Finally, a real application case is presented to test the method on a full-scale building structure.

减轻工程设计中的噪声和振动干扰需要在设计和项目前期阶段高效快速的仿真工具。在土木工程中，由板、梁和桩组成的骨架框架组成的梁柱结构，可以采用半解析方法，如流行的动刚度法。该方法避免了大量的计算资源，但依赖于一定的限制性建模简化。相比之下，三维有限元法允许考虑结构的所有几何细节，从而产生精确可靠的结果。在本文中，我们提出了一种将两种方法结合起来的混合方法，以提高动刚度法的精度，同时保持快速的计算时间。其中，梁柱之间的连接采用有限元法进行精确建模，其余部分采用动刚度法进行建模。一个额外的Craig-Bampton约简步骤也被应用到有限元域，以进一步减小问题的大小。本文详细介绍了混合动力刚度方法的发展，通过与完整有限元模型的比较，证明了其在精度和模型尺寸减小方面的有效性。最后，给出了一个实际应用案例，对该方法在全尺寸建筑结构上进行了验证。

{"title":"Development of a hybrid dynamic stiffness method adapted to Framed structures","authors":"A. Brisson , E. Perrey-Debain , J.D. Chazot , C. Guigou-Carter , C. Fillol , P. Jean","doi":"10.1016/j.compstruc.2025.108088","DOIUrl":"10.1016/j.compstruc.2025.108088","url":null,"abstract":"<div><div>Mitigating noise and vibration nuisances in engineering design requires efficient and rapid simulation tools during the design and pre-project phases. In civil engineering, post-and-beam structures which consist of skeletal frameworks composed of slabs, beams, and posts, semi-analytical methods such as the popular Dynamic Stiffness Method can be employed. The method avoids heavy computational resources but it relies on certain restrictive modeling simplifications. In contrast, the 3D Finite Element Method allows for taking into account all geometrical details of the structure, resulting in precise and reliable outcomes. In this paper, we propose a hybrid approach that combines both methods to enhance the accuracy of the Dynamic Stiffness Method while maintaining fast computation times. Specifically, the connections between posts and beams are modeled precisely with the FEM while the rest of the building structure is modeled with the Dynamic Stiffness Method. An additional Craig-Bampton reduction step is also applied to the finite element domain to further decrease the problem size. This article details the development of the hybrid dynamic stiffness method, demonstrating its effectiveness in terms of accuracy and model size reduction through comparisons with a full finite element model. Finally, a real application case is presented to test the method on a full-scale building structure.</div></div>","PeriodicalId":50626,"journal":{"name":"Computers & Structures","volume":"321 ","pages":"Article 108088"},"PeriodicalIF":4.8,"publicationDate":"2026-01-15","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145884461","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}

引用次数: 0

Model calibration and uncertainty quantification of cellular automata-based pitting corrosion model 基于元胞自动机的点蚀模型标定与不确定度量化

IF 4.8 2区工程技术 Q1 COMPUTER SCIENCE, INTERDISCIPLINARY APPLICATIONS

Computers & Structures

Pub Date : 2026-01-15 Epub Date: 2026-01-20 DOI: 10.1016/j.compstruc.2026.108103

J.Ramesh Babu, Pranav M. Karve, Sankaran Mahadevan

This study presents a Bayesian framework for calibrating cellular automata (CA) models of pitting corrosion. A high-fidelity two-dimensional CA model is used to simulate pit morphology evolution under coupled electrochemical and mass transport processes, incorporating eight uncertain model parameters related to metal dissolution, hydrolysis, diffusion, and sedimentation. Key geometric indicators—maximum pit depth and pit aspect ratio—are extracted from the simulation outputs and reduced to interpretable scalar features using power-law fitting. Correlation analysis is conducted to assess linear relationships between model parameters and power-law coefficients. Gaussian Process Regression (GPR) surrogate models are constructed using Latin Hypercube Sampling (LHS) to emulate the high-fidelity CA model. Global sensitivity analysis (GSA) using Sobol’ indices is performed by utilizing the trained surrogate models; it identifies sedimentation and ionic diffusion as the dominant contributors to output feature variability, with significant nonlinear interactions. The trained surrogates are integrated into a Bayesian inference framework using Metropolis-Hastings Markov Chain Monte Carlo (MH-MCMC) sampling to infer posterior distributions of the uncertain model parameters. The posterior samples are further propagated through the surrogates to quantify the output uncertainty, and the prediction is evaluated using a distance-based probabilistic model validation metric. Two numerical examples related to API-5L X65 steel pipelines are presented, demonstrating that incorporating multiple geometric features—both pit depth and aspect ratio—improves predictive accuracy. The proposed framework supports uncertainty-aware modeling and decision-making for structural health assessment and maintenance planning in corrosion-critical infrastructure.

本研究提出了一个贝叶斯框架来校准点蚀的元胞自动机（CA）模型。采用高保真二维CA模型，结合金属溶解、水解、扩散和沉降等8个不确定模型参数，模拟电化学和质量输运耦合过程下的矿坑形态演化。从模拟输出中提取关键几何指标-最大坑深和坑长比，并使用幂律拟合将其简化为可解释的标量特征。通过相关分析来评估模型参数与幂律系数之间的线性关系。利用拉丁超立方采样（LHS）构造高斯过程回归（GPR）代理模型来模拟高保真CA模型。使用Sobol指数进行全局敏感性分析（GSA），利用训练好的代理模型；它确定沉降和离子扩散是输出特征变异性的主要贡献者，具有显著的非线性相互作用。利用Metropolis-Hastings马尔可夫链蒙特卡罗（MH-MCMC）抽样方法，将训练好的替代模型整合到贝叶斯推理框架中，推断不确定模型参数的后验分布。后验样本通过代理进一步传播以量化输出不确定性，并使用基于距离的概率模型验证度量来评估预测。给出了两个与API-5L X65钢管道相关的数值实例，表明结合多种几何特征（坑深和纵横比）可以提高预测精度。提出的框架支持不确定性感知建模和决策，用于关键腐蚀基础设施的结构健康评估和维护规划。

{"title":"Model calibration and uncertainty quantification of cellular automata-based pitting corrosion model","authors":"J.Ramesh Babu, Pranav M. Karve, Sankaran Mahadevan","doi":"10.1016/j.compstruc.2026.108103","DOIUrl":"10.1016/j.compstruc.2026.108103","url":null,"abstract":"<div><div>This study presents a Bayesian framework for calibrating cellular automata (CA) models of pitting corrosion. A high-fidelity two-dimensional CA model is used to simulate pit morphology evolution under coupled electrochemical and mass transport processes, incorporating eight uncertain model parameters related to metal dissolution, hydrolysis, diffusion, and sedimentation. Key geometric indicators—maximum pit depth and pit aspect ratio—are extracted from the simulation outputs and reduced to interpretable scalar features using power-law fitting. Correlation analysis is conducted to assess linear relationships between model parameters and power-law coefficients. Gaussian Process Regression (GPR) surrogate models are constructed using Latin Hypercube Sampling (LHS) to emulate the high-fidelity CA model. Global sensitivity analysis (GSA) using Sobol’ indices is performed by utilizing the trained surrogate models; it identifies sedimentation and ionic diffusion as the dominant contributors to output feature variability, with significant nonlinear interactions. The trained surrogates are integrated into a Bayesian inference framework using Metropolis-Hastings Markov Chain Monte Carlo (MH-MCMC) sampling to infer posterior distributions of the uncertain model parameters. The posterior samples are further propagated through the surrogates to quantify the output uncertainty, and the prediction is evaluated using a distance-based probabilistic model validation metric. Two numerical examples related to API-5L X65 steel pipelines are presented, demonstrating that incorporating multiple geometric features—both pit depth and aspect ratio—improves predictive accuracy. The proposed framework supports uncertainty-aware modeling and decision-making for structural health assessment and maintenance planning in corrosion-critical infrastructure.</div></div>","PeriodicalId":50626,"journal":{"name":"Computers & Structures","volume":"321 ","pages":"Article 108103"},"PeriodicalIF":4.8,"publicationDate":"2026-01-15","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"146014886","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}

引用次数: 0

High-order time-spectral BEM for efficient elastodynamic analysis 高效弹性动力学分析的高阶时谱边界元

IF 4.8 2区工程技术 Q1 COMPUTER SCIENCE, INTERDISCIPLINARY APPLICATIONS

Computers & Structures

Pub Date : 2026-01-15 Epub Date: 2025-12-11 DOI: 10.1016/j.compstruc.2025.108076

Yan Gu , Wenzhen Qu , Chuanzeng Zhang , Vladimir Babeshko , Yuri V. Petrov

This study presents a novel boundary element method (BEM) framework for the accurate and efficient numerical solution of elastodynamic problems. By reformulating the time-derivative terms as equivalent body forces, the method enables the use of static fundamental solutions for dynamic analysis, thereby eliminating the need for frequency-domain transformations or the construction of complex time-dependent Green’s functions. In the temporal domain, instead of directly approximating the time-differentiation operators, a more stable time-spectral integration technique based on orthogonal polynomial expansions is introduced. This scheme can, in principle, achieve arbitrary-order of accuracy in the temporal space and eliminate the strict time-step limitations inherent in conventional finite-difference-based schemes. Moreover, the resulting coefficient matrix is time-independent and therefore needs to be computed only once for the entire time-marching process. To evaluate domain integrals, discontinuous triangular elements are employed for spatial discretization, and a scaled coordinate transformation (SCT) technique is employed to address singularities arising from source-field point coincidences. Preliminary numerical experiments in elastodynamic analysis demonstrate that the proposed framework is robust and flexible for long-time dynamic simulations, particularly in problems involving rapid transients or complex geometries.

本文提出了一种新的边界元法框架，用于精确、高效地求解弹性动力学问题。通过将时间导数项重新表述为等效体力，该方法可以使用静态基本解进行动态分析，从而消除了频域变换或构建复杂的时间相关格林函数的需要。在时域中，引入了一种基于正交多项式展开的更稳定的时谱积分技术，而不是直接逼近微分算子。该方案原则上可以在时间空间上实现任意阶的精度，消除了传统有限差分方案固有的严格的时间步长限制。此外，得到的系数矩阵是时间无关的，因此在整个时间推进过程中只需要计算一次。在计算域积分时，采用不连续三角形单元进行空间离散，并采用缩放坐标变换（SCT）技术处理源场点重合引起的奇异点。弹性动力学分析的初步数值实验表明，所提出的框架对于长时间的动态模拟具有鲁棒性和灵活性，特别是在涉及快速瞬态或复杂几何形状的问题中。

{"title":"High-order time-spectral BEM for efficient elastodynamic analysis","authors":"Yan Gu , Wenzhen Qu , Chuanzeng Zhang , Vladimir Babeshko , Yuri V. Petrov","doi":"10.1016/j.compstruc.2025.108076","DOIUrl":"10.1016/j.compstruc.2025.108076","url":null,"abstract":"<div><div>This study presents a novel boundary element method (BEM) framework for the accurate and efficient numerical solution of elastodynamic problems. By reformulating the time-derivative terms as equivalent body forces, the method enables the use of static fundamental solutions for dynamic analysis, thereby eliminating the need for frequency-domain transformations or the construction of complex time-dependent Green’s functions. In the temporal domain, instead of directly approximating the time-differentiation operators, a more stable time-spectral integration technique based on orthogonal polynomial expansions is introduced. This scheme can, in principle, achieve arbitrary-order of accuracy in the temporal space and eliminate the strict time-step limitations inherent in conventional finite-difference-based schemes. Moreover, the resulting coefficient matrix is time-independent and therefore needs to be computed only once for the entire time-marching process. To evaluate domain integrals, discontinuous triangular elements are employed for spatial discretization, and a scaled coordinate transformation (SCT) technique is employed to address singularities arising from source-field point coincidences. Preliminary numerical experiments in elastodynamic analysis demonstrate that the proposed framework is robust and flexible for long-time dynamic simulations, particularly in problems involving rapid transients or complex geometries.</div></div>","PeriodicalId":50626,"journal":{"name":"Computers & Structures","volume":"321 ","pages":"Article 108076"},"PeriodicalIF":4.8,"publicationDate":"2026-01-15","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145731206","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}

引用次数: 0

Multi-patch isogeometric shape optimization of complex free-form surfaces with buckling constraints 具有屈曲约束的复杂自由曲面的多斑块等几何形状优化

IF 4.8 2区工程技术 Q1 COMPUTER SCIENCE, INTERDISCIPLINARY APPLICATIONS

Computers & Structures

Pub Date : 2026-01-15 Epub Date: 2026-01-02 DOI: 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.

结构形状优化对于找到美观且力学性能合理的结构形状设计至关重要。等几何分析以其统一设计模型和分析模型的优点，为此类优化提供了一种很有前途的方法。本文提出了一种在等几何分析框架下的自由曲面的综合形状优化方法，解决了柔度和屈曲问题。导出了Kirchhoff-Love壳的解析解，实现了基于梯度的高效优化。对于包含多个非均匀有理b样条（NURBS）斑块的复杂自由曲面，采用无梯度优化策略保证鲁棒性。通过Nitsche的方法，实现了跨多补丁接口的连续性约束。通过几个数值算例验证了该方法的有效性，表明该方法在结构柔度和屈曲稳定性方面同时得到了改善。

引用次数: 0