Ball and Spring Model: Difference between revisions

From Physics Book
Jump to navigation Jump to search
 
(22 intermediate revisions by 3 users not shown)
Line 1: Line 1:
Page created by Ashley Fleck.  
Page created by Ashley Fleck.  


The Interactions of atoms can be modeled using balls to represent the atoms and a spring to represent the chemical bond between them.  
The interactions of atoms can be modeled using balls to represent the atoms and a spring to represent the chemical bond between them.  


==Introduction==
scene.title = "Ball and Spring Model of Matter\n"
scene.width = 900
scene.height = 600
scene.background = color.white


:Everything that has mass (solids, liquids, gases, etc.) is made of matter, which is made up of atoms. Each atom is composed of a dense positively charged nucleus and a much less dense negatively charged electron cloud (1). [Insert picture of a breakdown of an atom]. The charges of the atom affect how they interact with other atoms. Simply, the positive charge of one positively charged nucleus of one atom is attracted to the negatively charged electron cloud of the other atom, which bounds the two atoms together. This interaction is known as the electric force. However, this attraction only occurs up to a specific point. Once the atoms near too close to one another, the attractive forces of the nucleus and electron cloud are outweighed by the repulsive forces of the two atoms’ electron clouds (while opposite charges attract, negative charges repel). Therefore, if the atoms are pushed too close to one another, they will repel each other rather than attract. (1).
scene.append_to_caption("""
<h1>Ball and Spring Model of Matter</h1>
<h2>Jayani Mannam — Spring 2026</h2>


==The Main Idea==
<hr>
:These interactions can be modeled using a familiar object: a spring! Every spring has a relaxed length. Once the spring is stretched past this relaxed length, the spring exerts a force inward in order to restore it to its original, relaxed length. If the spring is compressed smaller than the relaxed length, then it exerts a force outward in order to restore it to its original, relaxed length. Atoms can be modeled in the same way. Suppose two balls represent two atoms. These two atoms are connected by a spring, which represents the chemical bond (a bond between atoms that is a result of electric forces that bind protons and electrons to each other) between the two atoms (5).


:This simple model of just two atoms connected by a spring can then be extended to represent an entire solid! If each atom was connected to six atoms (each at a 90 degree angle to the others), it would form a cube, or what’s known as a simple cubic lattice. A simple cubic lattice is the simplest model of a solid. A model of this can be seen below. [[File: Ball Spring Trim 1.gif| thumb | 250px | right |Example of Coded Simple Cubic Lattice]] If one of the atoms is moved in any direction, it has a ripple effect with the other atoms (causing some to move away and some to come closer). Because they are all connected by these springs (which again, represents the electric forces between the atoms through a chemical bond), the movement of one atom affects all the other atoms (4).
<h2>1. Introduction</h2>


===Macroscopic Stretch===
The ball and spring model of matter is a simplified model used to understand how atoms in a solid interact with each other.  
:That force exerted by a spring, mentioned earlier, in order to restore it back to its relaxed length after it has been compressed or stretched can be calculated. Every spring has a specific stiffness that is unique to that particular spring known as the spring constant (stiffness), k. It has units N/m. The stretch or compression of the spring is modeled by the difference in the final length minus the initial (relaxed length) and is denoted as ∆L = L – L0, where L is the final length and L0 is the relaxed length. Change in L is referred to as the macroscopic stretch. The force the spring exerts is a direct relationship between the spring stiffness and the macroscopic stretch, and is therefore:
In this model, atoms are represented as balls, and the bonds between atoms are represented as springs.
:::Fspring = k*∆L


===Microscopic Stretch===
Although real atoms are not literally connected by tiny mechanical springs, this model is useful because it captures an important physical idea:
:The ball-and-spring model of a solid consists of balls that represent atoms held together by springs representing chemical bonds as a result of electric forces. As with a spring, a relaxed length is needed as a reference point to determine stretch or compression; the relaxed length of the spring between the two atoms is the center-to-center distance between the atoms, d. This distance would simply be twice the radius since the atoms would fill the space up entirely, minimizing the space between atoms. Since these bonds are modeled as springs, they also maintain a spring stiffness, k, specific to the spring and how much it will stretch. This stretch is defined as s, and is referred to as the microscopic stretch.  
when atoms are displaced from their equilibrium positions, forces act to restore them.


===Young’s modulus===
This model is especially helpful for understanding vibration, elasticity, wave motion, and energy storage in matter.
:Recall that the stiffness of a wire depends on how long the wire is and the thickness of that wire only; it does not depend on the material of wire, so different wires from the same material could still have different stiffnesses.  
It also connects microscopic motion to macroscopic material properties.
:Calculating Young’s modulus is a way to measure the “springiness” of a material and factor out the size and shape of the particular wire. Young’s modulus is the ratio of the force exerted per square meter of area (F/A) and the stretch of the wire (∆L/L). It is a way to measure the “springiness” of a material.


===A Mathematical Model===
At the microscopic level, atoms in solids vibrate around stable positions.
At the macroscopic level, those vibrations help explain material stiffness, sound propagation, heat transfer, and elastic deformation.


What are the mathematical equations that allow us to model this topic.  For example <math>{\frac{d\vec{p}}{dt}}_{system} = \vec{F}_{net}</math> where '''p''' is the momentum of the system and '''F''' is the net force from the surroundings.
<hr>


===A Computational Model===


This model of a solid can be simulated using VPython (written by Bruce Sherwood and licensed under Creative Commons 4.0). 
scene.append_to_caption("""


<code>
<h2>2. Hooke's Law</h2>
    N = 3
    k = 1
    m = 1
    spacing = 1
    atom_radius = 0.3*spacing
    L0 = spacing-1.8*atom_radius
    V0 = pi*(0.5*atom_radius)**2*L0
    scene.background = color.white
    scene.center = 0.5*(N-1)*vector(1,1,1)
    dt = 0.04*(2*pi*sqrt(m/k))
    axes = [vector(1,0,0), vector(0,1,0), vector(0,0,1)]
   
    class crystal:
        def atomAt(self, np):
            if (np.x>=0 and np.y>=0 and np.z>=0 and np.x<N and np.y<N and np.z<N):
            return self.atoms[int(np.x + np.y*N + np.z*N*N)]
            w = box()
            w.visible = False 
            w.radius = atom_radius
            w.pos = np*spacing
            w.momentum = vector(0,0,0)
            return w
       
        def __init__(self,  N, atom_radius, spacing, momentumRange ):
            self.atoms = []
            self.springs = []


            for z in range(N):
The main equation behind the model is:
                for y in range(N):
                    for x in range(N):
                        atom = sphere()
                        atom.pos = vector(x,y,z)*spacing
                        atom.radius = atom_radius
                        atom.color = vector(0,0.58,0.69)
                        px = 2*random()-1 # ranges from -1 to +1
                        py = 2*random()-1
                        pz = 2*random()-1
                        atom.momentum = momentumRange*vector(px,py,pz)
                        self.atoms.append( atom )       


            for d in range(3):
<br><br>
                for z in range(-1,N):
<b>F_spring = -k Δs</b>
                    for y in range(-1,N):
<br><br>
                        for x in range(-1,N):
                            atom = self.atomAt(vector(x,y,z))
                            neighbor = self.atomAt(vector(x,y,z)+axes[d])
                            if (atom.visible or neighbor.visible):
                                spring = helix()
                                spring.visible = atom.visible and neighbor.visible
                                spring.thickness = 0.05
                                spring.radius = 0.5*atom_radius
                                spring.length = spacing
                                spring.up = vector(1,1,1)
                                spring.atoms = [ atom, neighbor ]
                                spring.color = vector(1,0.5,0)
                                self.springs.append(spring)


[[File: Imageedit_3_2948105260.gif| thumb | 200px | right |Example of Coded Simple Cubic Lattice]]
where:


    c = crystal(N, atom_radius, spacing, 0.1*spacing*sqrt(k/m))
<ul>
    while True:
<li><b>F_spring</b> = restoring force</li>
        rate(30)
<li><b>k</b> = spring constant</li>
        for atom in c.atoms:
<li><b>Δs</b> = displacement from equilibrium</li>
            atom.pos = atom.pos + atom.momentum/m*dt
</ul>
        for spring in c.springs:
            spring.axis = spring.atoms[1].pos - spring.atoms[0].pos
            L = mag(spring.axis)
            spring.axis = spring.axis.norm()
            spring.pos = spring.atoms[0].pos+0.5*atom_radius*spring.axis
            Ls = L-1*atom_radius
            spring.length = Ls
            Fdt = spring.axis * (k*dt * (1-spacing/L))
            spring.atoms[0].momentum = spring.atoms[0].momentum + Fdt
            spring.atoms[1].momentum = spring.atoms[1].momentum - Fdt


==Examples==
The negative sign means the force acts opposite the displacement.


Be sure to show all steps in your solution and include diagrams whenever possible
If the atom moves right, the spring pulls left.
If the atom moves downward, the spring pulls upward.


===Simple===
For a vertical spring:
[[File:simple ex.png| 200px |right]]


Young’s modulus for brass is 8.96*10^11 Pa. A 100N weight is attached to a 5m length of brass wire. Find the increase in length of the wire. The diameter is 1.5mm.
<br><br>
<b>F_spring = -k(L - L0)L_hat</b>
<br><br>


[[File:Simple exam.png]]
where:
<ul>
<li><b>L</b> = current spring length</li>
<li><b>L0</b> = natural spring length</li>
<li><b>L_hat</b> = unit vector along spring</li>
</ul>


===Middling===
<hr>
You stack 1000 pennies, face to face, and apply a force of 25,000N to the top penny. While the force is applied, you find the thickness of the stack decrease by 1mm. The diameter of a penny is 1.905e-2m. The thickness of a penny is 1.3*10^-3m).


