تحلیل ارتعاش آزاد تیر رایلی با میرایی ویسکوالاستیک غیرمحلی به روش گالرکین، فضای پوچ ماتریس و بسط سری نیومن

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

free vibration analysis of viscoelastic non-locally damped rayleigh beam using galerkin method, null space of the matrix & neumann series expansion

 while studying large-scale systems, non-local damping definition can beneficially model contact shear and long-range forces resulting from adjacent and non-adjacent elements in the set of interconnected dampers, damping patches and foundations which are modeled as non-local domains. if two or three dimensional systems are considered one dimensional (e.g. analyzing three dimensional beams based on euler-bernoulli, rayleigh or timoshenko theories which simplifies the behavior of structures), the concept of non-local damping models arises to improve the accuracy of numerical results. even though defining dissipative forces which are dependent on more parameters and quantities helpfully boost the validity of results compared to experimental cases in labs and three dimensional numerical analysis done by software, many researchers have widely employed viscous damping model to demonstrate damped behavior of the structures in the sake of simplicity. actually viscous damping does not model accurately the dissipative behavior of real systems and practical structures often demonstrate some kind of viscoelasticity while vibration.  in the recent study, external damping force at any point in the domain is influenced by the past history of velocity and long-range interactions through convolution integral over proper kernel functions. as a consequence of applying laplace transformation and using galerkin method, the integro-partial differential equation of rayleigh beam as a distributed parameter system turns to an ordinary differential equation governing a discrete system with finite degrees of freedom. galerkin method is established based on error minimization of assumed mode method and despite rayleigh and rayleigh-ritz methods can suitably analyze nonconservative systems including damped beams. corresponding undamped mode shapes of rayleigh beam which satisfy essential boundary conditions are chosen as the best admissible functions to expand the trial response of equation of motion. in order to get cm-1 continuous functions and finite weighted residual integral, the equation is presented in the weak form. afterwards stiffness, damping and two types of mass matrices are determined with respect to generalized coordinates. to get scaled mode shapes the mass change method is considered to evaluate scaling factor, then the results are mass normalized. by equating dynamic stiffness matrix to zero and solving the resulting algebraic equation, complex and real eigenvalues are obtained which are respectively elastic and non-viscous modes in stable systems. it’s noteworthy to announce that contrasted with the viscously damped systems the degree of algebraic equation in the case of viscoelastic non-locally damped rayleigh beam would generally be more than 2n (which n stands for total degrees of freedom). accordingly eigenvectors are found based on gaussian elimination method and null space of the dynamic stiffness matrix. additionally, neumann series expansion is developed to determine eigenvectors of vibrating systems effected by inertia. finally, numerical results of rayleigh and euler-bernoulli beams are compared. very thin rayleigh beams are verified and show great accuracy but when the thickness increases and inertial effect becomes bold, results of rayleigh and euler-bernoulli beams are not matched anymore. furthermore it’s understood from the graphs and tables that when the characteristic parameter of nonlocality and relaxation constant increase and tend to infinity, the nonlocal and non-viscous effects decline.