Joe,
To answer your questions:
1) MS is just a method by which a MDOF system is uncoupled into individual modes, which allows us to calculate the response for each mode individually, and then combines the individual response back into the overall behavior. MS does not allow for a solution on its own, it must be used in conjunction with other methods such as RSA, Duhamel's integral, or time-step integration. MS is just a tool that is used to simplify analysis for MDOF systems. RSA is a method that gives you the maximum response for a given mode. If you have a SDOF system RSA can be used on its own to determine the response (base shear, acceleration...). However if you have a MDOF system then you need MS to analyze the contribution of each mode, if you want an accurate answer.
2) A single response spectra is derived from a single input acceleration record, such as a recorded earthquake. The frequency content of the input acceleration is analyzed to give maximum responses for varying natural frequencies of structures. A single response spectra is not as smooth as you see in the codes, since a single input acceleration has varying frequency content. The way a design spectrum is derived is by combining many different response spectrums for different input accelerations (that are characteristic for that location), in order to envelope the potential responses. It is then smoothed to give you a nice equation that you see in the code. Response spectrum analysis can use either the design spectrum or the exact response spectrum, but in most cases the design spectrum is used due to uncertainty and conservatism.
3) Here is brief walk-through of RSA and MS:
MS:
- define mass [M] and stiffness [K] matrices
- use eigenvalue analysis to derive natural frequencies and mode shapes: det([K] - w2[M])[A] = 0
- calculate generalized masses: Mi = {Ai}T[M]{Ai]}
- calculate generalized loads: Pi = -{Ai}T[M]{r}a(t)
- Use RSA/Duhamels Integral/Time integration.... to calculate the response of each mode {Ui(t)}
- Combine responses into overall behavior: {x(t)} = [A]{Ui(t)}
- Calculate equivalent static loads: {Q(t)} = [k]{x(t)}
- Base shear: V(t) = (Q(t)}
The parts in bold are part of MS, everything else is algebra, dynamics.
RSA:
- define mass [M] and stiffness [K] matrices
- use eigenvalue analysis to derive natural frequencies and mode shapes: det([K] - w^2[M])[A] = 0
- calculate generalized masses: Mi = {Ai}T[M]{Ai]}
- calculate modal participation factors: ai={Ai}T[M]{r}/Mi
- Given the natural frequencies of each mode get the maximum response acceleration from a relevant response spectra - Sa
- Base shear: Vmax = Miai2Sa
Part in bold is RSA