Способ определения расходной характеристики гидравлического тракта для области перехода от турбулентного к ламинарному режиму истечения

Изобретение относится к области измерительной техники и может быть использовано для расчета пропускной способности проектируемых гидравлических трактов транспортных и дозирующих систем в химической, нефтехимической, авиационной, текстильной, лакокрасочной и других отраслях промышленности, в частности, - узлов транспортирования и дозирования клеевых составов при сборке малогабаритных изделий. Способ характеризуется тем, что для области перехода от турбулентного к ламинарному режиму истечения модель тракта с проектными конструктивами местных сопротивлений и участков соединительных трубопроводов длиной li≤Li с внутренним диаметром di с концевой раздачей через запорное устройство с цилиндрическим сливным наконечником с внутренним диаметром dн проливают водопроводной водой посредством ее многократного порционного дозирования в мерную тару при различных величинах действующего напора Hj, ступенчато изменяемых от максимального Hj,max, соответствующего турбулентному режиму истечения, до минимального Hj,min, соответствующего началу перехода от струйного к капельно-струйному режиму истечения, значений, которые во время каждой j-й операции проливки поддерживают на заданном постоянном уровне, где j=1, 2 …, n=31 - число проливок в каждом цикле измерений, и оперативно, с помощью электронных таблиц программы «Microsoft Excel», для каждой j-й операции проливки последовательно вычисляют объемный расход - Qj [см3/c]=Vj/tj, соответствующий напору Hj, где объем дозы Vj [см3] регистрируют с помощью взвешивания на электронных весах, а время дозирования tj[c] задают и отмеривают с помощью электронного таймера; скорость истечения воды из сливного наконечника - vн,j[см/c]=Qj/fн, потерю напора - hн,j [см]=α⋅vн,j2/2g и число Рейнольдса потока - Rej=vн,jdнв на его выходе, где коэффициент Кориолиса α принимают равным 1,05, как при турбулентном режиме истечения, а кинематическую вязкость воды νв рассчитывают по ее температуре в данном цикле проливок; скорости течения - vi,j[см/c]=vн,j(dн/di)2, числа Рейнольдса - Rei,j=vi,jdiв, коэффициенты трения - λi,j=64/Rei,j и потери напора на трение - hi,j [см]=λi,j(li/di)(dн/di)4 на отдельных участках трубопровода; суммарные потери напора на трение - hт,j [см]=Σi[hi,j], на местных сопротивлениях в тракте - hм,j[см]=Нj-hн,j-hт,j и суммарные коэффициенты потерь напора на местных сопротивлениях - ξj=2ghм,j/vн,j2. Далее полученную в табличном виде экспериментальную расходную характеристику тракта Qj=F(Hj) делят на зоны турбулентного, условно переходного и ламинарного режимов истечения; в программе «Advanced Grapher» для каждой выделенной зоны строят графики функций ξj,s=f(Rej,s), где индексы s=1, 2, 3 относятся к турбулентной, условно переходной и ламинарной зонам, соответственно, проводят их регрессионный анализ и получают для каждой из выделенных зон уравнения аппроксимирующих функций: ξ1=f1(Re), ξ2=f2(Re), ξ3=f3(Re), при этом для каждого s выбирают степенную или гиперболическую, по минимальному значению стандартного отклонения, аппроксимацию; строят графики аппроксимирующих функций и определяют точки Re1 и Re2 пересечения графиков: Re1 (ξ1 и ξ2) и Re2 (ξ2 и ξ3), фиксирующих, соответственно, верхнюю (число Re1) и нижнюю (число Re2) границы фактической переходной зоны прогнозируемой расходной характеристики натурной жидкости с вязкостью ν; для последней задают ряд значений расхода Qk=Qj, скоростей истечения vн,k=vн,j из сливного наконечника и потерь напора на его выходе - hн,k=hн,j, совпадающих с полученными на воде значениями; для каждого значения k=j и для каждой из выделенных зон режимов истечения натурной жидкости вычисляют соответствующие этим зонам числа Рейнольдса: Rek,ssvн,k2/2g, при этом числа Rek,1>Re1 относятся к зоне турбулентного режима истечения, числа Rek,3<Re2 - к зоне ламинарного режима, а числа Re2<Rek,2<Re1 - к фактической переходной зоне прогнозируемой расходной характеристики; далее по полученным ранее уравнениям ξs=fs(Re) аппроксимирующих функций вычисляют суммарные коэффициенты потерь напора на местных сопротивлениях ξk=f(Rek), скорости течения - vi,k=vн,k(dн/di)2, числа Рейнольдса - Rei„k=vi,kdi/ν, коэффициенты трения - λi,k=64/Rei,k и потери напора на трение - hi,ki,k(Li/di)(dн/di)4 на отдельных участках трубопровода, суммарные потери напора на трение - hт,k[см]=Σi[hi,k] и на местных сопротивлениях в тракте - hм,kk(vн,k2/2g) и, наконец, - расчетные, для каждой из выделенных зон, значения располагаемого напора - Нрk=hн,k+hт,k+hм,k. 3 з.п. ф-лы, 1 табл., 1 ил.

 

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

