Applied Mathematical Modelling最新文献_第9页

A new method of specifying parameter bounds for optimization of reduced-order models and application in design of controllers 一种确定降阶模型优化参数界的新方法及其在控制器设计中的应用

IF 4.4 2区工程技术 Q1 ENGINEERING, MULTIDISCIPLINARY

Applied Mathematical Modelling

Pub Date : 2026-01-10 DOI: 10.1016/j.apm.2026.116751

Anuj Goel , Amit Kumar Manocha , Parveen Bajaj , George Uwadiegwu Alaneme

In the present work, a new framework is proposed to design controllers using model order reduction techniques for linear time invariant complex engineering systems. The proposed model order reduction methodology employs optimization-based techniques namely ant lion optimization and moth flame optimization for which boundary conditions are systematically procured from an interim model derived using balancing free square-root algorithm. An area control coefficient is introduced to adjust the exploration range of the optimization process around the coefficients of the interim reduced-order model. The numerator as well as denominator coefficients of the desired reduced-order models are optimized to retain the performance characteristics of the original high-order systems. The effectiveness of the proposed approach is assessed based on different error metrics and unit step response plots. To validate the performance, seven benchmark systems of different pole configurations have been considered from the literature. It has been found that proposed approach provides reduced-systems with significant improvement in error and transient performance when compared to the literature work. The suggested model order reduction approach is further extended to design proportional-integral-derivative controller and fractional-order proportional-integral-derivative controller for an 84th-order benchmark system and a mechanical ventilator system respectively. The results demonstrate that the proposed model order reduction-based controller design approach achieves high-performance control with lesser steady-state error, improved time-domain specifications and robust disturbance rejection capability.

本文提出了一种利用模型降阶技术设计线性时不变复杂工程系统控制器的新框架。所提出的模型降阶方法采用基于优化的技术，即蚁狮优化和蛾焰优化，其边界条件是系统地从使用平衡自由平方根算法导出的临时模型中获得的。在中间降阶模型系数周围引入面积控制系数来调整优化过程的勘探范围。对期望的降阶模型的分子和分母系数进行了优化，以保持原高阶系统的性能特征。基于不同的误差度量和单位阶跃响应图对所提方法的有效性进行了评估。为了验证性能，从文献中考虑了七个不同极点配置的基准系统。研究发现，与文献相比，本文提出的方法在误差和瞬态性能方面提供了显著改善的简化系统。将模型降阶方法进一步推广到84阶基准系统和机械通风机系统的比例-积分-导数控制器和分数阶比例-积分-导数控制器的设计中。结果表明，所提出的基于模型阶数约简的控制器设计方法具有较小的稳态误差、较好的时域参数和鲁棒抗扰能力，实现了高性能的控制。

