IVP’s With Step Functions — Question 5

PDF ↗

Question 5

Let Ha(t)=0H_a(t)=0 for t<at<a and Ha(t)=1H_a(t)=1 for t≥at\ge a. Use one-sided Laplace transforms. Solutions are continuous (and have continuous first derivative for second-order equations); satisfy the equation away from switches and use one-sided derivatives there. Isolated input values do not change the solution.

A resonant input is windowed without restarting its phase: y″+y=sin⁡t[Ha(t)−Hb(t)],y(0)=y′(0)=0,a=π2,b=3π2.y''+y=\sin t\,[H_a(t)-H_b(t)],\quad y(0)=y'(0)=0,\qquad a=\frac\pi 2,\quad b=\frac{3\pi}2.

Tasks

  1. Rewrite the forcing in delayed coordinates and derive its transform. Pay particular attention to the sign at the second switch.

  2. Find and invert Y(s)Y(s), deriving the required repeated-quadratic transform pair.

  3. Find the exact motion after bb and its displacement, velocity and energy at shutoff.

  4. Determine the acceleration jumps and the exact real transform domain. Explain why a zero displacement at shutoff does not mean rest.

Original worksheet page 1: question and worked solution for 4-7-005
Show solutionHide solution

Question 5 – Solution

Strategy. The two endpoints have opposite sine phases. Rewriting both phases first prevents a sign error in the delayed inverse.

Step 1: Shift the phases. For delayed time τ\tau, sin⁡(τ+a)=cos⁡τ\sin(\tau+a)=\cos\tau and sin⁡(τ+b)=−cos⁡τ\sin(\tau+b)=-\cos\tau. Therefore f=Ha(t)cos⁡(t−a)+Hb(t)cos⁡(t−b),F=s(e−as+e−bs)s2+1.f=H_a(t)\cos(t-a)+H_b(t)\cos(t-b),\qquad F=\frac{s(e^{-as}+e^{-bs})}{s^2+1}. The plus sign at bb cancels the continuing negative-phase input; it does not add a second physical forcing window.

Step 2: Derive and invert the response pair. The zero-data equation gives Y=s(e−as+e−bs)(s2+1)2,s>0.Y=\frac{s(e^{-as}+e^{-bs})}{(s^2+1)^2},\qquad s>0. Differentiating ℒ{sin⁡t}=1/(s2+1)\mathcal L\{\sin t\}=1/(s^2+1) with respect to ss gives ℒ{tsin⁡t/2}=s/(s2+1)2\mathcal L\{t\sin t/2\}=s/(s^2+1)^2. With q(τ)=τsin⁡τ/2q(\tau)=\tau\sin\tau/2, y=Ha(t)q(t−a)+Hb(t)q(t−b).\boxed{y=H_a(t) q(t-a)+H_b(t) q(t-b).} Here q(0)=q′(0)=0q(0)=q'(0)=0 and q″+q=cos⁡τq''+q=\cos\tau, directly verifying the equation, the zero initial data and both matching conditions.

Step 3: Evaluate the residual motion. For t=b+vt=b+v, v≥0v\ge 0, the two ages are v+πv+\pi and vv. Hence y(b+v)=q(v+π)+q(v)=−π2sin⁡v.y(b+v)=q(v+\pi)+q(v)=\boxed{-\tfrac\pi 2\sin v}. In particular y(b)=0y(b)=0, y′(b)=−π/2y'(b)=-\pi/2, and the post-shutoff energy is E=12(y′2+y2)=π28.\boxed{E=\tfrac 12(y'^2+y^2)=\frac{\pi^2}{8}.} The pulse ends precisely at a zero crossing with nonzero velocity, so the system continues oscillating with amplitude π/2\pi/2.

Step 4: Inspect the switches and convergence. The input jumps from 00 to 11 at aa, and from −1-1 to 00 at bb. Continuity of yy gives [y″]a=[y″]b=1\boxed{[y'']_a=[y'']_b=1}. There are no jumps in yy or y′y'. The nonzero sinusoidal tail has exact real transform domain s>0\boxed{s>0}: its undamped primitive fails to converge at zero, and exponential growth of fixed-sign tail integrals rules out negative ss.

See the diagram in the original worksheet below.

Original worksheet page 2: question and worked solution for 4-7-005

Original worksheet layout. Use Enlarge or open the PDF for a closer view.