Подать статью
Стать рецензентом
Том 247
Страницы:
122-131
Скачать том:

Исследование токов привода подъема мостовых кранов металлургических предприятий для раннего диагностирования превышения массы груза

Авторы:
И. Ю. Семыкина1
М. В. Кипервассер2
А. В. Герасимук3
Об авторах
Дата отправки:
2020-05-16
Дата принятия:
2021-03-02
Дата публикации:
2021-04-21

Аннотация

В статье рассматривается основанный на анализе токов приводного электродвигателя подход к созданию дополнительного средства защиты от аварийных ситуаций при эксплуатации мостовых кранов, связанных с подъемом груза массой, превышающей допустимую. Описана математическая модель привода мостового крана, а также результаты компьютерного моделирования. Показано, что в процессе подъема до начала этапа отрыва груза ток статора приводного электродвигателя от массы груза не зависит, но при отрыве груза уже за единицы периодов сетевого напряжения после натяжения каната при превышении массы гру за фиксируется поддающееся измерению превышение амплитудного значения тока. Данная закономерность подтверждена для ряда кранов различной грузоподъемности, применяемых на металлургических предприятиях. Продемонстрирована возможность диагностирования превышения массы поднимаемого груза с более высоким быстродействием, чем у существующих механических способов контроля перегрузки, при этом не требуется вносить изменения в конструктивные элементы мостовых кранов.

Ключевые слова:
привод мостового крана трехмассовая механическая система асинхронный двигатель динамическая модель в дискретной форме этап отрыва груза амплитуда тока статора
10.31897/PMI.2021.1.13
Перейти к тому 247

Введение. В число основных видов технологического и ремонтного оборудования на металлургических предприятиях входят грузоподъемные механизмы: мостовые краны и кран-балки различного назначения. Удельный вес грузоподъемных механизмов в составе технологического и ремонтного оборудования цехов металлургических предприятий может достигать 10-13 %. Аварии, связанные с использованием грузоподъемной техники, зачастую приводят к травматизму людей, оказавшихся в опасной зоне [3]. Кроме того, от надежной работы грузоподъемных механизмов зависит ритмичность технологических процессов и высокие экономические показатели предприятий. Поэтому обеспечение грузоподъемных механизмов современными быстродействующими средствами защиты от аварийных ситуаций и средствами диагностики является актуальной задачей [5, 14, 15].

Для мостовых кранов металлургических предприятий стандартными являются следующие аварийные ситуации: переподъем грузозахватного органа, нарушение целостности и геометрии элементов конструкции крана, обрыв подъемных канатов, перегрев и выход из строя электродвигателя подъема. К перечисленным аварийным ситуациям может привести увеличение нагрузки крана при подъеме груза недопустимой массы.

Для защиты от аварий, вызванных подъемом грузов недопустимой массы, используются ограничители грузоподъемности. В настоящее время применяются несколько способов контроля перегрузки, среди которых устройства с применением пружинно-рычажных систем [6] и устройства с использованием тензодатчиков [8, 9], которые, однако, не обладают достаточными точностью, чувствительностью и быстродействием. Также разработаны схемы защиты с применением дополнительных силовых элементов, однако они усложняют конструкцию грузоподъемных механизмов, снижая этим надежность защиты [7].

С учетом изложенного, а также в условиях сокращения персонала ремонтных и эксплуатирующих подразделений, для распознавания аварийного превышения массы груза необходимо использовать дополнительные возможности, в том числе методы косвенного контроля. Для косвенного контроля оборудования горно-металлургических предприятий используются методы, основанные на анализе изменения во времени параметров рабочего цикла [1], в том числе анализе динамического состояния приводных электродвигателей [11, 12].

Существуют методы косвенного контроля массы груза для управления и защиты кранового оборудования в процессе подъема для кранов иного назначения. Так в источниках [21, 23, 32] предлагаются системы мониторинга и диагностики для мобильного, башенного и мостового кранов соответственно, предоставляющие дополнительную информацию о работе крана его оператору. Однако все эти системы для определения массы груза предполагают использование в конструкции крана дополнительных датчиков – энкодеров, тензодатчиков. В статье [24] описана система управления для портового крана со встроенным контуром оценивания массы груза, для функционирования которого используются дополнительные датчики смещения груза относительно подъемного блока стрелы, однако в ней не предполагается возможность наличия слабины каната. В статье [28] показана система управления для подъемного крана буровой платформы, в которой масса груза вычисляется по данным тензодатчика подъемного каната и датчика гидродинамического давления. Работы [30, 31] рассматривают систему защиты от перегрузки при подъеме для кранов морского назначения в составе системы управления с упреждением, для функционирования которой используются совместно динамическая модель серводвигателя тяговой лебедки и подъемного механизма с грузом, однако отличия конструкции не позволяют переложить такой подход на мостовые краны металлургических предприятий.

