Редактирование: Мещерский 48.2

Мещерский Задача 48.2

Визуализация 3D-задачи по динамике на JavaScript

Исполнитель: [[Сызранцев Денис|Сызранцев Денис]]

Группа 23632/2 Кафедра Теоретической механики

==Условие задачи==

Барабан центрифуги приводится во вращение электродвигателем через двухступенчатый редуктор. Заданы момент инерции J0 электродвигателя, момент инерции J2 барабана, момент инерции J1 промежуточного вала редуктора, передаточные числа i01 и i12 ступеней редуктора. К ротору электродвигателя приложен вращающий момент M0 и момент сил сопротивления M0s, к валу редуктора и к барабану — моменты сил сопротивления M1s и M2s соответственно. Составить дифференциальное уравнение вращения барабана центрифуги.

==Решение задачи==

<math> d(δT/δw)/dt - δT/δφ = 0</math><br/>
<math>T_0 = -J_0 w_0^2/2</math><br/>
<math>T_1 = -J_1 w_1^2/2</math><br/>
<math>T_2 = -J_2 w_2^2/2</math><br/>
<math>i_(01) = w_0/w_1</math><br/>
<math>i_(12) = w_1/w_2</math><br/>
<math>d(δT/δw)/dt = J_0 i_01^2 i_(12)^2 + J_0 i_(12)^2 + J_2) b_0</math><br/>
<math>(J_0 i_(01)^2 i_(12)^2 + J_0 i_(12)^2 + J_2) b_0 = (M - M_(0s)) i_(01) i_(12) - M_(1s) i_(12) - M_(2s)</math><br/>

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