Основной задачей, которую приходится решать при проектировании тракта транспортирования жидких продуктов, является определение его пропускной способности (объемного расхода) Q при различных значениях действующего напора Н, т.е. его расходной характеристики - функции Q=F(H). Предлагаемый способ определения расходной характеристики основан на анализе уравнения Бернулли для установившегося течения несжимаемой жидкости. Это уравнение для гидравлического тракта, состоящего из нескольких местных сопротивлений, соединенных участками трубопроводов длиной Li с внутренним диаметром di, с концевой раздачей через запорное устройство с цилиндрическим сливным наконечником с внутренним диаметром dн имеет вид:

где fн=[3,14⋅dн2/4] - площадь проходного сечения наконечника, g - ускорение силы тяжести, α - коэффициент Кориолиса потока на выходе сливного наконечника, λi - коэффициент сопротивления трения на i-м участке трубопровода, ξ - суммарный коэффициент потерь напора на местных сопротивлениях в тракте, который определяется конструктивными параметрами местных сопротивлений и зависит от числа Рейнольдса потока (Re) на выходе сливного наконечника. С учетом равенства Q=vнfн, где vн - скорость истечения жидкости из сливного наконечника, уравнение (1) принимает вид:

При заданных значениях параметров vн, dн и di и известной кинематической вязкости жидкости ν по общепринятым в гидравлике соотношениям легко рассчитывается коэффициент трения λi на отдельных участках трубопровода. Для ламинарного режима течения λi=64/Rei, где число Рейнольдса Rei=vidi/ν, vi=vн(dн/di)2 - скорость течения на i-м участке трубопровода. Таким образом, для определения интересующей нас функции Q=F(H) необходимо определить входящие в уравнение (2) заранее неизвестные параметры α и ξ.

Известен способ определения коэффициента гидравлического сопротивления тракта (RU 2240525 С1, 20.11.2004).

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

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

Известен также способ определения расходной характеристики гидравлического тракта и устройство для его осуществления (RU 2680987 С1, 01.03.2019.), принятый за прототип.

Недостатком способа-прототипа является недостаточная точность и достоверность определения и прогнозирования расходной характеристики проектируемого тракта для области перехода от турбулентного к ламинарному режиму истечения из-за использования при расчетах суммарного коэффициента потерь напора на местных сопротивлениях ξ единой для всей исследуемой области аппроксимирующей функции ξ=f(Re)=a⋅Reb (где а и b - коэффициенты аппроксимации); при этом не учитывается различный характер указанной функциональной зависимости для зон турбулентного и ламинарного режимов истечения, а также для зоны перехода от турбулентного к ламинарному режиму.

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

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

