Движение частицы в полях(JavaScript)

Материал из Department of Theoretical and Applied Mechanics
Версия от 13:09, 2 июня 2020; Kosoegor (обсуждение | вклад) (Новая страница: «==Описание== Движение частицы в гравитационном и электромагнитном полях Исполнители: Г…»)

(разн.) ← Предыдущая | Текущая версия (разн.) | Следующая → (разн.)
Перейти к: навигация, поиск

Описание

Движение частицы в гравитационном и электромагнитном полях

Исполнители: Галов Владислав, Шилов Михаил, Борисенков Богдан

Визуализация

Код программы

Код программы на языке JavaScript:
  1 window.addEventListener("load", main, false);
  2 
  3 function main(){
  4 	var ctx = canv.getContext('2d');
  5 	var h = canv.height;
  6 	var w = canv.width;
  7 	var E = 10; var Mod = 100000000; var Ma1 = 100000000; var Ma2 = 50000000;
  8 	var radod = 100;
  9 	var rad1 = 100;
 10         var rad2 = 50;
 11 	var rad = 10;
 12 	var x; var y; var m;
 13 	var vx = 0; var vy = 0; var aex = 0; var aey = 0; var atsx = 0; var atsy = 0;
 14 	var axg = 0; var ayg = 0;
 15 	var Flx = 0; var Fly = 0; var xo; var yo;
 16 	var fps = 60;
 17 	var dt = 1/fps;
 18 	
 19 	canv.onclick = function(event){
 20 		dt = 1/fps;
 21 		var rect = canv.getBoundingClientRect();
 22 		x = event.clientX - rect.left;
 23 		y = event.clientY - rect.top;
 24 }
 25 
 26 	
 27 function phys1(){
 28 	m = document.getElementById('polz').value;
 29 	var radf = document.getElementsByName('field');
 30 	var radn = document.getElementsByName('kolvo');
 31 	var rads = document.getElementsByName('sput')
 32 	var radd = document.getElementsByName('direct');
 33 	var radc = document.getElementsByName('ch');
 34 	if (radf[0].checked==true){
 35 	if (rads[1].checked==true){
 36 		if (radn[0].checked==true){
 37 			var rod = Math.sqrt((x-250)*(x-250)+(y-250)*(y-250));
 38 			var cosod = (250-x)/rod;
 39 			var sinod = (250-y)/rod;
 40 			var fg = m*Mod/(rod*rod);
 41 			var fgx = fg*cosod; var fgy = fg*sinod;
 42 			axg = fgx/m; ayg = fgy/m;
 43 			x = x+3*axg*dt*dt/2;
 44 			y = y+3*ayg*dt*dt/2;
 45 		}
 46 		if (radn[1].checked==true){
 47 			var r1 = Math.sqrt((x-350)*(x-350)+(y-350)*(y-350));
 48 			var r2 = Math.sqrt((x-100)*(x-100)+(y-100)*(y-100));
 49 			var cos1 = (350-x)/r1;
 50 			var cos2 = (100-x)/r2;
 51 			var sin1 = (350-y)/r1;
 52 			var sin2 = (100-y)/r2;
 53 			var fg1 = m*Ma1/(r1*r1);
 54 			var fg2 = m*Ma2/(r2*r2);
 55 			var fg1x = fg1*cos1; var fg2x = fg2*cos2;
 56 			var fg1y = fg1*sin1; var fg2y = fg2*sin2;
 57 			axg = (fg1x+fg2x)/m; ayg = (fg1y+fg2y)/m;
 58 			x = x+3*axg*dt*dt/2;
 59 			y = y+3*ayg*dt*dt/2;
 60 			if (r1<rad1+rad)dt=0;
 61 			if (r2<rad2+rad)dt=0;
 62 			}
 63 	}
 64 	if (rads[0].checked==true){
 65 		if (radn[0].checked==true){
 66 			var rsod = Math.sqrt((x-250)*(x-250)+(y-250)*(y-250));
 67 			var cossod = (250-x)/rsod;
 68 			var sinsod = (250-y)/rsod;
 69 			var v = Math.sqrt(0.1*Mod/rsod);
 70 			vx = v*sinsod; vy = v*(-cossod);
 71 			ac = v*v/rsod;
 72 			acx = ac*cossod; acy = ac*sinsod;
 73 			x = x+vx*dt+acx*dt*dt/2;
 74 			y = y+vy*dt+acy*dt*dt/2;
 75 		}
 76 		if (radn[1].checked==true){
 77 			var rs1 = Math.sqrt((x-350)*(x-350)+(y-350)*(y-350));
 78 			var rs2 = Math.sqrt((x-100)*(x-100)+(y-100)*(y-100));
 79 			if (rs1<=rs2+70){
 80 				var coss1 = (350-x)/rs1;
 81 				var sins1 = (350-y)/rs1;
 82 				var v = Math.sqrt(0.1*Ma1/rs1);
 83 				vx = v*sins1; vy = v*(-coss1);
 84 				ac = v*v/rs1;
 85 				acx = ac*coss1; acy = ac*sins1;
 86 				x = x+vx*dt+acx*dt*dt/2;
 87 				y = y+vy*dt+acy*dt*dt/2;
 88 			}
 89 			else{
 90 				var coss2 = (100-x)/rs2;
 91 				var sins2 = (100-y)/rs2;
 92 				var v = Math.sqrt(0.1*Ma2/rs2);
 93 				vx = -v*sins2; vy = v*coss2;
 94 				ac = v*v/rs2;
 95 				acx = ac*coss2; acy = ac*sins2;
 96 				x = x+vx*dt+acx*dt*dt/2;
 97 				y = y+vy*dt+acy*dt*dt/2;
 98 			}
 99 		}
100 	}
101 	}
102 	if (radf[1].checked==true){
103 		if (radd[0].checked==true && radc[0].checked==true ||
104 		radd[1].checked==true && radc[1].checked==true){
105 			atsy = vx/m;
106 			atsx = (-vy)/m;
107 		}
108 		else{
109 			atsy = (-vx)/m;
110 			atsx = vy/m;
111 		}
112 		var radd2 = document.getElementsByName('direct2');
113 		if (radd2[0].checked==true){aex = E/m;}
114 		if (radd2[1].checked==true){aex = -E/m;}
115 		x = x+vx*dt+(aex+atsx)*dt*dt/2;
116 		y = y+vy*dt+atsy*dt*dt/2;
117 		vx = vx+(aex+atsx)*dt;
118 		vy = vy+atsy*dt;
119 		if (x<-10 || y<-10 || x>w+10 || y>h+10) {vx = 0; vy = 0; atsx = 0; atsy = 0;};
120 	}
121 }
122 
123 	
124 function draw1(){
125 	ctx.clearRect(0,0,500,500);
126 	ctx.lineWidth="1";
127 	var radf = document.getElementsByName('field');
128 	var radn = document.getElementsByName('kolvo');
129 	var radd = document.getElementsByName('direct');
130 	var radc = document.getElementsByName('ch');
131 	if (radf[0].checked==true){
132 		if (radn[0].checked==true){
133 			ctx.beginPath();
134 			ctx.arc(250,250,radod,0,2*Math.PI);
135 			ctx.stroke();
136 			ctx.fillStyle="grey";
137 			ctx.fill();
138 			var rod = Math.sqrt((x-250)*(x-250)+(y-250)*(y-250));
139 			if (rod<=radod+rad){
140 				var sinod = Math.abs(y-250)/rod;
141 				var cosod = Math.abs(x-250)/rod;
142 				if (x>250){x = 250+(radod+rad)*cosod;}
143 				else{x = 250-(radod+rad)*cosod;}
144 				if (y>250){y = 250+(radod+rad)*sinod;}
145 				else{y = 250-(radod+rad)*sinod;}
146 			}
147 			ctx.beginPath();
148 			ctx.arc(x,y,rad,0,2*Math.PI);
149 			ctx.stroke(); ctx.fillStyle="black";
150 			ctx.fill();
151 		}
152 		if (radn[1].checked==true){
153 			ctx.beginPath();
154 			ctx.arc(350,350,rad1,0,2*Math.PI);
155 			ctx.stroke();
156 			ctx.fillStyle="brown";
157 			ctx.fill();
158 			ctx.beginPath();
159 			ctx.arc(100,100,rad2,0,2*Math.PI);
160 			ctx.stroke();
161 			ctx.fillStyle="purple";
162 			ctx.fill();
163 			var r1 = Math.sqrt((x-350)*(x-350)+(y-350)*(y-350));
164 			var r2 = Math.sqrt((x-100)*(x-100)+(y-100)*(y-100));
165 			if (r1<=rad1+rad){
166 				var sin1 = Math.abs(y-350)/r1;
167 				var cos1 = Math.abs(x-350)/r1;
168 				if (x>350){x = 350+(rad1+rad)*cos1;}
169 				else{x = 350-(rad1+rad)*cos1;}
170 				if (y>350){y = 350+(rad1+rad)*sin1;}
171 				else{y = 350-(rad1+rad)*sin1;}
172 			}
173 			if (r2<=rad2+rad){
174 				var sin2 = Math.abs(y-100)/r2;
175 				var cos2 = Math.abs(x-100)/r2;
176 				if (x>100){x = 100+(rad2+rad)*cos2;}
177 				else{x = 100-(rad2+rad)*cos2;}
178 				if (y>100){y = 100+(rad2+rad)*sin2;}
179 				else{y = 100-(rad2+rad)*sin2;}
180 			}
181 			ctx.beginPath();
182 			ctx.arc(x,y,rad,0,2*Math.PI);
183 			ctx.stroke(); ctx.fillStyle="black";
184 			ctx.fill();
185 		}
186 	}
187 	if (radf[1].checked==true){
188 		ctx.beginPath();
189 		ctx.arc(x,y,rad,0,2*Math.PI);
190 		ctx.stroke(); 
191                 ctx.fillStyle="black";
192 		ctx.fill();
193 		ctx.beginPath();
194 		ctx.arc(30,30,15,0,2*Math.PI);
195 		ctx.strokeStyle="grey";
196 		ctx.stroke();
197 		if (radd[1].checked==true){
198 			ctx.beginPath();
199 			ctx.moveTo(19.4,19.4);
200 			ctx.lineTo(40.6,40.6);
201 			ctx.stroke();
202 			ctx.beginPath();
203 			ctx.moveTo(19.4,40.6);
204 			ctx.lineTo(40.6,19.4);
205 			ctx.stroke();
206 		}
207 		else{
208 			ctx.beginPath();
209 			ctx.arc(30,30,1,0,2*Math.PI);
210 			ctx.fillStyle="grey";
211 			ctx.fill();
212 			ctx.stroke();
213 		}
214 		ctx.beginPath();
215 		ctx.moveTo(x-5,y);
216 		ctx.lineTo(x+5,y);
217 		if (radc[0].checked==true){
218 			ctx.moveTo(x,y-5);
219 			ctx.lineTo(x,y+5);
220 		}
221 		ctx.strokeStyle="white";
222 		ctx.lineWidth="3";
223 		ctx.stroke();
224 	}
225 }
226 
227 	
228 function control(){
229 	phys1();
230 	draw1();
231 }
232 
233 	
234 var radf = document.getElementsByName('field');
235 if (radf[0].checked==true){setInterval(control, 100/fps);}
236 if (radf[1].checked==true){setInterval(control, 0.1/fps);}
237 }