Солнце, ветер и вода ver 0.2

в 14:56, , рубрики: аэродинамика, математическая модель, оцифровка яхты, транспорт, физика, физика движения

Сезон второй

  

Физика яхтинга или яхтинг для физиков

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

Что нужно чтобы оцифровать яхту?

  1. Буксировочная характеристика
  2. Диаграмма остойчивости
  3. Геометрия парусов — измеряю «в живую»
    высота топового угла грота $h_{гр топ}=11,2 м$;
    высота галсового угла грота $h_{гр глс}=1,985 м$;
    высота шкотового угла грота $h_{гр шкт}=2,179 м$;
    расстояние шкотового угла грота от мачты $l_{гр шкт}=2,96 м$;
    высота топового угла стакселя $h_{ст топ}=10 м$;
    высота галсового угла стакселя $h_{ст глс}=1,111 м$;
    высота шкотового угла стакселя $h_{ст шкт}=1,5 м$;
    площадь грота $S_{гр}=17,2 м^2$;
    площадь стакселя $S_{ст}=14 м^2$;
    площадь боковой парусности корпуса $S_{anf}=7,315 м^2$;
    площадь фронтальной парусности корпуса $S_{fr}=3 м^2$;
    боковая проекция площади подводной части корпуса $S_{дп} w=1,82 м^2$;
    площадь килей и скегов $S_{ks}=2,33 м^2$;
    минимально возможный угол установки стакселя $psi_{ст}=18^o$.

    Может это мне так не повезло, но ни один из производителей, во время моего поиска яхты, не согласился (не смог) предоставить эти данные для своей яхты. Я уверен, что у конструктора яхты вся эта информация есть, но получить ее почему то не получается. Буду добывать ее сам.
    Буксировочную характеристику я уже получил с помощью оцифровки яхты в программе freeship, в награду за упорство программа сразу выдает и диаграмму остойчивости и точки центра тяжести и центра бокового сопротивления и много другой полезной геометрической информации. Диаграмма остойчивости показывает какой крутящий момент надо приложить к яхте, чтобы накренить ее на определенный угол.
    Солнце, ветер и вода ver 0.2 - 15Солнце, ветер и вода ver 0.2 - 16

    Все необходимые аргументы собраны, приступим к расчетам.

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

    $R_{fw}=0,61*S_{fr}*v^2, Н,$

    где v — скорость судна.
    В то же время на яхту действует сила реального ветра Frw, дующего под некоторым углом rw. Оба этих ветра создают силу равную векторной сумме частичных сил. Так называемую силу вымпельного ветра — ветра, дующего на борту судна.
    Основная проблема заключается в теории относительности. Наблюдатель (и все измерительные приборы) находится на борту яхты и для него нет возможности измерить силу и направление реального ветра, зато он может измерить направление — $alpha_{vw}$ и скорость вымпельного ветра — $v_{vw}$, а также направление — $alpha$, измеренное бортовым компасом и скорость судна — $v$, измеренную бортовым лагом (спидометром).
    Нахожу параметры реального ветра на основе измерений бортовых приборов.
    $v_{rw}=sqrt{v^2+v_{vw}^2-2*v*v_{vw}*cos(alpha_{vw})}$ — скорость реального ветра;
    $alpha_{rw}=pi-arccos(frac{v_{rw}^2+v^2-v_{vw}^2}{2*v*v_{rw}})$ — угол реального ветра отложенный от направления «вперед».
    Солнце, ветер и вода ver 0.2 - 24
    Для маленьких покатушек все это не имеет смысла. Связь реального и вымпельного ветра необходима для планирования достаточно длительных путешествий (когда цель не находится в прямой видимости). Ведь планирование выполняется на на карте планеты с указанием реальных ветров (по прогнозу погоды), а движение яхты происходит от вымпельного ветра.

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

    1. режим аэродинамического крыла — используется при встречных ветрах,
    2. режим аэродинамического тормоза — используется при попутных ветрах.

    Профиль паруса шьется таким образом, чтобы, будучи наполненным встречным ветром, принимать форму крыла. Т.е., есть крыло, движущееся с некоторой скоростью относительно воздуха — следовательно в нем возникает подъемная сила $F_{air}$. По аналогии с крылом самолета попробую упрощенно оценить ее величину и направление, относительно яхты.
    Чтобы парус принял рабочую форму его необходимо немного повернуть относительно направления ветра. Угол поворота называется «угол атаки» $alpha_{a}$. Для простоты расчетов крыло сравнивают с плоской пластинкой, а различия представлены в виде таблицы аэродинамических коэффициентов, где Су — коэффициент отличия подъемной силы, а Сх — коэффициент отличия силы сопротивления. Подъемная сила направлена перпендикулярно пластинке, а сила сопротивления — параллельно. Геометрия пластинки задается коэффициентом аэродинамического удлинения $A_{y}=h^2/S$, где h — высота паруса; S  — площадь паруса.

    а Cy;Ay=6 Cx Cy;Ay=3 Cx
    0 0 0,09 0 0,09
    5 0,83 0,1 0,65 0,1
    10 1,19 0,15 0,9 0,15
    15 1,22 0,3 1,1 0,25
    20 1,14 0,4 1,18 0,4
    25 1,06 0,485 1,2 0,5
    30 0,98 0,57 1 0,55
    40 0,9 0,73 0,82 0,65
    45 0,75 0,88 0,7 0,85
    50 0,6 1 0,6 1
    70 0,2 1,15 0,2 1,15
    90 0 1,2 0 1,2

    Крыло, в отличии от пластинки, трехмерная конструкция, поэтому отдельная таблица показывает как влияет размер выпуклости «пуза» на коэффициенты отличия от пластинки.

    пузо Су3 Сх3
    0,0667 -0,2 0,1
    0,1 0 0
    0,2 0 0,1
    0,25 0,1 0,2

    Проекция силы паруса $F_{air}$ на продольное направление яхты — полезная сила:
    $F_{s}=0,61*v_{vw}^2*S*((Cy+Cy3)*sin(alpha_{vw})-(Cx+Cx3)*cos(alpha_{vw}))$.
    Проекция силы паруса $F_{air}$ на поперечное направление яхты — сила дрейфа:
    $F_{d}=0,61*v_{vw}^2*S*((Cy+Cy3)*cos(alpha_{vw})-(Cx+Cx3)*sin(alpha_{vw}))$.
    Солнце, ветер и вода ver 0.2 - 32
    Сила ветрового сопротивления корпуса будет замедлять движение:
    $R_{fw}=0,61*S_{fr}*v_{vw}^2*sin(alpha_{vw})$
    Получив проекции сил, могу найти вектора продольной и поперечной скорости яхты. Продольную составляющую скорости $v_{дп}$ нахожу по графику буксировочной характеристики.
    С поперечной составляющей все сложнее.
    Для начала надо найти высоту центра парусности паруса. Парус представляет собой треугольник у которого одна сторона — выпуклая дуга. Представлю его как два паруса: треугольный и серпообразный кусок. Площадь треугольника вычисляется как половина произведения ширины на высоту, а площадь «серпа» это разность между площадью паруса и площадью треугольника.

    $ S=(h_{гр топ}-h_{гр глс})*l_{гр шкт}/2,$

    где $S_{серп}=S_{гр}-S$
    Высота центра парусности паруса:

    $h_{гр цп}=frac{h_{гр Delta цп}*SDelta+((h_{гр топ}-h_{гр шкт})*S_{серп}/1,8)}{S_{гр}},$

    где $h_{гр Delta цп}=h_{гр шкт}+(h_{гр топ}+h_{гр глс}-2*h_{гр шкт})/3 ,$ высота центра парусности треугольника паруса без серпа.
    Обычно яхта имеет два паруса: основной — грот и передний — стаксель. Стаксель обычно выполняется без серпа и тогда его высота центра парусности:

    $h_{ст Delta цп}=h_{ст шкт}+(h_{ст топ}+h_{ст глс}-2*h_{ст шкт})/3$

    Суммарный центр парусности:

    $h_{Sigma цп}=frac{h_{гр цп}*S_{гр}+h_{ст Delta цп}*S_{ст}}{S_{гр}+S_{ст}}$

    Теперь могу определить кренящий момент, действующий на яхту:

    $M_{кр}=F_{d}*h_{Sigma цп}.$

    По диаграмме остойчивости нахожу угол крена яхты $phi$.
    Выражаю силу сопротивления дрейфу яхты:

    $R_{d}=9,8*102*(S_{ks}*atan(frac{v_{d}}{v_{дп}})*v_{дп}^2*(cos(phi-gamma)-cos(phi+gamma))*0,96+(S_{ks}+S_{дп w})*1,15*v_{d}^2*cos(phi-gamma)),$

    где — угол наклона килей двух-килевой яхты относительно вертикали.
    Сила ветрового давления на корпус придает дополнительный дрейф.

    $F_{anf}=0,61*v_{d}^2*S_{anf}*cos(alpha_{vw}).$

    Скорость $v_{d}$ находится подстановкой по методу половинного деления до тех пор, пока не наступит баланс сил $F_{d}+F_{anf}=R_{d}$.
    Теперь очевидно, что яхта движется " косо", т.е. имеет поперечную и продольную скорость перемещения. Это приносит еще один сюрприз. Направление — $alpha$, измеренное бортовым компасом и скорость судна — $v=v_{дп}$, измеренное бортовым лагом (спидометром) показывают только продольную составляющую.
    Надо привязать измерения к реальности. Реальный угол перемещения судна $alpha_{real}$ называется «путевой угол»:

    $alpha_{real}=atan(frac{v_{d}}{v_{дп}})+alpha$

    и может отличаться от измеренного на 10-20 градусов.
    Реальная скорость яхты может быть определена по теореме Пифагора:

    $v_{real}=sqrt{v_{дп}^2+v_{d}^2}.$

    В завершении приведу пример расчета для моей яхты самой интересной ситуации — движение против ветра. Понятно, что строго против ветра под парусами не пойдешь, но под некоторым острым углом можно. 
    Начать надо с переднего паруса — стаксель закреплен сверху на мачте, передним углом к носу яхты, а задним с помощью шкота (мягкий трос) через ролик к лебедке. Его натяжение регулирует угол установки паруса относительно яхты. Минимальный угол установки $psi_{ст мин}=18^o$ получается при максимальном натяжении шкота. Аэродинамическое удлинение стакселя $A_{y ст}=h_{ст}^2/S_{ст}=5,8$. Методом линейной интерполяции нахожу оптимальный угол атаки $alpha_{a}=10^o$ из таблицы аэродинамических коэффициентов. Критерий оптимальности — максимизация $F_{s ст}$  при этом Су=1,17; Сх=0,15. Установленный таким образом стаксель будет эффективно работать при направлении вымпельного ветра $alpha_{vw}=psi_{ст}+alpha_{a}=18+10=28^o$. Хорошая прогулочная яхтенная погода, когда скорость ветра в районе 5-7 м/с. Для красоты цифр возьму $v_{vw}=6,55 м/с$
    Проекция силы стакселя на продольное направление яхты — полезная сила:

    $F_{s ст}=0,61*6,55^2*14*((1,17+0)*sin(28*pi/180)-(0,15+0)*cos(28*pi/180))=153,5 H$

    Проекция силы стакселя на поперечное направление яхты — сила дрейфа:

    $F_{d ст}=0,61*6,55^2*14*((1,17+0)*cos(28*pi/180)-(0,15+0)*sin(28*pi/180))=405,8 H$

    Ветровое сопротивление корпуса препятствует движению вперед.

    $R_{fr w}=0,61*6,55^2*3*cos(28*pi/180)=69 H$

    По буксировочной характеристике определяем скорость $v_{дп}=2,75$ узла=1,4 м/с

    $h_{ст Delta цп}=1,5+(10+1,5-2*1,1)/3=4,6 м$

    $M_{кр ст}=405,8*4,6=1867 Н/м$, по диаграмме остойчивости находим угол крена $phi_{ст}=3,9^o$. Ну это мелочи, потому добавим еще один парус — грот!
    Грот работает как закрылок у крыла самолета и вращается вокруг мачты на угол $70^o$. Коэффициент аэродинамического удлинения грота $A_{y гр}=h_{гр}^2/S_{гр}=4,9$;
    оптимальный угол атаки $alpha_{a}=10^o$; Су=1,09; Сх=0,15.
    Проекция силы грота на продольное направление яхты — полезная сила:

    $F_{s гр}=0,61*6,55^2*17,2*((1,09+0)*sin(28*pi/180)-(0,15+0)*cos(28*pi/180))=170 H$

    Проекция силы грота на поперечное направление яхты — сила дрейфа:

    $F_{d гр}=0,61*6,55^2*17,2*((1,09+0)*cos(28*pi/180)-(0,15+0)*sin(28*pi/180)))=464 H$

    $h_{Sigma цп}=4,16 м$, $M_{кр}=3127 H$, а угол крена $phi_{Sigma}=8,5^o$.
    Суммарная сила движения вперед:
    $F_{s}=F_{s ст}+F_{s гр}-R_{fr w}=153,5+170-69=254 H$, скорость движения вперед $v_{дп}=3,6$ узла или 1,86 м/с.
    Если предположить, что реальный ветер дует строго северный, то компас на борту будет показывать угол продольной составляющей скорости яхты $alpha_{vw}=38,1^o$.
    Скорость дрейфа составит $v_{d}=0,418 м/с$.
    Теперь надо привести эти результаты к реальности. 
    Путевой угол реального движения составит: 
    $alpha_{real}=atan(frac{v_{d}}{v_{дп}})+alpha=atan(frac{0,418}{1,86})+38,1=51^o$, это и есть угол реального ветра $alpha_{rw}$ к вектору движения яхты. 
    А реальная скорость перемещения в пространстве: $v_{real}=sqrt{v_{дп}^2+v_{d}^2}$=1,92+0,4182=1,9 м/с. 
    Скорость реального ветра:
    $v_{rw}=sqrt{v^2+v_{vw}^2-2*v*v_{vw}*cos(alpha_{vw})}=5 м/с$.
    При таких условиях скорость продвижения строго против ветра составит 1,2 м/с или 2,35 узла, а двигаться придется по зигзагообразной траектории — галсами.Солнце, ветер и вода ver 0.2 - 81
    Проведя расчет для других возможных углов вымпельного ветра можно получить круговую лавировочную диаграмму зависимости реальной скорости яхты от реального ветра. Солнце, ветер и вода ver 0.2 - 82С помощью ее уже можно планировать маршруты на карте по прогнозу ветров. Дополнительно стало понятно, что минимально возможный угол $alpha_{real min}=43,5^o$, а максимальная скорость движения против ветра достигается с курсовым углом $50^o<alpha_{real}<53^o$ для ветра 5 м/с.

Автор: Василий Карпов

Источник

* - обязательные к заполнению поля


https://ajax.googleapis.com/ajax/libs/jquery/3.4.1/jquery.min.js