Технический результат предлагаемого способа определения расходной характеристики гидравлического тракта, характеризующегося тем, что для области перехода от турбулентного к ламинарному режиму истечения модель тракта с проектными конструктивами местных сопротивлений и участков соединительных трубопроводов длиной li≤Li с внутренним диаметром di с концевой раздачей через запорное устройство с цилиндрическим сливным наконечником с внутренним диаметром dн проливают водопроводной водой посредством ее многократного порционного дозирования в мерную тару при различных величинах действующего напора Hj, ступенчато изменяемых от максимального Hj,max, соответствующего турбулентному режиму истечения, до минимального Hj,min, соответствующего началу перехода от струйного к капельно-струйному режиму истечения, значений, которые во время каждой j-й операции проливки поддерживают на заданном постоянном уровне, где j=1, 2 …, n=31 - число проливок в каждом цикле измерений, и оперативно, с помощью электронных таблиц программы «Microsoft Excel», для каждой j-й операции проливки последовательно вычисляют объемный расход - Qj [см3/c]=Vj/tj, соответствующий напору Hj, где объем дозы Vj [см3] регистрируют с помощью взвешивания на электронных весах, а время дозирования tj[c] задают и отмеривают с помощью электронного таймера; скорость истечения воды из сливного наконечника -vн,j[см/с]=Qj/fн, потерю напора - hн,j [см]=α⋅vн,j2/2g и число Рейнольдса потока - Rej=vн,jdнв на его выходе, где коэффициент Кориолиса α принимают равным 1,05 как при турбулентном режиме истечения, а кинематическую вязкость воды νв рассчитывают по ее температуре в данном цикле проливок; скорости течения - vi,j[см/c]=vн,j(dн/di)2, числа Рейнольдса - Rei,j=vi,jdiв, коэффициенты трения - λi,j=64/Rei,j и потери напора на трение - hi,j [см]=λi,j(li/di)(dн/di)4 на отдельных участках трубопровода; суммарные потери напора на трение - hт,j [см]=Σi[hi,j], на местных сопротивлениях в тракте - hм,j[см]=Hj-hн,j-hт,j и суммарные коэффициенты потерь напора на местных сопротивлениях - ξj=2ghм,j/vн,j2; при этом полученную в табличном виде экспериментальную расходную характеристику тракта Qj=F(Hj) делят на зоны турбулентного, условно переходного и ламинарного режимов истечения; в программе «Advanced Grapher» для каждой выделенной зоны строят графики функций ξj,s=f(Rej,s), где индексы s=1, 2, 3 относятся к турбулентной, условно переходной и ламинарной зонам, соответственно, проводят их регрессионный анализ и получают для каждой из выделенных зон уравнения аппроксимирующих функций: ξ1=f1(Re), ξ2=f2(Re), ξ3=f3(Re), при этом для каждого s выбирают степенную или гиперболическую, по минимальному значению стандартного отклонения, аппроксимацию; строят графики аппроксимирующих функций и определяют точки Re1 и Re2 пересечения графиков: Re1 (ξ1 и ξ2) и Re2 (ξ2 и ξ3), фиксирующих, соответственно, верхнюю (число Re1) и нижнюю (число Re2) границы фактической переходной зоны прогнозируемой расходной характеристики натурной жидкости с вязкостью ν; для последней задают ряд значений расхода Qk=Qj, скоростей истечения vн,k=vн,j из сливного наконечника и потерь напора на его выходе - hн,k=hн,j, совпадающих с полученными на воде значениями; для каждого значения k=j и для каждой из выделенных зон режимов истечения натурной жидкости вычисляют соответствующие этим зонам числа Рейнольдса: Rek,ssvн,k2/2g, при этом числа Rek,1>Re1 относятся к зоне турбулентного режима истечения, числа Rek,3<Re2 - к зоне ламинарного режима, а числа Re2<Rek,2<Re1 - к фактической переходной зоне прогнозируемой расходной характеристики; далее по полученным ранее уравнениям ξs=fs(Re) аппроксимирующих функций вычисляют суммарные коэффициенты потерь напора на местных сопротивлениях - ξk=f(Rek), скорости течения - vi,k=vн,k(dн/di)2, числа Рейнольдса - Rei„k=vi,kdi/ν, коэффициенты трения - λi,k=64/Rei,k и потери напора на трение - hi,ki,k(Li/di)(dн/di)4 на отдельных участках трубопровода, суммарные потери напора на трение - hт,k[см]=Σi[hi,k] и на местных сопротивлениях в тракте - hм,kk(vн,k2/2g) и, наконец, - расчетные, для каждой из выделенных зон, значения располагаемого напора - Hpk=hн,k+hт,k+hм,k.

Значения располагаемого напора Hj ступенчато уменьшают на 10 см вод. ст., начиная с максимального значения Hmax=Hj=1=184 см, до Hj=18=34 см и далее - ступенчато уменьшают на 2 см вод. ст. вплоть до нижней границы диапазона измерений - Hmin31=4 см.

Операция деления полученной на воде экспериментальной расходной характеристики тракта Qj=F(Hj) на зоны реализуется в программе «Advanced Grapher» посредством последовательного сравнительного анализа числовых значений стандартного отклонения точек Qj=p=10, 9, … 1 функции Oj=F(Hj) в верхней области диапазона измерений и точек Qj=q=22, 23, …31 - в нижней области диапазона измерений от ее аппроксимирующей характеристики вида Q=aHb, где а и b - коэффициенты аппроксимации; сравнительный анализ значений стандартного отклонения производится путем последовательного уменьшения числа р на 1 в верхней области диапазона измерений и последовательного увеличения числа q на 1 в нижней области диапазона измерений; результаты измерений расхода Oj=F(Hj) с числами j=p=1, 2…6 и с числами j=q=26, 27…31 принимают за представительские группы расходных характеристик турбулентной и ламинарной зон, соответственно, и в анализе значений величин стандартного отклонения эти результаты измерений не участвуют; зоны турбулентного и ламинарного режимов истечения определяют количествами измерений в верхней и нижней областях диапазона измерений, соответственно, для которых имеют место минимальные значения стандартного отклонения, причем группа результатов измерений, занимающих промежуточное место между результатами измерений с минимальными значениями стандартного отклонения для верхней и нижней областей диапазона измерений, образуют зону условно переходного режима истечения.

По полученным в табличном виде для каждой из выделенных зон функциям Qs,k=F(Hs,pk) путем их аппроксимации степенными функциями вида Qs=c⋅Hsd, где с и d - коэффициенты аппроксимации, определяют искомые для каждой из выделенных зон расходные характеристики натурной жидкости.

На Фиг. 1 представлена функциональная схема измерительной установки для экспериментального определения функций ξs=f(Res). Установка представляет собой систему автоматизированного дозирования (САД) воды, аналогичную системе дозирования, подробное описание которой представлено в патенте-прототипе как устройства для реализации предложенного там способа определения расходной характеристики (РХ) гидравлического тракта при турбулентном режиме истечения.

