Момент, що виникає на кормовій частині машини, яка рухається по воді

Обґрунтування та здійснення вибору математичної моделі, яка описує величину гідродинамічного моменту, який діє на кормову частину плаваючої машини. Оцінка точності обчислення значень величини гідродинамічного моменту. Оцінка якості математичної моделі.

Рубрика Производство и технологии
Вид статья
Язык украинский
Дата добавления 30.09.2018
Размер файла 350,2 K

Отправить свою хорошую работу в базу знаний просто. Используйте форму, расположенную ниже

Студенты, аспиранты, молодые ученые, использующие базу знаний в своей учебе и работе, будут вам очень благодарны.

Размещено на http://www.allbest.ru/

Одеський національний морський університет

МОМЕНТ, ЩО ВИНИКАЄ НА КОРМОВІЙ ЧАСТИНІ МАШИНИ, ЯКА РУХАЄТЬСЯ ПО ВОДІ

Бугаєв С.В., к.т.н.

Анотація

У роботі наведено обґрунтування і вибір математичної моделі, яка описує величину гідродинамічного моменту, який діє на кормову частину плаваючої машини.

Annotation

Mathematical model is proposed in this article. It describes the part of swimming machine. Model determination is also grounded.

Постановка проблеми

Визначення величин гідродинамічних моментів, діючих на елементи конструкції плаваючої машини (ПМ), що рухається по воді, - складна науково-технічна проблема. Внаслідок взаємодії ПМ із навколишнім середовищем під час її руху, на границі розділу двох середовищ (повітря і води) виникають аеро- та гідродинамічні моменти, які зі зростанням швидкості руху нелінійно змінюються і обумовлюють зміну диференту ПМ, що призводить до негативних наслідків.

Зв'язок з важливими науковими програмами

Обґрунтування і розрахунок математичної моделі (ММ), що описує зміну величини моменту, який діє на корму ПМ під час руху по воді, є складовою частиною методу розрахунку ходових якостей ПМ з урахуванням впливу диферентуючих сил і моментів. Метод розрахунку є складовою частиною рішення проблеми наукового обґрунтування підвищення експлуатаційних характеристик цього класу машин, що безпосередньо випливає з загального положення Державної програми розвитку транспорту.

Аналіз досліджень і публікацій

Аналіз наукової літератури за даним напрямком [1,2,3] показав, що робіт з визначення гідродинамічного моменту, який діє на кормову частину ПМ при взаємодії її з водою, на теперішній час не існує. Проведені модельні випробування в аеродинамічній трубі Одеського національного морського університету [4] дозволили одержати чисельні значення розподілу тиску, який викликаний впливом потоку повітря, що набігає на елементи конструкції ПМ. Використання результатів [5] і застосування методів теорії подібності і розмірності [5] дало можливість встановити кількісні закономірності зміни величини гідродинамічного моменту, що діє на корму ПМ при русі по воді, від швидкості руху та конструкційних особливостей машини. У статті викладено порядок розрахунку за MM, що дозволяє в залежності від швидкості руху і конструкційних особливостей ПМ визначати величину гідродинамічного моменту, який діє на корму машини.

Постановка завдання

Метою цієї статті є за допомогою положень теорії подібності і розмірності [6], а також методів рішення багатокритеріальних задач [7] одержати й обґрунтувати ММ, що описує зміну величини гідродинамічного моменту, що діє на кормову частину ПМ.

Основний матеріал

Приймемо, що кормова частина корпуса ПМ є сталева вертикальна пластина, встановлена перпендикулярно напрямку руху машини. Обтікання кормової частини корпуса веде до виникнення гідродинамічного моменту (Мк), що впливає на зміну диференту ПМ. На величину Мк, мають вплив наступні чинники: vм - швидкість руху ПМ; Тм- осідання машини у воду; кути нахилу носової частини (бн) і днища корпусу (вд).

В якості чинників, що істотно впливають на значення Мк, приймемо: vм; Тм; бн і вд. У якості функції відклику приймемо - Мк:

Мк = f(vпм; ТМ; бн.вд).

Досліджувана залежність Мк (1) у загальному випадку може бути представленою [7] у вигляді чотирьохмірної поверхні. Діапазон зміни значень прийнятих чинників наведено у табл. 1.

Таблиця 1

Значення чинників, які впливають на момент Мк, що використовуються у ММ

Чинники

vпм(м/с) (Re)

Тм (м)

бд (П)

вд (П)

Кодовані значення

-1

