文摘

能量平衡法(实证)延长高阶非线性振子。为了说明方法的有效性,cubic-quintic杜芬振荡器被选中。的最大相对误差的频率振荡器读第一和二阶近似1.25%和0.6%,分别。三阶近似的精度高达0.008%。优秀的协议的近似与精确的频率和周期解是演示了几个振荡器的参数值。

1。介绍

大量的工作一直致力于非线性问题领域的应用数学,物理,工程科学。一般而言,分析近似给定的非线性问题的解决方案是困难的,有时是不可能的;大量的技术实现了基于数值方法。其中包括变分迭代法(1,2),谐波平衡法(3,4),和能量平衡方法(5- - - - - -8)开发来解决非线性微分方程。在这项研究中,我们调查了应用高阶能量平衡方法cubic-quintic杜芬振荡器。非线性频率计算第一、二、三阶循证医学的结果并与不同的技术。

2。能量平衡方法的基本思想

本节简要介绍了能量平衡方法提出他(5]。在这种方法中,建立了振动的变分原理,那么相应的角频率的哈密顿被认为是可以很容易地通过伽辽金方法。

让我们考虑一般振荡器的运动的初始条件的形式 + ( ) = 0 , ( 0 ) = , ( 0 ) = 0 , ( 2 1 ) 在哪里 是初始振幅。

其变分可以写成 ( ) = 0 / 4 1 2 2 + ( ) ( 2 2 ) 在这里 = 2 / 非线性振荡的周期和吗 ( ) = ( )

的哈密顿(2.1可以书面形式): 1 ( ) = 2 2 + ( ) ( 2 3 ) 在(2.2)的动能 ( ) 和潜在的能量 ( ) 可以分别表示为 = ( 1 / 2 ) 2 , = ( )

以来整个振动系统是保守的,在运动的总能量保持不变;振荡器的哈密顿成为一个恒定值, = + = ( ) ( 2 4 ) 对于一阶近似,以下测试函数可以认为: ( ) = c o 年代 ( 2 5 ) 用(2.5)(2.3)收益率以下剩余: ( ) = 2 2 2 年代 n 2 + ( c o 年代 ) ( ) ( 2 6 ) 残被迫零,在平均意义上,通过设置加权剩余的积分为零 0 / 4 ( ) = 0 , = 1 , 2 , , ( 2 7 ) 在哪里 是一组权重函数(或试验)。

有很多的加权函数,也就是说,加勒金,最小二乘法,搭配等等。在这项研究中,我们利用伽辽金方法作为权重函数。

3所示。高阶能量平衡方法

为了延长他的能量平衡方法,让我们假设的解决方案(2.1)可以表示为 = 1 c o 年代 + 2 c o 年代 3 + + c o 年代 ( 2 1 ) ( 3 1 ) 从初始条件,系数应满足以下约束: = 1 + 2 + + ( 3 2 一个 ) 可以选择其中一个参数作为参数的依赖。因此, = 1 2 1 ( 3 2 b ) 通过插入(3所示。1)(2.7),可以获得以下系统: 0 / 4 ( ) c o 年代 = 0 , 0 / 4 ( ) c o 年代 3 = 0 , 0 / 4 ( ) c o 年代 ( 2 1 ) = 0 ( 3 3 )

4所示。例子

一个cubic-quintic杜芬振荡器。在下面几节中,非线性频率将与不同的技术的结果说明了能量平衡方法的效率和精度。

