本发明涉及非线性系统振动响应分析，尤其涉及一种次频振动边界条件的计算方法。
背景技术：
1、工程领域中，复杂的结构形状与组成部件极易引起系统的非线性振动。当外激励频率达到一定值时，非线性系统可能出现1/3阶次谐振动，这一现象已被学者大量研究，但相关研究中，1/3阶次谐振动对应的频率边界尚未被理论研究。
2、目前，在非线性振动领域中，谐波平衡法是寻求系统稳态周期振动的一种方法。其大体思路为：将响应方程代入平衡方程，得到幅频特性方程，根据幅频特性方程确定频率边界。该方法已被成功应用于确定一阶谐波的解，但尚未被成功用于确定1/3阶次谐振动，原因在于利用谐波平衡法理论研究1/3阶次谐振动时，方程数目小于未知量数目，导致无法准确求解得到1/3阶次谐振动，也就无法确定其振动边界。
技术实现思路
1、本发明提供了一种次频振动边界条件的计算方法，能够解决现有技术中无法确定次频振动频率边界的技术问题。
2、根据本发明的一方面，提供了一种次频振动边界条件的计算方法，计算方法包括：
3、根据谐波平衡方程和次频振动响应方程得到基于次频振动参数的参数方程；
4、根据能量理论建立基于次频振动参数的补充方程；
5、根据参数方程和补充方程计算得到次频振动参数的多组解；
6、将次频振动参数的多组解分别带入次频振动响应方程得到多组对应的振动响应曲线；
7、根据多组振动响应曲线对应得到多个频幅响应曲线，并根据多个频幅响应曲线确定次频振动的边界条件。
8、进一步地，次频振动为1/3阶次谐振动，次频振动参数包括一阶谐波分量的幅值、1/3阶次谐波分量的幅值、一阶谐波分量的相位和1/3阶次谐波分量的相位。
9、进一步地，谐波平衡方程为：
10、
11、上式中，x表示系统位移响应，ω0表示系统固有频率，ξ表示阻尼系数，ε表示非线性系数，b表示激振力幅值系数，ω表示激振力频率，t表示振动时间。
12、进一步地，次频振动响应方程为：
13、
14、上式中，a1表示一阶谐波分量的幅值，a2表示1/3阶次谐波分量的幅值，θ1表示一阶谐波分量的相位，θ2表示1/3阶次谐波分量的相位。
15、进一步地，根据能量理论建立的基于次频振动参数的补充方程为：
16、
17、进一步地，根据参数方程和补充方程计算得到的次频振动参数的解的数量为三组。
18、进一步地，将次频振动参数的多组解分别带入次频振动响应方程得到多组对应的振动响应曲线包括：
19、将次频振动参数的三组解分别带入次频振动响应方程以对应得到三个次频振动响应子方程；
20、对于每个次频振动响应子方程，均通过改变激振力频率的值得到多个振动响应曲线。
21、进一步地，根据多组振动响应曲线对应得到多个频幅响应曲线包括：
22、确定第一个次频振动响应子方程的多个振动响应曲线中每个曲线的响应峰值，并将其依次连接以得到第一幅频响应曲线；
23、确定第二个次频振动响应子方程的多个振动响应曲线中每个曲线的响应峰值，并将其依次连接以得到第二幅频响应曲线；
24、确定第三个次频振动响应子方程的多个振动响应曲线中每个曲线的响应峰值，并将其依次连接以得到第二幅频响应曲线。
25、进一步地，根据多个频幅响应曲线确定次频振动的边界条件包括：将第一幅频响应曲线、第二幅频响应曲线和第二幅频响应曲线这三条频幅响应曲线的交点对应的激振力频率确定为次频振动的边界条件。
26、应用本发明的技术方案，提供了一种次频振动边界条件的计算方法，该方法通过能量法得到补充方程，与谐波平衡法得到的参数方程组成完备方程组，从而能够求解得到次频振动的幅频响应曲线，进而根据幅频响应曲线快速确定次频振动的边界条件。该方法能够避免重复的数值积分运算，提高设计效率，为工程减振设计提供理论与技术支撑。
1.一种次频振动边界条件的计算方法，其特征在于，所述计算方法包括：
2.根据权利要求1所述的计算方法，其特征在于，所述次频振动为1/3阶次谐振动，所述次频振动参数包括一阶谐波分量的幅值、1/3阶次谐波分量的幅值、一阶谐波分量的相位和1/3阶次谐波分量的相位。
3.根据权利要求2所述的计算方法，其特征在于，所述谐波平衡方程为：
4.根据权利要求3所述的计算方法，其特征在于，所述次频振动响应方程为：
5.根据权利要求4所述的计算方法，其特征在于，根据能量理论建立的基于所述次频振动参数的补充方程为：
6.根据权利要求5所述的计算方法，其特征在于，根据所述参数方程和所述补充方程计算得到的次频振动参数的解的数量为三组。
7.根据权利要求6所述的计算方法，其特征在于，将所述次频振动参数的多组解分别带入所述次频振动响应方程得到多组对应的振动响应曲线包括：
8.根据权利要求7所述的计算方法，其特征在于，根据多组所述振动响应曲线对应得到多个频幅响应曲线包括：
9.根据权利要求8所述的计算方法，其特征在于，根据多个所述频幅响应曲线确定次频振动的边界条件包括：将所述第一幅频响应曲线、所述第二幅频响应曲线和所述第二幅频响应曲线这三条频幅响应曲线的交点对应的激振力频率确定为次频振动的边界条件。