Как сбит MH17: гражданское расследование

Информация о пользователе

Привет, Гость! Войдите или зарегистрируйтесь.


Вы здесь » Как сбит MH17: гражданское расследование » Математическое и компьютерное моделирование » Математическое моделирование траектории падения центроплана


Математическое моделирование траектории падения центроплана

Сообщений 1 страница 30 из 2006

1

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

Если сбросить балванку весом 250 тонн с высоты 10000 метров  без крылошек под углом к горизонту 70 градусов. то максимальная скорость у земли будет 460 м/сек. http://forumupload.ru/uploads/0016/23/c6/51/804874.png
http://forumupload.ru/uploads/0016/23/c6/51/804874.png

0

2

...............Советую почитать отчет Алмаз-Антей
Исследование зависимости баллистических характеристик ракет от параметров атмосферы.
http://www.almaz-antey.ru/upload/iblock … be9db0.pdf

http://d.radikal.ru/d38/2004/df/115884f029e4.jpg

http://d.radikal.ru/d00/2004/bd/b6437b4a23bc.jpg
900 км\ч = 250 м\с
- скорость звука на высоте 0 м......340 м\с

Материалы всероссийской научно-практической конференции
https://iupr.ru/domains_data/files/Cborniki_konf/7-14 maya 2019.pdf

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

Отредактировано Ралив (2020-04-27 17:09:29)

+1

3

Ралив написал(а):

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

Сбросим с высоты 6400 метров при скорости бомбардировщика 900 км/час ФАБ-250-М54. Бомба приземлится со скоростью 270м/сек. http://forumupload.ru/uploads/0016/23/c6/51/915120.png
http://forumupload.ru/uploads/0016/23/c6/51/915120.png

0

4

oper написал(а):

Сбросим с высоты 6400 метров при скорости бомбардировщика 900 км/час ФАБ-250-М54. Бомба приземлится со скоростью 270м/сек. http://forumupload.ru/uploads/0016/23/c6/51/915120.png

Опер! Разгон боинга по траектории басара на своей адской машинке когда проверишь? :)

https://www.politforums.net/redir/forei … 1586628869

Отредактировано Спасский (2020-04-27 18:28:06)

0

5

Спасский написал(а):

Опер! Разгон боинга по траектории басара на своей адской машинке когда проверишь? :)

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

+1

6

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

0

7

Ралив написал(а):

Сх = 0,49 - по всей видимости шар.

Сфера: Cx=0,47
В сущности мидель надо изменить так как бомба падает не вертикально, а под углом и он будет больше чем поперечное сечение цилиндра бомбы.
Достоверно известно только высота,скорость, масса и размеры бомбы недастает только Сх бомбы.

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

0

8

Спасский написал(а):

взял  какую-то  условную болванку

А если крылышки приделать к болванки, то она еще медленнее будет падать. У вас видимо "МН17" падает в безвоздушном пространстве и скорости аховские.

0

9

Спасибо Простой, что открыли важную тему.

Спасибо Раливу за эту ссылку https://www.webcitation.org/65Tg5GmBM?u … listic.pdf.

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

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

+1

10

Как обещал публикую код программы и любой загрузив Матлаб может использовать этот код.

весь код записан в файл с расширением .m.

function varargout = ballistic(varargin)
% BALLISTIC Application M-file for ballistic.fig
%    FIG = BALLISTIC launch ballistic GUI.
%    BALLISTIC('callback_name', ...) invoke the named callback.

% Last Modified by GUIDE v2.0 27-Apr-2020 12:52:21

if nargin == 0  % LAUNCH GUI

fig = openfig(mfilename,'reuse');

% Use system color scheme for figure:
set(fig,'Color',get(0,'defaultUicontrolBackgroundColor'));

% Generate a structure of handles to pass to callbacks, and store it.
handles = guihandles(fig);
guidata(fig, handles);

if nargout > 0
    varargout{1} = fig;
end

elseif ischar(varargin{1}) % INVOKE NAMED SUBFUNCTION OR CALLBACK