На Фиг. 1 представлены: 1 - дозирующее устройство (ДУ); 2 - устройство управления (УУ); 3 - расходная емкость с водопроводной водой; 4 - расходный резервуар (РР) из нержавеющей стали полезным объемом 30 л (диаметр резервуара - 20 см, высота - 100 см) с герметичной крышкой; 5 - уровнемерная трубка; 6 - барботажная трубка; 7 - питающая трубка; 8 - трубопровод линии пополнения РР; 9 - запорный вентиль на линии пополнения РР; 10 - мерная емкость (ME); 11 - запорное устройство - пневмоклапан (шаровой кран с двухсторонним пневмоприводом); 12 - сливной наконечник; 13 - гибкий трубопровод на линии налива воды в ME; 14 - струйный эжектор; 15-17 - запорные вентили; 18 - регулятор давления П; 19 - пневмотумблер «КонтрольПз/КонтрольП»; 20 - пневмотумблер «Давление в РР»; 21 - задатчик давления Пз; 22, 23 - пневмореле; 24, 25 - регулируемые дроссели; 26 - пятимембранный элемент сравнения; 27 - блок усилителей мощности; 28 - сборный коллектор; 29 - образцовый манометр; 30 - пневмоэлектронное временное устройство; 31 - электронный таймер; 32 - электропневмопреобразователь; 33 - логический элемент «НЕ».

Испытуемый гидравлический тракт состоит из стальной питающей трубки (ПТ) 7, гибкого трубопровода 13 из ПВХ, запорного устройства 11 - двухпозиционного пневмоклапана в виде шарового крана с двухсторонним пневмоприводом, снабженного сливным наконечником (СН) 12.

РХ САД представляет собой функциональную зависимость объемного расхода Q на выходе СН 12 от величины располагаемо напора H=(p+γвH1)/γв-H2, где р - избыточное давление в РР 4, γв=1 г/см3 - удельный вес воды, H1 - текущее положение уровня жидкости в РР, Н2=66 см - высота установки СН относительно нижней торцевой кромки ПТ 7, p+γвH1=П - полное давление в РР, измеряемое БТ 6.

В состав УУ 2 входят пневматический пропорциональный регулятор давления П (РД) 18 и пневмоэлектронное временное устройство (ВУ) 30.

РД 18 содержит пневмотумблер 19 «КонтрольПз/КонтрольП» и пневмореле 20 для контроля по манометру 29 заданного (Пз) (при выключенном положении тумблера 19) или фактического (П) (при включенном положении тумблера 19) значений полного давления в РР, пневмотумблер 20 «Давление в РР» для включения подачи сжатого воздуха в РР, задатчик 23 давления Пз, пневмореле 22 и 23, регулируемые дроссели 24 и 25, пятимембранный элемент сравнения 26, и блок усилителей мощности 27. Пневмотумблеры 19 и 20 и задатчик 21 функционально образуют пульт управления. Функции собственно регулятора давления выполняет пятимембранный элемент сравнения 26.

РД 18 обеспечивает автоматическое поддержание заданной постоянной величины давления Пj=p+γвH1з,j=const в каждом j-ом цикле проливки (дозирования) путем подкачки сжатого воздуха в РР от усилителей мощности 27 через сборный коллектор 28 и запорный вентиль 16, обеспечивая постоянство величины располагаемого напора H,jз,jв-H2=const, а следовательно, и расхода Qj при открытом положении пневмоклапана 11.

ВУ 30 содержит электронный таймер 31 с цифровой индикацией времени дозирования, устанавливаемого посредством набора кнопками «SET», «↑» и «↓», электропневмопреобразователь (ЭПП) 32 и логический элемент «НЕ» 33. Запуск операции дозирования производится по команде от сенсорной кнопки «СТАРТ», при касании которой вырабатывается управляющий сигнал (24 в) на переключение пневмоконтактов ЭПП 32 и перевод пневмоклапана 11 в открытое положение на время дозирования воды в ME 10.

Для выявления закономерностей зависимости параметра ξ от величин конструктивных параметров гидравлического тракта ДУ измерения расхода воды проводились с использованием в качестве трубопровода 13 гладкостенных гибких трубопроводов из ПВХ трех типоразмеров внутреннего диаметра: d1=0,6; 0,8 и 1,0 см и длины L1=329; 433 и 277 см. Типоразмеры питающей трубки 7 (d2=1,5 см и L2=112 см) во всех циклах измерений оставались неизменными. В качестве СН 12 использовались металлические трубки общей длиной 20 мм с зенковкой входной кромки внутренним диаметром выходного тонкостенного канала dн=0,2; 0,3; 0,4; 0,5 и 0,6 см.

