91Ó°ÊÓ

A forced damped spring-mass system (Fig. \(\mathrm{P} 28.38\) ) has the following ordinary differential equation of motion: \\[ m \frac{d^{2} x}{d t^{2}}+a\left|\frac{d x}{d t}\right| \frac{d x}{d t}+k x=F_{o} \sin (\omega t) \\] where \(x=\) displacement from the equilibrium position, \(t=\) time, \(m=2 \mathrm{kg}\) mass \(, a=5 \mathrm{N} /(\mathrm{m} / \mathrm{s})^{2}\), and \(k=6 \mathrm{N} /(\mathrm{m} / \mathrm{s})^{2}\). The damping term is nonlinear and represents air damping. The forcing function \(F_{0} \sin (\omega t)\) has values of \(F_{o}=2.5 \mathrm{N}\) and \(\omega=0.5 \mathrm{rad} / \mathrm{sec} .\) The initial conditions are Initial velocity' \\[ \frac{d x}{d t}=0 \mathrm{m} / \mathrm{s} \\] Initial displacement \\[ x=1 \mathrm{m} \\] Solve this equation using a numerical method over the time period \(0 \leq t \leq 15\) s. Plot the displacement and velocity versus time, and plot the forcing function on the same curve. Also, develop a separate plot of velocity versus displacement.

Step by step solution

01

Convert the second-order ODE into two first-order ODEs

To do this, we let the velocity (first derivative of x with respect to t) be another function y. Then the system may be described by: \(y = \frac{dx}{dt}\) \(m \frac{d^{2} x}{d t^{2}} = m \frac{dy}{dt} \) Now, we rewrite the original ODE in terms of x, y, and their derivatives: \(m \frac{dy}{dt} + a|y| y + kx = F_0 \sin(\omega t)\)

02

Define the numerical method for solving the ODEs

Now we have two first-order ODEs which can be written as: \( \frac{dx}{dt} = y\) \( \frac{dy}{dt} = -\frac{a}{m} |y| y - \frac{k}{m}x + \frac{F_0}{m} \sin(\omega t)\) The commonly used numerical method is Euler's method. However, to obtain better results we will use the Runge-Kutta method of order 4. Create a function for solving the system of ODEs using the 4th order Runge-Kutta method.

03

Define the given parameters and initial conditions

Given parameters: \(m = 2kg, a = 5N/(m/s)^2, k = 6N/(m/s)^2, F_0 = 2.5N, \omega = 0.5rad/s\) Initial conditions: \(x(0) = 1m, y(0) = v(0) = 0m/s\) Define a time period for the numerical solution: \(0 \leq t \leq 15s\)

04

Solve the ODEs using the defined method

Using Runge-Kutta 4th order method with given parameters and initial conditions, solve for x(t) and y(t) (velocity). Choose a reasonable time step size for the numerical method, e.g., 0.01s.

05

Plot the results according to the requests

1. Create a plot for the displacement x(t) vs time t. 2. Create a plot for the velocity y(t) vs time t. On the same plot, also show the forcing function \(F_0sin(\omega t)\). 3. Create a third plot for the velocity y(t) vs displacement x(t). With these steps, the given forced damped spring-mass system would be solved numerically, and the results can be visualized and analyzed through the generated plots.

Unlock Step-by-Step Solutions & Ace Your Exams!

• Full Textbook Solutions

  Get detailed explanations and key concepts

• Unlimited Al creation

  Al flashcards, explanations, exams and more...

• Ads-free access

  To over 500 millions flashcards

• Money-back guarantee

  We refund you if you fail your exam.

Over 30 million students worldwide already upgrade their learning with 91Ó°ÊÓ!

Key Concepts

These are the key concepts you need to understand to accurately answer the question.

Runge-Kutta Method

The Runge-Kutta method is one of the most popular numerical techniques for solving ordinary differential equations (ODEs). It's particularly useful when dealing with complex systems like the forced damped spring-mass system. Unlike basic methods such as Euler's, Runge-Kutta provides higher accuracy by taking multiple intermediate calculations within a single step.

The 4th Order Runge-Kutta method, often just referred to as RK4, takes four approximations within a time step to estimate the next value. These intermediate steps allow it to handle nonlinear terms effectively, which is crucial for systems with nonlinear damping like ours. The RK4 formulas are:

• Calculate slope estimates: \(k_1 = f(x_n, y_n)\), \(k_2 = f(x_n + h/2, y_n + k_1/2)\), \(k_3 = f(x_n + h/2, y_n + k_2/2)\), \(k_4 = f(x_n + h, y_n + k_3)\)
• Update next value: \(y_{n+1} = y_n + (h/6)(k_1 + 2k_2 + 2k_3 + k_4)\)
• Here, \(h\) is the time step size.

This method is most accurate and stable when small step sizes are used, such as in our exercise where we utilized a step size of 0.01s.

Ordinary Differential Equations

Ordinary Differential Equations (ODEs) are equations involving derivatives of a function with respect to one independent variable, usually time. These equations are fundamental in depicting how physical processes evolve over time.

In our forced damped spring-mass system, the ODE describes the motion of a mass on a spring as influenced by damping and an external force. The original ODE provided:

• Second-order ODE: \(m \frac{d^{2} x}{d t^{2}} + a\left|\frac{d x}{d t}\right| \frac{d x}{d t} + k x = F_{o} \sin (\omega t)\)

To solve this, it’s converted to two first-order ODEs, making it easier to apply numerical methods. This conversion involves defining velocity as a separate function, adding dimensionality to the problem. For numerical computation, managing ODEs in first-order form provides more flexibility and simplification.

Forced Damped Spring-Mass System

The forced damped spring-mass system is a classic mechanical model used to understand the dynamics of mechanical oscillators. This involves a mass attached to a spring that experiences both damping forces and an external periodic force.

Key elements of the system:

• The **mass (m)**: Determines the system's inertia, affecting how it responds to forces.
• The **spring constant (k)**: Represents the spring's stiffness. A higher value means a stiffer spring that resists displacement.
• The **damping (a)**: In this case, nonlinear air damping, which dissipates energy and reduces amplitude over time.
• The **forcing function**: In our model, \(F_0 \sin (\omega t)\), applies an oscillating force to the system.

This system exhibits complex behavior due to the interaction between damping, stiffness, and the external force, making analytical solutions difficult but fascinating for numerical simulations.

Numerical Simulation

Numerical simulation is a powerful approach to solving differential equations that cannot be handled analytically. This involves using computing algorithms to approximate solutions over small time steps, allowing researchers and students to visualize and analyze systems without an exact solution.

For our spring-mass system, a numerical simulation provides insights into system behavior over time. By implementing the 4th order Runge-Kutta method, we resolve the motion equations into displacement and velocity profiles. The results include:

• **Displacement vs. Time**: Shows how far the mass moves from its equilibrium position.
• **Velocity vs. Time**: Indicates the rate of position change, helping identify peak speeds.
• **Velocity vs. Displacement (phase plot)**: Offers insights into system stability and energy dissipation.

These plots enable engineers and scientists to predict system behavior under various conditions, improving design and analytical capabilities in mechanical systems.