A) Calculate Young’s modulus of zinc.  
<h2>3. Harmonic Approximation</h2>


[[File:Part 1.png]]
The model works best for small displacements.


B) Calculate the interatomic spring stiffness for zinc. The diameter of a single zinc atom is 2.5 * 10^-10m.  
Real atomic forces are more complex, but near equilibrium the potential energy curve behaves like a parabola.


[[file:Part 2.png]]
That is why Hooke's Law gives a good approximation for vibrating atoms.


===Difficult===
<hr>
A wire made of an unknown alloy hangs from a support in the ceiling. You measure the relaxed length of the wire to be 1.8m long, and the radius of the wire to be 0.00025m. When you hang a 5kg mass from the wire, you measure that it stretches a distance of 4 x 10^-3m. the average bond length between atoms is 1.3 x 10^-10m for this alloy.
""")


A) If you treat the wire as a macroscopic spring, what is the overall spring stiffness of the wire?
= Mass-Spring System Simulation =


[[File:Part A.png]]
== Overview ==


B) What is the value of Young’s modulus for this alloy?
This project models the motion of a ball attached to a vertical spring. The top of the spring is fixed to a support, while the ball hangs from the lower end. When the ball is pulled downward and released, it oscillates because the spring pulls upward while gravity pulls downward.


[[File: New2.png]]
The simulation uses numerical integration to update the ball’s momentum and position over small time intervals. Instead of solving the motion analytically, the computer repeatedly applies Newton’s Second Law to predict the next state of the system.


C) What is the stiffness of a typical interatomic bond in the alloy?
This model demonstrates how forces create acceleration, how springs store energy, and how oscillatory motion can be represented computationally.


[[File: New3.png]]
== Physical Setup ==


D) You cut the wire into four pieces of equal length and hang a 5kg mass from one of the wires. How much does this wire stretch?
The system contains three visible objects:


[[File: New4.png]]
* '''Holder''' – a fixed support at the top.
* '''Spring''' – connects the holder to the ball.
* '''Ball''' – the moving mass.


E) You bundle (side by side) each of the four pieces of wire together and hang a 5kg mass from the bundle. How much does this bundle of wires stretch?
The values used in the simulation are:


[[File: New5.png]]
* Unstretched spring length: <math>L_0 = 0.10 \text{ m}</math>
* Spring constant: <math>k = 15 \text{ N/m}</math>
* Ball mass: <math>m = 0.10 \text{ kg}</math>
* Gravitational field: <math>\vec{g} = \langle 0,-9.8,0\rangle \text{ m/s}^2</math>


==Connectedness==
The ball begins slightly below equilibrium so the spring is stretched at the start.
1. How is this topic connected to something that you are interested in?
One of my favorite hobbies is baking. I bake something just about every day and will make anything from snickerdoodles to lemon bars to banana bread. An important factor that often determines the success of the final product is its consistency. A brownie that is too soft won't support its own geometry, yet one too hard will be difficult to chew. The stiffness and consistency of a baked good is determined by its Young's modulus, which is a function of the applied stress and the resulting strain of the material.


== Coordinate System ==


The simulation uses a 3D coordinate system:


2. How is it connected to your major?
* Positive x-direction = horizontal right
I am majoring in Biomedical Engineering, a major in which we use different biomaterials. When developing devices for medical purposes (crutches, prosthetics, wheelchairs, etc.) from biomaterials or using biomaterials to replace parts of the body (i.e. bone structure replacement), there are many factors to take in to consideration. One of these factors deals with Young's modulus, stress, and strain. Understanding the ball-and-spring model and subsequently Young's modulus is important to ensure that the proper materials are used for the appropriate devices and won't degrade, break, or deform when in use.
* Positive y-direction = upward
* Positive z-direction = outward/inward depth


Since the motion is vertical, most changes occur only in the y-direction.


== Model 1: Position Vector ==


3. Is there an interesting industrial application?
The position of the ball is tracked using a vector:
Yes! As mentioned before, one of the ways that the ball and spring model of a solid (and specifically Young's modulus, stress, and strain that stem from the model) related to the biomedical industry is through bone structure replacement. Using this model, it can be determined whether or not the biomaterial has similar deformable properties with the material it will replace. Typically, these materials need high Young's modulus' because they bear a high amount of force. Therefore, a selected biomaterial can be determined a good fit for replacement if it's Young's modulus is similar to bone.
 
<math>\vec{r}_{ball}</math>
 
The holder also has a fixed position:
 
<math>\vec{r}_{holder}</math>
 
The spring vector is the displacement from the holder to the ball:
 
<math>\vec{L} = \vec{r}_{ball} - \vec{r}_{holder}</math>
 
This gives both the direction and length of the spring at every moment.
 
== Model 2: Spring Length ==
 
The current spring length is the magnitude of the spring vector:
 
<math>|\vec{L}| = \sqrt{L_x^2 + L_y^2 + L_z^2}</math>
 
If the spring length equals <math>L_0</math>, the spring is unstretched.
 
If:
 
<math>|\vec{L}| > L_0</math>
 
the spring is stretched.
 
If:
 
<math>|\vec{L}| < L_0</math>
 
the spring is compressed.
 
== Model 3: Hooke's Law ==
 
The spring exerts a restoring force that tries to return the spring to its natural length.
 
The force is:
 
<math>\vec{F}_s = -k(|\vec{L}| - L_0)\hat{L}</math>
 
where:
 
<math>\hat{L} = \frac{\vec{L}}{|\vec{L}|}</math>
 
Explanation:
 
* <math>k</math> determines spring stiffness.
* <math>(|\vec{L}| - L_0)</math> is the stretch/compression amount.
* <math>\hat{L}</math> gives the direction.
* The negative sign means the force acts opposite displacement.
 
If the ball is pulled downward, the spring force points upward.
 
== Model 4: Gravity ==
 
Gravity constantly pulls the ball downward.
 
The gravitational force is:
 
<math>\vec{F}_g = m\vec{g}</math>
 
Substituting values:
 
<math>\vec{F}_g = (0.10)\langle0,-9.8,0\rangle</math>
 
<math>\vec{F}_g = \langle0,-0.98,0\rangle \text{ N}</math>
 
This force remains constant throughout the simulation.
 
== Model 5: Net Force ==
 
The total force on the ball is the sum of spring force and gravity:
 
<math>\vec{F}_{net} = \vec{F}_s + \vec{F}_g</math>
 
This determines how the ball accelerates.
 
Cases:
 
* If spring force is larger than gravity, the ball accelerates upward.
* If gravity is larger, the ball accelerates downward.
* If forces balance, acceleration is zero.
 
== Model 6: Newton's Second Law ==
 
Newton’s Second Law states:
 
<math>\vec{F}_{net} = m\vec{a}</math>
 
So acceleration is:
 
<math>\vec{a} = \frac{\vec{F}_{net}}{m}</math>
 
Rather than directly solving for acceleration each time, the program updates momentum.
 
== Model 7: Momentum Principle ==
 
Momentum is:
 
<math>\vec{p} = m\vec{v}</math>
 
The Momentum Principle gives:
 
<math>\Delta \vec{p} = \vec{F}_{net}\Delta t</math>
 
So each loop:
 
<math>\vec{p}_{new} = \vec{p}_{old} + \vec{F}_{net}\Delta t</math>
 
This is the main physics engine of the simulation.
 
== Model 8: Position Update ==
 
Velocity is found from momentum:
 
<math>\vec{v} = \frac{\vec{p}}{m}</math>
 
Then position updates as:
 
<math>\vec{r}_{new} = \vec{r}_{old} + \vec{v}\Delta t</math>
 
Substituting velocity:
 
<math>\vec{r}_{new} = \vec{r}_{old} + \left(\frac{\vec{p}}{m}\right)\Delta t</math>
 
This allows the ball to move frame by frame.
 
== Model 9: Numerical Integration ==
 
The simulation does not solve one large equation for all time. Instead, it uses many tiny time steps:
 
<math>\Delta t = 0.01 \text{ s}</math>
 
For each step:
 
# Calculate spring vector
# Calculate spring force
# Calculate gravity
# Find net force
# Update momentum
# Update position
# Redraw spring
 
This method is called Euler-Cromer style numerical integration.
 
Smaller time steps improve accuracy.
 
== Model 10: Oscillatory Motion ==
 
Because the spring force always points toward equilibrium, the ball repeatedly moves up and down.
 
Sequence:
 
# Ball released downward
# Spring pulls upward
# Ball speeds upward
# Ball passes equilibrium
# Gravity and spring slow it
# Ball reverses
# Motion repeats
 
This creates periodic motion.
 
== Model 11: Equilibrium Position ==
 
Equilibrium occurs when:
 
<math>\vec{F}_s + \vec{F}_g = 0</math>
 
In magnitude form:
 
<math>k\Delta L = mg</math>
 
So:
 
<math>\Delta L = \frac{mg}{k}</math>
 
Substitute values:
 
<math>\Delta L = \frac{0.10(9.8)}{15}</math>
 
<math>\Delta L = 0.0653 \text{ m}</math>
 
So the ball hangs about 6.53 cm below the unstretched position at equilibrium.
 
== Model 12: Energy in the System ==
 
Energy shifts between gravitational, elastic, and kinetic forms.
 
Spring potential energy:
 
<math>U_s = \frac{1}{2}k(\Delta L)^2</math>
 
Gravitational potential energy:
 
<math>U_g = mgy</math>
 
Kinetic energy:
 
<math>K = \frac{1}{2}mv^2</math>
 
As the ball moves:
 
* At highest point: mostly potential energy
* At equilibrium: maximum speed
* At lowest point: spring energy largest
 
== Model 13: Why the Motion Continues ==
 
The ball has inertia. Even when it reaches equilibrium, it does not stop instantly because it still has momentum.
 
That momentum carries it past equilibrium, stretching/compressing the spring in the opposite direction. Then the restoring force reverses the motion.
 
This repeated exchange creates oscillation.
 
== Model 14: Code Connection ==
 
Each major equation appears directly in the code:
 
Spring vector:
 
<pre>
L = ball.pos - holder.pos
</pre>
 
Spring force:
 
<pre>
Fs = -k * (mag(L) - L0) * norm(L)
</pre>
 
Gravity:
 
<pre>
Fg = ball.m * g
</pre>
 
Net force:
 
<pre>
Fnet = Fs + Fg
</pre>
 
Momentum update:
 
<pre>
ball.p = ball.p + Fnet * dt
</pre>
 
Position update:
 
<pre>
ball.pos = ball.pos + (ball.p / ball.m) * dt
</pre>
 
Spring redraw:
 
<pre>
spring.axis = ball.pos - holder.pos
</pre>
 
== Real-World Applications ==
 
Mass-spring systems appear in many engineering designs:
 
* Vehicle suspension systems
* Building vibration control
* Shock absorbers
* Scales
* Seismographs
* Robotics
* Mechanical sensors
 
The same principles also apply to atoms vibrating in solids.
 
== Theory of Simple Harmonic Motion ==
 
A mass-spring system is one of the most common examples of simple harmonic motion (SHM). SHM occurs when the restoring force is directly proportional to displacement from equilibrium and always points back toward equilibrium.
 
Mathematically:
 
<math>F = -kx</math>
 
where:
 
* <math>F</math> = restoring force 
* <math>k</math> = spring constant 
* <math>x</math> = displacement from equilibrium 
 
The negative sign shows the force acts opposite the displacement.
 
This relationship causes repeating motion.
 
== Period of a Mass-Spring System ==
 
The time for one complete oscillation is called the period.
 
For an ideal spring-mass system:
 
<math>T = 2\pi\sqrt{\frac{m}{k}}</math>
 
where:
 
* <math>T</math> = period (seconds) 
* <math>m</math> = mass 
* <math>k</math> = spring constant 
 
Key ideas:
 
* Larger mass increases the period.
* Stiffer springs decrease the period.
* Period does not depend on amplitude for ideal SHM.
 
== Frequency ==
 
Frequency is the number of cycles completed each second.
 
<math>f = \frac{1}{T}</math>
 
Unit:
 
<math>\text{Hz}</math> (hertz)
 
Higher frequency means faster oscillations.
 
== Angular Frequency ==
 
Oscillating systems are often described using angular frequency:
 
<math>\omega = \sqrt{\frac{k}{m}}</math>
 
Also:
 
<math>\omega = 2\pi f</math>
 
Angular frequency is useful in sinusoidal equations of motion.
 
== Position as a Function of Time ==
 
Ideal spring motion can be modeled by:
 
<math>x(t)=A\cos(\omega t+\phi)</math>
 
where:
 
* <math>A</math> = amplitude 
* <math>\omega</math> = angular frequency 
* <math>t</math> = time 
* <math>\phi</math> = phase constant 
 
This equation describes how position changes smoothly over time.
 
== Velocity in SHM ==
 
Velocity is the derivative of position:
 
<math>v(t) = -A\omega\sin(\omega t+\phi)</math>
 
Maximum velocity occurs at equilibrium because all available energy is kinetic there.
 
Velocity is zero at the turning points.
 
== Acceleration in SHM ==
 
Acceleration is:
 
<math>a(t) = -\omega^2 x</math>
 
This shows acceleration is always proportional to displacement and points toward equilibrium.
 
Maximum acceleration occurs at maximum displacement.
 
== Energy in Simple Harmonic Motion ==
 
Total mechanical energy remains constant in an ideal system.
 
Kinetic energy:
 
<math>K = \frac{1}{2}mv^2</math>
 
Spring potential energy:
 
<math>U = \frac{1}{2}kx^2</math>
 
Total energy:
 
<math>E = K + U = \frac{1}{2}kA^2</math>
 
As the object moves:
 
* At equilibrium: kinetic energy is maximum.
* At maximum displacement: potential energy is maximum.
* Total energy stays constant.
 
== Damping ==
 
Real systems lose energy due to friction or air resistance. This causes damping.
 
Effects of damping:
 
* Amplitude decreases over time.
* Oscillations eventually stop.
* Mechanical energy is converted to thermal energy.
 
Examples:
 
* Car suspension systems
* Door closers
* Shock absorbers
 
== Resonance ==
 
Resonance occurs when an external driving force matches the natural frequency of a system.
 
When this happens:
 
* Amplitude grows significantly.
* Energy transfer becomes highly efficient.
 
Examples:
 
* Playground swings
* Musical instruments
* Bridges and buildings under vibration
 
Resonance can be useful or dangerous depending on the system.
 
== Hooke's Law Limits ==
 
Hooke’s Law works only when the spring remains in its elastic range.
 
If stretched too far:
 
* The spring may deform permanently.
* Force may no longer be proportional to displacement.
* Motion no longer follows ideal SHM.
 
This is why real materials have operating limits.
 
== Relationship to Circular Motion ==
 
Simple harmonic motion can be viewed as the projection of uniform circular motion onto one axis.
 
If a point moves in a circle with constant angular speed, its horizontal or vertical shadow moves sinusoidally.
 
This is why sine and cosine functions describe oscillations.
 
== Real-World Applications ==
 
Mass-spring concepts are used in:
 
* Vehicle suspension
* Earthquake engineering
* Robotics
* Mechanical watches
* Vibration isolation systems
* Industrial sensors
* Microelectromechanical systems (MEMS)
* Seismographs
 
Understanding oscillation helps engineers design stable systems.
 
== Experimental Quantities ==
 
Common values measured in spring experiments:
 
* Amplitude – maximum displacement
* Period – time for one cycle
* Frequency – cycles per second
* Equilibrium position
* Maximum speed
* Maximum acceleration
* Energy loss due to damping
 
These measurements help compare theory with real motion.
 
== Common Physics 2 Connections ==
 
Mass-spring systems connect to many topics:
 
* Newton’s Laws
* Energy conservation
* Differential equations
* Waves
* Resonance
* Damping
* Numerical modeling
* Momentum principle
 
They are foundational systems in introductory physics.
 
== Summary ==
 
Spring motion is one of the clearest examples of how forces create predictable motion. A restoring force causes repeated oscillation, while energy continuously shifts between kinetic and potential forms.
 
Because of this, mass-spring systems are used throughout physics, engineering, and modern technology.


==History==
==History==


Put this idea in historical context. Give the reader the Who, What, When, Where, and Why.
In the 1600's, Robert Hooke developed the idea of Hooke's Law, which states that for relatively small deformations of an object, the deformation is proportional to the force that deformed it. This lead to the development of the equation to calculate force of the spring; where it is equal to the spring stiffness multiplied by the change in length. This fundamental concept lead to other developments in Molecular Mechanics. However, the specific history for the ball-and-spring model is unknown, although the major developments in this area were accomplished in and around the 1930s and 1940s (including developments by T.L. Hill, I. Dostrovsky, and F. H. Westheimer) (6).


== See also ==
== See also ==
Line 182: Line 633:


===External links===
===External links===
#http://www.nuffieldfoundation.org/practical-physics/model-vibrating-atoms-solid
#https://www.physics.ncsu.edu/clarke/teaching/class.html
#http://umdberg.pbworks.com/w/page/46030513/A%20simple%20model%20of%20solid%20matter


==References==
==References==
Line 188: Line 642:
#http://bulldog2.redlands.edu/facultyfolder/eric_hill/Phys231/Lecture/Lect%209.pdf
#http://bulldog2.redlands.edu/facultyfolder/eric_hill/Phys231/Lecture/Lect%209.pdf
#http://www.matterandinteractions.org/Content/Materials/materials.html
#http://www.matterandinteractions.org/Content/Materials/materials.html
#Chabay, Ruth W., and Bruce A. Sherwood. Matter & Interactions. 3rd ed.  
#Chabay, Ruth W., and Bruce A. Sherwood. Matter & Interactions. 3rd ed.
#http://www.sdsc.edu/~kimb/history.html






[[Category:Which Category did you place this in?]]
[[Category:Energy]]

Latest revision as of 20:34, 26 April 2026

Page created by Ashley Fleck.

The interactions of atoms can be modeled using balls to represent the atoms and a spring to represent the chemical bond between them.

scene.title = "Ball and Spring Model of Matter\n" scene.width = 900 scene.height = 600 scene.background = color.white

scene.append_to_caption("""

