There are two bodies colliding Simulate the collision phenomena by num

Click For Summary
SUMMARY

This discussion focuses on simulating the collision phenomena of two bodies using numerical integration techniques, specifically through the Runge-Kutta method in MATLAB. The primary objectives are to confirm the conservation of total kinetic energy in the absence of damping forces and to demonstrate the reduction of kinetic energy during collisions when damping is present. The governing equations include mx''(t) + kx(t) = 0 and x''(t) + (w_n)^2 x(t) = 0, with parameters such as mass (m = 100), spring constant (k = 20), and damping coefficient (c = 10) being utilized in the simulation.

PREREQUISITES
  • Understanding of numerical integration techniques, specifically Runge-Kutta methods.
  • Familiarity with MATLAB programming for implementing simulations.
  • Knowledge of classical mechanics, particularly conservation of kinetic energy principles.
  • Basic understanding of differential equations and their applications in physics.
NEXT STEPS
  • Explore advanced MATLAB functions for simulating dynamic systems.
  • Research the implications of damping forces in collision simulations.
  • Learn about the mathematical derivation of conservation laws in mechanics.
  • Investigate alternative numerical methods for solving differential equations, such as Euler's method or Verlet integration.
USEFUL FOR

Students in physics or engineering disciplines, particularly those studying dynamics and numerical methods, as well as educators looking for practical examples of collision simulations in MATLAB.

max1546
Messages
7
Reaction score
0

Homework Statement


There are two bodies colliding
Simulate the collision phenomena by numerical integration.
(1) Confirm that the total kinetic energy is conserved if you ignore the damping forces.
(2) Show that the total kinetic energy is reduced through collision if there exists any damping.
http://img208.imageshack.us/img208/6222/springzj.jpg

Homework Equations


mx''(t)+kx(t) = 0
x''(t) + (w_n)^2 x(t) = 0
x(t) = A_1 e^(l w_n t) + A_2 e^(-l w_n t)


The Attempt at a Solution


I assumed that the second mass does not have any spring and infinite mass(wall)
Then used Runge Kutta through Matlab
===========
function Y=Func(X)
% X=(x,y), Y=(Dx,Dy)

m = 100;
k = 20;
c = 10;

x=X(1);
y=X(2);


Dx = y;
Dy = (- k* x - c* y)/m ;

Y=[Dx, Dy];

end
========
I am not sure how to arrange differential equations for this problem case in general.
 
Last edited by a moderator:
Physics news on Phys.org
Also, I am not sure how to prove the conservation of kinetic energy and the reduction of it.Any help will be appreciated.
 

Similar threads

  • · Replies 1 ·
Replies
1
Views
2K
  • · Replies 1 ·
Replies
1
Views
2K
  • · Replies 5 ·
Replies
5
Views
4K
  • · Replies 20 ·
Replies
20
Views
3K
  • · Replies 1 ·
Replies
1
Views
6K
Replies
21
Views
4K
Replies
2
Views
4K
  • · Replies 2 ·
Replies
2
Views
2K
Replies
4
Views
9K
Replies
7
Views
2K