в этой функции все и расчитывается.
function dzdt = jump(t, w)
global m c s sw g v0 teta ro Vmax anglew azimyt Val Vmax_z Hmax_v n N wan flag MT MV SM0 TM VM

x = w(1);
y= w(2);
z = w(3);
vx = w(4);
vy = w(5);
vz= w(6);
xdot = vx;
ydot = vy;
zdot=vz;
ro=1.225*(1-w(3)/44300)^4.256; %плотность воздуха в зависимости от высоты
S=s;
Sw=sw;
Vxy=sqrt(vx^2+vy^2); % полная скорость в плоскости Х0У
if MV
if abs(Vxy) > VM
S=SM0;   
Sw=SM0;
end;
end
if Val
[Vw,azimyt,Hnext]=getparamwind(w(3)); % функция которая возвращает скорость ветра и азимут в зависимости от высоты(табличные данные)
wx=Vw*sin(azimyt/180*pi); % состовляющая скорости ветра по оси Х
wy=Vw*cos(azimyt/180*pi); % состовляющая скорости ветра по оси У
if (wx & wy)
bw=(c*ro*Sw)/(2*m); % лобовое сопротивление ветра
vxwdot= -bw * wx * Vw; % ветровой напор по оси Х
vywdot= -bw * wy * Vw; % ветровой напор по оси У
else
vxwdot=0;
vywdot=0;
end
else
vxwdot=0;
vywdot=0;   
end
b = (c*ro*S)/(2*m); % лобовое сопротивлени фрагмента
v = sqrt(w(4)^2 + w(5)^2 + w(6)^2); % общая скорость фрагмента
vxdot = -b * vx * v; % сила сопротивления по оси Х
vydot= -b * vy * v; % сила сопротивления по оси Y
vzdot = -g - b * vz * v; % сила сопротивления по оси Z
vxdot=vxdot+vxwdot; %сложение силы сопротивления фрагмента воздуху с силой сопротивления ветра по оси Х
vydot=vydot+vywdot; % сложение силы сопротивления фрагмента воздуху с силой сопротивления ветра по оси У
dzdt =[ xdot; ydot; zdot;  vxdot; vydot; vzdot]; % данные отправляютя в решатель (внутренняя функция матлаб) ode45.

Отредактировано oper (2020-05-04 16:35:23)

Давайте попробуем.

Укажите массу, мидель и Сх.

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

Если у вас есть алгоритм изменения миделя, то я могу его применить в программе. Просто вам надо описать как это будет происходить.

Я пока не знаю как это реализовать. Но можно подумать.

Наибольшая  вовлечённость (к=1) будет для объектов, у которых объёмная плотность (m/V) будет равна весу вытесненного этим телом воздуха.
Проверим: Шар диаметром 1 м имеет объём V= 1/6* Пи* D³ = 0,52 м³
Вес такого объёма воздуха составит на земле 0,52 * 1,24 кг/м³ = 0,65 кг - выталкивающая сила.

Т.е. если такой шар будет весить 0,65 кг, то коэффициент его вовлечённости в движение воздуха составит 1  (нулевая плавучесть).
Положительной плавучести у нас быть не может, только отрицательная.
Значит, градация по отрицательной плавучести и может дать градацию по коэффициенту вовлечённости.

Для объектов одинакового объёма разный коэффициент вовлечённости будет полностью зависеть от веса объектов.
Чем больше вес, тем меньше коэффициент.
Проще всего принять 10-ти кратную шкалу (логарифмическую) уменьшения коэффициента.
Если выталкивающая сила в 10 раз меньше веса тела  коэффициент = 0,9
в 100 раз = 0,8, в 1000 раз = 0,7, и т.д.

Вес мы более-менее знаем. Осталось определить... объём тела! И тогда мы, рассчитав выталкивающую силу, узнаем коэффициент вовлечённости  и скорректируем скорость сноса ветром. умножив скорость ветра на этот коэффициент.

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

Да, по идее так можно сделать, но это сложнее, как мне кажется. Там возникает тогда куча синусов-косинусов с направление этого корыта относительно ветра. А ведь он стабилизаторов не имеет, поэтому его будет болтать по потоку. И мы не сможем установить положение корыта с минимальным лобовым сопротивлением, на которое оно и будет настраиваться в падении.

Поэтому и предлагаю идти по другому пути - уменьшение скорости ветра за счёт учёта коэффициента вовлечённости.  При этом будем рассчитыватся  вертикальная скорость падения по прямому миделю (синус от всей нижней площади) и скорость горизонтального сноса = скорости ветра * коэффициент вовлечённости + скорость горизонтальная от начальной горизонтальной.

Откеда могут быть истинными, если в отчете DSB разброс параметров по размеру и по весу одного и того же "скорпиона". Размеры 2х3 метра и на втором снимке 2,5х2,5.
А вес вообше с супер точностью 10-20 и 20-30 кг.

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

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

Ты столько колесил на велике ( говорил Сепаратист), мог бы взвесить и замерить и даже снять шаблон со скорпиона. Ты же местный?

В качестве предложения (Матлаб загрузить не смог без регистрации и без какой-то фигни, которая вместе с бесплатным загрузился - кое как вычистил):

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

Например:
На высоте 7800 ветер имел направление 210 град, а на высоте 5400 он изменил направление на 134 град, то
- разницу направлений (210-134=76) разделим на разницу высот (7800-5400 =2400) и получим изменение направления ветра на каждый метр высоты (0,031 град/м высоты). Проводим вычисления для всей высоты (10000м) по реперным точкам смены направления ветра и формируем массив "высота - направление ветра".

Так же  можно создать и массив "высота - скорость ветра)

При расчёте скоростей и направления сноса ветром мы будем знать на какой высоте какое будет точное направление и скорость ветра.

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

Если бы удар (взрыв) был справа. Но его нет. Значит только два варианта:
- таким было направление судна (хорошо подходит ракета с Красного Октября)
- так странно подействовала перегрузка торможения от взрыва.

Первый вариант легко проверить на программе, задав направление не 117, а 104 градуса и снова скинуть

Такого не может быть. Значит что-то не в порядке с ветром на высоте. Начальные условия по ветру покажешь?

Пример, как влияет ветер взависимости от массы объекта.
масса 10000 кг

масса 8000 кг

масса 6000 кг

масса 4000 кг

масса 2000 кг

масса 1000 кг

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

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

https://youtu.be/Vo15qE0dZvs?t=163

А вот так изобразил ее падение (с вырванным с ниши  передним шасси)  Сепаратист. Он видел ее падение. Она вращалась..

Также женщина, проживавшая в Торезе недалеко от шахты Прогресс, южнее.
https://youtu.be/ia55FmKp76c?t=25
Описывает после взрыва в воздухе падающий левее шахты предмет треугольной формы.Она связывает его появление явно с последствием взрыва, разметавшего в ее представлении самолет на куски. То есть, ниже облаков этот предмет треугольной формы при падении, как для женщины, визуально смещался не на шахту Прогресс, но от неё.

Отредактировано Спасский (2020-05-09 10:34:48)

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

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

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

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