БАЛЛИСТИЧЕСКИЙ СПОСОБ ОПРЕДЕЛЕНИЯ КОЭФФИЦИЕНТА ЛОБОВОГО СОПРОТИВЛЕНИЯ ТЕЛА Российский патент 2000 года по МПК G01M9/00 

Описание патента на изобретение RU2145067C1

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

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

где сx - коэффициент лобового сопротивления исследуемого тела;
X - сила лобового сопротивления;
ρ - плотность набегающего потока газа;
u - скорость набегающего потока газа;
S - площадь миделевого сечения тела,
либо косвенным способом, при котором измеряют давления на поверхности тела и по распределению давления рассчитывают коэффициент лобового сопротивления [1, 2].

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

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

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

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

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

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

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

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

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


где xи1, xи2 - величины перемещения характерных точек исследуемого тела относительно неподвижной системы координат за последовательные дискретные промежутки времени t1 и t2 соответственно;
xк1, xк2 - величины перемещения характерных точек контрольного тела относительно неподвижной системы координат за последовательные дискретные промежутки времени t1 и t2 соответственно;
с - коэффициент лобового сопротивления исследуемого тела;
с - коэффициент лобового сопротивления контрольного тела;
- коэффициент подобия исследуемого тела,
коэффициент подобия контрольного тела,
ρ - плотность газовой среды,
Sи, Sк - площадь миделевого сечения исследуемого и контрольного тела и соответственно,
mи, mк - масса исследуемого и контрольного тела соответственно.

Полученное значение коэффициента лобового сопротивления исследуемого тела относят к скорости

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

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

Сущность изобретения поясняется на фигурах 1 - 3, где на фиг.1 дана схема реализации предлагаемого способа на измерительном участке траектории, а на фиг. 2 и 3 приведен пример реализации способа.

Схема реализации предлагаемого способа на измерительном участке траектории приведена на фиг. 1. Oxy - неподвижная система координат; 1 - исследуемое тело; 2 - контрольное тело; 3, 4, 5 - соответственно начальное, первое и второе сечения регистрации координат исследуемого и контрольного тел в дискретные промежутки времени.

Способ реализуется следующим образом. Исследуемое тело 1 метают с помощью метательной установки. Одновременно с метанием исследуемого тела 1 производится параллельное метание контрольного тела 2 с той же начальной скоростью V0 начальном сечении регистрации 3, что и у исследуемого тела 1. На измерительном участке траектории производят одновременно регистрацию координат их характерных точек в начальном 3, первом 4 и втором 5 сечениях регистрации относительно неподвижной системы координат Oxy в дискретные моменты времени t1 и t2, например, путем фотографирования. Здесь t1 - время между моментами регистрации координат в начальном 3 и первом 4 сечениях регистрации, t2 - время между моментами регистрации координат в начальном 3 и втором 5 сечениях регистрации. Без уменьшения общности считают, что в начальном сечении регистрации 3 t0 = 0, x0 = 0. Поэтому определяют величины перемещения характерных точек исследуемого 1 и контрольного 2 тел между начальным 3 и первым 4 сечениями регистрации равные xи1 и xк1 соответственно за время t1 и величины перемещения характерных точек исследуемого 1 и контрольного 2 тел между начальным 3 и вторым 5 сечениями регистрации, равные xи2 и xк2 соответственно за время t2.

Коэффициенты лобового сопротивления с исследуемого 1 и с контрольного 2 тел полагают постоянными на измерительном участке траектории [3]. Это предположение выполняется достаточно точно, так как измерительный участок баллистической трассы имеет небольшую протяженность, порядка нескольких метров, замедление тел на нем невелико и, следовательно, мало падение скорости на измерительном участке. Учитывая сказанное, дважды интегрируют уравнение движения тела

