ChainSimulation/main.cpp at master · RichtersFinger/ChainSimulation · GitHub

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
/* g++ ChainElement.cpp ChainLink.cpp Chain.cpp Simulation.cpp main.cpp -std=c++11 -o run -O3
 */

#include <vector>
#include <iostream>
#include <sstream>
#include <fstream>
#include <stdio.h>
#include <iomanip> // precision //
#include <string.h>
#include <math.h>
#include <cstdlib> // random numbers //
#include <time.h> // for 'random' seed //

#include "Vector.h"
#include "ChainLink.h"
#include "ChainElement.h"
#include "Wall.h"
#include "Chain.h"
#include "Simulation.h"

int main(int argc, char *argv[]){
	Simulation thisSimulation;

	// setup container
	double containerposx = 0.0;
	double containerposy = 0.0;
	double containerheight = 1.0;
	double containerwallwidth = 0.1;
	double containerwidth = 0.6;
	double dropheight = 25.0;
	thisSimulation.addTrailofWalls(std::vector<Vector>{Vector(containerposx, containerposy + containerheight), Vector(containerposx, containerposy), // left wall
								Vector(containerposx, containerposy), Vector(containerposx + containerwidth, containerposy), // floor
								Vector(containerposx + containerwidth, containerposy), Vector(containerposx + containerwidth, containerposy + containerheight), // right wall
								Vector(containerposx + containerwidth + containerwallwidth, containerposy + containerheight), Vector(containerposx + containerwidth + containerwallwidth, containerposy - containerwallwidth), // right wall outer

								Vector(-8.0, containerposy + containerheight - dropheight), Vector(8.0, containerposy + containerheight - dropheight)});
	// change this to account for fast impact of chain
	thisSimulation.walls.back()->thickness = 1.0;
	// rounded off tips for container walls
	thisSimulation.addSphericalBoundary(containerposx + containerwidth + 0.5*containerwallwidth, containerposy + containerheight, 0.5*containerwallwidth);
	thisSimulation.addSphericalBoundary(containerposx - 0.5*containerwallwidth, containerposy + containerheight, 0.5*containerwallwidth);
	thisSimulation.writeWallsToFile("walls.dat");

	// make chain
	thisSimulation.initChain(16*200, 0.04, 1.0e6, 0.5, 1.0e6);

	// arrange chain elements in a loop above the container
	thisSimulation.chain.quickArrangementEntwined(Vector(containerposx + 0.5 * containerwidth, containerposy + 0.5 * containerheight), 0.33 * containerwidth, 16*8, 0.1);

	// load file or integrate drop
	if (false) {
		thisSimulation.loadChainBackup("chainbackup.dat");
	} else {
		for (int i = 0; i < 150; i++) {
			thisSimulation.writeChainToFile(0);
			thisSimulation.step(floor(1.0/60.0/thisSimulation.dt));
			std::cout << i << " " << thisSimulation.time << std::endl;
		}
		thisSimulation.saveChainBackup("chainbackup.dat");
	}
	// pull first element out of container
	thisSimulation.chain.elements[0]->externalax = 300.0;
	thisSimulation.chain.elements[0]->externalay = 1000.0;
	std::cout << "start" << std::endl;
	for (int i = 0; i < 1200; i++) {
		std::cout << i << " " << thisSimulation.time << std::endl;
		// stop pulling at some point
		if (thisSimulation.chain.elements[0]->y > containerposy + containerheight + 2.0 && thisSimulation.chain.elements[0]->x > containerposx + containerwidth + containerwallwidth + 0.5) {
			thisSimulation.chain.elements[0]->externalax = 0.0;
			thisSimulation.chain.elements[0]->externalay = 0.0;
		}
		thisSimulation.writeChainToFile(0);
		thisSimulation.writeChainToFile(i+1);
		thisSimulation.step(floor(1.0/60.0/thisSimulation.dt));
	}

	return 0;
}