首页 > 最新文献

Archive of Applied Mechanics最新文献

英文 中文
Reducing the dynamic response of self-supporting lattice towers with optimized connection dampers 优化连接阻尼器降低自支撑格构塔动力响应
IF 2.5 3区 工程技术 Q2 MECHANICS Pub Date : 2025-11-09 DOI: 10.1007/s00419-025-02981-7
Luiz Guilherme Grotto, Letícia Fleck Fadel Miguel, João Kaminski Junior

The design of metallic support towers for energy transmission and telecommunication networks increasingly demands efficiency and reliability. As these towers become more slender, they are highly susceptible to wind-induced vibrations, which can compromise structural stability. Conventional passive dampers are effective in controlling vibrations but often increase project costs, limiting practical implementation. This study investigates an alternative approach using connection dampers based on rubber rings installed in bolted joints. These dampers reduce connection stiffness while increasing damping, enhancing the structure’s energy dissipation capacity. Experimental tests characterized the mechanical properties of the rubber rings, and we implemented the results in a numerical optimization framework applied to a three-dimensional lattice steel tower under synoptic wind. The Whale Optimization Algorithm determined optimal stiffness and damping parameters. Comparisons between rigid and semi-rigid connections equipped with the proposed dampers demonstrate notable improvements: the structural damping ratio increased by approximately 173%, and the maximum displacement at the tower top decreased by around 28%. These results confirm the effectiveness of the proposed methodology in mitigating wind-induced vibrations and improving the dynamic response of slender lattice towers.