Рабочий цикл измерений представлял собой последовательные замеры расхода воды в гидравлическом тракте ДУ 1 для каждой (из 3-х) выбранной пары типоразмеров [d1_ L1] гибкого трубопровода и сливного насадка (dн) весовым методом. При этом каждый из 15 (3×5) рабочих циклов измерений заключался в последовательных проливках тракта водопроводной водой при различных постоянных значениях располагаемого напора Hjз,jв2=const, где Пз,j - заданная величина полного давления Пj (j=1, 2…31); при этом H1=Hmax=184 см, Н31=Hmin=4 см. Значения давлений Пз,j устанавливали задатчиком давления 21 по образцовому манометру 29 (пределы измерения давления - 0-1 кг/см2; 1 дел. шкалы - 2 см вод. ст.; класс точности - 0,16) с дискретностью 10 см вод. ст. при Пз,j=250-100 см вод. ст. и далее (при Пз,j=100-70 см вод. ст.) - с дискретностью 2 см вод. ст.

В результате проведения операций проливки для каждой пары [d1_L1]_dн типоразмеров конструктивных параметров гидравлического тракта ДУ (например, [0,2_329]_0,6), в соответствии с описанным в п. 1 Формулы изобретения способом, последовательно вычислялись объемный расход - Qj [см3/c]=Vj/tj, соответствующий напору Hj (где Vj [см3] - регистрируемый с помощью взвешивания на электронных весах объем дозы, tj[c] - время дозирования); скорости истечения воды из сливного наконечника - vн,j[см/c]=Qj/fн, потери напора - hн,j [см]=1,05⋅vн,j2/2g и числа Рейнольдса потока - Rej=vн,jdнв на его выходе; скорости течения - vi,j[см/c]=vн,j(dн/di)2, числа Рейнольдса - Rei,j=vi,jdiв, коэффициенты трения - λi,j=64/Rei,j и потери напора на трение - hi,j [см]=λi,j(li/di)(dн/di)4 на отдельных участках трубопровода; суммарные потери напора на трение - hт,j [см]=Σi[hi,j], на местных сопротивлениях в тракте - hм,j[см]=Hj-hн,j-hт,j и суммарные коэффициенты потерь напора на местных сопротивлениях - ξj=2ghм,j/vн,j2.

На последующих этапах обработки полученных экспериментальных данных, в соответствии с предлагаемым способом, с помощью программы «Advanced Grapher» полученную в табличном виде экспериментальную РХ тракта, Oj=F(Hj), по процедуре, описанной в п. 3 Формулы изобретения делили на зоны турбулентного, условно переходного и ламинарного режимов истечения. Далее, для каждой выделенной зоны строились графики функций ξj,s=f(Rej,s) (где индексы s=1, 2, 3 относятся к турбулентной, условно переходной и ламинарной зонам, соответственно), проводился их регрессионный анализ, и для каждой из выделенных зон определялись уравнения аппроксимирующих функций: ξ1=f1(Re), ξ2=f2(Re), ξ3=f3(Re), строились их графики и определялись точки Re1 и Re2 пересечения графиков: Re1 (ξ1 и ξ2) и Re2 (ξ2 и ξ3), фиксирующих, соответственно, верхнюю (число Re1) и нижнюю (число Re2) границы фактической переходной зоны прогнозируемой РХ натурной жидкости с вязкостью v. Для последней задавались значения расхода Qk=Oj, скорости истечения vн,k=vн,j из СН и потери напора на его выходе - hн,k=hн,j, совпадающие с полученными на воде значениями; для каждого значения k=j и для каждой из выделенных зон режимов истечения натурной жидкости вычислялись соответствующие этим зонам числа Рейнольдса: Rek,ssvн,k2/2g, и по полученным уравнениям ξs=fs(Re) аппроксимирующих функций вычислялись суммарные коэффициенты потерь напора на местных сопротивлениях - ξk=f(Rek), скорости течения - vi,k=vн,k(dн/di)2, числа Рейнольдса - Rei,k=vi,kdi/ν, коэффициенты трения - λi,k=64/Rei,k и потери напора на трение - hi,ki,k(Li/di)(dн/di)4 на отдельных участках трубопровода, суммарные потери напора на трение - hт,k [см]=Σi[hi,k] и на местных сопротивлениях в тракте - hм,kk(vн,k2/2g) и, наконец, - расчетные (прогнозируемые для каждой из выделенных зон) значения располагаемого напора - Hpk=hн,k+hт,k+hм,k.

Окончательно, по полученным в табличном виде для каждой из выделенных зон функциям Qs,k=F(Hs,pk) путем их аппроксимации степенными функциями вида Qs=c⋅Hsd определялись искомые (для каждой из выделенных зон) расходные характеристики натурной жидкости.