Системы защиты от превышения массы груза с применением методов косвенного контроля для крановых механизмов применяются активно, однако эффективный подход для мостовых кранов, не требующий установки в механической подсистеме крана дополнительных датчиков, а использующий информацию о динамическом состоянии электродвигателей, в настоящее время отсутствует.

Постановка проблемы. Для создания дополнительного средства предотвращения аварийных ситуаций – подъема грузов недопустимой массы – с использованием контроля токов электродвигателя механизма подъема мостового крана необходимо проводить анализ его динамического состояния, а при анализе динамических режимов электроприводов горно-металлургического оборудования необходимо обязательно учитывать динамику механической подсистемы [4,]. Для механизмов подъема мостовых кранов разработан ряд математических моделей, описывающих динамические режимы и процессы в узлах и элементах оборудования [2, 19, 22, 27]. Однако эти модели очень сложны и поэтому не подходят для решения поставленной задачи, предполагая в последующем необходимость технической реализации в виде программно-аппаратного средства, для чего нужна относительно простая модель с сосредоточенными параметрами, но адекватно описывающая технологическую операцию подъема груза краном и позволяющая достоверно идентифицировать аварийную ситуацию превышения массы груза.

Операция подъема груза состоит из трех последовательных этапов, а именно выбора слабины каната до появления в нем усилия натяжения, натяжения каната до уравновешивания усилием в канате веса груза и отрыва груза с его последующим подъемом. Важной особенностью этапа выбора слабины каната является фактическое отсутствие взаимной связи между грузом и приводным электродвигателем до появления в канате усилия натяжения. Это значит, что диагностировать аварийное превышение массы груза косвенным методом возможно только в результате детального анализа математической модели электродвигателя привода подъема и кранового механизма при переходе от первого этапа ко второму и третьему, в ходе которого необходимо выявить диагностические признаки режима подъема груза чрезмерной массы.

Методология. Для математического описания механизма привода подъема мостовых кранов наиболее распространенной на практике является модель трехмассовой механической системы с упругими связями между массами [25], представленная на рис.1, где m1, m2, m3 – приведенные к канату массы вращающихся частей механизма подъема, груза и моста; x1, x2, x3 – перемещения масс; C1, C2 – жесткости каната и моста крана; l – слабина каната; P – пусковая сила электродвигателя, определяемая как произведение электромагнитного момента на валу двигателя M на коэффициент приведения Kp = i/Rбар; F – усилие в канате.

Рис.1. Трехмассовая механическая система механизма подъема

На основании данной модели составляются системы дифференциальных уравнений, описывающие систему на этапах выбора слабины каната, его натяжения и отрыва груза от опоры [25]. Последовательное моделирование каждого из этапов с учетом начальных условий при переходе от одного к другому позволяет адекватно описать динамическое состояние механизма привода подъема мостовых кранов, а условия начала и окончания этапов приведены в таблице.

Условия начала и окончания этапов движения механизма подъема

