Периодические граничные условия
Материал из Department of Theoretical and Applied Mechanics
Виртуальная лаборатория>Периодические граничные условия
Содержание
Курсовой проект по механике дискретных сред[править]
- разработчик Теницкая Татьяна
- руководитель Кузькин Виталий
Краткое описание[править]
- Метод периодических граничных условий был разработан для решения задач теории жидкостей и плотных газов. Он состоит в том,что вокруг расчетной области строятся ее «образы» с актуальным положением частиц. И частицы «реальной» области взаимодействуют с частицами в «образе», а если частица пересекает границу расчетной области, она появляется с другой стороны.
- В теореме Нетер утверждается, что каждой непрерывной симметрии физической системы соответствует некоторый закон сохранения:
- однородности времени соответствует закон сохранения энергии,
- однородности пространства соответствует закон сохранения импульса,
- изотропии пространства соответствует закон сохранения момента импульса,
- калибровочной симметрии соответствует закон сохранения электрического заряда и т. д.
- Но для классической системы частиц с периодическими условиями сохранение момента импульса нарушается. Этот эффект наглядно проиллюстрирован в данной курсовой работе.
Цель проекта[править]
- Визуализация системы частиц с периодическими граничными условиями.
- Построение графиков зависимости кинетического момента от времени для одной частицы, двух частиц, многих частиц.
Математическая модель[править]
Граничные условия:
если , то
если , то
если , то
если , то
Где x и у - это координаты частицы, а w и h - ширина и длина окна соответственно.
Кинетический момент вычисляется по формуле:
Результаты[править]
Текст программы на языке JavaScript (разработчик Теницкая Татьяна):
Файл "One_Ball.js"
1 function main_particle_1() {
2
3 var ctx = canvas_particle_1.getContext("2d");
4 var w_1 = canvas_particle_1.width;
5 var h_1 = canvas_particle_1.height;
6 var w = w_1/3;
7 var h = h_1/3;
8
9 var r = 7;
10 var N = 1;
11 var rx,ry; //Содержит координаты частиц по оси Х //Содержит координаты частиц по оси У
12
13 var rxI0,rxI1,rxI2,rxI3,rxI4,rxI5,rxI6,rxI7;
14 var ryI0,ryI1,ryI2,ryI3,ryI4,ryI5,ryI6,ryI7;
15
16 var vx,vy;
17 var L;
18
19 dt = 0.3;
20 var alfa = 45 / 180 * Math.PI;
21 //0.04380741654852822
22 var v0 = 5;
23 vx = v0 * Math.cos(alfa);
24 vy = v0 * Math.sin(alfa);
25
26
27 rx = Math.random() * (w - 10);
28 ry = Math.random() * (h - 10);
29
30 var steps = 0;
31
32
33 function step()
34 {
35 tick();
36 draw();
37 }
38 var vGraph = new TM_graph( // определить график
39 "#vGraph", // на html-элементе #vGraph
40 200, // сколько шагов по оси "x" отображается
41 -1, 1,0.2); // мин. значение оси Y, макс. значение оси Y, шаг по оси Y
42
43 function tick()
44 {
45
46 steps += 1;
47
48 rx += vx*dt;
49 ry += vy*dt;
50
51 if (rx >= w) {rx = rx - w;}
52 if (rx <= 0) {rx = rx + w;}
53 if (ry >= h) {ry = ry - h;}
54 if (ry <= 0) {ry = ry + h;}
55
56 rxI0 = rx + w;
57 ryI0 = ry;
58
59 rxI3 = rx + w;
60 ryI3 = ry + h;
61
62 rxI6 = rx + w;
63 ryI6 = ry + 2*h;
64
65 rxI1 = rx + 2*w;
66 ryI1 = ry;
67
68 rxI4 = rx + 2*w;
69 ryI4 = ry + h;
70
71 rxI7 = rx + 2*w;
72 ryI7 = ry + 2*h;
73
74 rxI2 = rx;
75 ryI2 = ry + h;
76
77 rxI5 = rx;
78 ryI5 = ry + 2*h;
79
80
81 L = (rx * vy - ry * vx)/500;
82 if (steps % 1 == 0) vGraph.graphIter(steps, L); // подать данные на график
83 console.log(L);
84
85 }
86
87
88 function draw()
89 {
90 ctx.clearRect(0, 0, w_1 , h_1); // очистить экран
91
92 var x = rx;
93 var y = ry;
94 ctx.beginPath();
95 ctx.fillStyle = "#00008B";
96 ctx.arc(x, y, r , 0, 2 * Math.PI, false);
97 ctx.closePath();
98 ctx.fill();
99
100 var xI0 = rxI0;
101 var yI0 = ryI0;
102 ctx.beginPath();
103 ctx.fillStyle = "#00008B";
104 ctx.arc(xI0, yI0, r , 0, 2 * Math.PI, false);
105 ctx.closePath();
106 ctx.fill();
107 var xI1 = rxI1;
108 var yI1 = ryI1;
109 ctx.beginPath();
110 ctx.fillStyle = "#00008B";
111 ctx.arc(xI1, yI1, r , 0, 2 * Math.PI, false);
112 ctx.closePath();
113 ctx.fill();
114 var xI2 = rxI2;
115 var yI2 = ryI2;
116 ctx.beginPath();
117 ctx.fillStyle = "#00008B";
118 ctx.arc(xI2, yI2, r , 0, 2 * Math.PI, false);
119 ctx.closePath();
120 ctx.fill();
121 var xI3 = rxI3;
122 var yI3 = ryI3;
123 ctx.beginPath();
124 ctx.fillStyle = "#00008B";
125 ctx.arc(xI3, yI3, r , 0, 2 * Math.PI, false);
126 ctx.closePath();
127 ctx.fill();
128 var xI4 = rxI4;
129 var yI4 = ryI4;
130 ctx.beginPath();
131 ctx.fillStyle = "#00008B";
132 ctx.arc(xI4, yI4, r , 0, 2 * Math.PI, false);
133 ctx.closePath();
134 ctx.fill();
135 var xI5 = rxI5;
136 var yI5 = ryI5;
137 ctx.beginPath();
138 ctx.fillStyle = "#00008B";
139 ctx.arc(xI5, yI5, r , 0, 2 * Math.PI, false);
140 ctx.closePath();
141 ctx.fill();
142 var xI6 = rxI6;
143 var yI6 = ryI6;
144 ctx.beginPath();
145 ctx.fillStyle = "#00008B";
146 ctx.arc(xI6, yI6, r , 0, 2 * Math.PI, false);
147 ctx.closePath();
148 ctx.fill();
149 var xI7 = rxI7;
150 var yI7 = ryI7;
151 ctx.beginPath();
152 ctx.fillStyle = "#00008B";
153 ctx.arc(xI7, yI7, r , 0, 2 * Math.PI, false);
154 ctx.closePath();
155 ctx.fill();
156
157
158 ctx.beginPath(); // начать рисование
159 ctx.fillStyle="#000000";
160 ctx.moveTo(w, 0); // переместить "карандаш" в точку
161 ctx.lineTo(w, h_1); // нарисовать "карандашом" линию до точки
162 ctx.stroke();
163
164 ctx.beginPath();
165 ctx.fillStyle="#000000";
166 ctx.moveTo(2 * w, 0);
167 ctx.lineTo(2 * w, h_1);
168 ctx.stroke();
169
170 ctx.beginPath();
171 ctx.fillStyle="#000000";
172 ctx.moveTo(0, h);
173 ctx.lineTo(w_1, h);
174 ctx.stroke();
175
176 ctx.beginPath();
177 ctx.fillStyle="#000000";
178 ctx.moveTo(0, 2 * h);
179 ctx.lineTo(w_1, 2 * h);
180 ctx.stroke();
181 }
182
183 setInterval(step, 1000/20);
184 }
Текст программы на языке JavaScript (разработчик Теницкая Татьяна):
Файл "TwoBalls.js"
1 function main_particle_1() {
2
3 var ctx_11 = canvas_particle_11.getContext("2d");
4 var ctx_12 = canvas_particle_12.getContext("2d");
5 var ctx_13 = canvas_particle_13.getContext("2d");
6 var ctx_21 = canvas_particle_21.getContext("2d");
7 var ctx_22 = canvas_particle_22.getContext("2d");
8 var ctx_23 = canvas_particle_23.getContext("2d");
9 var ctx_31 = canvas_particle_31.getContext("2d");
10 var ctx_32 = canvas_particle_32.getContext("2d");
11 var ctx_33 = canvas_particle_33.getContext("2d");
12
13
14 var width = canvas_particle_11.width;
15 var height = canvas_particle_11.height;
16
17 var rx_1, ry_1, rx_2, ry_2; // координаты шара
18 var v0 = 4, vx_1, vy_1, vx_2, vy_2;
19 var r = 7; // радиус шара
20 var dt_1 = 0; dt_2 = 0;
21 var x0_1 = width/2, x0_2 = width/2, y0_1 = 0, y0_2 = height;
22 var steps = 0;
23 var val = 1;
24 var alfa = 45/180*Math.PI;
25
26 function step()
27 {
28 tick();
29 draw();
30 }
31
32 var vGraph = new TM_graph( // определить график
33 "#vGraph", // на html-элементе #vGraph
34 250, // сколько шагов по оси "x" отображается
35 -1,1,0.2);
36
37 function tick()
38 { // вычисление новой позиции шар
39 vx_1 = v0 * Math.cos(alfa);
40 vy_1 = v0 * Math.sin(alfa);
41 vx_2 = -v0 * Math.cos(alfa);
42 vy_2 = -v0 * Math.sin(alfa);
43
44 steps+=1;
45 if (val==1)
46 {
47 rx_1 = vx_1*dt_1 + x0_1 + r - 2.7; //
48 ry_1 = vy_1*dt_1 + y0_1 + r - 2.7;
49 dt_1+= 1;
50 rx_2 = vx_2*dt_2 + x0_2 - r + 2.7; //
51 ry_2 = vy_2*dt_2 + y0_2 - r + 2.7;
52 dt_2+= 1;
53
54
55 }
56 else
57 {
58 vy_1 = -vy_1;
59 vx_1 = -vx_1;
60 rx_1 = vx_1*dt_1 + width - r + 2.7;
61 ry_1 = vy_1*dt_1 + x0_1 - r + 2.7;
62 dt_1+= 1;
63 vy_2 = -vy_2;
64 vx_2 = -vx_2;
65 rx_2 = vx_2*dt_2 - 2.7 + r; //
66 ry_2 = vy_2*dt_2 + x0_2 - 2.7 + r;
67 dt_2+= 1;
68 }
69
70 if (rx_2 < 0 || rx_1 - r> width)
71 {
72 //alert(rx_2);
73 val = 0;
74 dt_1 = 0;
75 dt_2 = 0;
76 }
77
78 if (ry_1 < 0 || ry_2 - r > height)
79 {
80 val = 1;
81 dt_1=0;
82 dt_2 = 0;
83 }
84
85
86 var L = ((rx_1 * vy_1 - ry_1 * vx_1)+(rx_2 * vy_2 - ry_2 * vx_2))/376.2;
87 console.log(L);
88
89 if (dt_1 % 0.5 == 0) vGraph.graphIter(steps, L); // подать данные на график
90
91
92
93 }
94
95 function draw() // рисование шара
96 {
97
98 ctx_11.fillStyle= "#DC143C" ; // цвет закраски
99 ctx_11.clearRect(0, 0, width, height); // очистить экран
100 ctx_11.arc(rx_1, ry_1, r, 0, 2 * Math.PI);
101 ctx_11.fill();
102 ctx_11.beginPath();
103
104 ctx_11.fillStyle = "#00008B";
105 ctx_11.arc(rx_2, ry_2, r, 0, 2 * Math.PI);
106 ctx_11.fill();
107 ctx_11.beginPath();
108
109 ctx_12.fillStyle = "#DC143C"; // цвет закраски
110 ctx_12.clearRect(0, 0, width, height); // очистить экран
111 ctx_12.arc(rx_1, ry_1, r, 0, 2 * Math.PI);
112 ctx_12.fill();
113 ctx_12.beginPath();
114
115 ctx_12.fillStyle = "#00008B";
116 ctx_12.arc(rx_2, ry_2, r, 0, 2 * Math.PI);
117 ctx_12.fill();
118 ctx_12.beginPath();
119
120 ctx_13.fillStyle = "#DC143C"; // цвет закраски
121 ctx_13.clearRect(0, 0, width, height); // очистить экран
122 ctx_13.arc(rx_1, ry_1, r, 0, 2 * Math.PI);
123 ctx_13.fill();
124 ctx_13.beginPath();
125
126 ctx_13.fillStyle = "#00008B";
127 ctx_13.arc(rx_2, ry_2, r, 0, 2 * Math.PI);
128 ctx_13.fill();
129 ctx_13.beginPath();
130
131 ctx_21.fillStyle = "#DC143C"; // цвет закраски
132 ctx_21.clearRect(0, 0, width, height); // очистить экран
133 ctx_21.arc(rx_1, ry_1, r, 0, 2 * Math.PI);
134 ctx_21.fill();
135 ctx_21.beginPath();
136
137 ctx_21.fillStyle = "#00008B";
138 ctx_21.arc(rx_2, ry_2, r, 0, 2 * Math.PI);
139 ctx_21.fill();
140 ctx_21.beginPath();
141
142 ctx_22.fillStyle = "#DC143C"; // цвет закраски
143 ctx_22.clearRect(0, 0, width, height); // очистить экран
144 ctx_22.arc(rx_1, ry_1, r, 0, 2 * Math.PI);
145 ctx_22.fill();
146 ctx_22.beginPath();
147
148 ctx_22.fillStyle = "#00008B";
149 ctx_22.arc(rx_2, ry_2, r, 0, 2 * Math.PI);
150 ctx_22.fill();
151 ctx_22.beginPath();
152
153 ctx_23.fillStyle = "#DC143C"; // цвет закраски
154 ctx_23.clearRect(0, 0, width, height); // очистить экран
155 ctx_23.arc(rx_1, ry_1, r, 0, 2 * Math.PI);
156 ctx_23.fill();
157 ctx_23.beginPath();
158
159 ctx_23.fillStyle = "#00008B";
160 ctx_23.arc(rx_2, ry_2, r, 0, 2 * Math.PI);
161 ctx_23.fill();
162 ctx_23.beginPath();
163
164 ctx_31.fillStyle = "#DC143C"; // цвет закраски
165 ctx_31.clearRect(0, 0, width, height); // очистить экран
166 ctx_31.arc(rx_1, ry_1, r, 0, 2 * Math.PI);
167 ctx_31.fill();
168 ctx_31.beginPath();
169
170 ctx_31.fillStyle = "#00008B";
171 ctx_31.arc(rx_2, ry_2, r, 0, 2 * Math.PI);
172 ctx_31.fill();
173 ctx_31.beginPath();
174
175 ctx_32.fillStyle = "#DC143C"; // цвет закраски
176 ctx_32.clearRect(0, 0, width, height); // очистить экран
177 ctx_32.arc(rx_1, ry_1, r, 0, 2 * Math.PI);
178 ctx_32.fill();
179 ctx_32.beginPath();
180
181 ctx_32.fillStyle = "#00008B";
182 ctx_32.arc(rx_2, ry_2, r, 0, 2 * Math.PI);
183 ctx_32.fill();
184 ctx_32.beginPath();
185
186 ctx_33.fillStyle = "#DC143C"; // цвет закраски
187 ctx_33.clearRect(0, 0, width, height); // очистить экран
188 ctx_33.arc(rx_1, ry_1, r, 0, 2 * Math.PI);
189 ctx_33.fill();
190 ctx_33.beginPath();
191
192 ctx_33.fillStyle = "#00008B";
193 ctx_33.arc(rx_2, ry_2, r, 0, 2 * Math.PI);
194 ctx_33.fill();
195 ctx_33.beginPath();
196 }
197
198 setInterval(step, 1000/30); // функция step будет запускаться 30 раз в секунду (30 раз / 1000 мс)
199
200 }
Текст программы на языке JavaScript (разработчик Теницкая Татьяна):
Файл "NNN.js"
1 function main_particle_1() {
2
3 var ctx = canvas_particle_1.getContext("2d");
4 var w_1 = canvas_particle_1.width;
5 var h_1 = canvas_particle_1.height;
6 var w = w_1/3;
7 var h = h_1/3;
8
9 var r = 7;
10 var N = 10;
11 rx = new Array(); //Содержит координаты частиц по оси Х
12 ry = new Array(); //Содержит координаты частиц по оси У
13 rxI0 = new Array();rxI1 = new Array();rxI2 = new Array();rxI3 = new Array();rxI4 = new Array();rxI5 = new Array();rxI6 = new Array();rxI7 = new Array();
14 ryI0 = new Array();ryI1= new Array();ryI2= new Array();ryI3= new Array();ryI4= new Array();ryI5= new Array();ryI6= new Array();ryI7= new Array();
15 vx = new Array(); //Скорость по Х
16 vy = new Array(); //Скорость по У
17 vx_p = new Array();vx_m = new Array();vy_p = new Array();vy_m = new Array();
18 var distR0,distR2,distR3,distR;
19
20 L = new Array();
21 var D = 5;
22 var a0 = w * 0.1;
23 var c = 72 * D/r;
24 var W = Math.sqrt(c);
25 var T = 2 * Math.PI/W;
26 dt = 0.03;
27 //0.04380741654852822
28
29 for (i = 0; i < N/2; ++i)
30 {
31 vx_p[i] = Math.random()*15;
32 vy_p[i] = Math.random()*15;
33 vx_m[i] = -vx_p[i];
34 vy_m[i] = -vy_p[i];
35 vx = vx_p.concat(vx_m);
36 vy = vy_p.concat(vy_m);
37 }
38
39 var kol = 1;
40 var prx = 0;
41 var pry = 0;
42 var dist = 0;
43 var f;
44 rx[0] = Math.random() * (w - 10);
45 ry[0] = Math.random() * (h - 10);
46 while (kol < N)
47 {
48 f = 0;
49 prx = Math.random() * (w - 10);
50 pry = Math.random() * (h - 10);
51 for (i = 0; i < kol; i++)
52 {
53 dist = Math.pow((rx[i] - prx),2)+Math.pow((ry[i] - pry),2);
54 if (dist < Math.pow(2 * a0,2))
55 {
56 f = f+1;
57 }
58 }
59 if (f==0)
60 {
61 rx[kol] = prx;
62 ry[kol] = pry;
63 kol++;
64 }
65 }
66 var steps = 0;
67
68
69 function step()
70 {
71 tick();
72 draw();
73 }
74 var vGraph = new TM_graph( // определить график
75 "#vGraph", // на html-элементе #vGraph
76 1000, // сколько шагов по оси "x" отображается
77 -1, 1,0.2); // мин. значение оси Y, макс. значение оси Y, шаг по оси Y
78
79 function tick()
80 {
81 for (i = 0; i < N ; i++)
82 {
83
84 steps += 1;
85
86 rx[i] += vx[i]*dt;
87 ry[i] += vy[i]*dt;
88
89 if (rx[i] >= w) {rx[i] = rx[i] - w;}
90 if (rx[i] <= 0) {rx[i] = rx[i] + w;}
91 if (ry[i] >= h) {ry[i] = ry[i] - h;}
92 if (ry[i] <= 0) {ry[i] = ry[i] + h;}
93
94 rxI0[i] = rx[i] + w;
95 ryI0[i] = ry[i];
96
97 rxI3[i] = rx[i] + w;
98 ryI3[i] = ry[i] + h;
99
100 rxI6[i] = rx[i] + w;
101 ryI6[i] = ry[i] + 2*h;
102
103 rxI1[i] = rx[i] + 2*w;
104 ryI1[i] = ry[i];
105
106 rxI4[i] = rx[i] + 2*w;
107 ryI4[i] = ry[i] + h;
108
109 rxI7[i] = rx[i] + 2*w;
110 ryI7[i] = ry[i] + 2*h;
111
112 rxI2[i] = rx[i];
113 ryI2[i] = ry[i] + h;
114
115 rxI5[i] = rx[i];
116 ryI5[i] = ry[i] + 2*h;
117
118
119
120 //раcсчет расстояний
121
122 for (j = 0; j < N ; j++)
123 {
124 if (i != j)
125 {
126 distR0 = Math.sqrt(Math.pow((rx[i]-rxI0[j]),2)+Math.pow((ry[i]-ryI0[j]),2));
127 distR2 = Math.sqrt(Math.pow((rx[i]-rxI2[j]),2)+Math.pow((ry[i]-ryI2[j]),2));
128 distR3 = Math.sqrt(Math.pow((rx[i]-rxI3[j]),2)+Math.pow((ry[i]-ryI3[j]),2));
129 distR = Math.sqrt(Math.pow((rx[i]-rx[j]),2)+Math.pow((ry[i]-ry[j]),2));
130
131 //раcсчет сил
132 a = 4 * r;
133 if (distR0 < 2*r)
134 {
135 var mod_fR0,fxR0,fyR0;
136
137 mod_fR0 = 12*D*((Math.pow(r/distR0),13)-Math.pow((r/distR0),7))/a;
138 //console.log(mod_fR0);
139 if (mod_fR0 >=0)
140 {
141 fxR0 = mod_fR0*(rx[i]-rxI0[j]);
142 fyR0 = mod_fR0*(ry[i]-ryI0[j]);
143 vx[i]+=fxR0*dt;
144 vx[j]-=fxR0*dt;
145 vy[i]+=fyR0*dt;
146 vy[j]-=fyR0*dt;
147 }
148 }
149 if (distR2 < 2*r)
150 {
151 var mod_fR2,fxR2,fyR2;
152
153 mod_fR2 = 12*D*((Math.pow(r/distR2),13)-Math.pow((r/distR2),7))/a;
154
155 if (mod_fR2 >=0)
156 {
157 fxR2 = mod_fR2*(rx[i]-rxI2[j]);
158 fyR2 = mod_fR2*(ry[i]-ryI2[j]);
159 vx[i]+=fxR2*dt;
160 vx[j]-=fxR2*dt;
161 vy[i]+=fyR2*dt;
162 vy[j]-=fyR2*dt;
163 }
164 }
165 if (distR3 < 2*r)
166 {
167 var mod_fR3,fxR3,fyR3;
168
169 mod_fR3 = 12*D*((Math.pow(r/distR3),13)-Math.pow((r/distR3),7))/a;
170
171 if (mod_fR3 >=0)
172 {
173 fxR3 = mod_fR3*(rx[i]-rxI3[j]);
174 fyR3 = mod_fR3*(ry[i]-ryI3[j]);
175 vx[i]+=fxR3*dt;
176 vx[j]-=fxR3*dt;
177 vy[i]+=fyR3*dt;
178 vy[j]-=fyR3*dt;
179 }
180 }
181
182 if (distR < 2*r)
183 {
184 var mod_fR,fxR,fyR;
185 mod_fR = 12*D*((Math.pow(r/distR),13)-Math.pow((r/distR),7))/a;
186 if (mod_fR >=0)
187 {
188 fxR = mod_fR*(rx[i]-rx[j]);
189 fyR = mod_fR*(ry[i]-ry[j]);
190 vx[i]+=fxR*dt;
191 vx[j]-=fxR*dt;
192 vy[i]+=fyR*dt;
193 vy[j]-=fyR*dt;
194 }
195 }
196
197 Array.prototype.sum = function()
198 {
199 for (var i=0, sum=0; i < this.length; sum += this[i++]);
200 return sum;
201 }
202
203
204 L[i] = (rx[i] * vy[i] - ry[i] * vx[i])/5000;
205 var sumL = L.sum();
206 //график
207 if (steps % 50 == 0) vGraph.graphIter(steps, sumL); // подать данные на график
208 //console.log(sumL);
209
210 }
211
212 }
213 }
214 }
215
216
217 function draw()
218 {
219 ctx.clearRect(0, 0, w_1 , h_1); // очистить экран
220 for (var i = 0; i < N; i++)
221
222 {
223 var x = rx[i];
224 var y = ry[i];
225 ctx.beginPath();
226 ctx.fillStyle = "#00008B";
227 ctx.arc(x, y, r , 0, 2 * Math.PI, false);
228 ctx.closePath();
229 ctx.fill();
230
231 var xI0 = rxI0[i];
232 var yI0 = ryI0[i];
233 ctx.beginPath();
234 ctx.fillStyle = "#00008B";
235 ctx.arc(xI0, yI0, r , 0, 2 * Math.PI, false);
236 ctx.closePath();
237 ctx.fill();
238 var xI1 = rxI1[i];
239 var yI1 = ryI1[i];
240 ctx.beginPath();
241 ctx.fillStyle = "#00008B";
242 ctx.arc(xI1, yI1, r , 0, 2 * Math.PI, false);
243 ctx.closePath();
244 ctx.fill();
245 var xI2 = rxI2[i];
246 var yI2 = ryI2[i];
247 ctx.beginPath();
248 ctx.fillStyle = "#00008B";
249 ctx.arc(xI2, yI2, r , 0, 2 * Math.PI, false);
250 ctx.closePath();
251 ctx.fill();
252 var xI3 = rxI3[i];
253 var yI3 = ryI3[i];
254 ctx.beginPath();
255 ctx.fillStyle = "#00008B";
256 ctx.arc(xI3, yI3, r , 0, 2 * Math.PI, false);
257 ctx.closePath();
258 ctx.fill();
259 var xI4 = rxI4[i];
260 var yI4 = ryI4[i];
261 ctx.beginPath();
262 ctx.fillStyle = "#00008B";
263 ctx.arc(xI4, yI4, r , 0, 2 * Math.PI, false);
264 ctx.closePath();
265 ctx.fill();
266 var xI5 = rxI5[i];
267 var yI5 = ryI5[i];
268 ctx.beginPath();
269 ctx.fillStyle = "#00008B";
270 ctx.arc(xI5, yI5, r , 0, 2 * Math.PI, false);
271 ctx.closePath();
272 ctx.fill();
273 var xI6 = rxI6[i];
274 var yI6 = ryI6[i];
275 ctx.beginPath();
276 ctx.fillStyle = "#00008B";
277 ctx.arc(xI6, yI6, r , 0, 2 * Math.PI, false);
278 ctx.closePath();
279 ctx.fill();
280 var xI7 = rxI7[i];
281 var yI7 = ryI7[i];
282 ctx.beginPath();
283 ctx.fillStyle = "#00008B";
284 ctx.arc(xI7, yI7, r , 0, 2 * Math.PI, false);
285 ctx.closePath();
286 ctx.fill();
287
288 }
289 ctx.beginPath(); // начать рисование
290 ctx.fillStyle="#000000";
291 ctx.moveTo(w, 0); // переместить "карандаш" в точку
292 ctx.lineTo(w, h_1); // нарисовать "карандашом" линию до точки
293 ctx.stroke();
294
295 ctx.beginPath();
296 ctx.fillStyle="#000000";
297 ctx.moveTo(2 * w, 0);
298 ctx.lineTo(2 * w, h_1);
299 ctx.stroke();
300
301 ctx.beginPath();
302 ctx.fillStyle="#000000";
303 ctx.moveTo(0, h);
304 ctx.lineTo(w_1, h);
305 ctx.stroke();
306
307 ctx.beginPath();
308 ctx.fillStyle="#000000";
309 ctx.moveTo(0, 2 * h);
310 ctx.lineTo(w_1, 2 * h);
311 ctx.stroke();
312 }
313
314 setInterval(step, 1000/120); // функция step будет запускаться 60 раз в секунду (60 раз / 1000 мс)
315
316 }