Для экспериментальной проверки предлагаемого способа (вычисления и анализа погрешностей прогнозирования РХ) циклы «проливки» гидравлического тракта ДУ (фиг. 1) и измерения расхода воды проводились с использованием указанных выше трех типоразмеров внутреннего диаметра (d1) гибкого трубопровода и пяти типоразмеров (dн) внутреннего диаметра СН при различной (условно «высокой» и условно «низкой») сезонной температуре, т.е. при ее различной вязкости. Наибольшая разность температуры воды, при этом, составила ≈14°С. В результате этого этапа исследования было получено 15 (3×5) массивов данных измерений с условно «высокой» («+») и столько же массивов данных измерений с условно «низкой» («-») температурой воды. Т.е., для каждого из трех типоразмеров [d1_ L1] гибкого трубопровода было получено 5 пар результатов измерений (для различных типоразмеров СН) с условно «высокой» («+») и условно «низкой» («-») температурой воды. При этом (для экспериментальной проверки предлагаемого способа) прогнозировались расчетные РХ - Qpk=F(Hpk), которые получились бы при условно «низкой» температуре воды, по экспериментальным и расчетным данным, полученным при ее условно «высокой» температуре. О достоверности получаемых результатов прогнозирования РХ можно было судить по результатам вычисления и анализа относительных погрешностей расхода Qpk при их сравнении с соответствующими экспериментальными значениями Qэk, полученными при условно «низкой» температуре воды.

Эти погрешности имеют максимальные значения в верхней и нижней точках диапазона изменения располагаемого напора, т.е. при Н=Hmax=184 см и H=Hmin=4 см. Их вычисление проводится по формулам:

δQmax[%]=100⋅[Qp(Hmax)-Qэ(Hmax)]/Qэ(Hmax)];

δQmin[%]=100⋅[QP(Hmin)-Qэ(Hmin)]/Qэ(Hmin)],

где Qp(Hmax) и Qp(Hmin) - расчетные, a Qэ(Hmax) и Qэ(Hmin) - экспериментальные значения расхода при Н=184 см и Н=4 см, соответственно. Эти значения определялись посредством использования операции «трассировки» в программе «Advanced Grapher» соответствующих графиков степенных аппроксимирующих функций: Q1,p=cHd и Q1,э=cHd - при вычислении δQmax и Q3,p=cHd, Q3,э=cHd - При вычислении δQmin.

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

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

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

1. Способ определения расходной характеристики гидравлического тракта, характеризующийся тем, что для области перехода от турбулентного к ламинарному режиму истечения модель тракта с проектными конструктивами местных сопротивлений и участков соединительных трубопроводов длиной li≤Li с внутренним диаметром di с концевой раздачей через запорное устройство с цилиндрическим сливным наконечником с внутренним диаметром dн проливают водопроводной водой посредством ее многократного порционного дозирования в мерную тару при различных величинах действующего напора Hj, ступенчато изменяемых от максимального Hj,max, соответствующего турбулентному режиму истечения, до минимального Hj,min, соответствующего началу перехода от струйного к капельно-струйному режиму истечения, значений, которые во время каждой j-й операции проливки поддерживают на заданном постоянном уровне, где j=1, 2…, n=31 - число проливок в каждом цикле измерений, и оперативно, с помощью электронных таблиц программы «Microsoft Excel», для каждой j-й операции проливки последовательно вычисляют объемный расход - Qj [см3/c]=Vj/tj, соответствующий напору Hj, где объем дозы Vj [см3] регистрируют с помощью взвешивания на электронных весах, а время дозирования tj[c] задают и отмеривают с помощью электронного таймера; скорость истечения воды из сливного наконечника - vн,j[см/c]=Qj/fн, потерю напора - hн,j [см]=α⋅vн,j2/2g и число Рейнольдса потока - Rej=vн,jdнв на его выходе, где коэффициент Кориолиса α принимают равным 1,05, как при турбулентном режиме истечения, а кинематическую вязкость воды νв рассчитывают по ее температуре в данном цикле проливок; скорости течения - vi,j[см/c]=vн,j(dн/di)2, числа Рейнольдса - Rei,j=vi,jdiв, коэффициенты трения - λi,j=64/Rei,j и потери напора на трение - hi,j [см]=λi,j(li/di)(dн/di)4 на отдельных участках трубопровода; суммарные потери напора на трение - hт,j [см]=Σi[hi,j], на местных сопротивлениях в тракте - hм,j[см]=Hj-hн,j-hт,j и суммарные коэффициенты потерь напора на местных сопротивлениях - ξj=2ghм,j/vн,j2; при этом полученную в табличном виде экспериментальную расходную характеристику тракта Qj=F(Hj) делят на зоны турбулентного, условно переходного и ламинарного режимов истечения; в программе «Advanced Grapher» для каждой выделенной зоны строят графики функций ξj,s=f(Rej,s), где индексы s=1, 2, 3 относятся к турбулентной, условно переходной и ламинарной зонам, соответственно, проводят их регрессионный анализ и получают для каждой из выделенных зон уравнения аппроксимирующих функций: ξ1=f1(Re), ξ2=f2(Re), ξ3=f3(Re), при этом для каждого s выбирают степенную или гиперболическую, по минимальному значению стандартного отклонения, аппроксимацию; строят графики аппроксимирующих функций и определяют точки Re1 и Re2 пересечения графиков: Re1 (ξ1 и ξ2) и Re2 (ξ2 и ξ3), фиксирующих, соответственно, верхнюю (число Re1) и нижнюю (число Re2) границы фактической переходной зоны прогнозируемой расходной характеристики натурной жидкости с вязкостью ν; для последней задают ряд значений расхода Qk=Qj, скоростей истечения vн,k=vн,j из сливного наконечника и потерь напора на его выходе - hн,k=hн,j, совпадающих с полученными на воде значениями; для каждого значения k=j и для каждой из выделенных зон режимов истечения натурной жидкости вычисляют соответствующие этим зонам числа Рейнольдса: Rek,ssvн,k2/2g, при этом числа Rek,1>Re1 относятся к зоне турбулентного режима истечения, числа Rek,3<Re2 - к зоне ламинарного режима, а числа Re2<Rek,2<Re1 - к фактической переходной зоне прогнозируемой расходной характеристики; далее по полученным ранее уравнениям ξs=fs(Re) аппроксимирующих функций вычисляют суммарные коэффициенты потерь напора на местных сопротивлениях - ξk=f(Rek), скорости течения - vi,k=vн,k(dн/di)2, числа Рейнольдса - Rei„k=vi,kdi/ν, коэффициенты трения - λi,k=64/Rei,k и потери напора на трение - hi,ki,k(Li/di)(dн/di)4 на отдельных участках трубопровода, суммарные потери напора на трение - hт,k[см]=Σi[hi,k] и на местных сопротивлениях в тракте - hм,kk(vн,k2/2g) и, наконец, - расчетные, для каждой из выделенных зон, значения располагаемого напора - Нpk=hн,k+hт,k+hм,k.

