本文以在各类工程中广泛应用的充液系统为研究对象,建立了柱形储液箱中液体晃动模态的非线性动力学方程组。通过积分,得到了该类系统的同宿轨道的解析表达式,并利用Melnikov方法对此系统的全局动力学行为进行分析,计算出系统在扰动条件下进入混沌的阈值,表明此类系统在受到较大扰动时,将发生全局分岔和混沌现象。同时,文中给出的算例证明了理论分析结果的正确性。 The nonlinear governing equations of the liquid sloshing modals in cylindrical storage tank are established. Through integration, the analytical expressions of the homoclinic orbits of this kind of system are got. Using the Melnikov method, the global dynamical behaviors are analyzed. The threshold of the chaotic motion is calculated, which shows that there exists global bifurcation and chaos when the system is subjected to greater disturbance. Meanwhile, numerical simulations are also given, which confirm the analytical results.
曹俊灵1,2,钟顺2
1天津大学力学系,天津
2天津市联大通讯发展有限公司,天津
收稿日期:2018年5月9日;录用日期:2018年5月22日;发布日期:2018年5月29日
本文以在各类工程中广泛应用的充液系统为研究对象,建立了柱形储液箱中液体晃动模态的非线性动力学方程组。通过积分,得到了该类系统的同宿轨道的解析表达式,并利用Melnikov方法对此系统的全局动力学行为进行分析,计算出系统在扰动条件下进入混沌的阈值,表明此类系统在受到较大扰动时,将发生全局分岔和混沌现象。同时,文中给出的算例证明了理论分析结果的正确性。
关键词 :充液系统,同宿轨道,混沌阈值
Copyright © 2018 by authors and Hans Publishers Inc.
This work is licensed under the Creative Commons Attribution International License (CC BY).
http://creativecommons.org/licenses/by/4.0/
充液系统的非线性晃动动力学作为动力学与控制学科及其交叉学科的重要研究领域之一,广泛应用于工业、土木工程、汽车工程、船舶动力学、以及航空航天工业中 [
为建立有效的数学模型以研究充液系统的晃动动力学行为,尹立中 [
本文利用Melnikov方法对液体晃动方程的全局动力学行为 [
系统模型如图1所示,圆柱形储液箱半径为a,内部液体深度为h,定义OXYZ为绝对坐标系,局部坐标系 o θ r z 为柱坐标系。设液体燃料无旋,且与储液箱体接触角不变,根据速度势理论可得到液体晃动方程 [
∂ Φ ∂ t + 1 2 ( ∇ Φ − U ˙ − Ω × r ) 2 − 1 2 ( Ω × r ) 2 + Π ∓ 2 σ H ¯ ρ = 0 (1)
其中, Φ ( r , θ , z , t ) 为无旋流体的速度势函数, U ˙ 为储液箱平动速度矢量, Ω 为储液箱转动角速度, r 为转动半径, Π 为势函数, σ 为黏滞系数, ρ 为液体密度, H ¯ 为液面平均主曲率采用接触角不变假设时的常数,t为时间, ∇ 为散度算子。
不考虑储液箱体的转动,对液体的速度势函数进行模态展开 [
η ¨ m + 2 ς m η ˙ m + υ m 2 η m + ∑ i = 1 N ∑ j = 1 N σ i j m η i η j + ∑ i = 1 N ∑ j = 1 N ∑ k = 1 N σ i j k m η i η j η k + α ˜ m X ¨ + e ˜ m Y ¨ = 0 , m = 1 , 2 , 3 , ⋯ (2)
图1. 系统模型
图2. 液体晃动模态示意图
其中, η m ( m = 1 , 2 , 3 , ⋯ )为各阶广义模态位移, ς m ( m = 1 , 2 , 3 , ⋯ )为各阶模态阻尼, υ m ( m = 1 , 2 , 3 , ⋯ )为各阶模态频率, σ i j m 、 σ i j k m 、 α ˜ m 、 e ˜ m ,( m = 1 , 2 , 3 , ⋯ )为模态耦合系数。 X ¨ 和 Y ¨ 分别为X和Y方向上的扰动,上标一点为相应变量的一阶导数,二点为二阶导数。
设系统的扰动形式为 X ¨ = S 1 cos ( Ω 1 t ) , Y ¨ = 0 ,取液体晃动方程组中的第一阶主模态作为简化,得到
η ¨ 1 + υ 1 2 η 1 = K 31 η ˙ 1 + K 32 η 1 2 + K 33 η 1 3 + K 34 cos ( Ω 1 t ) (3)
其中, η 1 为液体晃动的第一阶主模态, υ 1 为一阶模态频率, S 1 为X方向扰动的振幅, Ω 1 为X方向扰动的频率, K 31 = − 2 ς 1 , K 32 = − σ 111 / μ 11 , K 33 = − σ 1111 / μ 11 和 K 34 = α ˜ 1 S 1 为简化后的方程系数。
为了便于讨论,将(3)改写为
f ¨ + λ 1 f + λ 2 f 2 + λ 3 f 3 = ε ( g cos ω t − δ f ˙ ) (4)
其中, f = η 1 , λ 1 = υ 1 2 , λ 2 = − K 32 , λ 3 = − K 33 , g = K 34 ε , δ = − K 31 ε , ω = Ω 1 。研究表明,在部分充液条件下,当充液深度相对较小时,上述各参数均为大于零的量 [
在方程(4)中 ε 为一小量,其对应的无扰动力系统为
{ y = f ˙ y ˙ = − λ 1 f − λ 2 f 2 − λ 3 f 3 (5)
其对时间的一次积分为
H = 1 2 f ˙ 2 + 1 2 λ 1 f 2 + 1 3 λ 2 f 3 + 1 4 λ 3 f 4 (6)
对于不同的H值,它代表相平面上不同的轨道。
系统(5)有3个不动点,它们分别为: O ( 0 , 0 ) , A ( − λ 2 2 λ 3 − λ 2 2 − 4 λ 1 λ 3 2 λ 3 , 0 ) 和 B ( − λ 2 2 λ 3 + λ 2 2 − 4 λ 1 λ 3 2 λ 3 , 0 ) ,不难判定O、A点是中心,B点是鞍点 [
现在来讨论无扰动动力系统的同宿运动,当 H = 0 ,由(6)式,得到
f ˙ * = ± − λ 1 f * 2 − 2 3 λ 2 f * 3 − 1 2 λ 3 f * 4 (7)
积分即得同宿轨道为:
− 1 − λ 1 ln | − λ 1 − 2 3 λ 2 f * − 1 2 λ 3 f * 2 + − λ 1 f * − λ 2 3 − λ 1 | = C 1 ± t (8)
或
| − λ 1 − 2 3 λ 2 f * − 1 2 λ 3 f * 2 + − λ 1 f * − λ 2 3 − λ 1 | = C 2 e ± − λ 1 t (9)
其中 C 1 和 C 2 为积分常量。
当 H < 0 时,得绕O、A两点的两族闭轨道,代表在稳定平衡位置附近的振动,当H增大时,振动的周期增大,当 H > 0 时,代表另一类型的周期性振动。
对动力系统(4),定义Melnikov函数为:
M ( τ ) = ∫ − ∞ + ∞ f ˙ * ϕ [ f * , f ˙ * , f ω ( t + τ ) ] d t (10)
其中, ϕ = ϕ ( f , f ˙ , ω t ) = g cos ω t − δ f ˙ 。
故
M ( τ ) = − δ ∫ − ∞ + ∞ f ˙ * 2 d t + g ∫ − ∞ + ∞ f ˙ * cos ω ( t + τ ) d t (11)
式(11)中的第一个积分为
∫ − ∞ + ∞ f ˙ * 2 d t = 4 3 λ 3 [ ( − λ 1 ) 3 2 + λ 2 2 − λ 1 3 λ 3 + λ 2 2 λ 2 2 − 9 λ 1 λ 3 9 λ 3 2 λ 3 arcsin 2 λ 2 2 λ 2 2 − 9 λ 1 λ 3 ] (12)
考虑式(11)中的第二个积分
∫ − ∞ + ∞ f ˙ * cos ω ( t + τ ) d t = ∫ − ∞ + ∞ cos ω ( t + τ ) d f * = 2 C 2 λ 1 ω ∫ − ∞ + ∞ e λ 1 t sin ω ( t + τ ) ( C 2 e λ 1 t + λ 2 3 λ 1 ) 2 + 1 2 λ 3 d t (13)
其中,
∫ − ∞ + ∞ e λ 1 t sin ω ( t + τ ) ( C 2 e λ 1 t + λ 2 3 λ 1 ) 2 + 1 2 λ 3 d t = 2 π ⋅ sin ω ( ξ + τ ) ( c h ω ψ 2 − c h ω ψ 1 ) C 2 2 λ 1 λ 3 ( 1 − c h 2 π ω λ 1 ) (14)
{ ξ = 1 λ 1 ln [ 1 3 C 2 2 λ 2 2 + 9 λ 1 λ 3 2 λ 1 ] ψ 1 = 1 λ 1 [ π + arctg ( 3 λ 2 λ 1 λ 3 2 ) ] ψ 2 = 1 λ 1 [ π − arctg ( 3 λ 2 λ 1 λ 3 2 ) ] (15)
将(12)和(13)代入式(11)中,得到Melnikov函数
M ( τ ) = − 4 δ 3 λ 3 [ λ 1 3 2 + λ 2 2 λ 1 3 λ 3 + λ 2 ⋅ 2 λ 2 2 + 9 λ 1 λ 3 9 λ 3 2 λ 3 arcsin 2 λ 2 2 λ 2 2 + 9 λ 1 λ 3 ] + 4 π g ω sin ω ( ξ + τ ) ( c h ω ψ 1 − c h ω ψ 2 ) 2 λ 3 ( c h 2 π ω λ 1 − 1 ) (16)
Melnikov函数如果有限,则稳定流形和不稳定流形相交,Poincare映射有马蹄,从而存在着奇怪的不变集合,对我们研究的耗散系统,就有可能出现混沌,由 M ( τ ) = 0 得到
sin ω ( ξ + τ ) = δ ( c h 2 π ω λ 1 − 1 ) 6 2 g ω λ 3 ( c h ω ψ 1 − c h ω ψ 2 ) [ λ 1 3 2 + λ 2 2 λ 1 3 λ 3 + λ 2 ⋅ 2 λ 2 2 + 9 λ 1 λ 3 9 λ 3 2 λ 3 arcsin 2 λ 2 2 λ 2 2 + 9 λ 1 λ 3 ] (17)
为使方程(17)有根,必须有 sin ω ( ξ + τ ) ≤ 1 ,故此,发生混沌的临界条件为
g δ ≥ ( c h 2 π ω λ 1 − 1 ) 6 2 ω λ 3 ( c h ω ψ 1 − c h ω ψ 2 ) [ λ 1 3 2 + λ 2 2 λ 1 3 λ 3 + λ 2 ⋅ 2 λ 2 2 + 9 λ 1 λ 3 9 λ 3 2 λ 3 arcsin 2 λ 2 2 λ 2 2 + 9 λ 1 λ 3 ] = R ( ω ) (18)
对非线性动力方程(4)进行数值计算,取 λ 1 = 40 ,
本文通过对液体晃动模态方程进行解析求解,得到了系统的同宿轨道的解析表达式,并利用Melnikov函数,确定了系统在受扰情况下发生混沌运动的阈值,表明系统在受到一定程度的扰动后,会发生同宿分岔和混沌运动。通过本文的结果,可以看到,全局动力学行为在液体晃动动力学分析中的重要性,局部分析结果往往只能对液体小幅晃动的情况进行预测和控制,当充液系统工作在较恶劣工况下,液体非小幅晃动时,系统的动力学行为更值得关注,本文的研究结果表明,当扰动较大时,系统发生混沌运动
图3. 分岔图
图4. g = 280时的各种曲线 (a) 时程曲线;(b) 相图;(c) Poincare 映射
的机理是同宿分岔,液体晃动虽非小幅,但为有限幅,在工程应用中,这一情况更容易得到控制,而不会对系统造成破坏性的影响,因此,此结果可为工程优化设计和控制提供理论依据。
国家自然科学基金青年项目(批准号:11502161)、天津市自然科学基金青年项目(合同号:13JCQNJC04300)。
曹俊灵,钟 顺. 充液系统中液体晃动的全局分岔和混沌 Global Bifurcation and Chaos of Liquid Sloshing in Liquid-Filled System[J]. 声学与振动, 2018, 06(02): 29-35. https://doi.org/10.12677/OJAV.2018.62004