Аңшылықтың тербелісі - Hunting oscillation

Теміржол доңғалақ дөңгелектеріндегі аң аулау тербелісі

Аңшылықтың тербелісі Бұл өздігінен тербеліс, әдетте қажет емес, туралы тепе-теңдік.[1] Өрнек 19 ғасырда қолданысқа еніп, тепе-теңдік үшін жүйенің қалай «аң аулайтынын» сипаттайды.[1] Өрнек электроника, авиация, биология және теміржол техникасы сияқты әртүрлі салалардағы құбылыстарды сипаттау үшін қолданылады.[1]

Теміржол доңғалақтары

Классикалық аң аулау тербелісі - а-ның тербелісі теміржол көлік құралы (жиі аталады жүк көлігін аулау) туындаған конинг бағытталатын әрекет тұрақтылық туралы адгезиялық теміржол байланысты. Бұл өзара әрекеттесуінен туындайды адгезия күштер және инерциялық күштер. Төмен жылдамдықта адгезия басым болады, бірақ жылдамдық өскен сайын адгезия күштері мен инерция күштері шамасы мен салыстырмалы болады тербеліс сыни жылдамдықпен басталады. Бұл жылдамдықтан жоғары қозғалыс зорлық-зомбылыққа, жол мен доңғалаққа зақым келтіруі мүмкін және мүмкін рельстен шығу. Мәселе а бар жүйелерде пайда болмайды дифференциалды өйткені әрекет а-ның екі дөңгелегіне де байланысты доңғалақ дөңгелегі бірдей бұрыштық жылдамдықпен айналады, дегенмен дифференциалдар сирек кездеседі, ал әдеттегі пойыздардың орнына дөңгелектері осьтерге жұппен бекітілген. Кейбір пойыздар, сияқты Тальго 350, дифференциалды емес, бірақ оларға көбінесе аң аулау тербелісі әсер етпейді, өйткені олардың дөңгелектерінің көпшілігі бір-бірінен тәуелсіз айналады. Қуатты автокөліктің дөңгелектеріне аң аулау тербелісі әсер етуі мүмкін, өйткені электр машинасының дөңгелектері кәдімгі боттардағыдай осьтерге екі-екіден бекітілген. Бір-бірінен тәуелсіз айналатын және білікке жұппен бекітілмеген тәуелсіз дөңгелектермен жабдықталған конустық доңғалақтар мен доңғалақтар пойыз боттары үшін қолайлы дифференциалға қарағанда арзанырақ.[2]

Мәселе алғаш рет 19 ғасырдың аяғында байқалды, сол кезде пойыздардың жылдамдығы оған тап болатындай жоғары болды. Бұған қарсы тұру үшін ауыр күштер 1930 жылдары басталды, бұл ұзартылған жүк көліктері мен бүйірлік демпферді тудырды. ілулі ілгіш жүк көлігі. Жапондықтардың дамуында Шинкансен жүк көлігінің жылдамдығын 225 км / сағ-тан (140 миль / сағ) ұзарту үшін конустық емес дөңгелектер және басқа дизайн өзгерістері қолданылды. Еуропадағы және Жапониядағы ғылыми-зерттеу және тәжірибелік-конструкторлық жұмыстарға негізделген дөңгелектер мен жүк машиналарын жобалаудағы жетістіктер болат доңғалақ жүйелерінің жылдамдығын түпнұсқаға қол жеткізгеннен әлдеқайда ұзартты Шинкансен, артқы үйлесімділіктің артықшылығы мұндай технологияны, мысалы, альтернативаға қарағанда үстем етіп ұстайды автокөлік және маглев жүйелер. Болат доңғалақты пойыздардың жылдамдық рекордын француздар ұстайды TGV, 574,9 км / сағ (357 миль).

Кинематикалық талдау

Теміржол доңғалағының кониндік әрекеті

Сапалы сипаттама құбылысты біраз түсінуді қамтамасыз етсе, тереңірек түсіну үшін көлік құралын математикалық талдау қажет динамика. Онда да нәтижелер тек шамамен болуы мүмкін.

