二阶时滞方程的调和解稳定性分析.docVIP

二阶时滞方程的调和解稳定性分析 　　引 言 　　对于含有参数激励的非线性动力系统是振动理论中的一类经典系统，对它的研究是极其丰富而又复杂的动力学行为，对此已经有了大量的研究[9].然而，对于含有参数激励的非线性时滞动力系统的研究，相关的文献较少，仅限于文献[10][1]等少量工作. 　　本文用某些承受周期激励的地震波和具有时滞弹性地基作用的结构为背景[2]，因为地基具有一定的长度，当弹性力在特定的时刻作用在物体上时没有马上导致物体活动状态，而是必须经过固定的时间距离，在物体的加速度达到均衡状态的时候，物体初步活动.在这个阶段一定包含着时间滞后，同时，在处理相关实际难点的体系开展探索时，滞量是必须重视的.所以，去掉滞量就不能使精确度准确，严重的能引起体系紊乱，总而言之，最后的数学模型为： 　　x??-（α-βx?2）x?+x（t-τ）+γx3（t-τ）+（KcosΩt）x=0（1） 　　其中K=εL，L=1，2，3，…，ε是一个小参数.Ω≈2ω，t为时间，τ为时滞量，t0，τ0. 　　当γ=ε=0时，系统（1）的Hopf分支以及稳定性分析已经在文献[8]中进行了研究，本文主要是研究在该系统存在Hopf分支的情况下，外加周期参数扰动时，调和解是否稳定的问题.本文主要运用多尺度分析的方法讨论（1）中介绍的地震波方程在周期参数扰动下，仍然存在稳定的周期解的问题，并给出解的近似表达式. 　　一、调和解的存在性 　　对于系统（1）存在调和解： 　　定理1.1 　　??.当τγ1-β10，ατ-20，1-2σω0 时，系统（1.1）存在调和解， 　　此时a=a01=16?-（τω2+α）ω2[（β1-τγ1）ω2+（β1ω2+αγ1）（1-ατ）]0，b=b′010. 　　??.当ατ-1=0且τγ1-β1=0，且（1-2σω）（τα-2）同号时，系统（1.1）存在调和解， 　　此时a=a02=（τα-2）ω-σ[（τω2+α）2+ω2（τα-2）2]6[τβ1ω5+（2γ1+αβ1）ω3+α2γ1ω]0，b=b′020. 　　??.当1-2σω0，τγ1-β10时，系统（1.1）存在调和解， 　　此时a=a03=a010，b=b03=b′010. 　　二、调和解的稳定性定理 　　记f（a，b）=ma+（n+σ2）b+6（-mQ-nP）a2b+6（mP-nQ）ab2+6（mP-nQ）a3+6（-nP-mQ）b3， 　　g（a，b）=（n-σ2）a-mb+6（mP-nQ）a2b+6（mQ+nP）ab2+6（nP+mQ）a3+6（mP-nQ）b3. 　　则有： 　　f（a，b）a=m+12（-mQ-nP）ab+6（mP-nQ）b2+18（mP-nQ）a2， 　　f（a，b）b=（n+σ2）+6（-mQ-nP）a2+12（mP-nQ）ab+18（-nP-mQ）b2， 　　g（a，b）a=（n-σ2）+12（mP-nQ）ab+6（mQ+nP）b2+18（nP+mQ）a2， 　　g（a，b）b=-m+6（mP-nQ）a2+12（mQ+nP）ab+18（mP-nQ）b2. 　　当a=a0，b=b=0时， 　　f（a，b）aa=a* 　　b=b*=m+18（mP-nQ）a2， 　　f（a，b）ba=a* 　　b=b*=（n+σ2）+6（-mQ-nP）a2， 　　g（a，b）aa=a* 　　b=b*=（n-σ2）+18（nP+mQ）a2， 　　g（a，b）ba=a* 　　b=b*=-m+6（mP-nQ）a2. 　　得到雅可比阵记为A=m+18（mP-nQ）a2 （n+σ2）+6（-mQ-nP）a2（n-σ2）+18（nP+mQ）a2 -m+6（mP-nQ）a2 　　其特征方程为A-λE=m+18（mP-nQ）a2-λ （n+σ2）+6（-mQ-nP）a2（n-σ2）+18（nP+mQ）a2 -m+6（mP-nQ）a2-λ=[m+18（mP-nQ）a2-λ][-m+6（mP-nQ）a2-λ]-[（n+σ2）+6（-mQ-nP）a2][（n-σ2）+18（nP+mQ）a2]=λ2-[24（mP-nQ）a2]λ+[12n2P-12m2P+36mnQ+12σnP+6σmQ]a2+108（mP-nQ）a4+108（-mQ-nP）（nP+mQ）a4-m2-n2-σ24. 　　由一元二次方程求根公式λ=-b±b2-4ac2a， 　　其中a=1， 　　b=-24（mP-nQ）a2， 　　c=[12n2P-12m2P+36mnQ+12σnP+6σmQ]a2+108（mP-nQ）a4+108（-mQ-nP）（nP+mQ）a4-m2-n2-σ24. 　　当b=-24（mP-nQ）a20，即mP