Диссертация (1024744), страница 39
Текст из файла (страница 39)
Внешний скручивающий момент может быть определен, как:M КР = qв ∆ϕ(П.2)7oПримем за расчетные значения qв = 3,8775 ⋅10 Нм, ∆ϕ = 0,066 , тогдаM КР = qв ∆ϕ = 44665,7 НмПримем, что на муфту гарантировано передается от 0,3 до 0,7 Мкр.Мощности 3,616 кВт соответствует значение крутящего момента,равного 11510 Нм.Силасдвига,возникающаяпринагружениимуфтывнешнимскручивающим моментом будет определяться выражением:FСДВИГ =M КРrШ n ,(П.3)где rШ - радиус расположения шпилек относительно оси валавозбудителя;n – количество шпилек муфты.Таким образом, с учетом полученной измерительной информациифазохронометрической системой и системой регистрации электрическихпараметров могут быть рассчитаны и оценены значения крутящихмоментов, воздействующих на наиболее опасное сечение валопровода.Статический расчет муфты РВ-РГ.302Поперечная нагрузка болта воспринимается силами трения.
Болтрассчитывается на силу затяжки [257]:2πdN= 1 [σ p ] ,Q0 =f ТР4(П.4)где N - поперечная нагрузка на одну шпильку;fТР – коэффициент трения на стыках деталей;[σ p ] - предел выносливости, принимается равным 0,6 σ Т ;d1 - диаметр шпильки в месте утонения, 20 мм.Тогда,Q0 =πd124[σ p ] =π (0,02) 24[0,6 ⋅ 650 ⋅ 10 6 ] = 122,522 кНПоперечная нагрузка на одну шпильку составляет:N = Q0 ⋅ f ТР = 18,378 кНТогда, передаваемая мощность через 12 шпилек муфты РВ-РГсоответствует максимальному значению передаваемой мощностиPmax = 12 ⋅ N ⋅ 2πν = 7,96 МВтРасчет фланцевых и фрикционных муфт на трениеМуфты турбоагрегатов ТВВ-200-2-К-200-130 могут быть отнесены ккатегории фланцевых муфт.Фланцевые муфты (в соответствии с ГОСТ 20761) из стали и чугунаобщемашиностроительного применения, предназначенные для соединениясоосных цилиндрических валов и передачи крутящего момента безуменьшения динамических нагрузок: от 16 до 40000 Нм при окружнойскорости на наружном диаметре муфт до 70 м/с - стальными муфтами, от 8до 20000 Нм при окружной скорости до 35 м/с - чугунными муфтами.
[258]В соответствии с ГОСТ 20761 приведены значения номинальныхкрутящих моментов для муфт, изготавливаемых из стали 40 (ГОСТ 1050-74).303Если нагрузка является переменной и периодически может достигатьдвукратного увеличения, значения номинального момента вращения должныбыть уменьшены в 1,4 раза. Муфте РВ-РГ соответствует номинальныймомент 6300 Нм, с учетом переменной нагрузки 4500 Нм ( 4,5 ⋅ 10 6 Н мм).Допускаемое соотношение крутящих моментов срабатывающих муфтвозьмем с учетом ГОСТ 15622, применяемое для расчета фрикционныхпредохранительных муфт:М max≤ 1,3М номгдеМ max -максимальный(П.5)допускаемыйкрутящиймомент,обеспечивающий работу муфты на трение;М ном-наибольшийдопускаемыйрабочийкрутящиймомент,передаваемый силами трения без относительного скольжения.В соответствии с [257] наиболее распространенными дефектамимонтажа жестких муфт являются:- коленчатость соединения в жесткой муфте, при монтаже или впроцессе эксплуатации по тем или иным причинам центры соединяемыхвалов не точно совпадают друг с другом;- угловой излом в жесткой муфте, приводящий к излому вала (напрактике это "маятник" или "раскрытие муфты");- различные дефекты посадки полумуфт на соединяемые валы,выражающиеся в отклонениях от условий идеального монтажа.Сила трения в муфте зависит от крутящего момента, передаваемогомуфтой, коэффициента трения материалов полумуфт и габаритных размеров.[257,258]FТР =PМ ⋅ f ТР2πωrМ ,где PМ - передаваемая мощность, Вт;fТР - коэффициент трения сталь по стали;(П.6)304rМ - радиус муфты РВ-РГ.FТР =PМ ⋅ f ТР 3616000 ⋅ 0,15== 12332 Н2πωrМ2π ⋅ 50 ⋅ 0,14Силу сжатия дисков можно определить следующим выражением1:M ⋅kнFсж =,r ⋅fм ТР(П.7)M ⋅ k 11,51 ⋅1,3 ⋅10 6нFсж === 1335 кНr ⋅f140 ⋅ 0,08м ТРРасчетмуфтывозбудителяприциклическихизмененияхкасательных напряженийАнализхронограмм(экспериментальныхданных)показал,чтовалопровод вследствие влияния внешней электрической сети постояннонаходится в режиме крутильных колебаний, при этом период вращения(частота) изменяется в определенных пределах.
Такое изменение частотывращения неизбежно приводит к изменению момента на валу роторагенератора и, как следствие, изменение момента передается на муфтывозбудителя и РНД, соединяющие генератор с соседними секциямивалопровода турбоагрегата. Момент изменяется циклически, поэтому расчетшпилек муфты возбудителя следует также проводить по циклическимизменениям касательных напряжений.Период вращения при грозовой активности может изменяться впределах ± 120 мкс, следовательно частота вращения изменяется в пределахω = (50,0 ± 0,3) Гц.Момент сопротивления всех болтов относительно центральной осиопределяется выражением:W1xвс.б.=nπd 2б ⋅y ,i4Анурьев В.И.
Справочник конструктора-машиностроителя. В трех томах. Том 2. С. 364.(П.8)305где n – количество шпилек на муфте;d - диаметр шпильки, мм;бy i - расстояние от центральной оси до шпильки.W= 12xвв.б.Wxвв.б.π ⋅ (24 ⋅ 10 − 3 ) 2 230 ⋅ 10 − 3⋅4=62π ⋅ (30 ⋅ 10 − 3 ) 2 230 ⋅ 10 − 3⋅42= 623980,8 ⋅ 10 −9 м 3= 487485 ⋅ 10 −9 м 3Так как моменты направлены в разные стороны, то на муфте возникаетразностный момент М ф = М Т − М Г , который передаётся на вал возбудителячерез муфту.Характеристики цикла соответственно равны [257,258]:τa =τm =τa =Mτm =τa =MMτm =max2 ⋅WMmax2 ⋅Wmax2 ⋅W−Mminxвс.б.+Mmin(П.9)xвс.б.3min = (11,51 − (−1,464)) ⋅10 = 10,4 МПа2 ⋅ 623980,8 ⋅10 − 9xвв.б.−Mmax2 ⋅Wmax2 ⋅WMM3min = (11,51 −1,464) ⋅10 = 8,05 МПа2 ⋅ 623980,8 ⋅10 − 9xвв.б.+M3min = (11,51 − (−1,464)) ⋅10 = 13,3 МПа2 ⋅ 487485 ⋅10 − 9xвв.б.−Mmax2 ⋅W3min = (11,51 − 1,464) ⋅ 10 = 10,3 МПа2 ⋅ 487485 ⋅10 − 9xвв.б.+MПри работе в условиях циклического нагружения коэффициент запасарассчитывается следующим выражением [257,258]:τ−1n =R K ⋅τ + ψ ⋅ττ aτ mгде τ−1- предел выносливости на кручение, Н/мм2;(П.10)306Kτ – коэффициент снижения предела выносливости при симметричномцикле;ψτ - коэффициент влияния асимметрии цикла на величину пределавыносливости.Тогда для муфты возбудителя при кручении коэффициент запаса поусталостному разрушению для материала шпилек равен:n=Rτ95−1== 4,0K ⋅τ +ψ ⋅τ1.5 ⋅13,3 + 0,05 ⋅10.3τ aτ mОбщий коэффициент безопасности с учетом изгибной жесткости равен:n =R⋅n4,0 ⋅1,73Rτ Rσ== 1,582222n Rτ + n Rσ4 + 1,73nПроведенные расчеты на прочность узла сочленения РВ-РГ показалинадежность применяемой конструкции для обеспечения неразъёмностисоединения и работы муфты РВ-РГ на трение со следующими результатами:- сборка муфты в соответствии с КД и технологическим процессом исоблюдение требований по вытяжке и расцентровкам;- по циклическому нагружению без учета изгибных колебаний(колебаний фундамента, коленчатость) шпильки проходят с хорошимкоэффициентом запаса, однако если учесть коэффициент запаса при ВКЗ наизгиб - значение 1,6 (то есть находимся на грани допустимого коэффициентазапаса), а при расчётахза основу взят метод Серенсена - Кисаношвили, атакже учтена теория Гудмана и Гербера;Расчетынациклическуюпрочностьвыполненыдляпределавыносливости из расчета 10 000 000 циклов, однако с учетом результатованализа в МГТУ им.
НЭ. Баумана и ЗАО «Уралэнерго-Союз», количествоциклов нагружения за 2 с может составить около 10 циклов, т.е. с учетомрасчетов за год (31 556 925 с) такую интенсивность можно выдержать тольков пределах от 25 до 30 дней. При этом грозовой активности региона (только307грозовые дни в Ханты-Мансийском округе) соответствует количествогрозовых дней в году в среднем 20 без учета КЗ на линиях, потребителей ит.д.308П.5. ЛИСТИНГ ПРОГРАММЫ РАСЧЕТА МАТЕМАТИЧЕСКОЙМОДЕЛИ ТУРБОАГРЕГАТА В ФАЗОХРОНОМЕТРИЧЕСКОМ ВИДЕС УЧЕТОМ ВОЗДЕЙСТВИЯ СИСТЕМЫ ВОЗБУЖДЕНИЯ ИВНЕШНЕЙ ЭЛЕКТРИЧЕСКОЙ СЕТИclear;%---------Расчет номинального режима работы генератора ТВВ-200-2---------%%---------Исходные данные - тип генератора ТВВ-200-2АУЗ------------------%P0=200000000;%----Номинальная активная мощность турбины, [Вт]-----%U0=15750/(sqrt(3)); %----Номинальное напряжение при соединении полюсов звездой, [В]--%I0=8620/(sqrt(3)); %----Номинальный ток, [А]--------%S0=235300000;%-----Номинальная полная мощность, [МВА]-----------%ra=0.00188;%-----Активное сопотивление обмотки статора, [Ом]------%rr=0.11;%-----Активное сопотивление обмотки ротора, [Ом]------%xs=0.2;%-----Индуктивное сопротивление рассеяния обмотки статора xs, [о.е.]----%xd=1.97;%-----Синхроное индуктивное сопротивление по продольной оси xd, [о.е.]---%xq=1.88;%-----Синхронное индуктивное сопротивление по поперечной оси xq,[о.е.]----%xdd=0.29;%-----Переходное сопротивление по продольной оси xd', [о.е.]----%xddd=0.2;%-----Сверхпереходное сопротивление по продольной оси xd'', [о.е.]----%xqqq=0.28;%-----Сверхпереходное сопротивление по поперечной оси xq'', [о.е.]----%x2=0.24;%-----Cопротивление обратной последовательности x2, [о.е.]----%xb=0.08;%-----Индуктивное сопротивление обмотки возбуждения xf, [о.е.]----%Tdd=0.12;%-----Сверхпереходная постоянная времени в продольной оси Td'', [c]---%Tqq=0.17;%-----Сверхпереходная постоянная времени в поперечной оси Tq'', [c]---%Ta3=0.03;%-----Постоянная времени затухания апериодической составляющей токатрехфазного короткого замыкания Ta3, [c]-%Td0=6.4;%-----Постоянная времени обмотки возбуждения при прочих разомкнутыхконтурах Td0, [c]---%Ub=300;Ib=2540;%-----Напряжение возбуждения номинальное, [В]-----%%-----Ток возбуждения номинальный, [А]-----%Ubb=89;Ibb=200;%-----Напряжение обмотки возбуждения номинальное, [В]-----%%-----Ток обмотки возбуждения номинальный, [А]-----%fi=acos(0.85);%w0=2*pi*50;%-----Угол fi, определяемый из коэффициента мощности cos(fi)=0.85----%-----Угловая скорость вращения, [рад/с]------%309%-----------------Перевод величин в систему СИ--------------------------%r0=U0/I0;%------Базовое значение сопротивления, [Ом]-----%xs_SI=xs*r0;%-----Индуктивное сопротивление рассеяния обмотки статора xs, [Ом]---%xd_SI=xd*r0;%-----Синхроное индуктивное сопротивление по продольной оси xd,[Ом]----%xq_SI=xq*r0;%-----Синхронное индуктивное сопротивление по поперечной оси xq,[Ом]----%xdd_SI=xdd*r0; %-----Переходное сопротивление по продольной оси xd', [Ом]----%xddd_SI=xddd*r0; %-----Сверхпереходное сопротивление по продольной оси xd'', [Ом]---%xqqq_SI=xqqq*r0; %-----Сверхпереходное сопротивление по поперечной оси xq'', [Ом]---%x2_SI=x2*r0;%-----Cопротивление обратной последовательности x2, [Ом]----%xb_SI=xb*r0;%-----Индуктивное сопротивление обмотки возбуждения xf, [Ом]----%Qi=atan((xq_SI*cos(fi)-ra*sin(fi))/(U0/I0+xq_SI*sin(fi)+ra*cos(fi)));%-------Внутреннийугол генератора для установившегося режима работы, [рад]-----%Qi_SI=Qi*180/pi;%-------Внутренний угол генератора для установившегося режимаработы, [град]-----%rb=Ubb/Ibb;Lb=rb*Td0;%----Активное сопротивление обмотки возбуждения, [Ом]------%%----Индуктивность обмотки возбуждения, [Гн]---------------%Id0=I0*sin(Qi+fi); %----Ток по оси d, [А]------%Iq0=I0*cos(Qi+fi); %----Ток по оси q, [А]------%Ud0=U0*sin(Qi+fi); %----Напряжение по оси d, [В]------%Uq0=U0*cos(Qi+fi); %----Напряжение по оси q, [В]------%Eq=U0*cos(Qi)+xdd_SI*Id0+ra*Iq0;%--------Синхронная ЭДС машины в поперечнойоси за переходным сопротивлением-----%kh=Eq/(w0*Ib);%--------Коэффициент взаимоиндуктивности, [Гн]--------------------%Mds=(kh*Ib*Iq0-(xd_SI-xq_SI)*Id0*Iq0/w0);%------ Момент сопротивления-----------%l1=10.205;l2=17.264;l3=3.84;J1=5600.1474;J2=7071.48;J3=1536;J4=663.22;J5=306.334;%-------Длина генератора, [м]------%%-------Длина турбины, [м]------%%-------Длина возбудителя, [м]------%%-------Момент инерции генератора, [кг*м^2]------%%-------Момент инерции ЦНД, [кг*м^2]------%%-------Момент инерции ЦСД, [кг*м^2]------%%-------Момент инерции ЦВД, [кг*м^2]------%%-------Момент инерции возбудителя, [кг*м^2]------%q1=1.5823e+008;%-------Жесткость генератора, [Нм]------%q2=6.5736e+009;%-------Жесткость ЦНД, [Нм]------%q3=2.123312111218046e+09;%-------Жесткость ЦСД, [Нм]------%q4=1.205875445746307e+09;%-------Жесткость ЦВД, [Нм]------%310q5=3.8775e+007;qM=166800000;%-------Жесткость возбудителя, [Нм]------%%-------Жесткость муфты, [Нм]------%q12=((q1*qM)/(q1+qM))*q2/(((q1*qM)/(q1+qM))+q2); %-------Жесткость, [Нм]------%q23=(q2*q3)/(q2+q3);%-------Жесткость, [Нм]------%q34=(q3*q4)/(q3+q4);%-------Жесткость, [Нм]------%q15=(q1*q5)/(q1+q5);%-------Жесткость, [Нм]------%k12=1000;k23=1000;k34=1000;k15=1000;e2=0.0001;e3=0.0001;e4=0.0001;%-------Вязкость секций, []------%%-------Вязкость секций, []------%%-------Вязкость секций, []------%%-------Вязкость секций, []------%%-------Эксцентриситет (радиальные биения) ЦНД, [м]------%%-------Эксцентриситет (радиальные биения) ЦСД, [м]------%%-------Эксцентриситет (радиальные биения) ЦВД, [м]------%%----------------Коэффициенты трения при изгибных колебаниях-------------%c12=500000;c13=500000;c14=500000;c2=10000;c3=10000;c4=10000;m2=26980;m3=14000;m4=6466;%-----Коэффициент внешнего трения ЦНД ----------%%-----Коэффициент внешнего трения ЦСД ----------%%-----Коэффициент внешнего трения ЦВД ----------%%-----Коэффициент внутреннего трения ЦНД ----------%%-----Коэффициент внутреннего трения ЦСД ----------%%-----Коэффициент внутреннего трения ЦВД ----------%%-----Масса ЦНД, [кг] ----------%%-----Масса ЦСД, [кг] ----------%%-----Масса ЦВД, [кг] ----------%%----------------Определение ускорения свободного падения точно, [м/с^2]-------------%ff=61.25*pi/180; %------------- Географическая широта места расположения турбоагрегата,[рад]--------%hh=40;%------------- Высота над уровнем моря места расположения турбоагрегата,[м]--------%g=9.780327*(1+0.0053024*(sin(ff))^2-0.0000058*(sin(2*ff))^2)-3.086*hh/1000000; %------------ Ускорение свободного падения, [м/с^2]-------------%%----------------Расчет номинального режима работы генератора-------------%kH=sqrt(Lb*((xd-xdd)/w0));%-----Взаимная индуктивность обмоток ротора и статора----%Ld=xd_SI/w0;%----Индуктивность по оси d, [Гн]---------------%Lq=xq_SI/w0;%----Индуктивность по оси q, [Гн]---------------%Mt=P0/w0;%-----Номинальный момент турбины, [Нм]--------------%%----------Решение системы СЛУ---------------%A=[q1 -q1 0 0 0;-q1 q1+q2 -q2 0 0;0 -q2 q2+q3 -q3 0;0 0 -q3 q3+q4 -q4;0 0 0 -q4 q4+q5];311B=[Mt-kH*Iq0*Ib+(Ld-Lq)*Id0*Iq0;0;0;0;0];Q0=inv(A)*B;C=rref([A B]); %Приведение расширенной матрицы к треугольному видуQ0=C(1:5,6:6)%Выделение последнего столбца из матрицы% A=[q1 -q1 0 0 0;1 0 0 0 0;1 0 0 0 0;-q1 q1+q2 -q2 0 0;0 -q2 q2+q3 -q3 0;0 0 -q3 q3+q4 -q4;00 0 -q4 q4];%B=[Mt-kH*Iq0*Ib+(Ld-Lq)*Id0*Iq0;asin((ra*Id0-w0*Lq*Iq0)/(U0));acos((w0*Ld*Id0+ra*Iq0-kH*w0*Ib)/(-U0));0;0;0;0];%----Задаем начальные условия--------%Q10=Q0(1);Q20=Q0(2);Q30=Q0(3);Q40=Q0(4);Q50=Q0(5);dQ10=0;dQ11=1;dQ20=0;dQ21=1;dQ30=0;dQ31=1;dQ40=0;dQ41=1;dQ50=0;dQ51=1;dId10=1;dIq10=1;dIb10=1;Xx20=0;Xx21=0.0001;Xx30=0;Xx31=0.0001;Xx40=0;Xx41=0.0001;Yy20=0;Yy21=0.0001;Yy30=0;Yy31=0.0001;Yy40=0;Yy41=0.0001;dQ10=0;dQ11=Q0(1)/1000;dQ20=0;dQ21=Q0(2)/1000;312dQ30=0;dQ31=Q0(3)/1000;dQ40=0;dQ41=Q0(4)/1000;dQ50=0;dQ51=Q0(5)/1000;dId10=1;dIq10=1;dIb10=1;% Вектор решений TRS(dQ1;dQ2;dQ3;dQ4;dQ5;dId;dIq;dIb);%------------------Кратность форсировки по току-----------------%KFIb=1;%------------------Кратность форсировки по напряжению-----------------%KFUb=1;KL=2000;Ub1=0;Ub5=0;Ibb1=0;Ibb2=0;Fbb1=0;Fbb2=0;FFb2=0;FFb3=0;U0=KFUb*Ub; %----Номинальное напряжение ротора при соединении полюсов звездой,[В]---%U01=0;U05=0;Ub=KFUb*Ubb;%------------------Постоянная времени интегрирования в контуре регулированиянапряжения статора-----------------%%------------------Тиристорный преобразоатель напряжения статора-----------------%for n=1:KLdt=0.001 ;%------Задаем шаг интегрирования, [с]---------%t=dt*n;Fb=314+rand(1);if n>=200 dU(n)=100000*sin((100*t)).*exp(-5*(t));%------------------Регулирование по каналу напряжения генератора -----------------%t=dt*n+0.005;Uab(n)=(U0)*sin(Fb*t);Ubc(n)=(U0)*sin(Fb*t+2*pi/3);Uca(n)=(U0)*sin(Fb*t+4*pi/3);Uabc(n)=(abs(Uab(n))+abs(Ubc(n))+abs(Uca(n)))/2;UG=Uabc(n);U0G=(UG)*sin(Fb*t);dU01=((1/0.005)*dt*(U0G-U01))+U01;313U0=((1/0.003)*dt*(dU01-U05))+U05;U01=dU01;U05=U0;%------------------Регулирование по каналу напряжения возбудителя -----------------%t=dt*n;Ubab(n)=(Ub)*sin(Fb*t);Ubbc(n)=(Ub)*sin(Fb*t+2*pi/3);Ubca(n)=(Ub)*sin(Fb*t+4*pi/3);Ubabc(n)=(abs(Ubab(n))+abs(Ubbc(n))+abs(Ubca(n)))/2;UbbG=Ubabc(n);UbG=(UbbG)*sin(Fb*t);dUb1=((1/0.005)*dt*(UbG-Ub1))+Ub1;dUb=((1/0.003)*dt*(dUb1-Ub5))+Ub5;Ub1=dUb1;Ub5=dUb;%------------------Регулирование по каналу производной тока генератора -----------------%Ib1(n)=(1/0.1)*dt*(Ib-Ibb1)-Ibb1;Ibb1=Ib1(n);Ib2(n)=(1/0.10)*dt*(Ibb1-Ibb2)-Ibb2;Ibb2=Ib2(n);Ib=Ib+Ibb2;%------------------Регулирование по каналу производной частоты вращения генератора ----------------%Fb1(n)=(1/0.05)*dt*(Fb-Fbb1)-Fbb1;Fbb1=Fb1(n);Fb2(n)=(1/0.08)*dt*(Fbb1-Fbb2)-Fbb2;Fbb2=Fb2(n);Fb=Fb+Fbb2;%------------------Регулирование по каналу частоты вращения генератора -----------------%FFb3(n)=(1/2.5)*((2/5)*Fbb1-FFb3)*dt-FFb3;FFb3=FFb3(n);FFb(n)=(1/0.1)*dt*(FFb3-FFb2)-FFb3;FFb2=FFb(n);Fb=Fb+FFb2;Chast(n)=Fb;t=dt;%------------------Решение системы дифференциальных уравнений модели -----------------%K1_dQ1(n)=1/t^2+k12/(t*J1)+k15/(t*J1); %----Коэффициент при dQ1 в 1 уравнении------%K1_dQ2(n)=-k12/(t*J1);%----Коэффициент при dQ2 в 1 уравнении------%K1_dQ3(n)=0;%----Коэффициент при dQ3 в 1 уравнении------%K1_dQ4(n)=0;%----Коэффициент при dQ4 в 1 уравнении------%K1_dQ5(n)=-k15/(t*J1);%----Коэффициент при dQ5 в 1 уравнении------%K1_dId(n)=(xd_SI-xq_SI)*Iq0/w0;%----Коэффициент при dId в 1 уравнении------%K1_dIq(n)=kH*Ib-(xd_SI-xq_SI)*Id0/w0; %----Коэффициент при dIq в 1 уравнении------%K1_dIb(n)=kH*Iq0;%----Коэффициент при dIb в 1 уравнении------%K1_Xx2(n)=0;%----Коэффициент при Xx2 в 1 уравнении------%314K1_Xx3(n)=0;K1_Xx4(n)=0;K1_Yy2(n)=0;K1_Yy3(n)=0;K1_Yy4(n)=0;%----Коэффициент при Xx3 в 1 уравнении------%%----Коэффициент при Xx4 в 1 уравнении------%%----Коэффициент при Yy2 в 1 уравнении------%%----Коэффициент при Yy3 в 1 уравнении------%%----Коэффициент при Yy4 в 1 уравнении------%K2_dQ1(n)=-k12/(t*J2);%----Коэффициент при dQ1 во 2 уравнении------%K2_dQ2(n)=1/t^2+k12/(t*J2)+k23/(t*J2); %----Коэффициент при dQ2 во 2 уравнении-----%K2_dQ3(n)=-k23/(t*J2);%----Коэффициент при dQ3 во 2 уравнении------%K2_dQ4(n)=0;%----Коэффициент при dQ4 во 2 уравнении------%K2_dQ5(n)=0;%----Коэффициент при dQ5 во 2 уравнении------%K2_dId(n)=0;%----Коэффициент при dId во 2 уравнении------%K2_dIq(n)=0;%----Коэффициент при dIq во 2 уравнении------%K2_dIb(n)=0;%----Коэффициент при dIb во 2 уравнении------%K2_Xx2(n)=0;%----Коэффициент при Xx2 во 2 уравнении------%K2_Xx3(n)=0;%----Коэффициент при Xx3 во 2 уравнении------%K2_Xx4(n)=0;%----Коэффициент при Xx4 во 2 уравнении------%K2_Yy2(n)=0;%----Коэффициент при Yy2 во 2 уравнении------%K2_Yy3(n)=0;%----Коэффициент при Yy3 во 2 уравнении------%K2_Yy4(n)=0;%----Коэффициент при Yy4 во 2 уравнении------%K3_dQ1(n)=0;%----Коэффициент при dQ1 в 3 уравнении------%K3_dQ2(n)=-k23/(t*J3);%----Коэффициент при dQ2 в 3 уравнении------%K3_dQ3(n)=1/t^2+k23/(t*J3)+k34/(t*J3); %----Коэффициент при dQ3 в 3 уравнении------%K3_dQ4(n)=-k34/(t*J3);%----Коэффициент при dQ4 в 3 уравнении------%K3_dQ5(n)=0;%----Коэффициент при dQ5 в 3 уравнении------%K3_dId(n)=0;%----Коэффициент при dId в 3 уравнении------%K3_dIq(n)=0;%----Коэффициент при dIq в 3 уравнении------%K3_dIb(n)=0;%----Коэффициент при dIb в 3 уравнении------%K3_Xx2(n)=0;%----Коэффициент при Xx2 в 3 уравнении------%K3_Xx3(n)=0;%----Коэффициент при Xx3 в 3 уравнении------%K3_Xx4(n)=0;%----Коэффициент при Xx4 в 3 уравнении------%K3_Yy2(n)=0;%----Коэффициент при Yy2 в 3 уравнении------%K3_Yy3(n)=0;%----Коэффициент при Yy3 в 3 уравнении------%K3_Yy4(n)=0;%----Коэффициент при Yy4 в 3 уравнении------%K4_dQ1(n)=0;%----Коэффициент при dQ1 в 4 уравнении------%K4_dQ2(n)=0;%----Коэффициент при dQ2 в 4 уравнении------%K4_dQ3(n)=-k34/(t*J4);%----Коэффициент при dQ3 в 4 уравнении------%K4_dQ4(n)=1/(t^2)+k34/(t*J4);%----Коэффициент при dQ4 в 4 уравнении------%K4_dQ5(n)=0;%----Коэффициент при dQ5 в 4 уравнении------%K4_dId(n)=0;%----Коэффициент при dId в 4 уравнении------%K4_dIq(n)=0;%----Коэффициент при dIq в 4 уравнении------%K4_dIb(n)=0;%----Коэффициент при dIb в 4 уравнении------%K4_Xx2(n)=0;%----Коэффициент при Xx2 в 4 уравнении------%K4_Xx3(n)=0;%----Коэффициент при Xx3 в 4 уравнении------%K4_Xx4(n)=0;%----Коэффициент при Xx4 в 4 уравнении------%K4_Yy2(n)=0;%----Коэффициент при Yy2 в 4 уравнении------%K4_Yy3(n)=0;%----Коэффициент при Yy3 в 4 уравнении------%K4_Yy4(n)=0;%----Коэффициент при Yy4 в 4 уравнении------%315K5_dQ1(n)=-k15/(t*J5);%----Коэффициент при dQ1 в 5 уравнении------%K5_dQ2(n)=0;%----Коэффициент при dQ2 в 5 уравнении------%K5_dQ3(n)=0;%----Коэффициент при dQ3 в 5 уравнении------%K5_dQ4(n)=0;%----Коэффициент при dQ4 в 5 уравнении------%K5_dQ5(n)=1/(t^2)+k15/(t*J5);%----Коэффициент при dQ5 в 5 уравнении------%K5_dId(n)=0;%----Коэффициент при dId в 5 уравнении------%K5_dIq(n)=0;%----Коэффициент при dIq в 5 уравнении------%K5_dIb(n)=0;%----Коэффициент при dIb в 5 уравнении------%K5_Xx2(n)=0;%----Коэффициент при Xx2 в 5 уравнении------%K5_Xx3(n)=0;%----Коэффициент при Xx3 в 5 уравнении------%K5_Xx4(n)=0;%----Коэффициент при Xx4 в 5 уравнении------%K5_Yy2(n)=0;%----Коэффициент при Yy2 в 5 уравнении------%K5_Yy3(n)=0;%----Коэффициент при Yy3 в 5 уравнении------%K5_Yy4(n)=0;%----Коэффициент при Yy4 в 5 уравнении------%K6_dQ1(n)=-xq_SI*Iq0/(w0*t);%----Коэффициент при dQ1 в 6 уравнении------%K6_dQ2(n)=0;%----Коэффициент при dQ2 в 6 уравнении------%K6_dQ3(n)=0;%----Коэффициент при dQ3 в 6 уравнении------%K6_dQ4(n)=0;%----Коэффициент при dQ4 в 6 уравнении------%K6_dQ5(n)=0;%----Коэффициент при dQ5 в 6 уравнении------%K6_dId(n)=ra+xd_SI/(w0*t);%----Коэффициент при dId в 6 уравнении------%K6_dIq(n)=-xq_SI;%----Коэффициент при dIq в 6 уравнении------%K6_dIb(n)=-kH/t;%----Коэффициент при dIb в 6 уравнении------%K6_Xx2(n)=0;%----Коэффициент при Xx2 в 6 уравнении------%K6_Xx3(n)=0;%----Коэффициент при Xx3 в 6 уравнении------%K6_Xx4(n)=0;%----Коэффициент при Xx4 в 6 уравнении------%K6_Yy2(n)=0;%----Коэффициент при Yy2 в 6 уравнении------%K6_Yy3(n)=0;%----Коэффициент при Yy3 в 6 уравнении------%K6_Yy4(n)=0;%----Коэффициент при Yy4 в 6 уравнении------%K7_dQ1(n)=-kH*Ib/t+xd_SI*Id0/(w0*t); %----Коэффициент при dQ1 в 7 уравнении------%K7_dQ2(n)=0;%----Коэффициент при dQ2 в 7 уравнении------%K7_dQ3(n)=0;%----Коэффициент при dQ3 в 7 уравнении------%K7_dQ4(n)=0;%----Коэффициент при dQ4 в 7 уравнении------%K7_dQ5(n)=0;%----Коэффициент при dQ5 в 7 уравнении------%K7_dId(n)=xd_SI;%----Коэффициент при dId в 7 уравнении------%K7_dIq(n)=ra+xq_SI/(w0*t);%----Коэффициент при dIq в 7 уравнении------%K7_dIb(n)=-kH*w0;%----Коэффициент при dIb в 7 уравнении------%K7_Xx2(n)=0;%----Коэффициент при Xx2 в 7 уравнении------%K7_Xx3(n)=0;%----Коэффициент при Xx3 в 7 уравнении------%K7_Xx4(n)=0;%----Коэффициент при Xx4 в 7 уравнении------%K7_Yy2(n)=0;%----Коэффициент при Yy2 в 7 уравнении------%K7_Yy3(n)=0;%----Коэффициент при Yy3 в 7 уравнении------%K7_Yy4(n)=0;%----Коэффициент при Yy4 в 7 уравнении------%K8_dQ1(n)=0;K8_dQ2(n)=0;K8_dQ3(n)=0;K8_dQ4(n)=0;K8_dQ5(n)=0;K8_dId(n)=kH/t;K8_dIq(n)=0;%----Коэффициент при dQ1 в 8 уравнении------%%----Коэффициент при dQ2 в 8 уравнении------%%----Коэффициент при dQ3 в 8 уравнении------%%----Коэффициент при dQ4 в 8 уравнении------%%----Коэффициент при dQ5 в 8 уравнении------%%----Коэффициент при dId в 8 уравнении------%%----Коэффициент при dIq в 8 уравнении------%316K8_dIb(n)=rb+Lb/t;K8_Xx2(n)=0;K8_Xx3(n)=0;K8_Xx4(n)=0;K8_Yy2(n)=0;K8_Yy3(n)=0;K8_Yy4(n)=0;%----Коэффициент при dIb в 8 уравнении------%%----Коэффициент при Xx2 в 8 уравнении------%%----Коэффициент при Xx3 в 8 уравнении------%%----Коэффициент при Xx4 в 8 уравнении------%%----Коэффициент при Yy2 в 8 уравнении------%%----Коэффициент при Yy3 в 8 уравнении------%%----Коэффициент при Yy4 в 8 уравнении------%K9_dQ1(n)=0;%----Коэффициент при dQ1 в 9 уравнении------%K9_dQ2(n)=c2*Yy21/t;%----Коэффициент при dQ2 в 9 уравнении------%K9_dQ3(n)=0;%----Коэффициент при dQ3 в 9 уравнении------%K9_dQ4(n)=0;%----Коэффициент при dQ4 в 9 уравнении------%K9_dQ5(n)=0;%----Коэффициент при dQ5 в 9 уравнении------%K9_dId(n)=0;%----Коэффициент при dId в 9 уравнении------%K9_dIq(n)=0;%----Коэффициент при dIq в 9 уравнении------%K9_dIb(n)=0;%----Коэффициент при dIb в 9 уравнении------%K9_Xx2(n)=m2/t^2+c12/t+c2/t;%----Коэффициент при Xx2 в 9 уравнении------%K9_Xx3(n)=0;%----Коэффициент при Xx3 в 9 уравнении------%K9_Xx4(n)=0;%----Коэффициент при Xx4 в 9 уравнении------%K9_Yy2(n)=0;%----Коэффициент при Yy2 в 9 уравнении------%K9_Yy3(n)=0;%----Коэффициент при Yy3 в 9 уравнении------%K9_Yy4(n)=0;%----Коэффициент при Yy4 в 9 уравнении------%K10_dQ1(n)=0;%----Коэффициент при dQ1 в 10 уравнении------%K10_dQ2(n)=0;%----Коэффициент при dQ2 в 10 уравнении------%K10_dQ3(n)=c3*Yy31/t;%----Коэффициент при dQ3 в 10 уравнении------%K10_dQ4(n)=0;%----Коэффициент при dQ4 в 10 уравнении------%K10_dQ5(n)=0;%----Коэффициент при dQ5 в 10 уравнении------%K10_dId(n)=0;%----Коэффициент при dId в 10 уравнении------%K10_dIq(n)=0;%----Коэффициент при dIq в 10 уравнении------%K10_dIb(n)=0;%----Коэффициент при dIb в 10 уравнении------%K10_Xx2(n)=0;%----Коэффициент при Xx2 в 10 уравнении------%K10_Xx3(n)=m3/t^2+c13/t+c3/t;%----Коэффициент при Xx3 в 10 уравнении------%K10_Xx4(n)=0;%----Коэффициент при Xx4 в 10 уравнении------%K10_Yy2(n)=0;%----Коэффициент при Yy2 в 10 уравнении------%K10_Yy3(n)=0;%----Коэффициент при Yy3 в 10 уравнении------%K10_Yy4(n)=0;%----Коэффициент при Yy4 в 10 уравнении------%K11_dQ1(n)=0;%----Коэффициент при dQ1 в 11 уравнении------%K11_dQ2(n)=0;%----Коэффициент при dQ2 в 11 уравнении------%K11_dQ3(n)=0;%----Коэффициент при dQ3 в 11 уравнении------%K11_dQ4(n)=c4*Yy41/t;%----Коэффициент при dQ4 в 11 уравнении------%K11_dQ5(n)=0;%----Коэффициент при dQ5 в 11 уравнении------%K11_dId(n)=0;%----Коэффициент при dId в 11 уравнении------%K11_dIq(n)=0;%----Коэффициент при dIq в 11 уравнении------%K11_dIb(n)=0;%----Коэффициент при dIb в 11 уравнении------%K11_Xx2(n)=0;%----Коэффициент при Xx2 в 11 уравнении------%K11_Xx3(n)=0;%----Коэффициент при Xx3 в 11 уравнении------%K11_Xx4(n)=m4/t^2+c14/t+c4/t;%----Коэффициент при Xx4 в 11 уравнении------%K11_Yy2(n)=0;%----Коэффициент при Yy2 в 11 уравнении------%K11_Yy3(n)=0;%----Коэффициент при Yy3 в 11 уравнении------%K11_Yy4(n)=0;%----Коэффициент при Yy4 в 11 уравнении------%317K12_dQ1(n)=0;%----Коэффициент при dQ1 в 12 уравнении------%K12_dQ2(n)=-c2*Xx21/t;%----Коэффициент при dQ2 в 12 уравнении------%K12_dQ3(n)=0;%----Коэффициент при dQ3 в 12 уравнении------%K12_dQ4(n)=0;%----Коэффициент при dQ4 в 12 уравнении------%K12_dQ5(n)=0;%----Коэффициент при dQ5 в 12 уравнении------%K12_dId(n)=0;%----Коэффициент при dId в 12 уравнении------%K12_dIq(n)=0;%----Коэффициент при dIq в 12 уравнении------%K12_dIb(n)=0;%----Коэффициент при dIb в 12 уравнении------%K12_Xx2(n)=0;%----Коэффициент при Xx2 в 12 уравнении------%K12_Xx3(n)=0;%----Коэффициент при Xx3 в 12 уравнении------%K12_Xx4(n)=0;%----Коэффициент при Xx4 в 12 уравнении------%K12_Yy2(n)=m2/t^2+c12/t+c2/t;%----Коэффициент при Yy2 в 12 уравнении------%K12_Yy3(n)=0;%----Коэффициент при Yy3 в 12 уравнении------%K12_Yy4(n)=0;%----Коэффициент при Yy4 в 12 уравнении------%K13_dQ1(n)=0;%----Коэффициент при dQ1 в 13 уравнении------%K13_dQ2(n)=0;%----Коэффициент при dQ2 в 13 уравнении------%K13_dQ3(n)=-c3*Xx31/t;%----Коэффициент при dQ3 в 13 уравнении------%K13_dQ4(n)=0;%----Коэффициент при dQ4 в 13 уравнении------%K13_dQ5(n)=0;%----Коэффициент при dQ5 в 13 уравнении------%K13_dId(n)=0;%----Коэффициент при dId в 13 уравнении------%K13_dIq(n)=0;%----Коэффициент при dIq в 13 уравнении------%K13_dIb(n)=0;%----Коэффициент при dIb в 13 уравнении------%K13_Xx2(n)=0;%----Коэффициент при Xx2 в 13 уравнении------%K13_Xx3(n)=0;%----Коэффициент при Xx3 в 13 уравнении------%K13_Xx4(n)=0;%----Коэффициент при Xx4 в 13 уравнении------%K13_Yy2(n)=0;%----Коэффициент при Yy2 в 13 уравнении------%K13_Yy3(n)=m3/t^2+c13/t+c3/t;%----Коэффициент при Yy3 в 13 уравнении------%K13_Yy4(n)=0;%----Коэффициент при Yy4 в 13 уравнении------%K14_dQ1(n)=0;%----Коэффициент при dQ1 в 14 уравнении------%K14_dQ2(n)=0;%----Коэффициент при dQ2 в 14 уравнении------%K14_dQ3(n)=0;%----Коэффициент при dQ3 в 14 уравнении------%K14_dQ4(n)=-c4*Xx41/t;%----Коэффициент при dQ4 в 14 уравнении------%K14_dQ5(n)=0;%----Коэффициент при dQ5 в 14 уравнении------%K14_dId(n)=0;%----Коэффициент при dId в 14 уравнении------%K14_dIq(n)=0;%----Коэффициент при dIq в 14 уравнении------%K14_dIb(n)=0;%----Коэффициент при dIb в 14 уравнении------%K14_Xx2(n)=0;%----Коэффициент при Xx2 в 14 уравнении------%K14_Xx3(n)=0;%----Коэффициент при Xx3 в 14 уравнении------%K14_Xx4(n)=0;%----Коэффициент при Xx4 в 14 уравнении------%K14_Yy2(n)=0;%----Коэффициент при Yy2 в 14 уравнении------%K14_Yy3(n)=0;%----Коэффициент при Yy3 в 14 уравнении------%K14_Yy4(n)=m4/t^2+c14/t+c4/t;%----Коэффициент при Yy4 в 14 уравнении------%P1(n)=Mt-((-2*dQ11+dQ10)/t^2+(k12/(t*J1))*(-dQ11+dQ21)+(k15/(t*J1))*(dQ11+dQ51)+(q12/J1)*(dQ11-dQ21)+(q15/J1)*(dQ11-dQ51)); %----Правая часть уравнения1-------%P2(n)=-((-2*dQ21+dQ20)/t^2+(k12/(t*J2))*(-dQ21+dQ11)+(k23/(t*J2))*(dQ21+dQ31)+(q12/J2)*(dQ21-dQ11)+(q23/J2)*(dQ21-dQ31)); %----Правая часть уравнения2-------%318P3(n)=-((-2*dQ31+dQ30)/t^2+(k23/(t*J3))*(-dQ31+dQ21)+(k34/(t*J3))*(dQ31+dQ41)+(q23/J3)*(dQ31-dQ21)+(q34/J3)*(dQ31-dQ41)); %----Правая часть уравнения3-------%P4(n)=-((-2*dQ41+dQ40)/t^2+(k34/(t*J4))*(-dQ41+dQ31)+(q34/J4)*(dQ41-dQ31));%---Правая часть уравнения 4-------%P5(n)=-((-2*dQ51+dQ50)/t^2+(k15/(t*J5))*(-dQ51+dQ11)+(q15/J5)*(dQ51-dQ11));%---Правая часть уравнения 5-------%P6(n)=-((xq_SI*Iq0/(w0*t))*dQ11+kH*dIb10/t+U0*dQ11*cos(Q10)+dU(n)*sin(Q10)); %---Правая часть уравнения 6-------%P7(n)=-((-xd_SI*Id0/(w0*t))*dQ11+kH*Ib*dQ11/t-xq_SI*dIq10/(w0*t)U0*dQ11*sin(Q10)+dU(n)*cos(Q10)); %----Правая часть уравнения 7-------%P8(n)=-(-Lb*dIb10/t-kH*dId10/t-dUb); %----Правая часть уравнения 8-------%P9(n)=q2*e2*cos(dQ21)-m2*(-2*Xx21+Xx20)/t^2-c12*(-Xx21/t)-q2*Xx21-c2*(-Xx21/tdQ21*Yy21/t); %----Правая часть уравнения 9 -------%P10(n)=q3*e3*cos(dQ31)-m3*(-2*Xx31+Xx30)/t^2-c13*(-Xx31/t)-q3*Xx31-c3*(-Xx31/tdQ31*Yy31/t); %----Правая часть уравнения 10 -------%P11(n)=q4*e4*cos(dQ41)-m4*(-2*Xx41+Xx40)/t^2-c14*(-Xx41/t)-q4*Xx41-c4*(-Xx41/tdQ41*Yy41/t); %----Правая часть уравнения 11 -------%P12(n)=q2*e2*sin(dQ21)-m2*g-m2*(-2*Yy21+Yy20)/t^2-c12*(-Yy21/t)-q2*Yy21-c2*(Yy21/t+dQ21*Xx21/t); %----Правая часть уравнения 12 -------%P13(n)=q3*e3*sin(dQ31)-m3*g-m3*(-2*Yy31+Yy30)/t^2-c13*(-Yy31/t)-q3*Yy31-c3*(Yy31/t+dQ31*Xx31/t); %----Правая часть уравнения 13 -------%P14(n)=q4*e4*sin(dQ41)-m4*g-m4*(-2*Yy41+Yy40)/t^2-c14*(-Yy41/t)-q4*Yy41-c4*(Yy41/t+dQ41*Xx41/t); %----Правая часть уравнения 14 -------%AA=[K1_dQ1(n) K1_dQ2(n) K1_dQ3(n) K1_dQ4(n) K1_dQ5(n) K1_dId(n) K1_dIq(n)K1_dIb(n) K1_Xx2(n) K1_Xx3(n) K1_Xx4(n) K1_Yy2(n) K1_Yy3(n) K1_Yy4(n);K2_dQ1(n) K2_dQ2(n) K2_dQ3(n) K2_dQ4(n) K2_dQ5(n) K2_dId(n) K2_dIq(n) K2_dIb(n)K2_Xx2(n) K2_Xx3(n) K2_Xx4(n) K2_Yy2(n) K2_Yy3(n) K2_Yy4(n);K3_dQ1(n) K3_dQ2(n) K3_dQ3(n) K3_dQ4(n) K3_dQ5(n) K3_dId(n) K3_dIq(n) K3_dIb(n)K3_Xx2(n) K3_Xx3(n) K3_Xx4(n) K3_Yy2(n) K3_Yy3(n) K3_Yy4(n);K4_dQ1(n) K4_dQ2(n) K4_dQ3(n) K4_dQ4(n) K4_dQ5(n) K4_dId(n) K4_dIq(n) K4_dIb(n)K4_Xx2(n) K4_Xx3(n) K4_Xx4(n) K4_Yy2(n) K4_Yy3(n) K4_Yy4(n);K5_dQ1(n) K5_dQ2(n) K5_dQ3(n) K5_dQ4(n) K5_dQ5(n) K5_dId(n) K5_dIq(n) K5_dIb(n)K5_Xx2(n) K5_Xx3(n) K5_Xx4(n) K5_Yy2(n) K5_Yy3(n) K5_Yy4(n);K6_dQ1(n) K6_dQ2(n) K6_dQ3(n) K6_dQ4(n) K6_dQ5(n) K6_dId(n) K6_dIq(n) K6_dIb(n)K6_Xx2(n) K6_Xx3(n) K6_Xx4(n) K6_Yy2(n) K6_Yy3(n) K6_Yy4(n);K7_dQ1(n) K7_dQ2(n) K7_dQ3(n) K7_dQ4(n) K7_dQ5(n) K7_dId(n) K7_dIq(n) K7_dIb(n)K7_Xx2(n) K7_Xx3(n) K7_Xx4(n) K7_Yy2(n) K7_Yy3(n) K7_Yy4(n);K8_dQ1(n) K8_dQ2(n) K8_dQ3(n) K8_dQ4(n) K8_dQ5(n) K8_dId(n) K8_dIq(n) K8_dIb(n)K8_Xx2(n) K8_Xx3(n) K8_Xx4(n) K8_Yy2(n) K8_Yy3(n) K8_Yy4(n);K9_dQ1(n) K9_dQ2(n) K9_dQ3(n) K9_dQ4(n) K9_dQ5(n) K9_dId(n) K9_dIq(n) K9_dIb(n)K9_Xx2(n) K9_Xx3(n) K9_Xx4(n) K9_Yy2(n) K9_Yy3(n) K9_Yy4(n);K10_dQ1(n) K10_dQ2(n) K10_dQ3(n) K10_dQ4(n) K10_dQ5(n) K10_dId(n) K10_dIq(n)K10_dIb(n) K10_Xx2(n) K10_Xx3(n) K10_Xx4(n) K10_Yy2(n) K10_Yy3(n) K10_Yy4(n);K11_dQ1(n) K11_dQ2(n) K11_dQ3(n) K11_dQ4(n) K11_dQ5(n) K11_dId(n) K11_dIq(n)K11_dIb(n) K11_Xx2(n) K11_Xx3(n) K11_Xx4(n) K11_Yy2(n) K11_Yy3(n) K11_Yy4(n);K12_dQ1(n) K12_dQ2(n) K12_dQ3(n) K12_dQ4(n) K12_dQ5(n) K12_dId(n) K12_dIq(n)K12_dIb(n) K12_Xx2(n) K12_Xx3(n) K12_Xx4(n) K12_Yy2(n) K12_Yy3(n) K12_Yy4(n);319K13_dQ1(n) K13_dQ2(n) K13_dQ3(n) K13_dQ4(n) K13_dQ5(n) K13_dId(n) K13_dIq(n)K13_dIb(n) K13_Xx2(n) K13_Xx3(n) K13_Xx4(n) K13_Yy2(n) K13_Yy3(n) K13_Yy4(n);K14_dQ1(n) K14_dQ2(n) K14_dQ3(n) K14_dQ4(n) K14_dQ5(n) K14_dId(n) K14_dIq(n)K14_dIb(n) K14_Xx2(n) K14_Xx3(n) K14_Xx4(n) K14_Yy2(n) K14_Yy3(n) K14_Yy4(n)];BB=[P1(n);P2(n);P3(n);P4(n);P5(n);P6(n);P7(n);P8(n);P9(n);P10(n);P11(n);P12(n);P13(n);P14(n)];TRS=inv(AA)*BB;dQ1(n)=TRS(1);dQ2(n)=TRS(2);dQ3(n)=TRS(3);dQ4(n)=TRS(4);dQ5(n)=TRS(5);dId(n)=TRS(6);dIq(n)=TRS(7);dIb(n)=TRS(8);Xx2(n)=TRS(9);Xx3(n)=TRS(10);Xx4(n)=TRS(11);Yy2(n)=TRS(12);Yy3(n)=TRS(13);Yy4(n)=TRS(14);dQ10=dQ11;dQ11=TRS(1);dQ20=dQ21;dQ21=TRS(2);dQ30=dQ31;dQ31=TRS(3);dQ40=dQ41;dQ41=TRS(4);dQ50=dQ51;dQ51=TRS(5);dId10=TRS(6);dIq10=TRS(7);dIb10=TRS(8);Xx20=Xx21;Xx21=TRS(9);Xx30=Xx31;Xx31=TRS(10);Xx40=Xx41;Xx41=TRS(11);Yy20=Yy21;Yy21=TRS(12);Yy30=Yy31;Yy31=TRS(13);Yy40=Yy41;Yy41=TRS(14);end;320end;for n=1:KL;x(n)=n*dt;end;figureplot(x,dU);title('Напряжение на обмотках статора');ylabel('Напряжение, В');xlabel('Время, сек');grid on;figureplot(x,(dQ1)*180*60/pi);title('Изменение угла качания генератора ');xlabel('Время, сек');ylabel('Амплитуда угла качания, угл.
![](https://s.studizba.com/z.php?f=/templates/si/i/noavatar.png&w=100&h=100&t=1"){kind=avatar}
