Делфи, моделирование физических процессов. Исходник.

[Алексей_2012]

Ребят, всем привет, есть у кого-нибудь любые программы, относящиеся к предмету "физика" с исходным кодом?

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

Просто любая программа, сам искал - не нашел

[Алексей_2012]

Вот нашел моделирование броска мяча, только не понял что это за параметры...

 IsRunning: boolean;
  Radius:    Cardinal=5;  
  k:         single=0.5; 
  R:         single=0.97; 
  m:         single=100;    
  g:         single=9.8; // я написал, но все ускорилось в разы, было 1.3
  fx,fy:     single;      
  x:         single=20;   
  y:         single=300;

вот полный код, достаточно маленький, и какая строчка отвечает за расчет и физику?

unit Unit1;
 
interface
 
uses
  Windows, Messages, Classes, Controls, Forms, Math;
 
type
  TForm1 = class(TForm)
    procedure FormMouseMove(Sender: TObject; Shift: TShiftState; X, Y: Integer);
    procedure FormMouseDown(Sender: TObject; Button: TMouseButton; Shift: TShiftState; X, Y: Integer);
  end;
 
var
  Form1:     TForm1;
  IsRunning: boolean;
  Radius:    Cardinal=5;  //ГђГ*äèóñ
  k:         single=0.5;  //ÊîýôôèöåГ*ГІ óïðóãîñòè Г±ГІГҐГ*
  R:         single=0.97; //ÑîïðîòèâëåГ*ГЁГҐ
  m:         single=100;    //Г¬Г*Г±Г±Г*
  g:         single=9.8;  //óñêîðåГ*ГЁГҐ ñâîáîäГ*îãî ГЇГ*äåГ*ГЁГї
  fx,fy:     single;      //Âåêòîð ñèëû
  x:         single=20;   //ÊîîðäèГ*Г*ГІГ» ГІГҐГ«Г*
  y:         single=300;
 
implementation
 
{$R *.dfm}
 
procedure TimerProc(hwnd,uMsg,idEvent,dwTime: LongInt);stdcall;
begin
  BitBlt(Form1.Canvas.Handle,0,0,Form1.Width,Form1.Height,0,0,0,WHITENESS);
  x :=x+fx;
  y :=y+fy;
  fx:=fx*R;
  fy:=fy*R+m*g;
  if Cardinal(round(x)-Radius)>Cardinal(Form1.ClientWidth-Radius*2) then
    fx:=-fx*k;
  if Cardinal(round(y)-Radius)>Cardinal(Form1.ClientHeight-Radius*2) then
    fy:=-fy*k;
  x:=min(max(x,Radius),Form1.ClientWidth-Radius);
  y:=min(max(y,Radius),Form1.ClientHeight-Radius);
  Ellipse(Form1.Canvas.Handle,round(x-Radius),round(y-Radius),round(x+Radius),round(y+Radius));
end;
 
procedure TForm1.FormMouseMove(Sender: TObject; Shift: TShiftState; X, Y: Integer);
begin
  if not IsRunning then
  begin
    BitBlt(Form1.Canvas.Handle,0,0,Form1.Width,Form1.Height,0,0,0,WHITENESS);
    fx:=(x-unit1.x)*0.1;
    fy:=(y-unit1.y)*0.1;
    MoveToEx(Form1.Canvas.Handle,round(unit1.x),round(unit1.y),0);
    LineTo(Form1.Canvas.Handle,round(x),round(y));
  end;
end;
 
procedure TForm1.FormMouseDown(Sender: TObject; Button: TMouseButton; Shift: TShiftState; X, Y: Integer);
begin
  if IsRunning then
  begin
    KillTimer(Form1.Handle,0);
    unit1.x:=20;
    unit1.y:=300;
  end
  else
    SetTimer(Form1.Handle,0,50,@TimerProc);
  IsRunning:=not IsRunning;
end;
 
end.

[phomm]

А что тут думать ? Таймерная процедура почти вся отвечает за расчёт, кроме первой и последней строки, это прорисовка. 2 ифки это очень похоже на отражение мяча от стенок, ещё там 2 строки, по сути тоже ифки на функциях min и max - некие ограничения, наверное невылета мяча за форму, следовательно весь расчёт вообще только в 4 строках в начале - приращение координат на текущую дельту и модификация дельты под действием тяжести и сопротивления воздуха (как я понял R). В маусдауне запуск мяча с фикс координат, в маусмуве задание направления (как я понял). Прогу не смотрел, только читнул код из поста.

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