CN103942374B - 一种刻画振动系统周期振幅与系统参数的关系的方法 - Google Patents
一种刻画振动系统周期振幅与系统参数的关系的方法 Download PDFInfo
- Publication number
- CN103942374B CN103942374B CN201410144087.XA CN201410144087A CN103942374B CN 103942374 B CN103942374 B CN 103942374B CN 201410144087 A CN201410144087 A CN 201410144087A CN 103942374 B CN103942374 B CN 103942374B
- Authority
- CN
- China
- Prior art keywords
- variable
- vibrational system
- solution
- relation
- parameter
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Fee Related
Links
Landscapes
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
Abstract
本发明公开了一种刻画振动系统周期振幅与系统参数的关系的方法。许多自然现象都具有周期性，大到天体运行，小到电子旋转。社会现象也其周期特性，如社会经济的涨跌，甚至民族的兴衰。研究这些周期现象可以预测至控制其运动演化规律。本发明给出的方法可以用于含有任意可数维数参数变量的任意自由度振动系统，可以计算参数变量的任意给定区域，可以求出具有多个周期解的振动系统的所有周期解，为振动系统的研究提供了有效的手段，为周期现象的研究提供了可靠的方法。
Description
技术领域
本发明公开了一种刻画振动系统周期振幅与系统参数的关系的方法，属于振动控制领域。
背景技术
狭义的振动又称机械振动，是指物体在平衡位置附近作的往复运动。广义的振动是指描述系统的状态变量在其基准值上下交替变化的过程。振动可分为线性振动和非线性振动。线性振动通常是对非线性振动简化和近似，自然界中真实的振动都是非线性的。周期振动是常见的振动形式。许多自然现象都具有周期性，大到天体运行，小到电子旋转。社会现象也其周期特性，如社会经济的涨跌，甚至民族的兴衰。研究这些周期现象可以预测至控制其运动演化规律。
振动系统通常用微分方程来描述，形式如下：
其中，t表示时间，x通常是向量，用来描述系统的状态变量，如物体的振动位移、电流的大小、声音的强弱，种群的数量，经济的总量等，表示状态变量对时间的导数，通常描述系统状态变量的变化速率，ε通常是向量，用来描述系统的参数变量，如惯性、弹性、激励的幅值或频率。g(t,x,ε)是关于时间、状态变量以及参数变量的非线性向量函数。
利用摄动法得到原系统的一阶近似解的调制方程具有如下形式：
X′＝G(X,ε), (2)
其中X描述原系统状态变量的一阶近似解的幅值和相位，X′表示X对一阶时间尺度的导数，G(X,ε)是非线性向量函数。
方程(2)的右端等于零时，调制方程有平衡解，即原系统状态变量的一阶近似解的幅值和相位都是常量，这对应原系统出现周期解。求解方程(2)的平衡解，等价于求解
G(X,ε)＝0. (3)
方程(3)称为调制方程的平衡解等价方程，是原系统的周期振动幅值和相位X关于系统参数变量ε的隐式函数。由于G(X,ε)通常是非线性函数，X很难表示成ε的显式函数。因此，用数值方法求解方程(3)成为得到原系统周期振动幅值和参数变量的关系的重要手段。有些软件含有求解非线性方程组或带参数非线性方程组的函数，如MATLAB的fsolve函数。该函数需要事先给出初始值，而且只能得到由该初始值求出的单个参数点处的一个解，但是振动系统周期振幅是参数变量的多值函数，而且刻画周期振幅与参数变量的关系需要在参数变量的整个取值范围进行。
根据一些科技文献，国外有些机构研发了一些用于非线性分析的软件，可以用来计算微分振动系统的周期解，如AUTO-07P，AnT，XPP-AUT等。这些软件通常运行在Linux环境下，或者需要在Windows系统中虚拟出Linux环境，而且其安装运行都很繁杂。国内有些文献中出现过类似的计算，但是其中的计算方法都是针对简单的系统，例如单自由度系统，或者单参数系统。
本发明针对以上情况给出一种刻画振动系统周期振幅与系统参数的关系的方法。本发明所述的方法利用嵌套循环结构遍历所有参数变量的取值范围，在每个参数点处计算周期振幅，从而刻画振动系统周期振幅与系统参数的关系。本发明可以在参数变量的取值范围内求出振动系统的多个周期解，并假设每个解在各自参数点处的邻域是连续的，由此以数值解的收敛性作为退出条件向该邻域的上下边界扩展，从而提高计算效率。本发明可以应用于任意可数维参数变量和任意自由度的振动系统，为振动系统的研究提供了有效的手段，为周期现象的研究提供了可靠的方法。
发明内容
本发明的目的在于，为克服现有方法的不足，给出一种刻画振动系统周期振幅与系统参数的关系的方法。针对具体振动系统的调制方程X′＝G(X,ε)，X＝[XM]表示M维状态变量，ε＝[εN]表示N维参数变量，本发明可以应用于任意可数维参数变量和任意可数维状态变量的调制方程。
1.本发明首先录入具体振动系统的状态变量的一阶近似解的调制方程的平衡解等价方程，即方程(3)。为了控制程序的运行，本发明定义了一些计算所需要的全局变量：记录解的最大个数的变量，限定最大求解数的变量，存储计算结果的变量，记录平衡解数据大小的变量，判断相同解的容差变量，限定振动系统维数的变量，划分参数变量网格数的变量，限定参数变量范围的变量。本发明给出嵌套循环结构遍历参数变量的整个取值范围，此嵌套循环结构的每一层循环结构对应一维参数变量εi(i＝1,......,n)。在嵌套循环结构的最底层，记ε＝ε0，即每一维参数变量都取一个具体数值。在嵌套循环结构的最底层定义局部变量：记录所给参数点处平衡解的变量，记录所给参数点处平衡解个数的变量，记录试探次数的变量，限定最大试探次数的变量。在嵌套循环结构的最底层，再以试探次数不超过最大限定试探次数和求得的解的个数不超过限定最大求解数作为退出条件给出一个循环结构。在此循环结构的内部，在求解范围内随机给出初始值X＝X0，并用迭代方法X＝X(X0)求出所有的解X＝ξ，使得调制方程的平衡解等价方程满足G(ξ,ε0)＝0。如果所求得的平衡解不与已经求解并记录的平衡解重复，利用局部变量记录此求得的平衡解。假设ξ在ε0处连续，即ξ+Δξ＝ξ(ε0+Δε)，以收敛性为退出条件向该平衡解所在的参数点的邻域的上下边界扩展即令ε＝ε0+Δε，以X＝ξ为初值迭代求解方程(3)，若能得到收敛的数值解ξ+Δξ，满足方程G(ξ+Δξ,ε0+Δε)＝0，则继续向边界扩展，否则就停止扩展。利用全局变量记录所有扩展后的收敛的平衡解，画出参数变量与平衡解的关系图，刻画出振动系统周期振幅与系统参数的关系X＝X(ε)。
X是向量，其元素个数可以是任意可数个。因此，本发明所给出的方法可用于任意可数维状态变量的调制方程，即用于任意自由度的振动系统。
ε是向量，其元素个数可以是任意可数个。因此，本发明所给出的方法可以用于含有任意可数维数参数变量的振动系统。
所定义的限定最大求解数的变量可以设定为任意整数，因此，本发明所给出的方法可以求出具有多个周期解的振动系统的所有周期解。
本发明所给出的方法可以计算参数变量的任意取值区域。
本发明的有益效果是：1)本发明给出的方法可以用于任意自由度的振动系统；2)本发明给出的方法可以用于含有任意可数维数的参数变量的振动系统；3)本发明给出的方法可以计算参数变量的任意给定区域；4)本发明所给出的方法可以求出具有多个周期解的振动系统的所有周期解。
附图说明
图1是本发明所述方法的计算程序设计流程图。
图2是本发明所述方法实施所得一维参数变量四维调制方程的平衡解，描述二自由度振动系统的周期振幅与一维参数变量外激励频率的关系以及周期振幅与一维参数变量外激励幅值的关系。
图3是本发明所述方法实施所得二维参数变量四维调制变量的平衡解，描述二自由度振动系统的周期振幅与二维参数变量外激励幅值和频率的关系。
具体实施方式
本发明可以用任意编程语言实现。MATLAB是常用的数学计算软件，本发明利用MATLAB软件通过具体实施例并结合附图作进一步详细的描述。
实施例
本发明利用MATLAB编写程序计算文献《Three-to-One Internal Resonances inParametrically Excited Hinged-Clamped Beams》中调制方程(31-34)：
令[q1′,p1′,q2′,p2′]＝0，方程(4)有平衡解，其对应的平衡解等价方程具有方程(3)的形式，
G(X,ε)＝0, (5)
其中
X＝[q1,p1,q2,p2], (6a)
ε＝[F,σ2], (6b)
方程(5)的解对应原系统的周期解，即文献《Three-to-One Internal ResonancesinParametrically ExcitedHinged-ClampedBeams》中方程(30)中Ak为常数。方程(30)如下
A1和A2分别是文献《Three-to-One Internal Resonances in ParametricallyExcitedHinged-Clamped Beams》中方程(23)的第一阶和第二阶周期振幅。方程(23)如下
其中，w1是文献中梁结构的横向振动位移。
方程(5)是含有两个可变参数的四维调制方程，两个可变参数分别是外激励的幅值F，外激励频率与结构第二阶自然频率的差值σ2，四维状态变量是[q1,p1,q2,p2]。方程(5)中其他参数均为已知。
1.图1给出N维参数变量任意可数维调制方程的平衡解的计算程序设计流程。为与文献中的计算结果进行对比，本实施例限定N＝2，以F和σ2作为参数变量，即本实施例计算二维参数变量四维调制方程的平衡解。程序首先录入方程(5)，然后定义全局变量：记录解的最大个数的变量，限定最大求解数的变量，存储计算结果的变量，记录平衡解数据大小的变量，判断相同解的容差变量，限定振动系统维数的变量，划分参数变量网格数的变量，限定参数变量范围的变量。
2.如图1所示，给出嵌套循环结构，即第1层至第N层循环，每一层循环对应一维参数变量。内层循环遍历完本层对应参数范围后返回其最近外层继续遍历外层所对应参数范围，直到遍历完参数变量所给定范围，退出嵌套循环结构，结束计算。本实施例给出两层嵌套循环结构，第一层遍历参数F的取值范围，第二层遍历参数σ2的取值范围。
3.如图1所示，在嵌套循环结构的最底层，即第N＝2层，定义局部变量：记录所给参数点处平衡解的变量，记录所给参数点处平衡解个数的变量，记录试探次数的变量，限定最大试探次数的变量。
4.如图1所示，嵌套循环结构的最底层再以试探次数不超过最大限定试探次数和求得的解的个数不超过限定最大求解数作为退出条件给出一个循环结构。在此循环结构的内部，在求解范围内随机给出初始值，并利用fsolve函数求解方程。如果所求得的平衡解符合收敛条件且不重复，即满足收敛条件且在全参数范围和该参数点处都不与已经求得的解相同，利用局部变量记录此求得的平衡解，直到试探次数大于限定最大试探次数的变量或求得的解的个数大于解的最大个数。
5.如图1所示，以收敛性为退出条件向记录所给参数点处平衡解的变量记录的平衡解所在的参数点的邻域的上下边界扩展，利用存储计算结果的变量记录所有扩展后的收敛的平衡解。
6.文献中的结论与本实施例的计算结果对比如图2所示。左边的曲线图分别是文献《Three-to-One Internal Resonances in Parametrically Excited Hinged-ClampedBeams》中的Figure 2和Figure 3，α1与σ2的曲线图表示周期振动一阶分量的幅值与外激励频率的关系，α2与σ2的曲线图表示周期振动二阶分量的幅值与外激励频率的关系，α1与F的曲线图表示周期振动一阶分量与外激励幅值的关系，α2与F的曲线图表示周期振动二阶分量与外激励幅值的关系。右边曲线图分别与其左边曲线图对应，是本实施例的计算结果。由图2可以看出，本发明所述方法能计算出一维参数变量四维调制方程的平衡解，刻画出了二自由度振动系统周期振幅与一维参数变量的关系。
7.本实施例的一个计算结果如图3所示。图3是文献《Three-to-One InternalResonances in Parametrically Excited Hinged-Clamped Beams》中的Figure 2和Figure 3综合。图3给出了周期振动一阶分量的幅值α1以及周期振动二阶分量的幅值α2分别在二维参数变量σ2和F中的取值情况。文献《Three-to-One Internal Resonances inParametrically Excited Hinged-Clamped Beams》没有给出与图3对应的情况。基于查阅大量文献，没有发现文献给出过二维以及高于二维参数变量四维调制方程的平衡解。由图3可以看出，本发明所述方法能计算出二维参数变量四维调制方程的平衡解。X和ε的元素个数都可以是任意可数个，因此，本发明所给出的方法能够刻画出含有任意可数维参数变量的任意自由度振动系统的周期振幅与参数变量的关系。
Claims (2)
1.一种刻画振动系统周期振幅与系统参数的关系的方法，其特征在于：首先录入具体振动系统的状态变量的一阶近似解的调制方程的平衡解等价方程；定义了一些计算所需要的全局变量：记录解的最大个数的变量，限定最大求解数的变量，存储计算结果的变量，记录平衡解数据大小的变量，判断相同解的容差变量，限定振动系统维数的变量，划分参数变量网格数的变量，限定参数变量范围的变量；给出嵌套循环结构遍历参数变量的整个取值范围，此嵌套循环结构的每一层循环结构对应一维参数变量εi，其中，i＝1,......,n；在嵌套循环结构的最底层，记ε＝ε0，即每一维参数变量都取一个具体数值；在嵌套循环结构的最底层定义局部变量：记录所给参数点处平衡解的变量，记录所给参数点处平衡解个数的变量，记录试探次数的变量，限定最大试探次数的变量；在嵌套循环结构的最底层，再以试探次数不超过最大限定试探次数和求得的解的个数不超过限定最大求解数作为退出条件给出一个循环结构；在此循环结构的内部，在求解范围内随机给出初值X＝X0，并用迭代方法X＝X(X0)求出ξ使得调制方程的平衡解等价方程满足G(ξ,ε0)＝0；如果所求得的平衡解不与已经求解并记录的平衡解重复，利用局部变量记录此求得的平衡解；假设ξ在ε0处连续，即ξ+Δξ＝ξ(ε0+Δε)，以收敛性为退出条件向该平衡解所在的参数点的邻域的上下边界扩展，即令ε＝ε0+Δε，以X＝ξ为初值迭代求解振动系统的一阶近似解的调制方程的平衡解等价方程,若能得到收敛的数值解ξ+Δξ，则继续向边界扩展，否则就停止扩展；利用全局变量记录所有扩展后的收敛的平衡解，画出参数变量与平衡解的关系图，刻画出振动系统周期振幅与系统参数的关系X＝X(ε)；
所定义的限定调制方程维数的变量，特征在于一种刻画振动系统周期振幅与系统参数的关系的方法适用于任意自由度的振动系统；
所述的嵌套循环结构，特征在于一种刻画振动系统周期振幅与系统参数的关系的方法可以用于含有任意可数维数参数变量的振动系统；
所定义的限定最大求解数的变量，特征在于一种刻画振动系统周期振幅与系统参数的关系的方法可以求出具有多个周期解的振动系统的所有周期解；
限定参数变量取值范围的变量，特征在于一种刻画振动系统周期振幅与系统参数的关系的方法可以计算参数变量的任意取值区域。
2.根据权利要求1所述一种刻画振动系统周期振幅与系统参数的关系的方法，其特征在于：所述状态变量的一阶近似解的调制方程描述的是周期振幅与系统参数的关系；其中，状态变量是振动系统的位移或速度可以任意随时间往复变化的特征量；系统参数是指振动系统的惯性、弹性以及系统受到的外激励的幅值及频率可以人为调节的固有的系统性质。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201410144087.XA CN103942374B (zh) | 2014-04-10 | 2014-04-10 | 一种刻画振动系统周期振幅与系统参数的关系的方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201410144087.XA CN103942374B (zh) | 2014-04-10 | 2014-04-10 | 一种刻画振动系统周期振幅与系统参数的关系的方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN103942374A CN103942374A (zh) | 2014-07-23 |
| CN103942374B true CN103942374B (zh) | 2017-04-05 |
Family
ID=51190042
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201410144087.XA Expired - Fee Related CN103942374B (zh) | 2014-04-10 | 2014-04-10 | 一种刻画振动系统周期振幅与系统参数的关系的方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN103942374B (zh) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110989497B (zh) * | 2019-12-10 | 2020-11-06 | 上海航天精密机械研究所 | 基于迭代优化的多轴多激励振动控制组合选择方法及系统 |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103034254A (zh) * | 2012-12-18 | 2013-04-10 | 重庆邮电大学 | 基于磁流变脂隔振器的地铁浮置板三向振动控制方法 |
| CN103559398A (zh) * | 2013-11-04 | 2014-02-05 | 湖南大学 | 汽车盘式制动器系统振动稳定性的不确定性分析方法 |
-
2014
- 2014-04-10 CN CN201410144087.XA patent/CN103942374B/zh not_active Expired - Fee Related
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103034254A (zh) * | 2012-12-18 | 2013-04-10 | 重庆邮电大学 | 基于磁流变脂隔振器的地铁浮置板三向振动控制方法 |
| CN103559398A (zh) * | 2013-11-04 | 2014-02-05 | 湖南大学 | 汽车盘式制动器系统振动稳定性的不确定性分析方法 |
Non-Patent Citations (1)
| Title |
|---|
| 复合材料角铺设层合板非线性振动特性研究;郭翔鹰等;《振动与冲击》;20121231;第31卷(第19期);第174-179页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN103942374A (zh) | 2014-07-23 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Wang et al. | Incorporating symmetry into deep dynamics models for improved generalization | |
| Dumbser et al. | ADER-WENO finite volume schemes with space–time adaptive mesh refinement | |
| Oseledets | Constructive representation of functions in low-rank tensor formats | |
| CN102648425B (zh) | 用于储集层仿真的多尺度有限体积方法 | |
| Kühnlein et al. | Modelling atmospheric flows with adaptive moving meshes | |
| Wang | Optimal point-wise error estimate of a compact difference scheme for the coupled Gross–Pitaevskii equations in one dimension | |
| Kočvara et al. | Free material optimization: recent progress | |
| Yang et al. | Simple positivity-preserving nonlinear finite volume scheme for subdiffusion equations on general non-conforming distorted meshes | |
| Antonietti et al. | A uniform additive Schwarz preconditioner for high-order discontinuous Galerkin approximations of elliptic problems | |
| Liu et al. | Unconditionally energy stable DG schemes for the Swift–Hohenberg equation | |
| Gao et al. | Aggregation-based method for computing absolute Boltzmann entropy of landscape gradient with full thermodynamic consistency | |
| Pautz et al. | Parallel deterministic transport sweeps of structured and unstructured meshes with overloaded mesh decompositions | |
| Degond et al. | Alignment of self-propelled rigid bodies: from particle systems to macroscopic equations | |
| Carrillo et al. | Compact approximate Taylor methods for systems of conservation laws | |
| Li et al. | Bifurcation and number of subharmonic solutions of a 4D non-autonomous slow–fast system and its application | |
| Wang et al. | The Magnus expansion for stochastic differential equations | |
| Lourenço et al. | An extension of Chubanov’s algorithm to symmetric cones | |
| Helzel et al. | A high-order unstaggered constrained-transport method for the three-dimensional ideal magnetohydrodynamic equations based on the method of lines | |
| Chen et al. | A dimensional splitting exponential time differencing scheme for multidimensional fractional Allen-Cahn equations | |
| Liu et al. | High resolution, entropy-consistent scheme using flux limiter for hyperbolic systems of conservation laws | |
| Axenides et al. | Strange attractors in dissipative Nambu mechanics: classical and quantum aspects | |
| Kanno | A fast first-order optimization approach to elastoplastic analysis of skeletal structures | |
| CN103942374B (zh) | 一种刻画振动系统周期振幅与系统参数的关系的方法 | |
| Behrens et al. | Efficiency considerations in triangular adaptive mesh refinement | |
| Dumbser et al. | On numerical methods for hyperbolic PDE with curl involutions |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20170405 Termination date: 20200410 |
|
| CF01 | Termination of patent right due to non-payment of annual fee |