0 (0)

0,75

34

+6 на ніс

0

1,666668 (1,1E-05)

1,1533

40

0

1

3,333336 (2,77E-05)

1,45

46

-6 на ніс

У табл. 2 наведені значення моменту Мк, що діє на корму ПМ, у залежності від швидкості руху (vпм), осідання ПМ (Тм), кутів нахилу носової частини (бн) і днища (вд).

Таблиця 2

Значення моменту (Мк), що діє на кормову частину ПМ, у залежності від: v; бн ;вд.

v

(чинник)

Величина кутів нахилу носової частини (бн) і днища (вд)

Тк=0,75м(-1);

б=46пн

вд=+6д

б=40пн

вд=+6д

б=34пн

вд=+6д

б=46пн

вд=0д

б=40пн вд=0д

б=34пн вд=0д

б=46пн

вд=-6д

б=40пн

вд=-6д

б=34пн

вд=-6д

-1

0

0

0

0

0

0

0

0

0

0

0,0029

0,0039

0,0048

0,0056

0,00648

0,0074

0,0083

0,009

0,012

1

0,025

0,0269

0,0288

0,0306

0,0325

0,0344

0,0363

0,038

0,04

Тк=1,1533м(0);

0

0,003

0,003875

0,00475

0,005625

0,0065

0,007375

0,00825

0,009125

0,01

1

0,02

0,02375

0,0275

0,03125

0,035

0,03875

0,0425

0,04625

0,05

Тк=1,45м (1)

0

0,0014

0,001713

0,002025

0,002338

0,00265

0,002963

0,003275

0,003588

0,0039

1

0,0056

0,006113

0,006625

0,007138

0,00765

0,008163

0,008675

0,009188

0,0097

Для більшої наочності, залежності Мк від v; LПМ; бн і вд наведені на рис.1.

Рис. 1 Значення моменту Мк в залежності від швидкості руху (vпм), кутів нахилу носу (бн) та днища ПМ (вд) при фіксованому значенні:

а) Тк=0,75 м (-1);

б) Тк=1,1533 м (0);

в) Тк=1,45 м (-1).

Таблиця 3

План проведення експерименту щодо визначення величини моменту Мк X1 - Швидкість руху ПМ, (у числах Рейнольдса Re) Х2 - Осідання ПМ (Тм) м Хз - Кут нахилу носової частини ПМ (бн) Х4 - Кут нахилу днища ПМ (вд) п щодо поверхні води п

Досліджувані в ході експерименту чинники

РІВНІ

X1

Х2

Хз

Х4

Основний

9,16x10-6 Re

1,1533 м

40п

0 п

Інтервал варіювання

9,16x10-6 Re

0,4033 м

6 п

6 п

Верхня межа

1,83x10-5 Re

1,45 м

46п

+6 п

Нижня межа

0 Re

0,75 м

34п

-6 п

Кодуємі значення чинників

Результат

№ досліду

Хо

XI

Х2

Хз

Х4

XI

Х2

XI

Хз

Х2

Хз

XI

Х4

X2

Х4

Хз

Х4

XI

Х2

Хз

XI

Х2

Х4

X2

Х3

Х4

XI Х2

Хз Х4

Y (функція віклику Мк)

Ядро

Чинники

Ефекти взаємодії

1

+1

-1

-1

+1

+1

+1

-1

-1

-1

-1

+1

+1

+1

-1

+1

0

2

+1

-1

+1

-1

+1

-1

+1

-1

-1

+1

-1

+1

-1

-1

+1

0

3

+1

+1

-1

-1

+1

-1

-1

+1

+1

-1

-1

+1

+1

+1

+1

0,0288

4

+1

+1

+1

+1

+1

+1

+1

+1

+1

+1

+1

+1

+1

+1

+1

0,007138

5

+1

-1

-1

-1

+1

+1

+1

+1

-1

-1

-1

-1

+1

+1

-1

0

6

+1

-1

+1

+1

+1

-1

-1

+1

-1

+1

+1

-1

-1

+1

-1

0

7

+1

+1

-1

+1

+1

-1

+1

-1

+1

-1

+1

-1

-1

-1

-1

0,025

8

+1

+1

+1

-1

+1

+1

-1

-1

+1

+1

-1

-1

+1

-1

-1

0,006625

9

+1

-1

-1

+1

-1

+1

-1

-1

+1

+1

-1

+1

-1

+1

-1

0

10

+1

-1

+1

-1

-1

-1

+1

-1

+1

