Centrifugal Pump

Центробежный насос с выбором опций параметризации

  • Библиотека:
  • Simscape/Жидкости/Гидравлика (изотермическая )/Насосы и двигатели

Описание

Блок Centrifugal Pump представляет центробежный насос любого типа как модель , основанная на таблице данных. В зависимости от данных, перечисленных в каталоге производителя или табличных данных для вашего конкретного насоса, можно выбрать одну из следующих опций модели параметризации:

  • By approximating polynomial - Задайте значения полиномиальных коэффициентов. Эти значения могут быть определены аналитически или экспериментально, в зависимости от имеющихся данных. Это метод по умолчанию.

  • By two 1D characteristics: P-Q and N-Q - Предоставить табличные данные характеристик Q перепада давления P и мощности сопротивления N от производительности насоса. Перепад давления и приводная мощность определяются одномерным поиском таблицы. У вас есть выбор из двух методов интерполяции и двух методов экстраполяции.

  • By two 2D characteristics: P-Q-W and N-Q-W - Предоставить табличные данные характеристик перепадов давления P и приводной мощности N в зависимости от производительности насоса Q при различных скоростях вращения W. Перепад давления и приводная мощность определяются двумерным поиском таблицы. У вас есть выбор из двух методов интерполяции и двух методов экстраполяции.

Эти опции параметризации дополнительно описаны более подробно:

Связи P и T являются гидравлическими портами, сопоставленными с выходным и входным отверстиями насоса, соответственно. Связь S является портом механического вращений, сопоставленным с ведущим валом насоса. Блок имеет положительное направление от порта T до порта P. Это означает, что насос перекачивает жидкость от T до P, когда его ведущий вал S вращается в глобально присвоенном положительном направлении.

Параметризация насоса путем аппроксимации полинома

Если вы устанавливаете параметр Model parameterization равным By approximating polynomialнасос параметризован полиномом, коэффициенты которого определяются аналитически или экспериментально для выбранной угловой скорости в зависимости от имеющихся данных. Характеристики насоса при других скоростях вращения определяются с помощью законов подобия.

Аппроксимирующий полином получен из уравнения момента импульса Эйлера, уравнения 1 и 2, которое для данного насоса, скорости вращения и жидкости может быть представлено как следующее:

pref=kpEpHLpD(1)

где

p refПерепад давления на насосе для режима ссылки, характеризующийся ссылкой скоростью вращения и плотностью
kКоэффициент коррекции. Коэффициент вводится для расчета размерных колебаний, несоответствия блейдов, объемов блейдов, внутреннего трения жидкости и так далее. Коэффициент должен быть установлен равным 1, если аппроксимационные коэффициенты определены экспериментально.
p EДавление Эйлера
p HLПадение давления из-за гидравлических потерь в полостях насоса
p DПадение давления, вызванное отклонениями производительности насоса от его номинального (номинального) значения

Давление Эйлера, pE, определяется уравнением Эйлера для центробежных машин в уравнениях 1 и 2 на основе известных размеров насоса. Для существующего насоса, работающего с постоянной скоростью вращения и заданной жидкостью, давление Эйлера может быть аппроксимировано уравнением

pE=ρref(c0c1qref)

где

ρrefПлотность жидкости
<reservedrangesplaceholder1> 0, <reservedrangesplaceholder0> 1Аппроксимация коэффициентов. Они могут быть определены или аналитически из уравнения Эйлера (уравнения 1 и 2) или экспериментально.
q refОбъемная подача насоса при исходном режиме

Падение давления из-за гидравлических потерь в полостях насоса, p HL, аппроксимируется уравнением

pHL=ρrefc2qref2

где

ρrefПлотность жидкости
<reservedrangesplaceholder0> 2Аппроксимационный коэффициент
q refОбъемная подача насоса при исходном режиме

Профиль лопасти определяется для скорости заданной жидкости, и отклонение от этой скорости результатов в потере давления из-за несоответствия между скоростью жидкости и скоростью профиля лопасти. Это падение давления, p D, оценивается уравнением

pD=ρrefc3(qDqref)2

где

ρrefПлотность жидкости
<reservedrangesplaceholder0> 3Аппроксимационный коэффициент
q refОбъемная подача насоса при исходном режиме
q DРасчетная поставка насоса (номинальная поставка)

Результат аппроксимации полинома принимает форму:

pref=ρref(k(c0c1qref)c2qref2c3(qDqref)2)(2)

Характеристики насоса, аппроксимированные четырьмя коэффициентами c 0, c 1, c 2 и c 3, определяются для конкретной жидкости и определенной угловой скорости ведущего вала насоса. Эти два параметра соответствуют, соответственно, Reference density и Reference angular velocity параметрам в диалоговом окне блока. Чтобы применить характеристики для другой скорости и плотности, используются законы подобия. С помощью этих законов доставка в ссылочном режиме, которая соответствует заданной производительности насоса и скорости вращения, вычисляется выражением