Ball and Spring Model of Matter

Jayani Mannam — Spring 2026


1. Introduction

The ball and spring model of matter is a simplified model used to understand how atoms in a solid interact with each other. In this model, atoms are represented as balls, and the bonds between atoms are represented as springs.

Although real atoms are not literally connected by tiny mechanical springs, this model is useful because it captures an important physical idea: when atoms are displaced from their equilibrium positions, forces act to restore them.

This model is especially helpful for understanding vibration, elasticity, wave motion, and energy storage in matter. It also connects microscopic motion to macroscopic material properties.

At the microscopic level, atoms in solids vibrate around stable positions. At the macroscopic level, those vibrations help explain material stiffness, sound propagation, heat transfer, and elastic deformation.



scene.append_to_caption("""

2. Hooke's Law

The main equation behind the model is:



F_spring = -k Δs

where:

  • F_spring = restoring force
  • k = spring constant
  • Δs = displacement from equilibrium

The negative sign means the force acts opposite the displacement.

If the atom moves right, the spring pulls left. If the atom moves downward, the spring pulls upward.

For a vertical spring:



F_spring = -k(L - L0)L_hat

where:

  • L = current spring length
  • L0 = natural spring length
  • L_hat = unit vector along spring

3. Harmonic Approximation

The model works best for small displacements.

Real atomic forces are more complex, but near equilibrium the potential energy curve behaves like a parabola.

That is why Hooke's Law gives a good approximation for vibrating atoms.


""")