2. Способ по п. 1, отличающийся тем, что значения располагаемого напора Hj ступенчато уменьшают на 10 см вод.ст., начиная с максимального значения Hmax=Hj=1=184 см, до Hj=18=34 см и далее - ступенчато уменьшают на 2 см вод.ст. вплоть до нижней границы диапазона измерений - Hmin31=4 см.

3. Способ по п. 1, отличающийся тем, что операция деления полученной на воде экспериментальной расходной характеристики тракта Qj=F(Hj) на зоны реализуется в программе «Advanced Grapher» посредством последовательного сравнительного анализа числовых значений стандартного отклонения точек Qj=p=10, 9, … 1 функции Qj=F(Hj) в верхней области диапазона измерений и точек Qj=q=22, 23, … 31 - в нижней области диапазона измерений от ее аппроксимирующей характеристики вида Q=aHb, где а и b - коэффициенты аппроксимации; сравнительный анализ значений стандартного отклонения производится путем последовательного уменьшения числа р на 1 в верхней области диапазона измерений и последовательного увеличения числа q на 1 в нижней области диапазона измерений; результаты измерений расхода Qj=F(Hj) с числами j=p=1, 2…6 и с числами j=q=26, 27…31 принимают за представительские группы расходных характеристик турбулентной и ламинарной зон, соответственно, и в анализе значений величин стандартного отклонения эти результаты измерений не участвуют; зоны турбулентного и ламинарного режимов истечения определяют количествами измерений в верхней и нижней областях диапазона измерений, соответственно, для которых имеют место минимальные значения стандартного отклонения, причем группа результатов измерений, занимающих промежуточное место между результатами измерений с минимальными значениями стандартного отклонения для верхней и нижней областей диапазона измерений, образуют зону условно переходного режима истечения.

4. Способ по п. 1, отличающийся тем, что по полученным в табличном виде для каждой из выделенных зон функциям Qs,k=F(Hs,pk) путем их аппроксимации степенными функциями вида Qs=c⋅Hsd, где c и d - коэффициенты аппроксимации, определяют искомые для каждой из выделенных зон расходные характеристики натурной жидкости.



 

Похожие патенты:

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

Устройство для моделирования содержит плоский прозрачный контейнер с узкой внутренней полостью, образованной двумя параллельными стенками контейнера, герметично соединенными с боковыми стенками, частично заполненной испытуемой жидкостью, поворотное основание для установки контейнера и механизм разворота основания относительно горизонтальной плоскости.
Предлагаемое изобретение относится к машиностроению, а именно к области испытания различных изделий, работающих при высоком внутреннем давлении. Перед испытанием изделие не менее чем на 90-95% заполняют частицами твердого сыпучего материала с малой сжимаемостью.

Стенд относится к оборудованию для гидравлических испытаний давлением преимущественно двухгорловых баллонов вместимостью до 210 литров на рабочее давление до 39,2 МПа.

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

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

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

Изобретение относится к области экспериментальной техники и может быть использовано для опытного определения динамических характеристик моделей подводных аппаратов.

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

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

Способ управления потоком текучей среды через клапан (100), расположенный в системе (400) для текучей среды, которая имеет расположенный в ней датчик (103) потока. Измеряют с помощью датчика (103) потока расход Qм текучей среды и сравнивают его с предварительно заданным пороговым значением Qт расхода текучей среды.

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

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

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

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

Изобретение относится к средствам для испытания или калибровки аппаратуры для измерения уровня и может использоваться в измерительной технике для поверки, в частности, радарных (радиолокационных) уровнемеров.

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

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

