پژوهشگاه بین المللی زلزله شناسی و مهندسی زلزلهفصلنامه علوم و مهندسی زلزله2476-60977420201221Evaluation of Dynamic and Kinematic Parameters of November 12, 2017 Sarpol-e Zahab Earthquake in Zagros Regionبررسی ویژگی های دینامیکی و سینماتیکی زمین لرزه 12 نوامبر 2017 سرپل ذهاب در استان کرمانشاه114243550FAناهیدخودی آغمیونیپژوهشکده زلزله شناسی، پژوهشگاه بینالمللی زلزلهشناسی و مهندسی زلزله، تهران، ایرانمهردادمصطفی زادهپژوهشکده زلزله شناسی، پژوهشگاه بینالمللی زلزلهشناسی و مهندسی زلزله، تهران، ایرانJournal Article20181009The Zagros fold Thrust belt is approximately 1,500 to 1,600 km long and 200 to 300 km wide, formed from western Iran to the Strait of Hormuz under the influence of the convergence of the Arabian and Eurasian plateaus. The faults in this belt are mainly parallel to the folds and are of the thrust type with slopes to the north to the northeast. The Zagros zone is divided into five units based on morpho-tectonic, seismic, and structural data, and it is separated by deep faults called the High Zagros Fault (HZF), Mountain Front Fault (MFF), Zagros Foredeep Fault (ZFF), Main Zagros Reverse Fault (MZRF) and Kazerun Fault (KF). The earthquake of November 12, 2017, which occurred in the Zagros region near to Iran-Iraq border (MW = 7.3), is the only earthquake with a magnitude of more than 7 recorded in the northwestern part of the Zagros Thrust Zone (ZTZ). The largest earthquake in the region near the source of the recent earthquake was the January 11, 1967 earthquake (M = 6.1), located about 100 km south of the epicenter of the 2017 earthquake. The closest city to the epicenter of the recent earthquake is the city of Azgeleh with a distance of approximately 5 km. The depth of the earthquake was also reported to be about 11 km. In this study, using seismic data recorded in the Global Digital Seismic Network (GDSN), the seismic moment of this earthquake was calculated using the Kelis-Borok method and its average value obtained was 1.45E+27 dyne-cm. According to the source time function of this earthquake and obtaining the peak value of the rise time, duration of the most fracturing and fall time, the total rupture time is calculated that is 4.998 seconds.زمینلرزه 12 نوامبر 2017 که در مرز شمالی عراق و غرب ایران (MW = 7.3) به وقوع پیوست، تنها زمینلرزهی دستگاهی است که با بزرگای بیش از 7 در ناحیه شمال غربی زون تراستی زاگرس (ZTZ) ثبت گردیده است. نزدیکترین شهر به کانون زمینلرزه، شهر ازگله با فاصله تقریبی 5 کیلومتر میباشد. همچنین عمق زمینلرزه حدود 11 کیلومتر گزارش شده است. این زمینلرزه موجب کشته شدن صدها نفر و مجروح شدن هزاران نفر گردیده و خسارات اقتصادی زیادی را در استان کرمانشاه به همراه داشته است. در این مطالعه با استفاده از اطلاعات لرزهای دورلرز ثبت شده در شبکه لرزهای جهانی (GDSN) ، گشتاور لرزهای این زمینلرزه با استفاده از روش کیلیس- بورک (1960) (Keilis-Borok) محاسبه و مقدار متوسط آن برابر 1.45E+27 دین- سانتیمتر تعیین گردیده است. با توجه به تابع زمانی چشمه و به دست آوردن اوج تابع زمان چشمه (δt1= 1.935) و گستره زمانی که بیشترین شکستگی طی آن ایجاد میشود (𝛿𝑡2= 3.225) و زمان شروع افت شکستگی بر روی سطح گسل (𝛿𝑡3= 1.612)، مدتزمان گسیختگی این زمینلرزه4.998 ثانیه تعیین گردیده است.https://www.bese.ir/article_243550_7f6846a98533d66b9c2f195336a3749b.pdfپژوهشگاه بین المللی زلزله شناسی و مهندسی زلزلهفصلنامه علوم و مهندسی زلزله2476-60977420201221Extracting Rayleigh Waves of Iran Earthquakes from Three-Component Signals in Time-Frequency Domainاستخراج امواج ریلی زلزله های ایران با استفاده از سه مولفه شتاب نگاشت در حوزه زمان-فرکانس1532243551FAپویانقشینگروه مهندسی عمران، دانشکده فنی و مهندسی، دانشگاه ارومیه، ارومیه، ایرانهادیبهادریگروه مهندسی عمران، دانشکده فنی و مهندسی، دانشگاه ارومیه، ارومیه، ایرانعباساسلامی حقیقتگروه مهندسی عمران، دانشکده فنی و مهندسی، دانشگاه ارومیه، ارومیه، ایرانJournal Article20190716Rayleigh waves can be more destructive in comparison with other types of waves during earthquakes. Identification and extraction of Rayleigh waves from seismic records can be more exactly performed in time-frequency domain. In the present study, Continuous Wavelet Transform and Stockwell Transform have been used to transfer three-components of seismic signals to time-frequency domain. Rayleigh waves can be extracted based on elliptical characteristics of these waves. On the other hand, Retrograde and Prograde Rayleigh waves and their propagation azimuth can be separately extracted with phase transferring of vertical component as well as normalized inner product technique. In the present study, the mentioned two algorithms have been programmed. Accuracy of the methods has been verified with using synthetic and Chi-Chi earthquake data. Finally, Rayleigh waves and their seismic parameters have been extracted from Tabas, Bam and Manjil earthquake records. Rayleigh waves of the earthquakes have been analyzed and compared with using different seismic parameters. According to the results, Retrograde Rayleigh waves are more predominant in comparison with Prograde Rayleigh waves in earthquake signals.<br />Introduction<br />Extraction of Rayleigh waves is one of the important tasks in seismology and applied geophysics due to their destructive effects. A lot of study has been conducted in order to develop an effective method for the extraction of Rayleigh waves. Most of these methods have been presented based on the ellipticity characteristics of Rayleigh waves. Each of them has their advantages and disadvantages. In this study, two of them have been investigated, programmed and compared with each other.<br />Methodology and Approaches<br />Rayleigh waves have been extracted from three-component signals in time-frequency domain in both of the investigated methods. In algorithm 1, continuous wavelet transform has been used to transfer three components of signals from time to time-frequency domain. However, Stockwell transform has been used to perform this approach in algorithm 2. Stockwell transform retain the absolute phase of each localized frequency component, the characteristic of which is required in algorithm 2 for time-frequency transform. On the other hand, Rayleigh waves have been extracted based on the instantaneous reciprocal ellipticity and phase difference in algorithm 1 and 2 respectively. The most important advantages of algorithm 2 is that retrograde and prograde Rayleigh waves and their propagation azimuth can be separately extracted while it is not possible in algorithm 2.<br />In this study, both of these algorithms have been clearly illustrated by flowcharts and programmed in Matlab. Synthetic signals and Chi-Chi earthquake signals have been used for verification. After verification, results of both algorithms have been compared, and finally they have been employed in extracting Rayleigh waves of Tabas, Bam and Manjil earthquakes. Finally, maximum displacements of Rayleigh waves and total earthquake signals have been investigated and compared in these three earthquakes. In addition, energy discharge of extracted Rayleigh waves during the earthquakes have been investigated. <br />Results and Conclusions<br />Results of the study showed that algorithm 1 can be effectively used if there is not required to extract retrograde and prograde Rayleigh waves and their propagation azimuth separately because the algorithm is simpler and faster in comparison with algorithm 2. However, algorithm 2 should be employed on the condition that the extraction of retrograde and prograde Rayleigh waves and their propagation azimuth have been required. In addition, it has been concluded that Rayleigh waves have more energy in earthquake signals compared to other types of waves. Consequently, elliptical polarization is more than linear polarization in earthquakes. In comparison of retrograde and prograde Rayleigh waves, it has been concluded that retrograde Rayleigh waves have more energy and therefore they can be more destructive. Finally, it has been resulted that ratio of Rayleigh waves are different in earthquakes. Induced displacement by Rayleigh waves in Manjil earthquake was more than Bam and Tabas earthquakes.<br /><br />امواج ریلی نسبت به سایر امواج، اثرات مخربتری بر روی سازههای ساخت بشر دارند. شناسایی و استخراج امواج ریلی از شتابنگاشتها در حوزه زمان- فرکانس از دقت بالایی برخوردار است. در تحقیق حاضر، از تبدیل موجک پیوسته و تبدیل استوکول برای انتقال سه مؤلفه شتابنگاشت هر ایستگاه زلزله به حوزه زمان- فرکانس و استخراج امواج ریلی استفاده شده است. با استفاده از مشخصات حرکت بیضوی، امواج ریلی از شتابنگاشتها در حوزه زمان- فرکانس قابل استخراج میباشند. همچنین با انتقال فاز حرکت قائم و با استفاده از ضرب داخلی نرمالایز شده، میتوان امواج ریلی پسرونده، پیشرونده و نیز زاویه انتشار این امواج را در حوزه زمان- فرکانس بهطور مجزا استخراج نمود و درنهایت پاسخ را در فضای زمان ارائه داد. در مطالعه حاضر، الگوریتمهای مذکور در زبان برنامهنویسی متلب کدنویسی شده است. صحت تحلیلهای صورت گرفته، با استفاده از امواج ریلی استخراج شده از سیگنالهای مصنوعی و نیز دادههای زلزله چیچی، ارزیابی شده است. در نهایت زلزلههای ایران شامل طبس، بم و منجیل مورد بررسی قرار گرفته و امواج ریلی آنها استخراج و پارامترهای لرزهای آنها محاسبه شده است. نتایج حاصله نشان میدهد که در شتابنگاشتها سهم امواج ریلی پسرونده بیشتر از سهم امواج ریلی پیشرونده و سایر امواج است.https://www.bese.ir/article_243551_0e93126f716d5d138775b090ec680d68.pdfپژوهشگاه بین المللی زلزله شناسی و مهندسی زلزلهفصلنامه علوم و مهندسی زلزله2476-60977420201221Predicting the Seismic Response of Embankments and Earth Dams with Artificial Neural Networkپیشبینی پاسخ لرزهای خاکریزها و سدهای خاکی با شبکه عصبی مصنوعی3350243552FAامینرضایاندانشکده عمران، واحد علوم و تحقیقات تهران، دانشگاه آزاد اسلامی، تهران، ایرانمحمدداودیپژوهشکده مهندسی ژئوتکنیک، پژوهشگاه بینالمللی زلزلهشناسی و مهندسی زلزله، تهران، ایرانمحمد کاظمجعفریپژوهشکده مدیریت خطرپذیری و بحران، پژوهشگاه بینالمللی زلزلهشناسی و مهندسی زلزله، تهران، ایرانJournal Article20181106The artificial neural network (ANN), one of the most powerful tools of artificial intelligence, was used to control the seismic responses of the cross section, and to predict earth dam seismic responses rather than utilizing dynamic time-based analyses. In terms of the artificial neural network, apart from the aforementioned ANN-based applications in various engineering problems, an increasing number of articles have been published over the last decade where the efficient implementation of ANNs in geotechnical earthquake engineering is presented. Most of these studies focus on liquefaction potential under seismic excitations, which is an extremely computationally intensive issue and therefore suitable for ANNs. In some of the studies in this field, the applicability of ANNs in soil dynamic analysis was examined. In addition to the generation of spectrum compatible accelerograms, prediction of earthquake parameters, wave propagation approximations, estimation of peak ground acceleration using microtremors, and the detection of earthquake electric field patterns. Not much has been done into predicting the seismic responses of embankments and earth dams using artificial neural networks tools. The first study on the domain was carried out by Tsompanakis et al. [1] in 2009 in which using the artificial neural network, a seismic response (maximum horizontal acceleration) of homogeneous and symmetric embankment in different locations was predicted. Although being a case study on a simple embankment, earthquake records (accelerograms) were used, which contained a variety of frequency concepts, simplicity or complexity, and various PGA levels, in addition to a general neural network that could be used in all PGA levels the seismic response was not accurately predicted. The objective of the present study was to investigate the ability of the artificial neural network to reach the dynamic response of the embankment and earth dams under seismic loads and different levels of intensity. In other words, it was planned to reduce the heavy cost of calculating the applicable issue in seismic geotechnical engineering. To fulfill the aim, a comprehensive study was carried out and various parameters were evaluated. According to the results, the artificial neural network could reach a good approximation of the seismic response of embankments and earth dams. Although comparing the results of artificial neural networks segmented from generic artificial neural network showed the accuracy of the estimated artificial neural network structured for each PGA level, especially in the linear region, the study was successful in eliminating the defects of previous research, making a perfect arrangement for a generic artificial neural network to be used for this type of neural network to achieve a favorable estimation of all PGA levels. Finally, it can be concluded that considering the complexities of the problem and the performance of the artificial neural networks in solving them, as well as the history of successful performance reports on this tool in other similar applications, the application of Artificial Neural Networks proved useful for estimating and predicting seismic responses of embankments and earth dams. Using artificial neural networks to meet this purpose would reduce the cost of calculating the seismic response assessment of embankments and earth dams. Therefore, applying this tool in the complicated field of seismic geotechnical engineering may be highly promising. Overall, future developments in this research field and promotions in artificial neural network training in other situations by changing the dimensions and properties of materials and load conditions would realize this idea and create a more general neural network for this application. <br />Reference<br />1. Tsompanakis, Y., Lagaros, N.D., and Stavroulakis, G.E. (2008) Soft computing techniques in parameter identification and probabilistic seismic analysis of structures. Advances in Engineering Software, 39(7), 612–624.مهندسی ژئوتکنیک لرزهای غالباً بهعنوان یک شاخه علمی غیردقیق به حساب میآید، زیرا در طی مراحل طراحی سازههای ژئوتکنیکی با عدم قطعیتهای اجتنابناپذیر و سادهسازیهایی مواجه هستیم که ناچاریم آنها را بپذیریم. بنابراین پیشبینیهای نسبتاً دقیق با استفاده از تکنیکهای محاسبات نرم (SC) پیشرفته میتواند دشواریهای کار در روشهای حل متعارف را هموار کند. شبکههای عصبی مصنوعی (ANN) یکی از مشهورترین تکنیکهای محاسبات نرم است که در زمینههای مختلف علمی و فناوری استفاده میشود. کاربرد این تکنیک در زمینههای کاربردی مهندسی زلزله در سازهها نیز همزمان با توسعه آن در سایر زمینههای علمی افزایش یافته است. این مقاله روی کاربرد شبکه عصبی مصنوعی بر روی شبیهسازی پاسخ لرزهای خاکریزها و سدهای خاکی متمرکز شده است. پاسخ دینامیکی خاکریزها و سدهای خاکی با استفاده از روش اجزای محدود و با استفاده از مدل معادل خطی ارزیابی شده است. در پژوهش حاضر، این فرآیند نسبتاً زمانبر با پیشبینیهای سریع شبکههای عصبی مصنوعی که بهطور صحیح آموزش دیده است جایگزین شده است. در اینجا ورودیهای مدل شبکه عصبی مصنوعی پارامترهای لرزهای تحریکهای زلزله وارد بر خاکریز یا سد خاکی بوده و خروجی آن در این پژوهش شتاب افقی حداکثر تاج سد است. بررسیهای انجام شده در این پژوهش نشان داد که مدل شبکه عصبی مصنوعی کلی ارائه شده جهت پیشبینی پاسخ لرزهای خاکریزها و سدهای خاکی میتواند نسبت به مدل جزئی پیشنهاد شده در پژوهشهای پیشین کاربردیتر باشد.https://www.bese.ir/article_243552_d317e55f830af13a116768848a2d73f7.pdfپژوهشگاه بین المللی زلزله شناسی و مهندسی زلزلهفصلنامه علوم و مهندسی زلزله2476-60977420201221Numerical Modeling of the Interaction of Micropile with Reverse Faultingمدلسازی عددی اندرکنش ریزشمع ها با گسلش سطحی معکوس5163243553FAمجتبیعلیزادهگروه مهندسی عمران، دانشکده فنی و مهندسی، دانشگاه قم، قم، ایرانمهدیخداپرستگروه مهندسی عمران، دانشکده فنی و مهندسی، دانشگاه قم، قم، ایرانعلی محمدرجبیدانشکده زمین شناسی، پردیس علوم، دانشگاه تهران، تهران، ایرانJournal Article20200214وقوع گسیختگیهای ناشی از گسلش سطحی زلزلههای اخیر، ضرورت بررسی اندرکنش گسل- سازه - فونداسیون را در طراحی سازههای مهندسی و شریانهای حیاتی واقع در حریم گسل، بیشازپیش میطلبد. در این مقاله، اندرکنش گسلش سطحی حاصل از عملکرد یک گسل معکوس و مجموعه پی شامل؛ فونداسیون سطحی و ریزشمع ها و همچنین اثر آن بر جابهجاییهای افقی و قائم فونداسیون، با استفاده از مدلسازی عددی مورد بررسی قرارگرفته است. مهمترین فاکتور مورد استفاده در این مطالعه، نسبت فاصله محل برخورد گسلش به زیر پی به عرض پی (S/B) میباشد. نتایج نشان میدهد جابهجایی افقی فونداسیون به حضور و موقعیت قرارگیری ریزشمعها حساسیت اندکی نشان میدهد درحالیکه جابهجایی قائم و چرخش فونداسیون با وجود ریز شمعها به موقعیت قرارگیری مجموعه پی نسبت به گسلش حساس بوده و وجود ریز شمعها همواره باعث کاهش مقادیر جابهجایی قائم فونداسیون میگردد. نتایج همچنین نشان میدهد بهترین حالت برای کنترل چرخش فونداسیون، حالت S/Bhttps://www.bese.ir/article_243553_b89b7576b0cc175679789aadec0e47d5.pdfپژوهشگاه بین المللی زلزله شناسی و مهندسی زلزلهفصلنامه علوم و مهندسی زلزله2476-60977420201221Investigating the Seismic Behaviour of the RC Wall-Framed System Considering Soil-Structure Interactionبررسی رفتار لرزهای سیستم دوگانه دیوار برشی– قاب خمشی با در نظر گرفتن اندرکنش بین خاک و سازه6581243554FAمحمودقطرنجیگروه مهندسی عمران، دانشکده فنی و مهندسی، دانشگاه بینالمللی امام خمینی (ره)، قزوین، ایرانمحمد مهدیمعمارپورگروه مهندسی عمران، دانشکده فنی و مهندسی، دانشگاه بینالمللی امام خمینی (ره)، قزوین، ایرانمنصوریخچالیانگروه مهندسی عمران، دانشکده فنی و مهندسی، دانشگاه بینالمللی امام خمینی (ره)، قزوین، ایران0000-0003-3077-5213Journal Article20171107In most cases of design or assessment of the buildings, the fixed foundation assumption is used and the effect of soil-structure interaction (SSI) is neglected. In fact, considering the SSI effects in analysis leads to substantial changes in the dynamic properties and the dynamic responses of the structure. On the other hand, the shear wall-frame lateral load resisting system is one of the most popular systems for resisting the earthquake and wind loads because of its ability to provide high energy dissipation, high lateral stiffness and high lateral capacity in comparison with shear wall or bare frame only. In this paper, a comprehensive study of the seismic behavior of the RC wall-frame system with and without SSI effects is performed. For this aim, three 5-, 10- and 15-story buildings were simulated in OpenSees software so that the wall-frame dual system for each of them is located at the perimeter of the structure. The shear wall modeling was performed using the shear flexure interaction (SFI) model, which was developed by (Kolozvari et al., 2015). For nonlinear modeling of beams and columns, the beamWithHinges and nonlinearBeamColumn elements<br />were used, respectively. Also, by using the beam on nonlinear Winkler foundation (BNWF) theory the SSI phenomena was simulated. Modal and pushover analyses were carried out in two cases, namely, fixed and flexible bases. The lateral load capacity, structure period, base shear contribution ratios of the shear wall and moment frame and the story drifts were studied in details. Also, for different performance levels, the base shear, story shear force and the development of the plastic hinges were investigated. The results of the modal analyses show that the SSI effects increase the first mode period of the structure and the percentages of increase in the first mode period for the 5-, 10- and 15-story models are 75.2, 23.8 and 13.2%, respectively. In addition, considering the SSI effects in eigenvalue analysis leads to change in the first mode shape of the 5-story model, so that it will be close to linear case. Also, the SSI effects lead to decrease in the base shear capacity of the 5-story model and by increasing the height of structure its effects on the base shear capacity are negligible. The results indicate that the foundation of the shear wall under SSI effects has a rotational movement that decreases the shear wall contribution in bearing the lateral loads and increases the frame base shear contribution at the linear stage. Besides, by increasing the lateral load level, the shear wall base shear contribution ratio increases and more nonlinearity will take place at this stage. In case of flexible base, the story drifts increase, especially in the lower stories and at the high damage stage the base shear contribution ratios of the shear wall and frame will be equal. Also, the SSI effects lead to increase in the seismic demands of the structural components in the lower stories and decrease them at the higher stories.در اکثر موارد در طراحی و ارزیابی سازهها، تکیهگاه سازه بهصورت گیردار فرض شده و اثرات اندرکنش خاک و سازه لحاظ نمیشود؛ درصورتیکه در نظر گرفتن اثرات اندرکنش خاک و سازه ویژگیهای دینامیکی سیستم را عوض میکند. در این مطالعه رفتار لرزهای سیستم دوگانه دیوار برشی- قاب خمشی با در نظر گرفتن اثرات اندرکنش خاک و سازه بررسی شده است. بدینمنظور سه سازه با تعداد طبقات 5، 10 و 15 در محیط نرمافزار OpenSees مدلسازی شده و تحلیلهای استاتیکی غیرخطی در دو حالت پایه گیردار و پایه انعطافپذیر انجام شده است. نتایج نشان میدهد که در نظر گرفتن SSI باعث افزایش زمان تناوب سازه و کاهش برش طبقات میشود. همچنین با در نظر گرفتن اثرات اندرکنش خاک و سازه ظرفیت سازه کوتاه کاهش پیدا میکند درحالیکه با افزایش ارتفاع سازه مقدار ظرفیت آن تفاوت چندانی ندارد. از سوی دیگر در اثر اندرکنش خاک و سازه فونداسیون دیوار برشی حرکت دورانی دارد و این باعث میشود که دیوار برشی در آغاز بارگذاری مقدار کمی از بار جانبی را تحمل نماید و هرچقدر بار جانبی بیشتر میشود سهم دیوار برشی در تحمل بار جانبی بیشتر شود.https://www.bese.ir/article_243554_d53f695083aa2b5b4170baa69b2107c2.pdfپژوهشگاه بین المللی زلزله شناسی و مهندسی زلزلهفصلنامه علوم و مهندسی زلزله2476-60977420201221Damage Detection by Energy and Displacement Damage Index on the Analytical Phase of ASCE Benchmarkتشخیص آسیب با استفاده از شاخص آسیب انرژی و جابهجایی در فاز تحلیلی مدل ASCE8399243558FAمحمدجوادخسرویانیگروه مهندسی عمران، واحد قزوین، دانشگاه آزاد اسلامی، قزوین، ایرانامیدبهارپژوهشکده مهندسی سازه، پژوهشگاه بینالمللی زلزلهشناسی و مهندسی زلزله، تهران، ایران0000-0002-9655-4512سید هومنقاسمیگروه مهندسی عمران، واحد قزوین، دانشگاه آزاد اسلامی، قزوین، ایرانJournal Article20181124Identification of the modal parameters of the damaged structure by signal processing of vibration based on changes in dynamic properties such as frequency, stiffness, and mode shape of structures. Some of these approaches fail when applied to civill engineering structures, the main reasons are the low sensitivity of the structural response to the damage location, or the low accuracy of structural response obtained by installed sensors. However, due to the rigorous evaluation and low cost of signal processing methods, this method has made a great progress. Signal processing methods have been extensively employed to examine the measured system responses and determine system variations. These methods include Fourier analysis, wavelet transform and Hilbert-Huang transform.<br />During the last decades, the number of vibration-based damage detection methods has been greatly developed and has influenced much of the research. The purpose of these methods is to determine the resulting changes in the modal characteristics of the structure including its natural frequency, mode shape, and damping ratio. For example, the basis of the Fourier transform method is to determine the structural modal parameters from the random vibration data in the frequency range. However, time-frequency analysis is introduced to overcome the limitations of the Fourier method, the most important of which is not providing a frequency-time range of a signal. The first case of frequency-time analysis was the short-time Fourier transform method based on the Fourier transform of the data divided by the time window function. According to this method, the interaction between time and frequency is difficult due to the existence of the time window function. If the windows are smaller in the time segment, its accuracy increases, and in the frequency domain it becomes less accurate.<br />Several damage detection methods have been proposed based on the vibration signal of structures. In most of them, a damage index has been described as the difference between damage and undamaged structure. This paper intends to propose a damage detection method based on the amplitude coefficient correlation of damaged and undamaged responses of structures, while a signal decomposes to IMFs and the changes appear in the first IMF. Therefore, every change on the original signal can be revealed on IMFs, since the original signal depends on IMFs. Also, these changes have an effect on the analytical signal and the Hilbert transform. The instantaneous frequency is measured joints on the structure is calculated by the Hilbert transform of the first IMF of response. Then, by introducing the instantaneous frequency energy (EDI), the location of damages are detected. To assess the feasibility and reliability of the proposed method, the ASCE benchmark problem has been used. To consider the robustness of the proposed method, contamination of signals during the data acquisition process is investigated. The ASCE benchmark study is carried out by the International Association (IASC) ASCE Structural Health Monitoring Task Group. The dynamic responses of the structure have been obtained by numerical analysis under random vibration loading. To evaluate the HHT method, it is required to attain the damaged responses of the ASCE benchmark. The first five damage patterns of the ASCE benchmark building is used. Then the damage on the structure is detected with a comparison of damage and undamaged dynamic responses. According to the measured noise levels, dynamic responses to noise values have been contaminated and the results have been evaluated. According to the results, this method can trace the location of the damage by the energy of instantaneous frequency. <br />Therefore, the locations of damages in different scenarios were located with the EDI index and velocity vector. The results show that the proposed method determined the location of damage with the acceptable accuracy for low and moderate damage in damage scenarios.تشخیص آسیب یکی از ابزارهای مهم پایش سلامت سازه برای ارزیابی بهتر سازهها در طول عمر آنها میباشد. بسیاری از مطالعات به ارائه روشهایی برای تعیین محل آسیب با استفاده از مدلهای تحلیلی و آزمایشگاهی مانند مدلهای شاخص پرداختهاند. هدف اصلی این مقاله ارائه روش جدید تشخیص محل آسیب ترکیبی برای شناسایی مکانهای آسیب با استفاده همزمان از شاخصهای آسیبپذیری انرژی و جابهجایی میباشد. در بخش اول از طریق شاخص انرژی فرکانس آنی EDI و پاسخهای شتاب سازه به تعیین الگوهای آسیب پرداخته شده است. در بخش دوم بهمنظور ارزیابی روش اول و همچنین ارائه روشی سریع برای ارزیابی آسیب از طریق شاخص آسیب جابهجایی که متشکل از شاخص قابلیت اعتماد خطا β و شاخص تابع چگالی احتمال نرمال G(x)با استفاده از پاسخهای نسبی جابهجایی سازه ASCE ارائه گردیده است. نوآوری این روش استفاده همزمان از پاسخ شتاب- جابهجایی در طی یک فرایند است که در ارزیابی سریع الگوهای آسیب مؤثرتر است. برای صحتسنجی روشهای ارائه شده، علاوه بر الگوهای آسیب موجود در مسئله شاخص، آسیب جدید دیگر مورد بررسی قرار گرفته است. تجزیهوتحلیل گسترده نشان میدهد که روش پیشنهادی، محل دقیق آسیب وارده به سازه را با دقت کافی و سرعت مناسب تعیین مینماید.https://www.bese.ir/article_243558_c878d0540f3c48772ed55d273489bf2f.pdfپژوهشگاه بین المللی زلزله شناسی و مهندسی زلزلهفصلنامه علوم و مهندسی زلزله2476-60977420201221Investigating the Effects of Underground Structures on the Scattering of Seismic Waves Reaching the Ground Surfaceبررسی اثرات سازه زیر زمینی بر روی تفرق امواج لرزه ای رسیده به سطح زمین101112243555FAمحسنایثاریدانشکده مهندسی عمران، دانشگاه تبریز، تبریز، ایرانرضاتاری نژاددانشکده مهندسی عمران، دانشگاه تبریز، تبریز، ایرانسید کاظمرضویدانشکده مهندسی عمران، دانشگاه تبریز، تبریز، ایرانJournal Article20190713In the present study, the FLAC Software is applied to investigate the different models of Daikai Subway Station, which was heavily damaged in the 1995 Kobe earthquake. The responses of the ground surface were analyzed with and without the presence of the structure to examine the phenomenon of amplification. The maximum amplification of 2.2 was obtained in the middle point of the structure, while the maximum acceleration was obtained in the point adjacent to the edge of the structure. The points located on the contact surface of the structure were experienced the highest amplification in comparison with the case of without structure. In the next step, different models of the structure with the central column thicknesses of 10-70 cm were developed to investigate the lateral displacements of the structure. The results indicated that the maximum lateral displacement of the structure between the maximum and minimum thicknesses was smaller than 10 mm.<br />1. Introduction<br />Today, the importance of underground structures such as subways is more than ever because of their no doubt critical roles in solving traffic problems in metropolises. Since subways are typically constructed in populated cities and pass beneath major commercial and economic centers, they can impose considerable casualties and economic losses if damaged or destroyed [1]. For many years, the performance of underground structures was believed to be better than that of ground structures when subjected to the earthquakes. However, the experience of some earthquakes demonstrated that underground structures might undergo large deformations or even major failures [2-3]. The present study models a real underground structure, which experienced serious damages during the earthquakes, by the finite difference method (FDM) using the FLAC Software. The present study aims to investigate the effects of the structure on the scattering of the seismic waves reaching the ground surface. The quantity and pattern of relative displacements for different sizes of the central column (destroyed by the earthquake), the relationship between the recorded accelerations on the ground surface and the effects of the presence of the structure on the acceleration response are investigated.<br />2. Numerical Modeling<br />The numerical modeling was performed by the Finite Difference Method using FLAC 3D software. The static mode of the gravity loading was employed to define stresses induced by the soil weight. Roller supports were used for vertical boundaries, while pinned supports were applied to the bottom of the model. These boundaries reflect waves into the model and thus they cannot be applied to the dynamic mode.<br />Conclusion<br />The present study employed the FDM to investigate the seismic responses of a real subway station. The obtained results are as follows: <br />1. Comparison of the amplification induced by the effects of the structure on the scattering of seismic waves indicated that the scattering effects increased near the edges of the structure, leading to the large amplification factors at the end area.<br />2. The maximum amplification of 2.2, and a 25 percent increase was obtained in the maximum acceleration of Daikai under the 1995 Kobe earthquake.<br />References<br />1. Dowding, C.H. and Rozan, A. (1978) Damage to rock tunnels from earthquake shaking. ASCE J. Geotech. Eng. Div., 104, 175-91.<br />2. Jing-Ming, W. and Litehiser, J.J. (1985) The distribution of earthquake damage to underground facilities during the 1976 Tang-Shan earthquake. Earthquake Spectra, 1, 741-57.<br />3. Sharma, S., and Judd, W.R. (1991) Underground opening damage from earthquakes. Eng. Geol., 30, 263-76.شناخت و درک واقعبینانه از تفرق امواج لرزهای ناشی از سازههای زیرزمینی که تحت عنوان پدیده بزرگنمایی در سطح زمین شناخته میشود، میتواند منجر به طراحی ایمن و مقاوم سازههای ساخته شده در سطح زمین شود. در پژوهش حاضر، با استفاده از نرمافزار تفاضل محدود FLAC مدلهای مختلفی از ایستگاه مترو دایکای که تحت زلزله کوبه در سال 1995 بهشدت آسیب دید، مورد بررسی قرار گرفت. برای مطالعه پدیده بزرگنمایی، پاسخهای سطح زمین در دو حالت بدون سازه و با وجود سازه بررسی شد. نتایج نشان میدهد که بیشترین مقدار بزرگنمایی به مقدار 2/2 در نقطه میانی سازه و بیشترین پیک شتاب در نقطه مجاور لبه سازه ثبت شد. در بین نقاط مختلف روی لایههای خاک، نقاطی که در تماس با سازه قرار داشتند بیشترین مقدار افزایش نسبت به حالت بدون سازه را تجربه کردهاند. در ادامه مدلهای مختلفی از سازه در ضخامتهای 10 تا 70 سانتیمتری ستون مرکزی تهیه و تغییر مکانهای جانبی سازه مورد بررسی قرار گرفت. نتایج نشان میدهد که حداکثر تغییر مکان جانبی سازه بین کمترین و بیشترین ضخامت کمتر از 10 میلیمتر میباشد.https://www.bese.ir/article_243555_c51cb512ecd182ea9ee61ff25d502747.pdfپژوهشگاه بین المللی زلزله شناسی و مهندسی زلزلهفصلنامه علوم و مهندسی زلزله2476-60977420201221Seismic Reliability of Steel Frames Systems Equipped with Viscose Dampersقابلیت اعتماد لرزهای و منحنی شکست سیستمهای سازهای قاب خمشی فولادی مجهز به میراگر ویسکوز113127243556FAفؤادکندریدانشکده فنی و مهندسی، دانشگاه خوارزمی، تهران، ایرانغلامرضانوریدانشکده فنی و مهندسی، دانشگاه خوارزمی، تهران، ایران0000-0002-2936-599Xپیمانهمامیدانشکده فنی و مهندسی، دانشگاه خوارزمی، تهران، ایرانJournal Article20190915Structures may undergo unexpected changes during the construction and assembly of materials and their sections, or may be affected by earthquakes of various magnitudes and mechanisms during their life time. Each of these changes may affect the possibility of structural failure. These variable parameters are examined as reliability. By the theory of reliability, uncertainties of the structural and non-structural parameters such as material properties, geometrical dimensions, and seismic inputs, etc. are taken into account. The theory of reliability with fragility curves shows which part of the model is more affected by uncertainties. These curves show the probability of passing a certain level of damage against seismic parameters.<br />This paper presents a reliability analysis of the steel moment resistance structure equipped with viscous dampers applied by uncertainties in seismic input and structural properties such as ductility and fragility of beam and column, damping coefficient of dampers by obtaining fragility curves. Given the intensity and frequency of near-field earthquakes that cause large displacements in the structure and it is possible that the braces on which the dampers are mounted enter the nonlinear domain before activating or collapse of viscous damper, the analysis are examined only on the basis of far field earthquakes. Nonlinear incremental dynamic analysis (IDA) is applied to the 5 and 10-story frames modeled in OpenSees software under a set of 22 pairs of acceleration records. Also, to reduce the time and number of analyzes, the response level method has been used in this study. The number of scenarios required for the response level method is obtained from the box-Wilson method. According to the results obtained from the analyzes (mean and deviation from the failure level criterion), a second-degree surface procedure for mean and standard deviation has been internalized. The Monte Carlo method also simulated 10,000 fragility curves to determine the final fragility curve. The results show that increasing structural story has a significant effect on the uncertainty in the response of structures. Based on the results, considering uncertainties in 5 and 10-story structure reduce failure spectral acceleration (Sa) by 7.1% and 9.3% respectively. Besides, in spectral values corresponding to the first mode of 5 and 10-story structures, with considering the uncertainty, probability of failure increase about 52.4% and 74.7%, respectively. The results showed that considering the uncertainty in some parameters increases and others reduce the capacity of the structure; however, by considering uncertainty in all mentioned parameters simultaneously, it will reduce the capacity and increase the probability of collapse.تئوری قابلیت اعتماد در کنار منحنیهای شکنندگی میزان تأثیرپذیری سازه از عدم قطعیتهای مدلسازی را نشان میدهد. این نمودارها احتمال فراگذشت از یک سطح خرابی مشخص را در مقابل پارامترهای لرزهای مشخص میکنند. در این مقاله تحلیل قابلیت اعتماد سیستم قاب خمشی تقویتشده با میراگر ویسکوز با اعمال عدم قطعیت در تحریک لرزهای ورودی و مشخصات سازهای و به دست آوردن منحنیهای شکنندگی انجام شده است. برای این منظور تحلیل دینامیکی افزایشی غیرخطی بر روی قاب 5 و10طبقه، مدلسازی شده در نرمافزار OpenSees تحت مجموعه 22 زوج شتابنگاشت رکوردهای حوزه دور اعمال شده است. همچنین به کمک روش مونتکارلو 10000 منحنی شکنندگی جهت تعیین منحنی شکست نهایی شبیهسازی شد. نتایج نشان میدهد با افزایش طبقات سازه تأثیر لحاظ نمودن عدم قطعیت در پاسخ سازهها قابلتوجه میباشد. بر اساس نتایج با در نظر گرفتن عدم قطعیت در قابهای 5 و 10 طبقه در احتمال فرو ریزش50 درصد مقدار شتاب طیفی (Sa) متناظر بهترتیب 1/7 درصد و3/9 درصد نسبت به زمانی که عدم قطعیت در مشخصات سازهای در نظر گرفته نشود، کمتر میگردد. همچنین در مقادیر شتاب طیفی متناظر با پریود اصلی سازه برای قاب 5 طبقه و 10طبقه با در نظر گرفتن عدم قطعیت درصد احتمال شکست بهترتیب حدود 4/52 درصد و 7/74 درصد بیشتر میشود.https://www.bese.ir/article_243556_aff636516d430724f52e12d8404ad937.pdfپژوهشگاه بین المللی زلزله شناسی و مهندسی زلزلهفصلنامه علوم و مهندسی زلزله2476-60977420201221Determination of the Modal Information of Structures under Impact Loads Using Proper Orthogonal Decompositionاستفاده از روش تجزیه متعامد بهینه جهت استخراج اطلاعات مودال سازهها تحت اثر بارهای ضربهای129145243557FAامیرزایری بغلانی نژادگروه سازه، دانشکده مهندسی عمران، دانشگاه تربیت دبیر شهید رجایی، تهران، ایرانموسیمحمودی صاحبیگروه سازه، دانشکده مهندسی عمران، دانشگاه تربیت دبیر شهید رجایی، تهران، ایرانJournal Article20190916Determining the modal characteristics of structures such as natural frequencies and damping ratios is one of the most important issues in structural engineering. In this regard, providing a low-cost and robust experimental method against all types of noises is very important. In the present paper, a new algorithm for determining the natural frequencies and damping ratios of structures under impact loads using the proper orthogonal decomposition technique is presented. This method uses the vibrational response of the structure to impact loads, without the need to calculate the impact magnitude. One of the strength points of the proposed methodology is the accumulation of laboratory noises in the latest modes. In other words, in the process of calculating the frequencies related to the first few modes, laboratory noise does not enter the calculations and will be aggregated in higher modes that are less important. The feasibility and efficiency of the new method was evaluated using numerical simulations as well as laboratory validation. In this research, four pure numerical models were used to evaluate and validate the accuracy of the proposed algorithm. These models include a simply supported beam, a two-dimensional portal frame, a three-dimensional truss and a clambed-clambed beam. The natural frequencies and damping ratios of the mentioned models were calculated using the new method, then the results were compared with those that obtained from the finite element method. Very good agreement was observed between the results of the two methods. For further investigation, various states such as the effect of the noise on the results, the effect of multiple impact loads and the effect of the number of sensors were also studied. The results of these numerical studies showed that the proposed method is very stable and robust against laboratory noises and has small errors. Acceptable results can be obtained in using a few numbers of sensors (even one sensor) and repetition of the test. Also, the results for multiple impact loads are almost similar to the results obtained from the excited state with a simple impulse load. In this study, in order to further investigate the efficiency of the POD based method, a small-scale laboratory model developed at Shahid Rajaee Teacher Training University and whose modal information has already been calculated using the closed image processing method, was studied. Good agreement was observed between the results of the two methods, so it can be another reason for the efficiency and capability of the new method.<br />Based on various studies, it was concluded that for structures with low damping, the proposed method has acceptable accuracy and with increasing damping, the accuracy of the results gradually decreases. Since conventional structures have a relatively low attenuation, the proposed algorithm is very suitable for determining their modal information. The proposed method due to the availability, cheapness and no need for complex experimental tools can be used as a useful algorithm to determine the modal information of a structure as well as control the results obtained from other experimental methods.تعیین مشخصات مودال سازهها نظیر فرکانسهای طبیعی و نسبتهای میرایی، از جمله موضوعات پر اهمیت در مهندسی سازه به شمار میرود. در این راستا ارائه یک روش آزمایشگاهی کمهزینه و پایدار در برابر انواع نویزها همواره حائز اهمیت است. در مقاله حاضر روش جدیدی جهت تعیین فرکانسهای طبیعی و نسبتهای میرایی سازهها با استفاده از تکنیک ریاضی تجزیه متعامد بهینه ارائه شده است. در این روش از پاسخ ارتعاشی سازه به بارهای ضربهای، بدون نیاز به محاسبه اندازه ضربه استفاده میشود. یکی از نقاط قوت روش پیشنهادی، تجمیع نویز آزمایشگاهی در مودهای بیاهمیت بالا میباشد. بهعبارتدیگر در فرایند محاسبه فرکانسهای مربوط به چند مود اول، نویز آزمایشگاهی وارد محاسبات نمیشود و در مودهای بالاتر که از اهمیت کمتری برخوردارند، تجمیع خواهد شد. کارایی روش جدید با استفاده از شبیهسازی عددی و نیز صحتسنجی آزمایشگاهی مورد ارزیابی قرار گرفت. روش پیشنهادی به دلیل سادگی، کمهزینه بودن و عدم نیاز به تجهیزات پیشرفته آزمایشگاهی میتواند بهعنوان ابزاری مفید جهت تعیین اطلاعات مودال یک سازه و همچنین کنترل نتایج بهدستآمده از روشهای آزمایشگاهی دیگر بهکار رود.https://www.bese.ir/article_243557_53d2d950d3385007e8be949ea45e7da7.pdfپژوهشگاه بین المللی زلزله شناسی و مهندسی زلزلهفصلنامه علوم و مهندسی زلزله2476-60977420201221Literature Review of Estimating the Bearing Capacity of Rough Footings by the Stress Characteristic Lines Methodتکنیکهای ارزیابی ظرفیت باربری پی سطحی زبر با روش خطوط مشخصهی تنش147161243559FAشرویناحمدیگروه مهندسی عمران، واحد تهران جنوب، دانشگاه آزاد اسلامی، تهران، ایرانمحسنکمالیانپژوهشکده مهندسی ژئوتکنیک، پژوهشگاه بین المللی زلزله شناسی و مهندسی زلزله، تهران، ایران0000-0001-5523-3520فرجالهعسکریپژوهشکده مهندسی ژئوتکنیک، پژوهشگاه بین المللی زلزله شناسی و مهندسی زلزله، تهران، ایران0000-0001-9157-2693Journal Article20210417The footing of all buildings will be located on the ground. Thus, civil engineers are always faced with the problem of estimating the bearing capacity of shallow foundations. So far, various analytical and numerical methods such as stress characteristic lines method, limit equilibrium method, lower and upper boundary limit analysis methods, as well as finite difference and finite element methods, have been utilized to evaluate the bearing capacity of shallow foundations.<br />In the meantime, the stress characteristic method, also known as the “slip-line method”, has been shown to be an efficient and useful technique for solving the bearing capacity problem due to its simplicity, high speed in calculations, and no need for neither meshing nor complex soil behavior models. These advantages attracted the attention of many researchers over the past three decades to the stress characteristic lines method.<br />Obviously, the bearing capacity of footings would be reduced, if they are subjected to seismic loadings or if they are located on slopes or adjacent to them. On the other hand, experimental as well as numerical studies reported in the literature revealed that the roughness of the footing would have an important effect on its bearing capacity. Hence, estimating the bearing capacity of a footing taking into consideration the roughness of its contact surface constitutes one of the important issues in foundation engineering, which of course has been studied only by a small number of researchers.<br />This paper intends to review the most important works reported in the literature concentrated on the estimation of the bearing capacity of rough shallow foundations by the stress characteristic lines method. We wanted to find the answer to the question that how different researchers considered the roughness factor in evaluating the bearing capacity of footings by the stress characteristic lines method? For this purpose, firstly the mathematical formulation of the stress characteristics method was briefly reviewed. After that, the numerical algorithm of implementing this method to estimate the bearing capacity of smooth as well as rough footings, the corresponding boundary conditions, and calculation procedure for determining the bearing capacity coefficients 𝑁𝑐, 𝑁𝑞, and 𝑁𝛾 were discussed. Various techniques for taking into consideration the roughness factor, especially their assumptions, were investigated and compared with the smooth footing case.<br />An important point to keep in mind when calculating the bearing capacity of rough footing is the non-uniformity of the contact pressure inclination. Various researchers employed different assumptions and obtained quite different results. Some researchers considered uniform roughness by regarding different inclination for the contact pressure along the soil-footing interface or by placing a non-plastic curved wedge immediately beneath the footing. Some other researchers adopted the non-uniform roughness and solved the roughness problem by assigning a non-plastic curved wedge under the footing, or by considering it as a line and two types of parabolic equations, or even without regarding any pre-assumption failure pattern. The influences of these assumptions and methods on estimating the bearing capacity of rough footings were investigated in detail and discussed. Finally, the most important challenges ahead in evaluating the bearing capacity coefficients of rough footing by the stress characteristic lines method were summarized and explained.برآورد ظرفیت باربری پیهای سطحی همواره مورد نیاز مهندسین بوده است. روش مشخصههای تنش برای حل مسائل ظرفیت باربری پیهای مستقر بر فضای نیمبینهایت به دلیل سادگی، سرعت در انجام محاسبات و عدم نیاز به مشبندی و مدلهای رفتاری پیچیدهی خاک، تکنیکی کارآمد و مؤثر است. در نوشتار حاضر، ابتدا فرمولبندی ریاضی روش خطوط مشخصهی تنش بهاختصار تشریح شده است. نحوهی بهکارگیری آن در ارزیابی ظرفیت باربری پیهای زبر و تکنیکهای مختلف حل مسئلهی زبری بررسی شده است. سپس، چگونگی در نظر گرفتن یکنواختی یا غیریکنواختی شیب فشار تماس وارد بر کف پی و تأثیر آن بر ضرایب ظرفیت باربری پیها مورد بحث و بررسی قرار گرفته است. اثرات ناشی از نحوهی ملحوظ نمودن زبری به روشهای مختلف بر این ضرایب سهگانه با یکدیگر مقایسه شده است. در نهایت، به چالشهای جدید و مهم پیش رو در استفاده از روش مشخصههای تنش در بررسی اثرات زبری بر ظرفیت باربری پی پرداخته شده است.https://www.bese.ir/article_243559_8b1154b0040c183c84871ae97624c0b1.pdf