{"title":"Radiation force and torque of elastic compressional Bessel waves on a solid sphere embedded in an unbounded elastic medium","authors":"Yuchen Zang , F.G. Mitri","doi":"10.1016/j.apm.2025.115958","DOIUrl":"10.1016/j.apm.2025.115958","url":null,"abstract":"<div><div>The time-averaged elastic radiation force and spin torque exerted on a solid viscoelastic sphere embedded in an unbounded elastic medium are considered, with the incident field composed of elastic compressional Bessel non-vortex or vortex progressive waves. Based on the multipole expansion method using spherical wave functions, partial-wave series expressions are derived for the elastic radiation force and torque through the integration of the elastodynamic Poynting vector as well as the cross product of the position vector and the time-averaged elastic radiation stress tensor. The dimensionless absorption, scattering and extinction efficiencies are also calculated. Numerical computations are performed for a brass sphere in a soft elastic gel matrix to illustrate the analysis with particular emphasis on varying the dimensionless size parameter of the sphere, the half-cone angle and order of the incident Bessel waves. The component related to mode preservation contributes dominantly to the total radiation force, while the component related to mode conversion alternates between positive, negative and neutral values. The elastic radiation torque components related to mode preservation and conversion are opposite in sign and approximately equal in magnitude. The results may find potential applications in the activation of implantable spherical devices, characterization of biological tissue, elastic wave scattering, non-destructive evaluation, and geophysical prospecting to name some examples.</div></div>","PeriodicalId":50980,"journal":{"name":"Applied Mathematical Modelling","volume":"142 ","pages":"Article 115958"},"PeriodicalIF":4.4,"publicationDate":"2025-01-19","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143174351","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}
Tengfei Tang , Hanliang Fang , Fufu Yang , Jun Zhang
{"title":"Adaptive Kriging elastodynamic modelling and analysis for a redundantly actuated parallel manipulator","authors":"Tengfei Tang , Hanliang Fang , Fufu Yang , Jun Zhang","doi":"10.1016/j.apm.2025.115961","DOIUrl":"10.1016/j.apm.2025.115961","url":null,"abstract":"<div><div>This paper presents a quantitative analysis comparing the elastodynamic performance of a redundantly actuated parallel manipulator to its non-redundantly actuated counterpart. A unified elastodynamic model is established to encompasses both actuation types, facilitating a comprehensive analysis of their performance characteristics. To assess the impact of redundant actuation, frequency variation indices are defined to measure the influence of the redundantly actuated limbs on the overall system dynamics. Two global sensitivity indices are introduced to quantify how variations in design variables affect elastodynamic performance. An adaptive Kriging-based algorithm is developed to predict elastodynamic behavior, and it is numerically validated for both accuracy and efficiency. Through a detailed elastodynamic comparison and sensitivity analysis, we identify the key effects of redundant actuation and pinpoint the most significant design variables. Our findings reveal that the enhancements in low-order natural frequencies due to redundant actuation are configuration- and order-dependent. The redundantly actuated parallel manipulator demonstrates heightened sensitivity to geometric variables in terms of elastodynamic performance, surpassing that of its non-redundantly actuated variant. The sensitivity analysis provides valuable insights, guiding future improvements in elastodynamic performance. Notably, the thickness of the moving platform emerges as a critical parameter that warrants optimization during the design phase. Therefore, this analysis is particularly pivotal in the initial evaluation stages of parallel manipulators, ensuring that designs are refined for enhanced dynamic responses. With necessary modifications, this work can be applied to other redundantly actuated parallel manipulators to contribute a deeper understanding of how redundancy influences their dynamic characteristics.</div></div>","PeriodicalId":50980,"journal":{"name":"Applied Mathematical Modelling","volume":"142 ","pages":"Article 115961"},"PeriodicalIF":4.4,"publicationDate":"2025-01-18","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143035294","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":"Analytical solution of shallow arbitrarily shaped tunnels in fractional viscoelastic transversely isotropic strata","authors":"Zhi Yong Ai, Lei Yang, Zi Kun Ye, Da Shan Wang","doi":"10.1016/j.apm.2025.115951","DOIUrl":"10.1016/j.apm.2025.115951","url":null,"abstract":"<div><div>An analytical method is proposed to obtain the stress and displacement around shallow arbitrarily shaped tunnels excavated in viscoelastic transversely isotropic strata. Firstly, combined with the fractional calculus theory, the fractional viscoelastic half-plane solution is obtained by displacement modification. Then, based on the solution of the half-plane, the Schwartz alternating method and the conformal mapping technique are introduced to obtain the solution of the fractional viscoelastic arbitrarily shaped shallow tunnels. According to the proposed theory, a MATLAB program is developed, and the accuracy of the theory is proved by degraded comparison with results of the existing literature and finite element software ABAQUS. Finally, numerical examples are designed to discuss the influence of tunnel shapes, transverse isotropy and viscoelastic parameters on the boundary stress and displacement around shallow arbitrarily shaped tunnels.</div></div>","PeriodicalId":50980,"journal":{"name":"Applied Mathematical Modelling","volume":"141 ","pages":"Article 115951"},"PeriodicalIF":4.4,"publicationDate":"2025-01-18","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143035297","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 theoretical model for a low-frequency two-stage hybrid vibration isolator with a nonlinear energy sink and a negative stiffness spring","authors":"Xingbao Huang , Biao Wang , Zhiwen Huang , Xugang Hua , Zhengqing Chen","doi":"10.1016/j.apm.2025.115948","DOIUrl":"10.1016/j.apm.2025.115948","url":null,"abstract":"<div><div>This paper proposes a new design philosophy of integrating a nonlinear energy sink (NES) with a negative stiffness spring (NSS) to construct a two-stage hybrid vibration isolator for low-frequency vibration attenuation. A theoretical model of the proposed two-stage hybrid vibration isolator is established to investigate the vibration isolation performance. The effect of a linear vibration absorber (LVA), the NES, and the NSS on the force transmissibility of the structure is compared, in which the parametric study including negative stiffness coefficient, mass ratio, frequency ratio, and barrier height is conducted. The nonlinear dynamic characteristics and energy transfer efficiency of the vibration isolation system are analyzed based on motion phase trajectories, Poincaré mapping and chaotic bifurcation analysis. The results demonstrate that the two-stage isolator with the NES and NSS architecture exhibits superior broadband vibration isolation compared to its counterpart with the LVA and NSS. The structure with the LVA and NSS shows a wider forbidden band range than the device without the LVA. Besides, the negative stiffness configuration between the vibration source and the protected base is a critical factor to influence vibration isolation. A larger negative stiffness coefficient results in lower force transmissibility and the onset frequency of the forbidden zone but the broader bandwidth of the forbidden zone. This work paves a new theoretical venue for low-frequency hybrid vibration isolators in precision instrument fields.</div></div>","PeriodicalId":50980,"journal":{"name":"Applied Mathematical Modelling","volume":"142 ","pages":"Article 115948"},"PeriodicalIF":4.4,"publicationDate":"2025-01-18","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143035296","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":"Dynamic stiffness method for exact modelling of acoustics black hole beams using Meijer G-functions","authors":"Le Chang, Li Cheng","doi":"10.1016/j.apm.2025.115945","DOIUrl":"10.1016/j.apm.2025.115945","url":null,"abstract":"<div><div>The acoustics black hole (ABH) effect shows promising potential for wave manipulation and vibration control. An ABH structure features a gradual reduction of the phase velocity of flexural waves alongside wave compression and energy accumulation when entering the tapered ABH portion where the thickness is tailored according to a power-law (with power index <em>m</em> no less than 2). The corresponding non-uniform wavelength distribution over the ABH structure poses great challenges to conventional modelling methods. To alleviate the problem, this paper proposes an exact dynamic stiffness method for modelling ABH beams with arbitrary exponent equal to or greater than 2 under the framework of Euler-Bernoulli beam theory. For ABH with <em>m</em> > 2, a change of variable and the Mellin integral transformation are conducted to derive the integral representations of the exact solution using Meijer G-functions. The solution for the case with <em>m</em> = 2 is also derived for completeness. Then the dynamic stiffness matrix is formulated through symbolic operation. The Wittrick-Williams (WW) algorithm is revamped to cope with the ABH-specific requirement. Numerical examples are given to validate the solution in integral form, the dynamic stiffness matrix, and the efficacy of the improved WW algorithm. The clear advantage of the accurate integral representations over series representations is justified in the higher frequency range. Covering all ABH-relevant scenarios (with <em>m</em> ≥ 2), the exact modelling framework established in this work offers a powerful tool for the modelling and investigation of more complex structures which are built upon ABH beam elements.</div></div>","PeriodicalId":50980,"journal":{"name":"Applied Mathematical Modelling","volume":"142 ","pages":"Article 115945"},"PeriodicalIF":4.4,"publicationDate":"2025-01-17","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143174366","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}
Leilei Chen , Haojie Lian , Chengmiao Liu , Yongsong Li , Sundararajan Natarajan
{"title":"Sensitivity analysis of transverse electric polarized electromagnetic scattering with isogeometric boundary elements accelerated by a fast multipole method","authors":"Leilei Chen , Haojie Lian , Chengmiao Liu , Yongsong Li , Sundararajan Natarajan","doi":"10.1016/j.apm.2025.115956","DOIUrl":"10.1016/j.apm.2025.115956","url":null,"abstract":"<div><div>With the capability of integrating computer-aided engineering and computer-aided design, isogeometric boundary element method has shown its great potential in electromagnetic scattering simulation. However, it has large memory requirements and computational costs because of its fully populated and unsymmetric system matrix. To address this issue, this paper derives fast multipole formulation for transverse electric polarized electromagnetic scattering sensitivity analysis with isogeometric boundary element methods, where the gradients of radar cross section are evaluated with respect to shape design variables and incident wave propagation angles. The effectiveness and reliability of the proposed algorithm are validated through numerical examples.</div></div>","PeriodicalId":50980,"journal":{"name":"Applied Mathematical Modelling","volume":"141 ","pages":"Article 115956"},"PeriodicalIF":4.4,"publicationDate":"2025-01-17","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143035299","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":"Numerical discretization errors of fluid-structure interaction model in pressurized pipeline systems","authors":"Hai Huang , Pengcheng Guo , Jianguo Yan","doi":"10.1016/j.apm.2025.115954","DOIUrl":"10.1016/j.apm.2025.115954","url":null,"abstract":"<div><div>The fluid-structure interaction (FSI) effects in pressurized pipeline systems may induce pipe vibrations, leading to more severe water hammer incidents. Therefore, establishing an FSI model for liquid-filled pipelines and developing efficient and accurate numerical simulation methods are crucial. These advancements enable precise prediction of FSI water hammer phenomena, which is essential for the design and protection of pressurized pipeline systems. When using the FSI four-equation axial model for water hammer calculations, techniques such as space-line interpolation (SLI), time-line interpolation (TLI), or wave speed adjustment (WSA) are commonly employed for grid processing. However, these techniques inevitably introduce additional numerical discretization errors. Existing literatures have primarily conducted quantitative analyses by comparing numerical results of different grid processing techniques with exact solutions or experimental data, lacking a theoretical analysis. To address this gap, this study establishes an equivalent hyperbolic differential equations (EHDE) approach for numerical discretization error analysis of the FSI four-equation axial model. EHDEs for SLI, TLI, and WSA are specifically derived to analyze how these techniques generate numerical discretization errors and the factors influencing these errors during the solution process. Theoretical findings are validated through numerical case studies. The results indicate that SLI introduces non-physical numerical dissipation terms into the fluid equations of the FSI four-equation axial model. TLI transforms a single fluid pressure wave into a superposition of two pressure waves with different wave speeds, causing additional numerical dissipation and dispersion due to this non-physical superposition. WSA introduces numerical errors when adjusting wave speeds artificially. Among the three grid processing techniques, SLI exhibits the highest sensitivity to grid resolution, followed by TLI, while WSA shows the least sensitivity. Under the identical spatial grid conditions, SLI yields the lowest computational accuracy, with the largest mean absolute error (MAE) compared to the recursive exact solution. TLI follows next, whereas the WSA demonstrates the smallest MAE. These findings enhance the theoretical analysis of numerical discretization errors in the FSI four-equation axial model, providing theoretical support for selecting grid processing techniques during the solution process.</div></div>","PeriodicalId":50980,"journal":{"name":"Applied Mathematical Modelling","volume":"142 ","pages":"Article 115954"},"PeriodicalIF":4.4,"publicationDate":"2025-01-17","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143175296","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":"Distributed optimization for penalized regression in massive compositional data","authors":"Yue Chao , Lei Huang , Xuejun Ma","doi":"10.1016/j.apm.2025.115950","DOIUrl":"10.1016/j.apm.2025.115950","url":null,"abstract":"<div><div>Compositional data have been widely used in various fields to analyze parts of a whole, providing insights into proportional relationships. With the increasing availability of extraordinarily large compositional datasets, addressing the challenges of distributed statistical methodologies and computations has become essential in the era of big data. This paper focuses on the optimization methodology and practical application of the distributed sparse penalized linear log-contrast model for massive compositional data, specifically in the context of medical insurance reimbursement ratio prediction. We propose two distributed optimization techniques tailored for centralized and decentralized topologies to effectively tackle the constrained convex optimization problems that arise in this application. Our algorithms are rooted in the frameworks of the alternating direction method of multipliers and the coordinate descent method of multipliers, making them available for distributed data scenarios. Notably, in the decentralized topology, we introduce a distributed coordinate-wise descent algorithm that employs a group alternating direction method of multipliers to achieve efficient distributed regularized estimation. We rigorously present convergence analysis for our decentralized algorithm, ensuring its reliability for practical applications. Through numerical experiments on both simulated datasets and a real-world medical insurance dataset, we evaluate the performance of our proposed algorithms.</div></div>","PeriodicalId":50980,"journal":{"name":"Applied Mathematical Modelling","volume":"141 ","pages":"Article 115950"},"PeriodicalIF":4.4,"publicationDate":"2025-01-16","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143035300","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}
Inhye Bang , Byung-In Kim , Jeongeol Park , Gwangsoo Kim
{"title":"Integrated cutting stock and multi-period inventory optimization considering raw material–product eligibility for steel-pipe manufacturers","authors":"Inhye Bang , Byung-In Kim , Jeongeol Park , Gwangsoo Kim","doi":"10.1016/j.apm.2025.115953","DOIUrl":"10.1016/j.apm.2025.115953","url":null,"abstract":"<div><div>This study addresses the complex cutting stock problem encountered by steel-pipe manufacturers. Key features of this problem include raw material–product eligibility, multi-period constraints, and inventory level restrictions. The primary objective is to allocate product orders to eligible raw materials while optimizing the production plan over multiple time periods. To tackle this challenge, we propose a comprehensive mathematical model with an objective function that minimizes raw-material costs, product inventory levels, the use of over-specified raw materials, unfulfilled orders, and excess inventory while adhering to eligibility constraints and yield limitations. To solve the model, we introduce a large neighborhood search algorithm integrated with a heuristic initial solution construction, destroy operators, and repair operators. A series of experiments on randomly generated instances and real-world data demonstrate the algorithm's effectiveness, outperforming conventional mathematical model-based approaches. For large-scale real-world problems, our method achieved a 24.8 % reduction in cost and a 14.0 % improvement in inventory optimization compared to current company practices.</div></div>","PeriodicalId":50980,"journal":{"name":"Applied Mathematical Modelling","volume":"142 ","pages":"Article 115953"},"PeriodicalIF":4.4,"publicationDate":"2025-01-16","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143035301","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":"Aeroelasticity of membrane airfoils and flexible-chord airfoils with permeable trailing sections","authors":"Omar S. Hussein","doi":"10.1016/j.apm.2025.115947","DOIUrl":"10.1016/j.apm.2025.115947","url":null,"abstract":"<div><div>This work addresses the aeroelasticity of airfoils with permeable trailing sections which is a topic that has not been studied in the literature before. There has been an increasing interest in the past few years regarding the analysis of permeable airfoils from the aerodynamics and aeroacoustics points of view. So, this work adds a new aspect to the analysis by considering the airfoil elasticity and flexibility effect. Two types of airfoils are considered: membrane airfoils and flexible-chord airfoils. Both static and dynamic aeroelastic studies of permeable airfoils are presented through a parametric analysis to investigate the effects of the permeable section length, permeability factor, boundary conditions represented by the membrane fixed or elastic supports and the elastic center position of flexible-chord airfoils, and the elastic properties of the airfoils like the bending rigidity of flexible-chord airfoils and the tension coefficient of membrane airfoils. The aeroelastic model of permeable airfoils is derived using the Hamilton's variational principle and the finite element method. The finite element system of equations is presented in a time domain nonlinear form to account for the nonlinear stiffening of membrane airfoils. The seepage flow through the permeable section is modeled using Darcy's law. The effects of permeability on the static aerodynamic coefficients, the dynamic response under a sinusoidal gust, and the aeroelastic instabilities (divergence and flutter) are investigated. The study aims at the aeroelastic performance of permeable airfoils subjected to incompressible airflows at small angles of attack and moderate to high Reynolds numbers, so the potential flow theory is adopted for the aerodynamic modeling. The analysis shows that permeability can have a significant effect on the aeroelastic responses of airfoils. The trailing permeable section changes the pressure distribution which changes the lift curve slope and the aerodynamic center position which in turn affects the aeroelastic instabilities and responses. This permeability effect strongly depends on the airfoil's stiffness represented by the bending rigidity of the flexible-chord airfoils and the tension in the membrane airfoils.</div></div>","PeriodicalId":50980,"journal":{"name":"Applied Mathematical Modelling","volume":"141 ","pages":"Article 115947"},"PeriodicalIF":4.4,"publicationDate":"2025-01-15","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143178391","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}