try
    if (nargout)
    [varargout{1:nargout}] = feval(varargin{:}); % FEVAL switchyard
    else
    feval(varargin{:}); % FEVAL switchyard
    end
catch
    disp(lasterr);
end

end
%---------------------------------------------------------------
function varargout = pushbutton1_Callback(h, eventdata, handles, varargin)
cla;
global m c s g v0 teta ro H Vmax L anglew azimyt Vw Val Vmax_z Hmax_v n N T wan
m=[]; c=[]; s=[]; g=[]; v0=[]; teta=[]; ro=[]; H=[]; Vmax=[]; L=[]; anglew=[]; azimyt=[]; Vw=[];
Lx=0;
Ly=0;
Vmax_z=0;
Hmax_v=0;
n=0;
N=0;
wan=1;
H=str2double(get(handles.edit1, 'String'));     
m = str2double(get(handles.edit3, 'String'));
c = str2double(get(handles.edit4, 'String'));
s = str2double(get(handles.edit6, 'String'));
v0 = str2double(get(handles.edit2, 'String'));
if get(handles.radiobutton1,'Value')
Val=1;
else
    Val=0;
end
angleG= str2double(get(handles.edit5, 'String'));
g = 9.81;
teta =-angleG/180*pi;
anglew=str2double(get(handles.edit7, 'String'));
kyrs=(anglew+180)/180*pi; % угол курса
[t w te ye ie ]=ODE(teta,kyrs,v0,H);

plot3(w(:,1), w(:,2), w(:,3), 'LineWidth', 2);
set(gca, 'YDir', 'reverse');
grid on;
set(handles.edit8,'String',num2str(floor(te)));
Lx=floor(ye(1));
set(handles.edit9,'String',num2str(Lx));
Ly=floor(ye(2));
set(handles.edit10,'String',num2str(Ly));
L=floor(sqrt(Lx^2+Ly^2));
set(handles.edit11,'String',num2str(L));
Vz=abs((ye(6)));
set(handles.edit12,'String',num2str(Vz));
set(handles.edit13,'String',num2str(abs(Vmax_z)));
%--------------------------------------------------------------
function varargout = edit1_Callback(h, eventdata, handles, varargin)

% --------------------------------------------------------------------
function varargout = edit2_Callback(h, eventdata, handles, varargin)

% --------------------------------------------------------------------
function varargout = edit3_Callback(h, eventdata, handles, varargin)

% --------------------------------------------------------------------
function varargout = edit4_Callback(h, eventdata, handles, varargin)

% --------------------------------------------------------------------
function varargout = edit5_Callback(h, eventdata, handles, varargin)

% --------------------------------------------------------------------
function varargout = edit6_Callback(h, eventdata, handles, varargin)

%-----------------------------------------------------------------
function dzdt = jump(t, w)
global m c s g v0 teta ro Vmax anglew azimyt Vw v Val
vx = w(4);
vy = w(5);
vz= w(6);
xdot = vx;
ydot = vy;
zdot=vz;
if Val
[Vw,azimyt]=getparamwind(w(3));
wx=Vw*sin((azimyt+180)/180*pi);
wy=Vw*cos((azimyt+180)/180*pi);
if (wx & wy)
vx=vx+wx;
vy=vy+wy;
v = sqrt((w(4)+wx)^2 + (w(5)+wy)^2 + w(6)^2);
end
else
v = sqrt(w(4)^2 + w(5)^2 + w(6)^2);
end
ro=1.225*(1-w(3)/44300)^4.256;
b = (c*ro*s)/(2*m);

vxdot = -b * vx * v;
vydot= -b * vy * v;
vzdot = -g - b * vz * v;
dzdt =[xdot; ydot; zdot; vxdot; vydot; vzdot];
%----------------------------------------------------------
function [value,isterminal,direction] = ex8events(t,w)
value = [w(3)];
isterminal = [0];
direction = [-1];
%-----------------------------------------------------------
function [wx,wy]=getspeed(w,kyrs,anglew)
wx=w*cos(anglew/180*pi);
wy=w*sin(anglew/180*pi);
if (kyrs>0&kyrs<90)
    if(anglew>90&anglew<180)
        wx=-wx;
    end
