- #1
polomeister
- 5
- 0
Hi All,
This is part solution, part request for peer review.
I wasn't able to find a model for the motion of a particle dropped onto an elastic band so I've created my own in Excel (attached).
Note that the model is primarily focussed on the motion of the particle, and not the elastic band, although in reality the two are dependent.
Key Results
Requests from the Forum
The attached file
Algorithm
There is a complication to the equations of motion in that, once the particle has hit the elastic band, the force exerted on the particle is a function of position.
Assumptions/Limitations
Improvements
Thanks & best wishes – I hope this is a worthy edition to the forum and perhaps it may even help someone out in future.
P
This is part solution, part request for peer review.
I wasn't able to find a model for the motion of a particle dropped onto an elastic band so I've created my own in Excel (attached).
Note that the model is primarily focussed on the motion of the particle, and not the elastic band, although in reality the two are dependent.
Key Results
- Once the particle has hit the band the vertical motion is like a sine wave.
- If the particle is dropped off-centre there is also horizontal motion.
- Movement towards the middle of the band might be expected – as this should minimise the potential energy.
- However, the particle seems to move briefly away from the centre before returning. (I found this interesting and counter-intuitive.)
- In general, the particle traverses a large proportion of the elastic band.
- For certain combinations of displacement, mass and elasticity, the particle bounces off the elastic band.
Requests from the Forum
- It would be great if someone could check over my working. This is one component of a larger simulation I'm working on.
- Tips for reducing the approximation errors in the working (see "Assumptions/Limitations" below).
- Experimental verification that a particle does initially move away from the centre. Might be a nice project for a physics student!
The attached file
- The format is the new Excel binary format .xlsb (even 1 row of simulation exceeds the 484k limit if saving as an .xlsx file). Send me a PM if you require the file in a different format.
- There are four tabs
- "Diagram" – Shows the setup when the particle has displaced the elastic band by a vertical distance
. Forces induced in the two components of the elastic band are assumed to act on the particle. - "Particle Dropped on Band" – This is the main working.
- Boxes highlighted in green are inputs
- I've had to curtail the number of rows to allow upload. I recommend at least 10,000.
- The time-step is 0.001s. This gives a relatively high error-rate (about 10% as measured by Energy). Reducing the time-step to 0.0005s or less gives much improved accuracy. (I was suprised that the algorithm required such small increments to calculate accurately.)
- See "Algorithm" below for basis of calculation
- There is a block on the right hand side that calculates the total Energy in the system, which I've used as a check.
- The last row will always have an error in it and should be excluded from summary calculations.
- "Graph of x, y and energy" – Graph showing
, and Energy with respect to time. Note: It won't look right until you complete enough rows. - "Forces on Elastic Band" – This was a bit of pre-working that I've left in. It calculates the forces in the elastic band just by pulling it down a vertical distance
at a horizontal displacement . Note is positive whereas in the main working it is negative.
- "Diagram" – Shows the setup when the particle has displaced the elastic band by a vertical distance
Algorithm
There is a complication to the equations of motion in that, once the particle has hit the elastic band, the force exerted on the particle is a function of position.
- Position, velocity and force are calculated component-wise. The formulae below apply for both components.
- Position:
- Change in position is an average of the starting and ending velocities. This improves the accuracy of the modelling.
- Velocity:
- Change in velocity is based on the force at the start of the time-step. Because force depends on position, it's not possible to take an average without generating circular calculations.
- Force:
- Force is calculated using the starting position of the particle for each time-step.
- The forces from the elastic band only operate in the region
. This is handled through the indicator function in Column U. - The elastic extension forces are
and , where-
and are the hypotenueses generated by extending the elastic band vertically down at horizontal position -
is the length of the elastic band -
is the elasticity coefficient
-
Assumptions/Limitations
- Ends of band are fixed
- Left peg of the band is the Origin
- Band is weightless
- Only two dimensions
- The particle is updated before the elastic band
- The elastic band is assumed to comprise of a left hand and right hand component with a vertex at the position of the particle.
- No friction
- Band is assumed to reset to a horizontal line once particle is higher than
– in reality it would continue to oscillate. This probably means the model is only valid for the period up to the particle achieving its maximum height for the first time. One might expect the band to oscillate in height as much as the maximum depth - Looks like a fixed amount of energy is being added per cycle
- I've only checked motion when the particle is dropped vertically
- Note: You can waste a lot of iterations in the 'freefall' phase of the model. A quick workaround is to work out the velocity the particle would have at
and start from there.
Improvements
- Sort out the Energy issue (but reducing the size of the time-steps mitigates this)
- Write a macro to animate the movement of the particle
- Improve the modelling of the elastic band:
- First step would be to treat the components as two elastic rods, with a mass per unit length
- Next step would be to further segment the elastic band into
such rods
- Please suggest more!
Thanks & best wishes – I hope this is a worthy edition to the forum and perhaps it may even help someone out in future.
P