Mass-Spring System Simulation

Overview

This project models the motion of a ball attached to a vertical spring. The top of the spring is fixed to a support, while the ball hangs from the lower end. When the ball is pulled downward and released, it oscillates because the spring pulls upward while gravity pulls downward.

The simulation uses numerical integration to update the ball’s momentum and position over small time intervals. Instead of solving the motion analytically, the computer repeatedly applies Newton’s Second Law to predict the next state of the system.

This model demonstrates how forces create acceleration, how springs store energy, and how oscillatory motion can be represented computationally.

Physical Setup

The system contains three visible objects:

  • Holder – a fixed support at the top.
  • Spring – connects the holder to the ball.
  • Ball – the moving mass.

The values used in the simulation are:

  • Unstretched spring length: [math]\displaystyle{ L_0 = 0.10 \text{ m} }[/math]
  • Spring constant: [math]\displaystyle{ k = 15 \text{ N/m} }[/math]
  • Ball mass: [math]\displaystyle{ m = 0.10 \text{ kg} }[/math]
  • Gravitational field: [math]\displaystyle{ \vec{g} = \langle 0,-9.8,0\rangle \text{ m/s}^2 }[/math]

The ball begins slightly below equilibrium so the spring is stretched at the start.

Coordinate System

The simulation uses a 3D coordinate system:

  • Positive x-direction = horizontal right
  • Positive y-direction = upward
  • Positive z-direction = outward/inward depth

