root/applications/dual/VYZ/vyz_text.lyx @ 1332

Revision 1332, 84.2 kB (checked in by vahalam, 13 years ago)
Line 
1#LyX 1.6.7 created this file. For more info see http://www.lyx.org/
2\lyxformat 345
3\begin_document
4\begin_header
5\textclass scrreprt
6\begin_preamble
7\usepackage[czech]{babel}
8\end_preamble
9\use_default_options true
10\language czech
11\inputencoding auto
12\font_roman default
13\font_sans default
14\font_typewriter default
15\font_default_family default
16\font_sc false
17\font_osf false
18\font_sf_scale 100
19\font_tt_scale 100
20
21\graphics default
22\paperfontsize default
23\spacing single
24\use_hyperref false
25\papersize default
26\use_geometry false
27\use_amsmath 1
28\use_esint 1
29\cite_engine basic
30\use_bibtopic false
31\paperorientation portrait
32\secnumdepth 2
33\tocdepth 2
34\paragraph_separation indent
35\defskip medskip
36\quotes_language german
37\papercolumns 1
38\papersides 1
39\paperpagestyle default
40\tracking_changes false
41\output_changes false
42\author ""
43\author ""
44\end_header
45
46\begin_body
47
48\begin_layout Standard
49\align left
50\begin_inset ERT
51status open
52
53\begin_layout Plain Layout
54
55
56\backslash
57thispagestyle{empty}
58\end_layout
59
60\end_inset
61
62
63\end_layout
64
65\begin_layout Standard
66\align center
67
68\size large
69České vysoké učení technické v Praze
70\end_layout
71
72\begin_layout Standard
73\align center
74
75\size large
76Fakulta jaderná a fyzikálně inženýrská
77\end_layout
78
79\begin_layout Standard
80\begin_inset VSpace bigskip
81\end_inset
82
83
84\end_layout
85
86\begin_layout Standard
87\align center
88Katedra matematiky
89\end_layout
90
91\begin_layout Standard
92\align center
93Obor: Inženýrská informatika
94\end_layout
95
96\begin_layout Standard
97\align center
98Zaměření: Softwarové inženýrství a matematická informatika
99\end_layout
100
101\begin_layout Standard
102\begin_inset VSpace bigskip
103\end_inset
104
105
106\end_layout
107
108\begin_layout Standard
109\align center
110\begin_inset Graphics
111        filename /home/michal/Dokumenty/Bakalarka/moje/komplet/logo_cvut.eps
112        lyxscale 20
113        scale 20
114
115\end_inset
116
117
118\end_layout
119
120\begin_layout Standard
121\begin_inset VSpace bigskip
122\end_inset
123
124
125\end_layout
126
127\begin_layout Standard
128\align center
129
130\size larger
131\color black
132Metody duálního řízení elektrických pohonů
133\end_layout
134
135\begin_layout Standard
136\begin_inset VSpace smallskip
137\end_inset
138
139
140\end_layout
141
142\begin_layout Standard
143\align center
144
145\size larger
146\color black
147Dual control methods for electrical drives
148\end_layout
149
150\begin_layout Standard
151\begin_inset VSpace bigskip
152\end_inset
153
154
155\end_layout
156
157\begin_layout Standard
158\align center
159
160\size largest
161\color black
162VÝZKUMNÝ ÚKOL
163\end_layout
164
165\begin_layout Standard
166\begin_inset VSpace vfill
167\end_inset
168
169
170\end_layout
171
172\begin_layout Standard
173\align center
174Vypracoval: Michal Vahala
175\end_layout
176
177\begin_layout Standard
178\align center
179Vedoucí práce: Ing.
180 Václav Šmídl, Ph.D.
181\end_layout
182
183\begin_layout Standard
184\align center
185Rok: 2011
186\end_layout
187
188\begin_layout Standard
189\begin_inset Newpage newpage
190\end_inset
191
192
193\end_layout
194
195\begin_layout Standard
196\begin_inset ERT
197status open
198
199\begin_layout Plain Layout
200
201
202\backslash
203thispagestyle{empty}~
204\end_layout
205
206\end_inset
207
208
209\end_layout
210
211\begin_layout Standard
212\begin_inset VSpace vfill
213\end_inset
214
215
216\end_layout
217
218\begin_layout Subsubsection*
219Prohlášení
220\end_layout
221
222\begin_layout Standard
223Prohlašuji, že jsem výzkumný úkol vypracoval samostatně a použil jsem pouze
224 podklady uvedené v přiloženém seznamu.
225\end_layout
226
227\begin_layout Standard
228\begin_inset VSpace bigskip
229\end_inset
230
231
232\end_layout
233
234\begin_layout Standard
235\noindent
236\align left
237V Praze dne \SpecialChar \ldots{}
238\SpecialChar \ldots{}
239\SpecialChar \ldots{}
240\SpecialChar \ldots{}
241\SpecialChar \ldots{}
242 
243\begin_inset space \hfill{}
244\end_inset
245
246\SpecialChar \ldots{}
247\SpecialChar \ldots{}
248\SpecialChar \ldots{}
249\SpecialChar \ldots{}
250\SpecialChar \ldots{}
251\SpecialChar \ldots{}
252
253\end_layout
254
255\begin_layout Standard
256\noindent
257\align block
258\begin_inset space \hfill{}
259\end_inset
260
261Michal Vahala
262\begin_inset ERT
263status open
264
265\begin_layout Plain Layout
266
267~~
268\end_layout
269
270\end_inset
271
272
273\end_layout
274
275\begin_layout Standard
276\begin_inset Newpage newpage
277\end_inset
278
279
280\end_layout
281
282\begin_layout Standard
283\begin_inset ERT
284status open
285
286\begin_layout Plain Layout
287
288
289\backslash
290thispagestyle{empty}~
291\end_layout
292
293\end_inset
294
295
296\end_layout
297
298\begin_layout Standard
299\begin_inset VSpace vfill
300\end_inset
301
302
303\end_layout
304
305\begin_layout Subsubsection*
306Poděkování
307\end_layout
308
309\begin_layout Standard
310Především bych chtěl poděkovat
311\begin_inset VSpace defskip
312\end_inset
313
314
315\end_layout
316
317\begin_layout Standard
318\begin_inset space \hfill{}
319\end_inset
320
321Michal Vahala
322\end_layout
323
324\begin_layout Standard
325\begin_inset Newpage newpage
326\end_inset
327
328
329\end_layout
330
331\begin_layout Standard
332\begin_inset ERT
333status open
334
335\begin_layout Plain Layout
336
337
338\backslash
339thispagestyle{empty}
340\end_layout
341
342\end_inset
343
344
345\end_layout
346
347\begin_layout Description
348
349\emph on
350Název
351\begin_inset space \space{}
352\end_inset
353
354práce:
355\emph default
356\color black
357
358\begin_inset ERT
359status open
360
361\begin_layout Plain Layout
362
363~
364\end_layout
365
366\end_inset
367
368
369\begin_inset Newline newline
370\end_inset
371
372Metody duálního řízení elektrických pohonů
373\end_layout
374
375\begin_layout Description
376\begin_inset VSpace defskip
377\end_inset
378
379
380\end_layout
381
382\begin_layout Description
383
384\emph on
385Autor:
386\emph default
387 Michal Vahala
388\end_layout
389
390\begin_layout Description
391
392\emph on
393Obor:
394\emph default
395 Inženýrská informatika
396\end_layout
397
398\begin_layout Description
399
400\emph on
401Druh
402\begin_inset space \space{}
403\end_inset
404
405práce:
406\emph default
407 Výzkumný úkol
408\end_layout
409
410\begin_layout Description
411
412\emph on
413Vedoucí
414\begin_inset space \space{}
415\end_inset
416
417práce:
418\emph default
419 Ing.
420 Václav Šmídl, Ph.D.
421\end_layout
422
423\begin_layout Description
424
425\emph on
426Abstrakt:
427\emph default
428 
429\end_layout
430
431\begin_layout Description
432
433\emph on
434Klíčová
435\begin_inset space \space{}
436\end_inset
437
438slova:
439\emph default
440 
441\end_layout
442
443\begin_layout Standard
444\begin_inset VSpace bigskip
445\end_inset
446
447
448\end_layout
449
450\begin_layout Description
451
452\emph on
453Title:
454\emph default
455\color black
456
457\begin_inset ERT
458status open
459
460\begin_layout Plain Layout
461
462~
463\end_layout
464
465\end_inset
466
467
468\begin_inset Newline newline
469\end_inset
470
471Dual control methods for eletrical drives
472\end_layout
473
474\begin_layout Description
475\begin_inset VSpace defskip
476\end_inset
477
478
479\end_layout
480
481\begin_layout Description
482
483\emph on
484Author:
485\emph default
486 Michal Vahala
487\end_layout
488
489\begin_layout Description
490
491\emph on
492Abstract:
493\emph default
494 
495\end_layout
496
497\begin_layout Description
498
499\emph on
500Key
501\begin_inset space \space{}
502\end_inset
503
504words:
505\emph default
506 
507\end_layout
508
509\begin_layout Standard
510\begin_inset Newpage newpage
511\end_inset
512
513
514\end_layout
515
516\begin_layout Standard
517\begin_inset ERT
518status open
519
520\begin_layout Plain Layout
521
522
523\backslash
524thispagestyle{empty}
525\end_layout
526
527\end_inset
528
529
530\end_layout
531
532\begin_layout Standard
533\begin_inset CommandInset toc
534LatexCommand tableofcontents
535
536\end_inset
537
538
539\end_layout
540
541\begin_layout Standard
542\begin_inset Newpage newpage
543\end_inset
544
545
546\end_layout
547
548\begin_layout Standard
549\begin_inset ERT
550status open
551
552\begin_layout Plain Layout
553
554
555\backslash
556thispagestyle{empty}
557\end_layout
558
559\end_inset
560
561
562\end_layout
563
564\begin_layout Chapter*
565Seznam použitého označení a zkratek
566\end_layout
567
568\begin_layout Subsubsection*
569Zkratky
570\end_layout
571
572\begin_layout Labeling
573\labelwidthstring 00.00.0000
574
575\emph on
576PMSM
577\emph default
578 synchronní stroj s permanentními magnety (
579\emph on
580Permanent Magnet Synchronous Machine
581\emph default
582)
583\end_layout
584
585\begin_layout Labeling
586\labelwidthstring 00.00.0000
587
588\emph on
589SMPMSM
590\emph default
591 PMSM s magnety na povrchu rotoru (
592\emph on
593Surface Mounted PMSM
594\emph default
595)
596\end_layout
597
598\begin_layout Labeling
599\labelwidthstring 00.00.0000
600
601\emph on
602IPMSM
603\emph default
604 PMSM s magnety uvnitř rotoru (
605\emph on
606Inner PMSM
607\emph default
608)
609\end_layout
610
611\begin_layout Labeling
612\labelwidthstring 00.00.0000
613
614\emph on
615LQG
616\emph default
617 lineárně kvadraticky gaussovské řízení (
618\emph on
619Linear-Quadratic-Gaussian
620\emph default
621)
622\end_layout
623
624\begin_layout Labeling
625\labelwidthstring 00.00.0000
626
627\emph on
628PI
629\emph default
630proporcionálně integrační regulátor
631\end_layout
632
633\begin_layout Labeling
634\labelwidthstring 00.00.0000
635
636\emph on
637EKF
638\emph default
639rozšířený Kalmanův filtr (
640\emph on
641Extended Kalman Filter
642\emph default
643)
644\end_layout
645
646\begin_layout Subsubsection*
647Označení
648\end_layout
649
650\begin_layout Labeling
651\labelwidthstring 00.00.0000
652\begin_inset Formula $\hat{a}$
653\end_inset
654
655 značí odhad veličiny
656\begin_inset Formula $a$
657\end_inset
658
659
660\end_layout
661
662\begin_layout Labeling
663\labelwidthstring 00.00.0000
664\begin_inset Formula $j$
665\end_inset
666
667 komplexní jednotka
668\end_layout
669
670\begin_layout Standard
671\begin_inset Newpage newpage
672\end_inset
673
674
675\end_layout
676
677\begin_layout Addchap
678Úvod
679\end_layout
680
681\begin_layout Standard
682Hlavní náplní této práce je řízení elektrických pohonů, konkrétně synchronního
683 motoru s permanentními magnety (v textu bude označován zkratkou PMSM z
684 anglického
685\emph on
686Permanent Magnet Synchronous Machine
687\emph default
688).
689 Jedná se o synchronní stroj, tedy rotor se otáčí současně (synchronně)
690 s točivým magnetickým polem statoru.
691 Na rotoru má ale místo budícího vinutí permanentní magnety.
692 Tato konstrukce nachází v poslední době stále větší uplatnění.
693 Je tomu tak především z důvodu snadnější dostupnosti kvalitních permanentních
694 magnetů, ale také díky možnosti využít stále výkonější polovodičová zařízení
695 pro řízení a napájení těchto strojů.
696\end_layout
697
698\begin_layout Standard
699Jak se ale ukazuje, řízení takovýchto strojů, zjeména pokud se jedná o takzvaný
700 bezsenzorový návrh je netriviální.
701 Je tedy třeba hledat vhodné řídící algoritmy, které zvládnou motor efektivně
702 řídit i v bezsenzorovém případě a umožní širší nasazení PMSM v praxi.
703\end_layout
704
705\begin_layout Standard
706V tomto textu je nejdříve stručně popsán samotný PMSM, následuje odvození
707 rovnic popisující tento stroj v nejčastěji používaných souřadných soustavách.
708 Dále je formulována problematika estimace a určovaní stavových veličin,
709 kdy je kladen důraz na bezsenzorový případ.
710 Následuje popis nejčastěji použavaných řídících technik, které jsou současně
711 dostatečně jednoduché, aby mohly být teoreticky nasazeny i pro případ řízení
712 v reálném čase.
713 Zvláštní pozornost je věnována řízení označovanému jako LQG.
714 
715\end_layout
716
717\begin_layout Standard
718Dále se text věnuje duálnímu řízení, které se zdá být vhodným kandidátem
719 na zvládnutí úlohy řízení PMSM.
720 Protože je však problém duálního řízení obecně velmi složitá úloha, zaměříme
721 se na jeho nejjednodušší případy, které by mohly být nasazeny i v reálném
722 čase.
723\end_layout
724
725\begin_layout Standard
726Na závěr jsou prezentovány výsledky simulací a jsou navrženy směry a metody,
727 které by mohly vést k úspěšnému řešení problému.
728\end_layout
729
730\begin_layout Paragraph
731Poznámka
732\end_layout
733
734\begin_layout Standard
735V celém textu bude
736\begin_inset Formula $j$
737\end_inset
738
739 označovat komplexní jednotku
740\begin_inset Formula $j=\sqrt{-1}$
741\end_inset
742
743.
744 Označení
745\begin_inset Formula $i$
746\end_inset
747
748 bude obvykle značit elektrický proud, komplexní jednotku však nikdy.
749\end_layout
750
751\begin_layout Chapter
752Popis PMSM
753\end_layout
754
755\begin_layout Section
756Vlastnosti
757\end_layout
758
759\begin_layout Subsection
760Permanentní magnety
761\end_layout
762
763\begin_layout Standard
764Jak již bylo řečeno pro PMSM mají velký význam kvalitní permanentní magnety.
765 Podle
766\begin_inset CommandInset citation
767LatexCommand cite
768key "cdern2010,novak2006"
769
770\end_inset
771
772 jsou magnety vhodné pro PMSM vyráběny ze speciálních slitin nejčastěji
773 na bázi prvků
774\begin_inset Formula $Sm-Co$
775\end_inset
776
777 nebo
778\begin_inset Formula $Nd-Fe-B$
779\end_inset
780
781.
782 Oproti klasickým feritovým magnetům se vyznačují velkou magnetickou indukcí
783 okolo
784\begin_inset Formula $1T$
785\end_inset
786
787 oproti přibližne
788\begin_inset Formula $0,3T$
789\end_inset
790
791 u feritových magnetů.
792\end_layout
793
794\begin_layout Standard
795Nevýhodou nejen těchto, ale permanentních magnetů obecně je změna jejich
796 magnetických vlastností s teplotou.
797 Jedná se především o hranici označovanou jako
798\emph on
799Courieův bod
800\emph default
801, kdy materiál přechází z feromagnetického stavu do paramagnetického a s
802 tím je spojen výrazný pokles magnetizmu.
803 Tato hodnota závisí na použítém materiálu a pohybuje se přibližně v rozmezí
804 
805\begin_inset Formula $200-1000^{\circ}C$
806\end_inset
807
808.
809 Z toho vyplývá, že je nutné udržovat motor na vhodné provozní teplotě a
810 tedy zajistit odpovídající chlazení.
811\end_layout
812
813\begin_layout Subsection
814Výhody a nevýhody PMSM
815\end_layout
816
817\begin_layout Standard
818Následující část popisující výhody a nevýhody čerpá především ze zdrojů
819 
820\begin_inset CommandInset citation
821LatexCommand cite
822key "cdern2010,novak2006"
823
824\end_inset
825
826
827\end_layout
828
829\begin_layout Subsubsection
830Výhody
831\end_layout
832
833\begin_layout Standard
834Proč se PMSM využívají a jaké mají výhody oproti jiným motorům.
835 Uveďme především:
836\end_layout
837
838\begin_layout Itemize
839rotor neobsahuje vinutí a tedy
840\end_layout
841
842\begin_deeper
843\begin_layout Itemize
844je možno jej konstruovat menší, což je velmi výhodné v aplikacích, kde záleží
845 na co nejmenší velikosti pohonu, příkladem mohou být dopravní prostředky,
846 kde lze ušetřené místo využít například pro cestující (nízkopodlažní tramvaj)
847\end_layout
848
849\begin_layout Itemize
850je možno jej konstruovat lehčí, což snižuje hmotnost celého zařízení
851\end_layout
852
853\begin_layout Itemize
854má menší moment setrvačnosti rotoru
855\end_layout
856
857\begin_layout Itemize
858není třeba
859\emph on
860 
861\emph default
862složitě přivádět
863\emph on
864 
865\emph default
866napájení
867\emph on
868 
869\emph default
870na rotor
871\end_layout
872
873\end_deeper
874\begin_layout Itemize
875není třeba motor před rozběhem budit a nepotřebuje zdroj budícího proudu
876\end_layout
877
878\begin_layout Itemize
879odpadá problém s přívodem proudu do buzení rotoru
880\end_layout
881
882\begin_layout Itemize
883vyšší účinnost -- nejsou jouleovy ztráty v rotoru (oproti asynchronnímu
884 stroji) popřipadě v buzení (oproti synchronnímu stroji s buzením)
885\end_layout
886
887\begin_layout Itemize
888momentová přetížitelnost
889\end_layout
890
891\begin_layout Itemize
892možnost konstrukce pomaluběžného stroje s dostatečným výkonem, který nepotřebuje
893 převedovku (výhody spojené s absencí převodovky)
894\end_layout
895
896\begin_layout Subsubsection
897Nevýhody
898\end_layout
899
900\begin_layout Standard
901Na druhou stranu toto řešení motoru má i své nevýhody, jedná se zejména
902 o:
903\end_layout
904
905\begin_layout Itemize
906technologicky složitější výroba -- připevnění permanentních magnetů na rotor
907 (nejčastěji lepení)
908\end_layout
909
910\begin_layout Itemize
911složitější opravy
912\end_layout
913
914\begin_layout Itemize
915vyšší cena (nezanetbatelné náklady na permanentní magnety)
916\end_layout
917
918\begin_layout Itemize
919menší robustnost
920\end_layout
921
922\begin_layout Itemize
923problematické odbuzování
924\end_layout
925
926\begin_layout Itemize
927nutnost dobrého chlazení -- závislot magnetických vlastností permanentních
928 magnetů na teplotě
929\end_layout
930
931\begin_layout Itemize
932problematika spojená s návrhem řízení těchto strojů (bude detailněji rozebrána
933 níže)
934\end_layout
935
936\begin_layout Section
937Konstrukce
938\end_layout
939
940\begin_layout Standard
941\begin_inset Float figure
942wide false
943sideways false
944status open
945
946\begin_layout Plain Layout
947\align center
948\begin_inset Tabular
949<lyxtabular version="3" rows="2" columns="2">
950<features>
951<column alignment="center" valignment="top" width="0">
952<column alignment="center" valignment="top" width="0">
953<row>
954<cell alignment="center" valignment="top" topline="true" bottomline="true" leftline="true" usebox="none">
955\begin_inset Text
956
957\begin_layout Plain Layout
958\begin_inset Graphics
959        filename pmsm_spec.eps
960        scale 35
961
962\end_inset
963
964
965\end_layout
966
967\end_inset
968</cell>
969<cell alignment="center" valignment="top" topline="true" bottomline="true" leftline="true" rightline="true" usebox="none">
970\begin_inset Text
971
972\begin_layout Plain Layout
973\begin_inset Graphics
974        filename pmsm_simple.eps
975        scale 60
976
977\end_inset
978
979
980\end_layout
981
982\end_inset
983</cell>
984</row>
985<row>
986<cell alignment="center" valignment="top" topline="true" bottomline="true" leftline="true" usebox="none">
987\begin_inset Text
988
989\begin_layout Plain Layout
990
991\emph on
992Ilustrativní obrázek konstrukce PMSM
993\end_layout
994
995\end_inset
996</cell>
997<cell alignment="center" valignment="top" topline="true" bottomline="true" leftline="true" rightline="true" usebox="none">
998\begin_inset Text
999
1000\begin_layout Plain Layout
1001
1002\emph on
1003Zjednodušený model PMSM
1004\end_layout
1005
1006\end_inset
1007</cell>
1008</row>
1009</lyxtabular>
1010
1011\end_inset
1012
1013
1014\end_layout
1015
1016\begin_layout Plain Layout
1017\begin_inset Caption
1018
1019\begin_layout Plain Layout
1020Konstrukce a model PMSM
1021\end_layout
1022
1023\end_inset
1024
1025
1026\begin_inset CommandInset label
1027LatexCommand label
1028name "Flo:obr1_ilupmsm"
1029
1030\end_inset
1031
1032
1033\begin_inset CommandInset label
1034LatexCommand label
1035name "Flo:obr2_simplepmsm"
1036
1037\end_inset
1038
1039
1040\end_layout
1041
1042\begin_layout Plain Layout
1043
1044\end_layout
1045
1046\end_inset
1047
1048
1049\end_layout
1050
1051\begin_layout Standard
1052Základní konstrukce PMSM je na obrázku
1053\begin_inset CommandInset ref
1054LatexCommand ref
1055reference "Flo:obr1_ilupmsm"
1056
1057\end_inset
1058
1059.
1060 Nákres je pouze ilustrativní, ale zobrazuje hlavní části PMSM: Vnější kruh
1061 představuje stator.
1062 Na něm jsou zuby, na kterých je navinuto statorové vinutí (v obrázku není
1063 zobrazeno).
1064 Vnitřní kruh je rotor, na jehož povrchu jsou umístěny právě permanentní
1065 magnety.
1066 U těchto magnetů je barevně rozlišen severní a jižní pól.
1067 
1068\end_layout
1069
1070\begin_layout Standard
1071Často se lze setkat i s opačnou konstrukcí, kdy je stator umístěn uvnitř
1072 a rotor s magnety se otáčí kolem něj.
1073 Tato konstrukce PMSM se využívá například k pohonu nejrůznějších vozidel,
1074 kdy je motor umístěn přímo v kole vozidla, nebo k pohonu bubnu automatické
1075 pračky.
1076 Existují i další konstrukce PMSM.
1077 Zajímavou je například verze, která má otočný stator i rotor a toto zařízení
1078 pak může sloužit jako dělič výkonu.
1079\end_layout
1080
1081\begin_layout Standard
1082Vyobrazená konstrukce je někdy také označováná jako SMPMSM (
1083\emph on
1084Surface Mounted PMSM
1085\emph default
1086), tedy PMSM s magnety na povrchu.
1087 Další častou konstrukcí je IPMSM (
1088\emph on
1089Inner PMSM
1090\emph default
1091), kde jsou permanentní magnety umístěny uvnitř rotoru.
1092 Tyto verze mají nepatrně odlišné vlastnosti, které ale mají významný vliv
1093 při návrhu řízení těchto strojů.
1094 Pod PMSM se ještě zahrnují reluktanční motory, které jsou založeny na poněkud
1095 odlišném principu a dále se jimi vůbec zabývat nebudeme.
1096\end_layout
1097
1098\begin_layout Standard
1099Pro představu a odvození základních rovnic však nepotřebujeme pracovat s
1100 příliš složitou konstrukcí a vystačíme si se zjednodušeným modelem, který
1101 je zobrazen na obrázku
1102\begin_inset CommandInset ref
1103LatexCommand ref
1104reference "Flo:obr2_simplepmsm"
1105
1106\end_inset
1107
1108.
1109 Na statoru jsou zde umístěny pouze tři cívky, které představují vinutí
1110 jednotlivých fází.
1111 Rotor je pak reprezentován jediným permanentním magnetem.
1112 Pro základní představu je tento model dostačující, dále ale bude třeba
1113 rozšířit model o více párů pólů.
1114 PMSM na nákresu (zjednodušený model) má 1 pár pólů, ale reálné motory jich
1115 mívají obvykle více.
1116\end_layout
1117
1118\begin_layout Section
1119Souřadné soustavy
1120\end_layout
1121
1122\begin_layout Standard
1123Pro popis a následné odvození rovnic se standartně používá několik souřadných
1124 systémů.
1125 
1126\end_layout
1127
1128\begin_layout Standard
1129\begin_inset Float figure
1130wide false
1131sideways false
1132status open
1133
1134\begin_layout Plain Layout
1135\align center
1136\begin_inset Tabular
1137<lyxtabular version="3" rows="2" columns="3">
1138<features>
1139<column alignment="center" valignment="top" width="0">
1140<column alignment="center" valignment="top" width="0">
1141<column alignment="center" valignment="top" width="0">
1142<row>
1143<cell alignment="center" valignment="top" topline="true" bottomline="true" leftline="true" usebox="none">
1144\begin_inset Text
1145
1146\begin_layout Plain Layout
1147\begin_inset Graphics
1148        filename pmsm_simple_abc.eps
1149        scale 35
1150
1151\end_inset
1152
1153
1154\end_layout
1155
1156\end_inset
1157</cell>
1158<cell alignment="center" valignment="top" topline="true" bottomline="true" leftline="true" usebox="none">
1159\begin_inset Text
1160
1161\begin_layout Plain Layout
1162\begin_inset Graphics
1163        filename pmsm_simple_albe.eps
1164        scale 35
1165
1166\end_inset
1167
1168
1169\end_layout
1170
1171\end_inset
1172</cell>
1173<cell alignment="center" valignment="top" topline="true" bottomline="true" leftline="true" rightline="true" usebox="none">
1174\begin_inset Text
1175
1176\begin_layout Plain Layout
1177\begin_inset Graphics
1178        filename pmsm_simple_dq.eps
1179        scale 35
1180
1181\end_inset
1182
1183
1184\end_layout
1185
1186\end_inset
1187</cell>
1188</row>
1189<row>
1190<cell alignment="center" valignment="top" topline="true" bottomline="true" leftline="true" usebox="none">
1191\begin_inset Text
1192
1193\begin_layout Plain Layout
1194
1195\emph on
1196Souřadný systém
1197\begin_inset Formula $a-b-c$
1198\end_inset
1199
1200
1201\end_layout
1202
1203\end_inset
1204</cell>
1205<cell alignment="center" valignment="top" topline="true" bottomline="true" leftline="true" usebox="none">
1206\begin_inset Text
1207
1208\begin_layout Plain Layout
1209
1210\emph on
1211Souřadný systém
1212\begin_inset Formula $\alpha$
1213\end_inset
1214
1215-
1216\begin_inset Formula $\beta$
1217\end_inset
1218
1219
1220\end_layout
1221
1222\end_inset
1223</cell>
1224<cell alignment="center" valignment="top" topline="true" bottomline="true" leftline="true" rightline="true" usebox="none">
1225\begin_inset Text
1226
1227\begin_layout Plain Layout
1228
1229\emph on
1230Souřadný systém
1231\begin_inset Formula $d-q$
1232\end_inset
1233
1234
1235\end_layout
1236
1237\end_inset
1238</cell>
1239</row>
1240</lyxtabular>
1241
1242\end_inset
1243
1244
1245\end_layout
1246
1247\begin_layout Plain Layout
1248\begin_inset Caption
1249
1250\begin_layout Plain Layout
1251Souřadné systémy
1252\end_layout
1253
1254\end_inset
1255
1256
1257\end_layout
1258
1259\begin_layout Plain Layout
1260\begin_inset CommandInset label
1261LatexCommand label
1262name "Flo:obr3_ssabc"
1263
1264\end_inset
1265
1266
1267\begin_inset CommandInset label
1268LatexCommand label
1269name "Flo:obr4_ssalbe"
1270
1271\end_inset
1272
1273
1274\begin_inset CommandInset label
1275LatexCommand label
1276name "Flo:obr5_ssdq"
1277
1278\end_inset
1279
1280
1281\end_layout
1282
1283\end_inset
1284
1285Prvním z nich je souřadný systém
1286\emph on
1287
1288\begin_inset Formula $a-b-c$
1289\end_inset
1290
1291
1292\emph default
1293 znázorněný na obrázku
1294\begin_inset CommandInset ref
1295LatexCommand ref
1296reference "Flo:obr3_ssabc"
1297
1298\end_inset
1299
1300.
1301 Jednotlivé osy tohoto souřadného systému (
1302\begin_inset Formula $a$
1303\end_inset
1304
1305,
1306\begin_inset Formula $b$
1307\end_inset
1308
1309,
1310\begin_inset Formula $c$
1311\end_inset
1312
1313) jsou směřují ve směru os vinutí jednotlivých fází a jsou tedy vzájemně
1314 pootočeny o
1315\begin_inset Formula $120^{\circ}$
1316\end_inset
1317
1318.
1319 
1320\end_layout
1321
1322\begin_layout Standard
1323Protože ale k popsaní polohy v rovině jsou tři souřadnice (v osách
1324\begin_inset Formula $a$
1325\end_inset
1326
1327,
1328\begin_inset Formula $b$
1329\end_inset
1330
1331,
1332\begin_inset Formula $c$
1333\end_inset
1334
1335) zbytečné a jedna z nich je vždy závislá, přecházíme k souřadnému systému
1336 
1337\begin_inset Formula $\alpha$
1338\end_inset
1339
1340-
1341\begin_inset Formula $\beta$
1342\end_inset
1343
1344, který je znázorněn na obrázku
1345\begin_inset CommandInset ref
1346LatexCommand ref
1347reference "Flo:obr4_ssalbe"
1348
1349\end_inset
1350
1351.
1352 Osa
1353\begin_inset Formula $\alpha$
1354\end_inset
1355
1356 se totožná s osou
1357\emph on
1358
1359\begin_inset Formula $a$
1360\end_inset
1361
1362
1363\emph default
1364 ze souřadného systému
1365\emph on
1366
1367\begin_inset Formula $a-b-c$
1368\end_inset
1369
1370
1371\emph default
1372, osa
1373\begin_inset Formula $\beta$
1374\end_inset
1375
1376 ja na ní pak kolmá.
1377 Osy
1378\begin_inset Formula $\alpha$
1379\end_inset
1380
1381-
1382\begin_inset Formula $\beta$
1383\end_inset
1384
1385 tedy tvoří ortogonální systém.
1386\end_layout
1387
1388\begin_layout Standard
1389Pro většinu aplikací se však ukazuje výhodným přejít do rotující soustavy
1390 
1391\emph on
1392
1393\begin_inset Formula $d-q$
1394\end_inset
1395
1396
1397\emph default
1398, která je svázána s rotorem.
1399 Její vyobrazení je na obrázku
1400\begin_inset CommandInset ref
1401LatexCommand ref
1402reference "Flo:obr5_ssdq"
1403
1404\end_inset
1405
1406.
1407 Opět se jedná o ortogonální systém, kdy osu
1408\emph on
1409
1410\begin_inset Formula $d$
1411\end_inset
1412
1413
1414\emph default
1415 orientujeme ve směru osy permanentního magnetu směřující k jeho severnímu
1416 pólu.
1417 Osa
1418\emph on
1419
1420\begin_inset Formula $q$
1421\end_inset
1422
1423
1424\emph default
1425 je pak na ní kolmá.
1426\end_layout
1427
1428\begin_layout Section
1429Transformace souřadnic
1430\end_layout
1431
1432\begin_layout Standard
1433Mezi výše zmíněnými souřadnými soustavami platí následující převodní vztahy.
1434\end_layout
1435
1436\begin_layout Subsection
1437Transformace
1438\begin_inset Formula $a-b-c\longleftrightarrow\alpha-\beta$
1439\end_inset
1440
1441
1442\end_layout
1443
1444\begin_layout Standard
1445Tato transformace se označuje také jako Clarkova transformace, rovnice lze
1446 nalézt například v
1447\begin_inset CommandInset citation
1448LatexCommand cite
1449key "fiser2006"
1450
1451\end_inset
1452
1453, nebo je možné je poměrně snadno odvodit.
1454\end_layout
1455
1456\begin_layout Subsubsection
1457Převod
1458\begin_inset Formula $a-b-c\rightarrow\alpha-\beta$
1459\end_inset
1460
1461
1462\end_layout
1463
1464\begin_layout Standard
1465Osa
1466\begin_inset Formula $\alpha$
1467\end_inset
1468
1469 je totožná s osou
1470\begin_inset Formula $a$
1471\end_inset
1472
1473 osy
1474\begin_inset Formula $b$
1475\end_inset
1476
1477 a
1478\begin_inset Formula $c$
1479\end_inset
1480
1481 jsou pak oproti ní otočeny o
1482\begin_inset Formula $120^{\circ}$
1483\end_inset
1484
1485 respektive
1486\begin_inset Formula $-120^{\circ}$
1487\end_inset
1488
1489.
1490 Tedy souřadnice v ose
1491\begin_inset Formula $\alpha$
1492\end_inset
1493
1494 získáme následujícím průmětem z os
1495\begin_inset Formula $a,\: b,\: c$
1496\end_inset
1497
1498:
1499\begin_inset Formula \[
1500\alpha=k\left(a+b\cdot\cos(120^{\circ})+c\cdot\cos(-120^{\circ})\right)=k\left(a-\frac{1}{2}b-\frac{1}{2}c\right),\]
1501
1502\end_inset
1503
1504kde
1505\begin_inset Formula $k$
1506\end_inset
1507
1508 značí konstantu
1509\begin_inset Formula $k=\frac{2}{3}$
1510\end_inset
1511
1512.
1513 Obdobně postupujeme v případě osy
1514\begin_inset Formula $\beta$
1515\end_inset
1516
1517.
1518 Osa
1519\begin_inset Formula $a$
1520\end_inset
1521
1522 je na ní kolmá a tedy její příspěvek je nulový.
1523 Osy
1524\begin_inset Formula $b$
1525\end_inset
1526
1527 a
1528\begin_inset Formula $c$
1529\end_inset
1530
1531 promítnutne do osy
1532\begin_inset Formula $\beta$
1533\end_inset
1534
1535 získáme vztah:
1536\begin_inset Formula \[
1537\beta=k\left(b\cdot\sin(120^{\circ})+c\cdot\sin(-120^{\circ})\right)=k\left(\frac{\sqrt{3}}{2}b-\frac{\sqrt{3}}{2}c\right).\]
1538
1539\end_inset
1540
1541Celkem tedy máme rovnice:
1542\begin_inset Formula \begin{eqnarray*}
1543\alpha & = & \frac{2}{3}\left(a-\frac{1}{2}b-\frac{1}{2}c\right),\\
1544\beta & = & \frac{\sqrt{3}}{3}\left(b-c\right).\end{eqnarray*}
1545
1546\end_inset
1547
1548
1549\end_layout
1550
1551\begin_layout Subsubsection
1552Převod
1553\begin_inset Formula $\alpha-\beta\rightarrow a-b-c$
1554\end_inset
1555
1556
1557\end_layout
1558
1559\begin_layout Standard
1560Pro inverzní transformaci platí následující vztahy:
1561\begin_inset Formula \begin{eqnarray*}
1562a & = & \alpha+\theta,\\
1563b & = & \left(-\frac{1}{2}\alpha+\frac{\sqrt{3}}{2}\beta\right)+\theta,\\
1564c & \text{=} & \left(-\frac{1}{2}\alpha-\frac{\sqrt{3}}{2}\beta\right)+\theta,\end{eqnarray*}
1565
1566\end_inset
1567
1568kde
1569\begin_inset Formula $\theta$
1570\end_inset
1571
1572 představuje takzvanou nulovou složku
1573\begin_inset Formula $\theta=\frac{1}{3}\left(a+b+c\right)$
1574\end_inset
1575
1576.
1577\end_layout
1578
1579\begin_layout Subsection
1580Transformace
1581\begin_inset Formula $\alpha-\beta\longleftrightarrow d-q$
1582\end_inset
1583
1584
1585\begin_inset CommandInset label
1586LatexCommand label
1587name "sub:Transformace_albe_dq"
1588
1589\end_inset
1590
1591
1592\end_layout
1593
1594\begin_layout Standard
1595Transformace je označována jako Parkova transformace a představuje přechod
1596 do rotujícího souřadného systému.
1597 Rovnice transformace lze najít opět například v
1598\begin_inset CommandInset citation
1599LatexCommand cite
1600key "fiser2006"
1601
1602\end_inset
1603
1604 nebo je možné je opět odvodit.
1605\end_layout
1606
1607\begin_layout Subsubsection
1608Převod
1609\begin_inset Formula $\alpha-\beta\rightarrow d-q$
1610\end_inset
1611
1612
1613\end_layout
1614
1615\begin_layout Standard
1616Předpokládáme otočení doustavy
1617\begin_inset Formula $d-q$
1618\end_inset
1619
1620 oproti
1621\begin_inset Formula $\alpha-\beta$
1622\end_inset
1623
1624 o úhel
1625\begin_inset Formula $\phi$
1626\end_inset
1627
1628 kolem společného počátku souřadných soustav a tedy:
1629\begin_inset Formula \begin{eqnarray*}
1630d & = & \alpha\cos\phi+\beta\sin\phi,\\
1631q & = & -\alpha\sin\phi+\beta\cos\phi.\end{eqnarray*}
1632
1633\end_inset
1634
1635
1636\end_layout
1637
1638\begin_layout Subsubsection
1639Převod
1640\begin_inset Formula $d-q\rightarrow\alpha-\beta$
1641\end_inset
1642
1643
1644\end_layout
1645
1646\begin_layout Standard
1647Inverzní transformaci provedeme pouze otočením na druhou stranu:
1648\begin_inset Formula \begin{eqnarray*}
1649\alpha & = & d\cos\phi-q\sin\phi,\\
1650\beta & = & d\sin\phi+q\cos\phi.\end{eqnarray*}
1651
1652\end_inset
1653
1654
1655\end_layout
1656
1657\begin_layout Section
1658Odvození rovnic
1659\end_layout
1660
1661\begin_layout Subsection
1662Odvození rovnic v
1663\begin_inset Formula $d-q$
1664\end_inset
1665
1666 soustavě
1667\begin_inset CommandInset label
1668LatexCommand label
1669name "sub:Odvození-rovnic-vdq"
1670
1671\end_inset
1672
1673
1674\end_layout
1675
1676\begin_layout Standard
1677Rovnice v
1678\begin_inset Formula $d-q$
1679\end_inset
1680
1681 soustavě lze odvodit buď přímo nebo transformací rovnic z jiné soustavy.
1682 Přímé odvození bude uvedeno počínaje následujícím odstavcem, transformace
1683 z jiné soustavy (konkrétně
1684\begin_inset Formula $\alpha-\beta$
1685\end_inset
1686
1687) bude pro srovnání a kontrolu uvedeno dále v textu.
1688\end_layout
1689
1690\begin_layout Standard
1691Rovnici pro napětí v obvodu statoru synchroního stroje lze zapsat jako
1692\begin_inset Formula \[
1693u_{s}=R_{s}i_{s}+u_{i},\]
1694
1695\end_inset
1696
1697tedy součet napětí v obvodu (Ohmův zákon) a indukovaného napětí, přičemž
1698 veličiny jsou uvažovány komplexní.
1699 Vyjáříme-li indukované napětí, jako změnu toku v čase (Faradayův zákon
1700 elektromagnetické indukce) přejde rovnice na tvar
1701\begin_inset Formula \[
1702u_{s}=R_{s}i_{s}+\frac{d\psi_{s}}{dt}.\]
1703
1704\end_inset
1705
1706Pro přechod do rotujícího souřadného systému předpokládáme obecně rotaci
1707 o úhel
1708\begin_inset Formula $\varepsilon$
1709\end_inset
1710
1711, kterou provedeme vynásobením všech veličin operátorem rotace v komplexních
1712 číslech
1713\begin_inset Formula $e^{j\varepsilon}$
1714\end_inset
1715
1716, kde
1717\begin_inset Formula $j$
1718\end_inset
1719
1720 značí komplexní jednotku.
1721 Tedy
1722\begin_inset Formula \begin{eqnarray*}
1723u_{s}e^{j\varepsilon} & = & R_{s}i_{s}e^{j\varepsilon}+\frac{d(\psi_{s}e^{j\varepsilon})}{dt},\\
1724u_{s}e^{j\varepsilon} & \text{=} & R_{s}i_{s}e^{j\varepsilon}+\frac{d\psi_{s}}{dt}e^{j\varepsilon}+\psi_{s}j\omega_{\varepsilon}e^{j\varepsilon},\\
1725u_{s} & \text{=} & R_{s}i_{s}+\frac{d\psi_{s}}{dt}+\psi_{s}j\omega_{\varepsilon},\end{eqnarray*}
1726
1727\end_inset
1728
1729kde symbol
1730\begin_inset Formula $\omega_{\varepsilon}$
1731\end_inset
1732
1733 označuje úhlovou rychlost -- změnu úhlu
1734\begin_inset Formula $\varepsilon$
1735\end_inset
1736
1737, jedná se tedy o derivaci
1738\begin_inset Formula $\omega_{\varepsilon}=\frac{d\varepsilon}{dt}$
1739\end_inset
1740
1741.
1742 Tato úhlová rychlost
1743\begin_inset Formula $\omega_{\varepsilon}$
1744\end_inset
1745
1746 odpovídá elektrickým otáčkám
1747\begin_inset Formula $\omega_{el}$
1748\end_inset
1749
1750 a lze ji přepočíst na mechanické otáčky pomocí vztahu
1751\begin_inset Formula $\omega_{el}=p_{p}\omega_{m},$
1752\end_inset
1753
1754 kde
1755\begin_inset Formula $p_{p}$
1756\end_inset
1757
1758 je počet párů polů rotoru a
1759\begin_inset Formula $\omega_{m}$
1760\end_inset
1761
1762 mechanické otáčky.
1763 Když předpokládáme počet párů polů roven 1, je
1764\begin_inset Formula $\omega_{e}=\omega_{m}$
1765\end_inset
1766
1767.
1768\end_layout
1769
1770\begin_layout Standard
1771Nyní můžeme přejít k rovnicím v souřadném systému
1772\begin_inset Formula $d-q$
1773\end_inset
1774
1775, který je natočen oproti souřadnému systému statoru (
1776\begin_inset Formula $\alpha-\beta$
1777\end_inset
1778
1779) o úhel
1780\begin_inset Formula $\varepsilon=\vartheta$
1781\end_inset
1782
1783 a otáčí se rychlostí
1784\begin_inset Formula $\omega$
1785\end_inset
1786
1787.
1788 Osa magnetického toku rotoru je osou
1789\begin_inset Formula $d$
1790\end_inset
1791
1792 a v tomto směru uvažujeme reálnou složku komplexních veličin, osa
1793\begin_inset Formula $q$
1794\end_inset
1795
1796 je pak na ní kolmá a bude reprezentovat složku imaginární.
1797 Dostáváme tedy
1798\begin_inset Formula \[
1799u_{d}+ju_{q}\text{=}R_{s}\left(i_{d}+ji_{q}\right)+\frac{d\left(\psi_{d}+j\psi_{q}\right)}{dt}+\left(\psi_{d}+j\psi_{q}\right)j\omega_{m},\]
1800
1801\end_inset
1802
1803což při rozepsání po složkách (reálná a imaginární) vede na rovnice
1804\end_layout
1805
1806\begin_layout Standard
1807\begin_inset Formula \begin{eqnarray*}
1808u_{d} & = & R_{s}i_{d}+\frac{d\psi_{d}}{dt}-\omega_{m}\psi_{q},\\
1809u_{q} & = & R_{s}i_{q}+\frac{d\psi_{q}}{dt}+\omega_{m}\psi_{d}.\end{eqnarray*}
1810
1811\end_inset
1812
1813Dále uvažujme vztahy pro magnetické toky
1814\begin_inset Formula \begin{eqnarray*}
1815\psi_{d} & = & L_{d}i_{d}+\psi_{pm},\\
1816\psi_{q} & = & L_{q}i_{q}.\end{eqnarray*}
1817
1818\end_inset
1819
1820Po dosazení získáme rovnice
1821\begin_inset Formula \begin{eqnarray*}
1822u_{d} & = & R_{s}i_{d}+L_{d}\frac{di_{d}}{dt}-\omega_{m}L_{q}i_{q},\\
1823u_{q} & = & R_{s}i_{q}+L_{q}\frac{di_{q}}{dt}+\omega_{m}L_{d}i_{d}+\omega_{m}\psi_{pm}.\end{eqnarray*}
1824
1825\end_inset
1826
1827Vydělením
1828\begin_inset Formula $L_{d}$
1829\end_inset
1830
1831 respektive
1832\begin_inset Formula $L_{q}$
1833\end_inset
1834
1835 získáme
1836\begin_inset CommandInset label
1837LatexCommand label
1838name "rovnice_i_dq_ruzneL"
1839
1840\end_inset
1841
1842
1843\begin_inset Formula \begin{eqnarray*}
1844\frac{di_{d}}{dt} & = & -\frac{R_{s}}{L_{d}}i_{d}+\frac{L_{q}}{L_{d}}\omega_{m}i_{q}+\frac{1}{L_{d}}u_{d},\\
1845\frac{di_{q}}{dt} & = & -\frac{R_{s}}{L_{q}}-\frac{\psi_{pm}}{L_{q}}\omega_{m}-\frac{L_{d}}{L_{q}}\omega_{m}i_{d}+\frac{1}{L_{q}}u_{q}.\end{eqnarray*}
1846
1847\end_inset
1848
1849Když ale položíme
1850\begin_inset Formula $L_{d}=L_{q}=L_{s}$
1851\end_inset
1852
1853 dostaneme rovnice
1854\begin_inset Formula \begin{eqnarray*}
1855u_{d} & = & R_{s}i_{d}+L_{s}\frac{di_{d}}{dt}-\omega_{m}L_{s}i_{q},\\
1856u_{q} & = & R_{s}i_{q}+L_{s}\frac{di_{q}}{dt}+\omega_{m}L_{s}i_{d}+\omega_{m}\psi_{pm}.\end{eqnarray*}
1857
1858\end_inset
1859
1860Vydělení
1861\begin_inset Formula $L_{s}$
1862\end_inset
1863
1864 pak vede na tvar
1865\begin_inset Formula \begin{eqnarray*}
1866\frac{di_{d}}{dt} & = & -\frac{R_{s}}{L_{s}}i_{d}+\omega_{m}i_{q}+\frac{u_{d}}{L_{s}},\\
1867\frac{di_{q}}{dt} & = & -\frac{R_{s}}{L_{s}}i_{q}-\frac{\psi_{pm}}{L_{s}}\omega_{m}-\omega_{m}i_{d}+\frac{u_{q}}{L_{s}}.\end{eqnarray*}
1868
1869\end_inset
1870
1871Toto vyjádření je shodné s tím, které dostaneme následně transformací z
1872 
1873\begin_inset Formula $\alpha-\beta$
1874\end_inset
1875
1876 souřadné soustavy.
1877\end_layout
1878
1879\begin_layout Subsection
1880Odvození rovnic v
1881\begin_inset Formula $\alpha-\beta$
1882\end_inset
1883
1884 soustavě
1885\end_layout
1886
1887\begin_layout Standard
1888I když se pro řízení ukazuje být lepší a v praxi více využíváné vyjádření
1889 v soustave
1890\begin_inset Formula $d-q$
1891\end_inset
1892
1893, rovnice v
1894\begin_inset Formula $\alpha-\beta$
1895\end_inset
1896
1897 jsou také důležité, protože představují přímý vztah mezi měřenými a řízenými
1898 veličinami.
1899 Mohou být využity například při návrhu rozšířeného Kalmanova filtru.
1900\end_layout
1901
1902\begin_layout Standard
1903Opět vyjdeme z rovnice
1904\begin_inset Formula \[
1905u_{s}=R_{s}i_{s}+\frac{d\psi_{s}}{dt}.\]
1906
1907\end_inset
1908
1909Magnetický tok
1910\begin_inset Formula $\psi_{s}$
1911\end_inset
1912
1913 vyjádříme jako tok vytvořený cívkami statoru a dále přičteme tok permanentních
1914 magnetů, je však třeba uvažovat, že rotor obsahující permanentní magnety
1915 je natočen obecně pod úhlem
1916\begin_inset Formula $\vartheta$
1917\end_inset
1918
1919.
1920 Tedy v komplexní rovině lze vyjádřit tok jako
1921\begin_inset Formula \[
1922\psi_{s}=L_{s}i_{s}+\psi_{pm}e^{j\vartheta}.\]
1923
1924\end_inset
1925
1926Dosadíme nyní do rovnice a rozepíšeme ji po složkách
1927\begin_inset Formula \begin{eqnarray*}
1928u_{s} & = & R_{s}i_{s}+\frac{d\left(L_{s}i_{s}+\psi_{pm}e^{j\vartheta}\right)}{dt},\\
1929u_{\alpha}+ju_{\beta} & \text{=} & R_{s}\left(i_{\alpha}+ji_{\beta}\right)+\frac{d}{dt}\left(L_{s}\left(i_{\alpha}+ji_{\beta}\right)+\psi_{pm}\left(\cos\vartheta+j\sin\vartheta\right)\right).\end{eqnarray*}
1930
1931\end_inset
1932
1933Rozepsaní na dvě rovnice je pak následující
1934\begin_inset Formula \begin{eqnarray*}
1935u_{\alpha} & \text{=} & R_{s}i_{\alpha}+L_{s}\frac{di_{\alpha}}{dt}-\frac{d\vartheta}{dt}\psi_{pm}\sin\vartheta,\\
1936u_{\beta} & = & R_{s}i_{\beta}+L_{s}\frac{di_{\beta}}{dt}+\frac{d\vartheta}{dt}\psi_{pm}\cos\vartheta.\end{eqnarray*}
1937
1938\end_inset
1939
1940Rovnice vydělíme indukčností
1941\begin_inset Formula $L_{s}$
1942\end_inset
1943
1944, vyjádříme z nich derivace proudů a derivace úhlu natočení
1945\family roman
1946\series medium
1947\shape up
1948\size normal
1949\emph off
1950\bar no
1951\noun off
1952\color none
1953označíme jako úhlovou rychlost
1954\begin_inset Formula $\frac{d\vartheta}{dt}$
1955\end_inset
1956
1957=
1958\begin_inset Formula $\omega$
1959\end_inset
1960
1961.
1962 Následně dostaneme rovnice v souřadné soustavě
1963\family default
1964\series default
1965\shape default
1966\size default
1967\emph default
1968\bar default
1969\noun default
1970\color inherit
1971 
1972\begin_inset Formula $\alpha-\beta$
1973\end_inset
1974
1975:
1976\begin_inset Formula \begin{eqnarray*}
1977\frac{di_{\alpha}}{dt} & = & -\frac{R_{s}}{L_{s}}i_{\alpha}+\frac{\psi_{pm}}{L_{s}}\omega\sin\vartheta+\frac{u_{\alpha}}{L_{s}},\\
1978\frac{di_{\beta}}{dt} & = & -\frac{R_{s}}{L_{s}}i_{\beta}-\frac{\psi_{pm}}{L_{s}}\omega\cos\vartheta+\frac{u_{\beta}}{L_{s}}.\end{eqnarray*}
1979
1980\end_inset
1981
1982
1983\end_layout
1984
1985\begin_layout Standard
1986Nyní je ještě třeba přidat další dvě diferenciální rovnice pro otáčky
1987\begin_inset Formula $\omega$
1988\end_inset
1989
1990 a polohu
1991\begin_inset Formula $\vartheta$
1992\end_inset
1993
1994.
1995 Rovnice pro
1996\begin_inset Formula $\vartheta$
1997\end_inset
1998
1999 je triviální a už byla užita, jedná se o
2000\begin_inset Formula \[
2001\frac{d\vartheta}{dt}=\omega.\]
2002
2003\end_inset
2004
2005
2006\end_layout
2007
2008\begin_layout Subsubsection
2009Rovnice pro
2010\begin_inset Formula $\omega$
2011\end_inset
2012
2013 
2014\end_layout
2015
2016\begin_layout Standard
2017Rovnice pro
2018\begin_inset Formula $\omega$
2019\end_inset
2020
2021 získáme následujícím postupem ze základních zákonů mechaniky: Pro točivý
2022 moment (speciální případ momentu síly pro silovou dvojici, kdy se vektory
2023 skládají na nulu, avšak mají točivý účinek, v anglické literatuře označeno
2024 jako
2025\emph on
2026torque
2027\emph default
2028)
2029\emph on
2030 
2031\emph default
2032platí obecně vztah
2033\begin_inset Formula \[
2034\tau=\frac{dL}{dt},\]
2035
2036\end_inset
2037
2038kde
2039\begin_inset Formula $L$
2040\end_inset
2041
2042 označuje moment hybnosti (
2043\emph on
2044angular momentum
2045\emph default
2046).
2047 Při uvažování působení více točivých momentu momentů pak
2048\begin_inset Formula \[
2049\tau_{1}+\ldots+\tau_{n}=\sum\tau=\frac{dL}{dt}.\]
2050
2051\end_inset
2052
2053Uvažujeme-li rotaci kolem pevné osy, lze moment hybnosti vyjádřit jako
2054\begin_inset Formula \[
2055L=J\omega_{m},\]
2056
2057\end_inset
2058
2059kde
2060\begin_inset Formula $J$
2061\end_inset
2062
2063 označuje moment setrvačnosti (
2064\emph on
2065moment of inertia
2066\emph default
2067) a
2068\begin_inset Formula $\omega_{m}$
2069\end_inset
2070
2071 je mechanická úhlová rychlost.
2072 Po dosazení tedy
2073\begin_inset Formula \[
2074\sum\tau=\frac{dL}{dt}=\frac{d(J\omega_{m})}{dt}=J\frac{d\omega_{m}}{dt}.\]
2075
2076\end_inset
2077
2078Točivé momenty
2079\begin_inset Formula $\sum\tau$
2080\end_inset
2081
2082 jsou:
2083\end_layout
2084
2085\begin_layout Itemize
2086moment získaný konverzním procesem elektrické energie, který vyjadřuje hlavní
2087 vlastnost točivého stroje, a to právě převod elektrické energie na mechanickou,
2088 tento mement označíme jako
2089\begin_inset Formula $T_{e}$
2090\end_inset
2091
2092
2093\end_layout
2094
2095\begin_layout Itemize
2096zátěžný moment reprezentující zatížení stroje, tedy v podstatě to, co je
2097 motorem poháněno, je však třeba uvažovat, že působí v opačném směru a stroj
2098 brzdí, označíme ho tedy
2099\begin_inset Formula $-T_{L}$
2100\end_inset
2101
2102
2103\end_layout
2104
2105\begin_layout Itemize
2106dále je ještě třeba uvažovat ztráty ve stroji v důsledku tření, tento moment
2107 opět působí v opačném směru a uvažujeme jej lineárně závislý na otáčkách
2108 
2109\begin_inset Formula $\omega_{m}$
2110\end_inset
2111
2112, tedy
2113\begin_inset Formula $-B\omega_{m}$
2114\end_inset
2115
2116, kde
2117\begin_inset Formula $B$
2118\end_inset
2119
2120 je koeficient viskozity (tření)
2121\end_layout
2122
2123\begin_layout Standard
2124Rovnice po dosazení tedy přejde na tvar
2125\begin_inset Formula \[
2126T_{e}-T_{L}-B\omega_{m}=J\frac{d\omega_{m}}{dt}.\]
2127
2128\end_inset
2129
2130Nyní je ještě třeba vyjádřit točívý moment
2131\begin_inset Formula $T_{e}$
2132\end_inset
2133
2134 na základě elektrických veličin.
2135 Toho lze dosáhnout výpočtem přes okamžitý elektrický výkon, pro trojfázový
2136 systém
2137\begin_inset Formula \[
2138P=u_{a}i_{a}+u_{b}i_{b}+u_{c}i_{c}.\]
2139
2140\end_inset
2141
2142Po transformaci do systému
2143\begin_inset Formula $\alpha-\beta$
2144\end_inset
2145
2146 získáme vyjádření
2147\begin_inset Formula \[
2148P=k_{p}\left(u_{\alpha}i_{\alpha}+u_{\beta}i_{\beta}\right),\]
2149
2150\end_inset
2151
2152kde
2153\begin_inset Formula $k_{p}$
2154\end_inset
2155
2156 označuje Parkovu konstantu s hodnotou
2157\begin_inset Formula $k_{p}=\frac{3}{2}$
2158\end_inset
2159
2160.
2161 Napětí je zde uvažováno indukované
2162\begin_inset Formula $u_{i}=\frac{d\psi_{s}}{dt}=\frac{d\left(L_{s}i_{s}+\psi_{pm}e^{j\vartheta}\right)}{dt}=L_{s}\frac{di_{s}}{dt}+j\omega\psi_{pm}e^{j\vartheta}$
2163\end_inset
2164
2165 a z něj využijeme pouze složku bez derivace proudu, protože ta slouží k
2166 tvorbě samotného magnetického pole stroje a nepodílí se na tvorbě výkonu,
2167 tedy
2168\begin_inset Formula $\omega\psi_{pm}j(\cos\vartheta+j\sin\vartheta)$
2169\end_inset
2170
2171.
2172 V systému
2173\begin_inset Formula $\alpha\beta$
2174\end_inset
2175
2176 získáme vyjádření
2177\begin_inset Formula \begin{eqnarray*}
2178u_{\alpha} & = & -\omega\psi_{pm}\sin\vartheta,\\
2179u_{\beta} & = & \omega\psi_{pm}\cos\vartheta,\end{eqnarray*}
2180
2181\end_inset
2182
2183tedy po dosazení
2184\begin_inset Formula \[
2185P=k_{p}\left(-i_{\alpha}\omega\psi_{pm}\sin\vartheta+i_{\beta}\omega\psi_{pm}\cos\vartheta\right).\]
2186
2187\end_inset
2188
2189Moment
2190\begin_inset Formula $T_{e}$
2191\end_inset
2192
2193 lze pak určit ze vztahu
2194\begin_inset Formula $P=\omega_{m}T_{e}$
2195\end_inset
2196
2197 a tedy
2198\begin_inset Formula \[
2199T_{e}=\frac{P}{\omega_{m}}=k_{p}\frac{i_{\beta}\omega\psi_{pm}\cos\vartheta-i_{\alpha}\omega\psi_{pm}\sin\vartheta}{\omega_{m}}=k_{p}p_{p}\psi_{pm}\left(i_{\beta}\cos\vartheta-i_{\alpha}\sin\vartheta\right),\]
2200
2201\end_inset
2202
2203kde jsme využili vztahu
2204\begin_inset Formula $\frac{\omega}{\omega_{m}}=p_{p}$
2205\end_inset
2206
2207.
2208\end_layout
2209
2210\begin_layout Standard
2211Dosazení do rovnice pro momenty pak vede na tvar
2212\begin_inset Formula \[
2213k_{p}p_{p}\psi_{pm}\left(i_{\beta}\cos\vartheta-i_{\alpha}\sin\vartheta\right)-T_{L}-B\omega_{m}=J\frac{d\omega_{m}}{dt}.\]
2214
2215\end_inset
2216
2217Ještě je třeba upravit rovnici tak, aby v ní nevystupovaly mechanické otáčky
2218 
2219\begin_inset Formula $\omega_{m}$
2220\end_inset
2221
2222, ale otáčky elektrické
2223\begin_inset Formula $\omega$
2224\end_inset
2225
2226.
2227 Toho je možno snadno dosáhnout násobením celé rovnice
2228\begin_inset Formula $p_{p}$
2229\end_inset
2230
2231.
2232 Rovnici ještě vydělíme momentem setrvačnosti
2233\begin_inset Formula $J$
2234\end_inset
2235
2236 a získáme tvar
2237\begin_inset Formula \[
2238\frac{d\omega}{dt}=\frac{k_{p}p_{p}^{2}\psi_{pm}}{J}\left(i_{\beta}\cos\vartheta-i_{\alpha}\sin\vartheta\right)-\frac{T_{L}p_{p}}{J}-\frac{B}{J}\omega.\]
2239
2240\end_inset
2241
2242Tedy máme poslední rovnici následující soustavy:
2243\end_layout
2244
2245\begin_layout Standard
2246\begin_inset Formula \begin{eqnarray*}
2247\frac{di_{\alpha}}{dt} & \text{=} & -\frac{R_{s}}{L_{s}}i_{\alpha}+\frac{\psi_{pm}}{L_{s}}\omega\sin\vartheta+\frac{u_{\alpha}}{L_{s}},\\
2248\frac{di_{\beta}}{dt} & \text{=} & -\frac{R_{s}}{L_{s}}i_{\beta}-\frac{\psi_{pm}}{L_{s}}\omega\cos\vartheta+\frac{u_{\beta}}{L_{s}},\\
2249\frac{d\omega}{dt} & \text{=} & \frac{k_{p}p_{p}^{2}\psi_{pm}}{J}\left(i_{\beta}\cos\vartheta-i_{\alpha}\sin\vartheta\right)-\frac{B}{J}\omega-\frac{p_{p}}{J}T_{L},\\
2250\frac{d\vartheta}{dt} & \text{=} & \omega.\end{eqnarray*}
2251
2252\end_inset
2253
2254
2255\end_layout
2256
2257\begin_layout Subsection
2258Odvození rovnice pro
2259\begin_inset Formula $\omega$
2260\end_inset
2261
2262 v
2263\begin_inset Formula $d-q$
2264\end_inset
2265
2266 soustavě pro různé indukčnosti
2267\end_layout
2268
2269\begin_layout Standard
2270Zatím jsme ve většině případů předchozího odvození učinili zjednodušující
2271 předpoklad stejných indukčností
2272\begin_inset Formula $L_{d}=L_{q}=L_{s}$
2273\end_inset
2274
2275.
2276 To relativně dobře platí pro případ SMPMSM.
2277 Pro IPMSM a přesnější model SMPMSM toto však neplatí a
2278\begin_inset Formula $L_{d}\neq L_{q}$
2279\end_inset
2280
2281.
2282 Tato vlastnost bude také velmi důležitá při užití estimačních technik označovan
2283ých jako
2284\emph on
2285injektáže
2286\emph default
2287 (detailněji dále v textu).
2288 Mít tedy k dispozici i rovnice pro různé indukčnosti je velmi žádoucí.
2289 Rovnice pro proudy v
2290\begin_inset Formula $d-q$
2291\end_inset
2292
2293 souřadnicích s různými indukčnostmi jsou již uvedeny v části
2294\begin_inset CommandInset ref
2295LatexCommand ref
2296reference "rovnice_i_dq_ruzneL"
2297
2298\end_inset
2299
2300.
2301 Rovnice pro
2302\begin_inset Formula $\omega$
2303\end_inset
2304
2305 bude odvozena nyní:
2306\end_layout
2307
2308\begin_layout Standard
2309Opět vyjdeme z analogických vztahů jako při předchozím odvození pro
2310\begin_inset Formula $\alpha-\beta$
2311\end_inset
2312
2313, tedy
2314\begin_inset Formula \[
2315T_{e}-T_{L}-B\omega_{m}=J\frac{d\omega_{m}}{dt},\]
2316
2317\end_inset
2318
2319kde vyjádříme
2320\begin_inset Formula $T_{e}$
2321\end_inset
2322
2323 ze vztahu
2324\begin_inset Formula \[
2325T_{e}=\frac{P}{\omega_{m}}.\]
2326
2327\end_inset
2328
2329Tedy transformujeme následující vyjádření pro výkond z
2330\begin_inset Formula $\alpha-\beta$
2331\end_inset
2332
2333 do
2334\begin_inset Formula $d-q$
2335\end_inset
2336
2337 
2338\begin_inset Formula \begin{eqnarray*}
2339P & = & k_{p}\left(u_{\alpha}i_{\alpha}+u_{\beta}i_{\beta}\right),\\
2340P & = & k_{p}\left(\left(u_{d}\cos\vartheta-u_{q}\sin\vartheta\right)\left(i_{d}\cos\vartheta-i_{q}\sin\vartheta\right)+\left(u_{q}\cos\vartheta+u_{d}\sin\vartheta\right)\left(i_{q}\cos\vartheta+i_{d}\sin\vartheta\right)\right),\\
2341P & \text{=} & k_{p}\left(u_{d}i_{d}+u_{q}i_{q}\right).\end{eqnarray*}
2342
2343\end_inset
2344
2345Opět dosadíme za
2346\begin_inset Formula $u_{d,q}$
2347\end_inset
2348
2349 složky indukovaného napětí bez derivace proudů
2350\begin_inset Formula \begin{eqnarray*}
2351u_{d} & = & -\omega L_{q}i_{q},\\
2352u_{q} & = & \omega L_{d}i_{d}+\omega\psi_{pm}.\end{eqnarray*}
2353
2354\end_inset
2355
2356To vede na
2357\begin_inset Formula \begin{eqnarray*}
2358P & = & k_{p}\left(-\omega L_{q}i_{q}i_{d}+\left(\omega L_{d}i_{d}+\omega\psi_{pm}\right)i_{q}\right),\\
2359P & = & k_{p}\omega\left(i_{d}i_{q}\left(L_{d}-L_{q}\right)+\psi_{pm}i_{q}\right).\end{eqnarray*}
2360
2361\end_inset
2362
2363A po dosazení získáme vyjádření pro moment
2364\begin_inset Formula $T_{e}$
2365\end_inset
2366
2367 ve tvaru
2368\begin_inset Formula \[
2369T_{e}=k_{p}p_{p}\left(i_{d}i_{q}\left(L_{d}-L_{q}\right)+\psi_{pm}i_{q}\right).\]
2370
2371\end_inset
2372
2373Rovnice
2374\begin_inset Formula $T_{e}-T_{L}-B\omega_{m}=J\frac{d\omega_{m}}{dt}$
2375\end_inset
2376
2377 pak po dosazení
2378\begin_inset Formula $T_{e}$
2379\end_inset
2380
2381, vydělení
2382\begin_inset Formula $J$
2383\end_inset
2384
2385 a násobení
2386\begin_inset Formula $p_{p}$
2387\end_inset
2388
2389 přejde na tvar
2390\begin_inset Formula \[
2391\frac{d\omega}{dt}=\frac{k_{p}p_{p}^{2}}{J}\left(\left(L_{d}-L_{q}\right)i_{d}i_{q}+\psi_{pm}i_{q}\right)-\frac{B}{J}\omega-\frac{p_{p}}{J}T_{L}.\]
2392
2393\end_inset
2394
2395
2396\end_layout
2397
2398\begin_layout Subsection
2399Diskretizace
2400\end_layout
2401
2402\begin_layout Standard
2403Výpočty jsou prováděny výhradně na počítači, simulace na PC a v případě
2404 řízení reálného stroje se obvykle užívá DSP.
2405 Je tedy třeba výše odvozené diferenciální rovnice diskretizovat a převést
2406 na rovnice diferenční.
2407 Diskretizaci je vhodné volit co možná nejjednodušší, aby se příliš nekomplikova
2408ly výsledné rovnice a aby bylo umožněno jejich případné rychlé zpracování
2409 v reálném čase.
2410\end_layout
2411
2412\begin_layout Standard
2413Diskretizací pomocí Eulerovy metody s časovým krokem
2414\begin_inset Formula $\Delta t$
2415\end_inset
2416
2417 získáme následující diskrétní diferenční rovnice:
2418\end_layout
2419
2420\begin_layout Standard
2421\begin_inset Formula \begin{eqnarray*}
2422i_{\alpha,t+1} & \text{=} & \left(1-\frac{R_{s}}{L_{s}}\Delta t\right)i_{\alpha,t}+\frac{\psi_{pm}}{L_{s}}\omega_{t}\sin\vartheta_{t}+\frac{u_{\alpha,t}}{L_{s}}\\
2423i_{\beta,t+1} & \text{=} & \left(1-\frac{R_{s}}{L_{s}}\Delta t\right)i_{\beta,t}-\frac{\psi_{pm}}{L_{s}}\omega_{t}\cos\vartheta_{t}+\frac{u_{\beta,t}}{L_{s}}\\
2424\omega_{t+1} & \text{=} & \left(1-\frac{B}{J}\Delta t\right)\omega_{t}+\frac{k_{p}p_{p}^{2}\psi_{pm}}{J}\Delta t\left(i_{\beta,t}\cos\vartheta_{t}-i_{\alpha,t}\sin\vartheta_{t}\right)-\frac{p_{p}}{J}T_{L}\Delta t\\
2425\vartheta_{t+1} & \text{=} & \vartheta_{t}+\omega_{t}\Delta t\end{eqnarray*}
2426
2427\end_inset
2428
2429
2430\end_layout
2431
2432\begin_layout Subsection
2433Rotace do
2434\begin_inset Formula $d-q$
2435\end_inset
2436
2437
2438\end_layout
2439
2440\begin_layout Standard
2441Nyní ještě provedeme rotaci rovnic ze souřadnic
2442\begin_inset Formula $\alpha-\beta$
2443\end_inset
2444
2445 do
2446\begin_inset Formula $d-q$
2447\end_inset
2448
2449.
2450 Jednak v diferenciálním případě, který bude následovat diskretizace, ale
2451 také v diskrétním případě diferenčních rovnic.
2452 Oba postupy pak budou srovnány.
2453\end_layout
2454
2455\begin_layout Standard
2456Převod do rotující souřadné soustavy
2457\begin_inset Formula $d-q$
2458\end_inset
2459
2460 pootočené o úhel
2461\begin_inset Formula $\vartheta$
2462\end_inset
2463
2464 a rotojící rychlostí
2465\begin_inset Formula $\omega$
2466\end_inset
2467
2468:
2469\begin_inset Formula \[
2470\left[\begin{array}{c}
2471x_{d}\\
2472x_{q}\end{array}\right]\text{=}\left[\begin{array}{cc}
2473\cos\vartheta & \sin\vartheta\\
2474-\sin\vartheta & \cos\vartheta\end{array}\right]\left[\begin{array}{c}
2475x_{\alpha}\\
2476x_{\beta}\end{array}\right],\]
2477
2478\end_inset
2479
2480viz
2481\begin_inset CommandInset ref
2482LatexCommand ref
2483reference "sub:Transformace_albe_dq"
2484
2485\end_inset
2486
2487 nebo stejného efektu lze dosáhnout i použítím komplexních souřadnic a zápisem
2488 
2489\begin_inset Formula $x_{dq}=e^{j\vartheta}x_{\alpha\beta}$
2490\end_inset
2491
2492, jako v části
2493\begin_inset CommandInset ref
2494LatexCommand ref
2495reference "sub:Odvození-rovnic-vdq"
2496
2497\end_inset
2498
2499.
2500\end_layout
2501
2502\begin_layout Standard
2503Následně tedy
2504\begin_inset Formula \begin{alignat*}{2}
2505i_{d} & = & i_{\alpha}\cos\vartheta+i_{\beta}\sin\vartheta,\\
2506i_{q} & = & i_{\beta}\cos\vartheta-i_{\alpha}\sin\vartheta,\end{alignat*}
2507
2508\end_inset
2509
2510
2511\end_layout
2512
2513\begin_layout Standard
2514a analogicky pro
2515\begin_inset Formula $u$
2516\end_inset
2517
2518.
2519 Naopak pro inverzní transformaci
2520\end_layout
2521
2522\begin_layout Standard
2523\begin_inset Formula \begin{alignat*}{2}
2524i_{\alpha} & = & i_{d}\cos\vartheta-i_{q}\sin\vartheta,\\
2525i_{\beta} & = & i_{q}\cos\vartheta+i_{d}\sin\vartheta,\end{alignat*}
2526
2527\end_inset
2528
2529
2530\end_layout
2531
2532\begin_layout Standard
2533a opět anoalogicky pro
2534\begin_inset Formula $u$
2535\end_inset
2536
2537.
2538 To po dosazení do původních diferenciálních rovnic vede na
2539\end_layout
2540
2541\begin_layout Standard
2542\begin_inset Formula \begin{eqnarray*}
2543\frac{d(i_{d}\cos\vartheta-i_{q}\sin\vartheta)}{dt} & \text{=} & -\frac{R_{s}}{L_{s}}(i_{d}\cos\vartheta-i_{q}\sin\vartheta)+\frac{\psi_{pm}}{L_{s}}\omega\sin\vartheta+\frac{u_{d}\cos\vartheta-u_{q}\sin\vartheta}{L_{s}},\\
2544\frac{d(i_{q}\cos\vartheta+i_{d}\sin\vartheta)}{dt} & \text{=} & -\frac{R_{s}}{L_{s}}(i_{q}\cos\vartheta+i_{d}\sin\vartheta)-\frac{\psi_{pm}}{L_{s}}\omega\cos\vartheta+\frac{u_{q}\cos\vartheta+u_{d}\sin\vartheta}{L_{s}},\\
2545\frac{d\omega}{dt} & \text{=} & \frac{k_{p}p_{p}^{2}\psi_{pm}}{J}\left(i_{q}\right)-\frac{B}{J}\omega-\frac{p_{p}}{J}T_{L},\\
2546\frac{d\vartheta}{dt} & \text{=} & \omega.\end{eqnarray*}
2547
2548\end_inset
2549
2550
2551\end_layout
2552
2553\begin_layout Standard
2554Ve třetí rovnici rovnou dosadíme
2555\begin_inset Formula $i_{q}$
2556\end_inset
2557
2558, čtvrtá se nemění a z prvních dvou vyjádříme rovnice pro proudy a napětí
2559 v
2560\begin_inset Formula $d$
2561\end_inset
2562
2563 a
2564\begin_inset Formula $q$
2565\end_inset
2566
2567, například tak, že první rovnici násobíme
2568\begin_inset Formula $\cos\vartheta$
2569\end_inset
2570
2571 a sečteme s druhou násobenou
2572\begin_inset Formula $\sin\vartheta$
2573\end_inset
2574
2575, dále pak první rovnici násobenou
2576\begin_inset Formula $-\sin\vartheta$
2577\end_inset
2578
2579 sečteme s druhou násobenou
2580\begin_inset Formula $\cos\vartheta$
2581\end_inset
2582
2583, tento postup vede na rovnice
2584\end_layout
2585
2586\begin_layout Standard
2587\begin_inset Formula \begin{eqnarray*}
2588\frac{di_{d}}{dt} & \text{=} & -\frac{R_{s}}{L_{s}}i_{d}-i_{q}\omega+\frac{u_{d}}{L_{s}},\\
2589\frac{di_{q}}{dt} & \text{=} & i_{d}\omega-\frac{R_{s}}{L_{s}}i_{q}-\frac{\psi_{pm}}{L_{s}}\omega+\frac{u_{q}}{L_{s}},\\
2590\frac{d\omega}{dt} & \text{=} & \frac{k_{p}p_{p}^{2}\psi_{pm}}{J}i_{q}-\frac{B}{J}\omega-\frac{p_{p}}{J}T_{L},\\
2591\frac{d\vartheta}{dt} & \text{=} & \omega.\end{eqnarray*}
2592
2593\end_inset
2594
2595
2596\end_layout
2597
2598\begin_layout Standard
2599Zde jsou zajímavé členy
2600\begin_inset Formula $-i_{q}\omega$
2601\end_inset
2602
2603 a
2604\begin_inset Formula $i_{d}\omega$
2605\end_inset
2606
2607 v první a druhé rovnici, protože když bychom nejdříve provedli diskretizaci
2608 a až následně převod do
2609\begin_inset Formula $d-q$
2610\end_inset
2611
2612 souřadnic, tyto členy zřejmě nevzniknou.
2613 Nevzniknou také, když soustavu
2614\begin_inset Formula $d-q$
2615\end_inset
2616
2617 definujeme ne jako pootočenou o
2618\begin_inset Formula $\vartheta$
2619\end_inset
2620
2621, ale jako soustavu pootočenou o nějaké konstantní
2622\begin_inset Formula $\varepsilon$
2623\end_inset
2624
2625.
2626 Z formálního hlediska se jeví jako nejvíce správné řešení zahrnující tyto
2627 členy.
2628 Pro praktické použití ale je vhodné otestovat, jaký je vliv těchto členů.
2629 Diskretizovaná verze rovnic v
2630\begin_inset Formula $d-q$
2631\end_inset
2632
2633 je tedy
2634\end_layout
2635
2636\begin_layout Standard
2637\begin_inset Formula \begin{eqnarray*}
2638i_{d,t+1}+\left|\overline{\underline{\left(-i_{q,t}\omega_{t}\right)}}\right| & \text{=} & \left(1-\frac{R_{s}}{L_{s}}\Delta t\right)i_{d,t}+\frac{u_{d,t}}{L_{s}},\\
2639i_{q,t+1}+\left|\underline{\overline{\left(+i_{d,t}\omega_{t}\right)}}\right| & \text{=} & \left(1-\frac{R_{s}}{L_{s}}\Delta t\right)i_{q,t}-\frac{\psi_{pm}}{L_{s}}\omega_{t}+\frac{u_{q,t}}{L_{s}},\\
2640\omega_{t+1} & \text{=} & \left(1-\frac{B}{J}\Delta t\right)\omega_{t}+\Delta t\frac{k_{p}p_{p}^{2}\psi_{pm}}{J}i_{q,t}-\frac{p_{p}}{J}T_{L}\Delta t,\\
2641\vartheta_{t+1} & \text{=} & \vartheta_{t}+\omega_{t}\Delta t,\end{eqnarray*}
2642
2643\end_inset
2644
2645kde
2646\begin_inset Quotes gld
2647\end_inset
2648
2649problematické
2650\begin_inset Quotes grd
2651\end_inset
2652
2653 členy jsou v rámečku.
2654\end_layout
2655
2656\begin_layout Section
2657Problematika modelu
2658\end_layout
2659
2660\begin_layout Standard
2661Dále budeme pracovat zpravidla převážně s rovnicemi odvozenými v předchozí
2662 části a skutečný stroj ustoupí do pozadí.
2663 Je však třeba mít na paměti, že za rovnicemi se skrývá fyzikální realita
2664 a mnoho jevů, které ji doprovází.
2665 Tyto jevy se totiž při aplikaci regulátoru na skutečném stroji projeví.
2666 Jedná se především o následující body:
2667\end_layout
2668
2669\begin_layout Itemize
2670
2671\series bold
2672nepřesnost modelu
2673\series default
2674 -- chyby způsobené zanedbáním nejrůznějších fyzikálních vlivů a důsledky
2675 zjednodušujících předpokladů, například závislosti některých veličin na
2676 teplotě, sycení magnetických obvodů, obecně nekonstantní parametry stroje
2677 atd.
2678\end_layout
2679
2680\begin_layout Itemize
2681
2682\series bold
2683nedokonalosti stroje
2684\series default
2685 -- žádný stroj nebude vyrobený přesně, aby odpovídal modelu, vyskytnou
2686 se různé nerovnosti, nesymetrie a podobně
2687\end_layout
2688
2689\begin_layout Itemize
2690
2691\series bold
2692diskretizační a zaokrouhlovací chyby
2693\series default
2694 -- řízení je navrhováno pro digitální počítač a tedy dříve nebo později
2695 je třeba provést diskretizaci a kvantizaci všech zpracovávaných veličin
2696\end_layout
2697
2698\begin_layout Itemize
2699
2700\series bold
2701chyby měření
2702\series default
2703 -- měřící přístroje a čidla, která získávají informace o motoru nejsou
2704 přesná, mají pouze určitou rozlišovací schopnost a také omezenou možnost
2705 předat informaci, zejména pokud se jedná o digitální zařízení
2706\end_layout
2707
2708\begin_layout Itemize
2709
2710\series bold
2711napájecí zdroj
2712\series default
2713 -- zařízení, které dodává regulátorem požadované napětí do stroje není
2714 ideální, naopak odpovídá ideálním požadavkům zpravidla velmi špatně, využívá
2715 pulzní šířkové modulace (PWM) a invertoru; tyto zařízení pak přinášejí
2716 množství negativních efektů
2717\end_layout
2718
2719\begin_layout Standard
2720Tyto jevy se velmi těžko popisují a jejich zachycení v modelu přináší mnoho
2721 komplikací.
2722 Většinu z nich ani nedokážeme popsat a předvídat.
2723 Proto se pokusíme co nejvíce z výše zmíněných problémů zahrnout pod pojem
2724 šum.
2725 Vzniká pak ale otázka, jak takový šum vhodně nastavit v modelu, aby alespoň
2726 přibližně odpovídal problematickým jevům.
2727 V rovnicích z předchozí části tedy budeme navíc ještě uvažovat jednoduchý
2728 model šumu a to aditivní bílý Gaussovský šum.
2729\end_layout
2730
2731\begin_layout Section
2732Estimace stavových veličin
2733\end_layout
2734
2735\begin_layout Subsubsection
2736Mechanické veličiny
2737\end_layout
2738
2739\begin_layout Standard
2740Pro řízení PMSM je důležité, že se jedná o synchronní stroj, kdy se rotor
2741 otáčí současně (synchronně) s točivým magnetickým polem vytvořeným cívkami
2742 statoru.
2743 Proto, když chceme navrhnout řízení takového stroje musíme nutně znát polohu
2744 rotoru
2745\begin_inset Formula $\vartheta$
2746\end_inset
2747
2748, a to s relativně velkou přesností.
2749 Dále, protože se v textu zaměřujeme na řízení rychlosti stroje (regulovanou
2750 veličinou jsou otáčky rotoru) potřebujeme znát i hodnotu otáček
2751\begin_inset Formula $\omega$
2752\end_inset
2753
2754.
2755 Problematika získání těchto hodnot se však ukazuje být netriviální.
2756 Obecně existuje několik přístupů, které budou detailněji rozebrány dále
2757 v textu.
2758\end_layout
2759
2760\begin_layout Paragraph
2761Poznámka:
2762\end_layout
2763
2764\begin_layout Standard
2765Zmiňované veličiny
2766\begin_inset Formula $\vartheta$
2767\end_inset
2768
2769 a
2770\begin_inset Formula $\omega$
2771\end_inset
2772
2773 jsou svázány jdenoduchým diferenciálním vztahem
2774\begin_inset Formula $\frac{d\vartheta}{dt}=\omega$
2775\end_inset
2776
2777.
2778 Při praktickém užití, kdy rovnice diskretizujeme, může být ale výpočet
2779 derivace popřípadě integrálu velmi nepřesný.
2780 Dáváme tedy přednost metodám estimace těchto veličin, které nám poskytují
2781 odhad obou.
2782\end_layout
2783
2784\begin_layout Subsubsection
2785Elektrické veličiny
2786\end_layout
2787
2788\begin_layout Standard
2789Co se týče dalších (elektrických) stavových veličin systému, ve výše uvedených
2790 rovnicích vystupují ještě proudy
2791\begin_inset Formula $i$
2792\end_inset
2793
2794 a napětí
2795\begin_inset Formula $u$
2796\end_inset
2797
2798.
2799 Proudy
2800\begin_inset Formula $i$
2801\end_inset
2802
2803 předpokládáme, že měříme, samozřejmě jen s určitou přesností.
2804 Napětí
2805\begin_inset Formula $u$
2806\end_inset
2807
2808 pak jsou vstupy, kterými řídíme systém.
2809 Ty navrhujeme a tedy je předpokládáme známé, je však třeba uvést, že řízením
2810 navržená napětí
2811\begin_inset Formula $u$
2812\end_inset
2813
2814 nejdou přímo do motoru, ale slouží pouze jako referenční hodnoty pro napájecí
2815 zdroj.
2816 Kontrolu nad napětím na vstupu do motoru tedy nemáme.
2817\end_layout
2818
2819\begin_layout Subsubsection
2820Bezsenzorové řízení
2821\end_layout
2822
2823\begin_layout Standard
2824Dále se v textu hovoří o
2825\emph on
2826bezsenzorovém řízení
2827\emph default
2828.
2829 Pod tímto pojmem je vždy bezvýhradně myšleno řízení, které nevyužívá senzorů
2830 k měření mechanických veličin.
2831 Elektrické veličiny jsou měřeny vždy.
2832\end_layout
2833
2834\begin_layout Subsection
2835Senzorové metody
2836\end_layout
2837
2838\begin_layout Subsubsection*
2839Senzory
2840\end_layout
2841
2842\begin_layout Standard
2843Nejpřímočařejším přístupem pro určování mechanických veličin je osazení
2844 stroje senzory.
2845 Často se může jednat o pulzní snímače na principu vhodného kódu
2846\begin_inset CommandInset citation
2847LatexCommand cite
2848key "novak2006"
2849
2850\end_inset
2851
2852.
2853 Další možností je využití Hallových senzorů
2854\begin_inset CommandInset citation
2855LatexCommand cite
2856key "PUK1"
2857
2858\end_inset
2859
2860.
2861 Využití senzorů přináší obecně mnoho nevýhod.
2862 Přidává do zařízení další části a tím zvyšuje jeho cenu i poruchovost.
2863 Je třeba řešit jeho připojení k motoru a vodiče pro sběr dat.
2864 Řízení využívající senzory je méně robustní a v případě selhání senzoru
2865 ztrácíme nad strojem kontrolu.
2866 To může být nežádoucí obvzláště, je-li motor využíván současně i jako brzda
2867 
2868\begin_inset CommandInset citation
2869LatexCommand cite
2870key "PCW1"
2871
2872\end_inset
2873
2874.
2875 
2876\end_layout
2877
2878\begin_layout Standard
2879Je tedy snaha se užití senzorů vyhnout a k určování polohy a otáček rotoru
2880 využít jiných,
2881\emph on
2882bezsenzorových
2883\emph default
2884, metod.
2885\end_layout
2886
2887\begin_layout Subsubsection*
2888Rezolvery
2889\end_layout
2890
2891\begin_layout Standard
2892Podle
2893\begin_inset CommandInset citation
2894LatexCommand cite
2895key "novak2006"
2896
2897\end_inset
2898
2899 a
2900\begin_inset CommandInset citation
2901LatexCommand cite
2902key "PAH1"
2903
2904\end_inset
2905
2906 se jedná o v praxi často využívaná zařízení k vyhodnocení úhlu natočení
2907 rotoru PMSM.
2908 Rezolver je speciální servomechanismus, v podstatě střídavý stroj.
2909 Pracuje na principu polohového transformátoru.
2910 Na rotoru má umístěné bezkontaktně napájené budící vinutí (primární vinutí
2911 transformátoru).
2912 Na statoru dvě vinutí posunutá o
2913\begin_inset Formula $90^{\circ}$
2914\end_inset
2915
2916 (představují sekundární vinutí).
2917 Zařízení je napájeno vysokofrekvenčním napětím okolo
2918\begin_inset Formula $5-10kHz$
2919\end_inset
2920
2921 o malé amplitudě cca
2922\begin_inset Formula $5V$
2923\end_inset
2924
2925.
2926 Velikosti napětí indukovaných ve statorovách vinutích jsou závislé na úhlovém
2927 natočení rotoru (
2928\begin_inset Formula $\sin$
2929\end_inset
2930
2931 a
2932\begin_inset Formula $\cos$
2933\end_inset
2934
2935).
2936 To následně může být získáno například pomocí fázového závěsu.
2937\end_layout
2938
2939\begin_layout Standard
2940Rezolvery jsou robustní a vyhodnocují přesně úhel natočení, toho se využívá
2941 například v robotice.
2942 Je však třeba složitějších obvodů, pro samotné vyhodnocení.
2943 Velkou nevýhodou ale je, že se jedná o přídavné zařízení a s tím jsou spojeny
2944 problémy již zmiňované u senzorů.
2945\end_layout
2946
2947\begin_layout Standard
2948Dále se ještě nabízí otázka, proč místo užití rezolvéru přímo nepoužít vysokofre
2949kvenční signál v samotném PMSM v rámci některé z injektážních metod.
2950\end_layout
2951
2952\begin_layout Subsection
2953Zpětné elektromotorické síly
2954\end_layout
2955
2956\begin_layout Standard
2957Využítí zpětné elektromotorické síly (
2958\emph on
2959back electromotiric force, back-EMF
2960\emph default
2961) je metoda, kdy informaci o úhlu natočení a otáčkách rotoru získáváme z
2962 indukovaného napětí.
2963 Princip je v podstatě velmi jednoduchý a nejlépe je vidět na rovnicích
2964 pro proudy v souřadnicích
2965\begin_inset Formula $\alpha-\beta$
2966\end_inset
2967
2968, které představují přímý vztah mezí řízením systému na vstupu a měřenými
2969 výstupu:
2970\begin_inset Formula \begin{eqnarray*}
2971\frac{di_{\alpha}}{dt} & = & -\frac{R_{s}}{L_{s}}i_{\alpha}\left|\underline{\overline{+\frac{\psi_{pm}}{L_{s}}\omega\sin\vartheta}}\right|+\frac{u_{\alpha}}{L_{s}},\\
2972\frac{di_{\beta}}{dt} & = & -\frac{R_{s}}{L_{s}}i_{\beta}\left|\underline{\overline{-\frac{\psi_{pm}}{L_{s}}\omega\cos\vartheta}}\right|+\frac{u_{\beta}}{L_{s}},\end{eqnarray*}
2973
2974\end_inset
2975
2976kde právě zarámované členy odpovídají indukovaným napětím a je z nich možno
2977 získat hodnoty
2978\begin_inset Formula $\vartheta$
2979\end_inset
2980
2981 a
2982\begin_inset Formula $\omega$
2983\end_inset
2984
2985.
2986 V ideálním případě by stačilo pouze členy extrahovat
2987\begin_inset Formula $e_{\alpha}=\frac{\psi_{pm}}{L_{s}}\omega\sin\vartheta$
2988\end_inset
2989
2990 a
2991\begin_inset Formula $e_{\beta}=-\frac{\psi_{pm}}{L_{s}}\omega\cos\vartheta$
2992\end_inset
2993
2994 a vypočítat
2995\begin_inset Formula \begin{eqnarray*}
2996\vartheta & = & \arctan\left(-\frac{e_{\alpha}}{e_{\beta}}\right),\\
2997\left|\omega\right| & = & \frac{L_{s}}{\psi_{pm}}\sqrt{e_{\alpha}^{2}+e_{\beta}^{2}}.\end{eqnarray*}
2998
2999\end_inset
3000
3001
3002\end_layout
3003
3004\begin_layout Subsubsection
3005Komplikace
3006\end_layout
3007
3008\begin_layout Standard
3009Ve skutečnosti ale postup není tak jednoduchý.
3010 Jednak je třeba ještě vyřešit problém se znaménkem
3011\begin_inset Formula $\mathrm{sign\,}\omega$
3012\end_inset
3013
3014, protože uvedené rovnice jsou symetrické na substituci
3015\begin_inset Formula $\left(\omega,\:\vartheta\right)\longleftrightarrow\left(-\omega,\:\vartheta+\pi\right)$
3016\end_inset
3017
3018.
3019 
3020\end_layout
3021
3022\begin_layout Standard
3023Dále do systému vstupuje šum a při malém odstupu signálu od šumu (
3024\emph on
3025S/N
3026\emph default
3027) bude výpočet výše značně nepřesný.
3028 To také souvisí dalším, největším, problémem tohoto přístupu.
3029 Zatímco amplitudu šumu uvažujeme neměnnou, amplituda indukovaných napětí
3030 je přímo závislá na otáčkách stroje
3031\begin_inset Formula $\omega$
3032\end_inset
3033
3034.
3035 A tedy při nízkých, nebo dokonce nulových, otáčkách tato metoda naprosto
3036 selhává.
3037 Tento případ je o to závažnější, že se s ním musíme vyrovnat při každém
3038 rozjezdu stroje.
3039 Úhel natočení
3040\begin_inset Formula $\vartheta$
3041\end_inset
3042
3043 je tedy v tomto případě nepozorovatelný stav.
3044 Navíc nemůžeme předpokládat žádnou počáteční hodnotu, protože nám s rotorem
3045 mohl otočit nějaký vnější zásah, popřípadě mohl oddriftovat.
3046 Je tedy vhodné předpokládat počáteční natočení
3047\begin_inset Formula $\vartheta_{0}$
3048\end_inset
3049
3050 jako náhodnou veličinu s rovnoměrným rozdělením v intervalu
3051\begin_inset Formula $\left(-\pi,\:\pi\right\rangle $
3052\end_inset
3053
3054.
3055 
3056\begin_inset Float figure
3057wide false
3058sideways false
3059status collapsed
3060
3061\begin_layout Plain Layout
3062\begin_inset Graphics
3063        filename badekfestim.eps
3064        scale 60
3065
3066\end_inset
3067
3068
3069\begin_inset Caption
3070
3071\begin_layout Plain Layout
3072Výsledek odhadování stavu pomocí EKF, který předpokládá počáteční hodnotu
3073 
3074\begin_inset Formula $\vartheta_{0}=0$
3075\end_inset
3076
3077, zatímco skutečná hodnota je jiná (viz legenda).
3078 Nahoře odhady otáček
3079\begin_inset Formula $\omega_{est}$
3080\end_inset
3081
3082 (červená přerušovaná čára značí referenční hodnotu
3083\begin_inset Formula $\overline{\omega}$
3084\end_inset
3085
3086 témeř přesně sledovanou systémem s řízením se znalostí stavu, tj.
3087 
3088\begin_inset Formula $\overline{\omega}\approx\omega_{sys}$
3089\end_inset
3090
3091).
3092 Vlevo dole skutečné hodnoty úhlu natočení
3093\begin_inset Formula $\vartheta_{sys}$
3094\end_inset
3095
3096 a vpravo dole estimované hodnoty
3097\begin_inset Formula $\vartheta_{est}$
3098\end_inset
3099
3100.
3101\end_layout
3102
3103\end_inset
3104
3105
3106\begin_inset CommandInset label
3107LatexCommand label
3108name "Flo:badekfestim"
3109
3110\end_inset
3111
3112
3113\end_layout
3114
3115\end_inset
3116
3117 Příkladem toho, jaké výsledky můžeme dosáhnout, když počítáme s očekávanou
3118 hodnotou
3119\begin_inset Formula $\vartheta_{0}=0$
3120\end_inset
3121
3122, zatímco skutečná hodnota je jiná zobrazuje obrázek
3123\begin_inset CommandInset ref
3124LatexCommand ref
3125reference "Flo:badekfestim"
3126
3127\end_inset
3128
3129.
3130 Jedná se o jednoduchý příklad odhadování stavu pomocí rozšířeného Kalmanova
3131 filtru, v tomto případě neuvažujeme šum.
3132 Zde však bylo použito odhadování stavových veličin již běžícího systému,
3133 který je řízen regulátorem využívajícím přesnou informaci o stavu.
3134 Získaný odhad se tedy nevyužíval pro řízení.
3135 Když bychom řídili na základě odhadu stavu, tj.
3136 přidali do systému zpětnou vazbu, výsledek by se nepatrně zlepšil viz obrázek
3137 
3138\begin_inset CommandInset ref
3139LatexCommand ref
3140reference "Flo:badekfestim2"
3141
3142\end_inset
3143
3144.
3145\begin_inset Float figure
3146wide false
3147sideways false
3148status collapsed
3149
3150\begin_layout Plain Layout
3151\begin_inset Graphics
3152        filename badekfestim2.eps
3153        scale 60
3154
3155\end_inset
3156
3157
3158\end_layout
3159
3160\begin_layout Plain Layout
3161\begin_inset Caption
3162
3163\begin_layout Plain Layout
3164Výsledek odhadování a řízení stavu pomocí EKF, který předpokládá počáteční
3165 hodnotu
3166\begin_inset Formula $\vartheta_{0}=0$
3167\end_inset
3168
3169, zatímco skutečná hodnota je jiná (viz legenda).
3170 Nahoře průběhy skutečných otáček systému
3171\begin_inset Formula $\omega_{sys}$
3172\end_inset
3173
3174 (červená přerušovaná čára značí referenční hodnotu
3175\begin_inset Formula $\overline{\omega}$
3176\end_inset
3177
3178 témeř přesně sledovanou řízením z estimátoru, tj.
3179 
3180\begin_inset Formula $\overline{\omega}\approx\omega_{est}$
3181\end_inset
3182
3183 pro všechny volby
3184\begin_inset Formula $\vartheta_{0}$
3185\end_inset
3186
3187).
3188 Vlevo dole skutečné hodnoty úhlu natočení
3189\begin_inset Formula $\vartheta_{sys}$
3190\end_inset
3191
3192 a vpravo dole estimované hodnoty
3193\begin_inset Formula $\vartheta_{est}$
3194\end_inset
3195
3196.
3197\end_layout
3198
3199\end_inset
3200
3201
3202\begin_inset CommandInset label
3203LatexCommand label
3204name "Flo:badekfestim2"
3205
3206\end_inset
3207
3208
3209\end_layout
3210
3211\begin_layout Plain Layout
3212
3213\end_layout
3214
3215\end_inset
3216
3217
3218\end_layout
3219
3220\begin_layout Subsubsection
3221Metody
3222\end_layout
3223
3224\begin_layout Standard
3225V praxi se pro určování parametrů z inukovaných napětí dle
3226\begin_inset CommandInset citation
3227LatexCommand cite
3228key "PAH1"
3229
3230\end_inset
3231
3232 nejčastěji používají nelineární pozorovatelé nebo adaptivní řízení s referenční
3233m modelem (MRAC).
3234 Nejčasteji užívaným nelineárním pozorovatelem je pak rozšířený Kalmanův
3235 filtr (
3236\emph on
3237EKF
3238\emph default
3239).
3240 Přístupy založené na EKF lze nalézt například v
3241\begin_inset CommandInset citation
3242LatexCommand cite
3243key "PEB1,PEB2,PSB1"
3244
3245\end_inset
3246
3247.
3248 
3249\end_layout
3250
3251\begin_layout Standard
3252V
3253\begin_inset CommandInset citation
3254LatexCommand cite
3255key "PSB1"
3256
3257\end_inset
3258
3259 představují bezsenzorové řízení založené na EKF estimátoru ve spojení s
3260 PI regulátory.
3261 To nepotřebuje znát počáteční natočení rotoru a zátěžný moment.
3262 PI regulátor napětí lze nastavit se zamčeným rotorem a je řešen i problém
3263 s rozpoznáním
3264\begin_inset Formula $\mathrm{sign}\,\omega$
3265\end_inset
3266
3267.
3268 Článek
3269\begin_inset CommandInset citation
3270LatexCommand cite
3271key "PEB2"
3272
3273\end_inset
3274
3275 je také zaměřen na využití EKF, nyní však v případe IPMSM.
3276 Návrh je komplikovanější v důsledku anizotropie stroje, autoři se ji však
3277 snaží využít k vylepšení výkonu systému.
3278 
3279\end_layout
3280
3281\begin_layout Standard
3282Dále kromě EKF je možno použít například klouzavého pozorovatele (
3283\emph on
3284sliding mode observer, SMO
3285\emph default
3286), jeho iterativní verzi využívají v
3287\begin_inset CommandInset citation
3288LatexCommand cite
3289key "PSK1"
3290
3291\end_inset
3292
3293.
3294 V
3295\begin_inset CommandInset citation
3296LatexCommand cite
3297key "PLU1"
3298
3299\end_inset
3300
3301 využívají také řízení na klouzavém pozorovateli, kde si ale navíc při nízkých
3302 otáčkách
3303\begin_inset Formula $\omega\thickapprox0$
3304\end_inset
3305
3306 pomáhají injektováním stejnosměrného proudu do
3307\begin_inset Formula $d$
3308\end_inset
3309
3310 osy.
3311 Nevyužívají však anizotropií ani nijak zvlášť neanalyzují injektovaný signál,
3312 tento přístup tedy v textu nezařadíme mezi injektáže.
3313\end_layout
3314
3315\begin_layout Standard
3316Pod metody využívající informaci ze zpětné elektromagnetické síly, můžeme
3317 zařadit ještě mnoho dalších, které možná na první pohled do této kategorie
3318 nespadají.
3319 Především se jedná o metody snažící se nějakým způsobem odstranit šum a
3320 tedy zvýšit rozlišovací schopnost indukovaných napětí.
3321 Opět zde narážíme na problém, že nefungují při
3322\begin_inset Formula $\omega\equiv0$
3323\end_inset
3324
3325.
3326 Jedná se o různé podoby odšumovacích filtrů, tedy filtrů typu dolní propusť
3327 (low-pass).
3328 V časové oblasti můžeme použít například klouzavé průměry (moving averages
3329 - MA) nebo jejich váženou verzi.
3330 Ve frekvenční oblasti lze užít (klouzavé) diskrétní Fouriefovy transformace,
3331 a buď odstranit vyšší frekvence, nebo si vybrat jen nějakou nízkou.
3332 Tím však nezískáváme o moc navíc, protože 0.
3333 harmonická odpovídá v podstatě průměru, další harmonické pak vhodně váženému
3334 průměru.
3335\end_layout
3336
3337\begin_layout Standard
3338Za zmínku ještě stojí další skupina metod využívající více paralelně běžících
3339 odhadů z nichž vybírá jeden, nějakým způsobem optimální.
3340 Takovou metodou je například sekvenční Monte Carlo metoda (Particle Filter).
3341 
3342\end_layout
3343
3344\begin_layout Standard
3345Dalším příkladem by mohlo být více paralelně běžících modelů, z nichž se
3346 vybere ten, jehož výstup nejlépe odpovídá výstupu skutečného systému.
3347 Nedostatkem těchto přístupů je poměrně velká výpočetní náročnost, přesto
3348 ale poskytují relativně dobré výsledky.
3349 Příkladem může být obrázek
3350\begin_inset CommandInset ref
3351LatexCommand ref
3352reference "Flo:multipleest"
3353
3354\end_inset
3355
3356 zachycujcí výsledek
3357\begin_inset Formula $12$
3358\end_inset
3359
3360 současně běžících modelů s různým počátečním odhadem
3361\begin_inset Formula $\vartheta_{0}$
3362\end_inset
3363
3364 rovnoměrně rozloženým v intervalu
3365\begin_inset Formula $\left(-\pi,\:\pi\right\rangle $
3366\end_inset
3367
3368.
3369 Skutečná počáteční hodnota natočení rotoru systému je
3370\begin_inset Formula $\vartheta_{0}=\frac{5}{12}\pi$
3371\end_inset
3372
3373.
3374\begin_inset Float figure
3375wide false
3376sideways false
3377status collapsed
3378
3379\begin_layout Plain Layout
3380\begin_inset Graphics
3381        filename multiest3.eps
3382        scale 60
3383
3384\end_inset
3385
3386
3387\end_layout
3388
3389\begin_layout Plain Layout
3390\begin_inset Caption
3391
3392\begin_layout Plain Layout
3393Grafy znázorňují průběh otáček
3394\begin_inset Formula $\omega$
3395\end_inset
3396
3397 a úhlů natočení systému
3398\begin_inset Formula $\theta$
3399\end_inset
3400
3401, kdy je k odhadování stavu použito
3402\begin_inset Formula $12$
3403\end_inset
3404
3405 současně běžících modelů, z nichž je v každém kroku vybírán nejlepší na
3406 základě shody s výstupem (měřené proudy) skutečného systému.
3407 Systém je řízen ze stavu, aby co nejlépe sledoval požadovanou hodnotu otáček,
3408 jeho počáteční úhel natočení je
3409\begin_inset Formula $\vartheta_{0}=\frac{5}{12}\pi$
3410\end_inset
3411
3412.
3413 Počáteční
3414\begin_inset Formula $\vartheta_{0}$
3415\end_inset
3416
3417 odhadovacích modelů jsou rovnoměrně rozloženy v intervalu
3418\begin_inset Formula $\left(-\pi,\:\pi\right\rangle $
3419\end_inset
3420
3421.
3422\end_layout
3423
3424\end_inset
3425
3426
3427\begin_inset CommandInset label
3428LatexCommand label
3429name "Flo:multipleest"
3430
3431\end_inset
3432
3433
3434\end_layout
3435
3436\begin_layout Plain Layout
3437
3438\end_layout
3439
3440\end_inset
3441
3442
3443\end_layout
3444
3445\begin_layout Subsubsection
3446Další vlastnosti
3447\end_layout
3448
3449\begin_layout Standard
3450Metody využívající zpětnou elektromotorickou sílu jsou obvykle založeny
3451 na modelu a je tedy důležitá znalost parametrů stroje.
3452 Bylo by tedy dobré najít přístupy, které na parametrech nezávisí, popřípadě
3453 které jsou odolné na jejich změnu.
3454 To se daří u mechanických parametrů stroje, jako je zátěžný moment například
3455 v
3456\begin_inset CommandInset citation
3457LatexCommand cite
3458key "PEB2,PSB1"
3459
3460\end_inset
3461
3462.
3463 
3464\end_layout
3465
3466\begin_layout Standard
3467Ve vyšších otáčkách poskytuje tento přístup dobré výsledky.
3468 Proto je součástí hybridních metod, které kombinují využití zpětné elektromotor
3469ické síly a injektáž.
3470\end_layout
3471
3472\begin_layout Subsection
3473Injektáže
3474\end_layout
3475
3476\begin_layout Standard
3477Injektážemi označujeme v textu metody, které využívají přídavného signálu
3478 k detekci anizotropií stroje a usnadňují určení jeho jinak obtížně pozorovateln
3479ých stavů, především úhlu natočení
3480\begin_inset Formula $\vartheta$
3481\end_inset
3482
3483.
3484 Anizotropie lze rodělit do dvou hlavních kategorií.
3485 První jsou vlastní magnetické výčnělky (
3486\emph on
3487saliency
3488\emph default
3489) rotoru, ty jsou charakteristické především pro IPMSM.
3490 Do druhé kategorie pak spadají lokální anizotropie vzniklé saturací magnetickým
3491 tokem, typické pro SMPMSM.
3492 Signál je přiváděn na vstup stroje spolu s řízením.
3493 Zpravidla je využíván vysokofrekvenční signál, aby docházelo k co možná
3494 nejmenšímu narušení průběhu samotného řízení.
3495 Tyto metody jsou ale ve většině případů založeny na nějakém
3496\begin_inset Quotes gld
3497\end_inset
3498
3499speciálním
3500\begin_inset Quotes grd
3501\end_inset
3502
3503 jevu (anizotropii), v tom smyslu, že jej v základních rovnicích nemáme.
3504 V reálném zařízení se samozřejmě vyskytují.
3505 
3506\end_layout
3507
3508\begin_layout Standard
3509Nejobvyklejším přístupem je, že anizotropie je v podstatě reprezentována
3510 rozdílnými indukčnostimi v osách
3511\begin_inset Formula $d$
3512\end_inset
3513
3514 a
3515\begin_inset Formula $q$
3516\end_inset
3517
3518, tedy
3519\begin_inset Formula $L_{d}\neq L_{q}$
3520\end_inset
3521
3522.
3523 Pro IPMSM s permanentními magnety uvnitř rotoru toto platí relativně velmi
3524 dobře.
3525 V případě SMPMSM je však situace horší, protože rozdíl
3526\begin_inset Formula $L_{q}-L_{d}$
3527\end_inset
3528
3529 je velmi malý, v krajním případě dokonce nulový.
3530 Za předpokladu
3531\begin_inset Formula $L_{d}\neq L_{q}$
3532\end_inset
3533
3534 lze této vlastnosti využít k určení polohy (úhlu natočení) rotoru
3535\begin_inset Formula $\vartheta$
3536\end_inset
3537
3538 pomocí injektování vhodného testovacího signálu do stroje.
3539 Obvykle se využívá vysokofrekvenčního signálu o frekvenci v řádu stovek
3540 
3541\begin_inset Formula $Hz$
3542\end_inset
3543
3544.
3545 Existují však i injektáže využívající nízkofrekvenční signály.
3546 
3547\end_layout
3548
3549\begin_layout Standard
3550Injektáž je aplikována jako vysokofrekvenční napěťový harmonický signál
3551 o frekvenci přibližně
3552\begin_inset Formula $500\, Hz$
3553\end_inset
3554
3555.
3556 Ten je injektovaný do estimované osy
3557\begin_inset Formula $d$
3558\end_inset
3559
3560 spolu s řídícím napětím.
3561 Následně je získána v
3562\begin_inset Formula $q$
3563\end_inset
3564
3565 složce proudu informace o úhlu natočení
3566\begin_inset Formula $\sin2\hat{\vartheta}$
3567\end_inset
3568
3569.
3570 Jedná se o obálku amplitudově modulovanou na nosné frekvenci.
3571 Demodulace je provedena vynásobením vysokofrekvenčním nosným signálem a
3572 následným užitím low-pass filtru.
3573 Je však třeba upozornit na nutnost předpokladu
3574\begin_inset Formula $L_{d}\neq L_{q}$
3575\end_inset
3576
3577, protože amplituda
3578\begin_inset Formula $\sin2\hat{\vartheta}$
3579\end_inset
3580
3581 závisí mimo jiné přímo úměrně na rozdílu indukčností
3582\begin_inset Formula $L_{q}-L_{d}$
3583\end_inset
3584
3585.
3586\end_layout
3587
3588\begin_layout Standard
3589Tato základní metoda je užívána například v
3590\begin_inset CommandInset citation
3591LatexCommand cite
3592key "PAH1,PAB1"
3593
3594\end_inset
3595
3596.
3597 Dále pak v
3598\begin_inset CommandInset citation
3599LatexCommand cite
3600key "PSJ1"
3601
3602\end_inset
3603
3604, kde se využívá principu, kdy v důsledku magnetického toku permanentních
3605 magnetů je syceno jádro vinutí kolem
3606\begin_inset Formula $q$
3607\end_inset
3608
3609 osy.
3610 To vytváří magnetickou nepravidelnost v motoru závislou na poloze rotoru.
3611 Tato nepravidelnost je následně detekována injektovaným vysokofrekvenčním
3612 napětím.
3613 Výhodou této metody je, že je přímo navrhována pro užití v SMPMSM.
3614 Vysokofrekvenční napěťový signál je opět injektován do estimované
3615\begin_inset Formula $d$
3616\end_inset
3617
3618 osy, informace o poloze rotoru je získána z proudu v ose
3619\begin_inset Formula $q$
3620\end_inset
3621
3622 násobením a low-pass filtrem.
3623 
3624\end_layout
3625
3626\begin_layout Standard
3627V článku
3628\begin_inset CommandInset citation
3629LatexCommand cite
3630key "PCB1"
3631
3632\end_inset
3633
3634 se zabývají srovnáním dvou metod injektáží.
3635 Zaměřují se jak na IPMSM, který má větší rozdíl indukčností
3636\begin_inset Formula $L_{q}-L_{d}$
3637\end_inset
3638
3639, tak i na SMPMSM.
3640 První metoda označovaná jako
3641\begin_inset Quotes gld
3642\end_inset
3643
3644pulzující napěťový vektor
3645\begin_inset Quotes grd
3646\end_inset
3647
3648 je v podstatě shodná s injektážní technikou z minulého odstavce.
3649 Oproti tomu druhý způsob, nazývaný jako
3650\begin_inset Quotes gld
3651\end_inset
3652
3653rotující napěťový vektor
3654\begin_inset Quotes grd
3655\end_inset
3656
3657, užívá injektáž v souřadnicích
3658\begin_inset Formula $\alpha-\beta$
3659\end_inset
3660
3661.
3662 Informaci o úhlu natočení, respektive chybě odhadu úhlu natočení je pak
3663 získána násobením a následnou aplikací high-pass filtru.
3664 Opět ale platí, že získaná informace je úměrná rozdílu indukčností
3665\begin_inset Formula $L_{q}-L_{d}$
3666\end_inset
3667
3668.
3669 Dále je v článku provedeno srovnání obou metoda na oba typy motorů, kdy
3670 je užit stejný stator a měněny rotory (SMPMSM a IPMSM).
3671 Mezi injektážními metodami nebyl shledán žádný zásadnější rozdíl.
3672 Rozdíly se projevily spíše při použití stejné metody na různé motory, to
3673 souvisí s jejich magnetickými vlastnostmi, v tomto textu se tímto však
3674 zabývat nebudeme.
3675 Srovnáním zmiňovaných dvou metod se zabývají i v
3676\begin_inset CommandInset citation
3677LatexCommand cite
3678key "PCK1"
3679
3680\end_inset
3681
3682, zaměřují se však na IPMSM.
3683\end_layout
3684
3685\begin_layout Standard
3686Články
3687\begin_inset CommandInset citation
3688LatexCommand cite
3689key "PSL1,PSL3"
3690
3691\end_inset
3692
3693 představují injektážní metodu k detekci anizotropií, která nepotřebuje
3694 znát parametry stroje.
3695 V případě
3696\begin_inset CommandInset citation
3697LatexCommand cite
3698key "PSL3"
3699
3700\end_inset
3701
3702 se navíc snaží kompenzovat negativní vliv invertoru, především jev označovaný
3703 jako
3704\emph on
3705dead-time effect
3706\emph default
3707.
3708 Díky tomu jsou schopni detekovat i malé nepravidelnosti typické pro SMPMSM.
3709 Je užíván vysokofrekvenční napěťový signál o frekvenci okolo
3710\begin_inset Formula $2\, kHz$
3711\end_inset
3712
3713.
3714 Injektovaný signál je složením dvou signálů rotojících proti sobě.
3715 V případě špatného odhadu úhlu
3716\begin_inset Formula $\hat{\vartheta}\neq\vartheta$
3717\end_inset
3718
3719 je vzniká aditivní vysokofrekvenční signál v proudech, ze kterého může
3720 být tato chyba získána pomocí pozorovatele
3721\emph on
3722(Tracking Observer
3723\emph default
3724).
3725\end_layout
3726
3727\begin_layout Standard
3728Zajímavou techniku představují v
3729\begin_inset CommandInset citation
3730LatexCommand cite
3731key "PIS1"
3732
3733\end_inset
3734
3735, kde vypočítají absolutní polohu rotoru v klidu.
3736 Metoda funguje i pro SMPMSM a je založena na injektování vhodných napěťových
3737 pulzů do vinutí každé z fází.
3738 Následně dochází k částečnému nasycení statoru, ze kterého je možno spočítat
3739 absolutní polohu rotoru i bez znalosti parametrů stroje.
3740\end_layout
3741
3742\begin_layout Standard
3743Další velmi zajímavý přístup je prezentován v
3744\begin_inset CommandInset citation
3745LatexCommand cite
3746key "PAP1"
3747
3748\end_inset
3749
3750.
3751 Tato metoda nevyužívá anizotropií rotoru, ani výčnělků, místo toho je založena
3752 na anizotropii samotných permanentních magnetů.
3753 Z tohoto důvodu může být dobře využita při estimaci PMSM, kde ostatní metody
3754 selhávají, například z důvodu
3755\begin_inset Formula $L_{q}=L_{d}$
3756\end_inset
3757
3758.
3759 K jejich detekci je ale třeba využít velmi vysokých frekvencí, řádově
3760\begin_inset Formula $100-500\, kHz$
3761\end_inset
3762
3763.
3764 Optimální hodnotu frekvence je navíc třeba naladit pro konkrétní typ magnetu.
3765 Tento přístup vypadá velmi slibně, ale jak autoři sami uvádějí, je tato
3766 metoda nová a vyvstává kolem ní ještě mnoho nezodpovězených otázek.
3767\end_layout
3768
3769\begin_layout Subsection
3770Hybridní metody
3771\end_layout
3772
3773\begin_layout Standard
3774Hybridními metodami v textu označujeme v podstatě vhodnou kombinaci předchozích
3775 dvou zmiňovaných přístupů.
3776 Techniky založené na zpětné elektromotorické síle fungují relativně velmi
3777 dobře, selhávají ale při nízkých a nulových otáčkách.
3778 Naopak užití injektáží je vhodné pro nízké a nulové otáčky, zatímco ve
3779 vyšších rychlostech způsobuje nežádoucí rušení.
3780 Z tohoto důvodu je snaha obě metody vhodným způsobem zkombinovat a využít
3781 předností obou.
3782 Základní idea je tedy jednoduchá.
3783 Dokud se pohybujeme v nízkých otáčkách, využíváme odhadů založených na
3784 injektáži, při vyšších otáčkách injektáž vypneme, aby nezpůsobovala nežádoucí
3785 zásahy a užíváme již jen odhadů získaných ze zpětné elektromotorické síly.
3786\end_layout
3787
3788\begin_layout Standard
3789Tento postup je použit například v
3790\begin_inset CommandInset citation
3791LatexCommand cite
3792key "PAP2"
3793
3794\end_inset
3795
3796, kdy jako estimátor používají adaptivního pozorovatele s referenčním modelem,
3797 který je pro nízké otáčky doplněn injektáží v podstatě v základním návrhu
3798 popsaném v předcházející části.
3799 Důležitou součástí těchto metod je způsob, jakým se vyřeší
3800\begin_inset Quotes gld
3801\end_inset
3802
3803bezproblémový
3804\begin_inset Quotes grd
3805\end_inset
3806
3807 přechod z jednoho estimátoru na jiný.
3808 V
3809\begin_inset CommandInset citation
3810LatexCommand cite
3811key "PHS1"
3812
3813\end_inset
3814
3815 je to například řešeno tak, že užívají estimátor rotorového toku založený
3816 na indukovaných napětích, který je funkční pořád.
3817 V nízkých otáčkách je pak doplňován injektáží, ta s rostoucími otáčkami
3818 postupně vymizí.
3819\end_layout
3820
3821\begin_layout Section
3822Řízení
3823\end_layout
3824
3825\begin_layout Standard
3826oddělit, potřeba dobrého odhadu, řízení v
3827\begin_inset Formula $\alpha-\beta$
3828\end_inset
3829
3830 oprodi
3831\begin_inset Formula $d-q$
3832\end_inset
3833
3834
3835\end_layout
3836
3837\begin_layout Subsection*
3838řídící strategii
3839\end_layout
3840
3841\begin_layout Standard
3842návrh standartně PI (vektorové), nebo přes LQ, zmínit DTC
3843\end_layout
3844
3845\begin_layout Subsection*
3846současný stav
3847\end_layout
3848
3849\begin_layout Standard
3850nejlepší je hybridní, ale třeba přepínat více modelů, řízení PI
3851\end_layout
3852
3853\begin_layout Subsection*
3854duální přístup
3855\end_layout
3856
3857\begin_layout Standard
3858výhody duálního přístupu, proč se na to laicky hodí, problém s reálným časem,
3859 jednoduché metody
3860\end_layout
3861
3862\begin_layout Section
3863Duální řízení
3864\end_layout
3865
3866\begin_layout Standard
3867stručně popis, proč jednoduché, jaké? - třeba filatov...
3868\end_layout
3869
3870\begin_layout Chapter
3871snaha o návrh
3872\end_layout
3873
3874\begin_layout Standard
3875injektáž-závěs-klaman-lq
3876\end_layout
3877
3878\begin_layout Chapter
3879vyhodnoncení a simulace
3880\end_layout
3881
3882\begin_layout Standard
3883možná něco vlastního v matlabu
3884\end_layout
3885
3886\begin_layout Standard
3887závěry ze simulátoru
3888\end_layout
3889
3890\begin_layout Standard
3891hlavně otestování toho
3892\begin_inset Quotes eld
3893\end_inset
3894
3895snaha o návrh
3896\begin_inset Quotes erd
3897\end_inset
3898
3899
3900\end_layout
3901
3902\begin_layout Standard
3903podloženo simulacemi i z těch předchozích sekcí
3904\end_layout
3905
3906\begin_layout Standard
3907\begin_inset Newpage newpage
3908\end_inset
3909
3910
3911\end_layout
3912
3913\begin_layout Addchap
3914Závěr
3915\end_layout
3916
3917\begin_layout Standard
3918\begin_inset Newpage newpage
3919\end_inset
3920
3921
3922\end_layout
3923
3924\begin_layout Standard
3925\begin_inset CommandInset bibtex
3926LatexCommand bibtex
3927bibfiles "vyz_clanky,vyz_texty"
3928options "bibtotoc,czechiso"
3929
3930\end_inset
3931
3932
3933\end_layout
3934
3935\end_body
3936\end_document
Note: See TracBrowser for help on using the browser.