一种半解析Legendre–Ritz方法用于动态分析阶梯功能梯度圆柱壳

《Vibration》：A Semi-Analytical Legendre–Ritz Method to Dynamically Analyze a Stepped Functionally Graded Cylindrical Shell

　　

本研究采用Legendre–Ritz方法引入了阶梯功能梯度(FG)圆柱壳在一般边界条件下的动态特性。计算模型基于一阶剪切变形理论(FSDT)和域分解方法(DDM)建立，并引入人工弹簧来模拟边界条件并确保段连续性。采用Legendre多项式和傅里叶级数构成容许位移函数。利用Rayleigh–Ritz方法确定了阶梯FG圆柱壳的自由振动和强迫振动特性。针对各种边界条件、材料参数和几何尺寸给出了结果，并与已发表研究进行了比较。该方法表现出良好的准确性，为分析阶梯FG圆柱壳的振动行为提供了基础。

复合材料圆柱壳作为基础承载构件，在航空航天、海洋工程及水下航行器、航天器、导弹系统等先进技术领域广泛应用。由于阶梯功能梯度(FG)圆柱壳兼具厚度不连续和材料梯度特性，其振动行为较均匀壳更为复杂，但现有研究多集中于各向同性或均匀FG圆柱壳，对阶梯FG圆柱壳在一般边界条件下的动态特性探索不足，尤其缺乏统一高效的半解析方法。针对这一空白，研究人员在《Vibration》上发表论文，提出一种基于Legendre–Ritz方法的半解析框架，系统研究阶梯FG圆柱壳的自由与强迫振动，为工程设计与优化提供理论指导。

研究人员采用一阶剪切变形理论(FSDT)建立位移场，基于域分解方法(DDM)将阶梯壳划分为多个薄段。在段间和边界处引入人工弹簧（刚度系数k_u, k_v, k_w, k_x, k_θ），统一模拟任意边界条件和几何连续性。位移函数在轴向采用Legendre正交多项式展开，周向采用傅里叶级数展开，以保证快速收敛和数值稳定性。利用Rayleigh–Ritz方法构建Lagrange能量函数并变分得到刚度矩阵K和质量矩阵M，求解特征值问题获得固有频率与振型。通过MATLAB自编程实现，并与文献对比验证精度（误差极小）。该方法无需生物或化学样本，主要基于数值仿真。

研究结果分为以下几个方面。

**模型验证与收敛性分析**：通过将计算结果与已发表文献中的频率参数对比（表1、表2），发现两者高度一致，验证了方法的准确性。截断系数研究表明，当轴向截断数M=8、周向最高波数N=8时，解达到稳定收敛，体现了Legendre多项式良好的正交性与数值稳定性。

**边界条件对动态特性的影响**：通过改变边界弹簧刚度从0到10¹⁵（对应自由到固支），发现增大边界约束可显著提高结构刚度，使固有频率增大且共振峰右移；约束减弱时共振峰数量增多，振动行为更复杂。

**长度-半径比对动态特性的影响**：保持半径不变，增大长度-半径比L/R导致整体刚度下降，固有频率降低，共振峰左移。同时给出了前五阶振型图，直观展示了长径比变化对变形形态的影响。

**厚度分布对动态特性的影响**：改变阶梯厚度比h₂/h₁，发现厚度增大时固有频率增大，共振峰右移，表明厚度分布可精确调节动态响应，是几何优化的直接依据。

**幂律指数对动态特性的影响**：随着幂律指数p增大（陶瓷成分减少、金属成分增加），材料刚度下降，固有频率逐渐降低，共振峰左移。材料梯度参数可作为核心设计变量用于振动特性优化。

讨论部分：该研究证实了边界约束、几何参数和材料梯度对阶梯FG圆柱壳动态特性的调控作用，但未涉及非线性、振动-声耦合等复杂工况。结论部分翻译如下：

本研究采用Legendre–Ritz方法研究了阶梯功能梯度(FG)圆柱壳的动态特性。结合一阶剪切变形理论、域分解方法和人工弹簧技术增强了当前方法的收敛性和有效性。利用Ritz方法可以评估自由振动和强迫振动行为。通过与已发表文献的比较验证了数值模型的准确性和可靠性。主要结论如下：

（1）边界约束强度和人工弹簧刚度是控制阶梯功能梯度圆柱壳振动特性的关键因素。加强边界约束可显著提高结构刚度和固有频率。

（2）长径比和阶梯厚度分布等几何参数对结构的动态特性具有显著的调控作用。长径比的增加会降低结构整体刚度，导致固有频率降低和共振峰左移。阶梯厚度比可以精确调节固有频率和共振响应，为几何优化提供了直接依据。

（3）随着幂律指数p值的增加，结构的固有频率逐渐降低。材料梯度参数可作为阶梯FG圆柱壳振动特性优化的核心设计变量。

本文提出的Legendre–Ritz半解析框架结合了域分解和非均匀结构容差建模技术，对几何不连续和材料非均匀的阶梯FG壳具有鲁棒性，可进一步扩展到涉及振动-声耦合的几何/材料非线性和复杂阶梯FG壳的动态分析。