Since the motion is vertical, most changes occur only in the y-direction.

Model 1: Position Vector

The position of the ball is tracked using a vector:

[math]\displaystyle{ \vec{r}_{ball} }[/math]

The holder also has a fixed position:

[math]\displaystyle{ \vec{r}_{holder} }[/math]

The spring vector is the displacement from the holder to the ball:

[math]\displaystyle{ \vec{L} = \vec{r}_{ball} - \vec{r}_{holder} }[/math]

This gives both the direction and length of the spring at every moment.

Model 2: Spring Length

The current spring length is the magnitude of the spring vector:

[math]\displaystyle{ |\vec{L}| = \sqrt{L_x^2 + L_y^2 + L_z^2} }[/math]

If the spring length equals [math]\displaystyle{ L_0 }[/math], the spring is unstretched.

If:

[math]\displaystyle{ |\vec{L}| \gt L_0 }[/math]

the spring is stretched.

If:

[math]\displaystyle{ |\vec{L}| \lt L_0 }[/math]

the spring is compressed.

Model 3: Hooke's Law

The spring exerts a restoring force that tries to return the spring to its natural length.

The force is:

[math]\displaystyle{ \vec{F}_s = -k(|\vec{L}| - L_0)\hat{L} }[/math]

where:

[math]\displaystyle{ \hat{L} = \frac{\vec{L}}{|\vec{L}|} }[/math]

