файл Unit1. h

// ---------------------------------------------------------------------------

#ifndef Unit1H

#define Unit1H

// ---------------------------------------------------------------------------

#include <Classes. hpp>

#include <Controls. hpp>

#include <StdCtrls. hpp>

#include <Forms. hpp>

#include <ComCtrls. hpp>

#include <gl\gl. h>

#include <gl\glu. h>

#include <gl\glaux. h>

// ---------------------------------------------------------------------------

struct tochka

{

GLfloat x,y; // текущие координаты точки

GLfloat x0,y0; // начальные координаты точки

double vx,vy; // проекции скоростей в текущее время

double m; // масса точки

double Fx,Fy; // проекции сил в текущее время

};

class TForm1: public TForm

{

__published: // IDE-managed Components

TMemo *Memo1;

TButton *Button1;

TButton *Button2;

TEdit *Edit1;

TEdit *Edit3;

TEdit *Edit4;

TEdit *Edit5;

TEdit *Edit6;

TLabel *Label2;

TLabel *Label3;

TLabel *Label4;

TLabel *Label5;

TLabel *Label7;

TEdit *Edit2;

TLabel *Label6;

void __fastcall FormCreate (TObject *Sender); // создание формы

void __fastcall Start (TObject *Sender); // кнопка start

void __fastcall Stop (TObject *Sender); // кнопка stop

void __fastcall Edit1KeyPress (TObject *Sender,

char &key); // проверка ввода данных для вещественных чисел

void __fastcall Edit2KeyPress (TObject *Sender,

char &Key); // проверка ввода данных для целых чисел

private: // User declarations

int PixelFormat; //

HDC hdc; //

HGLRC hrc; // переменные для рисования

GLuint texture1; //

tochka *t; // массив точек

int dlina; // длина нити

int kt; // количество точек

float shag; // шаг разбиения

float massa; // масса нити

float dconst; // константа демпфирования

float h; // шаг дискретизации

float ugol; // угол наклона

public: // User declarations

__fastcall TForm1 (TComponent* Owner);

void __fastcall SetPixelFormatDescriptor (); // установка дескрипта

void __fastcall DrawFunc (TObject*, bool&); // функция рисования

void Sila_RS (); // сила растяжения - сжатия

void Sila_Tiagesty (); // сила тяжести

void Sila_Soprotivlenia (); // сила сопротивления среды

void Sili (); // силы в целом и закреление точек

void Massa_tochek (); // вычисление массы каждой точки

void Pereshagivanie (); // нахождение новых координат и скоростей

};

// ---------------------------------------------------------------------------

extern PACKAGE TForm1 *Form1;

// ---------------------------------------------------------------------------

#endif

файл Unit1. cpp

// ---------------------------------------------------------------------------

#include <vcl. h>

#pragma hdrstop

#include "Unit1. h"

#include <math. h>

// ---------------------------------------------------------------------------

#pragma package (smart_init)

#pragma resource "*. dfm"

TForm1 *Form1;

bool stop = true;

// ---------------------------------------------------------------------------

__fastcall TForm1:: TForm1 (TComponent* Owner)

: TForm (Owner)

{

Application->OnIdle = DrawFunc; // по-умолчанию

}

// ---------------------------------------------------------------------------

void __fastcall TForm1:: DrawFunc (TObject *ob, bool &done) // функция рисования

{

int i;

GLfloat x,y;

done = false;

// gluLookAt (1,0,1,0,0,0,0,1,0);

if (! stop)

{

Sili ();

New_point ();

glClear (GL_COLOR_BUFFER_BIT|GL_DEPTH_BUFFER_BIT);

// рисование

glBegin (GL_LINE_STRIP);

glColor3d (0,0,0);

for (i = 0; i < kt; i++)

{

x = t [i]. x;

y = t [i]. y;

glVertex2f (x,y);

}

glEnd ();

glFlush ();

SwapBuffers (hdc);

}

}

// --------------------------------------------------------------------------

void __fastcall TForm1:: SetPixelFormatDescriptor () // установка дескрипта

{

PIXELFORMATDESCRIPTOR pfd = {

sizeof (PIXELFORMATDESCRIPTOR),

1,PFD_DRAW_TO_WINDOW | PFD_SUPPORT_OPENGL | PFD_DOUBLEBUFFER,

PFD_TYPE_RGBA,

24,0,0,0,0,0,0,0,0,0,0,0,0,0,32,0,0,PFD_MAIN_PLANE,

0,0,0,};

PixelFormat = ChoosePixelFormat (hdc, &pfd);

SetPixelFormat (hdc, PixelFormat, &pfd);

}

// -----------------------------------------------------------------------------

void __fastcall TForm1:: FormCreate (TObject *Sender) // создание формы