Изобретение может быть использовано для поверки средств измерений объема или массы жидкости в условиях эксплуатации. Передвижная поверочная установка содержит установленные на опорно-рамной конструкции эталонный мерник, тензометрические датчики и электронасос.

Изобретение относится к области измерительной техники и может быть использовано для расчета пропускной способности проектируемых гидравлических трактов транспортных и дозирующих систем в химической, нефтехимической, авиационной, текстильной, лакокрасочной и других отраслях промышленности, в частности, - узлов транспортирования и дозирования клеевых составов при сборке малогабаритных изделий. Способ характеризуется тем, что для области перехода от турбулентного к ламинарному режиму истечения модель тракта с проектными конструктивами местных сопротивлений и участков соединительных трубопроводов длиной li≤Li с внутренним диаметром di с концевой раздачей через запорное устройство с цилиндрическим сливным наконечником с внутренним диаметром dн проливают водопроводной водой посредством ее многократного порционного дозирования в мерную тару при различных величинах действующего напора Hj, ступенчато изменяемых от максимального Hj,max, соответствующего турбулентному режиму истечения, до минимального Hj,min, соответствующего началу перехода от струйного к капельно-струйному режиму истечения, значений, которые во время каждой j-й операции проливки поддерживают на заданном постоянном уровне, где j1, 2 …, n31 - число проливок в каждом цикле измерений, и оперативно, с помощью электронных таблиц программы «Microsoft Excel», для каждой j-й операции проливки последовательно вычисляют объемный расход - Qj [см3c]Vjtj, соответствующий напору Hj, где объем дозы Vj [см3] регистрируют с помощью взвешивания на электронных весах, а время дозирования tj[c] задают и отмеривают с помощью электронного таймера; скорость истечения воды из сливного наконечника - vн,j[смc]Qjfн, потерю напора - hн,j [см]α⋅vн,j22g и число Рейнольдса потока - Rejvн,jdнνв на его выходе, где коэффициент Кориолиса α принимают равным 1,05, как при турбулентном режиме истечения, а кинематическую вязкость воды νв рассчитывают по ее температуре в данном цикле проливок; скорости течения - vi,j[смc]vн,j2, числа Рейнольдса - Rei,jvi,jdiνв, коэффициенты трения - λi,j64Rei,j и потери напора на трение - hi,j [см]λi,j4 на отдельных участках трубопровода; суммарные потери напора на трение - hт,j [см]Σi[hi,j], на местных сопротивлениях в тракте - hм,j[см]Нj-hн,j-hт,j и суммарные коэффициенты потерь напора на местных сопротивлениях - ξj2ghм,jvн,j2. Далее полученную в табличном виде экспериментальную расходную характеристику тракта QjF делят на зоны турбулентного, условно переходного и ламинарного режимов истечения; в программе «Advanced Grapher» для каждой выделенной зоны строят графики функций ξj,sf, где индексы s1, 2, 3 относятся к турбулентной, условно переходной и ламинарной зонам, соответственно, проводят их регрессионный анализ и получают для каждой из выделенных зон уравнения аппроксимирующих функций: ξ1f1, ξ2f2, ξ3f3, при этом для каждого s выбирают степенную или гиперболическую, по минимальному значению стандартного отклонения, аппроксимацию; строят графики аппроксимирующих функций и определяют точки Re1 и Re2 пересечения графиков: Re1 и Re2, фиксирующих, соответственно, верхнюю и нижнюю границы фактической переходной зоны прогнозируемой расходной характеристики натурной жидкости с вязкостью ν; для последней задают ряд значений расхода QkQj, скоростей истечения vн,kvн,j из сливного наконечника и потерь напора на его выходе - hн,khн,j, совпадающих с полученными на воде значениями; для каждого значения kj и для каждой из выделенных зон режимов истечения натурной жидкости вычисляют соответствующие этим зонам числа Рейнольдса: Rek,sξsvн,k22g, при этом числа Rek,1>Re1 относятся к зоне турбулентного режима истечения, числа Rek,3<Re2 - к зоне ламинарного режима, а числа Re2<Rek,2<Re1 - к фактической переходной зоне прогнозируемой расходной характеристики; далее по полученным ранее уравнениям ξsfs аппроксимирующих функций вычисляют суммарные коэффициенты потерь напора на местных сопротивлениях ξkf, скорости течения - vi,kvн,k2, числа Рейнольдса - Rei„kvi,kdiν, коэффициенты трения - λi,k64Rei,k и потери напора на трение - hi,kλi,k4 на отдельных участках трубопровода, суммарные потери напора на трение - hт,k[см]Σi[hi,k] и на местных сопротивлениях в тракте - hм,kξk и, наконец, - расчетные, для каждой из выделенных зон, значения располагаемого напора - Нрkhн,k+hт,k+hм,k. 3 з.п. ф-лы, 1 табл., 1 ил.

Наверх