Explanation:

  • [math]\displaystyle{ k }[/math] determines spring stiffness.
  • [math]\displaystyle{ (|\vec{L}| - L_0) }[/math] is the stretch/compression amount.
  • [math]\displaystyle{ \hat{L} }[/math] gives the direction.
  • The negative sign means the force acts opposite displacement.

If the ball is pulled downward, the spring force points upward.

Model 4: Gravity

Gravity constantly pulls the ball downward.

The gravitational force is:

[math]\displaystyle{ \vec{F}_g = m\vec{g} }[/math]

Substituting values:

[math]\displaystyle{ \vec{F}_g = (0.10)\langle0,-9.8,0\rangle }[/math]

[math]\displaystyle{ \vec{F}_g = \langle0,-0.98,0\rangle \text{ N} }[/math]

This force remains constant throughout the simulation.

Model 5: Net Force

The total force on the ball is the sum of spring force and gravity:

[math]\displaystyle{ \vec{F}_{net} = \vec{F}_s + \vec{F}_g }[/math]

This determines how the ball accelerates.

Cases:

  • If spring force is larger than gravity, the ball accelerates upward.
  • If gravity is larger, the ball accelerates downward.
  • If forces balance, acceleration is zero.

Model 6: Newton's Second Law

Newton’s Second Law states:

[math]\displaystyle{ \vec{F}_{net} = m\vec{a} }[/math]

So acceleration is:

[math]\displaystyle{ \vec{a} = \frac{\vec{F}_{net}}{m} }[/math]

Rather than directly solving for acceleration each time, the program updates momentum.

