بررسی کمّی میزان دقت روش شتاب مرتبه 3 نسبت به روش‌های دیگر عددی برای بارگذاری هارمونیک؛ مطالعه موردی: یک سازه سخت‌شونده یک درجه آزاد

در این مقاله، پاسخ جابه جایی یک سیستم یک درجه آزاد سخت‌شونده غیرخطی، تحت طیفی از بارگذاری‌های سینوسی با فرکانس های مختلف و دامنه نیروهای متفاوت مورد بررسی قرار گرفته است. تغییرات سختی این سیستم، رابطه درجه 2 با مقدار جابه جایی داشت. از 9 روش عددی که شامل روش شتاب مرتبه 3، روش شتاب مرتبه 2، روش نیومارک (روش شتاب خطی)، روش نیومارک (روش شتاب متوسط)، روش ویلسون-تتا، روش تفاضل مرکزی، روش جنینگز، روش جنینگز ارتقاء یافته و روش ابداعی بود، برای محاسبه پاسخ جابه جایی این سیستم استفاده گردید. با استفاده از روش ابداعی و با استفاده از t∆ برابر با 0/001 ثانیه، پاسخ جابه جایی سیستم محاسبه شد و به‌عنوان پاسخ دقیق سیستم قلمداد گردید. سپس، برای t∆های 0/001، 0/002، 0/004، 0/006، 0/008 و 0/01 ثانیه و با استفاده 9 روش عددی معرفی شده در بالا پاسخ جابه جایی این سیستم به ازای دامنه بارگذاری های برابر با 200، 400، 600، 800 و 1000 نیوتن محاسبه گردید و مقادیر جذر میانگین مربعات خطا و ضریب تغییرات خطای این 9 روش به ازای t∆های تعریف شده در بالا محاسبه شد که نتایج کار به‌صورت نمودارهایی میله ای در مقاله آورده شده است. در مجموع، روش ابداعی، روش جنینگز ارتقا یافته و روش جنینگز، به ترتیب صعودی کمترین میزان خطا را داشتند و بیشترین خطا مربوط به روش ویلسون بود. از لحاظ ضریب تغییرات، روش های ویلسون، نیومارک شتاب متوسط و تفاضل مرکزی، به ترتیب صعودی، کمترین ضریب تغییرات را دارا بودند و روش جنینگز ارتقا یافته، بیشترین ضریب تغییرات را داشت.

quantitative evaluation of the accuracy of the third-order acceleration method compared to other numerical methods for a case study of a stiffening single-degree-of-freedom structure under harmonic loading

this paper investigates the displacement response of a nonlinear, stiffening single-degree-of-freedom system subjected to a range of sinusoidal loadings with varying frequencies and force amplitudes. the stiffness variation in this system follows a quadratic relationship with displacement. a total of nine numerical methods were utilized to compute the displacement response of the system. these methods include the third-order acceleration method, the second-order acceleration method, the newmark’s linear acceleration method, the newmark’s average acceleration method, the wilson-theta method, the central difference method, the jennings method, the improved jennings method, and a newly developed method. using the developed method and a time step ∆t of 0.001 seconds, the system’s displacement response was computed and considered as the reference solution. subsequently, for time steps of 0.001, 0.002, 0.004, 0.006, 0.008, and 0.01 seconds, the displacement response of the system was calculated using the aforementioned nine numerical methods under loading amplitudes of 200, 400, 600, 800, and 1000 newtons. the root mean square error and the cv of the error for each method were calculated for the defined time steps. results are presented as bar charts within the paper.this paper investigates a nonlinear stiffening sdof system without damping, with a mass of 1000 kg and an initial stiffness of 100,000 n/m. the mass and stiffness specifications of this system were taken from chang (2009, pp. 289-297), where a two-degree-of-freedom shear system without damping is presented with a mass of 1000 kg for each floor.in this study, the displacement response of the system under consideration, calculated with a time step of 0.001 seconds, is regarded as the accurate response of the system. to justify this, the studied nonlinear stiffening sdof system was subjected to a range of sinusoidal loads with frequencies ranging from 1 hz to 30 hz, increasing by 0.1 hz (291 frequencies), with a loading amplitude of 1000 newtons (the maximum loading amplitude), and the system response was calculated with a maximum time step of 0.01 seconds. then, using three methods (the innovative method, newmark’s average acceleration method, and the jennings method) and with a time step of 0.001 seconds, the system’s response was calculated and considered as the exact response. by reducing the number of points and changing the time interval from 0.001 seconds to 0.01 seconds, the maximum error and average errors of the nine methods discussed in this study were computed. examining the results reveals that there is essentially no significant difference in the error values obtained by the three mentioned methods, especially for the average error values. however, since the innovative method assumes a constant stiffness value within each time step of 0.001 seconds, this method appears to offer a higher degree of accuracy.overall, the newly developed method, the improved jennings method, and the original jennings method, respectively, exhibited the lowest error, while the wilson-theta method showed the highest error. in terms of the coefficient of variation, wilson-theta, newmark’s average acceleration, and the central difference method, respectively, had the lowest coefficients of variation in ascending order, with the improved jennings method having the highest coefficient of variation