qref=qωrefω(3)

где q и, - мгновенные значения производительности насоса и скорости вращения. Затем перепадом давления p ссылку на исходный режим, вычисленный уравнением 2 и преобразованный в перепады давления, p при текущей скорости вращения и плотности

p=pref(ωωref)2ρρref

Уравнение 2 описывает характеристику насоса для В > 0 и q > = 0. Вне этой области значений характеристика аппроксимируется следующими отношениями:

p={kleakqдля ω<=0pmaxkleakqдля ω>0,q<0kleak(qqmax)для ω>0,q>qmax(4)

qmax=b+b2+4ac2a

a=(c2+c3)α2

b=(kc12c3qD)α

c=kc0c3qD2

α=ωωref

qmax=ρ1α2(kc0c3qD2)

где

k утечкиКоэффициент сопротивления утечкам
q maxМаксимальная производительность насоса при заданной скорости вращения. Доставка определяется из уравнения 2 при p = 0.
p maxМаксимальное давление насоса при заданной скорости вращения. Давление определяется из уравнения 2 при q = 0.
kКоэффициент коррекции, как описано в уравнении 1.

Гидравлическая степень на выходе насоса при начальных условиях

Nhyd=prefqref

Гидравлический выход степени при произвольных скоростях вращения и плотности определяется законами подобия

N=Nref(ωωref)ρρref

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

Nhyd0=pErefqref(ωωref)3

где

<reservedrangesplaceholder0> hyd0Теоретическая гидравлическая степень насоса
p ЭрефДавление Эйлера. Теоретическое давление, развиваемое насосом до потерь, связанных с гидравлическими потерями и отклонением от производительности проекта.

Потери на трение аппроксимируются отношением:

Nfr=(T0+kpp)ω

где

NfrКоэффициент трения степени
T0Постоянный крутящий момент на ведущем валу, связанный с подшипниками вала, трением уплотнения и так далее
kpСоотношение крутящий момент и давление, которое характеризует влияние давления на крутящий момент ведущего вала

Мощность и крутящий момент на ведущем вале насоса (мощность тормоза N mech и крутящий момент привода T)

Nmech=Nhyd0+Nfr

T=Nmechω

Общий КПД η насоса вычисляется как

η=NhydNmech

Параметризация насоса перепадом давления и приводной мощностью в зависимости от производительности насоса

Если вы устанавливаете параметр Model parameterization равным By two 1D characteristics: P-Q and N-Qхарактеристики насоса вычисляются с помощью двух одномерных интерполяций таблицы: для перепада давления на основе производительности насоса и для приводной мощности насоса на основе производительности насоса. Обе характеристики заданы при одной и той же угловой скорости в ref (Reference angular velocity) и той же плотности жидкости в ref (Reference density ).

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

qref=qωrefω

где q - текущая производительность насоса. Затем перепадом давления на насосе при токе скорости вращения

p=pref(ωωref)2ρρref

где pref - перепад давления, определяемый из характеристики P-Q на qref производительности насоса.

Приводная мощность определяется уравнением

N=Nref(ωωref)3ρρref

где Nref - базовая мощность сопротивления, полученная из N-Q характеристики при производительности насоса qref.

Крутящий момент на приводном валу насоса вычисляется уравнением T = N/

Параметризация насоса перепадом давления и Приводной мощности в зависимости от Производительности насоса при различных Скоростях вращения

Если вы устанавливаете параметр Model parameterization равным By two 2D characteristics: P-Q-W and N-Q-Wхарактеристики насоса считываются из двухмерных интерполяций таблицы: для перепада давления на основе производительности насоса и скорости вращения и для приводной мощности насоса на основе производительности насоса и скорости вращения.

И перепадом давления, и тормозной мощностью масштабируются, если плотность жидкости

p=prefρρref

N=Nrefρρref

где pref и Nref - перепад давления и приводной мощности, полученные на графиках.

Основные допущения и ограничения

  • Сжимаемостью жидкости пренебрегают.

  • Насос вращается в положительном направлении со скоростью, которая больше или равна нулю.

  • Обратное течение через насос допускается только при неподвижном валу.

Порты

Сохранение

расширить все

Гидравлический порт сопоставлен с всасывающим или входным отверстием насоса.

Гидравлический порт сопоставлен с выходным отверстием насоса.

Механический вращательный порт сопоставлен с ведущим валом насоса.

Параметры

расширить все

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

  • By approximating polynomial - Задайте значения полиномиальных коэффициентов. Эти значения могут быть определены аналитически или экспериментально, в зависимости от имеющихся данных. Зависимость между характеристиками насоса и скоростью вращения определяется из законов подобия.

  • By two 1D characteristics: P-Q and N-Q - Предоставить табличные данные перепада давления и приводной мощности в зависимости от характеристик производительности насоса. Перепад давления и приводная мощность определяются одномерным поиском таблицы. У вас есть выбор из двух методов интерполяции и двух методов экстраполяции. Зависимость между характеристиками насоса и скоростью вращения определяется из законов подобия.

  • By two 2D characteristics: P-Q-W and N-Q-W - Предоставить табличные данные характеристик перепада давления и приводной мощности в зависимости от производительности насоса при различных скоростях вращения. Перепад давления и приводная мощность определяются двумерным поиском таблицы. У вас есть выбор из двух методов интерполяции и двух методов экстраполяции.

Аппроксимация коэффициента, c0 в описании блока, предшествующем. Этот параметр используется, если Model parameterization установлено на By approximating polynomial.

Аппроксимация коэффициента, c1 в описании блока, предшествующем. Этот параметр используется, если Model parameterization установлено на By approximating polynomial.

Аппроксимация коэффициента, c2 в описании блока, предшествующем. Этот коэффициент учитывает гидравлические потери в насосе. Этот параметр используется, если Model parameterization установлено на By approximating polynomial.

Аппроксимация коэффициента, c3 в описании блока, предшествующем. Этот коэффициент учитывает дополнительные гидравлические потери, вызванные отклонением от номинальной производительности. Этот параметр используется, если Model parameterization установлено на By approximating polynomial.

Коэффициент, обозначенный как k в описании блока, предшествующем, учитывает размерные колебания, несоответствие блейдов, объемы блейдов, внутреннее трение жидкости и другие факторы, которые уменьшают теоретическое давление Эйлера. Этот параметр используется, если Model parameterization установлено на By approximating polynomial.

Номинальная производительность насоса. Профиль блейдов, входное отверстие насоса и выход насоса имеют форму для этой конкретной производительности. Отклонение от этой производительности приводит к увеличению гидравлических потерь. Этот параметр используется, если Model parameterization установлено на By approximating polynomial.

Скорость вращения ведущего вала, при которой определяются характеристики насоса. Этот параметр используется, если Model parameterization установлено на By approximating polynomial или By two 1D characteristics: P-Q and N-Q.

Плотность жидкости, при которой определяются характеристики насоса.

Коэффициент сопротивления утечкам (см. Уравнение 4). Этот параметр используется, если Model parameterization установлено на By approximating polynomial.

Крутящий момент трения на валу при нулевой скорости. Этот параметр используется, если Model parameterization установлено на By approximating polynomial.

Коэффициент, который обеспечивает зависимость между крутящим моментом и давлением насоса. Значение по умолчанию 1e-6 Н * м/Па. Этот параметр используется, если Model parameterization установлено на By approximating polynomial.

Укажите вектор производительностей насоса, как одномерный массив, который будет использоваться вместе с вектором перепадов давления для определения насосной характеристики P-Q. Векторные значения должны быть строго увеличены. Значения могут быть неоднородно разнесены. Минимальное количество значений зависит от метода интерполяции: необходимо задать как минимум два значения для линейной интерполяции, как минимум три значения для сплайна-интерполяции. Этот параметр используется, если Model parameterization установлено на By two 1D characteristics: P-Q and N-Q.

Задайте вектор перепада давления на насосе как одномерный массив. Вектор будет использоваться вместе с вектором производительности насоса, чтобы задать насосную характеристику P-Q. Вектор должен быть того же размера, что и вектор производительности насоса для таблицы P-Q. Этот параметр используется, если Model parameterization установлено на By two 1D characteristics: P-Q and N-Q.

Укажите вектор производительностей насоса, как одномерный массив, который будет использоваться вместе с вектором приводной мощности насоса, чтобы задать характеристику насоса N-Q. Векторные значения должны быть строго увеличены. Значения могут быть неоднородно разнесены. Минимальное количество значений зависит от метода интерполяции: необходимо задать как минимум два значения для линейной интерполяции, как минимум три значения для сплайна-интерполяции. Этот параметр используется, если Model parameterization установлено на By two 1D characteristics: P-Q and N-Q.

Задайте вектор приводной мощности насоса как одномерный массив. Вектор будет использоваться вместе с вектором производительности насоса, чтобы задать характеристику насоса N-Q. Вектор должен быть того же размера, что и вектор производительности насоса для таблицы N-Q. Этот параметр используется, если Model parameterization установлено на By two 1D characteristics: P-Q and N-Q.

Укажите вектор скоростей вращения как одномерный массив, который будет использоваться для вычисления характеристик насоса P-Q-W и N-Q-W. Векторные значения должны быть строго увеличены. Значения могут быть неоднородно разнесены. Минимальное количество значений зависит от метода интерполяции: необходимо задать как минимум два значения для линейной интерполяции, как минимум три значения для сплайна-интерполяции. Этот параметр используется, если Model parameterization установлено на By two 2D characteristics: P-Q-W and N-Q-W.

Задайте вектор производительностей насоса, как одномерный массив, который будет использоваться вместе с вектором скоростей вращения и матрицей перепада давления, чтобы задать характеристику насоса P-Q-W. Векторные значения должны быть строго увеличены. Значения могут быть неоднородно разнесены. Минимальное количество значений зависит от метода интерполяции: необходимо задать как минимум два значения для линейной интерполяции, как минимум три значения для сплайна-интерполяции. Этот параметр используется, если Model parameterization установлено на By two 2D characteristics: P-Q-W and N-Q-W.

Задайте перепады давления на насосе как m-by- n матрица, где m - количество значений производительности насоса P-Q-W и n количество скоростей вращения. Эта матрица будет определять характеристику насос P-Q-W вместе с векторами производительности насоса и скорости вращения. Каждое значение в матрице задает перепад давления для определенной комбинации производительности насоса и скорости вращения. Размер матрицы должен совпадать с размерностями, заданными векторами производительности насоса и скорости вращения. Этот параметр используется, если Model parameterization установлено на By two 2D characteristics: P-Q-W and N-Q-W.

Задайте вектор производительностей насоса, как одномерный массив, который будет использоваться вместе с вектором скоростей вращения и приводной мощности матрицей, чтобы задать характеристику насоса N-Q-W. Векторные значения должны быть строго увеличены. Значения могут быть неоднородно разнесены. Минимальное количество значений зависит от метода интерполяции: необходимо задать как минимум два значения для линейной интерполяции, как минимум три значения для сплайна-интерполяции. Этот параметр используется, если Model parameterization установлено на By two 2D characteristics: P-Q-W and N-Q-W.

Задайте приводную мощность насоса как m-by- n матрица, где m - количество значений производительности насоса N-Q-W и n количество скоростей вращения. Эта матрица будет определять характеристику насоса N-Q-W вместе с векторами производительности насоса и скорости вращения. Каждое значение в матрице задает приводную мощность для определенной комбинации производительности насоса и скорости вращения. Размер матрицы должен совпадать с размерностями, заданными векторами производительности насоса и скорости вращения. Этот параметр используется, если Model parameterization установлено на By two 2D characteristics: P-Q-W and N-Q-W.

Угловая скорость вала, которая указывает порог перехода между прямыми и обратными течениями. Переходная область задается около 0 рад/с между положительным и отрицательным значениями порога скорости вращения. В пределах этой переходной области вычисленные уровни утечек и крутящий момент трения регулируются согласно α переходного периода, чтобы гарантировать плавность перехода от одного режима к другому.

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

  • Linear - Выберите эту опцию, чтобы получить лучшую эффективность.

  • Smooth - Выберите эту опцию, чтобы создать непрерывную кривую или поверхность с непрерывными производными первого порядка.

Этот параметр используется, если Model parameterization установлено на By By two 1D characteristics: P-Q and N-Q или By two By two 2D characteristics: P-Q-W and N-Q-W. Для получения дополнительной информации об алгоритмах интерполяции смотрите страницы PS Lookup Table (1D) и PS Lookup Table (2D) блочных ссылок.

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

  • Linear - Выберите эту опцию, чтобы создать кривую или поверхность с непрерывными производными первого порядка в области экстраполяции и на контуре с областью интерполяции.

  • Nearest - Выберите эту опцию, чтобы создать экстраполяцию, которая не идет выше высшей точки в данных или ниже самой нижней точки в данных.

Этот параметр используется, если Model parameterization установлено на By By two 1D characteristics: P-Q and N-Q или By two By two 2D characteristics: P-Q-W and N-Q-W. Для получения дополнительной информации об алгоритмах экстраполяции смотрите страницы PS Lookup Table (1D) и PS Lookup Table (2D) блочных ссылок.

Ссылки

[1] T.G. Hicks, T.W. Edwards, Pump Application Engineering, McGraw-Hill, NY, 1971

[2] I.J. Karassic, J.P. Messina, P. Cooper, C.C. Heald, Pump Handbook, Third edition, McGraw-Hill, NY, 2001

Расширенные возможности

Генерация кода C/C + +
Сгенерируйте код C и C++ с помощью Coder™ Simulink ®

.
Введенный в R2007a
Для просмотра документации необходимо авторизоваться на сайте