IVP’s With Step Functions — Question 3

PDF ↗

Question 3

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.

An undamped oscillator receives a rectangular force of selectable duration: y″+ω2y=F[1−HT(t)],y(0)=y′(0)=0,F>0,ω>0,T>0.y''+\omega^2y=F[1-H_T(t)],\qquad y(0)=y'(0)=0, \qquad F>0,\ \omega>0,\ T>0. Use the energy E=(y′2+ω2y2)/2E=(y'^2+\omega^2y^2)/2.

Tasks

  1. Derive Y(s)Y(s) and the delayed step-response formula for yy.

  2. Find the displacement and velocity at shutoff, and the residual oscillation amplitude.

  3. Classify every duration TT for which the system remains at rest after shutoff; explain why y(T)=0y(T)=0 alone is generally insufficient.

  4. Compute the residual energy and classify the exact real transform domain, including the rest-producing durations.

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

Question 3 – Solution

Strategy. Removing the force subtracts a delayed step response. The two surviving state coordinates determine whether motion persists.

Step 1: Transform and invert. For s>0s>0, Y=F(1−e−Ts)s(s2+ω2),y=Fω2[1−cosωt−HT(t)(1−cosω(t−T))].Y=\frac{F(1-e^{-Ts})}{s(s^2+\omega^2)},\qquad \boxed{y=\frac F{\omega^2}\left[1-\cos\omega t -H_T(t)(1-\cos\omega(t-T))\right].} Before TT this has zero initial data and satisfies the forced equation. The subtracted response and its first derivative vanish at activation, so y,y′y,y' match at TT; after TT the equation is homogeneous.

Step 2: Recover the shutoff state and amplitude. The matching state is y(T)=Fω2(1−cos⁡ωT),y′(T)=Fωsin⁡ωT.y(T)=\frac F{\omega^2}(1-\cos\omega T),\qquad y'(T)=\frac F\omega\sin\omega T. For t≥Tt\ge T, subtraction and the cosine difference identity give y(t)=2Fω2sin⁡(ωT/2)sin⁡(ω(t−T/2)).y(t)=\frac{2F}{\omega^2}\sin(\omega T/2) \sin(\omega(t-T/2)). Its amplitude is therefore R=2F|sin⁡(ωT/2)|/ω2\boxed{R=2F|\sin(\omega T/2)|/\omega^2}.

Step 3: Classify complete return to rest. A homogeneous oscillator remains zero exactly when both state coordinates are zero. Here that happens precisely for T=2πnω,n=1,2,….\boxed{T=\frac{2\pi n}{\omega},\qquad n=1,2,\ldots.} In this particular pulse family, y(T)=0y(T)=0 forces cos⁡ωT=1\cos\omega T=1 and hence also y′(T)=0y'(T)=0. For a general forcing history, zero displacement alone is insufficient: a nonzero velocity launches a new oscillation.

Step 4: Compute energy and convergence. After shutoff, E=12ω2R2=2F2ω2sin⁡2(ωT/2).\boxed{E=\tfrac 12\omega^2R^2 =\frac{2F^2}{\omega^2}\sin^2(\omega T/2).} If R>0R>0, the oscillatory tail has an ordinary Laplace transform exactly for s>0s>0: at zero its primitive does not converge, and negative ss gives growing integrals over fixed sign intervals. If R=0R=0, the output has compact support, so its transform exists for every real ss. Apparent poles in the displayed formula then cancel. The graph compares F=ω=1F=\omega=1 at two durations.

See the diagram in the original worksheet below.

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

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