A кинематикалық сипаттамасында геометрия сілтемесіз қозғалыс күштер оны тудырады, сондықтан талдау тікелей жолда жүретін дөңгелектер жиынтығының геометриясын сипаттаудан басталады. Бастап Ньютонның екінші заңы күштерді байланыстырады үдеу денелер, әсер ететін күштер кинематикадан компоненттердің үдеуін есептеу арқылы алынуы мүмкін. Алайда, егер бұл күштер кинематикалық сипаттаманы өзгертсе (бұл жағдайда олар солай болса), онда нәтижелер тек дұрыс болуы мүмкін.

Болжамдар және математикалық емес сипаттама

Бұл кинематикалық сипаттама күштерді елемейтіндіктен бірқатар жеңілдетілген болжамдар жасайды. Біреу үшін бұл деп санайды жылжымалы кедергі нөлге тең. Доңғалақ дөңгелегі (а-ға бекітілмеген) пойыз немесе жүк көлігі ), түзу және түзу жолда алға итеру беріледі. Доңғалақ дөңгелегі жағалаудан басталады және ешқашан бәсеңдемейді, өйткені күш жоқ (доңғалақ дөңгелегіндегі жолға жабысып, сырғып кетпеуі үшін төмен бағытталған күштерден басқа). Егер бастапқыда доңғалақ дөңгелегі теміржол бойында орналасқан болса, онда әр дөңгелектің тиімді диаметрлері бірдей және дөңгелектер тізбегі жол бойымен мәңгілікке түзу айналады. Егер дөңгелектер жиынтығы тиімді диаметрлер (немесе радиустары) әр түрлі болатындай етіп центрден сәл алшақ болса, онда доңғалақ жиынтығы R радиусының қисығында қозғала бастайды (осы дөңгелектер жиынтығының радиусына байланысты және т.б.; Кейінірек шығарылуы керек) . Мәселе кинематикалық пайымдауды табу үшін қолданылады траектория дөңгелектер жиынтығының, дәлірек айтқанда траектория трассаның ортасына тігінен проекцияланған доңғалақ дөңгелегі центрінің. Бұл жер деңгейінің траекториясы және xy графикалық сызбасында кескінделген, мұндағы х - теміржол бойындағы қашықтық, ал y - «қадағалау қателігі», доңғалақ жиынтығы центрінің түзу сызығынан ауытқуы. жолдың ортасынан өтетін теміржол (екі рельстің ортасында).

Дөңгелектер траекториясының қисық жолмен жүретіндігін көрсету үшін тырнақты немесе бұранданы тегіс үстелдің үстіне қойып, итеріп жіберуге болады. Ол қисық шеңбер бойымен айналады, өйткені тырнақ немесе бұранда диаметрі әр түрлі дөңгелектері бар дөңгелектерге ұқсайды. Басы үлкен диаметрлі дөңгелекке ұқсайды, ал ұшы ұсақ диаметрлі дөңгелекке ұқсайды. Шегенің немесе бұранданың шеңбер бойымен айналуы кезінде (және одан да көп) теміржол доңғалақтары өзгеше жүреді, өйткені ол қисыққа бұрыла бастағанда, тиімді диаметрлер жолдың қисаюын төмендететін етіп өзгереді. «Радиус» пен «қисықтық» теміржол қисықтығына емес, дөңгелектер траекториясының қисықтығына сілтеме жасайтынын ескеріңіз, өйткені бұл өте дұрыс жол. Доңғалақ дөңгелегі дөңгелектелгенде, дөңгелектер олардың тиімді диаметрлері тең болғанша және жол қисық болмайынша қисықтық азаяды. Бірақ траекторияның осы сәтте көлбеуі болады (бұл жолдың орта сызығы бойынша диагональ бойынша өтетін түзу сызық), ол трассаның орталық сызығын артқа айналдырып, тиімді диаметрлерді кері айналдырады (бұрынғы кіші диаметрлі дөңгелек үлкен диаметрге айналады және керісінше). Бұл дөңгелектің қисық бойымен қарама-қарсы бағытта қозғалуына әкеледі. Ол қайтадан орталық сызықты асып түсіреді және бұл құбылыс доңғалақ дөңгелегі бір жағынан екінші жағынан тербелісімен шексіз жалғасады. Доңғалақ екенін ескеріңіз фланец рельспен ешқашан байланысқа түспейді. Бұл модельде рельстер әрдайым доңғалақ протекторымен рельстің басындағы сол сызық бойымен жанасады деп есептеледі, ол рельстер пышақтың шетінен болады және тек сызық бойымен дөңгелек протекторымен байланыс жасайды (ені нөлге тең).