-1

+1

+1

-1

+1

-1

0

11

+1

+1

-1

-1

-1

-1

-1

+1

-1

+1

+1

+1

+1

-1

-1

0,04

12

+1

+1

+1

+1

-1

+1

+1

+1

-1

-1

-1

+1

-1

-1

-1

0,008675

13

+1

-1

-1

-1

-1

+1

+1

+1

+1

+1

+1

-1

-1

-1

-1

0

14

+1

-1

+1

+1

-1

-1

-1

+1

+1

-1

-1

-1

+1

-1

+1

0

15

+1

+1

-1

+1

-1

-1

+1

-1

+1

+1

-1

-1

+1

-1

+1

0,0363

16

+1

+1

+1

-1

-1

+1

-1

-1

-1

-1

+1

-1

-1

+1

-1

0,0097

17 Центр

0

0

0

0

0

0

0

0

0

0

0

0

0

0

0

0,0065

Таблиця 4

Коефіцієнти регресійних конкуруючих математичних моделей

Коеф. регр. ММ (1): y = ao +a1x1 +a2x2 +a3x3 +a4x4

похибка значень: 0,133361827; R=,98493 дисперсія результатів: 97,010%

Const.а0

a1

a2

a3

a4

0,009926

0,01014

-0,00612

-0,0005

-0,00169

Коеф. регрес. ММ (2) y = ao +a1x1 +a2x2 +a3x3 +a4x4 +a5x1 x2 +a6 x1 x3 +a7 x1 x4 + a8x2 x3 +a9x2x4 +a10x3x4

похибка значень:,133361827 R=,98493 дисперсія результатів: 97,010%

Const.а0

a1

a2

a3

a4

a5

a6

a7

0,009896

0,010179

-0,00614

-0,0005

-0,0017

-0,00614

-0,0005

-0,0017

a8

a9

a10

0,00044

0,001124

9,18E-05

Коеф. регрес. ММ (3): y = ao +a1x1 +a2x2 +a3x3 +a4x4 +a5x1 x2 +a6 x1 x3 +a7 x1 x4 + a8x2 x3 +a9x2x4 +a10x3x4 +a11 x1x2x3+a12 x1x3x4 + a13x1x3x4+ a14x2 x3x4+ a15 x1x2 x3x4

похибка значень:0,127111765 R=,98565 дисперсія результатів: 97,150%

Const.а0

a1

a2

a3

a4

a5

a6

a7

0,009926

0,01014

-0,00612

-0,0005

-0,00169

-0,00612

-0,0005

-0,00169

a8

a9

a10

a11

a12

a13

a14

a15

0,000437

0,001118

8,99E-05

0,000437

0,001118

8,99E-05

0,000102

0,000102

Коеф. регрес. ММ (4): y = ao +a1x1 +a2x2 +a3x3 +a4x4 +a5x1 x2 +a6 x1 x3 +a7 x1 x4 + a8x2 x3 +a9x2x4 +a10x3x4 +a11 x1x2x3+a12 x1x3x4 + a13x1x3x4+ a14x2 x3x4+ a15 x1x2 x3x4. + a16x12 +a17x22 +a18 x32 + a19x42

похибка значень:0, 000001220; R=1,0000 дисперсія результатів: 100,00%

Const.а0

a1

a2

a3

a4

a5

a6

a7

0,0065

0,010151

-0,00613

-0,0005

-0,00169

-0,00613

-0,0005

-0,0017

a8

a9

a10

a11

a12

a13

a14

a15

0,000439

0,00112

9,07E-05

0,000222

0,000532

5,82E-05

5,42E-05

5,75E-05

a16

a17

a18

a19

0,000455

0,000455

0,000455

0,000455

Таблиця 5

Оцінка точності обчислення значень величини моменту Мк

Точне

значення

№ оцінюваної математичної моделі (див. табл. 4)

ММ №1

ММ №2

ММ №3

ММ №4

Значення

Похибка

Значення

Похибка

Значення

Похибка

Значення

Похибка

Y

Yмм

ДY

Yмм

ДY

Yмм

ДY

Yмм

ДY

0

0,003713

-0,00371

-0,00176

0,001756

0,000214

-0,00021

-8E-05

8,01E-05

0

-0,00753

0,00753

0,000308

-0,00031

0,000214

-0,00021

8,24E-05

-8,2E-05

0,0288

0,024995

0,003805

0,029191

-0,00039

0,000214

0,028586

0,028876

-7,6E-05

0,0071

