The equation of motion of the car is\[ my^{\prime \prime }+cy^{\prime }+ky=cy_{g}^{\prime }+ky_{g}\]
Let \(y-y_{g}=u\) which is the distance between \(m~\) and the ground. Hence the equation of motion now becomes\begin {align*} m\left ( u^{\prime \prime }+y_{g}^{\prime \prime }\right ) +c\left ( u^{\prime }+y_{g}^{\prime }\right ) +k\left ( u+y_{g}\right ) & =cy_{g}^{\prime }+ky_{g}\\ mu^{\prime \prime }+cu^{\prime }+ku & =-my_{g}^{\prime \prime } \end {align*}
| imperial | SI | |
| span length \(\lambda \) | \(70\) ft | \(70\times 0.3048=\) \(21.336\) meter |
| upward camber \(\Delta \) | \(2.5^{\prime \prime }=\) \(0.208\,33\)ft | \(2.5\times 0.0254=\) \(0.0635\) meter |
| imperial | |
| mass of car | \(\frac {1500}{32.2}=\) \(46.584\) \(\frac {lb.s^{2}}{ft}\) |
| speed of car | \(30\) mile/hr\(=44.0\) ft/sec |
| critical damping ratio is \(\zeta \) | \(0.75\) |
| spring constant \(k\) | \(2400\) lb/ft |
| natural frequency \(\omega _{n}=\sqrt {\frac {k}{m}}\) | \(\sqrt {\frac {2400}{46.\,\allowbreak 584}}=7.177\,7\) rad/sec |
| natural frequency \(f_{n}=\frac {\omega _{n}}{2\pi }\) | \(\frac {7.\,177\,7}{2\pi }=\) \(1.142\,4\) Hz |
| natural period \(T_{n}=\frac {1}{f_{n}}\) | \(\frac {1}{1.142\,4}=\) \(0.875\,35\) sec |
| natural damped frequency \(\omega _{d}=\omega _{n}\sqrt {1-\zeta ^{2}}\) | \(7.177\,7\sqrt {1-0.75^{2}}=\) \(4.747\,6\) rad/sec |
| natural damped frequency \(f_{d}=\frac {\omega _{d}}{2\pi }\) | \(\frac {4.\,747\,6}{2\pi }=0.755\,6\) Hz |
| \(T_{p}\) time to driver over one span \(=\frac {\lambda }{v}\) | \(1.591\,\)sec |
| \(T_{a}\) time to cross the bridge (duration of loading) | \(7\times 1.591=11.137\) sec |
\(a_{n}\) values found for up to \(n=10\)
Peak relative displacement of the driver
Maximum relative displacement was \(0.24\) inch and it occurred during transient phase.
Peak total displacement of the driver
\(0.165\) inch + \(2.5\) inch =\(2.665\) inch and it occurred during steady state phase at multiples of half the period \(T_{p}\) while on the bridge.
Number of \(a_{n}\) terms used
In addition to \(a_{0}\) term, the next \(5\) terms are used for a total of \(6\) terms.
The equation of the ground is shown in this diagram
Therefore, the equation of span is \[ y_{g}\left ( x\right ) =\begin {array} [c]{lll}\Delta \sin \left ( \frac {\pi }{\lambda }x\right ) & & 0\leq x\leq \lambda \end {array} \]
Hence, we convert it to be a function of time using \(x=vt\), hence\begin {align*} y_{g}\left ( t\right ) & =\begin {array} [c]{lll}\Delta \sin \left ( \frac {\pi }{\lambda }vt\right ) & & 0\leq t\leq \frac {\lambda }{v}\end {array} \\ & =\begin {array} [c]{lll}\Delta \sin \left ( \varpi t\right ) & & 0\leq t\leq T_{p}\end {array} \\ & =\begin {array} [c]{lll}\Delta \sin \left ( \frac {\pi }{T_{p}}t\right ) & & 0\leq t\leq T_{p}\end {array} \end {align*}
Where in the above \(\varpi =\frac {\pi }{T_{p}}\) is the fundamental frequency of the ground motion. Hence \(y_{g}^{\prime }\left ( t\right ) =\Delta \frac {\pi }{T_{p}}\cos \left ( \frac {\pi }{T_{p}}t\right ) \) and \[ y_{g}^{\prime \prime }\left ( t\right ) =-\Delta \left ( \frac {\pi }{T_{p}}\right ) ^{2}\sin \left ( \frac {\pi }{T_{p}}t\right ) \]
And \[ \beta =\frac {T_{n}}{T_{p}}=\frac {0.875\,35}{1.591}=0.550\,19 \]
Then load in one span \(0<t<T_{p}\) is \[ P_{a}\left ( t\right ) =m\Delta \left ( \frac {\pi }{T_{p}}\right ) ^{2}\sin \left ( \frac {\pi }{T_{p}}t\right ) \]
Let\begin {align*} P_{o} & =m\Delta \left ( \frac {\pi }{T_{p}}\right ) ^{2}=\left ( 46.584\right ) \left ( 0.208\,33\right ) \left ( \frac {\pi }{1.591}\right ) ^{2}\\ & =37.\,840\text { lb} \end {align*}
Then the load becomes\begin {equation} P_{a}\left ( t\right ) =P_{o}\sin \left ( \frac {\pi }{T_{p}}t\right ) \label {eq:A} \end {equation}
Now we need to convert Eq ?? to Fourier series1. Let \(\tilde {P}_{a}\left ( t\right ) \) be the Fourier series approximation to \(P_{a}\left ( t\right ) \), hence\begin {align*} \tilde {P}_{a}\left ( t\right ) & =a_{0}+{\displaystyle \sum \limits _{n=1}^{\infty }} a_{n}\cos \left ( n\frac {2\pi }{T_{p}}t\right ) +{\displaystyle \sum \limits _{n=1}^{\infty }} b_{n}\sin \left ( n\frac {2\pi }{T_{p}}t\right ) \\ a_{0} & =\frac {1}{T_{p}}{\displaystyle \int \limits _{0}^{T_{p}}} P_{a}\left ( t\right ) dt\\ a_{n} & =\frac {2}{T_{p}}{\displaystyle \int \limits _{0}^{T_{p}}} P_{a}\left ( t\right ) \cos \left ( n\frac {2\pi }{T_{p}}t\right ) dt\\ b_{n} & =\frac {2}{T_{p}}{\displaystyle \int \limits _{0}^{T_{p}}} P_{a}\left ( t\right ) \sin \left ( n\frac {2\pi }{T_{p}}t\right ) dt \end {align*}
Hence \begin {align*} a_{0} & =\frac {1}{T_{p}}{\displaystyle \int \limits _{0}^{T_{p}}} P_{a}\left ( t\right ) dt=\frac {P_{o}}{T_{p}}{\displaystyle \int \limits _{0}^{T_{p}}} \sin \left ( \frac {\pi }{T_{p}}t\right ) dt=\frac {P_{o}}{T_{p}}\left ( \frac {-\cos \left ( \frac {\pi }{T_{p}}t\right ) }{\frac {\pi }{T_{p}}}\right ) _{0}^{T_{p}}=-\frac {P_{o}}{\pi }\left ( \cos \left ( \pi \right ) -1\right ) \\ & =\frac {2P_{o}}{\pi }=\frac {2\left ( 37.840\right ) }{\pi }\\ & =24.090\text { lb} \end {align*}
And\begin {align*} a_{n} & =\frac {2}{T_{p}}{\displaystyle \int \limits _{0}^{T_{p}}} P_{a}\left ( t\right ) \cos \left ( n\frac {2\pi }{T_{p}}t\right ) dt\\ & =\frac {2P_{o}}{T_{p}}{\displaystyle \int \limits _{0}^{T_{p}}} \sin \left ( \frac {\pi }{T_{p}}t\right ) \cos \left ( n\frac {2\pi }{T_{p}}t\right ) dt\\ & =\frac {4P_{o}}{\pi -4n^{2}\pi }\cos \left ( n\pi \right ) ^{2} \end {align*}
But \(\cos \left ( n\pi \right ) ^{2}=1\) Hence\[ a_{n}=\frac {4P_{o}}{\pi -4n^{2}\pi }\]
and\begin {align*} b_{n} & =\frac {1}{T_{p}}{\displaystyle \int \limits _{0}^{T_{p}}} P_{a}\left ( t\right ) \sin \left ( 2\pi n\frac {t}{T_{p}}\right ) dt\\ & =\frac {2P_{o}}{T_{p}}{\displaystyle \int \limits _{0}^{T_{p}}} \sin \left ( \frac {\pi }{T_{p}}t\right ) \cos \left ( n\frac {2\pi }{T_{p}}t\right ) dt\\ & =\frac {2P_{o}}{\pi -4n^{2}\pi }\sin \left ( 2n\pi \right ) \end {align*}
But \(\sin \left ( 2n\pi \right ) =0\) for all integer \(n\), hence \(b_{n}=0.\)Therefore\begin {align*} \tilde {P}_{a}\left ( t\right ) & =a_{0}+{\displaystyle \sum \limits _{n=1}^{\infty }} a_{n}\cos \left ( 2\pi n\frac {t}{T_{p}}\right ) \\ & =a_{0}+{\displaystyle \sum \limits _{n=1}^{\infty }} \frac {4P_{o}}{\pi -4n^{2}\pi }\cos \left ( 2\pi n\frac {t}{T_{p}}\right ) \end {align*}
Using the numerical values found, we obtain\[ \tilde {P}_{a}\left ( t\right ) =24.0897+{\displaystyle \sum \limits _{n=1}^{\infty }} \frac {4\left ( 37.840\right ) }{\pi -4n^{2}\pi }\cos \left ( 2\pi n\frac {t}{1.591}\right ) \]
This plot below shows \(P_{a}\left ( t\right ) \) and its Fourier series approximation \(\tilde {P}_{a}\left ( t\right ) \) as more terms are added. This was plotted for \(t=0\cdots 5\sec \). This was done to verify that the Fourier series approximation is correct before going to the next stage of the analysis. The actual calculations used the first 6 terms of \(a_{n}\).
The equation of motion of the car is\[ my^{\prime \prime }+cy^{\prime }+ky=cy_{g}^{\prime }+ky_{g}\]
Let \(y-y_{g}=u\) which is the distance between \(m~\) and the ground. Hence the equation of motion now becomes\begin {align} m\left ( u^{\prime \prime }+y_{g}^{\prime \prime }\right ) +c\left ( u^{\prime }+y_{g}^{\prime }\right ) +k\left ( u+y_{g}\right ) & =cy_{g}^{\prime }+ky_{g}\nonumber \\ mu^{\prime \prime }+cu^{\prime }+ku & =-my_{g}^{\prime \prime } \label {eq:1.4} \end {align}
Hence Eq ?? becomes\begin {align*} mu^{\prime \prime }+cu^{\prime }+ku & =m\Delta \left ( \frac {\pi }{T_{p}}\right ) ^{2}\sin \left ( \frac {\pi }{T_{p}}t\right ) \\ & =P_{a}\left ( t\right ) \\ & ={\displaystyle \sum \limits _{n=0}^{\infty }} a_{n}\cos \left ( 2\pi n\frac {t}{T_{p}}\right ) \\ & =\operatorname {Re}\left \{ {\displaystyle \sum \limits _{n=0}^{\infty }} a_{n}e^{in\varpi t}\right \} \end {align*}
Where \(\varpi =\frac {2\pi }{T_{p}}\) is the fundamental loading harmonic. Let \(u_{ss}\left ( n\right ) =\operatorname {Re}\left \{ U_{n}e^{in\varpi t}\right \} \) be the response due to the \(n\) term in the loading function. Hence the equation of motion now becomes\begin {align*} m\operatorname {Re}\left \{ {\displaystyle \sum \limits _{n=0}^{\infty }} -n^{2}\varpi ^{2}U_{n}e^{in\varpi t}\right \} +c\operatorname {Re}\left \{ {\displaystyle \sum \limits _{n=0}^{\infty }} i\varpi nU_{n}e^{in\varpi t}\right \} +k\operatorname {Re}\left \{ {\displaystyle \sum \limits _{n=0}^{\infty }} U_{n}e^{in\varpi t}\right \} & =\operatorname {Re}\left \{ {\displaystyle \sum \limits _{n=0}^{\infty }} a_{n}e^{in\varpi t}\right \} \\ \left ( -n^{2}\varpi ^{2}m+cin\varpi +k\right ) U_{n} & =a_{n}\\ U_{n} & =\frac {a_{n}}{-n^{2}\varpi ^{2}m+cin\varpi +k}\\ & =\frac {a_{n}}{k}\frac {1}{\left ( 1-n^{2}r^{2}\right ) +2i\zeta nr} \end {align*}
Hence the transfer function is
\begin {align*} \left ( -n^{2}\varpi ^{2}m+cin\varpi +k\right ) U_{n} & =a_{n}\\ U_{n} & =\frac {a_{n}}{-n^{2}\varpi ^{2}m+cin\varpi +k}\\ & =\frac {a_{n}}{k}\frac {1}{\left ( 1-n^{2}r^{2}\right ) +2i\zeta nr} \end {align*}
Therefore, steady state response is\begin {align} y_{ss}\left ( t\right ) & =\operatorname {Re}\left \{ {\displaystyle \sum \limits _{n=0}^{\infty }} U_{n}e^{in\varpi t}\right \} \nonumber \\ & =\operatorname {Re}\left \{ {\displaystyle \sum \limits _{n=0}^{\infty }} \frac {a_{n}}{k}\frac {1}{\left ( 1-n^{2}r^{2}\right ) +2i\zeta nr}e^{in\varpi t}\right \} \nonumber \\ & =\operatorname {Re}\left \{ {\displaystyle \sum \limits _{n=0}^{\infty }} \overset {U_{n}}{\overbrace {\frac {a_{n}}{k}D\left ( \zeta ,r,n\right ) }}e^{in\varpi t}\right \} \label {eq:B} \end {align}
Where \(D\left ( \zeta ,r,n\right ) \) is the \(n^{th}\) harmonic dynamic magnification factor \[ \left \vert D\left ( \zeta ,r,n\right ) \right \vert =\sqrt {\frac {1}{\left ( 1-n^{2}r^{2}\right ) ^{2}+\left ( 2\zeta nr\right ) ^{2}}}\]
and\[ \arg D\left ( \zeta ,r,n\right ) =-\tan ^{-1}\left ( \frac {2\zeta nr}{1-n^{2}r^{2}}\right ) \]
In the above, \[ r=\frac {\varpi }{\omega _{nat}}=\frac {\frac {2\pi }{T_{p}}}{\omega _{nat}}=\frac {\frac {2\pi }{1.591\,}}{7.1777}=\fbox {$0.550\,2$}\]
This is a list of the magnitude of \(U_{n}\) for different \(n\) value to examine the contribution of each harmonic to the steady state response.
From the steady state solution \(u_{ss}\left ( t\right ) \) we found above, we now find \(u_{ss}\left ( 0\right ) \) and \(u_{ss}^{\prime }\left ( 0\right ) \) these are the initial conditions, but in opposite sign, that the transient solution have to satisfy.
From above, we found the steady state solution to be\[ y_{ss}\left ( t\right ) =\operatorname {Re}\left \{ {\displaystyle \sum \limits _{n=0}^{\infty }} U_{n}e^{in\varpi t}\right \} \]
Hence\[ y_{ss}^{\prime }\left ( t\right ) =\operatorname {Re}\left \{ {\displaystyle \sum \limits _{n=0}^{\infty }} in\varpi U_{n}e^{in\varpi t}\right \} \]
At time \(t=0\) the above becomes\begin {align*} y_{ss}\left ( 0\right ) & =\operatorname {Re}\left \{ {\displaystyle \sum \limits _{n=0}^{\infty }} U_{n}\right \} \\ y_{ss}^{\prime }\left ( 0\right ) & =\operatorname {Re}\left \{ {\displaystyle \sum \limits _{n=0}^{\infty }} in\varpi U_{n}\right \} \end {align*}
Now we need to decide on how many harmonics to use in order to determine \(y_{ss}\left ( 0\right ) \) and \(y_{ss}^{\prime }\left ( 0\right ) \). From above we see that after \(n=5\) then \(a_{n}\) became very small. Hence we will use up to \(n=5\) to find the initial conditions from the above 2 equations.\begin {align*} y_{ss}\left ( 0\right ) & =\operatorname {Re}\left \{ {\displaystyle \sum \limits _{n=0}^{5}} U_{n}\right \} =\operatorname {Re}\left \{ {\displaystyle \sum \limits _{n=0}^{5}} \frac {a_{n}}{k}\frac {1}{\left ( 1-n^{2}r^{2}\right ) +2i\zeta nr}\right \} \\ & =0.0084435\text { ft}=0.101322\text { inch} \end {align*}
and for the initial velocity we obtain\begin {align*} y_{ss}^{\prime }\left ( 0\right ) & =\operatorname {Re}\left \{ {\displaystyle \sum \limits _{n=0}^{\infty }} in\varpi U_{n}\right \} \\ & =\operatorname {Re}\left \{ {\displaystyle \sum \limits _{n=0}^{\infty }} in\varpi \frac {a_{n}}{k}\frac {1}{\left ( 1-n^{2}r^{2}\right ) +2i\zeta nr}\right \} \\ & =-0.020207\text { ft/sec}=-0.242484\text { inch/sec} \end {align*}
Now the transient solution for damped system is given by\[ u_{tr}\left ( t\right ) =e^{-\zeta \omega _{n}t}\left ( A\cos \omega _{d}t+B\sin \omega _{d}t\right ) \]
with \begin {align*} u_{tr}\left ( 0\right ) & =-0.0084435\\ u_{tr}^{\prime }\left ( 0\right ) & =+0.020207 \end {align*}
Hence \[ A=u_{tr}\left ( 0\right ) =\fbox {$-0.0084435$}\]
Taking derivative of \(u_{tr}\left ( t\right ) \) gives\[ u_{tr}^{\prime }\left ( t\right ) =\zeta \omega _{n}e^{-\zeta \omega _{n}t}\left ( A\cos \omega _{d}t+B\sin \omega _{d}t\right ) +e^{-\zeta \omega _{n}t}\left ( -A\omega _{d}\sin \omega _{d}t+B\omega _{d}\cos \omega _{d}t\right ) \]
Hence at \(t=0\) we obtain\begin {align*} u_{tr}^{\prime }\left ( 0\right ) & =\zeta \omega _{n}A+B\omega _{d}\\ B & =\frac {u_{tr}^{\prime }\left ( 0\right ) -\zeta \omega _{n}A}{\omega _{d}} \end {align*}
But \(u_{tr}^{\prime }\left ( 0\right ) =+0.020207\) ft/sec\(,A=-0.0084435\) ft\(,\zeta =0.75\),\(\omega _{d}=4.\,\allowbreak 747\,6\) rad/sec,\(\omega _{n}=7.1777\) rad/sec, hence\begin {align*} B & =\frac {0.020207-0.75\times 7.1777\times \left ( -0.0084435\right ) }{4.7476}\\ & =0.01383 \end {align*}
Therefore \begin {align*} u_{tr}\left ( t\right ) & =e^{-\zeta \omega _{n}t}\left ( -0.0084435\cos \omega _{d}t+0.01383\sin \omega _{d}t\right ) \\ & =e^{-0.75\left ( 7.1777\right ) t}\left ( -0.0084435\cos \left ( 4.7476t\right ) +0.01383\sin 4.7476t\right ) \end {align*}
This solution is now added to the steady state solution.
Zooming on the first \(1.8\) seconds shows more clearly the effect of transient solution
The transient solution effect vanishes after about \(1.5\) second.
To better see the solution obtained, we plot the relative displacement. This is the displacement felt by the passenger. First the solution is shown for the whole time to cross the bridge, then we zoom to the first 2 seconds to better see the transient solution
From the above we see that the maximum relative displacement is about \(0.24\) inch and it occurs during transient phase. During steady state, the maximum relative displacement is about \(0.165\) inch
The Fourier series approximation can also be found using the complex representation. This is the derivation using this method which gives the same result as was found earlier.\[ \tilde {P}_{a}\left ( t\right ) =\frac {1}{2}Y_{0}+\operatorname {Re}\left ( {\displaystyle \sum \limits _{n=1}^{\infty }} Y_{n}e^{in\varpi t}\right ) \]
Where \begin {align} Y_{n} & =\frac {2}{T_{p}}{\displaystyle \int \limits _{0}^{T_{p}}} P_{a}\left ( t\right ) e^{^{-in\varpi t}}dt\nonumber \\ & =\frac {2P_{o}}{T_{p}}{\displaystyle \int \limits _{0}^{T_{p}}} \sin \left ( \frac {\pi }{T_{p}}t\right ) e^{^{-in\varpi t}}dt \label {eq:1.1} \end {align}
Integration by parts, \(\int udv=uv-\int vdu\), let \(u=\sin \left ( \frac {\pi }{T_{p}}t\right ) \), hence \(du=\frac {\pi }{T_{p}}\cos \left ( \frac {\pi }{T_{p}}t\right ) \) and \(v=\frac {e^{^{-in\varpi t}}}{-in\varpi }\), therefore the above becomes\begin {align} Y_{n} & =\frac {2P_{o}}{T_{p}}\left ( \left [ \sin \left ( \frac {\pi }{T_{p}}t\right ) \frac {e^{^{-in\varpi t}}}{-in\varpi }\right ] _{0}^{T_{p}}-{\displaystyle \int \limits _{0}^{T_{p}}} \frac {\pi }{T_{p}}\cos \left ( \frac {\pi }{T_{p}}t\right ) \frac {e^{^{-in\varpi t}}}{-in\varpi }dt\right ) \nonumber \\ & =\frac {2P_{o}}{T_{p}}\left ( \sin \left ( \frac {\pi }{T_{p}}T_{p}\right ) \frac {ie^{^{-in\varpi T_{p}}}}{n\varpi }-\frac {i}{2n}{\displaystyle \int \limits _{0}^{T_{p}}} \cos \left ( \frac {\pi }{T_{p}}t\right ) e^{^{-in\varpi t}}dt\right ) \nonumber \\ & =-\frac {iP_{o}}{nT_{p}}{\displaystyle \int \limits _{0}^{T_{p}}} \cos \left ( \frac {\pi }{T_{p}}t\right ) e^{^{-in\varpi t}}dt \label {eq:1.2} \end {align}
Now integrate by parts again where now \(\int udv=uv-\int vdu\), let \(u=\cos \left ( \frac {\pi }{T_{p}}t\right ) \), hence \(du=-\frac {\pi }{T_{p}}\sin \left ( \frac {\pi }{T_{p}}t\right ) \) and \(v=\frac {e^{^{-in\varpi t}}}{-in\varpi }\), therefore Eq ?? becomes\begin {align} Y_{n} & =-\frac {iP_{o}}{nT_{p}}\left [ \left ( \cos \left ( \frac {\pi }{T_{p}}t\right ) \frac {e^{^{-in\varpi t}}}{-in\varpi }\right ) _{0}^{T_{p}}-{\displaystyle \int \limits _{0}^{T_{p}}} -\frac {\pi }{T_{p}}\sin \left ( \frac {\pi }{T_{p}}t\right ) \frac {e^{^{-in\varpi t}}}{-in\varpi }dt\right ] \nonumber \\ & =-\frac {iP_{o}}{nT_{p}}\left [ \left ( \cos \left ( \frac {\pi }{T_{p}}T_{p}\right ) \frac {e^{^{-in\frac {2\pi }{Tp}T_{p}}}}{-in\varpi }-\frac {1}{-in\varpi }\right ) +\frac {i}{n2}{\displaystyle \int \limits _{0}^{T_{p}}} \sin \left ( \frac {\pi }{T_{p}}t\right ) e^{^{-in\varpi t}}dt\right ] \nonumber \\ & =-\frac {iP_{o}}{nT_{p}}\left [ \left ( -i\frac {e^{^{-in2\pi }}}{n\varpi }-\frac {i}{n\varpi }\right ) +\frac {i}{n2}{\displaystyle \int \limits _{0}^{T_{p}}} \sin \left ( \frac {\pi }{T_{p}}t\right ) e^{^{-in\varpi t}}dt\right ] \nonumber \\ & =-\frac {P_{o}}{nT_{p}}\left ( \frac {e^{^{-in2\pi }}+1}{n\varpi }\right ) +\frac {\Delta }{2n^{2}T_{p}}{\displaystyle \int \limits _{0}^{T_{p}}} \sin \left ( \frac {\pi }{T_{p}}t\right ) e^{^{-in\varpi t}}dt \label {eq:1.3} \end {align}
Now we see that the term \({\displaystyle \int \limits _{0}^{T_{p}}} \sin \left ( \frac {\pi }{T_{p}}t\right ) e^{^{-in\varpi t}}dt\) has repeated again. This term is the same as what we started with in Eq ?? , therefore, we write\[{\displaystyle \int \limits _{0}^{T_{p}}} \sin \left ( \frac {\pi }{T_{p}}t\right ) e^{^{-in\varpi t}}dt=\frac {T_{p}}{2\Delta }Y_{n}\]
and replace this term back into Eq ?? , hence it becomes\begin {align*} Y_{n} & =-\frac {P_{o}}{nT_{p}}\left ( \frac {e^{^{-in2\pi }}+1}{n\varpi }\right ) +\frac {\Delta }{2n^{2}T_{p}}\frac {T_{p}}{2\Delta }Y_{n}\\ & =-\frac {P_{o}}{nT_{p}}\left ( \frac {e^{^{-in2\pi }}+1}{n\varpi }\right ) +\frac {1}{2^{2}n^{2}}Y_{n}\\ Y_{n}-\frac {1}{2^{2}n^{2}}Y_{n} & =-\frac {P_{o}}{nT_{p}}\left ( \frac {e^{^{-in2\pi }}+1}{n\varpi }\right ) \\ Y_{n}\left ( 1-\frac {1}{\left ( 2n\right ) ^{2}}\right ) & =-\frac {P_{o}}{nT_{p}}\left ( \frac {e^{^{-in2\pi }}+1}{n\varpi }\right ) \\ Y_{n} & =-\frac {2P_{o}\left ( e^{^{-in2\pi }}+1\right ) }{\pi \left ( \left ( 2n\right ) ^{2}-1\right ) }=\frac {2P_{o}\left ( e^{^{-in2\pi }}+1\right ) }{\pi -\pi \left ( 2n\right ) ^{2}}\\ & =\frac {4P_{o}}{\pi \left ( 1-4n^{2}\right ) } \end {align*}
And \begin {align*} Y_{0} & =\frac {2}{T_{p}}{\displaystyle \int \limits _{0}^{T_{p}}} P_{o}\sin \left ( \frac {\pi }{T_{p}}t\right ) dt=\frac {2P_{o}}{T_{p}}{\displaystyle \int \limits _{0}^{T_{p}}} \sin \left ( \frac {\pi }{T_{p}}t\right ) dt=\frac {2P_{o}}{T_{p}}\left ( \frac {-\cos \left ( \frac {\pi }{T_{p}}t\right ) }{\frac {\pi }{T_{p}}}\right ) _{0}^{T_{p}}=-\frac {2P_{o}}{\pi }\left ( \cos \left ( \frac {\pi }{T_{p}}T_{p}\right ) -1\right ) \\ & =-\frac {2P_{o}}{\pi }\left ( -1-1\right ) \\ & =\frac {4P_{o}}{\pi } \end {align*}
Therefore, the Fourier series approximation for ground motion is now\begin {align*} \tilde {P}_{a}\left ( t\right ) & =\frac {1}{2}Y_{0}+\operatorname {Re}\left ( {\displaystyle \sum \limits _{n=1}^{\infty }} Y_{n}e^{in\varpi t}\right ) \\ & =\frac {4P_{o}}{2\pi }+\operatorname {Re}\left ( {\displaystyle \sum \limits _{n=1}^{\infty }} \frac {4P_{o}}{\pi \left ( 1-4n^{2}\right ) }e^{in\varpi t}\right ) \\ & =\frac {2P_{o}}{\pi }+\operatorname {Re}\left ( {\displaystyle \sum \limits _{n=1}^{\infty }} \frac {4P_{o}}{\pi \left ( 1-4n^{2}\right ) }e^{in\varpi t}\right ) \end {align*}
We see that we obtained the same result using the classical Fourier series form.
1The Fourier series can also be found using complex form. This was done in the appendix.