这种振荡器的控制微分方程的形式 + + 3 + 5 = 0 , w h e r e 0 ( 4 1 一个 ) 与初始条件 ( 0 ) = , ( 0 ) = 0 ( 4 1 b ) 的哈密顿(4.1)给出如下: 1 ( ) = 2 2 + 1 2 2 + 1 4 4 + 1 6 6 = 1 2 2 + 1 4 4 + 1 6 6 ( 4 2 ) 一阶近似的假设 ( ) 在以下形式: ( ) = c o 年代 ( 4 3 ) 用第一个近似成(4所示。2)的收益率 1 1 ( ) = 2 2 2 年代 n 2 1 ( ) + 4 4 c o 年代 4 1 ( ) + 6 6 c o 年代 6 1 ( ) + 2 2 c o 年代 2 ( ) 4 4 2 2 6 4 ( 4 4 ) 可以通过设置一阶近似 0 / 4 1 ( ) c o 年代 ( ) = 0 , = 2 ( 4 5 ) 第一幅频关系近似得到 ( ) = 7 0 + 4 9 2 + 3 8 4 7 0 ( 4 6 ) 获得更准确的结果,让我们定义 如下: ( ) = 1 c o 年代 + 1 c o 年代 3 ( 4 7 ) 用(4所示。7)(4所示。2以下剩余)结果: 2 ( ) = 2 2 1 年代 n ( ) + 3 1 年代 n ( 3 ) 2 + 2 1 c o 年代 ( ) + 1 c o 年代 ( 3 ) 2 + 4 1 c o 年代 ( ) + 1 c o 年代 ( 3 ) 4 + 6 1 c o 年代 ( ) + 1 c o 年代 ( 3 ) 6 2 2 4 4 6 6 ( 4 8 ) 我们设置 0 / 4 2 ( ) c o 年代 ( ) = 6 5 5 0 4 1 0 6 1 4 0 9 6 1 2 8 5 1 + 4 8 8 9 8 3 2 7 + 5 9 1 5 9 0 4 0 2 1 1 2 8 3 1 3 1 4 6 1 1 + 9 7 9 5 2 0 3 1 + 3 2 2 1 8 3 1 4 0 2 + 7 2 3 5 2 0 2 1 + 9 4 6 1 7 6 4 1 + 3 1 8 7 0 4 1 2 4 8 1 4 6 1 8 9 0 + 1 0 3 3 6 0 0 2 1 + 2 1 5 0 4 0 0 4 1 + 4 5 7 2 7 1 1 2 + 2 2 7 0 6 6 9 1 7 + 2 2 6 9 7 8 5 6 2 1 + 7 7 2 7 1 0 4 0 4 1 + 6 3 6 0 2 2 5 3 2 = 0 , ( 4 9 一个 ) 0 / 4 2 ( ) c o 年代 ( 3 ) = 5 2 6 5 5 4 6 6 + 8 3 2 9 1 9 0 4 5 1 + 4 6 7 8 9 7 8 3 3 1 3 9 5 0 7 2 0 2 1 + 2 1 7 6 3 4 5 4 2 9 1 + 3 0 1 1 2 0 3 1 + 6 2 1 6 1 6 6 1 5 3 7 3 3 3 6 3 2 2 1 1 2 8 9 2 1 6 0 0 4 1 1 4 5 4 9 5 3 5 2 9 6 2 1 1 0 1 6 1 5 8 + 1 0 1 2 9 2 8 2 1 + 1 3 9 6 7 3 6 4 1 + 1 8 0 1 3 7 1 2 + 4 8 1 2 6 7 8 9 6 2 + 5 0 0 2 6 2 4 2 1 + 1 0 3 0 9 6 3 2 4 1 + 4 5 7 2 7 1 1 2 = 0 ( 4 9 b )

通过求解(4.9)- (4.9 b)的同时,可以获得二阶近似幅频关系。为不同的 2 值,给出了近似频率表1

此外,结果将进一步提高通过定义的准确性 在以下形式: ( ) = 1 c o 年代 + 2 c o 年代 3 + 1 2 c o 年代 5 ( 4 1 0 ) 用(4.10)(4所示。2),我们得到以下剩余的三阶近似: 3 1 ( ) = 2 2 1 年代 n ( t ) + 3 2 年代 n ( 3 ) + 5 1 2 年代 n ( 5 ) 2 + 1 2 1 c o 年代 ( ) + 2 c o 年代 ( 3 ) + 1 2 c o 年代 ( 5 ) 2 + 1 4 1 c o 年代 ( ) + 2 c o 年代 ( 3 ) + 1 2 c o 年代 ( 5 ) 4 + 1 6 1 c o 年代 ( ) + 2 c o 年代 ( 3 ) + 1 2 c o 年代 ( 5 ) 6 2 2 4 4 6 6 ( 4 1 1 ) 插入(4.11)(3所示。3) = 3 上面解释和使用相同的过程,我们得到三个加权积分。同时解决这三个方程,三阶近似的幅频关系。对于高阶近似,可以应用类似的程序,然而,三阶近似的准确性是合适的几个参数的值 , ,

在下面,cubic-quintic振荡器的非线性频率两种不同情况下的计算:(i) = 0 和(2) 0

第一种情况考虑对应于立方杜芬振荡器。为 = 0 非线性频率得到的结果,第一,二,三阶能量平衡方法并与结果(4,9]。在第二种情况下, 0 ,给出了非线性频率表2。此外,所有情况下的数值解是通过标准龙格-库塔法(实际)。

在表1一阶近似读,相对误差为1.25%,虽然这在第二近似误差降低到0.59%。我们观察到三阶之间的差异和精确的频率是足够小。

见表2,三阶近似的结果很好的数值结果之间的协议。近似的比较和数值解也可以在图中找到1。可以看出,第一,和二阶结果略有差别相比,数值解。然而,三阶近似数值解的重叠。

5。结论

在本文中,延长高阶能量平衡方法的解决方案。达芬的一阶近似频率振荡器提供1.25%的相对误差,而第二个,和三阶近似频率相对误差达到0.59%和0.008%,分别。此外,高阶能量平衡的相对误差减小到较小的值比全局误差最小化和谐波平衡方法。因此,我们可以,扩展方法非常有效且方便的cubic-quintic杜芬振荡器。