Текущая версия |
Ваш текст |
Строка 8: |
Строка 8: |
| ==Who we are== | | ==Who we are== |
| * [[Степанов Алексей]] (responsible for contact model) | | * [[Степанов Алексей]] (responsible for contact model) |
− | * [[Дзенушко Дайнис]] (responsible for factory and EDEM simulation) | + | * [[Дзенушко Дайнис]] (responsible for factory and EDEM simylation) |
− | | |
| ==Factory== | | ==Factory== |
| We create 2 particles on a distance 0.1m and rotated at an angle of 5-15 () degrees;<br> | | We create 2 particles on a distance 0.1m and rotated at an angle of 5-15 () degrees;<br> |
Строка 24: |
Строка 23: |
| | | |
| ==Contact Model== | | ==Contact Model== |
− | In our project we used simple contact model.<br>
| + | Alexey will soon write about it... |
− | All central forces we made zero. Also we made zero all torques except torque, which operates on a "small" particle. <br>
| |
− | calculatedElem2AdditionalTorqueX = 1.0 * (elem1Orientation[7] * elem2Orientation[4] + elem2Orientation[5] * elem1Orientation[4] ); <br>
| |
− | In this formula С = 1.0
| |
− | | |
| ==EDEM simulation== | | ==EDEM simulation== |
− | ===Globals:===
| |
− | Interaction: Particle to particle<br>
| |
− | Model: our contact model<br>
| |
− | No gravity<br>
| |
− | There are two materials "material" and "material_2" with different density for "material" 1000 for "material_2" 1.7e+05 <br>
| |
− | Restitution: 0.5<br>
| |
− | No static and rolling friction<br>
| |
− | ===Particles:===
| |
− | We create particles of 2 types;"big" with big moment of inertia (100kgm2 X-axis) and "small"(0.000285kgm2 X-axis);Both particles are made of 2 surfaces placed along Z-axis on a distance of 2 particle radius
| |
− | ==Measures==
| |
− | We measured the period of oscillation<br>
| |
− | ===Analytics===
| |
− | <math>T = 2\pi\sqrt{\frac{C}{\theta}} = 0.106</math>
| |
− | ===Integration===
| |
− | We measured period using the Graph of angular velocity and got the result <br>
| |
− | <math>T = 0.082</math>
| |
− | ==Results==
| |
− | So we got the result that is different to the analytic;<br>
| |
− | Then we wrote small programm on C# that integrates the equalation in a simple way and got the result<br>
| |
− | <math>T = 0.106</math><br>
| |
− | We checked all in EDEM simulation; Changing parameters of particles, roll stiffness didn't gave any difference the results stayed different to analytic.
| |
− | ==Additional Measures==
| |
− | We chacked translational kinetic and potentional energy and they really equals to zero<br>
| |
− | Graph of Rotational Kinetic energy
| |
− | <gallery widths=580px heights=350px perrow = 1>
| |
− | Файл:EDEM_Kin_Energy.jpg
| |
− | </gallery>
| |
− | As we can see the maximum of Kinetic energy is conserving. This graph was extended to 100 seconds(~1220 period)<br>
| |
− | and conservation of Rotational Kinetic energy was prooved. Time step ~9.75e-07s.
| |
− |
| |
− |
| |
− | [[Category: Студенческие проекты]]
| |
− | [[Category: Механика дискретных сред]]
| |
− | [[Category: Программирование]]
| |