if(anglew>180&anglew<270)
        wx=-wx;
        wy=-wy;
    end
if(anglew>270&anglew<360)
        wy=-wy;
    end   
end
if (kyrs>90&kyrs<180)
    if(anglew>0&anglew<90)
        wx=-wx;
    end
if(anglew>180&anglew<270)
        wy=-wy;
    end
if(anglew>270&anglew<360)
        wx=-wx;
        wy=-wy;
    end   
end
%----------------------------------------

Отредактировано oper (2020-12-08 20:34:31)

0

11

global m c s g v0 teta ro H Vmax L anglew azimyt Vw Val Vmax_z Hmax_v n N T wan
m=[]; c=[]; s=[]; g=[]; v0=[]; teta=[]; ro=[]; H=[]; Vmax=[]; L=[]; anglew=[]; azimyt=[]; Vw=[];

Хорошо бы дать расшифровку этих параметров.
Да и комментарии к блокам расчёта тоже.

0

12

Продолжение.

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

Данные взял у uschen с ексель файла программы.

function [w,anglew]=getparamwind(Hz)
global Vw azimyt
if(Hz>=0 & Hz<185)
    Vw=5; azimyt=260;
end
if(Hz>=185 & Hz<380)
    Vw=9; azimyt=285;
end
if(Hz>=380 & Hz<668)
    Vw=10; azimyt=277;
end
if(Hz>=668 & Hz<771)
    Vw=11; azimyt=265;
end
if(Hz>=771 & Hz<1165)
    Vw=11; azimyt=265;
end
if(Hz>=1165 & Hz<1382)
    Vw=11; azimyt=268;
end
if(Hz>=1382 & Hz<1510)
    Vw=11; azimyt=269;
end
if(Hz>=1510 & Hz<1934)
    Vw=11; azimyt=270;
end
if(Hz>=1934 & Hz<2087)
    Vw=12; azimyt=285;
end
if(Hz>=2087 & Hz<2206)
    Vw=13; azimyt=291;
end
if(Hz>=2206 & Hz<2935)
    Vw=11; azimyt=295;
end
if(Hz>=2935 & Hz<2957)
    Vw=6; azimyt=318;
end
if(Hz>=2957 & Hz<3001)
    Vw=6; azimyt=319;
end
if(Hz>=3001 & Hz<3142)
    Vw=5; azimyt=320;
end
if(Hz>=3142 & Hz<3303)
    Vw=4; azimyt=325;
end
if(Hz>=3303 & Hz<3423)
    Vw=3; azimyt=350;
end
if(Hz>=3423 & Hz<4039)
    Vw=3; azimyt=348;
end
if(Hz>=4039 & Hz<4427)
    Vw=3; azimyt=340;
end
if(Hz>=4427 & Hz<4682)
    Vw=3; azimyt=335;
end
if(Hz>=4682 & Hz<4776)
    Vw=3; azimyt=347;
end
if(Hz>=4776 & Hz<5271)
    Vw=3; azimyt=352;
end
if(Hz>=5271 & Hz<5298)
    Vw=2; azimyt=16;
end
if(Hz>=5298 & Hz<5445)
    Vw=2; azimyt=18;
end
if(Hz>=5445 & Hz<5606)
    Vw=2; azimyt=25;
end
if(Hz>=5606 & Hz<5820)
    Vw=3; azimyt=16;
end
if(Hz>=5820 & Hz<6020)
    Vw=4; azimyt=5;
end
if(Hz>=6020 & Hz<6622)
    Vw=5; azimyt=5;
end
if(Hz>=6622 & Hz<6855)
    Vw=8; azimyt=5;
end
if(Hz>=6855 & Hz<6939)
    Vw=9; azimyt=18;
end
if(Hz>=6939 & Hz<7273)
    Vw=9; azimyt=23;
end
if(Hz>=7273 & Hz<7490)
    Vw=10; azimyt=43;
end
if(Hz>=7490 & Hz<7531)
    Vw=10; azimyt=55;