0,011748

-0,00461

0,005042

0,002096

0,000214

0,006924

0,007134

3,8E-06

0

0,004715

-0,00471

-0,00106

0,001059

0,000214

-0,00021

-5E-05

5,03E-05

0

-0,00853

0,008532

0,001372

-0,00137

0,000214

-0,00021

1,03E-05

-1E-05

0,025

0,023993

0,001007

0,02649

-0,00149

0,000214

0,024786

0,025077

-7,7E-05

0,0066

0,012749

-0,00612

0,005983

0,000642

0,000214

0,006411

0,00651

0,000115

0

0,007102

-0,0071

0,000308

-0,00031

0,000214

-0,00021

6,68E-05

-6,7E-05

0

-0,00414

0,004141

-0,00176

0,001756

0,000214

-0,00021

-0,0001

0,0001

0,04

0,028384

0,011616

0,038419

0,001581

0,000214

0,039786

0,039993

6,99E-06

0,0087

0,015137

-0,00646

0,009408

-0,00073

0,000214

0,008461

0,008689

-1,4E-05

0

0,008104

-0,0081

0,001372

-0,00137

0,000214

-0,00021

1,96E-05

-2E-05

0

-0,00514

0,005143

-0,00106

0,001059

0,000214

-0,00021

-4,4E-05

4,36E-05

0,0363

0,027382

0,008918

0,035351

0,000949

0,000214

0,036086

0,036259

4,13E-05

0,0097

0,016138

-0,00644

0,010716

-0,00102

0,000214

0,009486

0,009787

-8,7E-05

0,0065

0,009926

-0,00343

0,009896

-0,0034

-0,00343

0,009926

0,0065

2,47E-07

Як видно з наведених на рис.1 залежностей, величина моменту Мк визначається складною закономірністю, що обумовлено виникненням інтенсивного вихроутворення за корпусом ПМ. Для попередньої оцінки впливу прийнятих чинників на функцію, що досліджується (Mк), побудовані графіки впливу чинників (рис 2).

Рис 2 Вплив чинників: v, Lпм; бн (на графіку <н); вд (на графіку < д) на функцію відклику: Мк = f(v, Тм; бн; вд)

Вибір регресійної моделі.

Попередній аналіз графічних залежностей, наведених на рис. 1 і 2, показав, що апроксимація може бути виконана повним поліномом [7] при чотирьох аргументах виду:

y = ao +a1x1 +a2x2 +a3x3 +a4x4 … +an-2x34++an-1x44 + anx1 x2 x3 x4

гідродинамічний момент кормовий плаваючий

Такий поліном містить 70 доданків, що створює дуже громіздку модель. Для зменшення кількості доданків застосований метод послідовного ускладнення математичної залежності. План проведення експериментів наведено у табл. 3. В якості регресійної моделі приймалися частинні поліноми від повного (1).

Критерієм оцінки якості ММ приймалась похибка обчислення функції відклику (точність >0,001% у всьому діапазоні значень). Коефіцієнти регресії визначалися ітераційним квазі Ньютонівським методом [8] з критерієм збіжності 0,0001.

Результати обчислень коефіцієнтів регресійних конкуруючих ММ наведено у табл.4.

Оцінка точності обчислення значень моменту Мк, що діє на корму ПМ, проводилась з використанням розглянутих ММ у точках плану експерименту, результати наведено у табл. 5.

Отримані результати показують, що перші три ММ не можуть бути використані для розрахунків (мах. похибка більше, ніж прийнята для оцінки якості ММ).

Обчислення значень моменту Мк для точок плану при використанні ММ №4 дає високу точність (мах по модулю похибка 8,7х10-5).

Для оцінки адекватності прийнятої моделі був побудований графік залишків у логарифмічному масштабі (рис. 3).

Рис. 3 Графік похибок величини Мк, отриманих при застосуванні ММ №4

Висновки

Наведені результати показують, що отримана в ході проведення обчислювального експерименту ММ із достатньою для інженерних задач точністю описує зміну величини Мк, у залежності від швидкості руху ПМ і її конструктивних особливостей. Модель можна рекомендувати для розрахунків, пов'язаних з визначенням ходових якостей ПМ.

Література

1. Степанов А.П. Конструкция и расчет плавающих машин. М.: Машиностроение, 1978. 232. С.

2. Степанов А.П. Методика расчета основных водоходных качеств плавающих машин. Навчальний посібник. М. АБТВ, 1971. 41 С.

