Überblick
机械振动学习笔记
第1章 概论
1.1 机械振动研究的基本问题
研究振动的三要素
系统:研究对象(机械或结构)
激励:初始干扰和外界对振动系统的作用
响应:振动系统在激励作用下的动态行为
机械振动研究的三个基本问题
根据已知对象的不同的不同研究,可以分为如下三种:
振动分析:
系统识别:
载荷识别/环境预测:
以上三类问题时机械振动研究的基本问题
1.2 机械振动系统的力学模型与数学模型
1.3 机械振动的分类
时变系统、时不变系统
自由振动、强迫振动/受迫振动
无阻尼振动、阻尼振动
简谐振动、非简谐周期振动、非周期性振动
确定性振动、随机振动
线振动、角振动
usw.
第2章 机械振动系统
2.1 构成力学模型的基本原件
惯性元件
在运动时产生与(角)加速度成线性关系的惯性力(矩)的元件。用参振质量来描述它的特性,简称(广义)质量。
一个刚体在运动时的动能由平动动能和转动动能组成。计算公式为:
E
k
=
1
2
m
v
c
2
+
1
2
J
c
ω
2
E_k={1\over2}mv^2_c+{1\over2}J_c\omega^2
Ek=21mvc2+21Jcω2
其中
v
c
v_c
vc为质心C的运动速率,
J
c
J_c
Jc为刚体关于质心C的转动惯量
弹性元件
元件在变形时将产生与变形相关、抵抗变形的弹性恢复力(矩)的元件。描述线性弹性元件特性的参数通常称为刚度系数,简称刚度。
单个线性弹性元件变形时储存的势能为:
E
p
=
1
2
k
Δ
2
E_p={1\over2}k\Delta^2
Ep=21kΔ2
其中
Δ
\Delta
Δ为线位移或角位移,当其为角位移时,通常将系数
k
k
k表示为
k
ϕ
k_\phi
kϕ
阻尼元件
元件在变形时将产生与变形速度相关、阻碍变形变化变化的阻尼力(矩)的元件。黏性阻尼元件就是线性阻尼元件。描述线性阻尼元件特性的参数通常称为阻尼系数,简称阻尼。
阻尼元件在变形变化时消耗机械能功率的一半称为散逸函数。单个线性阻尼元件变形变化时对应的散逸函数为:
E
d
=
1
2
c
v
2
E_d={1\over2}cv^2
Ed=21cv2
类似的,当
v
v
v表示角速度时,
c
c
c写作
c
ϕ
c_\phi
cϕ
2.2 连接、约束和激励
2.3 建立力学模型的基本原则
等效性、简易性和逐步逼近三个原则。
一个振动系统的自由度是指确定系统运动所需要的独立坐标的数目。由离散元件或集中参数元件构造的系统为有限自由度的系统,其自由度通常决定于被弹性元件隔离的惯性元件的数目。
2.3 机械振动系统力学模型的建立
第3章 单自由度系统的振动理论
3.1 单自由度振动系统的运动方程
一个单自由度系统。质量经由弹簧和阻尼器同时和基座相连。那么我们就可以得到他的运动方程:
泛定方程:
m
x
¨
(
t
)
+
c
x
˙
(
t
)
+
k
x
(
t
)
=
f
(
t
)
(3-1)
m\ddot x(t)+c\dot x(t)+kx(t)=f(t)\tag{3-1}
mx¨(t)+cx˙(t)+kx(t)=f(t)(3-1)
定解方程:
{
m
x
¨
(
t
)
+
c
x
˙
(
t
)
+
k
x
(
t
)
=
f
(
t
)
x
(
0
)
=
x
0
,
x
˙
(
0
)
=
x
˙
0
(3-2)
\begin{cases}m\ddot x(t)+c\dot x(t)+kx(t)=f(t)\\x(0)=x_0, \dot x(0)=\dot x_0\end{cases}\tag{3-2}
{mx¨(t)+cx˙(t)+kx(t)=f(t)x(0)=x0,x˙(0)=x˙0(3-2)
在列运动方程时,除了使用牛顿定律,还可以使用功能原理。
KaTeX parse error: \tag works only in display equations
如果我们把动能,散逸函数,势能分别表示成如下的形式:
E
k
=
1
2
M
x
˙
2
E_k={1\over 2}M \dot x^2
Ek=21Mx˙2
E
d
=
1
2
C
x
˙
2
E_d={1\over 2}C\dot x^2
Ed=21Cx˙2
E
p
=
1
2
K
x
2
+
f
0
x
E_p={1\over 2}Kx^2+f_0x
Ep=21Kx2+f0x
由此可以得到泛定方程的一般形式:
M
x
¨
(
t
)
+
C
x
˙
(
t
)
+
K
x
(
t
)
=
p
(
t
)
(3-4)
M\ddot x(t)+C\dot x(t)+Kx(t)=p(t)\tag{3-4}
Mx¨(t)+Cx˙(t)+Kx(t)=p(t)(3-4)
其中
p
(
t
)
=
−
f
0
+
d
W
d
t
p(t)=-f_0+{\mathrm{d}W\over\mathrm{d}t}
p(t)=−f0+dtdW。如果
p
(
t
)
=
0
p(t)=0
p(t)=0则对应的振动称为自由振动,否则就称为受迫振动。
【例3-1】
取x为广义坐标,则其系统动能和势能分别为:
E
k
=
1
2
m
x
˙
2
+
1
2
J
θ
˙
2
=
1
2
(
m
+
J
R
2
)
⋅
x
˙
2
E_k={1\over 2}m\dot x^2+{1\over 2}J \dot \theta^2={1\over 2}(m+{J\over R^2})\cdot \dot x^2
Ek=21mx˙2+21Jθ˙2=21(m+R2J)⋅x˙2
E
p
=
1
2
k
x
2
−
m
g
x
E_p={1\over 2}kx^2-mgx
Ep=21kx2−mgx
由此我们可以得到系统的泛定方程
(
m
+
J
R
2
)
x
¨
2
+
k
x
(
t
)
=
m
g
(m+{J\over R^2})\ddot x^2+kx(t)=mg
(m+R2J)x¨2+kx(t)=mg
3.2 单自由度系统的自由振动
下面我们研究以下齐次线性微分方程的定解问题:
{
M
x
¨
(
t
)
+
C
x
˙
(
t
)
+
K
x
(
t
)
=
0
x
(
0
)
=
x
0
,
x
˙
(
0
)
=
x
˙
0
(3-5)
\begin{cases}M\ddot x(t)+C\dot x(t)+Kx(t)=0\\x(0)=x_0, \dot x(0)=\dot x_0\end{cases}\tag{3-5}
{Mx¨(t)+Cx˙(t)+Kx(t)=0x(0)=x0,x˙(0)=x˙0(3-5)
3.2.1 无阻尼自由振动
研究的方程为:
M
x
¨
(
t
)
+
K
x
(
t
)
=
0
(3-6)
M \ddot x(t)+Kx(t)=0\tag{3-6}
Mx¨(t)+Kx(t)=0(3-6)
试凑法
我们假设这个方程的非零解为:
x
(
t
)
=
x
u
s
i
n
(
ω
n
t
+
θ
)
(3-7)
x(t)=x_u\mathrm{sin}(\omega_nt+\theta)\tag{3-7}
x(t)=xusin(ωnt+θ)(3-7)
我们将(3-7)代入方程(3-6)可得
(
−
M
ω
n
2
+
K
)
x
u
s
i
n
(
ω
n
t
+
θ
)
=
0
(3-8)
(-M\omega^2_n+K)x_u\mathrm{sin}(\omega_nt+\theta)=0\tag{3-8}
(−Mωn2+K)xusin(ωnt+θ)=0(3-8)
显然
(
−
M
ω
n
2
+
K
)
x
u
(3-9)
(-M\omega^2_n+K)x_u\tag{3-9}
(−Mωn2+K)xu(3-9)
对于非零解,有
x
u
≠
0
x_u\not=0
xu=0,因此我们可以得到系统自由振动的频率必须满足
ω
n
=
K
M
(3-10)
\omega_n=\sqrt{K\over M}\tag{3-10}
ωn=MK(3-10)
ω
n
\omega_n
ωn也被称为固有(圆)频率。
如果我们把初始条件代入他的通解当中。那么我们就可以得到对应的系统自由振动的特解:
x
(
t
)
=
x
0
c
o
s
ω
n
t
+
x
˙
0
ω
n
s
i
n
ω
n
t
=
x
u
s
i
n
(
ω
n
t
+
θ
)
(3-11)
x(t)=x_0\mathrm{cos}\omega_nt+{\dot x_0\over \omega_n}\mathrm{sin}\omega_nt=x_u\mathrm{sin}(\omega_nt+\theta)\tag{3-11}
x(t)=x0cosωnt+ωnx˙0sinωnt=xusin(ωnt+θ)(3-11)
进一步,我们可以得到自由振动的振幅和相位角的计算公式:
x
u
=
x
0
2
+
(
x
˙
0
ω
n
)
2
x_u=\sqrt{x^2_0+({\dot x_0 \over \omega_n})^2}
xu=x02+(ωnx˙0)2
θ
=
t
g
−
1
x
0
ω
n
x
˙
0
\theta=\mathrm{tg^{-1}}{x_0\omega_n\over \dot x_0}
θ=tg−1x˙0x0ωn
圆频率
ω
n
\omega_n
ωn,固有频率
f
n
f_n
fn以及周期
T
n
T_n
Tn可以满足以下关系:
f
n
=
1
T
n
(3-12)
f_n={1 \over T_n}\tag{3-12}
fn=Tn1(3-12)
f
n
=
ω
n
2
π
(3-13)
f_n={\omega_n\over 2\pi}\tag{3-13}
fn=2πωn(3-13)
T
n
=
2
π
ω
n
(3-14)
T_n={2 \pi \over \omega _n}\tag{3-14}
Tn=ωn2π(3-14)
特征方程法
泛定方程(3-6)的特征方程为:
M
λ
2
+
K
=
0
(3-15)
M\lambda^2+K=0\tag{3-15}
Mλ2+K=0(3-15)
由此解的原方程的特征值为:
λ
1
,
2
=
±
j
ω
n
(3-16)
\lambda_{1, 2}=\pm \mathrm{j}\omega_n\tag{3-16}
λ1,2=±jωn(3-16)
由此该方程的通解为:
x
(
t
)
=
c
1
e
λ
1
t
+
c
2
e
λ
2
t
=
c
1
e
j
ω
n
t
+
c
2
e
−
j
ω
n
t
(3-17)
x(t)=c_1\mathrm{e}^{\lambda_1t}+c_2\mathrm{e}^{\lambda_2t}=c_1\mathrm{e}^{\mathrm{j}\omega_nt}+c_2\mathrm{e}^{-\mathrm{j}\omega_nt}\tag{3-17}
x(t)=c1eλ1t+c2eλ2t=c1ejωnt+c2e−jωnt(3-17)
那么根据欧拉公式
e
j
x
=
c
o
s
x
+
j
s
i
n
x
\mathrm{e}^{jx}=\mathrm{cos}x+j\mathrm{sin}x
ejx=cosx+jsinx可以得到,当
c
1
c_1
c1和
c
2
c_2
c2为共轭复数时,该方程有实数解。此时可以将式(3-17)表示为:
x
(
t
)
=
x
u
s
i
n
(
ω
n
t
+
θ
)
(3-18)
x(t)=x_u\mathrm{sin}(\omega_nt+\theta)\tag{3-18}
x(t)=xusin(ωnt+θ)(3-18)
3.2.2 有阻尼自由振动
有阻尼自由振动泛定方程的一般形式为:
M
x
¨
(
t
)
+
C
x
˙
(
t
)
+
K
x
(
t
)
=
0
(3-20)
M \ddot x (t)+C \dot x(t)+K x(t) = 0 \tag{3-20}
Mx¨(t)+Cx˙(t)+Kx(t)=0(3-20)
那么他的特征方程就是:
M
λ
2
+
C
λ
+
K
=
0
(3-21)
M \lambda^2+C \lambda + K = 0\tag{3-21}
Mλ2+Cλ+K=0(3-21)
由此可以解得特征方程的根(即特征值)为:
λ
1
,
2
=
−
C
2
K
±
C
2
−
4
M
K
2
M
(3-22)
\lambda_{1, 2}=-{C \over 2K} \pm {\sqrt{C^2-4MK}\over 2M}\tag{3-22}
λ1,2=−2KC±2MC2−4MK(3-22)
根据式(3-22)中的根式为实数、零还是虚数,特征值将得到不同的情况,对应
c
c
2
−
4
M
K
=
0
c_c^2-4MK=0
cc2−4MK=0时的阻尼称为临界阻尼,有:
c
c
=
2
M
K
(3-23)
c_c=2\sqrt{MK}\tag{3-23}
cc=2MK(3-23)
从而式(3-22)化为:
λ
1
,
2
=
(
−
ζ
±
ζ
2
−
1
)
ω
n
(3-24)
\lambda_{1,2}=(-\zeta \pm \sqrt{\zeta^2-1})\omega_n\tag{3-24}
λ1,2=(−ζ±ζ2−1)ωn(3-24)
其中,
ω
n
=
K
M
\omega_n=\sqrt{K\over M}
ωn=MK,为
C
=
0
C=0
C=0时的固有圆频率,称为系统的无阻尼固有圆频率;
ζ
=
C
c
c
=
C
2
M
K
\zeta={C\over c_c}={C \over 2 \sqrt{MK}}
ζ=ccC=2MKC,为系统实际阻尼与临界阻尼的比值,称为系统的阻尼比。
式(3-24)表明,特征值的性质取决于
ζ
\zeta
ζ值的大小。即系统的自由振动或微分方程的通解的特性决定于
ζ
\zeta
ζ值的大小。我们可以将有阻尼系统进行如下的分类:
1. 欠阻尼系统
阻尼比
ζ
<
1
\zeta <1
ζ<1(即实际阻尼小于临界阻尼)的系统称为欠阻尼系统或弱阻尼系统。根据式(3-24),欠阻尼系统的特征值为一对共轭复数:
λ
1
,
2
=
(
−
ζ
±
j
1
−
ζ
2
)
ω
n
(3-25)
\lambda_{1,2}=(-\zeta \pm \mathrm{j} \sqrt{1-\zeta^2})\omega_n \tag{3-25}
λ1,2=(−ζ±j1−ζ2)ωn(3-25)
那么其对应的微分方程方程的通解为:
x
(
t
)
=
c
1
e
(
−
ζ
+
j
1
−
ζ
2
)
ω
n
+
c
2
e
(
−
ζ
−
j
1
−
ζ
2
)
ω
n
=
e
−
ζ
ω
n
t
(
c
1
e
j
1
−
ζ
2
ω
n
t
+
c
2
e
−
j
1
−
ζ
2
ω
n
t
)
(3-26)
\begin{aligned} x(t)&=c_1 \mathrm{e}^{(-\zeta + \mathrm{j} \sqrt{1-\zeta^2})\omega_n}+c_2\mathrm{e}^{(-\zeta - \mathrm{j} \sqrt{1-\zeta^2})\omega_n} \\&=\mathrm{e}^{-\zeta \omega_n t}(c_1\mathrm{e}^{\mathrm{j}\sqrt{1-\zeta^2}\omega_n t}+c_2\mathrm{e}^{-\mathrm{j}\sqrt{1-\zeta^2}\omega_n t}) \end{aligned} \tag{3-26}
x(t)=c1e(−ζ+j1−ζ2)ωn+c2e(−ζ−j1−ζ2)ωn=e−ζωnt(c1ej1−ζ2ωnt+c2e−j1−ζ2ωnt)(3-26)
其中
c
1
,
c
2
c_1, c_2
c1,c2为任意常数,在他们共轭时,式(3-26)对应为实数解,可以表示为以下形式:
x
(
t
)
=
e
−
ζ
ω
n
t
(
a
c
o
s
ω
d
t
+
b
s
i
n
ω
d
t
)
=
A
e
−
ζ
ω
n
t
s
i
n
(
ω
d
+
θ
)
(3-27)
\begin{aligned} x(t)&=\mathrm{e}^{-\zeta \omega_n t}(a\mathrm{cos}\omega_dt+b\mathrm{sin}\omega_dt) \\&=A \mathrm{e}^{-\zeta \omega_n t}\mathrm{sin}(\omega_d+\theta) \end{aligned} \tag{3-27}
x(t)=e−ζωnt(acosωdt+bsinωdt)=Ae−ζωntsin(ωd+θ)(3-27)
其中,
a
,
b
a,b
a,b或
A
,
θ
A,\theta
A,θ为通解的任意实常数,由初始条件确定;
ω
d
=
1
−
ζ
2
ω
n
\omega_d=\sqrt{1-\zeta^2}\omega_n
ωd=1−ζ2ωn为系统的固有圆频率。我们一般更有区分的,更准确地称之为系统的有阻尼固有圆频率。对于小阻尼
(
ζ
≪
1
)
(\zeta \ll 1)
(ζ≪1)振动系统,
ω
d
≈
ω
n
\omega_d \approx \omega_n
ωd≈ωn。
如果我们把振动方程的初始条件代入式(3-27)我们就可以得到对应的自由整栋特解为:
x
(
t
)
=
A
e
−
ζ
ω
n
t
s
i
n
(
ω
d
+
θ
)
=
e
−
ζ
ω
n
t
(
x
0
c
o
s
ω
d
t
+
x
˙
0
+
ζ
ω
n
x
0
ω
d
s
i
n
ω
d
t
)
(3-27)
\begin{aligned} x(t)&=A \mathrm{e}^{-\zeta \omega_n t}\mathrm{sin}(\omega_d+\theta) \\&=\mathrm{e}^{-\zeta \omega_n t}(x_0\mathrm{cos}\omega_dt+{\dot x_0 + \zeta \omega_n x_0 \over \omega_d}\mathrm{sin}\omega_dt) \end{aligned} \tag{3-27}
x(t)=Ae−ζωntsin(ωd+θ)=e−ζωnt(x0cosωdt+ωdx˙0+ζωnx0sinωdt)(3-27)
其中,
A
=
(
x
˙
0
+
ζ
ω
n
x
0
ω
d
)
2
+
x
0
2
θ
=
t
g
−
1
x
0
ω
d
x
˙
0
+
ζ
ω
n
x
0
\begin{aligned} A&=\sqrt{({\dot x_0 + \zeta \omega_n x_0 \over \omega_d})^2+x_0^2} \\ \theta&=\mathrm{tg^{-1}}{x_0 \omega_d \over \dot x_0 + \zeta \omega_n x_0} \end{aligned}
Aθ=(ωdx˙0+ζωnx0)2+x02=tg−1x˙0+ζωnx0x0ωd
由式(3-28)可知,欠阻尼系统自由振动的振幅可以看作是一随时间变化函数:
x
u
(
t
)
=
A
e
−
ζ
ω
n
t
(3-29)
x_u(t)=A\mathrm{e}^{-\zeta \omega_n t}\tag{3-29}
xu(t)=Ae−ζωnt(3-29)
因此其相邻两个振幅之间的比值为:
x
u
(
t
)
x
u
(
t
+
T
d
)
=
e
−
ζ
ω
n
t
e
−
ζ
ω
n
(
t
+
T
d
)
=
e
ζ
ω
n
T
d
=
e
x
p
(
2
π
ζ
1
−
ζ
2
)
(3-30)
{x_u(t)\over x_u(t+T_d)}={\mathrm{e}^{-\zeta \omega_n t}\over \mathrm{e}^{-\zeta \omega_n (t +T_d)}}=\mathrm{e}^{\zeta \omega_n T_d}=\mathrm{exp}({2\pi \zeta \over \sqrt{1-\zeta^2}})\tag{3-30}
xu(t+Td)xu(t)=e−ζωn(t+Td)e−ζωnt=eζωnTd=exp(1−ζ22πζ)(3-30)
考虑到欠阻尼系统自由振动的相邻两个振幅之间的比值为常数,且振幅比值只与系统阻尼比有关,我们称之为衰减比。为了方便起见,将衰减比的自然对数定义为对数衰减率
δ
\delta
δ。即:
δ
=
l
n
x
u
(
t
)
x
u
(
t
+
T
d
)
=
2
π
ζ
1
−
ζ
2
(3-31)
\delta = \mathrm{ln}{x_u(t)\over x_u(t+T_d)} = {2\pi \zeta \over \sqrt{1-\zeta^2}}\tag{3-31}
δ=lnxu(t+Td)xu(t)=1−ζ22πζ(3-31)
其中,
T
d
T_d
Td代表系统的有阻尼固有周期。
T
d
=
2
π
ω
d
T_d={2 \pi \over \omega_d}
Td=ωd2π。当相邻
n
n
n个振幅时,易证:
δ
=
1
n
l
n
x
u
(
t
)
x
u
(
t
+
n
T
d
)
=
2
π
ζ
1
−
ζ
2
(3-32)
\delta = {1\over n}\mathrm{ln}{x_u(t)\over x_u(t+nT_d)} = {2\pi \zeta \over \sqrt{1-\zeta^2}}\tag{3-32}
δ=n1lnxu(t+nTd)xu(t)=1−ζ22πζ(3-32)
在阻尼较小时,式(3-32)可近似为:
δ
≈
2
π
ζ
(3-33)
\delta \approx 2 \pi \zeta \tag{3-33}
δ≈2πζ(3-33)
由式(3-33)进一步可得,当阻尼比较小时,欠阻尼系统的阻尼比可表示为:
ζ
≈
δ
2
π
(3-34)
\zeta \approx {\delta \over 2 \pi} \tag{3-34}
ζ≈2πδ(3-34)
2. 临界阻尼系统
阻尼比
ζ
=
1
\zeta = 1
ζ=1(即实际阻尼等于临界阻尼)的系统称为临界阻尼系统。微分方程的通解为:
x
(
t
)
=
(
c
1
+
c
2
t
)
e
−
ω
n
t
(3-35)
x(t)=(c_1+c_2t)\mathrm{e}^{-\omega_nt}\tag{3-35}
x(t)=(c1+c2t)e−ωnt(3-35)
3. 过阻尼系统
阻尼比
ζ
>
1
\zeta>1
ζ>1的系统称为过阻尼系统或强阻尼系统。过阻尼系统的特征值为两个不相等的负实数:
λ
1
,
2
=
(
−
ζ
±
ζ
2
−
1
)
ω
n
(3-36)
\lambda_{1,2}=(-\zeta \pm \sqrt{\zeta^2-1})\omega_n\tag{3-36}
λ1,2=(−ζ±ζ2−1)ωn(3-36)
对应的通解为:
x
(
t
)
=
c
1
e
(
−
ζ
+
ζ
2
−
1
)
ω
n
t
+
c
2
e
(
−
ζ
−
ζ
2
−
1
)
ω
n
t
x(t)=c_1 \mathrm{e}^{(-\zeta+\sqrt{\zeta^2-1})\omega_n t}+c_2 \mathrm{e}^{(-\zeta-\sqrt{\zeta^2-1})\omega_n t}
x(t)=c1e(−ζ+ζ2−1)ωnt+c2e(−ζ−ζ2−1)ωnt
其中
c
1
,
c
2
c_1,c_2
c1,c2为任意实常数。
3.3 单自由度系统简谐激励下的强迫振动
3.3.1 强迫振动的解的结构
由前面分析可知,强迫振动(或受迫振动)对应于微分方程式(3=2)中
p
(
t
)
≠
0
p(t) \not = 0
p(t)=0,即对应以下非齐次的泛定方程:
M
x
¨
(
t
)
+
C
x
˙
(
t
)
+
K
x
(
t
)
=
p
(
t
)
(3-38)
M\ddot x(t)+C \dot x(t)+Kx(t)=p(t) \tag{3-38}
Mx¨(t)+Cx˙(t)+Kx(t)=p(t)(3-38)
上式该用系统的固有参数表示为:
x
¨
(
t
)
+
2
ζ
ω
n
x
˙
(
t
)
+
ω
n
2
x
(
t
)
=
p
(
t
)
M
(3-39)
\ddot x(t) + 2 \zeta \omega_n \dot x(t) + \omega^2_n x(t)={p(t)\over M} \tag{3-39}
x¨(t)+2ζωnx˙(t)+ωn2x(t)=Mp(t)(3-39)
强调一遍,这里的
ω
n
=
K
M
\omega_n=\sqrt{K\over M}
ωn=MK,这里的
ζ
=
C
2
M
K
\zeta={C\over 2 \sqrt{MK}}
ζ=2MKC。
设
x
1
(
t
)
x_1(t)
x1(t),
x
2
(
t
)
x_2(t)
x2(t)分别为式(3-38)或式(3-39)的两个不同的解,令:
x
h
(
t
)
=
x
1
(
t
)
−
x
2
(
t
)
x_h(t)=x_1(t)-x_2(t)
xh(t)=x1(t)−x2(t)
因为
x
1
(
t
)
x_1(t)
x1(t),
x
2
(
t
)
x_2(t)
x2(t)分别满足泛定方程,则可得:
M
x
¨
h
(
t
)
+
C
x
˙
h
(
t
)
+
K
x
h
(
t
)
=
0
M \ddot x_h(t) + C \dot x_h(t)+Kx_h(t)=0
Mx¨h(t)+Cx˙h(t)+Kxh(t)=0
因此可知,非齐次线性微分方程的通解具有以下形式:
x
(
t
)
=
x
s
(
t
)
+
x
h
(
t
)
(3-40)
x(t)=x_s(t)+x_h(t)\tag{3-40}
x(t)=xs(t)+xh(t)(3-40)
式中,
x
s
(
t
)
x_s(t)
xs(t)是满足非齐次方程是(3-38)的任意一个特解;
x
h
(
t
)
x_h(t)
xh(t)为非齐次方程对应齐次方程(即在原微分方程中令
p
(
t
)
=
0
p(t)=0
p(t)=0)的通解,即振动系统在没有外界激励时系统自由振动的通解。
由于对于有阻尼系统,
t
→
∞
t \to \infty
t→∞时,
x
h
(
t
)
→
0
x_h(t)\to 0
xh(t)→0。因此对于任意的初始条件,在经过足够长时间后,将趋向于相同的强迫振动。求解强迫振动的关键是求出其一个特解
x
s
(
t
)
x_s(t)
xs(t)。
下面先研究:
p
(
t
)
=
p
u
s
i
n
(
ω
T
+
θ
p
)
(3-41)
p(t)=p_u\mathrm{sin}(\omega T+\theta_p)\tag{3-41}
p(t)=pusin(ωT+θp)(3-41)
其中
ω
\omega
ω为与系统固有频率无关的激励圆频率。那么简谐激励强迫振动的泛定方程为
M
x
¨
(
t
)
+
C
x
˙
(
t
)
+
K
x
(
t
)
=
p
u
s
i
n
(
ω
t
+
θ
p
)
(3-42)
M \ddot x(t) + C \dot x(t) + Kx(t) = p_u\mathrm{sin}(\omega t+ \theta _p)\tag{3-42}
Mx¨(t)+Cx˙(t)+Kx(t)=pusin(ωt+θp)(3-42)
为了研究的方便,我们将简谐振动用复数来表示。
3.3.2 简谐振动的复数表示
根据Eular公式有
e
j
σ
=
c
o
s
α
+
j
s
i
n
α
e
−
j
σ
=
c
o
s
α
−
j
s
i
n
α
(3-43)
\begin{aligned} \mathrm{e}^{\mathrm{j}\sigma}&=\mathrm{cos}\alpha+\mathrm{jsin}\alpha \\ \mathrm{e}^{-\mathrm{j}\sigma}&=\mathrm{cos}\alpha-\mathrm{jsin}\alpha \end{aligned} \tag{3-43}
ejσe−jσ=cosα+jsinα=cosα−jsinα(3-43)
由此可得:
s
i
n
α
=
1
2
j
(
e
j
σ
−
e
−
j
σ
)
=
I
m
(
e
j
α
t
)
(3-44)
\mathrm{sin}\alpha={1\over 2 \mathrm{j}}(\mathrm{e}^{\mathrm{j}\sigma}-\mathrm{e}^{-\mathrm{j}\sigma})=\mathrm{Im}(\mathrm{e}^{\mathrm{j}\alpha t})\tag{3-44}
sinα=2j1(ejσ−e−jσ)=Im(ejαt)(3-44)
因此简谐振动可以表示为以下的复数形式:
x
(
t
)
=
x
u
s
i
n
(
ω
t
+
θ
)
=
I
m
[
x
u
e
j
(
ω
t
+
θ
)
]
=
I
m
[
x
‾
(
t
)
]
(3-45)
x(t)=x_u \mathrm{sin}(\omega t + \theta)=\mathrm{Im}[x_u\mathrm{e}^{\mathrm{j}(\omega t + \theta)}]=\mathrm{Im}[\overline x(t)]\tag{3-45}
x(t)=xusin(ωt+θ)=Im[xuej(ωt+θ)]=Im[x(t)](3-45)
这里的
x
‾
(
t
)
\overline x(t)
x(t)代表简谐函数
x
(
t
)
x(t)
x(t)的(单边)复数表示,为如下的复指数函数:
x
‾
(
t
)
=
x
u
e
j
(
ω
t
+
θ
)
=
x
u
e
j
θ
e
j
ω
t
=
X
e
j
ω
t
(3-46)
\overline x(t)=x_u\mathrm{e}^{\mathrm{j}(\omega t + \theta)}=x_u\mathrm{e}^{\mathrm{j}\theta}\mathrm{e}^{\mathrm{j}\omega t}=X\mathrm{e}^{\mathrm{j}\omega t}\tag{3-46}
x(t)=xuej(ωt+θ)=xuejθejωt=Xejωt(3-46)
其中,
X
X
X为复常数,称为
x
(
t
)
x(t)
x(t)的(单边)复振幅,表示为:
X
=
x
u
e
j
θ
=
a
+
j
b
(3-47)
X=x_u\mathrm{e}^{\mathrm{j}\theta}=a+\mathrm{j}b\tag{3-47}
X=xuejθ=a+jb(3-47)
显然,复振幅同时包含了简谐振动的振幅和相位信息。在已知振动频率的情况下,简谐振动与复振幅具有一一对应的关系,即
x
(
t
)
↔
X
(3-48)
x(t)\leftrightarrow X\tag{3-48}
x(t)↔X(3-48)
(复振幅
X
X
X的模就是原简谐振动的的振幅的大小,
X
X
X的幅角就是原简谐振动的相位的大小。)
进一步我们可以把两个频率相同单振幅、相位不同的
x
1
(
t
)
x_1(t)
x1(t),
x
2
(
t
)
x_2(t)
x2(t)简谐振动的线形符合,表示成他们的复振幅的线形组合的形式。
C
1
x
1
(
t
)
+
C
2
x
2
(
t
)
↔
C
1
X
1
+
C
2
X
2
(3-49)
C_1x_1(t)+C_2x_2(t)\leftrightarrow C_1X_1+C_2X_2\tag{3-49}
C1x1(t)+C2x2(t)↔C1X1+C2X2(3-49)
简谐振动的导数仍然是同频率的简谐函数,并且
有
x
˙
(
t
)
↔
j
ω
X
(3-50)
\dot x(t) \leftrightarrow \mathrm{j} \omega X \tag{3-50}
x˙(t)↔jωX(3-50)
由此我们可以导出同频简谐函数满足微分方程式时,对应复振幅应满足的对应方程。设同频简谐函数
x
(
t
)
x(t)
x(t),
p
(
t
)
p(t)
p(t)满足以下关系:
d
2
x
¨
(
t
)
+
d
1
x
˙
(
t
)
+
d
0
x
(
t
)
=
p
(
t
)
(3-51)
d_2 \ddot x(t) + d_1 \dot x(t) + d_0 x(t) = p(t) \tag{3-51}
d2x¨(t)+d1x˙(t)+d0x(t)=p(t)(3-51)
对应的复振幅
X
X
X,
P
P
P应该满足:
d
2
(
j
ω
)
2
X
+
d
1
ω
j
X
+
d
0
X
=
P
(3-52)
d_2(\mathrm{j}\omega)^2X+d_1\omega \mathrm{j} X+d_0X=P \tag{3-52}
d2(jω)2X+d1ωjX+d0X=P(3-52)
简化得:
G
(
ω
)
X
=
P
(3-53)
G(\omega)X=P\tag{3-53}
G(ω)X=P(3-53)
其中,
G
(
ω
)
=
(
−
d
2
ω
2
+
d
0
)
+
d
1
j
ω
G(\omega)=(-d_2\omega ^2+d_0)+d_1\mathrm{j}\omega
G(ω)=(−d2ω2+d0)+d1jω。
由此,在分析频率
ω
\omega
ω确定时,已知
x
(
t
)
x(t)
x(t)可以计算
p
(
t
)
p(t)
p(t)反之亦然。
3.3.3 简谐激励的稳态响应
根据微分方程理论,在简谐激励时,总是存在一个与激励频率相同的简谐振动特解。我们设这个特解为:
x
s
(
t
)
=
x
s
u
s
i
n
(
ω
t
+
θ
s
x
)
(3-54)
x_s(t)=x_{su}\mathrm{sin}(\omega t+\theta_{sx})\tag{3-54}
xs(t)=xsusin(ωt+θsx)(3-54)
在不至引起误解的情况下,我们省略
x
s
(
t
)
x_s(t)
xs(t)的下标
s
s
s,而将上式简记为:
x
(
t
)
=
x
u
s
i
n
(
ω
t
+
θ
x
)
(3-55)
x(t)=x_u\mathrm{sin}(\omega t+\theta_x)\tag{3-55}
x(t)=xusin(ωt+θx)(3-55)
因稳态振动的频率与激励频率相同,所以式(3-55)代入泛定方程(3-42)后,两边都是同频率的简谐函数。
方程(3-42)相应复数形式为
[
M
(
j
ω
)
2
+
C
j
ω
+
K
]
X
=
P
[M(\mathrm{j} \omega)^2+C\mathrm{j} \omega +K]X=P
[M(jω)2+Cjω+K]X=P
其中,
P
P
P,
X
X
X分别为简谐激励
p
(
t
)
p(t)
p(t)和系统的稳态响应x(t)的复振幅,我们可以把它改写为:
Z
(
ω
)
⋅
X
=
P
(3-56)
Z(\omega)\cdot X =P \tag{3-56}
Z(ω)⋅X=P(3-56)
其中,
Z
(
ω
)
Z(\omega)
Z(ω)称为系统阻抗,定义为
Z
(
ω
)
=
(
K
−
ω
2
M
)
+
j
ω
C
(3-57)
Z(\omega)=(K-\omega^2M)+\mathrm{j}\omega C \tag{3-57}
Z(ω)=(K−ω2M)+jωC(3-57)
由此得稳态响应的复振幅为
X
=
P
Z
(
ω
)
=
H
(
ω
)
⋅
P
(3-58)
X={P\over Z(\omega)}=H(\omega)\cdot P \tag{3-58}
X=Z(ω)P=H(ω)⋅P(3-58)
式中,
H
(
ω
)
H(\omega)
H(ω)称为系统的频率响应函数或频率响应特性,为:
H
(
ω
)
=
h
u
e
j
θ
h
=
1
Z
(
ω
)
=
1
(
K
−
ω
2
M
)
+
j
ω
C
(3-59)
H(\omega)=h_u \mathrm{e}^{\mathrm{j}\theta_h}={1\over Z(\omega)}={1\over (K-\omega^2M)+\mathrm{j}\omega C}\tag{3-59}
H(ω)=huejθh=Z(ω)1=(K−ω2M)+jωC1(3-59)
由此得系统稳态振动的振幅和相位角分别为:
x
u
=
h
u
⋅
p
u
(3-60)
x_u=h_u \cdot p_u\tag{3-60}
xu=hu⋅pu(3-60)
θ
x
=
θ
p
+
θ
h
(3-61)
\theta_x = \theta_p + \theta_h\tag{3-61}
θx=θp+θh(3-61)
式中,
h
u
=
h
u
(
ω
)
h_u=h_u(\omega)
hu=hu(ω)和
θ
h
=
θ
h
(
ω
)
\theta_h=\theta_h(\omega)
θh=θh(ω)分别为
H
(
ω
)
H(\omega)
H(ω)的幅值和相位角(称为幅频特性和相频特性/相频特性),分别由式(3-62)和式(3-63)计算:
h
u
=
1
(
K
−
ω
2
M
)
2
+
(
ω
C
)
2
(3-62)
h_u={1 \over \sqrt{(K-\omega^2M)^2+(\omega C)^2}}\tag{3-62}
hu=(K−ω2M)2+(ωC)21(3-62)
θ
h
=
−
t
g
−
1
ω
c
K
−
ω
2
M
(3-63)
\theta_h=-\mathrm{tg^{-1}}{\omega c \over K- \omega^2M}\tag{3-63}
θh=−tg−1K−ω2Mωc(3-63)
稳态响应的振幅还可以表示为:
x
u
=
h
u
⋅
p
u
=
R
d
y
x
s
t
(3-64)
x_u=h_u \cdot p_u = R_{dy}x_{st}\tag{3-64}
xu=hu⋅pu=Rdyxst(3-64)
其中
x
s
t
x_{st}
xst为系统受到大小为简谐激励幅值
p
u
p_u
pu的静态力时对应的系统静态位移;
R
d
y
R_{dy}
Rdy为稳态振动的振幅与静态位移的比值,称为动力放大系数。
R
d
y
=
x
u
x
s
t
=
h
u
K
=
k
(
K
−
ω
2
M
)
2
+
(
ω
C
)
2
(3-65)
R_{dy}={x_u \over x_{st}}=h_u K={k \over \sqrt{(K-\omega^2 M)^2 + (\omega C)^2}}\tag{3-65}
Rdy=xstxu=huK=(K−ω2M)2+(ωC)2k(3-65)
如果改用系统的固有参数代替式(3-59)(3-62)(3-63)(3-65)中的质量
m
m
m,刚度
k
k
k和阻尼
c
c
c,系统的频率响应函数及幅频特性与相频特性,动力放大系数可以分别表示为:
H
(
ω
)
=
1
K
⋅
ω
n
2
(
ω
n
2
−
ω
2
)
+
j
2
ζ
ω
ω
n
=
1
K
⋅
1
(
1
−
λ
2
)
+
j
2
ζ
λ
(3-66)
H(\omega)={1\over K}\cdot{\omega^2_n\over (\omega^2_n-\omega^2)+\mathrm{j}2\zeta\omega \omega_n}={1\over K}\cdot {1\over (1-\lambda^2)+\mathrm{j}2\zeta\lambda} \tag{3-66}
H(ω)=K1⋅(ωn2−ω2)+j2ζωωnωn2=K1⋅(1−λ2)+j2ζλ1(3-66)
h
u
=
1
K
⋅
1
(
1
−
λ
2
)
2
+
(
2
ζ
λ
)
2
(3-67)
h_u={1\over K}\cdot {1\over \sqrt{(1-\lambda^2)^2+(2\zeta\lambda)^2}} \tag{3-67}
hu=K1⋅(1−λ2)2+(2ζλ)21(3-67)
θ
h
=
−
t
g
−
1
2
ζ
λ
1
−
λ
2
(3-68)
\theta_h=-\mathrm{tg^{-1}}{2\zeta\lambda \over 1-\lambda^2} \tag{3-68}
θh=−tg−11−λ22ζλ(3-68)
R
d
y
=
1
(
1
−
λ
2
)
2
+
(
2
ζ
λ
)
2
(3-69)
R_{dy}={1\over \sqrt{(1-\lambda^2)^2+(2\zeta\lambda)^2}} \tag{3-69}
Rdy=(1−λ2)2+(2ζλ)21(3-69)
其中
ω
n
\omega_n
ωn为系统的无阻尼固有圆频率;
ζ
\zeta
ζ为系统阻尼比;
λ
=
ω
ω
n
\lambda={\omega \over \omega_n}
λ=ωnω为激励频率与系统无阻尼固有频率的比值,简称频率比。
研究单自由度同频率响应函数的幅频特性和相频特性可以得出如下这些结论:
(1)当阻尼比
ζ
<
2
2
\zeta<{\sqrt{2}\over 2}
ζ<22时,在系统固有频率附近将出现共振峰值;
(2)在系统固有频率(即
λ
≈
1
\lambda \approx 1
λ≈1)的共振区域(阻尼控制区),振幅的大小主要由系统阻尼比
ζ
\zeta
ζ的大小决定,激励力主要由阻尼力平衡,阻尼对共振峰值有明显的抑止作用。
(3)
λ
≪
1
\lambda \ll 1
λ≪1的区域(刚度控制区),振幅的大小主要由系统刚度
K
K
K的大小决定,激励力主要由系统弹性力平衡。
(4)
λ
≫
1
\lambda \gg 1
λ≫1的区域(质量控制区),振幅的大小主要由系统质量
M
M
M的大小决定,激励力主要由系统惯性力平衡。
(5)稳态振动的相位角总是迟后于激励的相位角,并且激励频率越高迟后的相位角就越大;
(6)
λ
=
1
\lambda =1
λ=1时,响应与激励的相位差
θ
=
−
π
2
\theta = -{\pi \over 2}
θ=−2π,而与阻尼比
ζ
\zeta
ζ无关,对应的实频特性曲线为0.
(7)当阻尼比
ζ
=
2
2
\zeta = {\sqrt{2}\over 2}
ζ=22时,两个特性曲线在低频处都最为平直,这种特性最有利于测量仪表测试结果的“波形不失真”。
求出系统的稳态振动后,进一步就可以求出系统运动方程式(3-42)的通解。对于欠阻尼系统,强迫振动的通解为:
x
(
t
)
=
x
s
(
t
)
+
x
h
(
t
)
=
R
d
y
x
s
t
c
o
s
(
ω
t
+
θ
f
+
θ
p
)
+
e
−
ζ
ω
n
t
(
a
c
o
s
ω
d
t
+
b
s
i
n
ω
d
t
)
(3-70)
\begin{aligned} x(t)&=x_s(t)+x_h(t) \\&=R_{dy}x_{st}\mathrm{cos}(\omega t+\theta_f+\theta_p)+\mathrm{e}^{-\zeta \omega_n t}(a\mathrm{cos}\omega_d t+b\mathrm{sin}\omega_d t) \end{aligned} \tag{3-70}
x(t)=xs(t)+xh(t)=Rdyxstcos(ωt+θf+θp)+e−ζωnt(acosωdt+bsinωdt)(3-70)
如果系统的初始条件已知时,就可以求出对应的强迫振动的特解为:
x
(
t
)
=
x
s
(
t
)
+
e
−
ζ
ω
n
t
[
(
x
0
−
x
s
0
)
c
o
s
ω
d
t
+
(
x
˙
0
−
x
˙
s
0
+
ζ
ω
n
(
x
0
−
x
s
0
)
)
ω
d
s
i
n
ω
d
t
]
=
x
s
(
t
)
+
e
−
ζ
ω
n
t
[
x
0
c
o
s
ω
d
t
+
x
˙
0
+
ζ
ω
n
x
0
ω
d
s
i
n
ω
d
t
]
−
e
−
ζ
ω
n
t
[
x
s
0
c
o
s
ω
d
t
+
x
˙
s
0
+
ζ
ω
n
x
s
0
ω
d
s
i
n
ω
d
t
]
(3-71)
\begin{aligned} x(t)&=x_s(t)+\mathrm{e}^{-\zeta \omega_n t}[(x_0-x_{s0})\mathrm{cos}\omega_d t+{(\dot x_0 -\dot x_{s0}+\zeta \omega_n (x_0 - x_{s0}))\over \omega_d}\mathrm{sin}\omega_d t] \\&=x_s(t)+\mathrm{e}^{-\zeta \omega_n t}[x_0 \mathrm{cos}\omega_d t+{\dot x_0+\zeta \omega_n x_0\over \omega_d}\mathrm{sin}\omega_d t]-\mathrm{e}^{-\zeta \omega_n t}[x_{s0} \mathrm{cos}\omega_d t+{\dot x_{s0}+\zeta \omega_n x_{s0}\over \omega_d}\mathrm{sin}\omega_d t] \end{aligned} \tag{3-71}
x(t)=xs(t)+e−ζωnt[(x0−xs0)cosωdt+ωd(x˙0−x˙s0+ζωn(x0−xs0))sinωdt]=xs(t)+e−ζωnt[x0cosωdt+ωdx˙0+ζωnx0sinωdt]−e−ζωnt[xs0cosωdt+ωdx˙s0+ζωnxs0sinωdt](3-71)
上式中的第三项是由于稳态振动的初始值不等于零所产生的,它伴随稳态振动而存在,所以称为伴随自由振动。
有时我们将系统稳态响应作为其强迫振动响应,而将考虑初始条件后得到的强迫振动的完整解式(3-71)称为瞬态响应。
3.4 等效阻尼
非线性阻尼求起来太麻烦了。所以我们可以通过等效阻尼的方法把它近似地线性化。具体原则是:一个非线性阻尼如果在一个周期内耗散的能量与一个线性阻尼耗散的能量相同时,则可认为这两个阻尼等效。
易求得等效黏性阻尼一个周期做功为:
W
e
q
=
π
c
e
q
ω
x
u
2
(3-72)
W_{eq}=\pi c_{eq} \omega x^2_u \tag{3-72}
Weq=πceqωxu2(3-72)
其中,
x
u
x_u
xuf为物体作简谐振动时的振幅。由此我们可以获得等效阻尼的计算公式:
c
e
q
=
W
π
ω
x
n
2
(3-73)
c_{eq}={W\over \pi \omega x^2_n}\tag{3-73}
ceq=πωxn2W(3-73)
下面研究一些常见非线性阻尼等效阻尼的计算方法。
干摩擦力
W
=
4
f
m
x
u
(3-74)
W=4f_mx_u\tag{3-74}
W=4fmxu(3-74)
由此等效黏性阻尼为:
c
e
q
=
4
f
m
π
ω
x
u
(3-75)
c_{eq}={4f_m\over \pi \omega x_u}\tag{3-75}
ceq=πωxu4fm(3-75)
紊流流体阻尼
阻尼力满足:
f
c
=
c
l
x
˙
2
f_c=c_l \dot x^2
fc=clx˙2。因此:
W
=
8
c
l
x
u
3
ω
2
3
(3-76)
W={8c_l x^3_u \omega^2\over 3}\tag{3-76}
W=38clxu3ω2(3-76)
由此等效黏性阻尼为:
c
e
q
=
8
c
l
ω
x
u
3
π
(3-77)
c_{eq}={8c_l \omega x_u \over 3 \pi}\tag{3-77}
ceq=3π8clωxu(3-77)
结构阻尼
实验表明:
W
=
a
x
u
2
(3-78)
W=ax_u^2\tag{3-78}
W=axu2(3-78)
由此等效黏性阻尼为:
c
e
q
=
a
π
ω
(3-79)
c_{eq}={a\over \pi \omega}\tag{3-79}
ceq=πωa(3-79)
3.5 非简谐激励下的强迫振动——频域分析法
对于非简谐激励,一般可分解为各个谐波成分的叠加。将激励各谐波对应的稳态响应叠加在一起就可以获得对应非简谐激励的稳态响应。
3.5.1周期激励的谐波分解
在复数形式的Fourier级数分解形式中,通常将谐波函数表示为余弦的形式。与式(3-45)对应的简谐函数复数表示形式为:
x
(
t
)
=
x
u
c
o
s
(
ω
t
+
θ
)
=
R
e
[
x
‾
(
t
)
]
=
R
e
(
X
e
j
ω
t
)
(3-80)
x(t)=x_u\mathrm{cos}(\omega t +\theta)=\mathrm{Re}[\overline x(t)]=\mathrm{Re}(X\mathrm{e}^{\mathrm{j}\omega t}) \tag{3-80}
x(t)=xucos(ωt+θ)=Re[x(t)]=Re(Xejωt)(3-80)
式中,
x
‾
(
t
)
\overline x(t)
x(t),
X
X
X分别为简谐函数
x
(
t
)
x(t)
x(t)在余弦函数表示时的(单边)复数表示和(单边)复振幅。
这样要取实部进行计算比较繁琐,因此我们考虑利用双边复数表示对余弦函数进行分解。
x
(
t
)
=
x
u
c
o
s
(
ω
t
+
θ
)
=
1
2
x
u
e
j
(
ω
t
+
θ
)
+
1
2
x
u
e
−
j
(
ω
t
+
θ
)
=
x
‾
B
+
(
t
)
+
x
‾
B
−
(
t
)
(3-81)
x(t)=x_u \mathrm{cos}(\omega t + \theta)={1\over 2}x_u \mathrm{e}^{\mathrm{j}(\omega t+\theta)}+{1\over 2}x_u\mathrm{e}^{-\mathrm{j}(\omega t +\theta)}=\overline x_{B+}(t)+\overline x_{B-}(t) \tag{3-81}
x(t)=xucos(ωt+θ)=21xuej(ωt+θ)+21xue−j(ωt+θ)=xB+(t)+xB−(t)(3-81)
式中,
x
‾
B
+
(
t
)
\overline x_{B+}(t)
xB+(t),
x
‾
B
−
(
t
)
\overline x_{B-}(t)
xB−(t)为简谐函数
x
(
t
)
x(t)
x(t)的双边复数表示,为复指数函数,由下式表示:
x
‾
B
+
(
t
)
=
1
2
x
u
e
j
(
ω
t
+
θ
)
=
1
2
x
u
e
j
θ
⋅
e
j
ω
t
=
X
B
+
e
j
ω
t
(3-82)
\overline x_{B+}(t)={1\over 2}x_u \mathrm{e}^{\mathrm{j}(\omega t+\theta)}={1\over 2}x_u \mathrm{e}^{\mathrm{j}\theta}\cdot \mathrm{e}^{\mathrm{j}\omega t}=X_{B+}\mathrm{e}^{\mathrm{j}\omega t} \tag{3-82}
xB+(t)=21xuej(ωt+θ)=21xuejθ⋅ejωt=XB+ejωt(3-82)
x
‾
B
−
(
t
)
=
1
2
x
u
e
−
j
(
ω
t
+
θ
)
=
1
2
x
u
e
−
j
θ
⋅
e
−
j
ω
t
=
X
B
−
e
−
j
ω
t
(3-83)
\overline x_{B-}(t)={1\over 2}x_u \mathrm{e}^{-\mathrm{j}(\omega t+\theta)}={1\over 2}x_u \mathrm{e}^{-\mathrm{j}\theta}\cdot \mathrm{e}^{-\mathrm{j}\omega t}=X_{B-}\mathrm{e}^{-\mathrm{j}\omega t} \tag{3-83}
xB−(t)=21xue−j(ωt+θ)=21xue−jθ⋅e−jωt=XB−e−jωt(3-83)
其中,
X
B
+
X_{B+}
XB+,
X
B
−
X_{B-}
XB−均为复常数,称为
x
(
t
)
x(t)
x(t)的双边复振幅,由下式表示:
X
B
+
=
1
2
x
u
e
j
θ
X
B
−
=
1
2
x
u
e
−
j
θ
(3-84)
\begin{aligned} X_{B+}&={1\over 2}x_u\mathrm{e}^{\mathrm{j}\theta} \\X_{B-}&={1\over 2}x_u\mathrm{e}^{-\mathrm{j}\theta} \end{aligned} \tag{3-84}
XB+XB−=21xuejθ=21xue−jθ(3-84)
由此可知,双边复振幅
X
B
+
X_{B+}
XB+和
X
B
−
X_{B-}
XB−互为共轭复数,即
X
B
+
=
(
X
B
−
)
∗
X_{B+}=(X_{B-})^*
XB+=(XB−)∗,并且有
X
B
+
=
1
2
X
X_{B+}={1\over 2}X
XB+=21X。
函数满足某些条件就可以展开成傅立叶级数,分解形式为:
x
(
t
)
=
1
2
a
0
+
∑
k
=
1
∞
(
a
k
c
o
s
ω
k
t
+
b
k
s
i
n
ω
k
t
)
=
1
2
x
u
0
+
∑
k
=
1
∞
x
u
k
c
o
s
(
ω
k
t
+
θ
k
)
(3-85)
x(t)={1\over 2}a_0+\sum ^\infty _{k=1}(a_k \mathrm{cos}\omega_k t+b_k \mathrm{sin}\omega_k t)={1\over 2}x_{u0}+\sum^\infty_{k=1}x_{uk}\mathrm{cos}(\omega_k t+\theta_k)\tag{3-85}
x(t)=21a0+k=1∑∞(akcosωkt+bksinωkt)=21xu0+k=1∑∞xukcos(ωkt+θk)(3-85)
其中:
ω
k
=
2
k
π
T
\omega_k={2k\pi \over T}
ωk=T2kπ
a
k
=
1
π
∫
π
π
x
(
t
)
c
o
s
k
t
d
t
(
n
=
0
,
1
,
2
,
⋅
⋅
⋅
)
a_k={1\over \pi} \int^\pi_\pi x(t)\mathrm{cos}kt \mathrm{d}t(n=0,1,2, \cdot\cdot\cdot)
ak=π1∫ππx(t)cosktdt(n=0,1,2,⋅⋅⋅)
b
k
=
1
π
∫
π
π
x
(
t
)
s
i
n
k
t
d
t
(
n
=
1
,
2
,
⋅
⋅
⋅
)
b_k={1\over \pi} \int^\pi_\pi x(t)\mathrm{sin}kt \mathrm{d}t(n=1,2, \cdot\cdot\cdot)
bk=π1∫ππx(t)sinktdt(n=1,2,⋅⋅⋅)
式中,
x
(
t
)
x(t)
x(t)被分解为一系列简谐振动分量的叠加,分解后的简谐振动分量称为周期函数
x
(
t
)
x(t)
x(t)的谐波;
a
k
a_k
ak,
b
k
b_k
bk为Fourier级数的系数;
x
u
k
x_uk
xuk,
θ
k
\theta_k
θk为周期函数
x
(
t
)
x(t)
x(t)各谐波的振幅和相位角;
ω
k
\omega_k
ωk为各谐波的振动圆频率。
如果将周期函数的每个谐波按式(3-80)转化为复数表示,则可得周期函数的复Fourier级数展开:
x
(
t
)
=
1
2
X
0
+
∑
k
=
1
∞
R
e
(
X
k
e
j
ω
k
t
)
=
∑
k
=
−
∞
∞
X
B
k
e
j
ω
k
t
(3-86)
x(t)={1\over 2}X_0+\sum^\infty_{k=1}\mathrm{Re}(X_k\mathrm{e}^{\mathrm{j}\omega_k t})=\sum^\infty_{k=-\infty}X_{Bk}\mathrm{e}^{\mathrm{j}\omega_k t}\tag{3-86}
x(t)=21X0+k=1∑∞Re(Xkejωkt)=k=−∞∑∞XBkejωkt(3-86)
式中,
X
k
X_k
Xk和
X
B
k
X_{Bk}
XBk分别为周期函数
x
(
t
)
x(t)
x(t)谐波的单边和双边的复振幅,也称为复Fourier级数展开的复系数,由下式计算:
X
B
k
=
1
T
∫
0
T
x
(
t
)
e
−
j
ω
k
t
d
t
,
k
=
⋅
⋅
⋅
,
−
2
,
−
1
,
0
,
1
,
2
,
⋅
⋅
⋅
(3-87)
X_{Bk}={1\over T}\int^T_0x(t)\mathrm{e}^{-\mathrm{j}\omega_k t}\mathrm{d}t,k=\cdot \cdot \cdot, -2, -1, 0, 1, 2, \cdot \cdot \cdot \tag{3-87}
XBk=T1∫0Tx(t)e−jωktdt,k=⋅⋅⋅,−2,−1,0,1,2,⋅⋅⋅(3-87)
X
B
k
=
2
X
B
k
=
2
T
∫
0
T
x
(
t
)
e
−
j
ω
k
t
d
t
,
k
=
0
,
1
,
2
,
⋅
⋅
⋅
(3-88)
X_{Bk}=2X_{Bk}={2\over T}\int^T_0x(t)\mathrm{e}^{-\mathrm{j}\omega_k t}\mathrm{d}t,k=0, 1, 2, \cdot \cdot \cdot \tag{3-88}
XBk=2XBk=T2∫0Tx(t)e−jωktdt,k=0,1,2,⋅⋅⋅(3-88)
通过周期函数
x
(
t
)
x(t)
x(t)的复Fourier级数展开,就建立了周期函数与复振幅之间的一一对应关系,即:
x
(
t
)
↔
X
k
,
k
=
0
,
1
,
2
,
3
⋅
⋅
⋅
或
x
(
t
)
↔
X
B
k
,
k
=
⋅
⋅
⋅
,
−
2
,
−
1
,
0
,
1
,
2
,
⋅
⋅
⋅
(3-89)
x(t)\leftrightarrow X_k, k=0,1,2,3 \cdot\cdot\cdot或x(t)\leftrightarrow X_{Bk},k=\cdot\cdot\cdot,-2,-1,0,1,2,\cdot\cdot\cdot \tag{3-89}
x(t)↔Xk,k=0,1,2,3⋅⋅⋅或x(t)↔XBk,k=⋅⋅⋅,−2,−1,0,1,2,⋅⋅⋅(3-89)
3.5.2 非周期激励的谐波分解
将Fourier级数推广为Fourier变换就可以用于非周期振动的谐波分析。通过Fourier变换,可以将一个非周期函数分解为在所有频率范围内连续分布的谐波分量,即
x
(
t
)
=
1
2
π
∫
0
∞
R
e
[
X
(
ω
)
e
j
ω
t
]
d
ω
=
1
2
π
∫
−
∞
∞
X
B
(
ω
)
e
j
ω
t
d
ω
(3-90)
x(t)={1\over 2\pi}\int^ \infty_0\mathrm{Re}[X(\omega)\mathrm{e}^{\mathrm{j}\omega t}]\mathrm{d}\omega={1\over 2 \pi}\int^{\infty}_{-\infty}X_B(\omega)\mathrm{e}^{\mathrm{j}\omega t}\mathrm{d}\omega \tag{3-90}
x(t)=2π1∫0∞Re[X(ω)ejωt]dω=2π1∫−∞∞XB(ω)ejωtdω(3-90)
X
(
ω
)
=
2
∫
−
∞
∞
x
(
t
)
e
−
j
ω
t
d
t
,
X
B
(
ω
)
=
∫
−
∞
∞
x
(
t
)
e
−
j
ω
t
d
t
(3-91)
X(\omega)=2\int^{\infty}_{-\infty}x(t)\mathrm{e}^{-\mathrm{j}\omega t}\mathrm{d}t, X_B(\omega)=\int^{\infty}_{-\infty}x(t)\mathrm{e}^{-\mathrm{j}\omega t}\mathrm{d}t \tag{3-91}
X(ω)=2∫−∞∞x(t)e−jωtdt,XB(ω)=∫−∞∞x(t)e−jωtdt(3-91)
其中,
X
(
ω
)
X(\omega)
X(ω),
X
B
(
ω
)
X_B(\omega)
XB(ω)均为在频率域内连续分布的复函数,分别称为单边(余弦表示形式时的)幅值谱密度和双边幅值谱密度,简称为幅值谱。通过非周期函数
x
(
t
)
x(t)
x(t)的复Fourier变换,建立了非周期振动与幅值谱之间的一一对应关系,即
x
(
t
)
↔
X
(
ω
)
,
ω
∈
[
0
,
+
∞
)
或
x
(
t
)
↔
X
B
(
ω
)
,
ω
∈
(
−
∞
,
+
∞
)
(3-92)
x(t)\leftrightarrow X_(\omega),\omega \in [0,+\infty )或x(t)\leftrightarrow X_B(\omega),\omega \in (-\infty, +\infty) \tag{3-92}
x(t)↔X(ω),ω∈[0,+∞)或x(t)↔XB(ω),ω∈(−∞,+∞)(3-92)
x
(
t
)
x(t)
x(t)和
X
(
ω
)
X(\omega)
X(ω)或
X
B
(
ω
)
X_B(\omega)
XB(ω)称为一个Fourier变换对,称
X
(
ω
)
X(\omega)
X(ω)或
X
B
(
ω
)
X_B(\omega)
XB(ω)为
x
(
t
)
x(t)
x(t)的复Fourier变换,
x
(
t
)
x(t)
x(t)为
X
(
ω
)
X(\omega)
X(ω)或
X
B
(
ω
)
X_B(\omega)
XB(ω)的复Fourier逆变换。
3.5.3 非简谐激励的稳态响应
设
x
t
(
t
)
x_t(t)
xt(t),
x
2
(
t
)
x_2(t)
x2(t)分别对应于激励函数
p
1
(
t
)
p_1(t)
p1(t),
p
2
(
t
)
p_2(t)
p2(t)的解,由此我们可以获得:
m
(
x
¨
1
+
x
¨
2
)
+
c
(
x
˙
1
+
x
˙
2
)
+
k
(
x
1
+
x
2
)
=
p
1
(
t
)
+
p
2
(
t
)
(3-93)
m(\ddot x_1+\ddot x_2)+c(\dot x_1 +\dot x_2)+k(x_1+x_2)=p_1(t)+p_2(t) \tag{3-93}
m(x¨1+x¨2)+c(x˙1+x˙2)+k(x1+x2)=p1(t)+p2(t)(3-93)
也就是
x
(
t
)
=
x
1
(
t
)
+
x
2
(
t
)
x(t)=x_1(t)+x_2(t)
x(t)=x1(t)+x2(t)对应于激励函数
p
(
t
)
=
p
1
(
t
)
+
p
2
(
t
)
p(t)=p_1(t)+p_2(t)
p(t)=p1(t)+p2(t)的解。也就是说这是一个线性系统,我们可以利用叠加原理。
由此周期为
T
T
T的周期激励
p
(
t
)
p(t)
p(t)由叠加原理可得非简谐激励强迫振动的解,系统强迫振动的稳态响应为:
x
(
t
)
=
∑
i
=
−
∞
+
∞
H
(
ω
i
)
⋅
P
B
i
e
j
ω
i
t
=
∑
i
=
−
∞
+
∞
P
B
i
k
⋅
1
(
1
−
λ
i
2
)
+
j
2
ζ
λ
i
⋅
e
j
ω
i
t
(3-95)
\begin{aligned} x(t)&=\sum^{+\infty}_{i=-\infty}H(\omega_i)\cdot P_{Bi}\mathrm{e}^{\mathrm{j}\omega_i t} \\&=\sum^{+\infty}_{i=-\infty}{P_{Bi}\over k}\cdot {1 \over (1-\lambda^2_i)+\mathrm{j}2\zeta\lambda_i}\cdot \mathrm{e}^{\mathrm{j}\omega_i t} \end{aligned} \tag{3-95}
x(t)=i=−∞∑+∞H(ωi)⋅PBiejωit=i=−∞∑+∞kPBi⋅(1−λi2)+j2ζλi1⋅ejωit(3-95)
当激励为非周期激励时,叠加计算相应转变为积分计算,与式(3-95)对应的有
x
(
t
)
=
1
2
π
∫
−
∞
∞
H
f
,
x
(
ω
)
⋅
P
B
(
ω
)
e
j
ω
t
d
ω
=
1
2
π
∫
−
∞
∞
P
B
(
ω
)
k
⋅
1
(
1
−
λ
2
)
+
j
2
ζ
λ
⋅
e
j
ω
t
d
ω
(3-96)
\begin{aligned} x(t)&={1\over 2 \pi} \int ^\infty _{-\infty}H_{f,x}(\omega)\cdot P_B(\omega)\mathrm{e}^{\mathrm{j}\omega t}\mathrm{d}\omega \\&={1\over 2 \pi} \int ^\infty _{-\infty}{P_B(\omega)\over k}\cdot {1 \over (1-\lambda^2)+\mathrm{j}2\zeta\lambda}\cdot \mathrm{e}^{\mathrm{j}\omega t}\mathrm{d}\omega \end{aligned} \tag{3-96}
x(t)=2π1∫−∞∞Hf,x(ω)⋅PB(ω)ejωtdω=2π1∫−∞∞kPB(ω)⋅(1−λ2)+j2ζλ1⋅ejωtdω(3-96)