Математикалық талдау

Пойыз доңғалақтың конустық формасының арқасында трассада қалады протекторлар. Егер доңғалақ дөңгелегі бір жағына «у» мөлшерімен ығыстырылса (іздеу қателігі), бір жағында рельспен жанасатын протектор радиусы азаяды, ал екінші жағында ол өседі. The бұрыштық жылдамдық екі доңғалақ үшін бірдей (олар a арқылы біріктірілген қатаң ось ), сондықтан үлкенірек диаметрі протектор жылдамдатады, ал кіші баяулайды. Дөңгелектер жиынтығы дөңгелектер мен дөңгелектер жиынтығының осіндегі жанасу нүктелері арқылы өтетін конустың генераторының қиылысуымен анықталған қисықтық центрінің айналасында жүреді. Қолдану ұқсас үшбұрыштар, бізде бұрылыс радиусы бар:

NewTurnRadius.png

мұндағы d - жол өлшеуіш, r түзу жүгірген кезде дөңгелектің радиусы, ал к - протектор конустық (бұл жолға перпендикуляр көлденең бағытта протектордың көлбеуі).

Түзу жолға қатысты доңғалақ жиынтығының жолы y (x) функциясымен анықталады, мұндағы х - жол бойындағы прогресс. Мұны кейде бақылау қателігі деп атайды.[3] Қозғалыс бағыты азды-көпті болған жағдайда параллель рельстерге, қисықтық жолдың екіншісіне қатысты болуы мүмкін туынды жол бойындағы қашықтыққа қатысты y шамасында [4]

Бұдан шығатыны траектория жол бойында теңдеу басқарылады:[5]

Бұл қарапайым гармоникалық қозғалыс толқын ұзындығы:

 Клингел формуласы ретінде белгілі (1883 ж. шыққан)[6]

Бұл кинематикалық талдау пойыздардың үнемі бір жақтан екінші жаққа ауытқып отыратынын білдіреді. Шын мәнінде, бұл тербеліс сөндірілген критикалық жылдамдықтан төмен және жүру сәйкесінше ыңғайлы. Кинематикалық нәтиже қозғалысты тудыратын күштерді елемейді. Оларды қолдану арқылы талдауға болады жорғалаушы ұғым (сызықтық емес), бірақ олардың мөлшерін анықтау оңай емес, өйткені олар серпімді бұрмалану байланыс аймақтарындағы доңғалақ пен рельстің. Бұл тақырып үйкеліс байланысының механикасы; Картер ұсынған аңшылық қозғалысын талдауда осы әсерлерді қамтитын ерте презентация.[7] Knothe қараңыз[8] тарихи шолу үшін.

Егер қозғалыс рельстермен айтарлықтай параллель болса, дөңгелектің бұрыштық орын ауыстыруы орнатылады береді:

Демек:

Бұрыштық ауытқу сонымен қатар циклдің төрттен бірінен екінші жағына қарай артта қалатын қарапайым гармоникалық қозғалыспен жүреді. Екі түрлі күйді қамтитын гармоникалық қозғалыспен сипатталатын көптеген жүйелерде (бұл жағдайда осьтің иілісі және бүйірлік ығысу) екі қозғалыс арасындағы ширек цикл кідірісі жүйені алға қарай қозғалудан энергия алу мүмкіндігімен қамтамасыз етеді. Бұл әсер «байқаладықыбырлау «ұшақтың қанаттары және»жылтыр «автомобильдер, сондай-ақ теміржол көлігін аулау. Жоғарыда алынған кинематикалық шешім критикалық жылдамдықтағы қозғалысты сипаттайды.

Іс жүзінде, критикалық жылдамдықтан төмен, екі қозғалыс арасындағы кідіріс тоқсандық циклден аз, сондықтан қозғалыс бәсеңдейді, ал критикалық жылдамдықтан жоғары болса, қозғалыс күшейетін етіп, тоқсандық циклдан үлкен болады.

Бағалау үшін инерциялық күштер, қашықтық туындыларын уақыт ретінде көрсету керек туындылар. Бұл U деп саналатын автомобильдің жылдамдығын пайдаланып жасалады:

Ойықтағы осьтің бұрыштық үдеуі:

Инерциялық момент (гироскопиялық әсерлерді ескермеу):

