بررسی عدم قطعیت نتایج آزمایش‌های ارتعاش محیطی در روش‌های زیرفضای تصادفی

در روش های زیرفضای تصادفی مهترین عامل تاثیرگذار بر مشخصات دینامیکی مستخرج، ابعاد ماتریس هانکل بوده که شامل تعداد سطرها و ستون های آن است. استفاده از ابعاد کوچک ماتریس احتمال عدم شناسائی قطب های پایدار را داشته و انتخاب ابعاد بسیار بزرگ علاوه برآنکه احتمال بوجود آمدن قطب های مجازی و بایاس را افزایش داده، باعث افزایش هزینه محاسباتی نیز می گردد. از اینرو در این تحقیق قصد براین است که تعداد سطر های بلوکی و ستونهای مطلوب ماتریس هانکل در روش زیرفضای تصادفی متعادل طوری محاسبه شود که ضمن پوشش قطب های موجود، هزینه محاسباتی حداقلی نیز داشته باشد. برای این منظور از معیار عدد شرطی ماتریس هانکل و اندیکاتور انرژی استفاده خواهد شد. مراحل کار بدین طریق است که با استفاده از اندیکاتور انرژی، حداکثر مرتبه موثر سیستم برای سطرهای بلوکی مختلف ماتریس هانکل استخراج گردیده و سپس نمودار ماکزیمم عدد شرطی برای آن ها رسم می شود. تعداد سطرهای بلوکی مطلوب ماتریس هانکل از همگرائی معیار عدد شرطی با توجه به تعریف این معیار استخراج می شود. به منظور صحت سنجی روش پیشنهادی، از نتایج آزمایش ارتعاش محیطی پل روگذر شهر نمین که توسط محققین این مقاله انجام گرفته، استفاده شده است. نتایج این تحقیق نشان دهنده این است که از بُعد 352،  نمودار عدد شرطی نسبت به تغییر سطرها بلوکی حساسیت خود را از دست داده که بعنوان بعد مطلوب انتخاب شده است. همچنین در بعد بهینه با استفاده از تغییرات انحنای نمودار اندیکاتور انرژی مرتبه  66 بعنوان مرز بد شرطیدگی ماتریس سیستم، محاسبه گردیده که پیدایش مودهای بایاس و محاسباتی در نمودار پایداری، موید مرتبه محاسبه شده است.

the uncertainty of ambient vibration test results in stochastic subspace methods

in stochastic subspace methods, the most important factor influencing the dynamic specifications is the dimensions of the hankel matrix include the number of rows and columns. using small matrix dimensions is unlikely to identify existing poles, and selecting very large dimensions not only increases the likelihood of virtual and bias poles but also increases computational costs. in this study, it is intended that the optimal dimensions of the hankel matrix in the balanced stochastic subspace method be calculated in such a way that in addition to covering the existing poles, it also has a minimum computational cost. for this purpose, the condition number of the hankel matrix and energy indicator is used in two steps. the steps are as follows: first, calculate the optimal order of each cycle, and then use the optimal order to draw the condition number of the system matrix for different dimensions and calculate the desired dimension from its convergence. to verify the accuracy of the proposed method, the ambient vibration test of the namin entrance bridge has been used. this bridge is located at the entrance of namin city, 25 km from the center of ardabil province, iran, which includes two spans of 27.10m with a concrete deck. the deck of the bridge is located on beams with i sections, which are 2.5m away from each other, and the whole set of beams and deck is located on a system of foundations and piles with a diameter of 120cm. this bridge being the only entrance to the city and is exposed to various traffic loads, it was necessary to monitor the dynamic characteristics of the bridge as modal frequencies and damping ratios to evaluate the performance and ensure the health of the bridge structure. according to the numerical analysis and the length of the data (12000), the minimum order and the maximum number of cycles are 22 and 55, respectively. by diverging the curvature of the energy indicator graph, the optimal order is determined in the initial 5 12% of the singular values of cycles. for example, the maximum order of the 6th cycles was obtained, 28 62. also, from the convergence of the maximum condition number of cycles from the 8th   cycle, the optimal dimension was selected 352. in a general summary, it can be said that the use of the energy indicator concept in finding the effective order of the stability diagram has a significant effect on reducing the uncertainty of the extracted results. so that from the three identified stable poles, two poles have been extracted in the effective order area. also, using the concept of conditional number to find the optimal dimension of the system was effective, so that by drawing a stability diagram for the 15th cycle, it was found that the identified modal characteristics were not significantly different from the results of the optimal cycle (8th). finally, the extracted modal properties have an acceptable agreement with the numerical model and frequency domain decomposition method (fdd). the modal frequencies of both methods (fddand b ssi) have a good correlation but the damping ratios were very different. in frequency domain methods the damping ratios being very sensitive to the quality of data collection, one can expect that the results of the subspace method are closer to reality.