热环境下石墨烯纳米片增强多孔板振动分析的半解析公式

《Results in Engineering》：A semi-analytical formulation for vibration analysis of graphene platelet reinforced porous plates under thermal environments

　　

功能梯度石墨烯纳米片增强金属泡沫板 日益被考虑用于在热和弹性地基环境中工作的轻质结构部件。然而，在一般边界约束下准确预测其动态特性对工程设计而言仍具挑战性。本研究开发了一个统一的分析框架，用于评估在热环境和任意边界条件下放置于双参数弹性地基上的GPLRMF矩形板的自由振动行为。基于Kirchhoff板理论和哈密顿原理，推导出控制方程并采用改进傅里叶级数法进行求解，其中引入人工弹簧刚度以统一模拟一般边界约束。所提公式能够高效、准确地预测多种边界条件下的固有频率和振型，而无需重新构建求解流程。研究人员进行了全面的参数化研究，以量化温升、地基刚度、孔隙率系数、石墨烯重量分数、长宽比和剪切层刚度对结构动态性能的影响。结果表明，热载荷与地基参数之间存在强烈的耦合效应，这显著影响了刚度退化和频率演变。通过将改进傅里叶级数法与人工弹簧边界技术相结合，开发了一个统一的分析框架，能够对具有热环境和弹性地基效应的任意弹性边界条件下的GPLRMF矩形板进行精确的振动分析。与传统公式不同，当考虑不同的边界约束时，该框架无需重新推导。该模型与现有结果进行了对比验证，证实了其鲁棒性和计算效率。本框架为复杂服役环境下先进轻质复合材料结构部件的振动分析和设计优化提供了实用且可扩展的工具。

随着材料科学与制造技术的飞速发展，工程材料面临着日益严苛的性能要求。传统金属材料难以同时满足这些需求，从而推动了功能梯度材料 的发展。功能梯度材料具有空间变化的材料属性，相比单一金属或陶瓷具有更高的比强度和比刚度，被广泛应用于航空航天、声学、土木和海洋结构中。在实际工程应用中，结构部件常处于弹性地基和高温等复杂服役环境中，这易引发疲劳失效和渐进性损伤，从而降低结构可靠性和使用寿命。石墨烯纳米片 因其极高的弹性模量和拉伸强度，被视为有效的增强剂，将其引入基体能显著提升板结构的力学性能。因此，研究石墨烯纳米片增强金属泡沫板的自由振动行为具有重要意义。尽管已有大量关于传统FGM矩形板在温度和经典边界条件下的振动研究，但针对在热环境和任意边界条件下置于弹性地基上的石墨烯纳米片增强金属泡沫复合矩形板的自由振动研究仍然有限。现有研究多集中于特定结构配置、简化边界条件或有限的地基模型。为解决这些局限性，研究人员开展了相关研究，开发了一种统一的分析框架，能够处理任意弹性边界条件，无需针对不同支撑配置重新推导。该研究将石墨烯增强、孔隙率分布、热效应和弹性地基纳入单一力学模型，以更真实地表示耦合服役环境下的实际工程结构。研究得出的结论表明，热载荷与地基参数之间存在强烈的耦合效应，显著影响刚度退化和频率演变。该统一框架为复杂服役环境下先进轻质复合材料结构部件的振动分析和设计优化提供了实用且可扩展的工具，具有重要的工程意义，该论文发表在《Results in Engineering》上。

研究人员基于Kirchhoff薄板理论（适用于薄板小挠度弯曲的经典理论）和哈密顿原理（推导系统动力学方程的变分原理），推导了系统的控制方程。在求解方法上，采用改进傅里叶级数法 表示位移场，通过引入辅助三角项确保位移及其导数在边缘的连续性。同时，引入人工弹簧边界技术，在板的四条边缘布置平动和转动弹簧，通过赋予不同的弹簧刚度系数来统一模拟任意弹性边界条件。结合改进的Halpin-Tsai微观力学模型（用于预测复合材料有效弹性模量的方法）和混合律扩展法则，计算了GPLRMF板的有效材料属性。通过求解广义特征值问题获得结构的固有频率和振型。

