Pub Date : 2026-02-02DOI: 10.1016/j.mechrescom.2026.104634
Isidora Rapajić , Srboljub Simić
This work proposes an ODE model for a capillary rise in pipes with variable cross section, and compares it to the lubrication theory model. Two key assumptions are made: (1) radius of the pipe varies with axial coordinate, and (2) pipe’s convergence angle is small. The model reduction process involves the identification of critical parameters and simplifies the governing equations by neglecting higher-order terms. Under appropriate scaling, it is shown that generalized Washburn’s equation for capillary rise in pipes with variable cross section reduces to the lubrication theory model (Figliuzzi and Buie, 2013; Gorce et al., 2016).
本文提出了变截面管道中毛细上升的ODE模型,并与润滑理论模型进行了比较。提出两个关键假设:(1)管道半径随轴向坐标变化,(2)管道收敛角较小。模型约简过程包括关键参数的识别,并通过忽略高阶项来简化控制方程。在适当的标度下,变截面管道中毛细上升的广义Washburn方程可以简化为润滑理论模型(Figliuzzi and Buie, 2013; Gorce et al., 2016)。
{"title":"Capillary rise in pipes with variable cross section","authors":"Isidora Rapajić , Srboljub Simić","doi":"10.1016/j.mechrescom.2026.104634","DOIUrl":"10.1016/j.mechrescom.2026.104634","url":null,"abstract":"<div><div>This work proposes an ODE model for a capillary rise in pipes with variable cross section, and compares it to the lubrication theory model. Two key assumptions are made: (1) radius of the pipe varies with axial coordinate, and (2) pipe’s convergence angle is small. The model reduction process involves the identification of critical parameters and simplifies the governing equations by neglecting higher-order terms. Under appropriate scaling, it is shown that generalized Washburn’s equation for capillary rise in pipes with variable cross section reduces to the lubrication theory model (Figliuzzi and Buie, 2013; Gorce et al., 2016).</div></div>","PeriodicalId":49846,"journal":{"name":"Mechanics Research Communications","volume":"153 ","pages":"Article 104634"},"PeriodicalIF":2.3,"publicationDate":"2026-02-02","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"146116727","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2026-01-20DOI: 10.1016/j.mechrescom.2026.104633
Na Liu , Guoquan Nie
This study focuses on the dispersion characteristics of shear horizontal (SH) waves in a sandwich plate that incorporates a piezoelectric (PE) core layer and two piezomagnetic (PM) surface layers. The main objective is to reveal the effects of constituent PM materials and magnetic boundary conditions at the plate surface on dispersion behavior. Two different material combinations of PE and PM are considered. A more realistic assumption of magnetoelectric (ME) interface condition is applied and two different magnetic boundary conditions at the plate surface are discussed. Six types of phase velocity dispersion of SH waves resulted from the competitive relation of stiffness of PE and PM constituents are shown. The derived solutions are verified through finite element simulation, and a good agreement is demonstrated. The dispersion characteristics of SH waves in two typical ME sandwich plates are illustrated through several numerical examples. The effects of PM material property and magnetic boundary conditions on dispersion behaviors are revealed. The dependence of phase velocity of SH wave on volume fraction of components are shown. Some important characteristics of SH wave dispersion are drawn. The obtained results may be useful for the analysis and design of ME-based energy conversion devices.
{"title":"The effect of constituent materials on dispersion behaviors of SH waves in magnetoelectric sandwich plates","authors":"Na Liu , Guoquan Nie","doi":"10.1016/j.mechrescom.2026.104633","DOIUrl":"10.1016/j.mechrescom.2026.104633","url":null,"abstract":"<div><div>This study focuses on the dispersion characteristics of shear horizontal (SH) waves in a sandwich plate that incorporates a piezoelectric (PE) core layer and two piezomagnetic (PM) surface layers. The main objective is to reveal the effects of constituent PM materials and magnetic boundary conditions at the plate surface on dispersion behavior. Two different material combinations of PE and PM are considered. A more realistic assumption of magnetoelectric (ME) interface condition is applied and two different magnetic boundary conditions at the plate surface are discussed. Six types of phase velocity dispersion of SH waves resulted from the competitive relation of stiffness of PE and PM constituents are shown. The derived solutions are verified through finite element simulation, and a good agreement is demonstrated. The dispersion characteristics of SH waves in two typical ME sandwich plates are illustrated through several numerical examples. The effects of PM material property and magnetic boundary conditions on dispersion behaviors are revealed. The dependence of phase velocity of SH wave on volume fraction of components are shown. Some important characteristics of SH wave dispersion are drawn. The obtained results may be useful for the analysis and design of ME-based energy conversion devices.</div></div>","PeriodicalId":49846,"journal":{"name":"Mechanics Research Communications","volume":"152 ","pages":"Article 104633"},"PeriodicalIF":2.3,"publicationDate":"2026-01-20","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"146078982","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2026-01-19DOI: 10.1016/j.mechrescom.2026.104632
Feng Guo , Tao Fu , Chengfei Li , Jianfeng Jiang , Jinxiang Fang , Sen Wang , Chao Yang
Three-dimensional (3D) zero Poisson’s ratio (ZPR) metamaterials remain hold significant research potential, despite extensive research on two-dimensional (2D) ZPR and negative Poisson’s ratio (NPR) metamaterials in recent years. These materials exhibit extraordinarily mechanical properties that provide substantial advantages and broader application prospects over traditional materials. Therefore, in this paper, an orthogonal 3D orthogonal U-shaped accordion (3D-OUA) metamaterial was proposed on the basis of the 2D classical re-entrant hexagon core (RHC) metamaterial. It was applied to the modeling of honeycomb structures, sandwich panels and energy absorbing boxes, and the dynamic response under mechanical compression and dynamic impact was investigated. Finite element models (FEMs) were developed using ABAQUS. Specimens were fabricated using a metal 3D printer. Quasi-static compression experiments and low-velocity impact experiments were conducted, and the dynamic Poisson’s ratio of the honeycomb structure was filmed by a high-speed camera and image processing techniques. Decreasing the amplitude of the cosine function of the unit cell, increasing the tensile thickness of the unit cell and employing multidirectional unit cells to construct the 3D-OUA honeycomb structure can enhance the energy absorption and plateau stress. Reducing the cosine function amplitude of the unit cell, increasing the tensile thickness of the unit cell and employing multidirectional unit cells can enhance the impact resistance of 3D-OUA sandwich panels. In addition, 3D-OUA has obvious advantages over the classical energy absorbing box filling structure and 3D-RHC structure is expected to be applied to the design of new energy absorbing boxes.
{"title":"Dynamic response of 3D orthogonal U-shaped accordion honeycomb structures with zero Poisson’s ratio","authors":"Feng Guo , Tao Fu , Chengfei Li , Jianfeng Jiang , Jinxiang Fang , Sen Wang , Chao Yang","doi":"10.1016/j.mechrescom.2026.104632","DOIUrl":"10.1016/j.mechrescom.2026.104632","url":null,"abstract":"<div><div>Three-dimensional (3D) zero Poisson’s ratio (ZPR) metamaterials remain hold significant research potential, despite extensive research on two-dimensional (2D) ZPR and negative Poisson’s ratio (NPR) metamaterials in recent years. These materials exhibit extraordinarily mechanical properties that provide substantial advantages and broader application prospects over traditional materials. Therefore, in this paper, an orthogonal 3D orthogonal U-shaped accordion (3D-OUA) metamaterial was proposed on the basis of the 2D classical re-entrant hexagon core (RHC) metamaterial. It was applied to the modeling of honeycomb structures, sandwich panels and energy absorbing boxes, and the dynamic response under mechanical compression and dynamic impact was investigated. Finite element models (FEMs) were developed using ABAQUS. Specimens were fabricated using a metal 3D printer. Quasi-static compression experiments and low-velocity impact experiments were conducted, and the dynamic Poisson’s ratio of the honeycomb structure was filmed by a high-speed camera and image processing techniques. Decreasing the amplitude of the cosine function of the unit cell, increasing the tensile thickness of the unit cell and employing multidirectional unit cells to construct the 3D-OUA honeycomb structure can enhance the energy absorption and plateau stress. Reducing the cosine function amplitude of the unit cell, increasing the tensile thickness of the unit cell and employing multidirectional unit cells can enhance the impact resistance of 3D-OUA sandwich panels. In addition, 3D-OUA has obvious advantages over the classical energy absorbing box filling structure and 3D-RHC structure is expected to be applied to the design of new energy absorbing boxes.</div></div>","PeriodicalId":49846,"journal":{"name":"Mechanics Research Communications","volume":"152 ","pages":"Article 104632"},"PeriodicalIF":2.3,"publicationDate":"2026-01-19","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"146038679","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2026-01-19DOI: 10.1016/j.mechrescom.2026.104630
Luca Placidi , Francesco Fabbrocino , Daniel George , Yeaseul Lee , Michele Tepedino
Orthodontic treatments apply controlled displacement or mechanical forces (quasi-static but persistent) to teeth via dental braces. These forces are transmitted through the periodontal ligament (PDL) to the surrounding alveolar bone, initiating a bone remodeling process that allows the teeth to gradually shift position. The PDL transmits load to the alveolar bone and initial stress peaks relax with time as bone resorbs in regions of compression and deposits in regions of tension. Bone tissue, although mostly elastic on short time scales, exhibit therefore adaptive behaviors over longer durations due to cellular remodeling. In this paper we use viscoelastic Maxwell models to capture such an adaptive behavior. Two models are considered. The first is a two-degrees-of-freedom (2DOF) model. The second investigates the continuum two-dimensional case. We derive them, first time in this context, with the use of a Hamilton–Rayleigh principle. For both cases, we show that the load free configuration changes according to the experienced load and the models quantify these concept. The models are general and can be applied to any load history. Besides, the second 2D model can also be applied to any geometric configuration. The examples are nevertheless sufficiently simple to be solved analytically. The target is to help to predict optimal loading schedules to avoid excessive stress or damage.
{"title":"Application of Maxwell models for bone remodeling under orthodontic loading","authors":"Luca Placidi , Francesco Fabbrocino , Daniel George , Yeaseul Lee , Michele Tepedino","doi":"10.1016/j.mechrescom.2026.104630","DOIUrl":"10.1016/j.mechrescom.2026.104630","url":null,"abstract":"<div><div>Orthodontic treatments apply controlled displacement or mechanical forces (quasi-static but persistent) to teeth via dental braces. These forces are transmitted through the periodontal ligament (PDL) to the surrounding alveolar bone, initiating a bone remodeling process that allows the teeth to gradually shift position. The PDL transmits load to the alveolar bone and initial stress peaks relax with time as bone resorbs in regions of compression and deposits in regions of tension. Bone tissue, although mostly elastic on short time scales, exhibit therefore adaptive behaviors over longer durations due to cellular remodeling. In this paper we use viscoelastic Maxwell models to capture such an adaptive behavior. Two models are considered. The first is a two-degrees-of-freedom (2DOF) model. The second investigates the continuum two-dimensional case. We derive them, first time in this context, with the use of a Hamilton–Rayleigh principle. For both cases, we show that the load free configuration changes according to the experienced load and the models quantify these concept. The models are general and can be applied to any load history. Besides, the second 2D model can also be applied to any geometric configuration. The examples are nevertheless sufficiently simple to be solved analytically. The target is to help to predict optimal loading schedules to avoid excessive stress or damage.</div></div>","PeriodicalId":49846,"journal":{"name":"Mechanics Research Communications","volume":"152 ","pages":"Article 104630"},"PeriodicalIF":2.3,"publicationDate":"2026-01-19","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"146038678","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2026-01-16DOI: 10.1016/j.mechrescom.2026.104629
Like Pan , Yin Zhang , Tong Xing , Yuan Yuan , Han Wu
In the differential governing equation of the free vibration of an Euler-Bernoulli beam, the presence of multiple concentrated masses is modeled with the Dirac delta function. A new solution to the beam free vibration, which is based on the generalized functions and their distributional derivative rule, is provided. In the new solution, the shear force discontinuities due to the inertial forces cause by the concentrated masses are incorporated. As the result, the eigenvalue problem formulation with this new solution only needs to enforce the beam boundary conditions. In comparison with the determinant of a matrix (: the number of concentrated masses) formulated by the analytical method, the eigenvalue problem formulated by the new solution form is the determinant of a matrix irrespective of the number of concentrated masses, which significantly reduces the computation. In comparison with the Laplace transform and the Galerkin method, the procedure of formulating an eigenvalue problem is significantly simplified with this new solution. Furthermore, compared with the piece-wise solution form of the analytical method, the one-piece form of this new solution is convenient for enforcing the beam boundary conditions and finding the eigenvectors: The new solution form is with a recursive relation, from which the concise and analytical forms of the characteristic equations are derived and the eigenvectors are readily obtained.
{"title":"Generalized function based solution to free vibration of Euler–Bernoulli beam with arbitrary number of concentrated masses","authors":"Like Pan , Yin Zhang , Tong Xing , Yuan Yuan , Han Wu","doi":"10.1016/j.mechrescom.2026.104629","DOIUrl":"10.1016/j.mechrescom.2026.104629","url":null,"abstract":"<div><div>In the differential governing equation of the free vibration of an Euler-Bernoulli beam, the presence of multiple concentrated masses is modeled with the Dirac delta function. A new solution to the beam free vibration, which is based on the generalized functions and their distributional derivative rule, is provided. In the new solution, the shear force discontinuities due to the inertial forces cause by the concentrated masses are incorporated. As the result, the eigenvalue problem formulation with this new solution only needs to enforce the beam boundary conditions. In comparison with the determinant of a <span><math><mrow><mn>4</mn><mrow><mo>(</mo><mi>N</mi><mo>+</mo><mn>1</mn><mo>)</mo></mrow><mo>×</mo><mn>4</mn><mrow><mo>(</mo><mi>N</mi><mo>+</mo><mn>1</mn><mo>)</mo></mrow></mrow></math></span> matrix (<span><math><mi>N</mi></math></span>: the number of concentrated masses) formulated by the analytical method, the eigenvalue problem formulated by the new solution form is the determinant of a <span><math><mrow><mn>4</mn><mo>×</mo><mn>4</mn></mrow></math></span> matrix irrespective of the number of concentrated masses, which significantly reduces the computation. In comparison with the Laplace transform and the Galerkin method, the procedure of formulating an eigenvalue problem is significantly simplified with this new solution. Furthermore, compared with the piece-wise solution form of the analytical method, the one-piece form of this new solution is convenient for enforcing the beam boundary conditions and finding the eigenvectors: The new solution form is with a recursive relation, from which the concise and analytical forms of the characteristic equations are derived and the eigenvectors are readily obtained.</div></div>","PeriodicalId":49846,"journal":{"name":"Mechanics Research Communications","volume":"152 ","pages":"Article 104629"},"PeriodicalIF":2.3,"publicationDate":"2026-01-16","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145981031","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2026-01-14DOI: 10.1016/j.mechrescom.2026.104626
Ani Luo , Yaming Feng , Heping Liu , Ziying Cao , Guangzhen Xing , Jinxin Lu
This study presents a methodology for analyzing ground contact forces in tensegrity structures under unilateral constraints. A force-law-based contact formulation is first introduced for tensegrity systems. The results indicate that the large penetration depths permitted by this formulation lead to pronounced oscillations in both contact velocity and contact force at the interacting nodes. To overcome this limitation, the normal and tangential contact conditions are reformulated as standard linear complementarity problems. Normal contact is described using a Newtonian impact law to account for kinetic energy dissipation, while tangential contact is governed by Coulomb friction to model the stick–slip transition. For three-dimensional cases, the friction cone is approximated using a polyhedral representation. The resulting dynamic system is solved using a hybrid time-stepping scheme that combines the Moreau midpoint method with a fourth-order Runge–Kutta integrator. Finally, two numerical examples are presented to demonstrate the advantages of the proposed method in terms of accuracy and stability.
{"title":"Study on the nodal contact forces in tensegrity structures","authors":"Ani Luo , Yaming Feng , Heping Liu , Ziying Cao , Guangzhen Xing , Jinxin Lu","doi":"10.1016/j.mechrescom.2026.104626","DOIUrl":"10.1016/j.mechrescom.2026.104626","url":null,"abstract":"<div><div>This study presents a methodology for analyzing ground contact forces in tensegrity structures under unilateral constraints. A force-law-based contact formulation is first introduced for tensegrity systems. The results indicate that the large penetration depths permitted by this formulation lead to pronounced oscillations in both contact velocity and contact force at the interacting nodes. To overcome this limitation, the normal and tangential contact conditions are reformulated as standard linear complementarity problems. Normal contact is described using a Newtonian impact law to account for kinetic energy dissipation, while tangential contact is governed by Coulomb friction to model the stick–slip transition. For three-dimensional cases, the friction cone is approximated using a polyhedral representation. The resulting dynamic system is solved using a hybrid time-stepping scheme that combines the Moreau midpoint method with a fourth-order Runge–Kutta integrator. Finally, two numerical examples are presented to demonstrate the advantages of the proposed method in terms of accuracy and stability.</div></div>","PeriodicalId":49846,"journal":{"name":"Mechanics Research Communications","volume":"152 ","pages":"Article 104626"},"PeriodicalIF":2.3,"publicationDate":"2026-01-14","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"146038698","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2026-01-13DOI: 10.1016/j.mechrescom.2026.104628
Sei-ichiro SAKATA , George STEFANOU , Kazuto SHIRAHAMA , Seiji ONO
This paper is focused on the strength analysis of composite materials and structures considering the heterogeneity of their microstructure. The apparent stiffness and strength properties of composites are significantly affected by the microscopic heterogeneity and randomness, and this influence on the apparent properties of composite structures should be evaluated. In particular, the geometrical randomness in microstructure has a significant influence on the microscopic stress and fracture state, and thus should be appropriately taken into account in the strength estimation of composite structures.
In this work, a multiscale numerical analysis considering microscopic heterogeneity, randomness, and their spatial distribution is discussed. The moving window method-based approach is employed for random field modeling, and finite element-based microscopic fracture analysis is performed for the local apparent strength estimation of composites considering the random microstructure morphologies. The identified probabilistic properties are introduced into the multiscale probabilistic strength analysis of composite structures, and in addition, the influence of correlation between the local apparent elastic property and strength on the estimated apparent strength of composite structures is investigated.
{"title":"Probabilistic strength estimation analysis of composites considering cross-correlated random fields of local strength and apparent elastic properties","authors":"Sei-ichiro SAKATA , George STEFANOU , Kazuto SHIRAHAMA , Seiji ONO","doi":"10.1016/j.mechrescom.2026.104628","DOIUrl":"10.1016/j.mechrescom.2026.104628","url":null,"abstract":"<div><div>This paper is focused on the strength analysis of composite materials and structures considering the heterogeneity of their microstructure. The apparent stiffness and strength properties of composites are significantly affected by the microscopic heterogeneity and randomness, and this influence on the apparent properties of composite structures should be evaluated. In particular, the geometrical randomness in microstructure has a significant influence on the microscopic stress and fracture state, and thus should be appropriately taken into account in the strength estimation of composite structures.</div><div>In this work, a multiscale numerical analysis considering microscopic heterogeneity, randomness, and their spatial distribution is discussed. The moving window method-based approach is employed for random field modeling, and finite element-based microscopic fracture analysis is performed for the local apparent strength estimation of composites considering the random microstructure morphologies. The identified probabilistic properties are introduced into the multiscale probabilistic strength analysis of composite structures, and in addition, the influence of correlation between the local apparent elastic property and strength on the estimated apparent strength of composite structures is investigated.</div></div>","PeriodicalId":49846,"journal":{"name":"Mechanics Research Communications","volume":"152 ","pages":"Article 104628"},"PeriodicalIF":2.3,"publicationDate":"2026-01-13","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"146038677","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2026-01-13DOI: 10.1016/j.mechrescom.2026.104627
R. Rodríguez-Ramos , Y. Espinosa-Almeyda , J.A. Otero , H. Camacho-Montes , P. Rodríguez-Bermúdez , B.J. Mederos
In this work, laminated and fiber-reinforced elastic micropolar multiphase composites are analyzed. The aim is to evaluate the effective stiffness and torque properties, assuming isotropic centrosymmetric constituents and ideal interfacial bonding. From the asymptotic homogenization method (AHM), the corresponding formulations of the local cell problems and the effective properties are derived. The local problems are solved using the methodology reported for laminated and fibrous media in previous works. Herein, for the first time, a direct and systematic quantitative comparison between the effective properties of laminated and fibrous micropolar media is presented. In addition, an extensive parametric analysis is carried out for some real material systems (SyF, PUF, PMIF-WF51, PMIF-WF110), generating results useful for engineering applications, design, and model validation. Moreover, an integration of both microstructural architectures within a micropolar framework is shown, enabling a consistent comparison of how geometry influences the effective responses. Numerical results illustrate the findings for different two-phase composites. They reveal the influence of microstructural geometry on the effective behaviors, highlighting significant differences between laminated and fibrous composites. These findings could be useful for applications and offer a solid reference framework for validating both numerical and experimental studies.
{"title":"Multiscale analysis of periodic micropolar elastic media","authors":"R. Rodríguez-Ramos , Y. Espinosa-Almeyda , J.A. Otero , H. Camacho-Montes , P. Rodríguez-Bermúdez , B.J. Mederos","doi":"10.1016/j.mechrescom.2026.104627","DOIUrl":"10.1016/j.mechrescom.2026.104627","url":null,"abstract":"<div><div>In this work, laminated and fiber-reinforced elastic micropolar multiphase composites are analyzed. The aim is to evaluate the effective stiffness and torque properties, assuming isotropic centrosymmetric constituents and ideal interfacial bonding. From the asymptotic homogenization method (AHM), the corresponding formulations of the local cell problems and the effective properties are derived. The local problems are solved using the methodology reported for laminated and fibrous media in previous works. Herein, for the first time, a direct and systematic quantitative comparison between the effective properties of laminated and fibrous micropolar media is presented. In addition, an extensive parametric analysis is carried out for some real material systems (SyF, PUF, PMIF-WF51, PMIF-WF110), generating results useful for engineering applications, design, and model validation. Moreover, an integration of both microstructural architectures within a micropolar framework is shown, enabling a consistent comparison of how geometry influences the effective responses. Numerical results illustrate the findings for different two-phase composites. They reveal the influence of microstructural geometry on the effective behaviors, highlighting significant differences between laminated and fibrous composites. These findings could be useful for applications and offer a solid reference framework for validating both numerical and experimental studies.</div></div>","PeriodicalId":49846,"journal":{"name":"Mechanics Research Communications","volume":"152 ","pages":"Article 104627"},"PeriodicalIF":2.3,"publicationDate":"2026-01-13","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145981029","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2026-01-12DOI: 10.1016/j.mechrescom.2026.104625
Leslie D. Pérez-Fernández , Julián Bravo-Castillero
Motivated by the usual occurrence of interlayer defects in materials and structures produced via additive manufacturing, a homogenization methodology is presented here, based on the scale invariance of the strain energy density and asymptotic homogenization, to obtain analytical formulas for the effective coefficients of periodic generalized Hookean composites with spring-type interfacial bonding and strain gradient global behavior in a two-dimensional setting. The first finding, differing from the perfect bonding case, is that the general fourth-rank strain gradient effective coefficients differ from the classical ones. Next, specialization to two-phase laminates with isotropic layers and sandwich-type periodicity cell allowed obtaining analytical formulas for all the strain gradient effective coefficients. Three more findings follow from these formulas: the imperfect bonding causes the appearance of new nonzero effective coefficients, which is in contrast with the perfect bonding case; proportionality relations among effective coefficients of the same rank are revealed, which are useful for simplifying and controlling their implementations and validating other computational or semi-analytical homogenization approaches conceived for more complex structures; and the higher-order effective coefficients vanish for the pseudo-homogeneous material with periodicity cell made of two imperfectly bonded layers of the same constitutive material, so having internal structure of a composite material is insufficient for producing strain gradient effective behavior. Numerical examples, carried out with experimental material and spring data, illustrate the theoretical findings. To the best of our knowledge, the results obtained her e are new and generalizable in a straightforward manner to the three-dimensional setting with any number of layers and anisotropy.
{"title":"Fully analytical calculation of the strain gradient effective coefficients of periodic bilayered metamaterials with spring-type interlayer defects","authors":"Leslie D. Pérez-Fernández , Julián Bravo-Castillero","doi":"10.1016/j.mechrescom.2026.104625","DOIUrl":"10.1016/j.mechrescom.2026.104625","url":null,"abstract":"<div><div>Motivated by the usual occurrence of interlayer defects in materials and structures produced via additive manufacturing, a homogenization methodology is presented here, based on the scale invariance of the strain energy density and asymptotic homogenization, to obtain analytical formulas for the effective coefficients of periodic generalized Hookean composites with spring-type interfacial bonding and strain gradient global behavior in a two-dimensional setting. The first finding, differing from the perfect bonding case, is that the general fourth-rank strain gradient effective coefficients differ from the classical ones. Next, specialization to two-phase laminates with isotropic layers and sandwich-type periodicity cell allowed obtaining analytical formulas for all the strain gradient effective coefficients. Three more findings follow from these formulas: the imperfect bonding causes the appearance of new nonzero effective coefficients, which is in contrast with the perfect bonding case; proportionality relations among effective coefficients of the same rank are revealed, which are useful for simplifying and controlling their implementations and validating other computational or semi-analytical homogenization approaches conceived for more complex structures; and the higher-order effective coefficients vanish for the pseudo-homogeneous material with periodicity cell made of two imperfectly bonded layers of the same constitutive material, so having internal structure of a composite material is insufficient for producing strain gradient effective behavior. Numerical examples, carried out with experimental material and spring data, illustrate the theoretical findings. To the best of our knowledge, the results obtained her e are new and generalizable in a straightforward manner to the three-dimensional setting with any number of layers and anisotropy.</div></div>","PeriodicalId":49846,"journal":{"name":"Mechanics Research Communications","volume":"152 ","pages":"Article 104625"},"PeriodicalIF":2.3,"publicationDate":"2026-01-12","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145981028","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2026-01-09DOI: 10.1016/j.mechrescom.2026.104621
Jimei Wu , Rui Li , Mingyue Shao , Chaoyue Lin , Qiumin Wu
In actual printing production, the temperature distribution within the hot-air drying environment of the oven is non-uniform. Therefore, this paper investigates the nonlinear vibration of a fractional-order viscoelastic PET film under a non-uniform temperature field. The temperature distribution within the hot-air drying environment of the oven in printing equipment is assumed to be non-uniform. The thermal viscoelasticity of the fractional-order PET film is described based on thermal strain and the fractional-order viscoelastic Kelvin-Voigt constitutive model. Based on von Kármán theory and D'Alembert's principle, establish a fractional-order nonlinear vibrational differential equation for the motion of PET film under non-uniform temperature fields during the roll-to-roll printing process.The nonlinear vibration equation is discretized using the Galerkin method, combined with finite difference methods and discrete Caputo fractional derivatives, to obtain the nonlinear vibration response of the fractional-order viscoelastic PET film under non-uniform temperature effects. Programming in MATLAB generates bifurcation diagrams, phase trajectories, time history plots, phase diagrams, and Poincaré sections. The study analyzes the influence of parameters such as fractional order, non-uniform hot air temperature, and external hot air excitation on the nonlinear vibration of the film. This study establishes a film vibration model coupling non-uniform temperature fields with fractional viscoelasticity, quantifies the control mechanisms of multiple parameters on nonlinear vibration, and provides a theoretical basis for vibration control and parameter optimization of PET films during printing processes.
{"title":"Nonlinear vibration analysis of fractional-order viscoelastic pet membrane under non-uniformly distributed temperature field","authors":"Jimei Wu , Rui Li , Mingyue Shao , Chaoyue Lin , Qiumin Wu","doi":"10.1016/j.mechrescom.2026.104621","DOIUrl":"10.1016/j.mechrescom.2026.104621","url":null,"abstract":"<div><div>In actual printing production, the temperature distribution within the hot-air drying environment of the oven is non-uniform. Therefore, this paper investigates the nonlinear vibration of a fractional-order viscoelastic PET film under a non-uniform temperature field. The temperature distribution within the hot-air drying environment of the oven in printing equipment is assumed to be non-uniform. The thermal viscoelasticity of the fractional-order PET film is described based on thermal strain and the fractional-order viscoelastic Kelvin-Voigt constitutive model. Based on von Kármán theory and D'Alembert's principle, establish a fractional-order nonlinear vibrational differential equation for the motion of PET film under non-uniform temperature fields during the roll-to-roll printing process.The nonlinear vibration equation is discretized using the Galerkin method, combined with finite difference methods and discrete Caputo fractional derivatives, to obtain the nonlinear vibration response of the fractional-order viscoelastic PET film under non-uniform temperature effects. Programming in MATLAB generates bifurcation diagrams, phase trajectories, time history plots, phase diagrams, and Poincaré sections. The study analyzes the influence of parameters such as fractional order, non-uniform hot air temperature, and external hot air excitation on the nonlinear vibration of the film. This study establishes a film vibration model coupling non-uniform temperature fields with fractional viscoelasticity, quantifies the control mechanisms of multiple parameters on nonlinear vibration, and provides a theoretical basis for vibration control and parameter optimization of PET films during printing processes.</div></div>","PeriodicalId":49846,"journal":{"name":"Mechanics Research Communications","volume":"152 ","pages":"Article 104621"},"PeriodicalIF":2.3,"publicationDate":"2026-01-09","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"146038699","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}