{

hdc = GetDC (Memo1->Handle); // рисуем на объекте Мэмо

SetPixelFormatDescriptor ();

hrc = wglCreateContext (hdc);

if (hrc == NULL)

ShowMessage (": - ) ~ hrc == NULL");

if (wglMakeCurrent (hdc, hrc) == false)

ShowMessage ("Could not MakeCurrent");

glClearColor (1,1,1,0);

}

// -----------------------------------------------------------------------------

void TForm1:: Sili () // силы

{

int i;

for (i = 0; i < kt; i++)

{

t [i]. Fx = 0;

t [i]. Fy = 0;

}

Sila_RS ();

Sila_Tiagesty ();

Sila_Soprotivlenia ();

t [0]. Fx = 0;

t [0]. Fy = 0;

t [kt]. Fx = 0;

t [kt]. Fy = 0;

}

// -------------------------------------------------------------------------

void TForm1:: Sila_RS () // силы растяжения-сжатия между соседними точками

{

int i,k,l,p,ks;

double sqr;

ks = StrToInt (Edit2->Text);

for (i = 0; i < kt; i++)

{ if (i == 0) l = 0;

else l = i - 1;

if (i == kt-1) p = kt - 1;

else p = i + 1;

for (k = l; k <= p; k++)

{

if (k == i) continue;

sqr = sqrt ( (t [k]. x0-t [i]. x0) * (t [k]. x0-t [i]. x0) + (t [k]. y0-t [i]. y0) * (t [k]. y0-t [i]. y0)) /sqrt ( (t [k]. x-t [i]. x) * (t [k]. x-t [i]. x) + (t [k]. y-t [i]. y) * (t [k]. y-t [i]. y));

t [i]. Fx += ks* ( (t [k]. x - t [i]. x) - sqr * (t [k]. x - t [i]. x));

t [i]. Fy += ks* ( (t [k]. y - t [i]. y) - sqr * (t [k]. y - t [i]. y));

}

}

}

// ------------------------------------------------------------------------

void TForm1:: Sila_Tiagesty () // сила тяжести

{

int i;

for (i = 0; i < kt; i++)

t [i]. Fy - = 9.81*t [i]. m;

}

// ------------------------------------------------------------------------

void TForm1:: Sila_Soprotivlenia () // сила сопротивления среды

{

int i;

for (i = 0; i < kt; i++)

{

t [i]. Fx - = dconst*t [i]. vx;

t [i]. Fy - = dconst*t [i]. vy;

}

}

// ------------------------------------------------------------------------

void TForm1:: New_point () // Метод Эйлера

{

int i;

for (i = 0; i < kt; i++)

{

t [i]. x += t [i]. vx*h;

t [i]. y += t [i]. vy*h;

t [i]. vx += t [i]. Fx/t [i]. m*h;

t [i]. vy += t [i]. Fy/t [i]. m*h;

}

}

// ---------------------------------------------------------------------------

void TForm1:: Massa_tochek () // перераспределение массы

{

int i;

for (i = 0; i < kt; i++)

t [i]. m = massa/kt;

}

// --------------------------------------------------------

void __fastcall TForm1:: Start (TObject *Sender) // Начало моделирования

{

// gluLookAt (0,0,0,0,0,1,0,0,0);

int i;

dlina = StrToInt (Edit3->Text);

massa = StrToFloat (Edit4->Text);

shag = StrToFloat (Edit5->Text);

dconst = StrToFloat (Edit6->Text);

h = StrToFloat (Edit1->Text);

ugol = 90;

kt = dlina/shag;

t = new tochka [kt];

stop = false;

for (i = 0; i < kt; i++)

{

t [i]. x0 = t [i]. x = 1.2* ( (GLfloat) i + (GLfloat) kt/2) / (2*kt);

t [i]. y0 = t [i]. y = 0.07;

t [i]. vx = 0;

t [i]. vy = 0;

}

Massa_tochek ();

glColor3d (1,1,1);

Button1->Enabled = false;

Button2->Enabled = true;

}

// -------------------------------------------------------------------------------

void __fastcall TForm1:: Stop (TObject *Sender) // Окончание моделирования

{

stop = true;

Button1->Enabled = true;

Button2->Enabled = false;

delete t;

}

// -----------------------------------------------------------------------

void __fastcall TForm1:: Edit1KeyPress (TObject *Sender, char &Key)

{ // проверка дробных значений

if ( (Key >= '0') && (Key <= '9')) return;

if (Key == DecimalSeparator) return;

if (Key == VK_BACK) // клавиша <Backspace>

return;

Key=0;

}

// ---------------------------------------------------------------------------

void __fastcall TForm1:: Edit2KeyPress (TObject *Sender, char &Key)

{ // проверка целых значений

if ( (Key >= '0') && (Key <= '9')) return;

if (Key == VK_BACK)

return;

Key=0;

}

// ---------------------------------------------------------------------------

Размещено на Allbest.ru