黏弹性 Pasternak 地基上修正 Timoshenko 梁 横向自振特性分析

柳 伟，汪过兵，赵志鹏，赵晓军; LIU Wei; WANG Guo-bin; ZHAO Zhi-peng; ZHAO Xiao-jun

使用Chrome浏览器效果最佳，继续浏览，你可能不会看到最佳的展示效果，

确定继续浏览么?

复制成功，请在其他浏览器进行阅读

黏弹性Pasternak地基上修正Timoshenko梁横向自振特性分析 PDF

- ORCID：
柳伟 ¹
✉
- ORCID：
汪过兵 ²
- ORCID：
赵志鹏 ¹
- ORCID：
赵晓军 ¹

1. 兰州信息科技学院土木工程学院，甘肃兰州 730300； 2. 西安理工大学岩土工程研究所，陕西西安 710048

中图分类号： TU471⁺.2； TU348

最近更新：2024-09-04

DOI：10.16385/j.cnki.issn.1004-4523.2024.08.007

摘要

基于修正Timoshenko梁理论，建立黏弹性Pasternak地基上修正Timoshenko梁的横向振动控制方程，运用回传射线矩阵法推导出黏弹性Pasternak地基上两端简支修正Timoshenko梁自振频率和衰减系数的解析解，结合二分法和黄金分割法计算了经典边界条件下黏弹性Pasternak地基上修正Timoshenko梁的自振特性，对比分析了考虑剪切变形引起的转动惯量、梁长及不同的边界条件对结构自振特性的影响。研究表明：黏弹性Pasternak地基上修正Timoshenko梁的各阶自振频率和衰减系数小于经典Timoshenko梁的各阶自振频率和衰减系数；梁越短，剪切变形引起的转动惯量对结构自振频率和衰减系数的影响越大，且对高阶的影响明显大于低阶；边界约束条件越强，振动能量衰减越明显。

关键词

黏弹性Pasternak地基; 修正Timoshenko梁; 回传射线矩阵法; 解析解; 边界条件

引言

黏弹性地基梁在土木工程领域应用广泛^［

1‑2］，如机场、铁路、公路、高层建筑基础中均有涉及。而结构的自振特性是反映结构动力特性的重要物理量，通过对其进行准确高效的研究不仅对结构设计和施工计算提供理论基础，而且对避免发生共振、疲劳问题和对结构的减振有重要的帮助，因此，地基梁的动力学特性在工程领域及学术界备受关注^{［参考文献 3‑5}3‑5］。

目前各种地基上经典的Timoshenko梁研究较为丰富，Malekzadeh等^［

6］采用微积分法研究了弹性地基上Timoshenko梁的自振特性的数值解；彭丽等^{［参考文献 7‑9}7‑9］运用复模态方法研究了Pasternak地基上Timoshenko梁的自振特性及任一初始激励条件下外激励的响应；余云燕等^{［参考文献 10

百度学术}10］求解了黏弹性Pasternak地基上Timoshenko梁在不同的约束条件下单跨及两跨连续地基梁的自振频率、衰减系数和模态。这些理论虽然有较大的研究空间，但梁在高跨比比较大、局部悬空、局部高度承载情况下，由于未考虑剪切变形引起转动惯量的影响，导致经典Timoshenko梁存在挠度关于时间的四阶导数项、第二频谱、物理意义不明确等问题，使计算结果有较大偏差。因此，陈镕等^{［参考文献 11

百度学术}11］对传统的Timoshenko梁的运动方程进行了修正，发现修正Timoshenko梁对高频段有较大影响，指出考虑梁剪切变形引起的转动惯量后，时间的四阶导数项自然会消失；夏桂云^{［参考文献 12

百度学术}12］利用固有频率和临界频率的关系论证了Timoshenko梁产生第二频谱的原因，通过实例验证Timoshenko梁第二频谱的存在，因此，准确、合理地对其进行计算具有重要的意义。吴晓等^{［参考文献 13

百度学术}13］应用Timoshenko梁修正理论推导了泡沫铝合金梁的自振频率表达式，并求解了在简谐荷载作用下强迫振动的解析解；王家乐等^{［参考文献 14

百度学术}14］基于修正Timoshenko梁理论，采用复模态分析法推导多种边界条件下弹性地基梁振动超越方程及模态函数；Li等^{［参考文献 15

百度学术}15］推导出分数阶标准固体黏弹性地基上修正Timoshenko梁的运动控制方程，得到了自振频率的解析解；徐梅玲等^{［参考文献 16

百度学术}16］采用分离变量法，给出Euler梁模型相对于修正Timoshenko梁模型的误差计算公式。但上述研究未充分考虑地基土体颗粒之间相互剪切的连续性，也未分析剪切变形引起的转动惯量、梁长及边界条件对黏弹性地基上修正Timoshenko梁自振频率、衰减系数和模态的影响。

本文将修正Timoshenko梁理论与黏弹性Pasternak地基进行组合，建立新的黏弹性地基梁振动控制方程，运用回传射线矩阵法解耦，得到黏弹性Pasternak地基中两端简支修正Timoshenko梁自振频率和衰减系数的解析解，对比分析了考虑剪切变形引起的转动惯量、梁长和不同的边界条件对结构自振特性的影响。从而为黏弹性地基梁振动分析和计算提供理论基础。