4.1. 孔隙率系数变化下的频率响应
研究人员通过参数化分析得出，随着孔隙率系数 e₁ 从0增加到0.5，由于孔隙生长引起的有效刚度逐渐退化，大多数分布组合的固有频率呈下降趋势。但在GPL-C分布与Porosity-I（孔隙分布模式一）组合中，固有频率出现非单调变化，这表明石墨烯纳米片的刚度增强与孔隙演化导致的刚度损失之间存在竞争效应。在三种石墨烯纳米片分布模式中，GPL-A分布（纳米片富集于板上下表面）对应的固有频率最高，因为表面聚集的石墨烯有效增强了抗弯刚度，缓解了孔隙的不利影响。而GPL-B分布（纳米片集中于中性面）的固有频率最低。

4.2. 石墨烯重量分数的影响
随着石墨烯重量分数 W_gpl 从0.002增加到0.01，所有分布组合的固有频率均呈现明显的上升趋势。这表明增加石墨烯纳米片的含量能够有效增强金属泡沫材料的刚度。GPL-A分布始终产生最高的固有频率，且与Porosity-I结合时频率变化最陡峭。GPL-B分布的频率最低，且当与Porosity-II（孔隙分布模式二）结合时，频率随 W_gpl 的变化相对平缓，说明集中于中性面的石墨烯增强效率有限。

4.3. 热机械耦合效应
随着温度 T₀ 从300 K升至800 K，所有GPL分布模式的固有频率单调增加。在 T₀ = 600 K 时，GPL-A在Porosity-I和Porosity-III（孔隙分布模式三）下表现出比GPL-C更高的一阶固有频率。这种频率随温度升高的现象表明，由石墨烯纳米片分布引起的温度相关刚化效应在所研究的温度范围内超过了热软化效应。此外，当温度增量 ΔT 从0增加到100 K时，固有频率同样呈增加趋势，这主要归因于GPL增强金属泡沫基体有效刚度的温度依赖性增强。

4.4. 几何参数和边界条件的影响
随着长宽比 λ 从0.2增加到1.2，所有分布模式的固有频率单调增加。在固定长宽比下，GPL-A与Porosity-I组合表现出最强的协同增强效应。此外，在不同的边界条件下（如 E₁E₁E₁E₁、SSSF、SSSS、SCSC和CCCC等），除SSSF（三边简支一边自由）边界条件外，Porosity-I与GPL-A组合的固有频率始终最高，而Porosity-II与GPL-B组合的固有频率最低。随着边界约束的增强，板的总体结构刚度提升，固有频率随之增加。

4.5. 结构刚度参数的敏感性分析
随着无量纲弹性刚度系数 K_w 从0增加到500，所有分布模式的固有频率单调增加。对于固定 K_w 值，GPL-A与Porosity-I组合表现出最强的协同增强效应。此外，随着无量纲剪切刚度系数 Q_w 的增加，固有频率同样单调增加。地基提供的横向约束抑制了板的挠度和剪切变形，增强了全局结构刚度。在所有考虑的工况中，Porosity-II与GPL-A组合在整个 Q_w 范围内始终产生最高的固有频率。

讨论部分系统总结分析了各种参数对GPLRMF矩形板自由振动特性的影响。研究结论可以概括为以下几点：第一，结合改进傅里叶级数和人工弹簧边界技术的统一分析框架被成功开发，消除了改变边界约束时重新推导的需要，能够高效准确分析任意弹性边界条件下GPLRMF矩形板的自由振动。第二，随着 T₀ 和 ΔT 的增加，由于复合板有效刚度的温度诱导改变，所有分布模式的固有频率均增加。第三，随着孔隙率系数 e₁ 的增加，在GPL-C和Porosity-I配置下固有频率几乎保持不变，而其他分布模式呈单调下降趋势；当孔隙率系数从0增加到0.5时，一阶固有频率的最大下降幅度约为8.6%。第四，本分析框架为预测在耦合热环境和弹性地基效应下具有任意弹性边界条件的GPLRMF矩形板的自由振动特性提供了高效准确的工具，为轻质复合材料结构的设计和优化提供了有益参考。