где m - масса тела;
сx - коэффициент лобового сопротивления тела;
S - площадь миделевого сечения тела;
V - скорость движения тела;
ρ - плотность газовой среды;
t - время движения тела при начальных условиях: t0 = 0, x0 = О, V = V0.

В результате получают зависимость для перемещения тела

или
xbcx = ln(bcxV0t+1). (6)
Запишем уравнение (6) для перемещений исследуемого тела - xи1 и контрольного тела xк1, за промежуток времени t1 между начальным 3 и первым 4 сечениями регистрации
xи1bиc = ln(bиcV0t1 + 1)
xк1bкc = ln(bкcV0t + 1). (7)
Система уравнений (7) соответствует сущности изобретения, когда одновременно с метанием исследуемого тела 1 производят параллельное метание с той же начальной скоростью V0 контрольного тела 2 и регистрируют величину их перемещений за одно и то же время t1.

Потенцируя систему (7), после преобразования получают следующее уравнение:

Аналогично (7) составляют систему уравнений для перемещений исследуемого тела xи2 и контрольного тела xк2 за время t2 между начальным 3 и вторым 5 сечениями регистрации, которая после преобразований имеет вид

Уравнения (8) и (9) образуют систему двух уравнений с двумя неизвестными коэффициентами лобового сопротивления исследуемого тела с и контрольного тела ски.

Величины, входящие в систему уравнений (8), (9) p, mи, mк, Sи, Sк определяются до опыта. Величины перемещений характерных точек исследуемого тела 1 и контрольного тела 2 xи1, xи2, xк1, xк2 за дискретные промежутки времени t1, и t2 получают в опыте с точностью, определяемой разрешающей способностью фотоматериалов и юстировкой баллистической трассы [7]. При этом методическая ошибка предлагаемого способа обусловлена предположением равенства начальной скорости V0 исследуемого тела 1 и контрольного тела 2 в начальном сечении регистрации и предположением о постоянстве коэффициентов с и с на измерительном участке траектории.

В уравнения (8) и (9) не входят значения начальной скорости V0 и времени t1 и t2. Значение t1 необходимо только для расчета скорости V0 по (2).

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

Реализация предлагаемого способа, когда в качестве контрольного тела выбирают тело с произвольной аэродинамической формы описана выше. В случае, когда в качестве контрольного тела выбирают тело, имеющее ту же аэродинамическую форму, что и исследуемое тело, уменьшается число сечений регистрации до двух, отпадает необходимость в определении величин перемещений характерных точек исследуемого и контрольного тел xи2, xк2. При этом повышается точность определения с, а система уравнений (8) и (9) преобразуется в уравнение (3).

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

Оценку относительной ошибки определения коэффициента лобового сопротивления исследуемого тела с производят традиционным способом [8], дифференцируя (1), (3) или (4). Тогда, например, для (4) получают
δc = δρ+2(δxи1+δmи+δsи), (10)
где через δ обозначены относительные погрешности измерении исходных параметров, а множитель 2 перед скобкой указывает на то, что производится метание двух тел.

Предлагаемый способ был реализован на баллистическом стенде. В качестве исследуемого тела было выбрано тело из алюминия сферической формы с Sи = 12 мм2 имело массу mи = 4.05 г. Контрольное тело имело ту же аэродинамическую форму с Sк = 11,5 мм2, но было изготовлено из стали и имело массу mк = 11.7 г. Метание исследуемого 1 и контрольного 2 тел производилось одновременно с помощью специального поддона, что позволяло получать в начальном сечении регистрации 3 одинаковые начальные скорости для исследуемого и контрольного тел. В опыте расстояние между начальным 3 и первым 4 сечениями регистрации составляло 150 см. Определение координат характерных точек тел, расположенных на сферической части донного среза, производилось в начальном 3 (фиг. 2) и первом 4 (фиг. З) сечениях регистрации с помощью последовательного фотографирования в дискретные моменты времени. По полученным значениям координат определялись величины перемещения характерных точек исследуемого и контрольного тел, которые являлись исходными данными для расчета коэффициента лобового сопротивления исследуемого тела. Из соотношения (3) было получено с = с = 0,98. Рассчитанная относительная ошибка определения коэффициента лобового сопротивления составила 1,6%. Определение коэффициента лобового сопротивления по способу, выбранному в качестве прототипа, позволяет определить искомый коэффициент с точностью около 5%.

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

