一种刻画振动系统周期振幅与系统参数的关系的方法
【专利摘要】本发明公开了一种刻画振动系统周期振幅与系统参数的关系的方法。许多自然现象都具有周期性，大到天体运行，小到电子旋转。社会现象也其周期特性，如社会经济的涨跌，甚至民族的兴衰。研究这些周期现象可以预测至控制其运动演化规律。本发明给出的方法可以用于含有任意可数维数参数变量的任意自由度振动系统，可以计算参数变量的任意给定区域，可以求出具有多个周期解的振动系统的所有周期解，为振动系统的研究提供了有效的手段，为周期现象的研究提供了可靠的方法。
【专利说明】一种刻画振动系统周期振幅与系统参数的关系的方法
【技术领域】
[0001]本发明公开了一种刻画振动系统周期振幅与系统参数的关系的方法，属于振动控制领域。
【背景技术】
[0002]狭义的振动又称机械振动，是指物体在平衡位置附近作的往复运动。广义的振动是指描述系统的状态变量在其基准值上下交替变化的过程。振动可分为线性振动和非线性振动。线性振动通常是对非线性振动简化和近似，自然界中真实的振动都是非线性的。周期振动是常见的振动形式。许多自然现象都具有周期性，大到天体运行，小到电子旋转。社会现象也其周期特性，如社会经济的涨跌，甚至民族的兴衰。研究这些周期现象可以预测至控制其运动演化规律。
[0003]振动系统通常用微分方程来描述，形式如下:
[0004]X = g(t,x,s).(I)
[0005]其中，t是标量，表示时间。X通常是向量，用来描述系统的状态变量，如物体的位移、电流的大小、声音的强弱，种群的数量，经济的总量等。i表示状态变量对时间的导数，通常描述系统状态变量的变化速率。ε通常是向量，用来描述系统的参数变量，如惯性、弹性、激励的幅值或频率。g(t，x，ε )是关于时间、状态变量以及参数变量的非线性向量函数。
[0006]利用摄动法得到原系统的一阶近似解的调制方程具有如下形式:
[0007]X' =G(X, ε).(2)
[0008]其中X描述原系统状态变量的一阶近似解的幅值和相位，X'表示X对一阶时间尺度的导数，G(X，ε)是非线性向量函数。
[0009]方程(2)的右端等于零时，调制方程有平衡解，即原系统状态变量的一阶近似解的幅值和相位都是常量，这对应原系统出现周期解。求解方程(2)的平衡解，等价于求解
[0010]G(X, ε ) = 0.(3)
[0011]方程(3)称为调制方程的平衡解等价方程，是原系统的周期振动幅值和相位X关于系统参数变量ε的隐式函数。由于G(X，ε)通常是非线性函数，X很难表示成ε的显式函数。因此，数值求解变系数非线性方程(3)成为得到原系统周期振动幅值和参数变量的关系的重要手段。有些软件含有求解非线性方程组或带参数非线性方程组的函数，如MATLAB的fsolve函数。该函数需要事先给出初始值，而且只能得到由该初始值求出的单个参数点处的一个解，但是振动系统周期振幅是参数变量的多值函数，而且刻画周期振幅与参数变量的关系需要在参数变量的整个取值范围进行。
[0012]根据一些科技文献，国外有些机构研发了一些用于非线性分析的软件，可以用来计算微分振动系统的周期解，如AUT0-07P，AnT,XPP-AUT等。这些软件通常运行在Linux环境下，或者需要在Windows系统中虚拟出Linux环境，而且其安装运行都很繁杂。国内有些文献中出现过类似的计算，但是其中的计算方法都是针对简单的系统，例如单自由度系统，或者单参数系统。[0013]本发明针对以上情况给出一种刻画振动系统周期振幅与系统参数的关系的方法。本发明所述的方法利用嵌套循环结构遍历所有参数变量的取值范围，在每个参数点处计算周期振幅，从而刻画振动系统周期振幅与系统参数的关系。本发明可以在参数变量的取值范围内求出振动系统的多个周期解，并假设每个解在各自参数点处的邻域是连续的，由此以数值解的收敛性作为退出条件向该邻域的上下边界扩展，从而提高计算效率。本发明可以应用于任意可数维参数变量和任意自由度的振动系统，为振动系统的研究提供了有效的手段，为周期现象的研究提供了可靠的方法。
【发明内容】
[0014]本发明的目的在于，为克服现有方法的不足，给出一种刻画振动系统周期振幅与系统参数的关系的方法。针对具体振动系统的调制方程X' =G(X, ε),Χ= [Xm]表示m维状态变量，ε = [εη]表示η维参数变量，本发明可以应用于任意可数维参数变量和任意可数维状态变量的调制方程。
[0015]1.本发明首先录入具体振动系统的状态变量的一阶近似解的调制方程的平衡解等价方程，即方程(3)。为了控制程序的运行，本发明定义了一些计算所需要的全局变量:记录解的最大个数的变量，限定最大求解数的变量，存储计算结果的变量，记录平衡解数据大小的变量，判断相同解的容差变量，限定振动系统维数的变量，划分参数变量网格数的变量，限定参数变量范围的变量。本发明给出嵌套循环结构遍历参数变量的整个取值范围，此嵌套循环结构的每一层循环结构对应一维参数变量ε i (i = 1...η)。在嵌套循环结构的最底层，记ε = ε°，即每一维参数变量都取一个固定数值。在嵌套循环结构的最底层定义局部变量:记录所给参数点处平衡解的变量，记录所给参数点处平衡解个数的变量，记录试探次数的变量，限定最大试探次数的变量。在嵌套循环结构的最底层，再以试探次数不超过最大限定试探次数和求得的解的个数不超过限定最大求解数作为退出条件给出一个循环结构。在此循环结构的内部，在求解范围内随机给出初值X = Χ°，并用迭代方法X = Χ(Χ°)求出ξ使得调制方程的平衡解等价方程满足GU，ε°) =O0如果所求得的平衡解不与已经求解并记录的平衡解重复，利用局部变量记录此求得的平衡解。假设ξ在ε°处连续，即ξ+Λ ξ = ξ (ε°+Λ ε)，以收敛性为退出条件向该平衡解所在的参数点的邻域的上下边界扩展即令ε = ε°+Λ ε，以X = ξ为初值迭代求解方程(3)，若能得到收敛的数值解ξ+Λ ξ，则继续向边界扩展，否则就停止扩展。利用全局变量记录所有扩展后的收敛的平衡解，画出参数变量与平衡解的关系图，刻画出振动系统周期振幅与系统参数的关系X =Χ( ε )。
[0016]其中，状态变量可以是任意随时间往复变化的特征量，例如机械的振动、电磁的振荡、种群数目的变化、经济总量的起伏。系统参数是指可以人为调节的固有的系统性质，如惯性、弹性、外激励幅值及频率。
[0017]本发明所给出的方法可用于任意可数维状态变量的调制方程，即用于任意自由度的振动系统。
[0018]本发明所给出的方法可以用于含有任意可数维数参数变量的振动系统。
[0019 ]所定义的限定最大求解数的变量，本发明所给出的方法可以求出具有多个周期解的振动系统的所有周期解。[0020]本发明所给出的方法可以计算参数变量的任意取值区域。
[0021]所述的假设平衡解在其参数点的邻域内连续，以收敛性为退出条件向该平衡解所在的参数点的邻域的上下边界扩展，该方法提高了程序的计算效率。
[0022]本发明的有益效果是:1)本发明给出的方法可以用于任意自由度的振动系统；2)本发明给出的方法可以用于含有任意可数维数的参数变量的振动系统；3)本发明给出的方法可以计算参数变量的任意给定区域；4)本发明所给出的方法可以求出具有多个周期解的振动系统的所有周期解。
【专利附图】
【附图说明】
[0023]图1是本发明所述方法的计算程序设计流程图。
[0024]图2是本发明所述方法实施所得一维参数变量四维调制方程的平衡解，描述二自由度振动系统的周期振幅与一维参数变量外激励频率的关系以及周期振幅与一维参数变量外激励幅值的关系。
[0025]图3是本发明所述方法实施所得二维参数变量四维调制变量的平衡解，描述二自由度振动系统的周期振幅与二维参数变量外激励幅值和频率的关系。
【具体实施方式】
[0026]本发明可以用任意编程语言实现。MATLAB是常用的数学计算软件，本发明利用MATLAB软件通过具体实施例并结合附图作进一步详细的描述。
[0027]实施例
[0028]本发明利用MATLAB 编写程序计算文献《Three-to-One Internal Resonances inParametrically Excited Hinged-Clamped Beams))中调制方程(31-34):
【权利要求】
1.一种刻画振动系统周期振幅与系统参数的关系的方法，其特征在于:本发明首先录入具体振动系统的状态变量的一阶近似解的调制方程的平衡解等价方程；为了控制程序的运行，本发明定义了一些计算所需要的全局变量:记录解的最大个数的变量，限定最大求解数的变量，存储计算结果的变量，记录平衡解数据大小的变量，判断相同解的容差变量，限定振动系统维数的变量，划分参数变量网格数的变量，限定参数变量范围的变量；本发明给出嵌套循环结构遍历参数变量的整个取值范围，此嵌套循环结构的每一层循环结构对应一维参数变量SiG = 1...n);在嵌套循环结构的最底层，记ε = ε°，即每一维参数变量都取一个固定数值；在嵌套循环结构的最底层定义局部变量:记录所给参数点处平衡解的变量，记录所给参数点处平衡解个数的变量，记录试探次数的变量，限定最大试探次数的变量；在嵌套循环结构的最底层，再以试探次数不超过最大限定试探次数和求得的解的个数不超过限定最大求解数作为退出条件给出一个循环结构；在此循环结构的内部，在求解范围内随机给出初值X = Χ°，并用迭代方法X = Χ(Χ°)求出I使得调制方程的平衡解等价方程满足GU，ε°) = O;如果所求得的平衡解不与已经求解并记录的平衡解重复，利用局部变量记录此求得的平衡解；假设ξ在ε °处连续，即ξ + Λ ξ = ξ ( ε °+Λ ε )，以收敛性为退出条件向该平衡解所在的参数点的邻域的上下边界扩展，即令ε = ε °+Δ ε，以X= ξ为初值迭代求解方程(3)，若能得到收敛的数值解ξ+Λ ξ，则继续向边界扩展，否则就停止扩展；利用全局变量记录所有扩展后的收敛的平衡解，画出参数变量与平衡解的关系图，刻画出振动系统周期振幅与系统参数的关系X = Χ( ε )； 所定义的限定调制方程维数的变量，本发明所给出的方法可用于任意可数维状态变量的调制方程，即用于任意自由度的振动系统； 所述的嵌套循环 结构，本发明所给出的方法可以用于含有任意可数维数参数变量的振动系统； 所定义的限定最大求解数的变量，本发明所给出的方法可以求出具有多个周期解的振动系统的所有周期解； 限定参数变量取值范围的变量，本发明所给出的方法可以计算参数变量的任意取值区域； 所述的假设平衡解在其参数点的邻域内连续，以收敛性为退出条件向该平衡解所在的参数点的邻域的上下边界扩展，该方法提高了程序的计算效率。
2.根据权利I所述一种刻画振动系统周期振幅与系统参数的关系方法，其特征在于:所述状态变量的一阶近似解的调制方程描述的是周期振幅与系统参数的关系；其中，状态变量可以是任意随时间往复变化的特征量，例如机械的振动、电磁的振荡、种群数目的变化、经济总量的起伏；系统参数是指可以人为调节的固有的系统性质，如惯性、弹性、外激励幅值及频率。
【文档编号】G06F17/50GK103942374SQ201410144087
【公开日】2014年7月23日 申请日期:2014年4月10日 优先权日:2014年4月10日
【发明者】张伟, 胡文华, 曹东兴, 杨晓峰 申请人:北京工业大学