Model 7: Momentum Principle

Momentum is:

[math]\displaystyle{ \vec{p} = m\vec{v} }[/math]

The Momentum Principle gives:

[math]\displaystyle{ \Delta \vec{p} = \vec{F}_{net}\Delta t }[/math]

So each loop:

[math]\displaystyle{ \vec{p}_{new} = \vec{p}_{old} + \vec{F}_{net}\Delta t }[/math]

This is the main physics engine of the simulation.

Model 8: Position Update

Velocity is found from momentum:

[math]\displaystyle{ \vec{v} = \frac{\vec{p}}{m} }[/math]

Then position updates as:

[math]\displaystyle{ \vec{r}_{new} = \vec{r}_{old} + \vec{v}\Delta t }[/math]

Substituting velocity:

[math]\displaystyle{ \vec{r}_{new} = \vec{r}_{old} + \left(\frac{\vec{p}}{m}\right)\Delta t }[/math]

This allows the ball to move frame by frame.

Model 9: Numerical Integration

The simulation does not solve one large equation for all time. Instead, it uses many tiny time steps:

[math]\displaystyle{ \Delta t = 0.01 \text{ s} }[/math]

For each step:

  1. Calculate spring vector
  2. Calculate spring force
  3. Calculate gravity
  4. Find net force
  5. Update momentum
  6. Update position
  7. Redraw spring

This method is called Euler-Cromer style numerical integration.

Smaller time steps improve accuracy.

Model 10: Oscillatory Motion

Because the spring force always points toward equilibrium, the ball repeatedly moves up and down.

Sequence:

  1. Ball released downward
  2. Spring pulls upward
  3. Ball speeds upward
  4. Ball passes equilibrium
  5. Gravity and spring slow it
  6. Ball reverses
  7. Motion repeats

This creates periodic motion.

Model 11: Equilibrium Position

Equilibrium occurs when:

[math]\displaystyle{ \vec{F}_s + \vec{F}_g = 0 }[/math]

In magnitude form:

[math]\displaystyle{ k\Delta L = mg }[/math]

So:

[math]\displaystyle{ \Delta L = \frac{mg}{k} }[/math]

Substitute values:

[math]\displaystyle{ \Delta L = \frac{0.10(9.8)}{15} }[/math]

[math]\displaystyle{ \Delta L = 0.0653 \text{ m} }[/math]

So the ball hangs about 6.53 cm below the unstretched position at equilibrium.

Model 12: Energy in the System

Energy shifts between gravitational, elastic, and kinetic forms.

Spring potential energy:

[math]\displaystyle{ U_s = \frac{1}{2}k(\Delta L)^2 }[/math]

Gravitational potential energy:

[math]\displaystyle{ U_g = mgy }[/math]

Kinetic energy:

[math]\displaystyle{ K = \frac{1}{2}mv^2 }[/math]

As the ball moves:

  • At highest point: mostly potential energy
  • At equilibrium: maximum speed
  • At lowest point: spring energy largest

Model 13: Why the Motion Continues

The ball has inertia. Even when it reaches equilibrium, it does not stop instantly because it still has momentum.

That momentum carries it past equilibrium, stretching/compressing the spring in the opposite direction. Then the restoring force reverses the motion.

This repeated exchange creates oscillation.

Model 14: Code Connection

Each major equation appears directly in the code:

Spring vector:

L = ball.pos - holder.pos

Spring force:

Fs = -k * (mag(L) - L0) * norm(L)

Gravity:

Fg = ball.m * g

Net force:

Fnet = Fs + Fg

Momentum update:

ball.p = ball.p + Fnet * dt

Position update:

ball.pos = ball.pos + (ball.p / ball.m) * dt

Spring redraw:

spring.axis = ball.pos - holder.pos

Real-World Applications

Mass-spring systems appear in many engineering designs:

  • Vehicle suspension systems
  • Building vibration control
  • Shock absorbers
  • Scales
  • Seismographs
  • Robotics
  • Mechanical sensors

The same principles also apply to atoms vibrating in solids.

Theory of Simple Harmonic Motion

A mass-spring system is one of the most common examples of simple harmonic motion (SHM). SHM occurs when the restoring force is directly proportional to displacement from equilibrium and always points back toward equilibrium.

Mathematically:

[math]\displaystyle{ F = -kx }[/math]

where:

  • [math]\displaystyle{ F }[/math] = restoring force
  • [math]\displaystyle{ k }[/math] = spring constant
  • [math]\displaystyle{ x }[/math] = displacement from equilibrium

The negative sign shows the force acts opposite the displacement.

This relationship causes repeating motion.

Period of a Mass-Spring System

The time for one complete oscillation is called the period.

For an ideal spring-mass system:

[math]\displaystyle{ T = 2\pi\sqrt{\frac{m}{k}} }[/math]

where:

  • [math]\displaystyle{ T }[/math] = period (seconds)
  • [math]\displaystyle{ m }[/math] = mass
  • [math]\displaystyle{ k }[/math] = spring constant

Key ideas:

  • Larger mass increases the period.
  • Stiffer springs decrease the period.
  • Period does not depend on amplitude for ideal SHM.

Frequency

Frequency is the number of cycles completed each second.

[math]\displaystyle{ f = \frac{1}{T} }[/math]

Unit:

[math]\displaystyle{ \text{Hz} }[/math] (hertz)

Higher frequency means faster oscillations.

Angular Frequency

Oscillating systems are often described using angular frequency:

[math]\displaystyle{ \omega = \sqrt{\frac{k}{m}} }[/math]

Also:

[math]\displaystyle{ \omega = 2\pi f }[/math]

Angular frequency is useful in sinusoidal equations of motion.