end
if(Hz>=7857 & Hz<8286)
    Vw=9; azimyt=60;
end
if(Hz>=8286 & Hz<9102)
    Vw=9; azimyt=35;
end
if(Hz>=9102 & Hz<9530)
    Vw=11; azimyt=40;
end
if(Hz>=9530 & Hz<10060)
    Vw=15; azimyt=40;
end
if(Hz>=10060)
    Vw=20; azimyt=45;
end
w=Vw;
anglew=azimyt;
%---------------------------------------------------------

function varargout = edit7_Callback(h, eventdata, handles, varargin)

% --------------------------------------------------------------------
function varargout = pushbutton2_Callback(h, eventdata, handles, varargin)
cla;
global m c s sw g v0 teta ro H Vmax L anglew azimyt Vw
m=[]; c=[]; s=[]; g=[]; v0=[]; teta=[]; ro=[]; H=[]; Vmax=[]; L=[]; anglew=[]; azimyt=[]; Vw=[];
Lx=0;
Ly=0;
H=str2double(get(handles.edit1, 'String'));
m = str2double(get(handles.edit3, 'String'));
c = str2double(get(handles.edit4, 'String'));
s = str2double(get(handles.edit6, 'String'));
sw = str2double(get(handles.edit14, 'String'));
v0 = str2double(get(handles.edit2, 'String'));
angleG= str2double(get(handles.edit5, 'String'));
g = 9.81;
teta =-angleG/180*pi;
anglew=str2double(get(handles.edit7, 'String'));
kyrs=anglew/180*pi; % угол курса 
[t w te ye ie ]=ODE(teta,kyrs,v0,H);
plot(w(:,1), w(:,2), 'LineWidth', 3);
grid on;
te
Lx=floor(ye(1))
Ly=floor(ye(2))
L=floor(sqrt(Lx^2+Ly^2))
Vz=abs(floor(ye(6)))
% --------------------------------------------------------------------
function varargout = edit8_Callback(h, eventdata, handles, varargin)

% --------------------------------------------------------------------
function varargout = edit9_Callback(h, eventdata, handles, varargin)

% --------------------------------------------------------------------
function varargout = edit10_Callback(h, eventdata, handles, varargin)

% --------------------------------------------------------------------
function varargout = edit11_Callback(h, eventdata, handles, varargin)

% --------------------------------------------------------------------
function varargout = edit12_Callback(h, eventdata, handles, varargin)

% --------------------------------------------------------------------
function varargout = radiobutton1_Callback(h, eventdata, handles, varargin)

% --------------------------------------------------------------------
function varargout = edit13_Callback(h, eventdata, handles, varargin)

%----------------------------------------------------------------------
function [T, W, TE, YE, IE]=ODE(teta, kyrs, v0,H)
tmax = 10000;
dt = [0 tmax];
w0=[];
w0 = [0 0 H   v0*cos(teta)*sin(kyrs) v0*cos(teta)*cos(kyrs) v0*sin(teta)];
opt = odeset('Events',@ex8events,'RelTol', 1e-8, 'AbsTol', 1e-8);
[t w te ye ie ] = ode45(@jump, dt, w0, opt);%=================================решатель диффиренциальных уравнений

if te>0
tmax=floor(te);
end
dt = [0 tmax];
t=[];
w=[];
Hmax_v=[0 0];
N=0;
wan=0;
w0 = [0 0 H   v0*cos(teta)*sin(kyrs) v0*cos(teta)*cos(kyrs) v0*sin(teta)];
opt = odeset('RelTol', 1e-8, 'AbsTol', 1e-8);
[t w ] = ode45(@jump, dt, w0, opt);%=================================решатель диффиренциальных уравнений
T=t;
W=w;
TE=te;
YE=ye
IE=ie;
%----------------------------------------------------
function varargout = edit14_Callback(h, eventdata, handles, varargin)