Источники информации
1. Прикладная аэродинамика. /Под редакцией Н.Ф.Краснова. - М.: Высшая школа. - 1974. - c. 7-23, c. 250-268.

2. Горлин С. М. Экспериментальная аэродинамика. - М.: Высшая школа. - 1970. - c. 189- 261.

3. Баллистические установки и их применение в экспериментальных исследованиях. / Под ред. Н.А. Златина и Г.И. Мишина. - М.: Наука. - 1974. - с. 7-11, с. 241-249, с. 278-279, с. 285- 289.

4. Физико-газодинамические баллистические исследования. /Под редакцией Г.И.Мишина. - Л.: Наука. - 1980. - с.9-24.

5. Аэрофизические исследования сверхзвуковых течений. /Под редакцией Ю. А.Дунаева. М. - Л.: Наука. - 1967. - с. 163-169, с. 190-192.

6. Дмитриевский А.А. Внешняя баллистика. - М.: Машиностроение. - 1972. - с. 541-543.

7. Оптические методы исследований в баллистическом эксперименте. /Под ред. Г.И. Мишина - М.: Наука. - 1979. - с. 21-31.

8. Фихтенгольц Г.М. Основы математического анализа. - М.: Физматгиз. - 1960. - 440 с.

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

название год авторы номер документа
СПОСОБ ОПРЕДЕЛЕНИЯ ИЗМЕНЕНИЯ ДАВЛЕНИЯ АТМОСФЕРЫ ПО ЕЕ ПРОФИЛЮ 2018
  • Казаков Геннадий Викторович
  • Мочалов Виктор Васильевич
  • Пономарёв Владимир Анатольевич
  • Терентьев Олег Сергеевич
  • Кунавин Владимир Иванович
RU2691921C1
Измерительно-регистрирующий комплекс для определения внутрибаллистических и траекторных параметров метаемого тела 2020
  • Липанов Алексей Матвеевич
  • Шрагер Эрнест Рафаилович
  • Ищенко Александр Николаевич
  • Буркин Виктор Владимирович
  • Майстренко Иван Викторович
  • Дьячковский Алексей Сергеевич
  • Корольков Леонид Валерьевич
  • Саммель Антон Юрьевич
RU2752131C1
СПОСОБ ОПРЕДЕЛЕНИЯ ДАВЛЕНИЯ АТМОСФЕРЫ ПО ЕЁ ПРОФИЛЮ 2021
  • Казаков Геннадий Викторович
  • Мочалов Виктор Васильевич
  • Терентьев Олег Сергеевич
  • Алпатов Евгений Вячеславович
  • Кунавин Владимир Иванович
RU2759851C1
СПОСОБ ПОСТРОЕНИЯ ОПТИМАЛЬНОЙ АЭРОДИНАМИЧЕСКОЙ ПОВЕРХНОСТИ ГИПЕРЗВУКОВОГО ЛЕТАТЕЛЬНОГО АППАРАТА 2013
  • Галактионов Алексей Юрьевич
  • Лихачева Екатерина Сергеевна
RU2558498C2
УСТРОЙСТВО ДЛЯ ОПРЕДЕЛЕНИЯ КОЭФФИЦИЕНТА ЛОБОВОГО СОПРОТИВЛЕНИЯ ПЕРЕМЕЩАЮЩЕГОСЯ ТЕЛА 2010
  • Копылов Геннадий Алексеевич
  • Ковалёв Вячеслав Данилович
  • Лопатин Андрей Сергеевич