{{#widget:Iframe |url=http://tm.spbstu.ru/htmlets/Syzrantsev_DA/index.html |width=850 |height=400 |border=0 }}

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

<syntaxhighlight lang="javascript" line start="1" enclose="div">
<!DOCTYPE html>
<html lang="en">
	<head>
		<title>48.2</title>
		<meta charset="utf-8">
		<meta name="viewport" content="width=device-width, user-scalable=no, minimum-scale=1.0, maximum-scale=1.0">
		<style>
			body {
				font-family: Monospace;
				background-color: #f0f0f0;
				margin: 0px;
				overflow: hidden;
			}
		</style>
	</head>
	<body>

		<script src = "http://tm.spbstu.ru/htmlets/Belousova/Statika/three.js">
		</script>
		<script src ="http://tm.spbstu.ru/htmlets/Belousova/Statika/stats.min.js">
		</script>
		<script src ="http://tm.spbstu.ru/htmlets/Belousova/Statika/OrbitControls.js">
		</script>
		<script src ="http://tm.spbstu.ru/htmlets/Belousova/Statika/dat.gui.js">
		</script>

		<script>

			var container, stats;

			var camera, controls, scene, renderer;

			var start = Date.now();
			
			var A2 , w0, w1, w2;

			
			var c = new function() {
				this.i01 = 2;
				this.i12 = 3;
				this.J0 = 80;
				this.J1 = 40;
				this.J2 = 80;
				this.M0 = 3;
				this.M0s = 2;
				this.M1s = 1.6;
				this.M2s = 1;
				this.stop = false;
				this.restart = function(){
				start = Date.now()
				ED.rotation.x = 0;
				rot.rotation.x =  0;
				B.rotation.x = 0;
				}
				
				this.w2 = 0.0;
				this.A2 = ((this.M0-this.M0s)*this.i01*this.i12-this.M1s*this.i12-this.M2s)/(this.J0*(this.i01*this.i12)^2+this.J1*this.i12^2+this.J2a)
			}
			
			
			var width = window.innerWidth;
			var height = window.innerHeight;
			var le1 = 7;
			var le2 = 3;
			var re1 = 2;
			var re2 = 0.3;
			var lb1 = 9;
			var rb1 = 3;
			var rr1 = re1/c.i01;
			var rb3 = rb1/c.i12;
			var rr3 = rb1;
			
			var texture = new THREE.TextureLoader().load( "texture.jpg" );
				texture.wrapS = THREE.RepeatWrapping;
				texture.wrapT = THREE.RepeatWrapping;
				texture.repeat.set( 3, 3 );
			var material = new THREE.MeshLambertMaterial({color: 0xFFC673,map:texture});

			function physics(t,c) {
			
				rr1 = re1/c.i01;
				rb3 = rb1/c.i12;
				
				ED.position.y = -re1-rr1;
				B.position.y = rb1+rr3;
				rot.remove(rot1);
				rot1 = new THREE.Mesh( new THREE.CylinderGeometry( rr1, rr1, 0.01, 64 ), material );
				rot1.rotation.z = - Math.PI / 2;
				rot1.position.x = -le2/2-5;
				rot1.castShadow = true;	
				rot.add( rot1 );
				B.remove(B3);
				B3 = new THREE.Mesh( new THREE.CylinderGeometry( rb3, rb3, 0.01, 64 ), material );
				B3.rotation.z = - Math.PI / 2;
				B3.position.x = le2/2-5;
				B3.castShadow = true;
				B.add( B3 );
				B.position.y = rb3+rr3;
				
				if (c.stop == false){
					var w1,w0;
					c.A2 = ((c.M0-c.M0s)*c.i01*c.i12-c.M1s*c.i12-c.M2s)/(c.J0*(c.i01*c.i12)^2+c.J1*c.i12^2+c.J2a);
					if ((c.A2<0)&&(c.w2 <= 0)){
					c.A2 = 0;
					c.w2 = 0;
					} else {
					c.w2 = c.A2*t;
					w1 = c.w2/c.i12;
					w0 = w1/c.i01;
					
				
					ED.rotation.x = t * w0*Math.pow(10, -3);
					rot.rotation.x =  - t * w1*Math.pow(10, -3);
					B.rotation.x = t * c.w2*Math.pow(10, -3);
					}
				}
				
				

			}
			
			function init() {
				
				
				container = document.createElement( 'div' );
				document.body.appendChild( container );
				
				scene = new THREE.Scene();
				scene.background = new THREE.Color( 0xbfd1e5 );

				camera = new THREE.PerspectiveCamera( 45, width / height, 0.1, 1000 );
				camera.position.x = -20;
				camera.position.y = 5;
				camera.position.z = 20;
				camera.lookAt(0,0,0)

				

				

				
				
				var ambient = new THREE.AmbientLight( 0xffffff, 0.3 );
				scene.add( ambient );
				
				spotLight = new THREE.SpotLight( 0xffffff, 2 );
				spotLight.position.set( -30, 45, 35 );
				spotLight.angle = Math.PI / 3;
				spotLight.penumbra = 0.05;
				spotLight.decay = 2;
				spotLight.distance = 200;
				spotLight.castShadow = true;
				spotLight.shadow.mapSize.width = 1024;
				spotLight.shadow.mapSize.height = 1024;
				spotLight.shadow.camera.near = 10;
				spotLight.shadow.camera.far = 200;
				scene.add( spotLight );
				
				
				ED = new THREE.Object3D;
				ED1 = new THREE.Mesh( new THREE.CylinderGeometry( re1, re1, le1, 64 ), material );
				ED1.rotation.z = - Math.PI / 2;
				ED1.position.x = -le1/2-3*le2/2-5;
				ED1.castShadow = true;
				ED.add(ED1);
				ED2 = new THREE.Mesh( new THREE.CylinderGeometry( re2, re2, le2, 64 ), material );
				ED2.rotation.z = - Math.PI / 2;
				ED2.position.x = -le2-5;
				ED2.castShadow = true;
				ED.add(ED2);
				ED3 = new THREE.Mesh( new THREE.CylinderGeometry( re1, re1, 0.01, 64 ), material );
				ED3.rotation.z = - Math.PI / 2;
				ED3.position.x = -le2/2-5;
				ED3.castShadow = true;
				ED.add( ED3 );
				ED.position.y = -re1-rr1;
				scene.add( ED );
				
				rot = new THREE.Object3D;
				rot1 = new THREE.Mesh( new THREE.CylinderGeometry( rr1, rr1, 0.01, 64 ), material );
				rot1.rotation.z = - Math.PI / 2;
				rot1.position.x = -le2/2-5;
				rot1.castShadow = true;				
				rot.add( rot1 );
				rot2 = new THREE.Mesh( new THREE.CylinderGeometry( re2, re2, le2, 64 ), material );
				rot2.rotation.z = - Math.PI / 2;
				rot2.position.x = -5;
				rot2.castShadow = true;	
				rot.add( rot2 );
				rot3 = new THREE.Mesh( new THREE.CylinderGeometry( rr3, rr3, 0.01, 64 ), material );
				rot3.rotation.z = - Math.PI / 2;
				rot3.position.x = le2/2-5;
				rot3.castShadow = true;	
				rot.add( rot3 );
				scene.add( rot );
				
				B = new THREE.Object3D;
				B1 = new THREE.Mesh( new THREE.CylinderGeometry( rb1, rb1, lb1, 64 ), material );
				B1.rotation.z = - Math.PI / 2;
				B1.position.x = lb1/2+3*le2/2-5;
				B1.castShadow = true;	
				B.add(B1);
				B2 = new THREE.Mesh( new THREE.CylinderGeometry( re2, re2, le2, 64 ), material );
				B2.rotation.z = - Math.PI / 2;
				B2.position.x = le2-5;
				B2.castShadow = true;
				B.add(B2);
				B3 = new THREE.Mesh( new THREE.CylinderGeometry( rb3, rb3, 0.01, 64 ), material );
				B3.rotation.z = - Math.PI / 2;
				B3.position.x = le2/2-5;
				B3.castShadow = true;
				B.add( B3 );
				B.position.y = rb3+rr3;
				scene.add( B );

				plane = new THREE.Mesh( new THREE.BoxGeometry( 2000, 1, 2000 ), new THREE.MeshPhongMaterial( { color: 0x679ED2, dithering: true } ) );
				plane.position.y = - 15;
				plane.receiveShadow = true;
				scene.add( plane );
				
				renderer = new THREE.WebGLRenderer();
				renderer.setSize( width, height );
				renderer.shadowMap.enabled = true;
				renderer.shadowMap.type = THREE.PCFSoftShadowMap;
				container.appendChild( renderer.domElement );
				
				controls = new THREE.OrbitControls( camera,renderer.domElement );
				
				var gui = new dat.GUI();
					gui.add(c, 'i01',1,5);
					gui.add(c, 'i12',1,5);
					gui.add(c, 'J0',1,1000);
					gui.add(c, 'J1',1,1000);
					gui.add(c, 'J2',1,1000);
					gui.add(c, 'M0',0,10);
					gui.add(c, 'M0s',0,10);
					gui.add(c, 'M1s',0,10);
					gui.add(c, 'M2s',0,10);
					gui.add(c, 'stop');
					gui.add(c, 'restart');
					gui.add(c, 'w2').listen();
					gui.add(c, 'A2').listen();
			}

			function onWindowResize() {

				camera.aspect = window.innerWidth / window.innerHeight;
				camera.updateProjectionMatrix();

				renderer.setSize( window.innerWidth, window.innerHeight );
				

			}

			//

			function animate() {
				
					requestAnimationFrame( animate );
					render();
			}

			function render() {

				var t = Date.now() - start;
				
				physics(t,c)
				
				ED.update;
				B.update;
				rot.update;
					
				controls.update();

				renderer.render( scene, camera );

			}
			
			init();
			animate();
			

		</script>
	</body>
</html>
</syntaxhighlight>