Presuming that you are running a 2D plane strain analysis of the soil-tunnel system subjected to vertically propagated shear waves, the acceleration is imposed along the base of the numerical grid as acceleration or displacement time history.
If your contact elements act as a suspension / elastic foundation of your tunnel in regards to the soil, then you might apply the displacements to the soil nodes. Be aware of damping in this elasto-dynamic foundation. The displacements can be evaluated once you work on the periodgram of the probable earthquake in the area of your tunnel.
with periodgram I intend a statistical distribution of the energy/acceleration vs oscillation periods in the earthquake, thus depending from studies in a specific region (e.g., for my country https://ingvterremoti.wordpress.com/la-pericolosita-sismica/ or http://esse1-gis.mi.ingv.it/s1_en.php). The periodgrams give amplitude, e.g. acceleration [g] vs period [s], but not phase relationship, also depending on the long time prediction base.
If you are analyzing an actual case study, where you have the acceleration time history at surface, then you should run a deconvolution analysis, using a simplified 1D code (e.g. CyberQuake) to estimate the acceleration time history at the bedrock (e.g. base of your numerical model). You may check this indicative link for further information regarding this procedure: http://geoserver.ing.puc.cl/info/IMG/pdf/UsersGuide2-01.pdf (search for deconvolution).
If you are running a generic analysis you may select real records from recent earthequakes that are closer to the needs of your study (e.g. in terms of frequency contain, amplitude etc). In this case, you should pick acceleration time histories recorded preferably in rock or in quite stiff soil.
In both cases, you may introduce the acceleration either directly on the nodes along the base of your numerical model (assuming a rigid bedrock) or through proper dashpots at the model base to account for the finite rigidity of the bedrock (e.g. Lysmer and Kuhlemeyer analog).
I assumed my rigid bed rock by applying fixed boundary condition to the model base. and time history of earthquake was applied to the nodes of model base. my database error because of to many boundary condition to mentioned nodes. this is reasonable of course. do you have any suggestion to improve this problem?
If you want to simulate a rigid bedrock, you should fix the vertical degree of freedom of the nodes at the base of your numerical model. You will leave the horizontal degree of freedom free, and you will introduce either an acceleration or displacement time history as an acceleration or displacement boundary condition in the horizontal direction. To do so, you should set the horizontal degree of freedom equal to unity (e.g. U1 or A1) and then select an amplitude function in the Amplitude tab. This function will actually be the acceleration or the displacement time history. In any case be careful with the units you use (e.g. displacement in meters, acceleration in m/s2).