% --------------------------------------------------------------------
function varargout = pushbutton3_Callback(h, eventdata, handles, varargin)
cla;
global m c s sw g v0 teta ro H Vmax L anglew azimyt Vw
m=[]; c=[]; s=[]; g=[]; v0=[]; teta=[]; ro=[]; H=[]; Vmax=[]; L=[]; anglew=[]; azimyt=[]; Vw=[];
Lx=0;
Ly=0;
H=str2double(get(handles.edit1, 'String'));
m = str2double(get(handles.edit3, 'String'));
c = str2double(get(handles.edit4, 'String'));
s = str2double(get(handles.edit6, 'String'));
sw = str2double(get(handles.edit14, 'String'));
v0 = str2double(get(handles.edit2, 'String'));
angleG= str2double(get(handles.edit5, 'String'));
g = 9.81;
teta =-angleG/180*pi;
anglew=str2double(get(handles.edit7, 'String'));
kyrs=anglew/180*pi; % угол курса 
[t w te ye ie ]=ODE(teta,kyrs,v0,H);
plot(w(:,2), w(:,3), 'LineWidth', 3);
set(gca, 'XDir', 'reverse');

grid on;
te
Lx=floor(ye(1))
Ly=floor(ye(2))
L=floor(sqrt(Lx^2+Ly^2))
Vz=abs(floor(ye(6)))
%-----------------------------------------------------------
function varargout = pushbutton4_Callback(h, eventdata, handles, varargin)
cla;
global m c s sw g v0 teta ro H Vmax L anglew azimyt Vw
m=[]; c=[]; s=[]; g=[]; v0=[]; teta=[]; ro=[]; H=[]; Vmax=[]; L=[]; anglew=[]; azimyt=[]; Vw=[];
Lx=0;
Ly=0;
H=str2double(get(handles.edit1, 'String'));
m = str2double(get(handles.edit3, 'String'));
c = str2double(get(handles.edit4, 'String'));
s = str2double(get(handles.edit6, 'String'));
sw = str2double(get(handles.edit14, 'String'));
v0 = str2double(get(handles.edit2, 'String'));
angleG= str2double(get(handles.edit5, 'String'));
g = 9.81;
teta =-angleG/180*pi;
anglew=str2double(get(handles.edit7, 'String'));
kyrs=anglew/180*pi; % угол курса 
[t w te ye ie ]=ODE(teta,kyrs,v0,H);
plot(w(:,1), w(:,3), 'LineWidth', 3);
grid on;
te
Lx=floor(ye(1))
Ly=floor(ye(2))
L=floor(sqrt(Lx^2+Ly^2))
Vz=abs(floor(ye(6)))

Отредактировано oper (2020-05-02 18:13:58)

0

13

Простой написал(а):

global m c s g v0 teta ro H Vmax L anglew azimyt Vw Val Vmax_z Hmax_v n N T wan
    m=[]; c=[]; s=[]; g=[]; v0=[]; teta=[]; ro=[]; H=[]; Vmax=[]; L=[]; anglew=[]; azimyt=[]; Vw=[];

Хорошо бы дать расшифровку этих параметров.
Да и комментарии к блокам расчёта тоже.

m - масса
с - коэффициент лобового сопротивления (Сх)
s - мидель
g - ускорение свободного падения,  константа, 9,81
v0 - начальная скорость
teta - угол к горизонту в радианах
ro - плотность воздуха
H - начальная высота
L - расстояние по прямой от проекции на горизонтальную плоскость начала траектории до точки падения на горизонтальную плоскость
Vmax - зарезервировано (гдето хотел применить так и осталась)
anglew - направление - курс в градусах
azimyt - метерологический азимут ветра
Vw - скорость ветра
Val - зарезервировано (гдето хотел применить так и осталась)
Vmax_z - максимальная скорость
Hmax_v - зарезервировано (гдето хотел применить так и осталась)
n - колличество циклов в момент достижения максимальной скорости
N - общее колличество циклов 
T - резерв
wan - флаг принимает два значения 1 и 0

m=[]; c=[]; s=[]; g=[]; v0=[]; teta=[]; ro=[]; H=[]; Vmax=[]; L=[]; anglew=[]; azimyt=[]; Vw=[]; все данные стираются.

+1

14

Замечательно!
А текущие значения в циклах