{"title":"A new method of specifying parameter bounds for optimization of reduced-order models and application in design of controllers","authors":"Anuj Goel , Amit Kumar Manocha , Parveen Bajaj , George Uwadiegwu Alaneme","doi":"10.1016/j.apm.2026.116751","DOIUrl":"10.1016/j.apm.2026.116751","url":null,"abstract":"<div><div>In the present work, a new framework is proposed to design controllers using model order reduction techniques for linear time invariant complex engineering systems. The proposed model order reduction methodology employs optimization-based techniques namely ant lion optimization and moth flame optimization for which boundary conditions are systematically procured from an interim model derived using balancing free square-root algorithm. An area control coefficient is introduced to adjust the exploration range of the optimization process around the coefficients of the interim reduced-order model. The numerator as well as denominator coefficients of the desired reduced-order models are optimized to retain the performance characteristics of the original high-order systems. The effectiveness of the proposed approach is assessed based on different error metrics and unit step response plots. To validate the performance, seven benchmark systems of different pole configurations have been considered from the literature. It has been found that proposed approach provides reduced-systems with significant improvement in error and transient performance when compared to the literature work. The suggested model order reduction approach is further extended to design proportional-integral-derivative controller and fractional-order proportional-integral-derivative controller for an 84th-order benchmark system and a mechanical ventilator system respectively. The results demonstrate that the proposed model order reduction-based controller design approach achieves high-performance control with lesser steady-state error, improved time-domain specifications and robust disturbance rejection capability.</div></div>","PeriodicalId":50980,"journal":{"name":"Applied Mathematical Modelling","volume":"155 ","pages":"Article 116751"},"PeriodicalIF":4.4,"publicationDate":"2026-01-10","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145956844","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

Efficient simulation method of non-stationary non-Gaussian stochastic ground motions based on adaptive interpolation strategy 基于自适应插值策略的非平稳非高斯随机地震动高效模拟方法

IF 4.4 2区工程技术 Q1 ENGINEERING, MULTIDISCIPLINARY

Applied Mathematical Modelling

Pub Date : 2026-01-10 DOI: 10.1016/j.apm.2026.116746

Xiangqian Sheng , Kuahai Yu , Wenliang Fan , Baiwan Su

The translation process-based spectral representation method is widely used to simulate the non-stationary non-Gaussian stochastic ground motions. However, the computation of the underlying evolutionary power spectral density matrix and its decomposition require substantial computational effort at discrete time-frequency points. To address this problem, this paper proposes an adaptive interpolation strategy for selecting the time-frequency interpolation points to improve the simulation efficiency. Firstly, the correlation function equations between non-Gaussian stochastic processes and the underlying Gaussian stochastic processes are constructed using Mehler's formula. A fast calculation method for the evolutionary power spectral density of the underlying Gaussian processes is introduced based on the interpolation technique. Secondly, the discrete time-frequency interpolation points are determined based on the amplitude information of the evolutionary power spectral density of the underlying Gaussian stochastic processes. The evolutionary power spectral density matrix is decomposed at these time-frequency interpolation points. The decomposed spectrum is then expressed as a sum of products of various time and frequency components. Spline interpolation is applied to these components at the discrete time-frequency points to approximate the matrix decomposition required by the spectral representation method, improving the efficiency of the decomposition. Additionally, the Fast Fourier Transform further accelerates simulation efficiency. Finally, the accuracy and efficiency of the proposed method for simulating non-stationary non-Gaussian stochastic ground motions are verified by considering the real ground motion record, stochastic vector processes, different probability distribution types, different power spectrum density types, and the number of variates.

基于平移过程的谱表示方法被广泛用于模拟非平稳非高斯随机地震动。然而，底层进化功率谱密度矩阵的计算及其分解需要在离散时频点进行大量的计算。针对这一问题，本文提出了一种选择时频插值点的自适应插值策略，以提高仿真效率。首先，利用梅勒公式建立了非高斯随机过程与底层高斯随机过程的相关函数方程；介绍了一种基于插值技术的高斯过程演化功率谱密度的快速计算方法。其次，根据高斯随机过程演化功率谱密度的幅值信息确定离散时频插值点；在这些时频插值点处对演化功率谱密度矩阵进行分解。然后将分解后的频谱表示为各种时间和频率分量的乘积和。在离散时频点对这些分量进行样条插值，逼近谱表示法所需的矩阵分解，提高了分解效率。此外，快速傅里叶变换进一步提高了仿真效率。最后，通过考虑实际地震动记录、随机矢量过程、不同概率分布类型、不同功率谱密度类型和变量数量，验证了所提方法模拟非平稳非高斯随机地震动的准确性和效率。

{"title":"Efficient simulation method of non-stationary non-Gaussian stochastic ground motions based on adaptive interpolation strategy","authors":"Xiangqian Sheng , Kuahai Yu , Wenliang Fan , Baiwan Su","doi":"10.1016/j.apm.2026.116746","DOIUrl":"10.1016/j.apm.2026.116746","url":null,"abstract":"<div><div>The translation process-based spectral representation method is widely used to simulate the non-stationary non-Gaussian stochastic ground motions. However, the computation of the underlying evolutionary power spectral density matrix and its decomposition require substantial computational effort at discrete time-frequency points. To address this problem, this paper proposes an adaptive interpolation strategy for selecting the time-frequency interpolation points to improve the simulation efficiency. Firstly, the correlation function equations between non-Gaussian stochastic processes and the underlying Gaussian stochastic processes are constructed using Mehler's formula. A fast calculation method for the evolutionary power spectral density of the underlying Gaussian processes is introduced based on the interpolation technique. Secondly, the discrete time-frequency interpolation points are determined based on the amplitude information of the evolutionary power spectral density of the underlying Gaussian stochastic processes. The evolutionary power spectral density matrix is decomposed at these time-frequency interpolation points. The decomposed spectrum is then expressed as a sum of products of various time and frequency components. Spline interpolation is applied to these components at the discrete time-frequency points to approximate the matrix decomposition required by the spectral representation method, improving the efficiency of the decomposition. Additionally, the Fast Fourier Transform further accelerates simulation efficiency. Finally, the accuracy and efficiency of the proposed method for simulating non-stationary non-Gaussian stochastic ground motions are verified by considering the real ground motion record, stochastic vector processes, different probability distribution types, different power spectrum density types, and the number of variates.</div></div>","PeriodicalId":50980,"journal":{"name":"Applied Mathematical Modelling","volume":"155 ","pages":"Article 116746"},"PeriodicalIF":4.4,"publicationDate":"2026-01-10","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145956846","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 efficient adaptive time-integration method for second-order nonlinear dynamics 二阶非线性动力学的一种有效自适应时间积分方法

IF 4.4 2区工程技术 Q1 ENGINEERING, MULTIDISCIPLINARY

Applied Mathematical Modelling

Pub Date : 2026-01-09 DOI: 10.1016/j.apm.2026.116740

Zhi Duan , Xiaohui Liu , Hongbing Guo , Chuan Wu , Zhongfei Ye , Zhongbin Lu

Achieving an optimal balance among computational efficiency, robustness, and accuracy is a central challenge in simulating second-order nonlinear dynamical systems. While the existing parameterized two-sub-step composite integrator provides rigorous nonlinear stability and controllable dissipation, its fixed-step formulation limits efficiency in simulations with strongly varying dynamics. This paper presents a novel adaptive time integration method that augments the second-order base scheme with an explicit auxiliary stage for efficient error estimation. Its key innovation is a cost-free, direct error estimator, constructed by rigorously deriving the embedding coefficients via order conditions and analytically combining the implicit base stages with an extrapolated explicit stage to derive a local error estimate based on a third-order embedding without additional nonlinear iterations or matrix operations. Combined with a proportional–integral–derivative-like step-size controller, systematic numerical tests show that the proposed method achieves a significantly better computational cost-to-accuracy trade-off than high-order algebraically stable singly diagonally implicit Runge–Kutta methods. The algorithm demonstrates strong robustness in stiff and large-scale nonlinear problems while preserving the unconditional nonlinear stability and controllable dissipation of the base scheme. In summary, the proposed adaptive method offers an efficient, reliable, and self-starting tool for simulating large-scale, long-duration, strongly nonlinear systems.

实现计算效率、鲁棒性和准确性之间的最佳平衡是模拟二阶非线性动力系统的核心挑战。虽然现有的参数化两子步复合积分器提供了严格的非线性稳定性和可控耗散，但其定步式限制了其在强动态变化模拟中的效率。本文提出了一种新的自适应时间积分方法，该方法在二阶基格式的基础上增加了一个显式的辅助阶段，以实现有效的误差估计。它的关键创新是一种无成本的直接误差估计器，该估计器通过严格地根据阶条件推导嵌入系数，并将隐式基本阶段与外推的显式阶段解析地结合起来，从而基于三阶嵌入推导出局部误差估计，而无需额外的非线性迭代或矩阵运算。结合比例-积分-导数类步长控制器，系统数值试验表明，与高阶代数稳定的单对角隐式龙格-库塔方法相比，该方法具有更好的计算成本-精度权衡。该算法对刚性和大规模非线性问题具有较强的鲁棒性，同时保持了基本方案的无条件非线性稳定性和可控耗散。总之，所提出的自适应方法为模拟大规模、长时间、强非线性系统提供了一种高效、可靠和自启动的工具。

{"title":"An efficient adaptive time-integration method for second-order nonlinear dynamics","authors":"Zhi Duan , Xiaohui Liu , Hongbing Guo , Chuan Wu , Zhongfei Ye , Zhongbin Lu","doi":"10.1016/j.apm.2026.116740","DOIUrl":"10.1016/j.apm.2026.116740","url":null,"abstract":"<div><div>Achieving an optimal balance among computational efficiency, robustness, and accuracy is a central challenge in simulating second-order nonlinear dynamical systems. While the existing parameterized two-sub-step composite integrator provides rigorous nonlinear stability and controllable dissipation, its fixed-step formulation limits efficiency in simulations with strongly varying dynamics. This paper presents a novel adaptive time integration method that augments the second-order base scheme with an explicit auxiliary stage for efficient error estimation. Its key innovation is a cost-free, direct error estimator, constructed by rigorously deriving the embedding coefficients via order conditions and analytically combining the implicit base stages with an extrapolated explicit stage to derive a local error estimate based on a third-order embedding without additional nonlinear iterations or matrix operations. Combined with a proportional–integral–derivative-like step-size controller, systematic numerical tests show that the proposed method achieves a significantly better computational cost-to-accuracy trade-off than high-order algebraically stable singly diagonally implicit Runge–Kutta methods. The algorithm demonstrates strong robustness in stiff and large-scale nonlinear problems while preserving the unconditional nonlinear stability and controllable dissipation of the base scheme. In summary, the proposed adaptive method offers an efficient, reliable, and self-starting tool for simulating large-scale, long-duration, strongly nonlinear systems.</div></div>","PeriodicalId":50980,"journal":{"name":"Applied Mathematical Modelling","volume":"155 ","pages":"Article 116740"},"PeriodicalIF":4.4,"publicationDate":"2026-01-09","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145956848","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

Interval theory-embedded data-driven identification framework for uncertain thermo-elastic parameters 区间理论嵌入的数据驱动不确定热弹性参数识别框架

IF 4.4 2区工程技术 Q1 ENGINEERING, MULTIDISCIPLINARY

Applied Mathematical Modelling

Pub Date : 2026-01-09 DOI: 10.1016/j.apm.2026.116743

Xin Qiang , Chong Wang , Huanyu Zhang

Since system parameters can reflect fluctuations in structural performance, identifying the thermo-elastic parameters based on measured responses is becoming increasingly important for health monitoring of thermo-mechanical systems. To avoid the drawback of traditional probabilistic methods in handling limited experimental samples, this paper proposes a novel interval theory-integrated computational framework for efficient and robust identification of uncertain thermo-elastic parameters. For the coupled thermo-mechanical problem, the thermo-elastic governing equation is derived and the thermal stress effect is discussed. In view of the limitation of extremum searching in capturing potential supplementary data, a confidence-based unbiased interval estimation method is introduced to quantify experimental response bounds of limited experimental samples. Subsequently, a gene expression programming support vector regression (GEP-SVR) metamodel is constructed to replace the full-scale finite element simulations, thereby alleviating the computational burden of the nested dual-loop optimization in interval parameter identification. The effectiveness of the proposed framework is demonstrated through three case studies. Numerical results show that the proposed method achieves identification errors below 3.0% while improving computational efficiency by 87.08% compared to full-scale finite element simulation, providing a practical and efficient tool for uncertainty-aware parameter identification of thermo-mechanical systems.

由于系统参数可以反映结构性能的波动，因此基于测量响应识别热弹性参数对于热机械系统的健康监测变得越来越重要。为了避免传统概率方法在处理有限实验样本时的缺点，提出了一种新的区间理论集成计算框架，用于高效、鲁棒地识别不确定热弹性参数。对于热-力耦合问题，推导了热弹性控制方程，讨论了热应力效应。针对极值搜索在获取潜在补充数据方面的局限性，提出了一种基于置信度的无偏区间估计方法来量化有限实验样本的实验响应界。随后，构建基因表达编程支持向量回归（GEP-SVR）元模型取代全尺寸有限元模拟，减轻了区间参数辨识中嵌套双环优化的计算负担。通过三个案例研究证明了所提出框架的有效性。数值结果表明，与全尺寸有限元模拟相比，该方法辨识误差在3.0%以下，计算效率提高87.08%，为热机械系统的不确定性参数辨识提供了实用、高效的工具。

{"title":"Interval theory-embedded data-driven identification framework for uncertain thermo-elastic parameters","authors":"Xin Qiang , Chong Wang , Huanyu Zhang","doi":"10.1016/j.apm.2026.116743","DOIUrl":"10.1016/j.apm.2026.116743","url":null,"abstract":"<div><div>Since system parameters can reflect fluctuations in structural performance, identifying the thermo-elastic parameters based on measured responses is becoming increasingly important for health monitoring of thermo-mechanical systems. To avoid the drawback of traditional probabilistic methods in handling limited experimental samples, this paper proposes a novel interval theory-integrated computational framework for efficient and robust identification of uncertain thermo-elastic parameters. For the coupled thermo-mechanical problem, the thermo-elastic governing equation is derived and the thermal stress effect is discussed. In view of the limitation of extremum searching in capturing potential supplementary data, a confidence-based unbiased interval estimation method is introduced to quantify experimental response bounds of limited experimental samples. Subsequently, a gene expression programming support vector regression (GEP-SVR) metamodel is constructed to replace the full-scale finite element simulations, thereby alleviating the computational burden of the nested dual-loop optimization in interval parameter identification. The effectiveness of the proposed framework is demonstrated through three case studies. Numerical results show that the proposed method achieves identification errors below 3.0% while improving computational efficiency by 87.08% compared to full-scale finite element simulation, providing a practical and efficient tool for uncertainty-aware parameter identification of thermo-mechanical systems.</div></div>","PeriodicalId":50980,"journal":{"name":"Applied Mathematical Modelling","volume":"155 ","pages":"Article 116743"},"PeriodicalIF":4.4,"publicationDate":"2026-01-09","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145956847","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

Closed-form solutions for wave scattering coefficients of a Timoshenko beam with an edge crack considering axial–bending–shear coupling 考虑轴向-弯曲-剪切耦合的边裂纹Timoshenko梁波散射系数的封闭解

IF 4.4 2区工程技术 Q1 ENGINEERING, MULTIDISCIPLINARY

Applied Mathematical Modelling

Pub Date : 2026-01-09 DOI: 10.1016/j.apm.2026.116748

Taejeong Lim, Nur Indah Mukharromah, Hyun Woo Park

To date, no rigorous closed-form solutions can fully address the complexity of wave scattering phenomena, although such solutions are essential for fast and accurate assessment of edge cracks in beams. This study presents closed-form solutions for wave scattering coefficients when an incident bending wave encounters an edge crack in a Timoshenko beam, resulting in wave transmission, reflection, and mode conversion. Unlike its shallow counterpart, a deep edge crack produces complicated wave scattering phenomena due to axial–bending–shear coupling (ABSC) with increasing driving frequency. Herein, three compatibility equations for modes I and II crack in linear fracture mechanics are incorporated into the Timoshenko beam theory to account for ABSC. Closed-form solutions are derived by applying spectral solutions to the three compatibility equations and three equilibrium equations at the edge crack. The sole contribution of mode I crack and simultaneous contribution of modes I and II crack to wave scattering coefficients are thoroughly investigated with respect to the normalized frequency and crack depth ratio. Finally, the proposed closed-form solutions are validated via comparison with previously reported finite element analysis and experimental results. The findings of the study provide physical insights into complex wave scattering phenomena and are anticipated to enable faster inverse analysis in ultrasonic crack evaluation.

迄今为止，没有严格的封闭解可以完全解决波散射现象的复杂性，尽管这样的解对于快速准确地评估梁的边缘裂纹是必不可少的。本研究提出了当入射弯曲波遇到Timoshenko光束中的边缘裂纹，导致波的传输、反射和模式转换时，波散射系数的封闭解。与浅边裂纹不同，随着驱动频率的增加，深边裂纹由于轴-弯-剪切耦合而产生复杂的波散射现象。本文将线性断裂力学中I型和II型裂纹的三个相容方程纳入Timoshenko梁理论来解释ABSC。对三种相容方程和三种边缘裂纹处的平衡方程应用谱解，推导出闭型解。从归一化频率和裂纹深度比的角度深入研究了I型裂纹的唯一贡献和I、II型裂纹的同时贡献对波散射系数的影响。最后，通过与先前报道的有限元分析和实验结果的比较，验证了所提出的封闭解。研究结果为复杂的波散射现象提供了物理见解，并有望在超声裂纹评估中实现更快的逆分析。

{"title":"Closed-form solutions for wave scattering coefficients of a Timoshenko beam with an edge crack considering axial–bending–shear coupling","authors":"Taejeong Lim, Nur Indah Mukharromah, Hyun Woo Park","doi":"10.1016/j.apm.2026.116748","DOIUrl":"10.1016/j.apm.2026.116748","url":null,"abstract":"<div><div>To date, no rigorous closed-form solutions can fully address the complexity of wave scattering phenomena, although such solutions are essential for fast and accurate assessment of edge cracks in beams. This study presents closed-form solutions for wave scattering coefficients when an incident bending wave encounters an edge crack in a Timoshenko beam, resulting in wave transmission, reflection, and mode conversion. Unlike its shallow counterpart, a deep edge crack produces complicated wave scattering phenomena due to axial–bending–shear coupling (ABSC) with increasing driving frequency. Herein, three compatibility equations for modes I and II crack in linear fracture mechanics are incorporated into the Timoshenko beam theory to account for ABSC. Closed-form solutions are derived by applying spectral solutions to the three compatibility equations and three equilibrium equations at the edge crack. The sole contribution of mode I crack and simultaneous contribution of modes I and II crack to wave scattering coefficients are thoroughly investigated with respect to the normalized frequency and crack depth ratio. Finally, the proposed closed-form solutions are validated via comparison with previously reported finite element analysis and experimental results. The findings of the study provide physical insights into complex wave scattering phenomena and are anticipated to enable faster inverse analysis in ultrasonic crack evaluation.</div></div>","PeriodicalId":50980,"journal":{"name":"Applied Mathematical Modelling","volume":"155 ","pages":"Article 116748"},"PeriodicalIF":4.4,"publicationDate":"2026-01-09","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145956905","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

Complex-mode framework for sloshing damping in cylindrical tanks with porous baffles 多孔隔板圆柱罐晃动阻尼的复杂模态框架

IF 4.4 2区工程技术 Q1 ENGINEERING, MULTIDISCIPLINARY

Applied Mathematical Modelling

Pub Date : 2026-01-08 DOI: 10.1016/j.apm.2026.116742

Jiadong Wang, Xianqing Liu, Shiyu Zhao, Gang Liu

Liquid sloshing in partially filled cylindrical tanks can generate large hydrodynamic loads that threaten the safety of marine vessels, spacecraft, and storage tanks. This work introduces a complex-mode framework for liquid sloshing in cylindrical tanks equipped with annular porous baffles, based on linear potential-flow theory with a porous-jump condition. The governing equations are cast as a non-self-adjoint quadratic eigenvalue problem and solved via state-space linearization. Natural frequencies and damping ratios then follow directly from the complex eigenvalues, providing a physically based description of energy dissipation. An orthogonality relation for the complex sloshing modes is established, enabling modal decoupling of the dynamic response equations. Parametric analyses show that the sloshing characteristics are strongly influenced by the porous baffle elevation, inner radius, and porous-effect parameter. Raising the elevation or decreasing the inner radius lowers the fundamental sloshing frequency while increasing the damping ratio. For each elevation and inner radius, an optimal porous-effect parameter exists that maximizes the damping ratio of the fundamental mode. Frequency response analyses confirm that a properly tuned porous-effect parameter suppresses near-resonant sloshing responses effectively. These results provide a robust basis for sloshing mitigation in practical engineering applications.

液体在部分充满的圆柱形储罐中晃动会产生巨大的水动力载荷，威胁到船舶、航天器和储罐的安全。本文基于具有多孔跃变条件的线性势流理论，提出了一种复杂模态框架，用于研究装有环形多孔挡板的圆柱槽内液体晃动。将控制方程转化为一个非自伴随的二次特征值问题，并通过状态空间线性化进行求解。然后，固有频率和阻尼比直接遵循复特征值，提供基于物理的能量耗散描述。建立了复杂晃动模态的正交关系，实现了动力响应方程的模态解耦。参数分析表明，多孔挡板高度、内部半径和多孔效应参数对晃动特性有较大影响。提高俯仰或减小内半径降低了基频晃动频率，同时增加了阻尼比。对于每个高程和内半径，存在一个使基模态阻尼比最大化的最优多孔效应参数。频率响应分析证实，适当调整多孔效应参数可有效抑制近共振晃动响应。这些结果为在实际工程应用中减少晃动提供了坚实的基础。

{"title":"Complex-mode framework for sloshing damping in cylindrical tanks with porous baffles","authors":"Jiadong Wang, Xianqing Liu, Shiyu Zhao, Gang Liu","doi":"10.1016/j.apm.2026.116742","DOIUrl":"10.1016/j.apm.2026.116742","url":null,"abstract":"<div><div>Liquid sloshing in partially filled cylindrical tanks can generate large hydrodynamic loads that threaten the safety of marine vessels, spacecraft, and storage tanks. This work introduces a complex-mode framework for liquid sloshing in cylindrical tanks equipped with annular porous baffles, based on linear potential-flow theory with a porous-jump condition. The governing equations are cast as a non-self-adjoint quadratic eigenvalue problem and solved via state-space linearization. Natural frequencies and damping ratios then follow directly from the complex eigenvalues, providing a physically based description of energy dissipation. An orthogonality relation for the complex sloshing modes is established, enabling modal decoupling of the dynamic response equations. Parametric analyses show that the sloshing characteristics are strongly influenced by the porous baffle elevation, inner radius, and porous-effect parameter. Raising the elevation or decreasing the inner radius lowers the fundamental sloshing frequency while increasing the damping ratio. For each elevation and inner radius, an optimal porous-effect parameter exists that maximizes the damping ratio of the fundamental mode. Frequency response analyses confirm that a properly tuned porous-effect parameter suppresses near-resonant sloshing responses effectively. These results provide a robust basis for sloshing mitigation in practical engineering applications.</div></div>","PeriodicalId":50980,"journal":{"name":"Applied Mathematical Modelling","volume":"155 ","pages":"Article 116742"},"PeriodicalIF":4.4,"publicationDate":"2026-01-08","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145956849","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

Size-dependent damped free vibration of viscoelastic curved Timoshenko nanobeams resting on viscoelastic foundation through viscoelastic integral nonlocal strain gradient model 基于粘弹性积分非局部应变梯度模型的粘弹性弯曲Timoshenko纳米梁在粘弹性基础上的尺寸相关阻尼自由振动

IF 4.4 2区工程技术 Q1 ENGINEERING, MULTIDISCIPLINARY

Applied Mathematical Modelling

Pub Date : 2026-01-07 DOI: 10.1016/j.apm.2026.116739

Chang Li , Rongjun Chen , Limin Guo , Hai Qing

Viscoelastic nanobeams on viscoelastic foundations enable enhanced vibration suppression and novel functionalities in M/NEMS devices. However, accurately predicting their size-dependent, time-damped vibrations demands robust nonclassical theories. This study presents a precise theoretical framework for the free damping vibration characteristics of a curved viscoelastic Timoshenko nanobeam resting on a size-dependent viscoelastic foundation. We develop a mathematically well-posed viscoelastic integral nonlocal strain gradient theory (VINSGT) by integrating integral nonlocal strain gradient theory (INSGT with the Kelvin-Voigt model, while accounting for the size effect in the foundations’ reaction force. The integral constitutive equations are transformed into an equivalent differential form, incorporating essential constitutive boundary conditions (CBCs). The governing equations are discretized via the generalized differential quadrature method (GDQM), yielding a complex eigenvalue problem. A two-step numerical scheme resolves the vibration frequencies and establishes the relationship between damping and viscous coefficients. Numerical examples validate the VINSGT framework and systematically investigate size effects in the damped vibration behavior of viscoelastic curved Timoshenko nanobeams on size-dependent foundations. This work provides a reliable theoretical basis for designing and optimizing vibration control in advanced M/NEMS with viscoelastic nanobeam-foundation systems.

粘弹性纳米梁在粘弹性基础上增强了M/NEMS设备的振动抑制和新功能。然而，要准确预测它们的尺寸相关、时间阻尼的振动，需要强大的非经典理论。本研究提出了基于尺寸相关粘弹性基础的弯曲粘弹性Timoshenko纳米梁自由阻尼振动特性的精确理论框架。将积分非局部应变梯度理论（INSGT）与Kelvin-Voigt模型相结合，考虑地基反力中的尺寸效应，建立了数学上适定的粘弹性积分非局部应变梯度理论（VINSGT）。将积分本构方程转化为包含基本本构边界条件的等效微分形式。采用广义微分正交法（GDQM）对控制方程进行离散化，得到一个复特征值问题。采用两步数值格式求解振动频率，建立阻尼系数与粘性系数之间的关系。数值算例验证了VINSGT框架，并系统地研究了尺寸对粘弹性弯曲Timoshenko纳米梁在尺寸相关基础上的阻尼振动行为的影响。该工作为先进的粘弹性纳米梁基础M/NEMS振动控制设计和优化提供了可靠的理论依据。

{"title":"Size-dependent damped free vibration of viscoelastic curved Timoshenko nanobeams resting on viscoelastic foundation through viscoelastic integral nonlocal strain gradient model","authors":"Chang Li , Rongjun Chen , Limin Guo , Hai Qing","doi":"10.1016/j.apm.2026.116739","DOIUrl":"10.1016/j.apm.2026.116739","url":null,"abstract":"<div><div>Viscoelastic nanobeams on viscoelastic foundations enable enhanced vibration suppression and novel functionalities in M/NEMS devices. However, accurately predicting their size-dependent, time-damped vibrations demands robust nonclassical theories. This study presents a precise theoretical framework for the free damping vibration characteristics of a curved viscoelastic Timoshenko nanobeam resting on a size-dependent viscoelastic foundation. We develop a mathematically well-posed viscoelastic integral nonlocal strain gradient theory (VINSGT) by integrating integral nonlocal strain gradient theory (INSGT with the Kelvin-Voigt model, while accounting for the size effect in the foundations’ reaction force. The integral constitutive equations are transformed into an equivalent differential form, incorporating essential constitutive boundary conditions (CBCs). The governing equations are discretized via the generalized differential quadrature method (GDQM), yielding a complex eigenvalue problem. A two-step numerical scheme resolves the vibration frequencies and establishes the relationship between damping and viscous coefficients. Numerical examples validate the VINSGT framework and systematically investigate size effects in the damped vibration behavior of viscoelastic curved Timoshenko nanobeams on size-dependent foundations. This work provides a reliable theoretical basis for designing and optimizing vibration control in advanced M/NEMS with viscoelastic nanobeam-foundation systems.</div></div>","PeriodicalId":50980,"journal":{"name":"Applied Mathematical Modelling","volume":"155 ","pages":"Article 116739"},"PeriodicalIF":4.4,"publicationDate":"2026-01-07","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145956904","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

Dynamic modelling and vortex dynamics in elastic-joint caudal fin propulsor for efficient swimming 高效游泳弹性关节尾鳍推进器的动力学建模和涡动力学

IF 4.4 2区工程技术 Q1 ENGINEERING, MULTIDISCIPLINARY

Applied Mathematical Modelling

Pub Date : 2026-01-05 DOI: 10.1016/j.apm.2025.116737

Yiming Lu , Haicheng Zhang , Xubin Sun , Daolin Xu

Marine animals achieve remarkable swimming performance through precisely coordinated interactions among different body segments. Current bio-inspired swimmers generally fail to account for the distinct mechanical functions of various anatomical structures and often overlook the critical nonlinear fluid-structure interaction mechanisms. To bridge this gap, an elastic-joint rigid fin model is developed that conceptualizes propulsion as three synergistic phases: active muscular actuation, elastic peduncle energy transmission, and vortex-generating caudal fin motion. This model uniquely enables the quantitative dissection of individual component contributions to overall performance. To address the fluid-structure interaction problem in swimming, a Nonlinear Vortex Sheet Method (NVSM) is proposed. The method seamlessly integrates a vortex sheet model based flow reconstruction, a coupled system modeling for structural motion and vortex tracking, and the Broyden method based synchronized solving, achieving robust and efficient performance. Rigorously validated against computational fluid dynamics, the NVSM accurately captures the key physics of dynamic response and vortex shedding with significantly reduced computational cost.

Comparative analyses demonstrate a threefold enhancement in peak propulsion efficiency at optimal frequencies over rigid counterparts. A detailed examination of the fluid-structure interaction reveals the propulsion mechanisms. These are identified as hydrodynamic phase synchronization, suppression of energy-dissipating vortex structures, and superior fluid energy utilization efficiency. Furthermore, this study establishes a quantitative mapping between structural stiffness and key performance metrics, providing actionable design guidelines for next-generation biomimetic propulsion systems.

海洋动物通过不同身体部位之间精确协调的相互作用，获得卓越的游泳性能。目前的仿生游泳者普遍未能考虑到各种解剖结构的独特力学功能，往往忽视了关键的非线性流固相互作用机制。为了弥补这一差距，研究人员开发了一种弹性关节刚性鳍模型，该模型将推进概念定义为三个协同阶段：主动肌肉驱动、弹性肢节能量传输和尾鳍运动产生涡流。该模型独特地能够定量分析单个组件对整体性能的贡献。针对游泳过程中流固耦合问题，提出了一种非线性涡片法（NVSM）。该方法将基于旋涡片模型的流动重建、结构运动与旋涡跟踪的耦合系统建模和基于Broyden方法的同步求解无缝集成，实现了鲁棒性和高效率。NVSM经过了严格的计算流体动力学验证，能够准确地捕捉到动态响应和涡流脱落的关键物理特性，大大降低了计算成本。对比分析表明，在最佳频率下，峰值推进效率比刚性频率提高了三倍。对流固相互作用的详细研究揭示了推进机制。这些特性被确定为流体动力相位同步，抑制能量耗散涡结构，以及优越的流体能量利用效率。此外，该研究建立了结构刚度和关键性能指标之间的定量映射，为下一代仿生推进系统提供了可操作的设计指南。

{"title":"Dynamic modelling and vortex dynamics in elastic-joint caudal fin propulsor for efficient swimming","authors":"Yiming Lu , Haicheng Zhang , Xubin Sun , Daolin Xu","doi":"10.1016/j.apm.2025.116737","DOIUrl":"10.1016/j.apm.2025.116737","url":null,"abstract":"<div><div>Marine animals achieve remarkable swimming performance through precisely coordinated interactions among different body segments. Current bio-inspired swimmers generally fail to account for the distinct mechanical functions of various anatomical structures and often overlook the critical nonlinear fluid-structure interaction mechanisms. To bridge this gap, an elastic-joint rigid fin model is developed that conceptualizes propulsion as three synergistic phases: active muscular actuation, elastic peduncle energy transmission, and vortex-generating caudal fin motion. This model uniquely enables the quantitative dissection of individual component contributions to overall performance. To address the fluid-structure interaction problem in swimming, a Nonlinear Vortex Sheet Method (NVSM) is proposed. The method seamlessly integrates a vortex sheet model based flow reconstruction, a coupled system modeling for structural motion and vortex tracking, and the Broyden method based synchronized solving, achieving robust and efficient performance. Rigorously validated against computational fluid dynamics, the NVSM accurately captures the key physics of dynamic response and vortex shedding with significantly reduced computational cost.</div><div>Comparative analyses demonstrate a threefold enhancement in peak propulsion efficiency at optimal frequencies over rigid counterparts. A detailed examination of the fluid-structure interaction reveals the propulsion mechanisms. These are identified as hydrodynamic phase synchronization, suppression of energy-dissipating vortex structures, and superior fluid energy utilization efficiency. Furthermore, this study establishes a quantitative mapping between structural stiffness and key performance metrics, providing actionable design guidelines for next-generation biomimetic propulsion systems.</div></div>","PeriodicalId":50980,"journal":{"name":"Applied Mathematical Modelling","volume":"155 ","pages":"Article 116737"},"PeriodicalIF":4.4,"publicationDate":"2026-01-05","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145902634","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 new meshless ENO scheme based on FVS methods for compressible flows 基于FVS方法的可压缩流无网格ENO新方案

IF 4.4 2区工程技术 Q1 ENGINEERING, MULTIDISCIPLINARY

Applied Mathematical Modelling

Pub Date : 2026-01-03 DOI: 10.1016/j.apm.2025.116732

Tao Zhang , Xiaofeng Zhou , Jian Li

We propose a meshless essentially non-oscillatory (MENO) scheme based on flux vector splitting (FVS) methods that can be integrated with various meshless methods to solve compressible flows problems involving upwind bias and discontinuities. The central idea is to generate multiple generalized finite difference method (GFDM) approximations at each node by employing asymmetric or eccentric stars. The smoothness indicator is then applied to each candidate approximation to select the optimal one for suppressing oscillations at discontinuities, in a manner similar to the Essentially Non-Oscillatory (ENO) scheme. First, based on the meshless methods, the Steger-Warming flux splitting method for the Euler equations is developed; Second, a one-side upwind stencil is constructed; Third, two additional stencils are constructed along the upwind or downwind direction to identify discontinuities; Finally, the stencil with the smallest smoothness indicator is selected as the optimal one. It seems that this is the first attempt to develop a shock-capturing scheme in the meshless FVS methods. Compared to the meshless FDS methods, the proposed meshless FVS methods combined with MENO scheme is simpler and somewhat more efficient, as it identifies the upwind direction with less effort and does not require mid-point reconstruction. Several benchmark tests demonstrate that the proposed MENO scheme under meshless FVS methods effectively achieves second- or third-order accuracy in smooth regions while maintaining robust and non-oscillatory shock-capturing, representing a significant improvement over the upwind GFDM scheme.

我们提出了一种基于通量矢量分裂（FVS）方法的无网格基本无振荡（MENO）方案，该方案可以与各种无网格方法相结合，以解决涉及逆风偏置和不连续的可压缩流动问题。其核心思想是利用不对称星或偏心星在每个节点生成多个广义有限差分法（GFDM）逼近。然后将平滑指标应用于每个候选近似，以选择最优的近似来抑制不连续处的振荡，类似于本质非振荡（ENO）方案。首先，在无网格方法的基础上，提出了Euler方程的Steger-Warming通量分裂法；其次，构造单侧逆风模板；第三，沿逆风或下风方向构建两个附加模板以识别不连续面；最后，选取平滑度指标最小的模板作为最优模板。这似乎是在无网格FVS方法中开发冲击捕获方案的第一次尝试。与无网格FDS方法相比，本文提出的无网格FVS与MENO方案相结合的方法更简单，效率更高，因为它可以更轻松地识别逆风方向，并且不需要重建中点。几个基准测试表明，在无网格FVS方法下提出的MENO方案在光滑区域有效地实现了二阶或三阶精度，同时保持了鲁棒性和非振荡冲击捕获，比逆风GFDM方案有了显着改进。

{"title":"A new meshless ENO scheme based on FVS methods for compressible flows","authors":"Tao Zhang , Xiaofeng Zhou , Jian Li","doi":"10.1016/j.apm.2025.116732","DOIUrl":"10.1016/j.apm.2025.116732","url":null,"abstract":"<div><div>We propose a meshless essentially non-oscillatory (MENO) scheme based on flux vector splitting (FVS) methods that can be integrated with various meshless methods to solve compressible flows problems involving upwind bias and discontinuities. The central idea is to generate multiple generalized finite difference method (GFDM) approximations at each node by employing asymmetric or eccentric stars. The smoothness indicator is then applied to each candidate approximation to select the optimal one for suppressing oscillations at discontinuities, in a manner similar to the Essentially Non-Oscillatory (ENO) scheme. First, based on the meshless methods, the Steger-Warming flux splitting method for the Euler equations is developed; Second, a one-side upwind stencil is constructed; Third, two additional stencils are constructed along the upwind or downwind direction to identify discontinuities; Finally, the stencil with the smallest smoothness indicator is selected as the optimal one. It seems that this is the first attempt to develop a shock-capturing scheme in the meshless FVS methods. Compared to the meshless FDS methods, the proposed meshless FVS methods combined with MENO scheme is simpler and somewhat more efficient, as it identifies the upwind direction with less effort and does not require mid-point reconstruction. Several benchmark tests demonstrate that the proposed MENO scheme under meshless FVS methods effectively achieves second- or third-order accuracy in smooth regions while maintaining robust and non-oscillatory shock-capturing, representing a significant improvement over the upwind GFDM scheme.</div></div>","PeriodicalId":50980,"journal":{"name":"Applied Mathematical Modelling","volume":"155 ","pages":"Article 116732"},"PeriodicalIF":4.4,"publicationDate":"2026-01-03","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145894396","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

Hierarchical multi-agent reinforcement learning algorithm for multi-UAV roundup strategy 多无人机围捕策略的分层多智能体强化学习算法

IF 4.4 2区工程技术 Q1 ENGINEERING, MULTIDISCIPLINARY

Applied Mathematical Modelling

Pub Date : 2026-01-03 DOI: 10.1016/j.apm.2025.116728

Qichen Xu, Zhaohui Zhang, Jing Li, Xiaogang Qi

This paper addresses collision-avoidance multi-UAV roundup by introducing a hierarchical instructor-student reinforcement learning framework that decomposes high-dimensional decisions into discrete orders. An instructor coordinates global context, issues approach/stretch/contract commands, and is guided by a reward machine capturing task-stage progress, while students learn continuous control under these commands. Compared with a centralized training and decentralized execution baseline, the proposed scheme converges faster, reaches higher cumulative rewards, and improves success rate with no observed collisions in simulations. It also maintains robust roundup when obstacles move, indicating practical potential for swarm containment missions.

本文通过引入分层师生强化学习框架，将高维决策分解为离散顺序，解决了多无人机避撞问题。教练协调全局环境，发出接近/拉伸/收缩命令，并由奖励机器引导，捕捉任务阶段的进展，而学生在这些命令下学习连续控制。与集中式训练和分散式执行基线相比，该方案收敛速度更快，获得更高的累积奖励，并且在模拟中没有观察到碰撞的情况下提高了成功率。当障碍物移动时，它还能保持强大的围捕，这表明了蜂群控制任务的实际潜力。

引用次数: 0