мұндағы F - рельстер бойымен әсер ететін күш, C - инерция моменті дөңгелектер жиынтығы.

максимум үйкелісті доңғалақ пен рельстің арасындағы күш:

мұндағы W - осьтің жүктемесі және болып табылады үйкеліс коэффициенті. Жалпы сырғанау жылдамдық пен осьтің ауытқуымен жүреді:

бұл өрнек критикалық жылдамдықты едәуір асыра бағалайды, бірақ аң аулаудың пайда болуының физикалық себебін көрсетеді, яғни инерциялық күштер белгілі бір жылдамдықтан жоғары адгезия күштерімен салыстырылады. Үйкелісті шектеу - бұл жағдайда адгезия күшінің нашар көрінісі.

Нақты адгезия күштері жанасу аймағында протектор мен рельстің бұрмалануынан пайда болады. Жалпы сырғанау жоқ, тек серпімді бұрмалану және жергілікті сырғанау (сырғанау сырғуы). Қалыпты жұмыс кезінде бұл күштер үйкелудің шектеулі шекарасында болады. Толық талдау осы күштерді қолдана отырып ескереді домалақ байланыс механикасы теориялар.

Алайда кинематикалық талдау доңғалақ рельсті байланыста сырғанау мүлдем болмады деп болжады. Енді дөңгелектер жиынтығының есептелген синусоидалық траекториясын (Клингел формуласы бойынша) дәл емес етіп жасайтын сырғудың сырғуы бар екендігі анық.

Энергия балансы

Критикалық жылдамдықтың бағасын алу үшін біз осы кинематикалық шешім дұрыс болатын жағдай тор жоқ жағдайға сәйкес келетіндігін қолданамыз энергия қоршаған ортамен алмасу, сондықтан қарастыру арқылы кинетикалық және потенциалды энергия жүйенің, біз критикалық жылдамдықты алуымыз керек.

Келіңіздер:

Операторды пайдалану:

бұрыштық үдеу теңдеуі бұрыштағы жылдамдықпен өрнектелуі мүмкін:

интеграциялау:

сондықтан айналуға байланысты кинетикалық энергия:

NewWheelSetYaw.png

Ось серпілген кезде, түйісу нүктелері протекторларда сыртқа қарай қозғалады, осылайша осьтің биіктігі төмендетіледі. Тірек нүктелерінің арақашықтығы:

(кішігірім шамалардың екінші ретінен) .тіректердің центрлерінен тіреу нүктесінің орын ауыстыруы:

осьтің жүктемесі түсіп кетеді

Осьтің жүктемесін төмендету арқылы орындалатын жұмыс:

Бұл жүйеден жоғалған энергия, сондықтан қозғалысты жалғастыру үшін доңғалақ парағының алға жылжуынан энергияның тең мөлшерін алу керек.

Сыртқы дөңгелектің жылдамдығы:

Кинетикалық энергия:

ішкі дөңгелегі үшін бұл

мұндағы m масса екі дөңгелектің де

Ұлғаюы кинетикалық энергия бұл:

Қозғалыс тұрақты болғанша тұрақты амплитудада жалғасады энергия алға жылжытудан алынған және нөлдік иіліске орнатылған доңғалақтың өскен кинетикалық энергиясы ретінде көрінетін, потенциалды энергия максималды иық кезінде біліктің жүктемесін төмендету арқылы жоғалтты.

Енді кинематикадан:

бірақ

Трансляциялық кинетикалық энергия

Толық кинетикалық энергия:

Критикалық жылдамдық энергия теңгерімінен анықталады:

Осыдан критикалық жылдамдық беріледі

Бұл доңғалақ конустықынан тәуелсіз, бірақ осьтің қатынасына байланысты жүктеме массаға дөңгелекке. Егер протекторлар конустық пішінде болса, критикалық жылдамдық конустыққа тәуелсіз болар еді. Іс жүзінде дөңгелектің тозуы конустың протекторлық ені бойынша өзгеруіне әкеледі, сондықтан потенциалдық энергияны анықтауға арналған конустың мәні кинетикалық энергияны есептегеннен өзгеше болады. Алғашқысын а деп белгілеп, критикалық жылдамдық келесідей болады:

