1.一种基于声子晶体的汽车减震器的活塞杆减震方法，其特征在于：包括以下步骤：(1)确定汽车减震器的活塞杆的尺寸以及活塞杆工作时受到的振动频率区间；

(2)将活塞杆近似成有限自由度的活塞杆离散结构，建立所述活塞杆离散结构‑弹簧质子的振动方程；所述振动方程由下列公式推导得出：第j个质子质量mj表示为：

2

mj＝πρjdjr       (1)式中，π表示圆周率；ρj表示第j个质子的密度；dj表示每j个质子对应的活塞杆离散单元长度；r表示活塞杆的半径；

质子所受的力即活塞杆离散单元左侧的应力F‑，与质子所对应的活塞杆离散单元左侧应变△x‑，以及活塞杆离散单元右侧的应力F+，与质子所对应的活塞杆离散单元右侧应变△x+的关系如下：式中，E为质子对应的离散单元的杨氏模量；

由式(2)和(3)得第j个质子左侧的等效弹簧刚度kj,‑和右侧的等效弹簧刚度kj,+如下：由式(4)和(5)得第j个质子和第j+1个质子间的第j个弹簧的等效弹簧刚度kj为：式中，kj+1,‑为第j+1个质子左侧的等效弹簧刚度，kj,+为第j个质子右侧的等效弹簧刚度；Ej+1为活塞杆第j+1段离散单元的杨氏模量；Ej为活塞杆第j段离散单元的杨氏模量；

式(6)结合牛顿第二定律得活塞杆离散结构‑‑弹簧质子链的振动方程如下：式中，第j个质子的位移 Aj表示第j个质子振幅，i表示虚数的单位，q表示波矢，l表示杆的长度，n表示离散单元数，ω表示频率，t表示时间；xj‑1第j‑1个质子的位移；xj+1第j+1个质子的位移；kj‑1表示第j‑1个质子和第j个质子间的等效弹簧刚度；

(3)在振动方程的基础上结合Bloch定理建立活塞杆减震频率计算的弥散方程；推导过程如下：将式(1)和式(6)带入式(7)得：

在周期边界条件下不同元胞之间的参量的关系如下：式中，k表示等效弹簧刚度；m表示振子质量；d表示离散单元的轴向长度；A表示质子振幅；下标符号对应的是离散结构的序列数，如k1表示单个周期结构的离散结构‑‑质子链中第个1弹簧的等效刚度；n为活塞杆元胞离散结构的振子数；

将(9)式带入(8)式并展开为矩阵形式，表示为：2

[ωI‑X(q)]A＝0      (10)式中，I表示单位矩阵；

T

A＝[A1 A2 A3 A4 A5 … An‑3 An‑2 An‑1 An]      (12)式中，n为活塞杆元胞离散结构的振子数；

2

由矩阵X(q)‑ωI行列式为零，得活塞杆元胞的弥散方程如下：

2n 2(n‑1) 8 6 4 2

Bnω +Bn‑1ω +…+B4ω+B3ω+B2ω+B1ω+B0＝0     (13)其中，B0～Bn为与由q、r、E、ρ、l、n组成的一系列系数；

当方程(8)不等于0时，定义弥散方程如下：

2 2n 2(n‑1) 8 6 4 2f(ω)＝Bnω +Bn‑1ω +…+B4ω+B3ω+B2ω+B1ω+B0 (14)；

其中，ω为活塞杆的固有频率；

(4)基于弥散方程建立活塞杆材料参数满足减震要求的存在性依据；设定目标减震频率区间为 基于弥散方程(13)建立所述活塞杆材料参数满足减振要求的存在性依据如下：*

式中，波矢q＝q＝0或π/l；ωa表示初始频率，ωb表示截止频率；

(5)将活塞杆的已知参数和未知参数带入所述步骤(4)中的存在性依据中判定所述活塞杆的振动频率是否在振动频率区间内；所述已知参数包括结构参数，所述结构参数为半径、长度；所述未知参数包括弹性模量和密度。

2.根据权利要求1所述的基于声子晶体的汽车减震器的活塞杆减震方法，其特征在于：步骤(5)中，当活塞杆在振动频率区间外时，通过有限元分析判定所述活塞杆的振动频率在步骤(1)的振动频率区间内是否存在衰减。

3.根据权利要求1所述的基于声子晶体的汽车减震器的活塞杆减震方法，其特征在于：所述活塞杆为由不同声子晶体材料制成的活塞杆离散结构。

4.根据权利要求3所述的基于声子晶体的汽车减震器的活塞杆减震方法，其特征在于：所述活塞杆离散结构包括多个与质子对应的离散单元。

5.根据权利要求1所述的基于声子晶体的汽车减震器的活塞杆减震方法，其特征在于：步骤(2)中，所述活塞杆离散结构为弹簧质子结构。