Этап Начальные условия Условия окончания
1. Выбор слабины каната $$ \begin{array} [c]{ll} x_1=0; \ x_2=0; \ x_3=0; \\ x_1^{'}=0; x_2^{'}=0; x_3^{'}=0; t=0; \end{array} $$ $$ \begin{array} \biggl\ \begin{cases} x_{1}=l; \\ F=C_1(x_1-l)>0; \end{cases} \\ x_{11k}=x_1; \ x_{21k}=0; \ x_{31k}=x_3; \\ x_{11k}^{'}=x_1^{'}; \ x_{21k}^{'}=0; \ x_{31k}^{'}=x_3^{'}; \ t=t_{1k}.\end{array} $$
2. Натяжение каната до уравновешивания $$ \begin{array} [c]{ll} x_1=x_{11k}; \ x_2=x_{21k}; \ x_3=x_{31k}; \\ x_1^{'}=x_{11k}^{'}; x_2^{'}=x_{21k}^{'}; x_3^{'}=x_{31k}^{'}; t=t_{1k}; \end{array} $$ $$ \begin{array} [c]{ll} F=C_1(x_1-x_2)=m_2g; \\ x_{12k}=x_1; \ x_{22k}=0; \ x_{32k}=x_3; \\ x_{12k}^{'}=x_1^{'}; \ x_{22k}^{'}=0; \ x_{32k}^{'}=x_3^{'}; \ t=t_{2k}.\end{array} $$
3. Отрыв груза с последующим подъемом $$ \begin{array} [c]{ll} x_1=x_{12k}; \ x_2=x_{22k}; \ x_3=x_{32k}; \\ x_1^{'}=x_{12k}^{'}; x_2^{'}=x_{22k}^{'}; x_3^{'}=x_{32k}^{'}; t=t_{2k}; \end{array} $$ -

Движение масс трехмассовой механической системы после отрыва груза описывается следующей системой уравнений:

$$ \begin{gathered}\left\{ \begin{array}{ll} m_1\frac{d^2x_1}{dt}=k_pM-C_1[x_1-x_2]; \\ m_2\frac{d^2x_2}{dt}=C_1[x_1-x_2]-m_2g+m_2 \frac {d^2x_3}{dt}; \\ m_3\frac{d^2x_3}{dt}=m_3g+C_1[x_1-x_2]-C_2x_3. \end{array} \right.\end{gathered} \tag*{(1)}$$

Для последующего анализа систему уравнений удобно привести к линейному виду, выполнив замену переменных:

$$ \begin{gathered}\left\{ \begin{array}{ll} \frac{dw_1}{dt}=w_4; \ \frac{dw_4}{dt}=\frac{1}{m_1} \big [k_pM-C_1(w_1-w_2)\big]; \\ \frac{dw_2}{dt}=w_5; \ \frac{dw_5}{dt}=\frac{1}{m_2} \big [C_1(w_1-w_2)-m_1g\big]+\frac{dw_6}{dt}; \\ \frac{dw_3}{dt}=w_6; \ \frac{dw_6}{dt}=g+\frac{1}{m_3} \big [C_1(w_1-w_2)-C_2w_3\big]. \end{array} \right.\end{gathered} \tag*{(2)}$$

Данная система уравнений показывает, что электромагнитный момент двигателя M и масса груза m2 при переходе к этапам натяжения каната и отрыва груза оказываются динамически связанными, а значит анализ величин, влияющих на электромагнитный момент на этом этапе, способен выявить диагностические признаки режима подъема груза чрезмерной массы.

Систему уравнений (2) можно решить различными методами. Поскольку современные системы диагностики технически реализуются на базе микроконтроллерных устройств, для которых характерно квантование сигналов по времени, удобно привести систему уравнений (2) к дискретному виду. Для выполнения дискретизации используются последовательно преобразование Лапласа и билинейное преобразование (преобразование Тастина).

После перехода от пространства оригиналов к пространству комплексной переменной p, выполняя преобразование Лапласа, и далее от нее через билинейное преобразование к переменной z, получим следующую дискретную модель:

$$ \begin{cases} w_1(k)-w_1(k-1)=(w_4(k)+w_4(k-1))T/2; \\ w_2(k)-w_2(k-1)=(w_5(k)+w_5(k-1))T/2; \\ w_3(k)-w_3(k-1)=(w_6(k)+w_6(k-1))T/2; \\ w_4(k)-w_4(k-1)=\frac{T}{2m_1}\left( \begin{array}{ll}k_pM(k)+kM(k-1)-\\ -C_1\left(\begin{array}{ll}w_1(k)+w_1(k-1)-\\ -w_2(k)-w_2(k-1)\end{array} \right)\end{array}\right);\\ w_5(k)-w_5(k-1)=\frac{T}{2m_2} \begin{array}{ll} -\left( C_1\Biggl(\begin{array}{ll}w_1(k)+w_1(k-1)-\\ -w_2(k)-w_2(k-1)\end{array}\right) \end{array}-m_1g\Biggr)+\frac{2}{T}(w_6(k)-w_6(k-1));\\ w_6(k)-w_6(k-1)=\frac{T}{2}g+\frac{T}{2m_3}\left(C_1\left(\begin{array}{ll}w_1(k)+w_1(k-1)-\\ -w_2(k)-w_2(k-1)\end{array}\right)-C_2w_3(1+z^{-1})\right). \end{cases} \qquad\tag*{(3)} $$

Система уравнений (3) должна решаться совместно с математической моделью электропривода, включая мостовой кран в состав электромеханического комплекса [13]. В качестве электропривода рассматривается привод на базе асинхронного двигателя.

Для составления математического описания асинхронного двигателя может быть использовано несколько подходов [16, 17, 26, 29]. Одним из них является моделирование в фазных координатах, которое наиболее соответствует задачам диагностирования, поскольку является более подходящим в случае неисправностей самого двигателя.

Асинхронный двигатель в фазных координатах a, b, c описывается системой дифференциальных уравнений [20]:

$$ \begin{gathered}\left\{ \begin{array}{ll} \frac{d\Psi_A}{dt}+r_Ai_A=U_A; \ \frac{d\Psi_a}{dt}+r_ai_a=0; \ W_\mathrm{ЭМ}=\frac{1}{2}(\Psi_Ai_A+\Psi_Bi_B+\Psi_Ci_C+\Psi_ai_a+\Psi_bi_b+\Psi_ci_c); \\ \frac{d\Psi_B}{dt}+r_Bi_B=U_B; \ \frac{d\Psi_b}{dt}+r_bi_b=0; \ M=p\frac{\delta W_\mathrm{ЭМ}}{\delta \gamma} \\ \frac{d\Psi_C}{dt}+r_Ci_C=U_C; \ \frac{d\Psi_c}{dt}+r_ci_c=0; \ \frac{d \gamma}{dt}=\omega_r, \end{array} \right.\end{gathered} \tag*{(4)}$$

где UA, UB, UC – напряжения на зажимах обмотки статора; ΨA, ΨB, ΨC, iA, iB, iC – потокосцепления и токи обмоток фаз статора; rA, rB, rC – активные сопротивления обмоток фаз статора; Ψa, Ψb, Ψc, ia, ib, ic – потокосцепления и токи обмоток фаз ротора; ra, rb, rc – активные сопротивления обмоток фаз ротора; γ – угловое положение ротора; WЭМ – запас электромагнитной энергии машины; p – число пар полюсов; ωr – частота вращения ротора; ωs – синхронная частота вращения поля статора.

Потокосцепления статорных и роторных обмоток описываются следующими уравнениями:

$$ \begin{gathered}\left\{ \begin{array}{ll} \Psi_{A}=L_Ai_A+M_{AB}i_B+M_{AC}i_C+M_{Aa}i_a+M_{Ab}i_b+M_{Ac}i_c; \\ \Psi_{B}=L_Bi_B+M_{BC}i_C+M_{AB}i_A+M_{Ba}i_a+M_{Bb}i_b+M_{Bc}i_c;\\ \Psi_{C}=L_Ci_C+M_{AC}i_C+M_{BC}i_B+M_{Ca}i_a+M_{Cb}i_b+M_{Cc}i_c; \\ \Psi_{a}=L_ai_a+M_{ab}i_b+M_{ac}i_c+M_{aA}i_A+M_{aB}i_B+M_{aC}i_C;\\ \Psi_{b}=L_bi_b+M_{bc}i_c+M_{ab}i_a+M_{bA}i_A+M_{bB}i_B+M_{bC}i_C;\\ \Psi_{c}=L_ci_c+M_{ac}i_c+M_{bc}i_b+M_{cA}i_A+M_{cB}i_B+M_{cC}i_C, \end{array} \right.\end{gathered} \tag*{(5)}$$

где LA, LB, LC – индуктивности фаз статора; La, Lb, Lc – индуктивности фаз ротора; MAB, MBC, MAC – взаимные индуктивности между обмотками статора; Mab, Mbc, Mac – взаимные индуктивности между обмотками ротора.

Поскольку обмотки статора и ротора перемещаются относительно друг друга с частотой ωr, взаимные индуктивности между обмотками определяются гармоническими зависимостями [20], где MAam, MBbm, MCcm – максимальные значения взаимной индуктивности при прохождении обмотками ротора фаз a, b и с соответственно точно под обмотками статора фаз А, В и С. Например, для фазы А:

$$ \begin{gathered} M_{Aa}=M_{aA}=M_{Aam} \cos (\omega_rt);\\ M_{Ab}=M_{bA}=M_{Aam} \cos (\omega_rt+\frac{2\pi}{3}); \\ M_{Ac}=M_{cA}=M_{Aam} \cos (\omega_rt-\frac{2\pi}{3}). \end{gathered} \tag*{(6)}$$

В результате дифференциальная трехфазная модель асинхронного электродвигателя в естественных фазных координатах, полученная из системы уравнений (4), выглядит следующим образом:

$$ \begin{gathered}\left\{ \begin{array}{ll} i_A=\frac{1}{r_A}\Biggl(U_A-\Biggl(\frac{\delta \Psi_A}{\delta \gamma}\omega_r+L_A\frac{di_A}{dt}+M_{AB}\frac{di_B}{dt}+M_{AC}\frac{di_C}{dt}+M_{Aa}\frac{di_a}{dt}+M_{Ab}\frac{di_b}{dt}+M_{Ac}\frac{di_c}{dt}\Biggl)\Biggl); \\ i_B=\frac{1}{r_B}\Biggl(U_B-\Biggl(\frac{\delta \Psi_B}{\delta \gamma}\omega_r+M_{AB}\frac{di_A}{dt}+L_{B}\frac{di_B}{dt}+M_{BC}\frac{di_C}{dt}+M_{Ba}\frac{di_a}{dt}+M_{Bb}\frac{di_b}{dt}+M_{Bc}\frac{di_c}{dt}\Biggl)\Biggl); \\ i_C=\frac{1}{r_C}\Biggl(U_C-\Biggl(\frac{\delta \Psi_C}{\delta \gamma}\omega_r+M_{AC}\frac{di_A}{dt}+L_{C}\frac{di_C}{dt}+M_{BC}\frac{di_B}{dt}+M_{Ca}\frac{di_a}{dt}+M_{Cb}\frac{di_b}{dt}+M_{Cc}\frac{di_c}{dt}\Biggl)\Biggl); \\ i_a=-\frac{1}{r_a}\Biggl(\frac{\delta \Psi_a}{\delta \gamma}\omega_r+L_{a}\frac{di_a}{dt}+M_{ab}\frac{di_B}{dt}+M_{ac}\frac{di_c}{dt}+L_{a}\frac{di_A}{dt}+M_{aB}\frac{di_B}{dt}+M_{aC}\frac{di_B}{dt}\Biggl); \\ i_b=-\frac{1}{r_b}\Biggl(\frac{\delta \Psi_b}{\delta \gamma}\omega_r+L_{b}\frac{di_b}{dt}+M_{ab}\frac{di_a}{dt}+M_{bc}\frac{di_c}{dt}+M_{ba}\frac{di_A}{dt}+M_{bB}\frac{di_B}{dt}+M_{bC}\frac{di_C}{dt}\Biggl); \\ i_b=-\frac{1}{r_b}\Biggl(\frac{\delta \Psi_b}{\delta \gamma}\omega_r+L_{b}\frac{di_b}{dt}+M_{ab}\frac{di_a}{dt}+M_{bc}\frac{di_c}{dt}+M_{ba}\frac{di_A}{dt}+M_{bB}\frac{di_B}{dt}+M_{bC}\frac{di_C}{dt}\Biggl); \\ i_c=-\frac{1}{r_c}\Biggl(\frac{\delta \Psi_c}{\delta \gamma}\omega_r+L_{c}\frac{di_c}{dt}+M_{ac}\frac{di_a}{dt}+M_{bc}\frac{di_B}{dt}+M_{cA}\frac{di_A}{dt}+M_{cB}\frac{di_B}{dt}+M_{cC}\frac{di_C}{dt}\Biggl); \\ M=\frac{p}{2}\Biggl(\frac{\delta \Psi_A}{\delta \gamma}i_A+\frac{\delta \Psi_B}{\delta \gamma}i_B+\frac{\delta \Psi_C}{\delta \gamma}i_C+\frac{\delta \Psi_a}{\delta \gamma}i_a+\frac{\delta \Psi_b}{\delta \gamma}i_b+\frac{\delta \Psi_c}{\delta \gamma}i_c\Biggl);\ \frac{d \gamma}{dt}=\omega_r, \\ \end{array} \right.\end{gathered} \qquad\tag*{(7)}$$

где производные зависимых от угла γ потокосцеплений определяются путем дифференцирования системы уравнений (5) по времени как сложной функции с учетом гармонических зависимостей взаимной индуктивности (6) для фазы А и аналогичных им для фаз В и С.

Систему уравнений (7), аналогично модели трехмассовой механической системы, посредством билинейного преобразования [10] удобно представить в дискретной форме:

$$ i(k)=[A(k)+C(k)]^{-1}[U(k)+(-1)U(k-1)]+[A(k)+C(k)]^{-1}[(-1)B(k-1)+D(k-1)]i(k-1), \qquad\tag*{(8)}$$

где матрицы напряжений и токов:

$$ U(k)={\begin{bmatrix} U_A(k)\\ U_B(k)\\ U_C(k) \\ U_a(k)\\ U_b(k) \\ U_c(k) \end{bmatrix}}; \ U(k-1)={\begin{bmatrix} U_A(k-1)\\ U_B(k-1)\\ U_C(k-1) \\ U_a(k-1)\\ U_b(k-1) \\ U_c(k-1) \end{bmatrix}}; i(k)={\begin{bmatrix} i_A(k)\\ i_B(k)\\ i_C(k) \\ i_a(k)\\ i_b(k) \\ i_c(k) \end{bmatrix}}; i(k-1)={\begin{bmatrix} i_A(k-1)\\ i_B(k-1)\\ i_C(k-1) \\ i_a(k-1)\\ i_b(k-1) \\ i_c(k-1) \end{bmatrix}}; \qquad\tag*{(9)} $$
Рис.2. Блоки моделирования асинхронного двигателя (а) и механической подсистемы привода (б) в Matlab Simulink

Аналитическое решение задач Коши для этапов подъема груза с использованием полученного описания электродвигателя и трехмассовой механической системы представляется затруднительным ввиду наличия большого количества переменных, а также производной по углу поворота ротора, в связи с чем полученное математическое описание анализировалось методом компьютерного моделирования в программном комплексе Matlab Simulink [18]. Модель привода подъема мостового крана, соответствующая уравнениям (2)-(9), показана на рис.2.

Обсуждение. Покажем использование приведенной на рис.2 модели для расчетов токов статора двигателя механизма подъема ремонтного мостового крана грузоподъемностью 32 т. Система моделировалась в режимах подъема груза номинальной массы и подъема груза массой, превышающей допустимую на 3, 15 и 25 %, соответственно при проведении вычислительных экспериментов в качестве массы груза m2 в системе уравнений (3) были заданы следующие значения: mI = 32000 кг, mII = 33000 кг, mIII = 36000 кг и mIV = 40000 кг; величина слабины каната l принята одинаковой для всех грузов. В качестве результатов моделирования получены величины токов статора электродвигателя подъема при отрыве от поверхности грузов (рис.3).

Анализируя полученные результаты можно отметить, что в рассматриваемой системе до момента времени t = 4,8 с происходит выбор слабины каната, вследствие чего ток статора одинаков для всех грузов, и его амплитуда в течение всего этапа не меняется.

Рис.3. Колебания тока статора электродвигателя при подъеме грузов различной массы (а), на этапе отрыва груза (б) 1 – выбор слабины каната; 2 – натяжение каната; 3 – отрыв груза

После того как первый этап подъема груза заканчивается, в канате возникает усилие, а значение тока возрастает с ростом натяжения каната механизма подъема. Условием окончания второго этапа является равенство усилия в канате веса груза, вследствие чего время окончания второго этапа для каждого из грузов различно, а разделение на рис.3, а показано условно. Для грузов большей массы длительность второго этапа увеличивается.

На третьем этапе подъема груза амплитудные значения тока колеблются с частотой, определяемой соотношением приведенных масс и жесткостью связей. Кроме изменения времени начала этапа отрыва груза, отличаются и амплитудные значения тока для разных масс, а также наблюдается смещение колебаний по фазе для грузов различной массы.

Результаты моделирования показывают, что уже за первые пять периодов сетевого напряжения после отрыва груза от опорной поверхности обнаруживается поддающаяся фиксации разница амплитуды тока статора (рис.3, б), которая для mII составила 1,8 %, для mIII – 12,7 %, а для mIV – 23,6 %. Появление различий в величинах амплитуды тока статора свидетельствует о чувствительности модели к вариации массы поднимаемого груза и подтверждает возможность диагностирования перегрузки за единицы периодов сетевого напряжения.

Моделирование процессов подъема грузов различной массы, аналогичное показанному выше, было проделано для ряда кранов различной грузоподъемности, применяемых на металлургических предприятиях. В качестве основы рассматривался парк оборудования предприятия АО «ЕВРАЗ ЗСМК» (Новокузнецк), насчитывающий 354 единицы мостовых кранов. В результате его анализа в рассматриваемый ряд мостовых кранов вошла 51 единица типоразмеров с мощностью привода 5-160 кВт и массой поднимаемого груза 5-200 т.

С использованием систем уравнений (2)-(9) и модели, показанной на рис.2, были определены величины токов статора асинхронных двигателей механизмов подъема для всех мостовых кранов рассматриваемого ряда. Для каждого из типоразмеров было выполнено моделирование ситуации перегрузки механизма подъема массой груза, превышающей номинальную на 1 %. Определялись значения фазного тока статора в первом полупериоде напряжения питающей сети при отрыве груза, и во всех случаях фиксировалось увеличение его амплитуды относительно опыта подъема с грузом номинальной массы со средневзвешенным по ряду значением отклонения 0,579 %. Для полученных результатов моделирования были исследованы статистические зависимости между мощностью привода N, грузоподъемностью крана W, увеличением амплитуды тока I и получены уравнения регрессии с коэффициентом линейной корреляции 0,75:

$$ \begin{gathered} I=r_{NI}\frac{\sigma_{I}}{\sigma_N}(N- \bar N )+\bar I \ или \ I=0,014N+0,009; \\ I=r_{WI}\frac{\sigma_{I}}{\sigma_W}(W- \bar W )+\bar I \ или \ I=0,01W+0,36.\end{gathered} \qquad\tag*{(10)}$$

Дополнительные исследования показали, что увеличение амплитуды тока статора с первого по пятый период сетевого напряжения носит линейный характер, как и увеличение амплитуды тока статора при повышении массы от 1 до 10 %.

Таким образом, уравнения регрессии (10) позволяют прогнозировать увеличение амплитуды тока статора при превышении массы груза для мостовых кранов, не вошедших в рассмотренный перечень типоразмеров, и, следовательно, выступают основой для алгоритма раннего диагностирования превышения массы груза для любого мостового крана металлургических предприятий и могут быть в дальнейшем применены для определения уставки срабатывания программно-аппаратного средства защиты.

Заключение. Результаты всех вычислительных экспериментов подтвердили следующее:

• на этапе выбора слабины каната величина тока статора электродвигателя механизма подъема не зависит от массы груза;

• время до начала третьего этапа с ростом массы груза увеличивается;

• амплитудные значения тока на этапе отрыва груза возрастают с увеличением массы груза, а поддающаяся фиксации разница появляется уже после пяти периодов сетевого напряжения с начала этапа;

• частота колебаний тока после отрыва груза зависит от его массы, но изменяется в небольших пределах.

Объем полученных результатов позволил определить статистические зависимости между величинами мощности привода подъема, грузоподъемности мостового крана и отклонениями амплитуды тока электродвигателя, которые являются основой для алгоритма раннего диагностирования превышения массы груза. Таким образом, контроль тока электродвигателя подъема позволяет создать защиту грузоподъемного механизма от превышения массы груза, не требующую внесения изменений в конструктивные элементы мостовых кранов. Быстродействие предлагаемого подхода существенно выше, чем у существующих механических способов контроля перегрузки.

Литература

1. Афанасьев А.С. Методика диагностирования дизелей по временным параметрам рабочего цикла / А.С.Афанасьев, А.А.Третьякова // Записки Горного института. 2015. Т. 214. С. 51-56. Методика диагностирования дизелей по временным параметрам рабочего цикла

2. Визуальное моделирование двухбалочного мостового крана как сложной динамической системы / Л.Н.Ахтулова, А.Л.Ахтулов, О.М.Кирасиров, В.А.Машонский // Омский научный вестник. 2014. № 1(127). С. 147-152. Визуальное моделирование двухбалочного мостового крана как сложной динамической системы

3. Голдобина Л.А. Обеспечение безопасности строительно-монтажных работ при возведении зданий и сооружений / Л.А.Голдобина, П.А.Деменков, О.В.Трушко // Записки Горного института. 2019. Т. 239. С. 583-595. DOI: 10.31897/pmi.2019.5.583

4. Ещин Е.К. Расчеты динамических режимов работы электроприводов самоходных горных машин // Записки Горного института. 2018. Т.233. С.534-538. DOI: 10.31897/PMI.2018.5.534

5. Оценка рисков потенциально опасных металлургических объектов за пределами гарантийных сроков эксплуатации в менеджменте качества промышленного предприятия / Ю.А.Извеков, Е.М.Гугина, А.Л.Анисимов, В.В.Шеметова // Современные наукоемкие технологии. 2018. № 11 (Ч. 2). C. 179-182. Оценка рисков потенциально опасных металлургических объектов за пределами гарантийных сроков эксплуатации в менеджменте качества промышленного предприятия

6. Патент № 2144901 РФ. Устройство для измерения и ограничения грузоподъемности лебедки / Г.Г.Платонов, А.В.Кудрявцев, С.Н.Сологубов. Опубл. 27.01.2000. Бюлл. № 3.

7. Патент № 2381984 РФ. Устройство для определения веса груза мостового крана / Х.М.Кадыров. Опубл. 20.02.2010. Бюлл. № 5.

8. Патент № 2464220 РФ. Тензометрическая ось для измерения нагрузки на крюке грузоподъемного крана / В.А.Потапов, Ю.Ф.Тимин, М.В.Корников. Опубл. 20.10.2012. Бюлл. № 29.

9. Патент № 2483016 РФ. Ограничитель нагрузки мостового крана / В.А.Коровин, К.В.Коровин. Опубл. 27.05.2013. Бюлл. № 15.

10. Полищук В.И. Развитие теории построения систем диагностики синхронных машин: Автореф. дис. … д-ра техн. наук. Самара: Cамарский государственный технический университет, 2016. 36 с. Развитие теории построения систем диагностики синхронных машин

11. Савельев А.Н. Диагностика аварийных режимов ленточных транспортеров в горно-металлургической промышленности / А.Н.Савельев, М.В.Кипервассер, Д.С.Аниканов // Известия высших учебных заведений. Черная металлургия. 2015. Т. 58. № 12. C. 906-911. DOI: 10.17073/0368-0797-2015-12-906-911

12. Сафин Н.Р. Диагностика неисправностей вентиляторных установок с помощью спектрального анализа токов статора / Н.Р.Сафин, В.А.Прахт, В.А.Дмитриевский // Энергобезопасность и энергосбережение. 2016. № 4. C. 37-42. DOI: 10.18635/2071-2219-2016-4-37-42

13. Семыкина И.Ю. Исследование переходного процесса в цепи тока при дефекте мультипликатора турбокомпрессора и разработка структурной схемы системы обнаружения дефекта мультипликатора / И.Ю.Семыкина, М.В.Кипервассер, А.В.Герасимук // Вестник Чувашского университета. 2017. № 3. С. 138-144. Исследование переходного процесса в цепи тока при дефекте мультипликатора турбокомпрессора и разработка структурной схемы системы обнаружения дефекта мультипликатора

14. Турышева Е.С. Совершенствование процесса автоматической защиты гидрофицированного крана от перегрузки и опрокидывания: Автореф. дис. … канд. техн. наук. Красноярск: Сибирский федеральный университет, 2009. 18 с.

15. Abdel-Rahman E.M. Dynamics and Control of Cranes: A Review / E.M.Abdel-Rahman, A.H.Nayfeh, Z.N.Masoud // Journal of Vibration and Control. 2003. Vol. 9(7), P. 863-908. DOI: 10.1177/1077546303009007007

16. Ansari А. Mathematical Model of Asynchronous Machine in Matlab Simulink / A.Ansari, D.Deshpande // International Journal of Engineering Science and Technology. 2010. Vol. 2(5). P. 1260-1267. DOI: 10.1.1.165.5447

17. Arya M.K. Transient Analysis of Three Phase Squirrel Cage Induction Machine using Matlab / M.K.Arya, S.Wadhwani // International Journal of Engineering Research and Applications (IJERA). 2009. Vol. 1. Iss. 3. P. 918-922. DOI: 10.1.1.300.1599

18. Chee-Mun OngDynamic simulation of electric machinery using Matlab / Simulink. Prentice Hall PTR, 1997. 626 p.

19. Chmurawa M. Modeling of bridge cranes for dimensioning needs of their load-carrying structures / M.Chmurawa, D.Gąska // The International Journal of Ingenium. 2005. Vol. 4. P. 409-414.

20. Electric Motor Handbook / J.L.Kirtley, H.W.Beaty, N.K.Ghai, S.B.Leeb, R.H.Lyon. New York, USA: McGraw-Hill Education, 1998. 400 p. Electric Motor Handbook

21. Fang Y. A framework for real-time pro-active safety assistance for mobile crane lifting operations / Y.Fang, Y.K.Cho, J.Chen // Automation in Construction. 2016. Vol. 72. Part 3. P. 367-379. DOI: 10.1016/j.autcon.2016.08.025

22. Gąska D. Numeryczno-statystyczna metoda oceny nośności i stateczności stalowych ustrojów nośnych dźwignic: Praca doktorska. Katowice: Politechnika Śląska, 2007. 204 p.

23. Gu Liang-Yao.Research on intelligent monitoring and protection system of distributive multi-tower cranes / Gu Liang-Yao, Yang Yue, Chen Feng // International Conference on Computer Application and System Modeling, 2010, October 22-24, Taiyuan, China. IEEE, 2010. P. 540-543. DOI: 10.1109/ICCASM.2010.5623149

24. He Wei. Adaptive Control of a Flexible Crane System With the Boundary Output Constraint / He Wei, Shuang Zhang, Shuzhi Sam Ge // IEEE Transactions On Industrial Electronics. 2014. Vol. 61. Iss. 8. P. 4126-4133. DOI: 10.1109/TIE.2013.2288200

25. Hong K.-S. Lumped mass models of mobile cranes / K.-S.Hong, U.H. Shah // Dynamics and Control of Industrial Cranes: Advances in Industrial Control. Singapore: Springer Nature Singapore Pte Ltd, 2019. P. 49-65. DOI: 10.1007/978-981-13-5770-1_4

26. Krishnan R. Motor Drives: Modelling, Analysis and Control. New Jersy, USA: Upper Saddle River, 2001. 626 p.

27. Matyja T. Modeling of the Lift Crane Vibration Caused by the Lifting Loads / T.Matyja, A.Sładkowski // Zdvihací Zařízení v Teorii a Praxi, 2007, April 10-11, Brno, Czech Republic. VUT Brno, 2007. P. 98-105. Modeling of the Lift Crane Vibration Caused by the Lifting Loads

28. Messineo S. Offshore Crane Control Based on Adaptive External Models / S.Messineo, A.Serrani // Automatica. 2009. Vol. 45. Iss. 11. P. 2546-2556. DOI: 10.1016/j.automatica.2009.07.032

29. Parameter Imputation Methods and Typical Parameter Analysis of Load Transient Model of Induction Motor / T.Liang, N.Zhou, T.Lu, H.Wu, P.Ju // Dianli Xitong Zidonghua – Automation of Electric Power Systems. 2020. Vol. 44. № 1. P. 74-82. DOI: 10.7500/AEPS20190507002

30. Ren Z. A Crane Overload Protection Controller for Blade Lifting Operation Based on Model Predictive Control / Z.Ren, R.Skjetne, Z.Gao // Energies. 2019. Vol. 12(1). P. 1-22. DOI: 10.3390/en12010050

31. Ren Z. Modeling and Control of Crane Overload Protection During Marine Lifting Operation Based on Model Predictive Control / Z.Ren, R.Skjetne, Z.Gao // International Conference on Ocean, Offshore and Arctic Engineering, 2017, June 25-30, Trondheim, Norway. ASME, 2017. Vol. 57779. 9 p. DOI: 10.1115/OMAE2017-62003

32. Szpytko J. Control and Fault Diagnosis System of Crane Operation // IFAC Proceedings Volumes. 2001. Vol. 34. Iss. 17. P. 237-242. DOI: 10.1016/S1474-6670(17)33286-X

Похожие статьи

Автоматизированная оцифровка круговых диаграмм
2021 Н. В. Васильева, А. В. Бойков, О. О. Ерохина, А. Ю. Трифонов
Результаты комплексных геофизических исследований по поиску склепов на территории загородного некрополя Херсонеса Таврического в Карантинной балке
2021 В. В. Глазунов, А. С. Агеев, Г. Д. Горелик, Т. В. Сарапулкина
Прогнозирование горно-геологических процессов на основе анализа подземного пространства рудника Купол как многокомпонентной системы (Чукотский автономный округ, Анадырский р-н)
2021 Р. Э. Дашко, И. С. Романов
Проведение промышленных взрывов вблизи газопроводов
2021 С. В. Хохлов, С. Т. Соколов, Ю. И. Виноградов, И. Б. Френкель
Нефтегазоносность малоизученной части северо-запада Тимано-Печорской нефтегазоносной провинции по результатам бассейнового моделирования
2021 О. М. Прищепа, И. С. Боровиков, Е. И. Грохотов
Динамическое моделирование промышленного цикла кристаллизации гиббсита
2021 В. О. Голубев, Т. Е. Литвинова