能源传输和电信网络金属支撑塔的设计对效率和可靠性的要求越来越高。随着这些塔变得越来越细长,它们非常容易受到风引起的振动的影响,这可能会损害结构的稳定性。传统的被动阻尼器在控制振动方面是有效的,但往往会增加工程成本,限制了实际应用。本研究探讨了另一种方法,即在螺栓连接中安装基于橡胶环的连接阻尼器。这些阻尼器在减小连接刚度的同时增加了阻尼,增强了结构的耗能能力。实验测试表征了橡胶环的力学性能,并将结果应用于天气风作用下的三维晶格钢塔的数值优化框架中。鲸鱼优化算法确定了最优的刚度和阻尼参数。通过刚性连接与半刚性连接的比较,发现了显著的改善:结构阻尼比增加了约173%,塔顶最大位移减少了约28%。这些结果证实了所提出的方法在减轻风致振动和改善细长晶格塔的动力响应方面的有效性。
{"title":"Reducing the dynamic response of self-supporting lattice towers with optimized connection dampers","authors":"Luiz Guilherme Grotto,&nbsp;Letícia Fleck Fadel Miguel,&nbsp;João Kaminski Junior","doi":"10.1007/s00419-025-02981-7","DOIUrl":"10.1007/s00419-025-02981-7","url":null,"abstract":"<div><p>The design of metallic support towers for energy transmission and telecommunication networks increasingly demands efficiency and reliability. As these towers become more slender, they are highly susceptible to wind-induced vibrations, which can compromise structural stability. Conventional passive dampers are effective in controlling vibrations but often increase project costs, limiting practical implementation. This study investigates an alternative approach using connection dampers based on rubber rings installed in bolted joints. These dampers reduce connection stiffness while increasing damping, enhancing the structure’s energy dissipation capacity. Experimental tests characterized the mechanical properties of the rubber rings, and we implemented the results in a numerical optimization framework applied to a three-dimensional lattice steel tower under synoptic wind. The Whale Optimization Algorithm determined optimal stiffness and damping parameters. Comparisons between rigid and semi-rigid connections equipped with the proposed dampers demonstrate notable improvements: the structural damping ratio increased by approximately 173%, and the maximum displacement at the tower top decreased by around 28%. These results confirm the effectiveness of the proposed methodology in mitigating wind-induced vibrations and improving the dynamic response of slender lattice towers.</p></div>","PeriodicalId":477,"journal":{"name":"Archive of Applied Mechanics","volume":"95 12","pages":""},"PeriodicalIF":2.5,"publicationDate":"2025-11-09","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145469249","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
引用次数: 0
Structural responses of functionally graded graphene nanoplatelets-reinforced composite plates using inverse hyperbolic shear deformation theory 利用反双曲剪切变形理论研究功能梯度石墨烯纳米片增强复合材料板的结构响应
IF 2.5 3区 工程技术 Q2 MECHANICS Pub Date : 2025-11-09 DOI: 10.1007/s00419-025-02971-9
Sumit Kumar, Rosalin Sahoo, Rachit Panda

In this work, the inverse hyperbolic shear deformation theory (IHSDT) is used to examine the bending, buckling, and free vibration responses of functionally graded graphene nanoplatelets-reinforced composite (FG-GNPRC) plates. This theory ensures that there are no traction forces on both the top and bottom sides of the plate. It achieves this without needing to use a shear correction factor, and it allows for a nonlinear distribution of transverse shear stresses across the plate. The study utilizes a finite element approach incorporating a nonlinear function based on inverse hyperbolic sine function. A C0 finite element model is created for the GNPRC plate in the framework of IHSDT to determine the structural responses of the plate in MATLAB environment. It studies graphene nanoplatelets (GPLs) distributions, patterned as UD, FG-X, FG-O, and FG-A throughout the thickness in composite plate. The weight fraction (wt%) of GPLs varies along the thickness direction and is evenly distributed throughout the matrix of each layer, follows a specific distribution pattern. The Halpin–Tsai micromodel is used to estimate the effective Young’s modulus of the GNPRC plate, and the rule of mixtures is used to calculate the Poisson’s ratio and mass density. The plate domain is discretized using an eight-noded elements, each with 56 degrees of freedom. Further the analysis looks at the effects of a variety of factors, including the number of layers (NL) of GNPRC plate, length/thickness, and length/width ratio of GPLs, wt% of GPLs, and dispersion patterns of GPLs on the structural responses of FG-GNPRC plate. The numerical results demonstrate that the rigidity of plates can be significantly enhances by incorporating GPLs and the outcomes were compared with prior findings in order to evaluate the performance and effectiveness of the suggested mathematical approach.

在这项工作中,利用逆双曲剪切变形理论(IHSDT)来研究功能梯度石墨烯纳米片增强复合材料(FG-GNPRC)板的弯曲、屈曲和自由振动响应。这个理论保证了在板的上下两侧都没有牵引力。它不需要使用剪切校正系数就可以实现这一点,并且它允许横向剪切应力在板上的非线性分布。该研究采用了基于双曲正弦反函数的非线性函数的有限元方法。建立了在IHSDT框架下GNPRC板的C0有限元模型,在MATLAB环境下确定了该板的结构响应。它研究了石墨烯纳米片(gpl)在复合材料板厚度上的分布,图案为UD、FG-X、FG-O和FG-A。gpl的重量分数(wt%)沿厚度方向变化,并均匀分布在每层基质中,遵循特定的分布模式。采用Halpin-Tsai微观模型估计GNPRC板的有效杨氏模量,采用混合规律计算泊松比和质量密度。板域采用8结点单元离散化,每个单元56个自由度。进一步分析了多种因素对FG-GNPRC板结构响应的影响,包括GNPRC板的层数(NL)、gpl的长度/厚度、长宽比、gpl的wt%和gpl的色散模式。数值结果表明,加入gpl可以显著提高板的刚度,并将结果与先前的研究结果进行了比较,以评估所建议的数学方法的性能和有效性。
{"title":"Structural responses of functionally graded graphene nanoplatelets-reinforced composite plates using inverse hyperbolic shear deformation theory","authors":"Sumit Kumar,&nbsp;Rosalin Sahoo,&nbsp;Rachit Panda","doi":"10.1007/s00419-025-02971-9","DOIUrl":"10.1007/s00419-025-02971-9","url":null,"abstract":"<div><p>In this work, the inverse hyperbolic shear deformation theory (IHSDT) is used to examine the bending, buckling, and free vibration responses of functionally graded graphene nanoplatelets-reinforced composite (FG-GNPRC) plates. This theory ensures that there are no traction forces on both the top and bottom sides of the plate. It achieves this without needing to use a shear correction factor, and it allows for a nonlinear distribution of transverse shear stresses across the plate. The study utilizes a finite element approach incorporating a nonlinear function based on inverse hyperbolic sine function. A C<sup>0</sup> finite element model is created for the GNPRC plate in the framework of IHSDT to determine the structural responses of the plate in MATLAB environment. It studies graphene nanoplatelets (GPLs) distributions, patterned as UD, FG-X, FG-O, and FG-A throughout the thickness in composite plate. The weight fraction (wt%) of GPLs varies along the thickness direction and is evenly distributed throughout the matrix of each layer, follows a specific distribution pattern. The Halpin–Tsai micromodel is used to estimate the effective Young’s modulus of the GNPRC plate, and the rule of mixtures is used to calculate the Poisson’s ratio and mass density. The plate domain is discretized using an eight-noded elements, each with 56 degrees of freedom. Further the analysis looks at the effects of a variety of factors, including the number of layers (<i>N</i><sub><i>L</i></sub>) of GNPRC plate, length/thickness, and length/width ratio of GPLs, wt% of GPLs, and dispersion patterns of GPLs on the structural responses of FG-GNPRC plate. The numerical results demonstrate that the rigidity of plates can be significantly enhances by incorporating GPLs and the outcomes were compared with prior findings in order to evaluate the performance and effectiveness of the suggested mathematical approach.</p></div>","PeriodicalId":477,"journal":{"name":"Archive of Applied Mechanics","volume":"95 12","pages":""},"PeriodicalIF":2.5,"publicationDate":"2025-11-09","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145469295","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
引用次数: 0
Hybrid finite-discrete element simulation of the bending fracture process in sandstone with prefabricated cracks under different loading conditions 不同加载条件下预制裂缝砂岩弯曲断裂过程的有限-离散元混合模拟
IF 2.5 3区 工程技术 Q2 MECHANICS Pub Date : 2025-11-09 DOI: 10.1007/s00419-025-02978-2
Huaming An, Bingbing Yang, Hongyuan Liu, Xinghai Mu, Xin Zhang

Bending fracture is a critical failure mechanism in many engineering structures such as rock beams, slabs and bridge components. This is particularly the case when prefabricated cracks are present. In order to accurately simulate such failure processes, a method is required that can capture both continuous deformation and the transition to discontinuous fracturing. The hybrid finite-discrete element method (HFDEM) is particularly well-suited for this purpose, as it integrates the advantages of continuum-based finite element methods with those of discontinuum-based discrete element methods. However, conventional HFDEM approaches are computationally intensive and often impractical for large-scale problems. This paper introduces a general-purpose graphic processing unit-parallelized HFDEM that is general-purpose and achieves execution speeds exceeding 128 times those of sequential code. The model is employed to simulate the bending fracture processes of sandstone specimens with prefabricated cracks under varying test configurations, revealing the mechanism by which loading rate influences fracture mechanics behavior. The results demonstrate good agreement with existing research findings, thereby validating the model. The core contribution, however, lies in revealing the intrinsic mechanisms behind the strongly rate-dependent fracture behavior. The simulations elucidate how higher loading rates weaken the guiding effect of prefabricated cracks and promote a transition to complex mixed-mode fracture, thereby establishing that the apparent toughness increase is primarily driven by inertial effects and exhibits fundamental divergence between different fracture modes.

在岩石梁、板、桥梁构件等工程结构中,弯曲断裂是一种重要的破坏机制。当预制裂缝存在时,情况尤其如此。为了准确地模拟这种破坏过程,需要一种既能捕捉连续变形又能捕捉到向不连续压裂过渡的方法。混合有限-离散单元法(HFDEM)特别适合于这一目的,因为它结合了基于连续单元的有限元方法和基于非连续单元的离散单元方法的优点。然而,传统的HFDEM方法计算量大,对于大规模问题往往不切实际。本文介绍了一种通用图形处理单元——并行化HFDEM,它具有通用性,执行速度超过顺序码的128倍。利用该模型模拟了含预制裂缝砂岩试件在不同试验配置下的弯曲断裂过程,揭示了加载速率影响断裂力学行为的机理。结果与已有的研究结果吻合较好,从而验证了模型的有效性。然而,核心的贡献在于揭示了强烈依赖速率的破裂行为背后的内在机制。模拟结果表明,较高的加载速率削弱了预制裂纹的导向作用,促进了预制裂纹向复杂混合模式断裂的过渡,从而证明了表观韧性的增加主要是由惯性效应驱动的,不同断裂模式之间表现出根本的差异。
{"title":"Hybrid finite-discrete element simulation of the bending fracture process in sandstone with prefabricated cracks under different loading conditions","authors":"Huaming An,&nbsp;Bingbing Yang,&nbsp;Hongyuan Liu,&nbsp;Xinghai Mu,&nbsp;Xin Zhang","doi":"10.1007/s00419-025-02978-2","DOIUrl":"10.1007/s00419-025-02978-2","url":null,"abstract":"<div><p>Bending fracture is a critical failure mechanism in many engineering structures such as rock beams, slabs and bridge components. This is particularly the case when prefabricated cracks are present. In order to accurately simulate such failure processes, a method is required that can capture both continuous deformation and the transition to discontinuous fracturing. The hybrid finite-discrete element method (HFDEM) is particularly well-suited for this purpose, as it integrates the advantages of continuum-based finite element methods with those of discontinuum-based discrete element methods. However, conventional HFDEM approaches are computationally intensive and often impractical for large-scale problems. This paper introduces a general-purpose graphic processing unit-parallelized HFDEM that is general-purpose and achieves execution speeds exceeding 128 times those of sequential code. The model is employed to simulate the bending fracture processes of sandstone specimens with prefabricated cracks under varying test configurations, revealing the mechanism by which loading rate influences fracture mechanics behavior. The results demonstrate good agreement with existing research findings, thereby validating the model. The core contribution, however, lies in revealing the intrinsic mechanisms behind the strongly rate-dependent fracture behavior. The simulations elucidate how higher loading rates weaken the guiding effect of prefabricated cracks and promote a transition to complex mixed-mode fracture, thereby establishing that the apparent toughness increase is primarily driven by inertial effects and exhibits fundamental divergence between different fracture modes.</p></div>","PeriodicalId":477,"journal":{"name":"Archive of Applied Mechanics","volume":"95 12","pages":""},"PeriodicalIF":2.5,"publicationDate":"2025-11-09","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145469248","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
引用次数: 0
Photothermal response and wave propagation in size-dependent functionally graded semiconductor materials: a nonlocal modified thermoelastic analysis 尺寸相关功能梯度半导体材料中的光热响应和波传播:非局部修正热弹性分析
IF 2.5 3区 工程技术 Q2 MECHANICS Pub Date : 2025-11-09 DOI: 10.1007/s00419-025-02986-2
Abeer Alhashash, Ahmed E. Abouelregal

This study presents a novel analytical framework to investigate the photothermal response of functionally graded semiconductor materials under laser pulse excitation, utilizing a nonlocal thermoelastic model. By integrating the single-phase-lagging (SPL) concept with a nonlocal thermal length-scale, the model captures finite thermal wave speeds critical for micro- and nanoscale transient phenomena. The analysis examines coupled thermal, elastic, and plasma wave interactions in materials with spatially varying properties. Numerical results show that increasing the nonlocal thermal length-scale significantly reduces temperature and carrier density, while material gradation enhances near-surface temperature but markedly suppresses displacement. Longer laser pulse durations lower peak temperatures but increase stresses, highlighting key trade-offs. By overcoming limitations of classical Fourier heat conduction and incorporating size-dependent effects, the model provides insights into thermoelastic wave propagation in nanostructured materials. These findings emphasize the importance of material design and pulse optimization for improving energy transfer and stress management, with applications in semiconductor device engineering, thermal management, and laser-based manufacturing in nanotechnology and optoelectronics.

本研究提出了一种新的分析框架,利用非局部热弹性模型来研究激光脉冲激励下功能梯度半导体材料的光热响应。通过将单相滞后(SPL)概念与非局部热长度尺度相结合,该模型捕获了对微纳米瞬态现象至关重要的有限热波速度。分析检查耦合热,弹性和等离子体波相互作用的材料与空间变化的性质。数值结果表明,增加非局部热长度尺度可以显著降低温度和载流子密度,而材料级配可以提高近表面温度,但显著抑制位移。较长的激光脉冲持续时间降低峰值温度,但增加应力,突出关键的权衡。通过克服经典傅立叶热传导的局限性,并结合尺寸依赖效应,该模型为纳米结构材料中的热弹性波传播提供了见解。这些发现强调了材料设计和脉冲优化对于改善能量传递和应力管理的重要性,并在半导体器件工程、热管理以及纳米技术和光电子学中的激光制造中得到了应用。
{"title":"Photothermal response and wave propagation in size-dependent functionally graded semiconductor materials: a nonlocal modified thermoelastic analysis","authors":"Abeer Alhashash,&nbsp;Ahmed E. Abouelregal","doi":"10.1007/s00419-025-02986-2","DOIUrl":"10.1007/s00419-025-02986-2","url":null,"abstract":"<div><p>This study presents a novel analytical framework to investigate the photothermal response of functionally graded semiconductor materials under laser pulse excitation, utilizing a nonlocal thermoelastic model. By integrating the single-phase-lagging (SPL) concept with a nonlocal thermal length-scale, the model captures finite thermal wave speeds critical for micro- and nanoscale transient phenomena. The analysis examines coupled thermal, elastic, and plasma wave interactions in materials with spatially varying properties. Numerical results show that increasing the nonlocal thermal length-scale significantly reduces temperature and carrier density, while material gradation enhances near-surface temperature but markedly suppresses displacement. Longer laser pulse durations lower peak temperatures but increase stresses, highlighting key trade-offs. By overcoming limitations of classical Fourier heat conduction and incorporating size-dependent effects, the model provides insights into thermoelastic wave propagation in nanostructured materials. These findings emphasize the importance of material design and pulse optimization for improving energy transfer and stress management, with applications in semiconductor device engineering, thermal management, and laser-based manufacturing in nanotechnology and optoelectronics.</p></div>","PeriodicalId":477,"journal":{"name":"Archive of Applied Mechanics","volume":"95 12","pages":""},"PeriodicalIF":2.5,"publicationDate":"2025-11-09","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145469294","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
引用次数: 0
Dynamic response of the tunnel embedded in a saturated frozen half-space with wave functions method 基于波函数法的饱和冻结半空间隧道动力响应分析
IF 2.5 3区 工程技术 Q2 MECHANICS Pub Date : 2025-11-08 DOI: 10.1007/s00419-025-02982-6
Shuocheng Zhang, Menglei Ji, Yingzhe Zhang, Wenhua Chen, Ruichen Zheng

An important influence that distinguishes the dynamic response of tunnels embedded in saturated frozen soil from that of unfrozen soil is the freezing of liquid-phase water in the pores, and the formation of ice alters the propagation of loads through the frozen soil. The dynamic response of the tunnel embedded in saturated frozen half-space when loads are applied to the lining invert is addressed by the wave function method. The saturated frozen soil is regarded as a poroelastic medium, the lining is regarded as a hollow cylinder, and the total wave field in the surrounding medium consists of outgoing cylindrical waves and down-going plane waves. The plane wave function and cylindrical wave function could be converted to obtain the boundary conditions that are convenient for solving the ground surface and soil-lining interface in the rectangular and cylindrical coordinate systems. The parametric analysis demonstrates that the dynamic response of a tunnel embedded in a saturated frozen half-space and the surrounding medium could be influenced by the permeability, porosity, ice content of frozen soil and the change of the tunnel depth. This method provides a novel approach to the safety and stability of tunnels in cold regions and serves as a reference for predicting the vibration of tunnels embedded in saturated frozen half-spaces.

饱和冻土中埋深隧道动力响应与非冻土中埋深隧道动力响应的一个重要影响因素是孔隙中液相水的冻结,冰的形成改变了荷载在冻土中的传播。采用波函数法分析了饱和冻结半空间隧道衬砌仰拱受荷载作用时的动力响应。将饱和冻土视为多孔弹性介质,将衬砌视为空心圆柱体,周围介质中的总波场由向外的柱面波和向内的平面波组成。将平面波函数与柱面波函数进行转换,得到便于在直角坐标系和柱面坐标系下求解地表与衬砌界面的边界条件。参数分析表明,埋置在饱和冻结半空间及周围介质中的隧道动力响应受土体渗透性、孔隙度、冻土含冰量及隧道深度变化的影响。该方法为研究寒冷地区隧道的安全与稳定提供了一种新的途径,并为饱和冻结半空间隧道的振动预测提供了参考。
{"title":"Dynamic response of the tunnel embedded in a saturated frozen half-space with wave functions method","authors":"Shuocheng Zhang,&nbsp;Menglei Ji,&nbsp;Yingzhe Zhang,&nbsp;Wenhua Chen,&nbsp;Ruichen Zheng","doi":"10.1007/s00419-025-02982-6","DOIUrl":"10.1007/s00419-025-02982-6","url":null,"abstract":"<div><p>An important influence that distinguishes the dynamic response of tunnels embedded in saturated frozen soil from that of unfrozen soil is the freezing of liquid-phase water in the pores, and the formation of ice alters the propagation of loads through the frozen soil. The dynamic response of the tunnel embedded in saturated frozen half-space when loads are applied to the lining invert is addressed by the wave function method. The saturated frozen soil is regarded as a poroelastic medium, the lining is regarded as a hollow cylinder, and the total wave field in the surrounding medium consists of outgoing cylindrical waves and down-going plane waves. The plane wave function and cylindrical wave function could be converted to obtain the boundary conditions that are convenient for solving the ground surface and soil-lining interface in the rectangular and cylindrical coordinate systems. The parametric analysis demonstrates that the dynamic response of a tunnel embedded in a saturated frozen half-space and the surrounding medium could be influenced by the permeability, porosity, ice content of frozen soil and the change of the tunnel depth. This method provides a novel approach to the safety and stability of tunnels in cold regions and serves as a reference for predicting the vibration of tunnels embedded in saturated frozen half-spaces.</p></div>","PeriodicalId":477,"journal":{"name":"Archive of Applied Mechanics","volume":"95 11","pages":""},"PeriodicalIF":2.5,"publicationDate":"2025-11-08","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145510288","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
引用次数: 0
Dispersion characteristics of Lamb waves in porous functionally graded plate with free boundaries 自由边界多孔功能梯度板中Lamb波的色散特性
IF 2.5 3区 工程技术 Q2 MECHANICS Pub Date : 2025-11-06 DOI: 10.1007/s00419-025-02984-4
Jiezhi Wang, Jiahe Guo, Zhaochun Teng

This study investigates the dispersion characteristics of Lamb waves in porous functionally graded materials (FGM) with free boundaries, with a focus on the influence of porosity. The material properties of porous FGM are characterized by a mixed power-law model. Based on elastic wave theory, the wave equations for an infinite porous FGM plate are derived, and the dispersion equations of Lamb waves are obtained using the Wentzel–Kramers–Brillouin–Jeffreys (WKBJ) method. To validate the proposed approach, the porous FGM layer is reduced to an isotropic material, and the results are compared with existing literature, demonstrating the method’s accuracy. Numerical solutions of the dispersion equations are implemented via MATLAB to analyze the effects of porosity, layer thickness and gradient index on dispersion curves. The results indicate that porosity, layer thickness and gradient index significantly affect Lamb wave dispersion. This study provides both theoretical and numerical foundations for understanding and utilizing Lamb waves in porous FGM structures, with potential applications in structural health monitoring and nondestructive testing for aerospace and civil engineering.

本文研究了具有自由边界的多孔功能梯度材料(FGM)中Lamb波的频散特性,重点研究了孔隙率的影响。多孔FGM的材料性能采用混合幂律模型表征。基于弹性波理论,推导了无限多孔FGM板的波动方程,并采用WKBJ方法得到了Lamb波的频散方程。为了验证所提出的方法,将多孔FGM层简化为各向同性材料,并将结果与现有文献进行了比较,证明了该方法的准确性。通过MATLAB对色散方程进行数值求解,分析孔隙度、层厚和梯度指数对色散曲线的影响。结果表明,孔隙度、层厚和梯度指数对兰姆波色散有显著影响。该研究为理解和利用多孔FGM结构中的Lamb波提供了理论和数值基础,在航空航天和土木工程的结构健康监测和无损检测方面具有潜在的应用前景。
{"title":"Dispersion characteristics of Lamb waves in porous functionally graded plate with free boundaries","authors":"Jiezhi Wang,&nbsp;Jiahe Guo,&nbsp;Zhaochun Teng","doi":"10.1007/s00419-025-02984-4","DOIUrl":"10.1007/s00419-025-02984-4","url":null,"abstract":"<div><p>This study investigates the dispersion characteristics of Lamb waves in porous functionally graded materials (FGM) with free boundaries, with a focus on the influence of porosity. The material properties of porous FGM are characterized by a mixed power-law model. Based on elastic wave theory, the wave equations for an infinite porous FGM plate are derived, and the dispersion equations of Lamb waves are obtained using the Wentzel–Kramers–Brillouin–Jeffreys (WKBJ) method. To validate the proposed approach, the porous FGM layer is reduced to an isotropic material, and the results are compared with existing literature, demonstrating the method’s accuracy. Numerical solutions of the dispersion equations are implemented via MATLAB to analyze the effects of porosity, layer thickness and gradient index on dispersion curves. The results indicate that porosity, layer thickness and gradient index significantly affect Lamb wave dispersion. This study provides both theoretical and numerical foundations for understanding and utilizing Lamb waves in porous FGM structures, with potential applications in structural health monitoring and nondestructive testing for aerospace and civil engineering.</p></div>","PeriodicalId":477,"journal":{"name":"Archive of Applied Mechanics","volume":"95 11","pages":""},"PeriodicalIF":2.5,"publicationDate":"2025-11-06","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145456448","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
引用次数: 0
The influence of poling direction on an interface crack with contact zones in a piezoelectric bimaterial 极化方向对压电双材料接触面裂纹的影响
IF 2.5 3区 工程技术 Q2 MECHANICS Pub Date : 2025-11-04 DOI: 10.1007/s00419-025-02983-5
Volodymyr Govorukha, Tetyana Kagadiy, Marc Kamlah

The plane strain problem for an interface crack between two different piezoelectric materials with an arbitrary polarization direction under remote mechanical loading is analyzed. Mechanically frictionless and electrically permeable contact zones are assumed at the crack tips, and the remaining part of the crack is considered as electrically permeable. An analytical method is used based upon the presentations of the field variables via sectionally analytic vector functions with the resulting formulation and the solution of the combined Dirichlet–Riemann boundary value problem. Clear analytical expressions for crack opening and stresses along the bimaterial interface as well as for stress intensity factors are derived. Furthermore, a system of the transcendental equations for the determination of the contact zones lengths in the sense of Comninou is obtained. The numerical analysis performed for a certain piezoelectric bimaterial showed an essential influence of the changing poling direction on the stress intensity factors and the contact zones lengths as well as on the crack opening displacement and the variations in the stress at the bimaterial interface.

分析了任意极化方向的两种不同压电材料界面裂纹在远载荷作用下的平面应变问题。假定裂纹尖端处无机械摩擦和电渗透接触区,而裂纹的其余部分被认为是电渗透的。本文采用了一种解析方法,将场变量用截面解析向量函数表示出来,并给出了Dirichlet-Riemann联合边值问题的表达式和解。导出了裂纹开度、双材料界面应力及应力强度因子的解析表达式。进一步,得到了在Comninou意义下确定接触区长度的超越方程组。对某压电双材料的数值分析表明,极化方向的变化对材料的应力强度因子、接触区长度、裂纹张开位移和双材料界面应力变化有重要影响。
{"title":"The influence of poling direction on an interface crack with contact zones in a piezoelectric bimaterial","authors":"Volodymyr Govorukha,&nbsp;Tetyana Kagadiy,&nbsp;Marc Kamlah","doi":"10.1007/s00419-025-02983-5","DOIUrl":"10.1007/s00419-025-02983-5","url":null,"abstract":"<div><p>The plane strain problem for an interface crack between two different piezoelectric materials with an arbitrary polarization direction under remote mechanical loading is analyzed. Mechanically frictionless and electrically permeable contact zones are assumed at the crack tips, and the remaining part of the crack is considered as electrically permeable. An analytical method is used based upon the presentations of the field variables via sectionally analytic vector functions with the resulting formulation and the solution of the combined Dirichlet–Riemann boundary value problem. Clear analytical expressions for crack opening and stresses along the bimaterial interface as well as for stress intensity factors are derived. Furthermore, a system of the transcendental equations for the determination of the contact zones lengths in the sense of Comninou is obtained. The numerical analysis performed for a certain piezoelectric bimaterial showed an essential influence of the changing poling direction on the stress intensity factors and the contact zones lengths as well as on the crack opening displacement and the variations in the stress at the bimaterial interface.</p></div>","PeriodicalId":477,"journal":{"name":"Archive of Applied Mechanics","volume":"95 11","pages":""},"PeriodicalIF":2.5,"publicationDate":"2025-11-04","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145456344","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
引用次数: 0
Thickness prediction and uniformity optimization in AFP of composite conical shells via fiber accumulation modeling 基于纤维积累模型的复合材料锥形壳AFP厚度预测及均匀性优化
IF 2.5 3区 工程技术 Q2 MECHANICS Pub Date : 2025-11-03 DOI: 10.1007/s00419-025-02976-4
Zhiheng Wang, Haojie Xu, Jun Hu, Zhou Yu

This study tackles the problem of wall thickness nonuniformity in the automated fiber placement (AFP) of composite conical shells by introducing a layer-wise control strategy grounded in a fiber accumulation model. A mathematical model for the single-layer fiber placement path is constructed, and the equivalent area method is utilized to quantitatively describe the overlap ratio of prepreg tows, from which the axial ply thickness distribution function is derived. To achieve uniform thickness, a multi-layer placement strategy with decreasing placement heights is introduced, establishing a mapping relationship between layer termination positions and cumulative thickness. A parametric analysis reveals that achieving a target thickness of 4 mm necessitates 28 layers, with fiber placement lengths decreasing beyond the 15th layer. The optimized scheme results in a mean squared error (MSE) of 0.021 mm in the total thickness. Based on the Tsai-Wu failure criterion, the predicted ultimate internal pressure is 4.2 MPa, slightly lower than the pre-optimization value of 4.8 MPa. This indicates that the optimized conical shell experiences only a marginal decline in failure performance while achieving a 22.35% reduction in material usage and still satisfying the strength requirements. Finite element analysis using Abaqus shows that under an internal pressure load of 4 MPa, the maximum displacement magnitude is 4.1 × 10–2 mm. This method effectively alleviates fiber accumulation induced by curvature gradients and offers a theoretical foundation for the high-precision manufacturing of thin-walled composite structures in aerospace applications.

本文通过引入基于纤维积累模型的分层控制策略,解决了复合材料锥形壳体自动铺放纤维时壁厚不均匀的问题。建立了单层纤维铺放路径的数学模型,利用等效面积法定量描述了预浸料束的重叠率,并由此导出了轴向铺层厚度分布函数。为了实现厚度的均匀化,引入了一种逐渐降低放置高度的多层放置策略,建立了层终止位置与累积厚度之间的映射关系。参数分析表明,达到4毫米的目标厚度需要28层,纤维放置长度减少到第15层以上。优化方案得到的总厚度均方误差(MSE)为0.021 mm。基于Tsai-Wu破坏准则,预测极限内压为4.2 MPa,略低于预优化值4.8 MPa。这表明,优化后的锥形壳破坏性能仅略有下降,但材料用量减少了22.35%,仍然满足强度要求。利用Abaqus进行有限元分析表明,在4 MPa的内压荷载作用下,最大位移量为4.1 × 10 - 2mm。该方法有效地缓解了曲率梯度引起的纤维堆积,为航空航天应用中薄壁复合材料结构的高精度制造提供了理论基础。
{"title":"Thickness prediction and uniformity optimization in AFP of composite conical shells via fiber accumulation modeling","authors":"Zhiheng Wang,&nbsp;Haojie Xu,&nbsp;Jun Hu,&nbsp;Zhou Yu","doi":"10.1007/s00419-025-02976-4","DOIUrl":"10.1007/s00419-025-02976-4","url":null,"abstract":"<div><p>This study tackles the problem of wall thickness nonuniformity in the automated fiber placement (AFP) of composite conical shells by introducing a layer-wise control strategy grounded in a fiber accumulation model. A mathematical model for the single-layer fiber placement path is constructed, and the equivalent area method is utilized to quantitatively describe the overlap ratio of prepreg tows, from which the axial ply thickness distribution function is derived. To achieve uniform thickness, a multi-layer placement strategy with decreasing placement heights is introduced, establishing a mapping relationship between layer termination positions and cumulative thickness. A parametric analysis reveals that achieving a target thickness of 4 mm necessitates 28 layers, with fiber placement lengths decreasing beyond the 15th layer. The optimized scheme results in a mean squared error (MSE) of 0.021 mm in the total thickness. Based on the Tsai-Wu failure criterion, the predicted ultimate internal pressure is 4.2 MPa, slightly lower than the pre-optimization value of 4.8 MPa. This indicates that the optimized conical shell experiences only a marginal decline in failure performance while achieving a 22.35% reduction in material usage and still satisfying the strength requirements. Finite element analysis using Abaqus shows that under an internal pressure load of 4 MPa, the maximum displacement magnitude is 4.1 × 10<sup>–2</sup> mm. This method effectively alleviates fiber accumulation induced by curvature gradients and offers a theoretical foundation for the high-precision manufacturing of thin-walled composite structures in aerospace applications.</p></div>","PeriodicalId":477,"journal":{"name":"Archive of Applied Mechanics","volume":"95 11","pages":""},"PeriodicalIF":2.5,"publicationDate":"2025-11-03","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145456467","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
引用次数: 0
Bifurcations and chaotic geometrically nonlinear vibrations of parametrically excited beams with breathing cracks 带呼吸裂纹的参数激励梁的分岔和混沌几何非线性振动
IF 2.5 3区 工程技术 Q2 MECHANICS Pub Date : 2025-11-03 DOI: 10.1007/s00419-025-02969-3
S. Malyshev, K. Avramov

The parametric vibrations of the flexible beams with account of the geometrically nonlinear deformations and breathing crack are studied. The crack breathing is described by using contact parameter and crack function. The nonlinear integro-differential equation, which describes geometrically nonlinear parametric vibrations of the beams with breathing cracks, is derived from the variational principal of Hu–Washizu. The nonlinear integro-differential equation is transformed into the system of piecewise—nonlinear ordinary differential equations by using the weighted residuals method. Arc length continuation technique and shooting technique are used jointly to study numerically nonlinear oscillations, their stability and bifurcations. The Neimark–Sacker and the saddle-nodes bifurcations are observed due to numerical analysis. Quasi-periodic and chaotic oscillations are studied numerically.

研究了考虑几何非线性变形和呼吸裂纹的柔性梁的参数振动问题。用接触参数和裂纹函数来描述裂纹呼吸。根据胡和柱变分原理,导出了描述含呼吸裂纹梁几何非线性参数振动的非线性积分微分方程。利用加权残差法将非线性积分微分方程转化为分段非线性常微分方程系统。采用弧长连续技术和射击技术,研究了数值非线性振荡及其稳定性和分岔问题。通过数值分析,观察到neimmark - sacker分岔和鞍节点分岔。对拟周期振荡和混沌振荡进行了数值研究。
{"title":"Bifurcations and chaotic geometrically nonlinear vibrations of parametrically excited beams with breathing cracks","authors":"S. Malyshev,&nbsp;K. Avramov","doi":"10.1007/s00419-025-02969-3","DOIUrl":"10.1007/s00419-025-02969-3","url":null,"abstract":"<div><p>The parametric vibrations of the flexible beams with account of the geometrically nonlinear deformations and breathing crack are studied. The crack breathing is described by using contact parameter and crack function. The nonlinear integro-differential equation, which describes geometrically nonlinear parametric vibrations of the beams with breathing cracks, is derived from the variational principal of Hu–Washizu. The nonlinear integro-differential equation is transformed into the system of piecewise—nonlinear ordinary differential equations by using the weighted residuals method. Arc length continuation technique and shooting technique are used jointly to study numerically nonlinear oscillations, their stability and bifurcations. The Neimark–Sacker and the saddle-nodes bifurcations are observed due to numerical analysis. Quasi-periodic and chaotic oscillations are studied numerically.</p></div>","PeriodicalId":477,"journal":{"name":"Archive of Applied Mechanics","volume":"95 11","pages":""},"PeriodicalIF":2.5,"publicationDate":"2025-11-03","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145456468","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
引用次数: 0
Dynamics and stability in the perturbed photogravitational CR3BP with disc structures and continued fractional effects 具有盘状结构和持续分数效应的摄动光引力CR3BP的动力学和稳定性
IF 2.5 3区 工程技术 Q2 MECHANICS Pub Date : 2025-11-01 DOI: 10.1007/s00419-025-02973-7
Lata Kumari Bairwa, Ashok Kumar Pal, Bhupendra Jangid, Sergey Ershkov, Elbaz I. Abouelmagd

This study examines the perturbed photogravitational Circular Restricted Three-Body Problem (CR3BP) under the combined effects of radiation pressure, albedo, continued fractional perturbations, and a gravitational disc. For the Sun–Earth system with the asteroid belt and the Sun–Jupiter system with the Kuiper belt, we analyze how the modified potential influences equilibrium points, Zero-Velocity Curves (ZVCs), and stability. Numerical results show that as the continued fractional parameter (epsilon ) increases, new collinear ((L_1', L_2')) and non-collinear ((L_4', L_5')) points emerge, while disc mass and other perturbations merely shift their locations. The Jacobi constant decreases with increasing (epsilon ), enlarging the allowed motion regions. Stability analysis indicates that (L_1, L_3, L_4', L_5') are unstable, whereas (L_4, L_5) remain stable within certain ranges of the critical mass ratio (mu _c). The new point (L_1') can be stable for specific (epsilon ) values in both systems. Overall, continued fractional perturbations, disc effects, and surface properties reshape equilibrium configurations and modify classical stability boundaries, enriching celestial dynamics and space mission design.

本文研究了在辐射压力、反照率、连续分数扰动和引力盘的综合作用下的摄动光引力圆形受限三体问题(CR3BP)。对于带小行星带的太阳-地球系统和带柯伊伯带的太阳-木星系统,我们分析了修正势对平衡点、零速度曲线(zvc)和稳定性的影响。数值结果表明,随着分数参数(epsilon )的增加,新的共线点((L_1', L_2'))和非共线点((L_4', L_5'))出现,而圆盘质量和其他扰动只是移动了它们的位置。雅可比常数随(epsilon )的增大而减小,使得允许的运动区域增大。稳定性分析表明(L_1, L_3, L_4', L_5')是不稳定的,而(L_4, L_5)在一定的临界质量比(mu _c)范围内保持稳定。对于两个系统中的特定(epsilon )值,新的点(L_1')可以是稳定的。总的来说,持续的分数扰动、圆盘效应和表面性质重塑了平衡构型,修改了经典的稳定性边界,丰富了天体动力学和空间任务设计。
{"title":"Dynamics and stability in the perturbed photogravitational CR3BP with disc structures and continued fractional effects","authors":"Lata Kumari Bairwa,&nbsp;Ashok Kumar Pal,&nbsp;Bhupendra Jangid,&nbsp;Sergey Ershkov,&nbsp;Elbaz I. Abouelmagd","doi":"10.1007/s00419-025-02973-7","DOIUrl":"10.1007/s00419-025-02973-7","url":null,"abstract":"<div><p>This study examines the perturbed photogravitational Circular Restricted Three-Body Problem (CR3BP) under the combined effects of radiation pressure, albedo, continued fractional perturbations, and a gravitational disc. For the Sun–Earth system with the asteroid belt and the Sun–Jupiter system with the Kuiper belt, we analyze how the modified potential influences equilibrium points, Zero-Velocity Curves (ZVCs), and stability. Numerical results show that as the continued fractional parameter <span>(epsilon )</span> increases, new collinear (<span>(L_1', L_2')</span>) and non-collinear (<span>(L_4', L_5')</span>) points emerge, while disc mass and other perturbations merely shift their locations. The Jacobi constant decreases with increasing <span>(epsilon )</span>, enlarging the allowed motion regions. Stability analysis indicates that <span>(L_1, L_3, L_4', L_5')</span> are unstable, whereas <span>(L_4, L_5)</span> remain stable within certain ranges of the critical mass ratio <span>(mu _c)</span>. The new point <span>(L_1')</span> can be stable for specific <span>(epsilon )</span> values in both systems. Overall, continued fractional perturbations, disc effects, and surface properties reshape equilibrium configurations and modify classical stability boundaries, enriching celestial dynamics and space mission design.</p></div>","PeriodicalId":477,"journal":{"name":"Archive of Applied Mechanics","volume":"95 11","pages":""},"PeriodicalIF":2.5,"publicationDate":"2025-11-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145405605","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
引用次数: 0
期刊
Archive of Applied Mechanics
全部 Acc. Chem. Res. ACS Applied Bio Materials ACS Appl. Electron. Mater. ACS Appl. Energy Mater. ACS Appl. Mater. Interfaces ACS Appl. Nano Mater. ACS Appl. Polym. Mater. ACS BIOMATER-SCI ENG ACS Catal. ACS Cent. Sci. ACS Chem. Biol. ACS Chemical Health & Safety ACS Chem. Neurosci. ACS Comb. Sci. ACS Earth Space Chem. ACS Energy Lett. ACS Infect. Dis. ACS Macro Lett. ACS Mater. Lett. ACS Med. Chem. Lett. ACS Nano ACS Omega ACS Photonics ACS Sens. ACS Sustainable Chem. Eng. ACS Synth. Biol. Anal. Chem. BIOCHEMISTRY-US Bioconjugate Chem. BIOMACROMOLECULES Chem. Res. Toxicol. Chem. Rev. Chem. Mater. CRYST GROWTH DES ENERG FUEL Environ. Sci. Technol. Environ. Sci. Technol. Lett. Eur. J. Inorg. Chem. IND ENG CHEM RES Inorg. Chem. J. Agric. Food. Chem. J. Chem. Eng. Data J. Chem. Educ. J. Chem. Inf. Model. J. Chem. Theory Comput. J. Med. Chem. J. Nat. Prod. J PROTEOME RES J. Am. Chem. Soc. LANGMUIR MACROMOLECULES Mol. Pharmaceutics Nano Lett. Org. Lett. ORG PROCESS RES DEV ORGANOMETALLICS J. Org. Chem. J. Phys. Chem. J. Phys. Chem. A J. Phys. Chem. B J. Phys. Chem. C J. Phys. Chem. Lett. Analyst Anal. Methods Biomater. Sci. Catal. Sci. Technol. Chem. Commun. Chem. Soc. Rev. CHEM EDUC RES PRACT CRYSTENGCOMM Dalton Trans. Energy Environ. Sci. ENVIRON SCI-NANO ENVIRON SCI-PROC IMP ENVIRON SCI-WAT RES Faraday Discuss. Food Funct. Green Chem. Inorg. Chem. Front. Integr. Biol. J. Anal. At. Spectrom. J. Mater. Chem. A J. Mater. Chem. B J. Mater. Chem. C Lab Chip Mater. Chem. Front. Mater. Horiz. MEDCHEMCOMM Metallomics Mol. Biosyst. Mol. Syst. Des. Eng. Nanoscale Nanoscale Horiz. Nat. Prod. Rep. New J. Chem. Org. Biomol. Chem. Org. Chem. Front. PHOTOCH PHOTOBIO SCI PCCP Polym. Chem.
×
引用
GB/T 7714-2015
复制
MLA
复制
APA
复制
导出至
BibTeX EndNote RefMan NoteFirst NoteExpress
×
0
微信
客服QQ
Book学术公众号 扫码关注我们
反馈
×
意见反馈
请填写您的意见或建议
请填写您的手机或邮箱
×
提示
您的信息不完整,为了账户安全,请先补充。
现在去补充
×
提示
您因"违规操作"
具体请查看互助需知
我知道了
×
提示
现在去查看 取消
×
提示
确定
Book学术官方微信
Book学术文献互助
Book学术文献互助群
群 号:604180095
Book学术
文献互助 智能选刊 最新文献 互助须知 联系我们:info@booksci.cn
Book学术提供免费学术资源搜索服务,方便国内外学者检索中英文文献。致力于提供最便捷和优质的服务体验。
Copyright © 2023 Book学术 All rights reserved.
ghs 京公网安备 11010802042870号 京ICP备2023020795号-1