1 振动控制方程及方程的解

黏弹性地基上修正Timoshenko梁的力学模型如图1所示。地基与修正Timoshenko梁的相互作用采用考虑土体连续性的巴氏（Pasternak）模型，即在Winkler地基模型的基础上通过一层只能产生横向剪切变形而不可压缩的剪切层来实现。建立整体坐标系 $(x, v)$ ，引入2个对偶局部坐标系 $x^{12}$ 和 $x^{21}$ 。

图1 黏弹性Pasternak地基上修正Timoshenko梁力学模型

Fig.1 Mechanical model of modified Timoshenko beam on viscoelastic Pasternak foundation

对图1取微段隔离体，其受力情况如图2所示。

图2 微段隔离体受力分析图

Fig.2 Force analysis diagram of micro-segment isolator

根据达朗贝尔原理，对微段隔离体建立竖向力及力矩平衡略去高阶项，得：

\{\begin{array}{l} \frac{\partial V (x, t)}{\partial x} - q_{v} (x, t) = ρ A \frac{\partial^{2} v (x, t)}{\partial t^{2}} \\ \frac{\partial M (x, t)}{\partial x} + V (x, t) = ρ I_{z} \frac{\partial^{3} v (x, t)}{\partial x \partial t^{2}} \end{array}

（1）

式中 $ρ$ ， $I_{z}$ 分别为密度、横截面惯性矩； $V (x, t)$ 和 $M (x, t)$ 分别为截面剪力和弯矩， $V (x, t) =$ $k' A G [\frac{\partial v (x, t)}{\partial x} - ϕ (x, t)]$ ， $M (x, t) = - E I_{z} \frac{\partial ϕ (x, t)}{\partial x}$ ，其中 $v$ 和 $ϕ$ 分别为梁的总挠度和截面转角， $E$ ， $G$ ， $A$ 和 $k'$ 分别为弹性模量、剪切模量、横截面面积和截面剪切系数； $q_{v} (x, t)$ 为黏弹性地基梁的地基竖向反力， $q_{v} (x, t) = k_{v} v (x, t) + β_{v} \frac{\partial v (x, t)}{\partial t} - T_{v} \frac{\partial^{2} v (x, t)}{\partial x^{2}}$ ，其中 $k_{v}$ ， $β_{v}$ 和 $T_{v}$ 分别为土体弹簧系数、土体阻尼系数和地基剪切系数。

整理式（1），得到黏弹性Pasternak地基上修正Timoshenko梁的振动控制方程为：

\{\begin{array}{l} \begin{array}{l} k' A G [\frac{\partial^{2} v (x, t)}{\partial x^{2}} - \frac{\partial ϕ (x, t)}{\partial x}] = ρ A \frac{\partial^{2} v (x, t)}{\partial t^{2}} + \\ k_{v} v (x, t) + β_{v} \frac{\partial v (x, t)}{\partial t} - T_{v} \frac{\partial^{2} v (x, t)}{\partial x^{2}} \end{array} \\ \begin{array}{l} k' A G [\frac{\partial v (x, t)}{\partial x} - ϕ (x, t)] + E I_{z} \frac{\partial^{2} ϕ (x, t)}{\partial x^{2}} = \\ ρ I_{z} \frac{\partial^{3} v (x, t)}{\partial x \partial t^{2}} \end{array} \end{array}

（2）

对式（2）进行Fourier变换有：

\{\begin{array}{l} Q_{1} \hat{v} - {\bar{T}}_{v} \frac{d^{2} \hat{v}}{d {\bar{x}}^{2}} + \frac{d \hat{ϕ}}{d \bar{x}} - \frac{d^{2} \hat{v}}{d {\bar{x}}^{2}} = 0 \\ Q_{3} \frac{d \hat{v}}{d \bar{x}} + \frac{d^{2} \hat{ϕ}}{d {\bar{x}}^{2}} - Q_{2} (\hat{ϕ} - \frac{d \hat{v}}{d \bar{x}}) = 0 \end{array}

（3）

式中顶标“ $\hat{}$ ”和“-”表示频域中的变量； $i = \sqrt[]{- 1}$ ； $Q_{1} = {\bar{k}}_{v} + i ω {\bar{β}}_{v} - \frac{ω^{2}}{c_{1}^{2}}$ ， $Q_{2} = \frac{1}{α R_{z}^{2}}$ ， $Q_{3} = \frac{ω^{2}}{c_{0}^{2}}$ ， ${\bar{k}}_{v} = \frac{k_{v}}{k' A G}$ ， $ω = {\bar{ω}}_{n} + i δ_{n}$ 为圆频率， ${\bar{ω}}_{n}$ 和 $δ_{n}$ 分别为自振频率和衰减系数； ${\bar{β}}_{v} = \frac{β_{v}}{k' A G}$ ， ${\bar{T}}_{v} = \frac{T_{v}}{k' A G}$ ， $α = \frac{E}{k' G}$ ， $R_{z} = \sqrt[]{I_{z} / A}$ 为截面的回转半径， $c_{0} = \sqrt[]{E / ρ}$ 为纵波波速， $c_{1} = \sqrt[]{k' G / ρ}$ 为横波波速。

求解式（3）得：

\hat{v} (x, ω) = a_{1} (ω) e^{i k_{1} x} + d_{1} (ω) e^{- i k_{1} x} + a_{2} (ω) e^{i k_{2} x} + d_{2} (ω) e^{- i k_{2} x}

\hat{ϕ} (x, ω) = g_{1} a_{1} (ω) e^{i k_{1} x} - g_{1} d_{1} (ω) e^{- i k_{1} x} + g_{2} a_{2} (ω) e^{i k_{2} x} - g_{2} d_{2} (ω) e^{- i k_{2} x}

（4）

式中 $a_{1} (ω)$ ， $a_{2} (ω)$ 为待定的入射波波幅； $d_{1} (ω)$ ， $d_{2} (ω)$ 为待定的出射波波幅； $k_{1}$ ， $k_{2}$ 为波数，满足：

k_{j} (ω) = \sqrt[]{\frac{- b \pm {(- 1)}^{j + 1} \sqrt[]{b^{2} - 4 a c}}{2 a}}

，

j = 1,2

（5）

式中 $a = {\bar{T}}_{v} + 1$ ； $b = Q_{1} + ({\bar{T}}_{v} + 1) Q_{2} - (Q_{2} + Q_{3})$ ； $c = Q_{1} Q_{2}$ 。

对应于波数 $k_{1}$ 和 $k_{2}$ ， $\hat{v}$ 与 $\hat{ϕ}$ 的比值为：

g_{j} = - \frac{Q_{1} + ({\bar{T}}_{v} + 1) k_{j}^{2}}{i k_{j}^{2}} = \frac{i k_{_{j}} (Q_{2} + Q_{3})}{k_{j}^{2} + Q_{2}}

，

j = 1,2

（6）

弯矩和剪力在频域中的表达式为：

\{\begin{array}{l} \hat{M} (x, ω) = β_{1} a_{1} (ω) e^{i k_{1} x} + β_{1} d_{1} (ω) e^{- i k_{1} x} + \\ β_{2} a_{2} (ω) e^{i k_{2} x} + β_{2} d_{2} (ω) e^{- i k_{2} x} \\ \hat{V} (x, ω) = γ_{1} a_{1} (ω) e^{i k_{1} x} - γ_{1} d_{1} (ω) e^{- i k_{1} x} + \\ γ_{2} a_{2} (ω) e^{i k_{2} x} - γ_{2} d_{2} (ω) e^{- i k_{2} x} \end{array}

（7）

式中 $β_{j} = - i E I_{z} k_{j} g_{j}$ ， $γ_{j} = k' A G (i k_{j} - g_{j})$ ， $j = 1, 2$ 。

2 经典边界条件下的自振频率及模态求解

2.1 两端简支条件下黏弹性地基梁的自振频率方程

以两端简支的黏弹性修正Timoshenko梁为例，运用回传射线矩阵法推导自振频率方程，其自由振动的边界条件为：

\{\begin{matrix} {\hat{v}}^{12} (0, ω) = 0, & {\hat{v}}^{21} (0, ω) = 0 \\ {\hat{M}}^{12} (0, ω) = 0, & {\hat{M}}^{21} (0, ω) = 0 \end{matrix}

（8）

将式（4），（7）代入式（8），并整理成矩阵形式有：

\{\begin{matrix} d_{1}^{12} \\ d_{2}^{12} \\ d_{1}^{21} \\ d_{2}^{21} \end{matrix}\} = {[\begin{matrix} 1 & 1 & 0 & 0 \\ β_{1} & β_{2} & 0 & 0 \\ 0 & 0 & 1 & 1 \\ 0 & 0 & β_{1} & β_{2} \end{matrix}]}^{- 1} [\begin{matrix} - 1 & - 1 & 0 & 0 \\ - β_{1} & - β_{2} & 0 & 0 \\ 0 & 0 & - 1 & - 1 \\ 0 & 0 & - β_{1} & - β_{2} \end{matrix}] \{\begin{matrix} a_{1}^{12} \\ a_{2}^{12} \\ a_{1}^{21} \\ a_{2}^{21} \end{matrix}\}

（9）

式（9）可写为：

d = S a

（10）

式中 $d$ ， $a$ 分别为出射波和入射波波幅向量； $S$ 为散射矩阵。三者表达式如下：

d = {[d_{1}^{12} d_{2}^{12} d_{1}^{21} d_{2}^{21}]}^{T}

，

a = {[a_{1}^{12} a_{2}^{12} a_{1}^{21} a_{2}^{21}]}^{T}

，

S = {[\begin{matrix} 1 & 1 & 0 & 0 \\ β_{1} & β_{2} & 0 & 0 \\ 0 & 0 & 1 & 1 \\ 0 & 0 & β_{1} & β_{2} \end{matrix}]}^{- 1} [\begin{matrix} - 1 & - 1 & 0 & 0 \\ - β_{1} & - β_{2} & 0 & 0 \\ 0 & 0 & - 1 & - 1 \\ 0 & 0 & - β_{1} & - β_{2} \end{matrix}]

。

入射波波幅向量和出射波波幅向量之间的相位关系为：

a = P \tilde{d}

（11）

式中 $\tilde{d} = {[d_{1}^{21} d_{2}^{21} d_{1}^{12} d_{2}^{12}]}^{T}$ ； $P$ 为传播矩阵， $P = d i a g \{- e^{- i k_{1} l} - e^{- i k_{2} l} - e^{- i k_{1} l} - e^{- i k_{2} l}\}$ 。

$\tilde{d}$ 与 $d$ 中的元素完全相同，只是各元素的排列位置有所调整，引入置换矩阵 $U$ ：

\tilde{d} = U d

（12）

U = [\begin{matrix} 0 & 0 & 1 & 0 \\ 0 & 0 & 0 & 1 \\ 1 & 0 & 0 & 0 \\ 0 & 1 & 0 & 0 \end{matrix}]

（13）

将式（12）代入式（11）得：

a = P U d

（14）

将式（14）代入式（10）得：

d [I - R] = 0

（15）

式中 $R = S P U$ 为回传射线矩阵； $I$ 为单位矩阵。

要使式（15）中的 $d$ 有非零解，则其系数行列式为零，则：

|I - R| = (e^{2 i k_{1} l} - 1) (e^{2 i k_{2} l} - 1) = 0

（16）

利用指数函数与三角函数的关系，有：

\begin{array}{l} |I - R| = [c o s (2 k_{1} l) + i s i n (2 k_{1} l) - 1] [c o s (2 k_{2} l) + \\ i s i n (2 k_{2} l) - 1] = 0 \end{array}

（17）

忽略其奇异解，式（17）可解：

k = \frac{(n - 1) π}{l}

, n = 1,2, 3, \dots

（18）

式中 l为梁长。

将式（18）代入式（5）求解复系数一元二次方程，解得：

ω = \frac{{\bar{β}}_{v} (k^{2} + Q_{2})}{2 (\frac{k^{2} + Q_{2}}{c_{1}^{2}} + \frac{k^{2}}{c_{0}^{2}})} i + \frac{\sqrt[]{- {\bar{β}}_{v}^{2} (k^{2} + Q_{2})^{2} + 4 (\frac{k^{2} + Q_{2}}{c_{1}^{2}} + \frac{k^{2}}{c_{0}^{2}}) [({\bar{T}}_{v} + 1) (k^{4} + Q_{2} k^{2}) + {\bar{k}}_{v} (k^{2} + Q_{2}) - Q_{2} k^{2}]}}{2 (\frac{k^{2} + Q_{2}}{c_{1}^{2}} + \frac{k^{2}}{c_{0}^{2}})}

（19）

进一步得到两端简支边界条件下黏弹性Pasternak地基上修正Timoshenko梁的自振频率为：

{\bar{ω}}_{n} = \frac{\sqrt[]{- {\bar{β}}_{v}^{2} (k^{2} + Q_{2})^{2} + 4 (\frac{k^{2} + Q_{2}}{c_{1}^{2}} + \frac{k^{2}}{c_{0}^{2}}) [({\bar{T}}_{v} + 1) (k^{4} + Q_{2} k^{2}) + {\bar{k}}_{v} (k^{2} + Q_{2}) - Q_{2} k^{2}]}}{2 (\frac{k^{2} + Q_{2}}{c_{1}^{2}} + \frac{k^{2}}{c_{0}^{2}})}

（20）

衰减系数为：

δ_{n} = \frac{{\bar{β}}_{v} (k^{2} + Q_{2})}{2 (\frac{k^{2} + Q_{2}}{c_{1}^{2}} + \frac{k^{2}}{c_{0}^{2}})}

（21）

由线性代数，有：

[I - R] * a d j [I - R] = d e t [I - R] * I

（22）

式中 $a d j [I - R]$ 为矩阵 $[I - R]$ 的伴随矩阵； $d e t [I - R]$ 为矩阵 $[I - R]$ 的模。

当 $d e t [I - R]$ 中 $ω$ 的实部取自振频率 ${\bar{ω}}_{n}$ ，虚部取衰减系数 $δ_{n}$ 时，式（22）可表示为：

[I - R (ω_{k})] * a d j [I - R (ω_{k})] = 0

（23）

式中 $ω_{k}$ 为第k阶频率。

设 $N$ 是出射波波幅向量 $d_{k}$ 的维数， $T$ 为 $1 \leq T \leq N$ 的任意正整数，矩阵 $a d j [I - R]$ 的第 $T$ 列记为 $d_{k}$ ，则：

[I - R (ω_{k})] * d_{k} = 0

（24）

由式（24）可知， $a d j [I - R (ω_{k})]$ 的每一个非零列均可作为黏弹性地基梁在自由振动时非零出射波的波幅向量，求得 $d_{k}$ 以后，代入式（14）可求得 $a_{k}$ ，将 $a_{k}$ 和 $d_{k}$ 代入式（4）中，将各节点处的位移归一化处理后，得到黏弹性地基梁的模态曲线。

2.2 其他边界条件下黏弹性地基梁的自振频率方程

固定‑简支边界条件下黏弹性Pasternak地基上修正Timoshenko梁，其自由振动的边界条件为：

\{\begin{matrix} {\hat{v}}^{12} (0, ω) = 0, & {\hat{v}}^{21} (0, ω) = 0 \\ {\hat{ϕ}}^{12} (0, ω) = 0, & {\hat{M}}^{21} (0, ω) = 0 \end{matrix}

（25）

散射矩阵 $S$ 为：

S = {[\begin{matrix} 1 & 1 & 0 & 0 \\ g_{1} & g_{2} & 0 & 0 \\ 0 & 0 & 1 & 1 \\ 0 & 0 & β_{1} & β_{2} \end{matrix}]}^{- 1} [\begin{matrix} - 1 & - 1 & 0 & 0 \\ g_{1} & g_{2} & 0 & 0 \\ 0 & 0 & - 1 & - 1 \\ 0 & 0 & - β_{1} & - β_{2} \end{matrix}]

（26）

$[I - R]$ 为：

[I - R] = (g_{1} - g_{2} + g_{1} e^{2 i k_{1} l} - g_{1} e^{2 i k_{2} l} + g_{2} e^{2 i k_{1} l} - g_{2} e^{2 i k_{2} l} - g_{1} e^{2 i k_{1} l} e^{2 i k_{2} l} + g_{2} e^{2 i k_{1} l} e^{2 i k_{2} l}) = 0

（27）

式（27）是关于自振频率 ${\bar{ω}}_{n}$ 及衰减系数 $δ_{n}$ 的二维复数超越方程。经计算，其他边界条件下的自振频率方程也是隐式超越方程，在数学上此类问题只有数值解而没有解析解。考虑到以上因素，根据回传射线矩阵法的列式特点提出将二分法和黄金分割法（简称求根法）结合起来进行迭代求解复杂超越方程，即分别对 $ω$ 的实部 ${\bar{ω}}_{n}$ 和虚部 $δ_{n}$ 进行循环，当 $[I - R]$ 的模小于预先给定的误差时，取出 $ω$ 对应的实部 ${\bar{ω}}_{n}$ 和虚部 $δ_{n}$ ，得到其他边界条件下黏弹性地基上修正Timoshenko梁的自振频率和衰减系数。

3 算例分析

黏弹性Pasternak地基上修正Timoshenko梁计算模型如图1所示，修正Timoshenko梁的计算参数如表1所示，土体的计算参数参考文献［

17］中的数值，如表2所示。

表1 修正Timoshenko梁的计算参数

Tab.1 Calculation parameters of modified Timoshenko beam

弹性模量 $E / P a$	剪切模量 $G / P a$	截面尺寸 $A / m$ ²	梁长 $l / m$	惯性矩 $I_{z} / m^{4}$	密度 $ρ / (k g \cdot {m^{-}}^{3})$	截面剪切系数 $k'$
4.322×10¹⁰	1.751×10¹⁰	$1.5 \times 1.5$	6	0.422	2700	$π^{2} / 12$

表2 黏弹性地基的各项物理计算参数

Tab.2 Physical calculation parameters of viscoelastic foundation

弹簧系数

$k_{v} / (N \cdot {m^{-}}^{2})$

阻尼系数

$β_{v} / (N \cdot s \cdot {m^{-}}^{2})$

剪切系数

$T_{v} / N$

1×10⁶

1×10⁴

1×10⁷

3.1 不同梁理论对黏弹性Pasternak地基梁自振特性的影响

以两端简支的黏弹性Pasternak地基上修正Timoshenko梁为例，利用表1和2中的计算参数，改变梁长，得出不同梁理论下黏弹性Pasternak地基梁前8阶自振频率和衰减系数的精确解如表3和4所示。其中黏弹性Pasternak地基上经典Timoshenko梁（P‑T）和修正Timoshenko梁（P‑MT）理论的区别在于振动控制方程中是否考虑了剪切变形引起的转动惯量影响。表中误差为相对误差，其自振频率相对误差计算公式为 $ζ = ({\bar{ω}}_{P ‑ T} - {\bar{ω}}_{P ‑ M T}) / {\bar{ω}}_{P ‑ T}$ ，其衰减系数相对误差计算公式为 $ψ = (δ_{P ‑ T} - δ_{P ‑ M T}) / δ_{P ‑ T}$ 。

表3 两端简支黏弹性地基梁前8阶自振频率

Tab.3 The first eight order natural frequencies of simply supported viscoelastic foundation beam at both ends

阶数	l=30 m			l=6 m
阶数	${\bar{ω}}_{P ‑ T}$	${\bar{ω}}_{P ‑ M T}$	$ζ$ /%	${\bar{ω}}_{P ‑ T}$	${\bar{ω}}_{P ‑ M T}$	$ζ$ /%
1	12.804	12.804	0.00	12.804	12.804	0.00
2	23.228	23.187	0.18	434.390	432.699	0.39
3	76.323	76.145	0.23	1436.618	1407.195	2.05
4	166.013	165.557	0.27	2649.427	2531.576	4.45
5	286.854	286.000	0.30	3925.501	3668.667	6.54
6	434.390	432.699	0.39	5215.121	4791.268	8.13
7	604.326	600.823	0.58	6501.793	5898.138	9.28
8	792.713	785.937	0.85	7780.585	6992.438	10.13

由表3可知，黏弹性Pasternak地基上经典Timoshenko梁和修正Timoshenko梁理论第1阶自振频率相同，第2阶之后，随着阶数的增大，两者频率值的差距逐渐增大，并且修正Timoshenko梁的各阶自振频率小于经典Timoshenko梁的各阶自振频率。当梁长l=30 m时，高跨比为0.05，其7，8阶自振频率相对误差为0.58%和0.85%，当梁长l=6 m时，高跨比为0.25，其7，8阶自振频率相对误差分别为9.28%和10.13%。从相对误差的角度来看，在黏弹性Pasternak地基上考虑剪切变形引起转动惯量的修正Timoshenko梁与不考虑其影响的经典Timoshenko梁自振频率计算结果在高跨比比较大的高频段振动有较大偏差。所以在涉及黏弹性地基梁具体工程计算中，在梁对高跨比比较大的高频段计算结果有重大影响的分析中（诸如冲击等问题），应采用修正Timoshenko梁理论。

由表4可知，黏弹性Pasternak地基上经典Timoshenko梁和修正Timoshenko梁理论第1阶衰减系数相同，第2阶之后，随着阶数的增大，两者衰减系数值的差距逐渐增大，并且修正Timoshenko梁的各阶衰减系数小于经典Timoshenko梁的各阶衰减系数。梁长为6 m时，第2阶衰减系数相对误差值与梁长为30 m时第6阶衰减系数相对误差值相等，且梁长为6 m时，其7，8阶衰减系数相对误差分别高达41.97%和43.66%。所以在工程实践中涉及振动能量损耗控制时，应考虑剪切变形引起转动惯量的影响，否则会带来较大偏差。

表4 两端简支黏弹性地基梁前8阶衰减系数

Tab.4 The first eight order attenuation coefficients of simply supported viscoelastic foundation beam at both ends

阶数	l=30 m			l=6 m
阶数	$δ_{P ‑ T}$	$δ_{P ‑ M T}$	$ψ$ /%	$δ_{P ‑ T}$	$δ_{P ‑ M T}$	$ψ$ /%
1	0.8230	0.8230	0.00	0.8230	0.8230	0.00
2	0.8214	0.8180	0.41	0.7921	0.7260	8.34
3	0.8167	0.8037	1.59	0.7577	0.5958	21.38
4	0.8096	0.7819	3.42	0.7476	0.5205	30.38
5	0.8011	0.7552	5.73	0.7509	0.4809	35.96
6	0.7921	0.7260	8.34	0.7589	0.4588	39.55
7	0.7834	0.6965	11.09	0.7677	0.4455	41.97
8	0.7754	0.6680	13.85	0.7758	0.4371	43.66

图3和4为不同梁长情况下黏弹性Pasternak地基上经典Timoshenko梁与修正Timoshenko梁的前6阶模态。从图3和4可知，l=30 m时，黏弹性Pasternak地基上经典Timoshenko梁与修正Timoshenko梁的第1阶振型曲线之间的差距最大，振型峰值相同，从第2阶开始，随着阶数的增大，其振型曲线之间的差距逐渐减小，但振型峰值逐渐增大。l=6 m相对于l=30 m的短梁，黏弹性Pasternak地基上经典Timoshenko梁与修正Timoshenko梁振型峰值增大得越明显，表明梁越短，考虑剪切变形引起转动惯量的修正Timoshenko梁对振型峰值的影响越大。所以在工程结构抗震计算中，按照修正Timoshenko梁设计是偏于安全的。

图3 黏弹性Pasternak地基上经典Timoshenko梁与修正Timoshenko梁的前6阶模态（l=30 m）

Fig.3 The first six order modes of classical Timoshenko beam and modified Timoshenko beam on viscoelastic Pasternak foundation（l=30 m）

图4 黏弹性Pasternak地基上经典Timoshenko梁与修正Timoshenko梁的前6阶模态（l=6 m）

Fig.4 The first six order modes of classical Timoshenko beam and modified Timoshenko beam on viscoelastic Pasternak foundation（l=6 m）

3.2 边界条件对黏弹性Pasternak地基上修正Timoshenko梁自振特性的影响

黏弹性Pasternak地基上修正Timoshenko梁的各项物理参数采用表1和2中的数值，通过回传射线矩阵法，结合二分法和黄金分割法，计算了黏弹性Pasternak地基上修正Timoshenko梁在不同边界条件下前8阶自振频率和衰减系数的数值解如表5所示。

表5 不同边界条件下黏弹性Pasternak地基上修正Timoshenko梁的前8阶自振频率和衰减系数

Tab.5 The first eight order natural frequencies and attenuation coefficients of modified Timoshenko beam on viscoelastic Pasternak foundation under different boundary conditions

阶数	自由-自由		固定-自由		固定-固定		固定-简支		简支-简支
阶数	${\bar{ω}}_{n}$	$δ_{n}$	${\bar{ω}}_{n}$	$δ_{n}$	${\bar{ω}}_{n}$	$δ_{n}$	${\bar{ω}}_{n}$	$δ_{n}$	${\bar{ω}}_{n}$	$δ_{n}$
1	13.0102	0.82308	12.9521	0.82351	12.9432	0.82368	12.9325	0.82341	12.8036	0.82305
2	895.9527	0.70329	816.2489	0.78087	794.0807	0.77087	608.6540	0.77750	432.6988	0.72604
3	1954.0449	0.57165	1840.2387	0.72457	1723.0747	0.74967	1574.4856	0.73863	1407.1950	0.59576
4	3083.8719	0.51572	2918.6867	0.68364	2768.2405	0.71916	2653.7663	0.70508	2531.5759	0.52051
5	4190.3426	0.46777	4015.2194	0.65482	3839.4595	0.69345	3756.7773	0.67992	3668.6674	0.48086
6	5282.7712	0.44280	5100.7618	0.63782	4921.8337	0.67384	4857.7592	0.66270	4791.2676	0.45877
7	6355.1587	0.43669	6180.0322	0.62739	6001.2015	0.66009	5950.8560	0.65114	5898.1376	0.44552
8	7420.6872	0.41373	7249.4563	0.62143	7077.9599	0.65041	7035.9882	0.64326	6992.4382	0.43706

由表5可知，五种边界条件下黏弹性Pasternak地基上修正Timoshenko梁的前8阶自振频率由大到小依次排序为：自由‑自由>固定‑自由>固定‑固定>固定‑简支>简支‑简支；衰减系数由大到小依次排序为：固定‑固定>固定‑简支>固定‑自由>简支‑简支>自由‑自由，表明约束越强，振动能量衰减越明显。

图5为五种边界条件下黏弹性Pasternak地基上修正Timoshenko梁的前6阶模态。由图5可见，将各节点处的位移归一化处理后，其振型峰值变化并没有统一的规律，但不同边界条件下的振型曲线差异十分明显。

图5 五种边界条件下黏弹性Pasternak地基上修正Timoshenko梁的前6阶模态

Fig.5 The first six order modes of modified Timoshenko beam on viscoelastic Pasternak foundation under five boundary conditions

4 结论

将回传射线矩阵法推广至黏弹性地基梁的振动分析中，求解了黏弹性Pasternak地基中两端简支修正Timoshenko梁自振频率和衰减系数的解析解，对比分析了考虑剪切变形引起的转动惯量、梁长和不同的边界条件对结构自振特性的影响。得出以下结论：

（1）黏弹性Pasternak地基上修正Timoshenko梁的各阶自振频率和衰减系数小于经典Timoshenko梁的各阶自振频率和衰减系数；梁越短，对黏弹性Pasternak地基上修正Timoshenko梁自振频率和衰减系数的影响就越大，所以在黏弹性地基梁具体工程计算中涉及振动能量损耗控制时，或在梁对高跨比比较大的高频段计算结果有重大影响的分析中（诸如冲击等问题），应采用考虑剪切变形引起转动惯量的修正Timoshenko梁理论，否则会带来较大偏差。

（2）梁越短，考虑剪切变形引起转动惯量的修正Timoshenko梁对振型峰值的影响越大。所以在工程结构抗震计算中，按照修正Timoshenko梁设计是偏于安全的。

（3）不同边界条件下黏弹性Pasternak地基上修正Timoshenko梁的前8阶自振频率大小关系为：自由‑自由>固定‑自由>固定‑固定>固定‑简支>简支‑简支；其衰减系数大小关系为：固定‑固定>固定‑简支>固定‑自由>简支‑简支>自由‑自由，表明约束越强，振动能量衰减越明显。

（4）运用回传射线矩阵法可以求解埋置结构的高阶自振频率、模态等动力参数，列式统一，求解稳定，易于编程，且具有较高的计算效率和精度，从而为工程实践及工程设计提供理论基础。

参考文献

白海峰. 基于连续弹性地基梁的轨枕静力响应研究［J］. 铁道工程学报， 2007，24（5）： 22-27. [百度学术]

BAI Haifeng. Research on the static response of tie on continuous elastic grade beam［J］. Journal of Railway Engineering Society， 2007， 24（5）： 22-27． [百度学术]

刘飞禹，吴杰杰，陈江，等. 考虑水平摩阻的Pasternak路基基层及路面变形分析［J］. 中国公路学报， 2019， 32（5）： 38-46. [百度学术]

LIU Feiyu， WU Jiejie， CHEN Jiang， et al. Analysis of pavement and subbase on Pasternak foundation considering horizontal friction［J］. China Journal of Highway and Transport， 2019， 32（5）： 38-46. [百度学术]

Calim F F. Dynamic analysis of beams on viscoelastic foundation［J］. European Journal of Mechanics， A： Solid， 2009， 28（3）： 469-476． [百度学术]

Attarnejad R， Shahba A， Semnani S J. Application of differential transform in free vibration analysis of Timoshenko beams resting on two-parameter elastic foundation［J］. Arabian Journal for Science and Engineering， 2010， 35（2B）： 125-132. [百度学术]

Ding H， Chen L Q， Yang S P. Convergence of Galerkin truncation for dynamic response of finite beams on nonlinear foundations under a moving load［J］. Journal of Sound and Vibration， 2012， 331（10）： 2426-2442． [百度学术]

Malekzadeh P， Karami G， Farid M. DQEM for free vibration analysis of Timoshenko beams on elastic foundations［J］. Computational Mechanics， 2003， 31（3-4）： 219-228． [百度学术]

彭丽，丁虎，陈立群. 黏弹性三参数地基梁横向自由振动［J］. 振动与冲击， 2014， 33（1）： 101-105. [百度学术]

PENG Li， DING Hu， CHEN Liqun. Transverse free vibration of a beam resting on a three-parameter viscoelastic foundation［J］. Journal of Vibration and Shock， 2014， 33（1）： 101-105. [百度学术]

彭丽，丁虎，陈立群. 黏弹性Pasternak地基梁振动的复模态分析［J］. 振动与冲击， 2013， 32（2）： 143-146． [百度学术]

PENG Li， DING Hu， CHEN Liqun. Complex modal analysis for vibrations of a beam on a viscoelastic Pasternak foundation［J］. Journal of Vibration and Shock， 2013， 32（2）： 143-146． [百度学术]

彭丽，丁虎，陈立群.黏弹性三参数地基上Timoshenko梁横向自由振动［J］. 噪声与振动控制， 2013， 33（5）： 107-110． [百度学术]

PENG Li， DING Hu， CHEN Liqun. Transverse free vibration of a Timoshenko beam rested on transverse free vibration of a Timoshenko beam rested on three-parameter viscoelastic foundation［J］.Noise and Vibration Control， 2013， 33（5）：107-110． [百度学术]

余云燕，付艳艳，张伟. 黏弹性Pasternak地基上两跨连续Timoshenko梁横向自振特性分析［J］. 振动与冲击， 2023， 42（11）： 1-10. [百度学术]

YU Yunyan， FU Yanyan， ZHANG Wei. Analysis of transverse natural frequency of two-span continuous Timoshenko beam on viscoelastic Pasternak foundation［J］. Journal of Vibration and Shock， 2023， 42（11）： 1-10. [百度学术]

陈镕，万春风，薛松涛，等. Timoshenko梁运动方程的修正及其影响［J］. 同济大学学报， 2005， 33（6）： 711-715. [百度学术]

CHEN Rong， WAN Chunfeng， XUE Songtao， et al. Modification of motion equation of Timoshenko beam and its effect［J］. Journal of Tongji University， 2005， 33（6）： 711-715. [百度学术]

夏桂云. Timoshenko梁的第二频谱分析［J］. 湖南大学学报， 2021， 48（11）： 142-149. [百度学术]

XIA Guiyun. Analysis on the second frequency spectrum of Timoshenko beam［J］. Journal of Hunan University， 2021， 48（11）： 142-149. [百度学术]

吴晓，孙晋，黄翀，等. 用Timoshenko梁修正理论研究泡沫金属铝合金梁的动力响应［J］. 振动与冲击， 2011， 30（1）： 124-127. [百度学术]

WU Xiao， SUN Jin， HUANG Chong， et al. Dynamic responses of a metal aluminum alloy beam with Timoshenko beam theory［J］. Journal of Vibration and Shock， 2011， 30（1）： 124-127. [百度学术]

王家乐，夏桂云. Winkler地基上修正Timoshenko梁振动分析［J］. 振动与冲击， 2020， 39（3）： 30-37. [百度学术]

WANG Jiale， XIA Guiyun. Vibration analysis for a modified Timoshenko beam on Winkler elastic foundation［J］. Journal of Vibration and Shock， 2020， 39（3）： 30-37. [百度学术]

Li M L ， Wei P J ， Zhou X L. Wave propagation and free vibration of a Timoshenko beam mounted on the viscoelastic Pasternak foundation modeled by fraction-order derivatives［J］. Mechanics of Time-Dependent Materials， 2023， 27（4）： 1209-1223. [百度学术]

徐梅玲，叶茂，付明科，等. 修正Timoshenko 梁自由振动及Euler 梁误差分析［J］. 科学技术与工程， 2015， 15（15）： 88-94. [百度学术]

XU Meiling， YE Mao， FU Mingke， et al. Systematic modal analysis of simply supported beam based on the modified motion equation of Timoshenko［J］. Science Technology and Engineering， 2015， 15（15）： 88-94. [百度学术]

Kargarnovin M H， Younesian D， Thompson D J， et al. Response of beams on nonlinear viscoelastic foundations to harmonic moving loads［J］. Computers and Structures， 2005， 83（23-24）： 1865-1877. [百度学术]