3. Лактионов К.Н. Методика оценки сопротивления воды с целью поиска путей его уменьшения. // Сб. научных трудов МАДИ "Оптимизация при расчетах и конструировании специальных транспортных средств". 1986. С. 111-116.

4. Бугаев С.В. Аэродинамические исследования распределения давления по контуру погруженной части плавающей машины.// Вісник ОДМУ. 1999. №4. С. 85-93.

5. Бугаев С.В. Плавающая машина. Теория и эксперимент. Одеса Друк. ОДМУ., 2000. 253 С.

6. Методы подобия и размерности в механике. М.,: Наука, 1967. 148 С.

7. Адлер Ю. П., Маркова Е. В.,. Грановский Ю. В. Планирование эксперимента при поиске оптимальных условий М.: Наука, 1976. 69-154 С.

8. Боровиков В.П., Боровиков И.П. Statistica статистический анализ и обработка данных в среде Windows. М: «Филинъ», 1997. 583 с.

Размещено на Allbest.ru


Подобные документы

  • Обґрунтування найбільш раціонального типу вихоревої турбіни, що відповідає умовам роботи приводу гідродинамічного очисника. Параметри силової взаємодії потоку робочої рідини з лопатками робочого колеса вихоревої турбіни, розробка практичних рекомендацій.

    автореферат [444,2 K], добавлен 26.07.2009

  • Основні причини виникнення похибок. Їх класифікація і принципи оцінювання. Визначення відносної і приведеної похибок. Особливості математичної моделі їх визначення. Правила округлення значень і форми запису кінцевого результату. Критерії оцінки промахів.

    реферат [592,9 K], добавлен 23.08.2013

  • Положення розмірного аналізу конструкції. Основні методичні положення розмірного аналізу машини чи складальної одиниці. Порядок проведення розмірного аналізу конструкції машини чи складальної одиниці. Вибір методу досягнення точності замикальної ланки.

    реферат [448,3 K], добавлен 08.07.2011

  • Аналіз технологічного процесу складання заготовки і устаткування, яке використовується в діючому цеху. Аналіз якості взуття. Обґрунтування вибору моделі відповідно до напряму моди. Обґрунтування способу формування заготовки на колодці і методу кріплення.

    контрольная работа [51,8 K], добавлен 25.03.2014

  • Процеси, що протікають в посудомийних машинах. Шляхи поліпшення якості миття. Пристрої автоматизації миття посуду. Розробка лабораторного стенду для дослідження характеристик посудомийної машини. Опис гідравлічної принципової схеми, порядок роботи.

    курсовая работа [721,1 K], добавлен 20.06.2013

  • Розробка моделі зачіски відповідно до історичної епохи, типу обличчя і напрямків моди. Розробка технологічної послідовності виконання зачіски. Обґрунтування вибору, парфумерно-косметичних засобів, інструментів, обладнання, необхідних для виконання моделі.

    дипломная работа [1,2 M], добавлен 26.10.2012

  • Сутність клейового методу кріплення підошви до заготовки. Обґрунтування вибору колодки і матеріалів для взуття. Розмірно-повнотний асортимент для туфель. Проектування моделі методом копіювально-графічної системи. Технологічний процес складання заготовки.

    курсовая работа [412,9 K], добавлен 24.11.2015

  • Машинно-тракторний парк ТОВ "Агрофірма "Маяк". Призначення мельничного комплексу, його технічна характеристика. Будова та опис технологічного процесу млина. Підготовка мельниці до роботи. Призначення і будова оббивальної машини. Розрахунок проводу машини.

    дипломная работа [535,5 K], добавлен 07.06.2012

  • Вибір і обґрунтування моделі повсякденної сукні. Технічне завдання на її розробку. Основні матеріали для її виконання. Особливості розробки комплекту лекал, етапи раціональної технології виготовлення проектної моделі з врахуванням можливостей обладнання.

    курсовая работа [1,4 M], добавлен 19.02.2014

  • Розробка модельного ряду молодіжних жакетів. Обґрунтування вибору методу технічного моделювання та методики конструювання моделі молодіжного жакету. Розкладка деталей крою швейного виробу. Вивчення основних способів з’єднання деталей швейного виробу.

    курсовая работа [1,2 M], добавлен 31.10.2014

Работы в архивах красиво оформлены согласно требованиям ВУЗов и содержат рисунки, диаграммы, формулы и т.д.
PPT, PPTX и PDF-файлы представлены только в архивах.
Рекомендуем скачать работу.