мұндағы а - дөңгелекпен анықталатын формалық фактор кию. Бұл нәтиже алынған [9] стандартты қолдана отырып, жүйенің динамикасын талдаудан басқару инженері әдістер.

Оңайлатылған талдауды шектеу

Дөңгелектер жиынтығының қозғалысы осы талдау көрсеткеннен әлдеқайда күрделі. Автокөліктің тоқтата тұру күші бар қосымша күштер бар[10] және жоғары жылдамдықта дөңгелектер жиынтығы қосымша күш береді гироскопиялық критикалық жылдамдықтың бағасын өзгертетін моменттер. Әдетте теміржол көлігі төмен жылдамдықта тұрақты қозғалысқа ие, жоғары жылдамдыққа жеткенде тұрақтылық тұрақсыз қалыпқа ауысады. Теміржол көлігі жүйесінің динамикасын бейсызықтық талдаудың негізгі мақсаты - жанама жолдағы бифуркацияны, сызықтық емес бүйірлік тұрақтылықты және аңшылық мінез-құлықты аналитикалық зерттеу көзқарасын көрсету. Бұл зерттеуде талдау үшін Боголиубов әдісі бар[11]

Екі негізгі мәселе, яғни денені тұрақты тірек ретінде қабылдау және сызықтық емес элементтердің аң аулау жылдамдығын есептеу кезінде әсері көбінесе зерттеулерге бағытталған.[12] Нағыз теміржол көлігі көптеген еркіндік дәрежелеріне ие, демек, бірнеше жылдамдыққа ие болуы мүмкін; ең төменгісі доңғалақ қозғалысымен белгіленетіні сөзсіз.

Алайда, талдау ғибратты, өйткені аң аулаудың не үшін пайда болатынын көрсетеді. Жылдамдық өскен сайын инерциялық күштер адгезия күштерімен салыстырыла бастайды. Сондықтан критикалық жылдамдық осьтің жүктемесінің (адгезия күшін анықтайтын) дөңгелектер жиынтығының массасына (инерциялық күштерді анықтайтын) қатынасына байланысты болады.

Сонымен қатар, белгілі бір жылдамдықтан төмен, алға қозғалудан алынатын энергия осьтер мен қозғалыс демпстерін түсіру арқылы жоғалған энергияның орнын толтыру үшін жеткіліксіз; осы жылдамдықтан жоғары алынған энергия потенциалды энергияның жоғалуынан үлкен және амплитудасы жинақталады.

Максималды осьтің иілуіндегі потенциалдық энергияны серіппелі шектеуді енгізу арқылы арттыруға болады, осылайша серіппелі шиеленістен туындайтын үлес болады. Доңғалақ дөңгелектерінің қозғалысындағы шектеулерді арттыру үшін дөңгелектерді дөңгелектерге орналастыру және серіппелі шектеулерді креслоларға қолдану критикалық жылдамдықты жоғарылатады. Теңдеулерге серпімді күштерді енгізу суспензияға рұқсат береді, олар классикалық аң аулау емес, жалпы сырғанаудың басталуымен шектеледі. Аңшылықты виртуалды түрде жойғаны үшін төленетін айыппұл тікелей жол болып табылады, оған жолда жүру проблемасы және бұрынғы инфрақұрылымға сәйкес келмейді.

Аңшылық - бұл динамикалық мәселе, оны, кем дегенде, жолдың сапасына бейімделуі мүмкін кері байланысты бақылау арқылы шешуге болады. Алайда белсенді бақылауды енгізу сенімділік пен қауіпсіздік мәселелерін тудырады.

Аңшылық басталғаннан кейін көп ұзамай, жалпы тайғанақ пайда болады және доңғалақ фланецтері рельстерге әсер етіп, екеуіне де зақым келтіруі мүмкін.

Теміржол көлігі

Тәуелсіз рельсті доңғалақ осьтері автомобиль-рельсті көліктерде кең таралған

Көптеген теміржол көлігі рельстің әр дөңгелегіндегі тәуелсіз осьтер мен аспалы жүйелердің ерекшелігі. Бұл рельсте жол дөңгелектерінің болуымен үйлескенде жоғарыдағы формулаларды қолдану қиынға соғады. Тарихи тұрғыдан алғанда, автомобиль-теміржол көлігінің алдыңғы дөңгелектері сәл орнатылған саусақ көлік құралы рельсте жүргізіліп жатқан кезде аң аулауды азайтуға мүмкіндік беретіні анықталды.

