ارزیابی پایداری روش اجزای مرزی حوزه زمان در تحلیل لرزهای محیط های ناهمگن دو بعدی

روش‌های عددی یکی از ابزارهای مورد استفاده برای تحلیل پاسخ لرزه‌ای است، که در این میان روش اجزای مرزی جایگاه ویژه‌ای دارد. در این پژوهش، مطالعه جامعی بر چگونگی تاثیر ابعاد المان و طول گام زمانی بر پایداری و دقت نتایج روش اجزای مرزی حوزه زمان صورت پذیرفته است. به این منظور دو پارامتر β و l/λ که به‌طور گسترده در ادبیات فنی برای ارزیابی پایداری روش‌های عددی شناخته شده‌اند مورد استفاده قرار گرفته‌اند. سه محیط همگن، شبه‌همگن و ناهمگن با تحلیل مجموعا 280 مدل عددی مورد مطالعه قرار گرفته و مشخص شده است که پایداری و دقت نتایج در تحلیل محیط ناهمگن وابستگی اساسی به دو پارامتر مذکور داشته و نتایج قابل قبول تنها هنگامی حاصل می‌شود که موج در هر گام زمانی فاصله‌ای در حدود یک چهارم تا نصف طول المان را طی کرده  0/4 - 0/24=β و حداقل یک و نیم گره به ازای طول موج کمینه تعریف شود  (λ/l >1.0). همچنین مشخص شد در تحلیل محیط ناهمگن، ضریب β مربوط به محیطی که سرعت کمتری دارد، تعیین‌کننده پایداری و دقت نتایج تحلیل اجزای مرزی حوزه زمان خواهد بود.

evaluation of the stability of time domain boundary element method in seismic analysis of heterogeneous environments

numerical approaches are one of the best tools for seismic response analysis. in between,  the boundary element method (bem) has attracted special attention. in this paper, a comprehensive study has been performed to characterize the dependence of stability and accuracy of the time domain bem on the chosen time step duration and effective length of the elements. to this end, the two parameters β and λ/l, widely known and used in the literature for the investigation of numerical stability and accuracy, have been employed. three different environments as homogeneous, pseudo-homogeneous and non- homogeneous have been analyzed through total number of 280 numerical models. it is found that the stability and accuracy of the used algorithm is considerably influenced by the mentioned parameters, in a way that stable and accurate results will be achieved merely when the wave travels one-fourth to less than half the element size during each time step (0.24<β<0.4) and also when at least one and a half node is defined per the shortest wave-length (λ/l>1.5). it also became clear that in the modeling of non-homogeneous environments, the β value for the environment with the lowest wave velocity specifies the range of acceptable results.