x = w(1);
y= w(2);
z = w(3);
vx = w(4);
vy = w(5);
vz= w(6);
xdot = vx;
ydot = vy;
zdot=vz;

0

15

Простой написал(а):

Замечательно!
А текущие значения в циклах

    x = w(1);
    y= w(2);
    z = w(3);
    vx = w(4);
    vy = w(5);
    vz= w(6);
    xdot = vx;
    ydot = vy;
    zdot=vz;

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

W это и есть многомерная матрица с шестью столбцами. Функция jump имеет два аргумента t и w, первый время, второй матрица w в которые решателем ode45 (внутренняя функция Матлаб) при вызове записываются предыдушие значения, которые и подставляются в дифференциальные уравнения и так цикл за циклом.
x, y, z координаты в момент времени t, vx, vy, vz -  состовляюшие скорости.

И эта матрица w возвращается решателем [t w ] = ode45(@jump, dt, w0, opt); и она используется в конечном итоге функциями plot и plot3 для отображения графиков.

Отредактировано oper (2020-04-28 00:13:26)

0

16

Это к вопросу учёта площади "утюга"
Мидель - видимая площадь объекта перпендикулярная его движению.
Синий вектор - направление движения объекта
http://forumupload.ru/uploads/0016/23/c6/2/t493043.jpg

0

17

http://mh17.webtalk.ru/viewtopic.php?id … =18#p36405

Анна написал(а):
1. Правда, мне всё равно непонятны Ваши картинки.
2. Вот по Вашей картинке указан вес 0,40. Это сколько - 400 грамм или 40 кг? И то и другое на мой взгляд не катит. Обломок с двумя шпангоутами и дырками примерно размером метр на метр, весит около 5 кг. Всё сугубо субъективно.
3. Второй обломок - скальп, который тоже с дырками, размером примерно 2х1,5 метра. весом около 20 кг. Лежат на одном поле (раз уж никто не захотел сделать точную локацию, то пусть будут на одном поле, т.к. точные координаты есть только для скальпа).