Сондай-ақ қараңыз

Мәселенің осы класына қатысты жалпы әдістерді қараңыз

Пайдаланылған әдебиеттер

  1. ^ а б c Оксфорд ағылшын сөздігі (2-ші басылым). Оксфорд университетінің баспасы. 1989 ж. f. Машинаның, құралдың, жүйенің және т.б. аң аулайтын әрекеті (аң аулауды 7б. Қараңыз); тепе-теңдік жылдамдығы, жағдайы немесе күйі туралы жағымсыз тербеліс.
  2. ^ https://www.talgo.com/kz/rolling-stock/very-high-speed/350/
  3. ^ Егер дөңгелектердің жолы трассаның бойымен абсолютті түзу жүрсе және доңғалақ жұбы трек бойында орналасқан болса, «қадағалау қатесі» нөлге тең болады.
  4. ^ Қараңыз Қисықтық # Графиктің қисықтығы математикалық бөлшектер үшін. Шамамен теңдік теңдікке тең болады, егер бақылау қателігі, x-ге қатысты нөлдік көлбеу болса. Бақылау қателігі синус толқынына айналатындықтан, нөлдік көлбеу нүктелері максималды бақылау қателігінің нүктелерінде болады. Бірақ теңдік y шамасы төмен болған жағдайда шамамен дұрыс болады.
  5. ^ Ескертіп қой у оң және керісінше болған кезде теріс болады. R үшін басқа теңдеу, y теріс болған кезде дұрыс емес, өйткені R радиусы теріс болуына жол бермейді (математикалық анықтамаға сәйкес). Бірақ екі теңдеуді біріктіру арқылы R радиусы жойылғаннан кейін, алынған теңдеу екі жағдайды тексеру арқылы дұрыс болады: у теріс және у оң.
  6. ^ Ивички, 7-формула 2.1
  7. ^ Картер, Ф.В. (1928 ж. 25 шілде). «Локомотивтердің жұмысының тұрақтылығы туралы». Корольдік қоғамның еңбектері. А. 121 (788): 585–610. Бибкод:1928RSPSA.121..585C. дои:10.1098 / rspa.1928.0220.
  8. ^ Knothe, K. (2008). «Доңғалақ / рельсті байланыс механикасының тарихы: Редтенбахерден Калкерге дейін». Көлік жүйесінің динамикасы. 46 (1–2): 9–26. дои:10.1080/00423110701586469.
  9. ^ Уиккенс, А.Х. (1965–66). «Түзу жолдағы теміржол көлігінің динамикасы: бүйірлік тұрақтылықтың негізгі мәселелері». Proc. Инст. Мех. Eng.: 29–.
  10. ^ Уиккенс, А.Х .; Гилхрист А.О .; А.Е.В. Хоббс (1969–70). «Жоғары білікті екі осьті жүк көліктеріне арналған ілулі жобалау». Инженер-механиктер институтының материалдары: 22–.
  11. ^ Сераджиан, Реза (2013). «Богидің аңшылық мінез-құлқын сызықтық емес талдауға әр түрлі бүйірлік қаттылықпен параметрлердің өзгеретін әсері». Техника саласындағы өлшеу журналы: 195–206.
  12. ^ Сераджиан, Реза (2011). «Боги мен дене инерциясының hopf бифуркация теориясы мойындаған сызықтық емес дөңгелектегі аң аулауға әсері». Int J Auto Engng: 186–196.
  • Ивницки, Саймон (2006). Теміржол көлігі динамикасы туралы анықтама. CRC Press.
  • Шабана, Ахмед А .; т.б. (2008). Теміржол көлігінің динамикасы: есептеу тәсілі. CRC Press.
  • Уиккенс, A H (1 қаңтар, 2003). Теміржол көлігі динамикасының негіздері: бағыттаушы және көлденең тұрақтылық. Swets & Zeitlinger.
  • Сераджиан, Реза (2013). Богидің аң аулау мінез-құлқын сызықтық емес талдауға әр түрлі бүйірлік қаттылықпен параметрлердің өзгеретін әсері. CRC Press.
  • Сераджиан, Реза (2011). Боги мен дене инерциясының hopf бифуркация теориясы мойындаған сызықсыз дөңгелектегі аң аулауға әсері. CRC Press.