 |
 |
 |
 |
 |
 |
 |
|---|
Как сбит MH17: гражданское расследование |
Привет, Гость! Войдите или зарегистрируйтесь.
Вы здесь » Как сбит MH17: гражданское расследование » Математическое и компьютерное моделирование » Математическое моделирование траектории падения центроплана
На картинке в клеточку запредельные скорости, этого неможет быть.
Если сбросить балванку весом 250 тонн с высоты 10000 метров без крылошек под углом к горизонту 70 градусов. то максимальная скорость у земли будет 460 м/сек. http://forumupload.ru/uploads/0016/23/c6/51/804874.png
По всей видимости, наши форумские баллистики должны получить похожие результаты, затем экстраполировать на мн17.
Сбросим с высоты 6400 метров при скорости бомбардировщика 900 км/час ФАБ-250-М54. Бомба приземлится со скоростью 270м/сек. http://forumupload.ru/uploads/0016/23/c6/51/915120.png
Сбросим с высоты 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)
Опер! Разгон боинга по траектории басара на своей адской машинке когда проверишь?
Вначале поросите Простого открыть тему по компьютерному моделированию. В этой теме я опубликую код программы, чтоб все было по честному и потом уже можно будет с вашей помощью разогнать Боинг. Сделаем горку и вперед вниз... Можно будет добавить не только 3D график и вид сверху, но и вид двух сторон, вертикальная проекция.
Открывая отдел для переноса сюда всего, что сможет помочь понять, почему боинг падал по такой странной траектории.
Сх = 0,49 - по всей видимости шар.
Сфера: Cx=0,47
В сущности мидель надо изменить так как бомба падает не вертикально, а под углом и он будет больше чем поперечное сечение цилиндра бомбы.
Достоверно известно только высота,скорость, масса и размеры бомбы недастает только Сх бомбы.
Ралив если вы хотите протестировать программу, то я и предлагаю опубликовать открытый код программы и каждый может определить верно или нет рассчитываются траектории.
взял какую-то условную болванку
А если крылышки приделать к болванки, то она еще медленнее будет падать. У вас видимо "МН17" падает в безвоздушном пространстве и скорости аховские.
Спасибо Простой, что открыли важную тему.
Спасибо Раливу за эту ссылку https://www.webcitation.org/65Tg5GmBM?u … listic.pdf.
В разделе Appendix K окончательного отчета DSB не опубликованы данные баллистических траекторий падения обломков, которые обязаны были опубликовать, примерно такими же как китайцы в 2002 году.
Значит можно создать общую картину падения обломков, китайцы смогли сможем и мы.
Как обещал публикую код программы и любой загрузив Матлаб может использовать этот код.
весь код записан в файл с расширением .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;
endelseif 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);
endend
%---------------------------------------------------------------
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)
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=[];
Хорошо бы дать расшифровку этих параметров.
Да и комментарии к блокам расчёта тоже.
Продолжение.
Отдельно опубликую функцию которая в зависимости от высоты возвращает скорость ветра и метереологический азимут.
Данные взял у 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)
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=[]; все данные стираются.
Замечательно!
А текущие значения в циклах
x = w(1);
y= w(2);
z = w(3);
vx = w(4);
vy = w(5);
vz= w(6);
xdot = vx;
ydot = vy;
zdot=vz;
Замечательно!
А текущие значения в циклах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)
Это к вопросу учёта площади "утюга"
Мидель - видимая площадь объекта перпендикулярная его движению.
Синий вектор - направление движения объекта
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 метров. Есть ли у кого данные по этому обломку?
Хорошо бы послушать Сепаратиста
Не взвешивал 
.
Лежал на мусорнике. Значит кто то туда его принёс. Значит упал кому то во двор. Скорее всего по крайней улице от поля в Рассыпном.
Отредактировано Сепаратист (2020-05-02 02:24:00)
А можно координаты этого "мусорника" ?
В фильме разве не видно? От пятиэтажек в Рассыпном дорога на Петропавловку. На выезде когда то был ж.д. переезд. Теперь там стоит крест... и была небольшая незаконная .свалка. Мэр наверное кипишь и поднял чтобы его за этот мусорник не вздрючили.
Первым оторвался фрагмент под названием "Скорпион".
Анна говорит про фрагмент
2. Вот по Вашей картинке указан вес 0,40. Это сколько - 400 грамм или 40 кг? И то и другое на мой взгляд не катит. Обломок с двумя шпангоутами и дырками примерно размером метр на метр, весит около 5 кг. Всё сугубо субъективно.
По ее словам фрагмент метр на метр и весом около 5 кг. Это "скорпион"?
Добавил в программу возможность отключать влияние ветра, добавленны кнопки отображения 2D графика в плоскости ZОY и ZOX, вертикальные плоскость юг и запад и добавлена возможность вводить мидель фрагмента при учете влияния ветра.
Поэтому выше отредактировал листинг программы.
http://mh17.webtalk.ru/viewtopic.php?id … =15#p86641
Сепаратист написал
Бутблэк,правый это один из тех что Яна отдала в фильме в петропавловский сельсовет.По-видимому в облаке у Вас фотки с того склада. Левый - это скальп-скорпион,который исчез ,не дождавшись сбора обломков.
Как так исчез? Сепаратист вы успели его взвесить и снять размеры?
Вот и данные по скальпам, с большим разбросом по весу. Данные их окончательного очета, раздел - MH17 About the investigation.
Будем пробывать моделировать.
http://mh17.webtalk.ru/viewtopic.php?id … =15#p86641
Как так исчез? Сепаратист вы успели его взвесить и снять размеры?
Вот так - исчез! После падения самолёта до освобождения Дебальцево я в Петропавловке ни разу не был. Там где скальп- скорпион вообще территория бд. Ходить с фотоаппаратом среди позиций военных нужно быть либо чокнутым либо с лицензией на разрешение. Знал командира и его телефон .На просьбу поснимать обломки ( для первого фильма Ерлашовой) получил отказ. Мотив- после каждого приезда журналистов обстрелы учащались.
Vam написал:
это не скальпы, а скальп. Один и тот же фрагмент - "скорпион", снятый с разных ракурсов, но с разными размерами и весом по мнению профи:) Моделируйте... если найдете его по тем координатам, что указаны в отчете.
Это що получается? Один и тот же фрагмент имеет два размера и три веса -10-20-30 кг.
Это что издеваются? Формулы с ошибками украдены у китайцев, а рулетки и весов нема?
А координаты одни и те же.)) Кто стащил "скорпион"? Он точно знает вес и размеры.
Как моделировать?
Картинка "Скорпиона" из отчета раздела MH17 Crash.
По моему мнению приблизительно размеры 1,4 х 2,4 метра. У кого другое мнение?
Зная площадь фрагмента и толшину обшивки, а также примерный вес стрингеров, можно будет определить вес всего фрагмента.
Отредактировано oper (2020-05-03 14:52:59)
три веса -10-20-30 кг.
Меньше 20 кг быть не может, потому что одна обшивка из сплава Д16 будет весить больше 15 кг.
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 
Вид сверху
Вид с Запада
Вид с Юга
По моим приблезительным подсчетом вес фрагмента получился 24 кг.
получаем следующую траекторию.
3 D
Вид сверху
Вид с Запада
Вид с Юга
Вид с Запада - не так что-то! 
Вид с Запада - не так что-то!
Это вертикальная плоскость Z0Y.
"Стена" с юга на север, фрагмент сносит на север.
можно на нее смотреть и с востока и с запада (зеркальное отражение).
oper написал(а):
Вид сверху
Время падения 12 минут.
Траектория в точности повторяет траекторию падения в виде "клюшки" на калькуляторе uschen-а, а это значит , oper, что у тебя "скорпион" падает со скоростью РАВНОЙ скорости Ветра.Твоя программа не про то как падают обломки. Твоя программа про то КАК ПАДАЕТ Ветер.
У тебя есть "тяжёлый" Ветер, есть "средний" Ветер, и есть "лёгкий" Ветер...Но ты же понимаешь что это АБСУРД ? Ветер не может падать.
Высота от 0 до 10 км разделена НА СЛОИ. У каждого воздушного слоя своя скорость и направление.
А зачем так много слов?
Вам необходимо вместо слов написать формулу с учетом ветра.
Чуть выше есть опубликованный листинг программы (код) вот там и укажите мою ошибку.
Вы здесь » Как сбит MH17: гражданское расследование » Математическое и компьютерное моделирование » Математическое моделирование траектории падения центроплана