Так вот... Можно ли каким-либо образом  подсчитать их путь в обратную сторону, то есть наверх? Чтобы определить высоту и точку, с которой они начали свободное падение? Либо же скажите мне просто по-русски - с какой высоты они отвалились, т.к. по картинкам я нифига не понимаю.(((( Просто словами...

Хорошо бы послушать Сепаратиста. Действительно ли  "второй обломок - скальп, который тоже с дырками, размером примерно 2х1,5 метра. весом около 20 кг"?
Сепаратист измерял и взвешиваал "скальп сепаратиста"?

Это скальп первым оторвался и ушел в свободное падение и можно было бы проверить действительно ли он падал с высоты 10060 метров. Есть ли у кого данные по этому обломку?

0

18

oper написал(а):

Хорошо бы послушать Сепаратиста

Не взвешивал :dontknow: .
Лежал на мусорнике. Значит кто то туда его принёс. Значит упал кому то во двор. Скорее всего по крайней улице от поля в Рассыпном.

Отредактировано Сепаратист (2020-05-02 02:24:00)

0

19

ГРОЗ написал(а):

А можно координаты этого "мусорника" ?

В фильме разве  не видно? От пятиэтажек в Рассыпном дорога на Петропавловку. На выезде когда то был ж.д. переезд. Теперь там стоит крест... и была небольшая незаконная .свалка. Мэр наверное кипишь и поднял чтобы его за этот мусорник не вздрючили.

0

20

ГРОЗ написал(а):

Первым оторвался фрагмент под названием "Скорпион".

Анна говорит про фрагмент

2. Вот по Вашей картинке указан вес 0,40. Это сколько - 400 грамм или 40 кг? И то и другое на мой взгляд не катит. Обломок с двумя шпангоутами и дырками примерно размером метр на метр, весит около 5 кг. Всё сугубо субъективно.

По ее словам фрагмент метр на метр и весом около 5 кг. Это "скорпион"?

0

21

Добавил в программу возможность отключать влияние ветра, добавленны кнопки отображения 2D графика в плоскости ZОY и ZOX, вертикальные плоскость юг и запад и добавлена возможность вводить мидель фрагмента при учете влияния ветра.

Поэтому выше отредактировал листинг программы.

0

22

http://mh17.webtalk.ru/viewtopic.php?id … =15#p86641

Сепаратист написал

Бутблэк,правый это один из тех что Яна отдала в фильме в петропавловский сельсовет.По-видимому в облаке у Вас фотки с того склада. Левый - это скальп-скорпион,который исчез ,не дождавшись сбора обломков.

Как так исчез? Сепаратист вы успели его взвесить и снять размеры?

0

23

Вот и данные по скальпам, с большим разбросом по весу. Данные их окончательного очета, раздел - MH17 About the investigation.
http://forumupload.ru/uploads/0016/23/c6/51/576460.png

Будем пробывать моделировать.

0

24

oper написал(а):

http://mh17.webtalk.ru/viewtopic.php?id … =15#p86641
Как так исчез? Сепаратист вы успели его взвесить и снять размеры?

Вот так - исчез! После падения самолёта до освобождения Дебальцево я в Петропавловке ни разу не был. Там где скальп- скорпион вообще территория бд. Ходить с фотоаппаратом среди позиций военных нужно быть либо чокнутым либо с лицензией на разрешение. Знал командира и его телефон .На просьбу поснимать обломки ( для первого фильма Ерлашовой) получил отказ. Мотив- после каждого приезда журналистов обстрелы учащались.

0

25

Vam написал:
это не скальпы, а скальп. Один и тот же фрагмент - "скорпион", снятый с разных ракурсов, но с разными размерами и весом по мнению профи:) Моделируйте... если найдете его по тем координатам, что указаны в отчете.

Это що получается? Один и тот же фрагмент имеет два размера и три веса -10-20-30 кг.
Это что издеваются? Формулы с ошибками украдены у китайцев, а рулетки и весов нема?
А координаты одни и те же.)) Кто стащил "скорпион"? Он точно знает вес и размеры.

Как моделировать?

0

26

Картинка "Скорпиона" из отчета раздела MH17 Crash.
По моему мнению приблизительно размеры 1,4 х 2,4 метра. У кого другое мнение?

http://forumupload.ru/uploads/0016/23/c6/51/27951.png

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

Отредактировано oper (2020-05-03 14:52:59)

0

27

oper написал(а):

три веса -10-20-30 кг.

Меньше 20 кг быть не может, потому что одна обшивка из сплава Д16 будет весить больше 15 кг.

0

28

Простой написал(а):

oper написал(а):

    три веса -10-20-30 кг.

Меньше 20 кг быть не может, потому что одна обшивка из сплава Д16 будет весить больше 15 кг.

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

Площадь 1,4х2,4=3,36 или примерно 3 м^2.
Толшина листов 2 мм или 0.002 м.
Плотность сплава берем 2500 кг/м^3
3х0,002=0,006 м^3
0,006х2500=15 кг
стрингеры - общая длина 6 м, общая ширина 0,2 м, толшина 0,003 м
6х0,2х0,003х2500=9 кг.
итого 15+9=24 кг.

получаем следующую траекторию.

3 D
http://forumupload.ru/uploads/0016/23/c6/51/852395.png
Вид сверху
http://forumupload.ru/uploads/0016/23/c6/51/540680.png
Вид с Запада
http://forumupload.ru/uploads/0016/23/c6/51/596530.png
Вид с Юга
http://forumupload.ru/uploads/0016/23/c6/51/522532.png

0

29

oper написал(а):

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

получаем следующую траекторию.

3 D

Вид сверху

Вид с Запада

Вид с Юга

Вид с Запада - не так что-то! :(

0

30

Простой написал(а):

Вид с Запада - не так что-то! :(

Это вертикальная плоскость Z0Y.
"Стена" с юга на север, фрагмент сносит на север.
можно на нее смотреть и с востока и с запада (зеркальное отражение).

0


Вы здесь » Как сбит MH17: гражданское расследование » Математическое и компьютерное моделирование » Математическое моделирование траектории падения центроплана