路面激扰与电磁激扰下电动轮车辆系统非线性全局动态特性研究.pdf

第 卷第期 年月动 力 学 与 控 制 学 报J OUR NA LO FD YNAM I C SAN DC ON T R O LV o l N o J u n 文章编号:()D O I:/收到第稿,收到修改稿甘肃省教育厅 年度甘肃省高等学校创新基金项目(B )和甘肃省青年科技计划资助项目(J R R A ),G a n s uP r o v i n c i a lD e p a r t m e n to fE d u c a t i o n G a n s uP r o v i n c i a lH i g h e rE d u c a t i o nI n n o v a t i o nF u n dP r o j e c t(B )a n dG a n s uP r o v i n c i a lY o u t hS c i e n c ea n dT e c h n o l o g yP r o j e c t(J R RA )通信作者E m a i l:c l x c o m路面激扰与电磁激扰下电动轮车辆系统非线性全局动态特性研究陈丽霞 石建飞(兰州石化职业技术大学 汽车工程学院,兰州 )(兰州理工大学 机电工程学院,兰州 )摘要为研究电动轮车辆系统在路面电磁联合激励下的非线性振动特性,并挖掘参数与初始条件协同作用下系统全局动力学信息,建立了计及悬架系统和轮胎的非线性弹簧力和阻尼力的电动轮车辆系统两自由度/垂向振动非光滑非线性时变动力学模型,考虑外部激励的谐波性和随机性以及电磁激励的分段周期性,得到了含谐波性、随机性和周期性的复杂外激励模型,分析了驱动电机转速与多初值的关联性对系统动力学共存行为及其分岔的影响规律,研究了共存运动的多样性及其吸引域的拓扑构型,揭示了转速对系统全局动力学稳定性的影响机理结果表明:电机转速与多初值协同作用诱发系统出现多样性的共存运动,共存运动吸引域的拓扑结构复杂多样,部分共存运动的振动幅值可能较大亦可能很小,其全局动力学稳定性差异较大电机转速与初值的协同作用对电动轮车辆系统非线性振动和全局动态特性的影响不予忽视本文研究对电动轮车辆系统的动态性能改善和结构优化具有重要的理论意义和工程价值关键词电动轮车辆振动,路面电磁激励,非线性振动,吸引域,全局动力学特性中图分类号:O 文献标志码:AN o n l i n e a rG l o b a lD y n a m i cC h a r a c t e r i s t i c so fE l e c t r i cW h e e l e dV e h i c l eS y s t e m sw i t hP a v e m e n t a n dM a g n e t i cE x c i t a t i o n sC h e nL i x i a S h i J i a n f e i(S c h o o l o fA u t o m o t i v eE n g i n e e r i n g,L a n z h o uP e t r o c h e m i c a lU n i v e r s i t yo fV o c a t i o n a lT e c h n o l o g y,L a n z h o u ,C h i n a)(S c h o o l o fM e c h a n i c a l a n dE l e c t r o n i c a lE n g i n e e r i n g,L a n z h o uU n i v e r s i t yo fT e c h n o l o g y,L a n z h o u ,C h i n a)A b s t r a c t I no r d e rt os t u d yt h en o n l i n e a rv i b r a t i o no ft h ee l e c t r i cw h e e l e dv e h i c l es y s t e m u n d e rt h ep a v e m e n t m a g n e t i ce x c i t a t i o na n dt oe x c a v a t et h eg l o b a ld y n a m i ci n f o r m a t i o no ft h es y s t e mu n d e rt h es y n e r g i s t i ce f f e c to f p a r a m e t e r s a n d i n i t i a l c o n d i t i o n s,an e wn o n s m o o t hn o n l i n e a r t i m e v a r y i n gd y n a m i cm o d e l o f t h e t w o d e g r e e o f f r e e d o m/v e r t i c a l v i b r a t i o no f a ne l e c t r i cw h e e l e dv e h i c l e s y s t e mi n c l u d i n gt h en o n l i n e a rs p r i n gf o r c ea n dd a m p i n go f t h es u s p e n s i o ns y s t e ma n dt i r e sa r ee s t a b l i s h e d C o m p l e xe x t e r n a l e x c i t a t i o nm o d e l sw i t hh a r m o n i c i t y,r a n d o m n e s sa n dp e r i o d i c i t ya r eo b t a i n e db yc o n s i d e r i n gt h eh a r m o n i c i t ya n dr a n d o m n e s so fr o a de x c i t a t i o na n dt h ep i e c e w i s ep e r i o d i c i t yo fe l e c t r o m a g n e t i ce x c i t a t i o n T h e i n f l u e n c eo f t h ec o r r e l a t i o nb e t w e e nm o t o rs p e e da n dm u l t i p l e i n i t i a l v a l u e so nt h ec o e x i s t e n c eb e h a v i o ro f t h e s y s t e ma n d i t sb i f u r c a t i o n i s a n a l y z e d T h ed i v e r s i t yo f c o e x i s t i n gm o t i o n s a n d t h e t o p o l o g yo f t h e i rb a s i n s o f a t t r a c t i o na r e i n v e s t i g a t e d T h em e c h a n i s mo f t h e i n f l u e n c eo f r o t a t i o n a l s p e e do n t h eg l o b a ld y n a m i c s t a b i l i t yo f t h e s y s t e mi s r e v e a l e d T h e r e s u l t s s h o wt h a t t h e s y n e r g i s t i c e f f e c t o f t h em o t o rs p e e da n d t h em u l t i p l e i n i t i a l v a l u e s i n d u c e s av a r i e t yo f c o e x i s t i n gm o t i o n s i n t h e s y s t e m T h e t o p o l 第期陈丽霞等:路面激扰与电磁激扰下电动轮车辆系统非线性全局动态特性研究o g i c a l s t r u c t u r e so f t h eb a s i n so f a t t r a c t i o no f c o e x i s t i n gm o t i o n s a r e c o m p l e xa n dd i v e r s e,a n d t h ev i b r a t i o na m p l i t u d eo fs o m ec o e x i s t i n g m o t i o n sm a yb el a r g eo rs m a l l,a n dt h e i rg l o b a ld y n a m i cs t a b i l i t yc h a n g e sg r e a t l y T h e i m p a c to f t h es y n e r g i s t i ce f f e c to f t h em o t o r s p e e da n dt h e i n i t i a l v a l u eo nt h en o n l i n e a rv i b r a t i o na n dg l o b a l d y n a m i c so f t h es y s t e mc a n n o tb e i g n o r e d T h i s r e s e a r c hh a s i m p o r t a n t t h e o r e t i c a l s i g n i f i c a n c ea n de n g i n e e r i n gv a l u e f o r t h ed y n a m i cp e r f o r m a n c e i m p r o v e m e n ta n ds t r u c t u r a lo p t i m i z a t i o no f t h ee l e c t r i cw h e e l v e h i c l es y s t e mK e yw o r d s e l e c t r i cw h e e l e dv e h i c l ev i b r a t i o n,p a v e m e n t m a g n e t i c e x c i t a t i o n,n o n l i n e a r v i b r a t i o n,b a s i no f a t t r a c t i o n,g l o b a l d y n a m i c s引言电动轮汽车由于其结构灵活紧凑、动力源智能化、易操作控制等优点,已成为电动汽车的重要发展方向,轮毂电机驱动通常安装在车轮内部,在汽车行驶过程中会产生较大的电磁力,其在垂直方向的分量严重影响车辆的垂向振动,汽车悬架系统、轮胎等非线性使车辆系统呈现较强的非线性动力学行为,研究电动轮车辆系统非线性动力学特性、揭示其非线性振动机理,对其结构优化设计、动态性能改善、减振降噪及可靠性应用具有重要的理论意义,是电动轮汽车主动振动控制亟待解决的关键问题之一很多学者对汽车系统的动力学特性开展了大量的研究,取得了 很多有价值 的研究成果宋作军基于弹簧、减振器及轮胎的非线性方程,对汽车半主动悬架系统稳定性进行了分析,并研究了系统的H o p f分岔特性卢剑伟等考虑运动副间隙的影响,建立含间隙汽车摆振系统非线性动力学模型,研究了系统的非线性动态特性陈晓建等建立了考虑非线性弹簧力和非线性阻尼力的两自由度四分之一车辆动力学模型,分析了激励比和非线性刚度系数对悬架系统的动态特性的影响对于电动轮汽车而言,金智林等研究了轮毂电机驱动电动汽车的侧倾稳定性,并分析了不同路面激励下非簧载质量的影响近年来,杨蔚华等 通过推导电磁转矩和电机切向力的解析表达式,建立了路面和电磁双重激励下的/车辆动力学模型并分析了其平顺性,为路面电磁激励下车辆系统动力学研究奠定了基础李韶华及其团队 在路面电磁激励下电动轮系统的动力学建模、动态特性、平顺性评价指标等方面进行了开创性的研究,分析了电磁激励、路面二次激励、车速和车辆非线性对车辆平顺性和道路友好性的影响,研究了胎路多点接触、非线性地基、路面不平顺幅值等对车辆响应的影响规律,以及基于分形理论的轮胎与路面的三维动态接触行为随后,分析了刚柔耦合特种车辆通过凸台障碍路面时的动力学特性,并从车辆纵、横、垂向动力学三个方面,介绍了目前电动车辆动力学与控制领域的研究成果 张云等 建立了采用悬置轮毂式电机驱动的电动汽车四分之一主动悬架模型,并研究了其滑 模控制行为韩彦伟等 利用拉格朗日方程,研究了车辆行驶时的非线性动力学响应尽管考虑各种影响因素的电动轮汽车系统非线性动力学模型已被建立,其动态响应特性被广泛研究,但单参数影响下电动轮车辆系统非线性动力学响应的分岔与混沌特性却鲜有研究非线性系统的动力学响应不仅与参数相关,对初值亦具有较强的敏感依赖性,然而对于参数与初始条件协同作用下电动轮车辆系统的全局动力学特性的研究未见文献报道,系统非线性振动机理尚未清楚基于以上分析,本文考虑路面电磁联合激励的影响,考虑悬架系统和轮胎的非线性刚度和非线性阻尼,建立一个新的电动轮车辆系统两自由度/垂向振动非光滑非线性时变动力学模型,基于变步长四阶R u n g e K u t t a法求解模型,借助相图、P o i n c a r 映射图、多初值分岔图、最大L y a p u n o v指数(T o pL y a p u n o vE x p o n e n t,T L E)、胞映射等方法研究该系统在电机转速变化下的非线性全局动力学特性,为电动轮车辆系统动力学行为预测、评价及控制、减振降噪及可靠性应用奠定理论基础路面电磁激励下电动轮车辆非线性动力学建模考虑车辆垂向振动,假设车身结构和质量分配对称,路面电磁激励下/车辆的两自由度振动力学模型如图所示M和M为簧下和簧上质动力学与控制学报 年第 卷量,M为轮毂电机质量,MM为总非簧载质量,Fk()和Fc()为悬架系统的非线性弹簧力和阻尼力,表示时间,Fk()和Fc()为车辆轮胎的非线性弹簧力和阻尼力,Fd()为轮毂电机垂向激振力,X和X为簧下和簧上质点振动位移,X为谐波随机外激励轮胎的非线性弹簧力和非线性阻尼力表达式为:Fk()K(XX)K(XX)K(XX)()图谐波随机电磁激励下电动轮车辆力学模型F i g M e c h a n i c a lm o d e l o f e l e c t r i cw h e e l v e h i c l eu n d e rh a r m o n i c r a n d o m m a g n e t i ce x c i t a t i o nFc()C(XX)C(XX)(XX)C(XX)(XX)()式中,Ki和Ci(i,)分别表示轮胎的线性、平方和三次方非线性刚度和阻尼系数悬架系统的非线性弹簧力Fk()表示为:Fk()K(XX)K(XX)()其中,K为悬架系统刚度,为小参数,表示非线性程度悬架系统的非线性阻尼力Fc()由式()得到,Fc()Cs i g n(XX)(XX)()式中,C为悬架系统阻尼,为非对称性系数,s i g n()为符号函数,表征阻尼力的非光滑性,表达式为:s i g n(XX),XX,XX()考虑轮毂盘式永磁直流电机,其垂向激振力Fd()由式()计算得到 Fd()Te mRc o s(n),z T(z)T,(z)T(z)T()式中,R为电机定子的平均半径,n为电机转速,z为正整数,T/(n m)为激振力的作用周期,m为电机极数,Te m为电磁转矩,可由式()得到 Te m iBwN I(DoDi)()其中,i/为计算极弧系数,BT为气隙磁密幅值,w为绕组相数,N定子绕组每相线圈匝数,I为电流有效值,Do和Di为 永磁体的外径和内径考虑路面激励具有谐波性和随机性(高斯白噪声),则X表示为:Xc o s(rt)(i)()式中,为谐波激励幅值,r为谐波激励频率,为调节参数,(i)为随机变量,用M o n t e C a r l o数值模拟法表示为式()(i)/i I nis i n(i)()其中,H为白噪声的激励强度,i为计算步长,i和i为在,之间均匀变化的独立随机数根据牛顿运动定律,图所示的路面电磁联合激励下车辆系统两自由度/垂向振动非线性动力学模型如式()所示(MM)X Fk()Fc()Fk()Fc()Fd()MX Fk()Fc()()方程()()经合并、化简后带入上式,得(MM)XK(XX)K(XX)K(XX)C(XX)C(XX)(XX)C(XX)(XX)K(XX)K(XX)Cs i g n(XX)(XX)Fd()MXK(XX)K(XX)Cs i g n(XX)(XX)()第期陈丽霞等:路面激扰与电磁激扰下电动轮车辆系统非线性全局动态特性研究令初始位移D为长度尺度,时间尺度为T(MM)/K,定义以下无量纲参数:k K/(MM),k K D/(MM),Dk K D/(MM),C/(MM),t/T C D/(MM),/D,k K/(MM)C D/(MM),k KD/(MM)CD/(MM),fd(t)Fd()/D(MM)xX/D,xX/D,xX/D,nT/mM/(MM),/R D(MM),/D则方程()的无量纲形式为:x k(xx)k(xx)k(xx)(xx)(xx)(xx)(xx)(xx)k(xx)k(xx)s i g n(xx)(xx)fd(t)mx k(xx)k(xx)s i g n(xx)(xx)()式中,s i g n(xx),xx,xx()fd(t)iBwN I(DoDi)c o s(t),z Tt(z)T,(z)Tt(z)T()x(t)c o s(t)(i)()方程()是一个非光滑非线性系统,其中,x,fd(t)为系统输入激励,x,x,x,x为系统动态响应谐波随机电磁激励下电动轮车辆系统动力学模型是一个非光滑非线性系统,其实验验证非常困难本文模型是在文献 的基础上考虑系统外激励的谐波性和随机性,以及悬架系统的非线性弹性力和非线性阻尼力车辆系统外激励的谐波性和随机性是真实存在的,文献,分析了悬架系统的非线性弹簧力和非线性阻尼力,故本文建立的谐波随机电磁激励下车辆系统两自由度/垂向振动非线性动力学模型是正确的和有效的利用变步长四阶R u n g e K u t t a法数值求解方程()路面激励与电磁激励如图所示,无量纲谐波随机联合激励x(图 c)由谐波分量(图 b)和随机分量(图 a)构成系统量纲参数见表,令电磁驱动转速为 r/m i n,图为无量纲垂向电机电磁激励力fd(t)的时域图,fd(t)在两个时间程度上呈现周期性变化规律在小的时间尺度T 上,电磁力在fd(t)和fd(t)之间周期性切换;在大的时间尺度T上,电磁力在fd(t)和fd(t)之间周期性切换两时间尺度上电磁力的周期性变化使系统振动具有较强的时变特性图路面激励xF i g H a r m o n i c r a n d o me x c i t a t i o nx图电磁激励fd(t)F i g E l e c t r o m a g n e t i ce x c i t a t i o nfd(t)参数影响下系统多初值分岔特性方程()是一个含时变参数的非光滑非线性系统,其动力学特性不仅受系统参数的影响较大,而且对初始条件具有较强的敏感依赖性,为此研究动力学与控制学报 年第 卷多初值域内系统的全局动力学特性是非常有必要的考虑表系统量纲样本参数,经无量纲化后得系统无量纲参数为:k、k、k、k 、k、m,定义轮毂驱动电机无量纲转速为控制参数,考虑簧下质量MM的初始条件对系统全局动态特性的影响,令多初值域为(x,x,x,x)|x,x,xx,图为电机无量纲转速减小时系统簧下质量MM的振动位移x的多初值分岔图 图(a)和相应的T L E图 图(b),图(c)为图(a)中Z区域的局部放大图,相同参数下簧上质量的振动位移和速度的变化规律与图相似,仅是振动幅值有所差异,故本文只讨论簧下质量MM振动行为的多初值分岔及演化规律表系统样本参数T a b l e S y s t e ms a m p l ep a r a m e t e r sP a r a m e t e rV a l u eP a r a m e t e rV a l u ei/Ip ABTR mwT/nDo mM k gDi mM k gN M k gK N/mC Ns/mK N/mC Ns/mK N/mC Ns/mK N/mC Ns/m图(a)中,较大时(A点右侧),系统仅表现为稳定的周期运动P,对应T L E值小于零,见图(b),P 相轨迹如图(a),振动幅值较小,系统较稳定,图中符号“”表示P o i n c a r 映射图随着减小到A点,P 在部分初值下转迁为周期运动P,相应T L E值在零点附近发生跳跃,在另一部分初值下P 依然存在,故在A点左侧(A和B点之间)系统表现为P 和P 共存(双稳态),其相轨迹如图(b)所示,尽管P 与P 均为稳定周期运动,但P 的振动幅值远大于P,P 行为下系统动力学品质远高于P 因此,在相同参数域内可通过调节初始条件,使系统在稳定的动力学品质较好的P 上运行,以达到减振降噪、提高可靠性的目的随着进一步减小到B点,P 在部分初值域内再次经历极短暂的混沌行为后转迁为周期运动P,对应T L E值在B点附近出现大于零的情况,另一部分初值域内P 依然存在,故在B和C点之间,系统表现为P、P 和P 共存(三稳态),对应相轨迹如图(c)所示,图中符号“”表示P 的P o i n c a r 映射图相比之下,P 和P 的振动较大,而P 的振动较小,因此P 是系统期望的动力学行为在C点,P 在初值域内完全转迁为混沌行为CM,对应的T L E值大于零,随后在C和D点之间,系统表现为P、P 和混沌CM共存,相轨迹如图(d)所示,CM的振幅较大,P 的较小,而P 的最小,在图(d)近似为一个点(a)减小时系统多初值分岔图(a)M u l t i i n i t i a l b i f u r c a t i o nd i a g r a mo f t h es y s t e md e c r e a s i n g(b)与图(a)对应的T L E图(b)T L Em a pc o r r e s p o n d i n gt op l o t(a)(c)图(a)中Z区域局部放大图(c)E n l a r g e m e n to f t h eZr e g i o n i np l o t(a)图系统多初值分岔图与T L E图F i g B i f u r c a t i o nd i a g r a ma n dT L Ed i a g r a mo f t h es y s t e mw i t hm u l t i p l e i n i t i a l v a l u e s继续减小,在D点混沌运动CM在其初值域内完全退化为周期运动P,在D和E点之第期陈丽霞等:路面激扰与电磁激扰下电动轮车辆系统非线性全局动态特性研究间,系统表现为P、P 和P 共存(三稳态),相轨迹如图(e),P 的振动幅值最大,容易导致系统振动结构破坏,在工程中应尽可能避免在E点,P 完全转迁为P 运动(P 消失),在E和F点之间,系统表现为P 和P 共存,如图(f),P 的振动幅值依然较大在F点,P 在部分初值域内转迁为一个新的周期运动Q,在另一部分初值域内P 依然存在,故在F和G点之间系统表现为P、Q 和P 共存,相轨迹如图(g),P 和Q 相轨迹的拓扑结构不同,其中Q 的振动最大,P 次之,P 最小另外,随着的减小,P 的振动幅值不断地减小图多初值条件下系统相图和P o i n c a r 映射图,符号“”和“”表示P o i n c a r 映射点F i g P h a s ed i a g r a m sa n dP o i n c a r m a p p i n gd i a g r a m su n d e rm u l t i p l e i n i t i a l v a l u e s,t h es y m b o l sa n d i n d i c a t et h eP o i n c a r m a p p i n gp o i n t随着减小到G点,P 在初值域内经周期跳跃完全转迁为Q(P 消失),系统在G和H点之间表现为P 和Q 共存,相轨迹如图(h)所示,P 和Q 的振动随着的减小而减小在H点,P 退化为周期运动P,P 的振动明显增大,随后系统在H和I点之间表现为P 和Q 共存,但P 动力学与控制学报 年第 卷的振动幅值明显大于Q,如图(i)所示,此时Q 的动力学品质优于P 继续减小,如局部放大图(c),在I点Q 经周期跳跃转迁为较短暂的高周期运动,随后在J点进入混沌运动C,对应T L E值大于零,系统表为P 与混沌C 共存,如图(j)所示,尽管混沌C 不稳定,但其振动幅值较小,因此可根据实际需求,通过控制初始条件,使其在期望的轨道上运行在K点混沌C 退化为周期运动P,系统在K和L点之间表现为P 与Q 共存,相轨迹如图(k)所示在L点,Q 转迁为准周期运动P Z,其对应T L E值近似等于零,系统在L和M点之间表现为P 与准周期P Z共存,如图(l)在M点,准周期P Z退化为稳定的周期运动Q,在M和N点之间系统表现为P 和Q 共存,如图(m),Q 的振动相对较小在N点Q 转迁为混沌CM,随后在N和O点之间系统表现为P 和混沌CM共存,如图(n),相比之下,混沌行为的振幅较小,而P 的振幅较大在O点,混沌行为在初值域内完全退化为P,故系统在较小时(O点左侧),表现为一个稳定的周期运动P 综上,轮毂电机无量纲转速与多初值域的协同作用对系统多稳态行为的影响规律被完全揭示较大时系统表现为一个稳定的周期运动P,其振动幅值较小;随着减小,在初值域内,系统逐渐出现了多样性的共存运动,这些共存运动中仅有部分运动是稳定的、安全的共存运动诱发系统动力学行为在同一参数条件下,随着初值扰动发生转迁或者分岔,加剧系统振动甚至动态失稳在共存运动之间,一些运动振动较小,另一些运动振动较大或者不稳定,因此在不改变参数下可通过调节初始条件控制系统动力学行为,以期获得稳定、振动较小、动力学品质较优的动力学行为参数与多初值域的协同交互下,两种特殊的现象被揭示:)在G和A点之间,系统出现周期运动P 与多种类型的周期运动或混沌运动共存的情况 图(a),且P 运动始终是最稳定、振动最小、动力学品质最优的运动,可通过改变初始状态,使系统表现为期望的P 运动;)较小时(I点左侧),系统出现周期运动P 与多种类型的周期运动或混沌运动共存的情况,但P 的振动幅值较大,动力学品质较差,而共存的周期运动和混沌运动振动较小,动力学品质反而较好 图(c)另外,由于谐波随机分量的影响,系统运动轨迹在较小领域内具有发散性初值域内系统全局动力学行为及其演化规律由第节的分析可知,随着变化系统在多初值域内存在多样性的共存运动,这些共存运动的振动幅值或稳定性相差较大,为全面评价共存运动的多样性与全局特性,下面详细分析共存运动在初值域内吸引域的拓扑结构及其演化规律表不同参数范围内共存运动的数量及类型T a b l e N u m b e ra n dt y p e so f c o e x i s t i n gm o t i o n s i nd i f f e r e n tp a r a m e t e r r a n g e so fV a l u er a n g eo fC o e x i s t i n gm o t i o nn u m b e rC o e x i s t i n gm o t i o nt y p e sF o rF i g P P、P(a)P、P、P(b)P、P、CM(c)P、P、P(d)P、P(e)Q、P、P(f)Q、P(g)Q、P(h)P、C(i)P、Q(j)P、P Z(k)P、Q P、CM(l)P 图为变化时系统在内共存运动吸引域的拓扑结构及其演化过程,图为不同参数范围对应的多样性共存运动数量及其类型结合图,当较大时系统在多初值条件下仅表现为一个稳定的周期运动P,随着减小,系统出现稳定的P 与P 共存,其吸引域如图(a)所示,图中橙色区域为P 的初值域(吸引域),湛蓝色区域为P 的吸引域,相比之下,P 的吸引域的面积大于P 吸引域的面积,故P 的全局稳定性要强于P 随着减小,P 在其部分初值域内转迁为P(周期运动),随后出现P、P 和P 三稳态共存,其吸引域如图(b)所示,图中灰绿色表示P 的吸引域,第期陈丽霞等:路面激扰与电磁激扰下电动轮车辆系统非线性全局动态特性研究由于P 吸引域的面积最小,故P 的全局稳定性最弱,其初始条件的轻微扰动可能使系统运动状态由P 转迁为P 或P 随着进一步减小,P 在其初值域内完全转迁为混沌运动CM,随后系统出现P、P 和CM共存的情况,如图(c),尽管P 的吸引域面积有所增大,但其全局动力学稳定性依然是最弱的随着减小,混沌运动C M完全退化为周期运动P,故系统出现P、P 和P 共存情形,如图(d),随后P 完全转迁为P,系统表现为P 和P 共存,如图(e),由于P 吸引域面积较小,故P 的全局稳定性大于P,且随着减小,P 吸引域的面积不断缩小图变化时多初值域内系统共存运动的拓扑及其演化机制F i g T o p o l o g ya n de v o l u t i o no f c o e x i s t e n c em o t i o no f t h es y s t e mi nm u l t i p l e i n i t i a l v a l u ed o m a i n sw h e nc h a n g e s动力学与控制学报 年第 卷随着继续减小,P 在其部分初值域内转迁为另一种新的周期运动Q P 和Q 相轨迹的拓扑结构不同,见图(i),在内系统出现两个周期运动和一个周期运动共存,即P、Q 和P 共存,其吸引域如图(f)所示,图中天蓝色区域为Q 的吸引域,相比而言,P 吸引域的面积最大,Q 的次之,P 的最小,即P 的全局动力学稳定性最强,而P 的最弱随后,P 完全转迁为Q,系统表现为Q 和P 共存,且Q 的全局动力学稳定性较强,如图(g)所示随着减小,P 退化为P,系统表现为两个不同的周期运动共存(P 和Q 共存),如图(h)所示,粉色区域为P 的吸引域,且Q 的全局稳定在大于P 继续减小,P 运动持续存在,Q 运动先后转迁为混沌运动C、周期运动Q、准周期运动P Z以及混沌运动CM,故而系统在的不同范围内出现多样性的共存运动,如,P 和C 图(i)、P 和Q 图(j)、P 和P Z 图(k)以及P 和CM 图(l),且P 的吸引域面积(粉色区域)不断增大,故P 的全局稳定性增强,但P 的振动幅值较大对比图、图和图,在变化一定范围内,系统出现多样性的共存行为,当较大时,共存运动中P 持续存在且其吸引域面积较大,其全局动力学稳定性较强,由于P 的振动幅值较小,故系统表现为P 运动时,其全局动力学稳定性较好且振动小,是期望的动力学行为当较小时,P 运动吸引域面积不断增大,其全局动力学稳定性较强,故系统在大多初值下表现为P 运动,但其振动幅值较大,加剧系统振动,甚至引起系统结构破坏,在工程实际中应尽量避免此类运动可见,随着减小,多样性的共存运动吸引域的拓扑结构发生复杂的变化,吸引域面积或拓扑结构的演变引起系统共存运动类型的变化,某一共存运动吸引域面积的突然增大可能减小共存运动的数量,而某一共存运动吸引域面积的突然减少可能增加共存运动的数量考虑参数与多初值的协同作用,系统出现多样性的共存行为,其中有部分行为振动较小或全局动力学较稳定,从而为实现系统全局动力学行为的预测与评价以及期望运动的控制提供一定的理论和数据支撑结论建立了路面激扰与电磁激扰下电动轮车辆系统两自由度/垂向振动非线性动力学模型,采用定性与定量相结合的分析方法,如多初值分岔图、T L E、相图、P o i n c a r 映射图等,研究了多初值条件下系统动力学共存行为或多稳态随无量纲轮毂电机转速变化的演变过程,揭示了参数变化和初值扰动协同作用下动力学行为的演化规律,基于胞映射方法分析了多初值条件下共存行为的吸引域或初值域的拓扑特性,探讨了系统全局动力学稳定性研究结果为电动轮车辆的/垂向动力学分析提供理论参考价值具体结论如下:()建立的电动轮车辆系统垂向振动非线性动力学模型考虑了谐波、随机和电磁三重激励,电磁激励在时域内呈较强的非光滑非连续性,考虑了悬架系统和轮胎的非线性弹簧力和阻尼力,且非线性阻尼力具有分段特性因此,所建模型是一个含时变参数和随机参数的非光滑非线性动力学系统()电机转速与初值的协同作用对系统/垂向振动运动类型完整性及其分岔的影响不予忽视转速较大时,系统运动对初值扰动的影响不敏感,仅表现为一个稳定的周期运动且振幅较小;随着转速减小,系统运动受初值扰动的影响较大,出现了多样性的共存运动,且部分共存运动的振幅很大;转速较小时,系统运动不受初值扰动的影响,表现为稳定的周期运动且振幅较大()状态平面内系统共存运动的吸引域或初值域的拓扑结构复杂多样,多样性共存运动的全局动力学稳定性差异较大转速较大时,共存行为中周期运动的全局稳定性较大,且振动幅值较小,是车辆系统运行中期望的动力学行为;随着转速减小,另一个周期运动的全局动力学稳定性较好,但其振动幅值较大,为此应根据工程实际需求,通过合理地匹配转速与初值,使系统运行在期望的动态响应上本文的研究为车辆系统垂向振动设计、动力学行为预测与控制、动态性能改善提供理论依据参考文献Z HAN GXP,L I ANGYN,YUEH,e t a l R e v i e w第期陈丽霞等:路面激扰与电磁激扰下电动轮车辆系统非线性全局动态特性研究o f e l e c t r i cv e h i c l ep o l i c i e s i nC h i n a:c o n t e n t s u mm a r ya n de f f e c ta n a l y s i sJ R e n e w a b l ea n dS u s t a i n a b l eE n e r g yR e v i e w s,:李勇,徐兴,孙晓东,等轮毂电机驱动技术研究概况及发展综述 J电机与控制应用,():L IY,X U X,S UN X D,e ta l R e v i e wa n df u t u r ed e v e l o p m e n t o f i n w h e e lm o t o r d r i v e t e c h n o l o g yJE l e c t r i cM a c h i n e s&C o n t r o lA p p l i c a t i o n,():(i nC h i n e s e)钟银辉,李以农,杨超,等基于主动悬架控制轮边驱动电动车垂向振动研究 J振动与冲击,():Z HON GY H,L IY N,YAN G C,e ta l V e r t i c a lv i b r a t i o no f i n w h e e lm o t o r e