RU2442172C2
СПОСОБ ОПРЕДЕЛЕНИЯ БАЛЛИСТИЧЕСКОГО КОЭФФИЦИЕНТА ОБЪЕКТА 2015
  • Воропаев Анатолий Павлович
  • Захаров Владимир Николаевич
  • Пономарев Владимир Анатольевич
RU2600515C1
Снаряд для стрельбы в водной среде 2017
  • Архипов Владимир Афанасьевич
  • Коноваленко Алексей Иванович
  • Перфильева Ксения Григорьевна
RU2677506C1
СПОСОБ ПРЕДСКАЗАНИЯ ЗЕМЛЕТРЯСЕНИЙ 2000
  • Давыдов В.Ф.
  • Есаков В.А.
  • Новоселов О.Н.
  • Комаров Е.Г.
  • Галкин Ю.С.
RU2170447C1
СПОСОБ ОПРЕДЕЛЕНИЯ КОЭФФИЦИЕНТА АЭРОДИНАМИЧЕСКОГО СОПРОТИВЛЕНИЯ ДВИЖУЩЕГОСЯ ТЕЛА 2000
  • Павлюк А.С.
  • Павлюк С.А.
  • Баранов А.С.
  • Бизяев С.Н.
RU2172937C1
Способ определения температуропроводности и коэффициента теплопроводности 2022
  • Соколов Анатолий Константинович
  • Якубина Ольга Анатольевна
RU2785084C1

Иллюстрации к изобретению RU 2 145 067 C1

Реферат патента 2000 года БАЛЛИСТИЧЕСКИЙ СПОСОБ ОПРЕДЕЛЕНИЯ КОЭФФИЦИЕНТА ЛОБОВОГО СОПРОТИВЛЕНИЯ ТЕЛА

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

Формула изобретения RU 2 145 067 C1

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


где xu1, xu2 - величины перемещения характерных точек исследуемого тела относительно неподвижной системы координат за последовательные дискретные промежутки времени t1 и t2 соответственно;
xk1, xk2 - величины перемещения характерных точек контрольного тела относительно неподвижной системы координат за последовательные дискретные промежутки времени t1 и t2 соответственно;
cxu - коэффициент лобового сопротивления исследуемого тела;
cxk - коэффициент лобового сопротивления контрольного тела;
коэффициент подобия исследуемого тела;
коэффициент подобия контрольного тела;
ρ - плотность газовой среды;
su, sk - площадь миделевого сечения исследуемого и контрольного тела соответственно;
mu, mk - масса исследуемого и контрольного тела соответственно,
и относят полученные значения коэффициентов лобового сопротивления исследуемого тела и контрольного тела к скорости

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

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

Документы, цитированные в отчете о поиске Патент 2000 года RU2145067C1

Аэрофизические исследования сверхзвуковых течений/Под ред
Ю.А.Дунаева
М.-Л.: Наука, 1967, с.163 - 169, 190 - 192
СПОСОБ ОПРЕДЕЛЕНИЯ КОЭФФИЦИЕНТА ЛОБОВОГО СОПРОТИВЛЕНИЯ ТЕЛ 1991
  • Попов Валентин Николаевич
  • Кривоцюк Виктор Иванович
RU2103666C1
US 4532801 A, 06.09.85
УСТРОЙСТВО ДЛЯ МЕТАНИЯ МОДЕЛЕЙ 0
SU255078A1
US 3456503 A, 22.07.69
УСТРОЙСТВО для ОПРЕДЕЛЕНИЯ СОПРОТИВЛЕНИЯ СРЕДЫ ПРИ ПОЛЕТЕ ТВЕРДЫХ ЧАСТИЦ 0
SU211845A1

RU 2 145 067 C1

Авторы

Якушев В.К.

Биматов В.И.

Николаев Р.Г.

Даты

2000-01-27Публикация

1998-05-08Подача