Position as a Function of Time

Ideal spring motion can be modeled by:

[math]\displaystyle{ x(t)=A\cos(\omega t+\phi) }[/math]

where:

  • [math]\displaystyle{ A }[/math] = amplitude
  • [math]\displaystyle{ \omega }[/math] = angular frequency
  • [math]\displaystyle{ t }[/math] = time
  • [math]\displaystyle{ \phi }[/math] = phase constant

This equation describes how position changes smoothly over time.

Velocity in SHM

Velocity is the derivative of position:

[math]\displaystyle{ v(t) = -A\omega\sin(\omega t+\phi) }[/math]

Maximum velocity occurs at equilibrium because all available energy is kinetic there.

Velocity is zero at the turning points.

Acceleration in SHM

Acceleration is:

[math]\displaystyle{ a(t) = -\omega^2 x }[/math]

This shows acceleration is always proportional to displacement and points toward equilibrium.

Maximum acceleration occurs at maximum displacement.

Energy in Simple Harmonic Motion

Total mechanical energy remains constant in an ideal system.

Kinetic energy:

[math]\displaystyle{ K = \frac{1}{2}mv^2 }[/math]

Spring potential energy:

[math]\displaystyle{ U = \frac{1}{2}kx^2 }[/math]

Total energy:

[math]\displaystyle{ E = K + U = \frac{1}{2}kA^2 }[/math]

As the object moves:

  • At equilibrium: kinetic energy is maximum.
  • At maximum displacement: potential energy is maximum.
  • Total energy stays constant.

Damping

Real systems lose energy due to friction or air resistance. This causes damping.

Effects of damping:

  • Amplitude decreases over time.
  • Oscillations eventually stop.
  • Mechanical energy is converted to thermal energy.

Examples:

  • Car suspension systems
  • Door closers
  • Shock absorbers

Resonance

Resonance occurs when an external driving force matches the natural frequency of a system.

When this happens:

  • Amplitude grows significantly.
  • Energy transfer becomes highly efficient.

Examples:

  • Playground swings
  • Musical instruments
  • Bridges and buildings under vibration

Resonance can be useful or dangerous depending on the system.

Hooke's Law Limits

Hooke’s Law works only when the spring remains in its elastic range.

If stretched too far:

  • The spring may deform permanently.
  • Force may no longer be proportional to displacement.
  • Motion no longer follows ideal SHM.

This is why real materials have operating limits.

Relationship to Circular Motion

Simple harmonic motion can be viewed as the projection of uniform circular motion onto one axis.

If a point moves in a circle with constant angular speed, its horizontal or vertical shadow moves sinusoidally.

This is why sine and cosine functions describe oscillations.

Real-World Applications

Mass-spring concepts are used in:

  • Vehicle suspension
  • Earthquake engineering
  • Robotics
  • Mechanical watches
  • Vibration isolation systems
  • Industrial sensors
  • Microelectromechanical systems (MEMS)
  • Seismographs

Understanding oscillation helps engineers design stable systems.

Experimental Quantities

Common values measured in spring experiments:

  • Amplitude – maximum displacement
  • Period – time for one cycle
  • Frequency – cycles per second
  • Equilibrium position
  • Maximum speed
  • Maximum acceleration
  • Energy loss due to damping

These measurements help compare theory with real motion.

Common Physics 2 Connections

Mass-spring systems connect to many topics:

  • Newton’s Laws
  • Energy conservation
  • Differential equations
  • Waves
  • Resonance
  • Damping
  • Numerical modeling
  • Momentum principle

They are foundational systems in introductory physics.

Summary

Spring motion is one of the clearest examples of how forces create predictable motion. A restoring force causes repeated oscillation, while energy continuously shifts between kinetic and potential forms.

Because of this, mass-spring systems are used throughout physics, engineering, and modern technology.

History

In the 1600's, Robert Hooke developed the idea of Hooke's Law, which states that for relatively small deformations of an object, the deformation is proportional to the force that deformed it. This lead to the development of the equation to calculate force of the spring; where it is equal to the spring stiffness multiplied by the change in length. This fundamental concept lead to other developments in Molecular Mechanics. However, the specific history for the ball-and-spring model is unknown, although the major developments in this area were accomplished in and around the 1930s and 1940s (including developments by T.L. Hill, I. Dostrovsky, and F. H. Westheimer) (6).

See also

The Ball-and-Spring model deals primarily with Kinds of Matter, Young's Modulus, Hooke's Law, and the Length and Stiffness of an Interatomic Bond.

Further reading

The Kinetics of materials is a book written by Robert W. Balluffi that goes in-depth about many different concepts, including the ball-and-spring model and Young's modulus.

The Molecular Dynamics and Spectroscopy by Stimulated Emission Pumping, written by Hai-Lung Dai and Robert W. Field, looks at a different aspect of the ball-and-spring model and determines its role in quantum eigenstate spectra.

External links

  1. http://www.nuffieldfoundation.org/practical-physics/model-vibrating-atoms-solid
  2. https://www.physics.ncsu.edu/clarke/teaching/class.html
  3. http://umdberg.pbworks.com/w/page/46030513/A%20simple%20model%20of%20solid%20matter

References

  1. http://p3server.pa.msu.edu/coursewiki/doku.php?id=183_notes:model_of_solids
  2. http://webs.morningside.edu/slaven/Physics/atom/atom7.html
  3. http://bulldog2.redlands.edu/facultyfolder/eric_hill/Phys231/Lecture/Lect%209.pdf
  4. http://www.matterandinteractions.org/Content/Materials/materials.html
  5. Chabay, Ruth W., and Bruce A. Sherwood. Matter & Interactions. 3rd ed.
  6. http://www.sdsc.edu/~kimb/history.html