توضیحات

توجه : به همراه فایل word این محصول فایل پاورپوینت (PowerPoint) و اسلاید های آن به صورت هدیه ارائه خواهد شد

  مقاله بررسی بردارهای ریتز وابسته به بار و روش MPA در word دارای 173 صفحه می باشد و دارای تنظیمات و فهرست کامل در microsoft word می باشد و آماده پرینت یا چاپ است

فایل ورد مقاله بررسی بردارهای ریتز وابسته به بار و روش MPA در word  کاملا فرمت بندی و تنظیم شده در استاندارد دانشگاه  و مراکز دولتی می باشد.

توجه : توضیحات زیر بخشی از متن اصلی می باشد که بدون قالب و فرمت بندی کپی شده است

بخشی از فهرست مطالب پروژه مقاله بررسی بردارهای ریتز وابسته به بار و روش MPA در word

فصل اول: آنالیز دینامیکی با استفاده از بردارهای ریتز وابسته به بار
بخش اول: تحلیل دینامیکی
مقدمه
1-1- اصول اولیه تحلیل دینامیکی
2-1- تعادل دینامیکی
3-1- روش حل گام به گام
4-1- روش برهم نهی مدی
5-1- تحلیل طیف پاسخ
6-1- حل در حوزه فرکانس
7-1- حل معادلات خطی
بخش دوم: محاسبه بردارهای متعامد بر جرم و سختی
مقدمه
1-2- روش جستجوی دترمینانی
2-2- کنترل ترتیب استورم
3-2- متعامد سازی گرام اشمیت
4-2- تکرار زیر فضای بلوکی
5-2- حل سیستمهای منفرد
6-2- ایجاد بردارهای ریتز وابسته به بار
بخش سوم: کلیات روش LDR
1-3- روش جداسازی دو مرحله ای در تحلیل سازه ها
1-1-3- جداسازی مسائل خطی دینامیکی به وسیله برهم نهی مدی
2-3- استفاده از بردارهای ریتز در دینامیک سازه ها
1-2-3- روش ریلی برای سیستمهای تک درجه آزادی
3-3- تولید خودکار بردارهای ریتز وابسته به بار
4-3- تاثیر فرمول بندی اجزای محدود بر ایجاد بردارهای ریتز وابسته به بار
1-4-3- ماتریس جرم
2-4-3- بردار بارگذاری
1-2-4-3- محتوای فرکانسی
2-2-4-3- توزیع مکانی
بخش چهارم: ارتباط میان الگوریتم بردارهای ریتز وابسته به بار و روش Lanczos
1-4- روش Lanczos
عنوان صفحه
2-4- خواص اساس بردارهای ریتز وابسته به بار
3-4- نکاتی در مورد تعامد بردارهای پایه ریتز وابسته به بار
4-4- تحلیل سیستمهای با میرایی
1-4-4- روند حل برای میرایی متناسب (با ماتریس سختی)
2-4-4- روند حل برای میرایی غیر متناسب
5-4- فلسفه اساسی فراسوی بردارهای ریتز وابسته به بار
بخش پنجم: توسعه تخمین خطا برای بردارهای ریتز وابسته به بار
1-5- تخمین های خطای مکانی برای ارائه بارگذاری
2-5- ارائه بارگذاری به وسیله پایه بردارهای ریتز وابسته به بار
3-5- تخمین های خطا با استفاده از مجموع بارهای ارائه شده
4-5- تخمین خطا براساس معیار اقلیدسی بردار خطای نیرو
5-5- روشهای جمع بندی برای آنالیز برهم نهی مستقیم بردار
1-5-5- روش تصحیح استاتیکی
2-5-5- روش شتاب مدی
6-5- رابطه میان بردارهای ریتز وابسته به بار و حل مقدار ویژه دقیق
بخش ششم: الگوریتمی جدید برای ایجاد بردارهای ریتز وابسته به بار
1-6- استقلال خطی بردارهای ریتز وابسته به بار
1-1-6- روش Lanczos و مساله از دست دادن تعامد
2-1-6- بردارهای ریتز وابسته به بار و مساله از دست دادن تعامد
3-1-6- باز متعامد سازی انتخابی
4-1-6- کاربرد کامپیوتری متعامد سازی انتخابی
2-6- تنوع محاسباتی الگوریتم بردارهای ریتز وابسته به بار
1-2-6- بردارهای ریتز LWYD
2-2-6- کاربرد کامپیوتری با استفاده از فرم کاهش یافته سه قطری
3-6- کاربرد عددی روی سیستمهای ساده سازه‌ای
1-3-6- حل مثال با استفاده از برنامه CALSAP
2-3-6- توضیح مدل ریاضی
3-3-6- ارزیابی گونه های محاسباتی الگوریتم ریتز
بخش هفتم: تحلیل دینامیکی غیرخطی با برهم نهی مستقیم بردارهای ریتز
1-7- منبع و حد رفتار غیرخطی
2-7- تکنیک های راه حل برای تحلیل دینامیکی غیرخطی
3-7- روشهای انتگرال گیری مستقیم
عنوان صفحه
4-7- روشهای برهم نهی برداری
5-7- گزینش بردارهای انتقال برای روشهای برهم نهی
6-7- خط مشی های حل سیستمهای غیرخطی کلی
7-7- خط مشی های حل سیستمهای غیرخطی محلی
بخش هشتم: توصیف فیزیکی الگوریتم ریتز و ارائه چند مثال
1-8- مقایسه حل با استفاده از بردارهای ویژه و بردارهای ریتز
مثال
مثال
مثال
بخش نهم: تحلیل دینامیکی با استفاده از بردارهای ریتز
1-9- معادله حرکت کاهش یافته
نتیجه
مراجع فصل اول
ضمیمه
فصل دوم: آنالیز استاتیکی فزاینده غیرخطی مودال (MPA)
بخش اول: آنالیز استاتیکی فزاینده غیرخطی
1-1- روندهای تحلیلی
2-1- پیدایش روش غیرخطی استاتیکی
3-1- فرضیات اساسی
1-3-1- کنترل براساس نیرو یا تغییر مکان
2-3-1- الگوهای بارگذاری
3-3-1- تبدیل سازه MDF به SDF
4-3-1- تغییر مکان هدف
5-3-1- حداکثر شتاب زمین
4-1- روش آنالیز استاتیکی غیرخطی
5-1- روش گام به گام در محاسبه منحنی ظرفیت
1-5-1- روش گام به گام محاسبه منحنی ظرفیت
6-1- محدودیتهای POA
بخش دوم: MPA
1-2- معادلات حرکت
2-2- معرفی سیستمهای مورد بررسی و حرکت زمین
3-2- روند تقریبی تحلیل
1-3-2- بسط مدی نیروهای موثر
2-3-2- ایده اساسی
4-2- روشUMRHA
1-4-2- سیستمهای خطی
2-4-2- سیستمهای غیرخطی
5-2- MPA
1-5-2- سیستمهای الاستیک
2-5-2- سیستمهای غیرالاستیک
6-2- خلاصه MPA
7-2- برآورد روش

بخشی از منابع و مراجع پروژه مقاله بررسی بردارهای ریتز وابسته به بار و روش MPA در word

 1- Building Seismic Safety Council (1997) NEHRP Guidelines for the Seismic Rehabilitation of Buildings, FEMA-273, Federal Emergency Management Agency, Washington, D.C

2- American Society of Civil Engineers (2000). Prestandard and Commentary for the Seismic Rehabilitation of Buildings, FEMA-356, Federal Emergency Management Agency, Washington, D.C

3- Fajfar, P., and fischinger, M. (1988). N2- a method for nonlinear seismic analysis of regular structures, Proc., 9 th World Conf. Earthq. Engrg., 5:111-116, Tokyo-kyoto, Japan

4- Krawinkler, H,, and Seneviratna, G.D.P.K. (1998). Pros and cons of a pushover analysis of seismic performance evaluation. Engrg. Struc., 20(4-6):5452-

5- Kim. B., and D’Amore, E. (1999). Pushover analysis procedure in earthquake engineering, Earthq. Spectra, 13(2):417-

6- Gupta, B., and Kunnath, S.K. (2000). Adaptive spectra- based pushover procedure for seismic evaluation of structures, Earthq. Spectra, 16(2): 367-

7- Sasaki, k.k., Freeman, S.A., and Paret, T.F. (1998)Multimode pushover procedure (MMP) – A method to identify the effects of higher modes in a pushover analysis, Proc., 6 th U.S. Nat. Conf. Earthq. Engrg., Seattle, Washington

8- Chopra, A.K., and Goel, R.K. (2002). A modal pushover analysis procedure for estimating seismic demands for buildings, Earthq. Engrg. Struc. Dyn., 31(3): 561-

9- Chintanapakdee, C, and Chopra, A.K. (2003). Evaluation of modal pushover analysis using generic frames, Engrg. Struc. Dyn., 32(3): 417-

10- Kilar, V., and Fajfar,P. (1997). Simple push-over analysis of asymmetric buildings, Earthq. Engrg. Stuc. Dyn., 26(2):233-

11- Moghadam, A.S., and Tso, W.K. (1998). Pushover analysis for asymmetrical multistory buildings, Proc., 6th U.S. Nat. Conf. Eagrg., EERI, Oakland, Calif., 13 pgs

12- Chopra, A.K. (2001). Dynamics of Structures: Theory and Applications to Earthquake Engineering. 2nd Edition, New jersey: Prentice Hall

13- Chopra, A.k. and Chintanapakdee, C., (2003) Inelastic deformation ratios for design and evaluation of structures: Single-  degree- of- freedom bilinear systems, ASCE, J, Struc. Engrg., to appear

14- Chopra, A. K., Goel, R. K., and Chintanapakdee, C. (2003) Statistics of single- degree- of- freedom estimate of displacements for pushover analysis of buildings, ASCE, J. Struc. Engrg., 129:1-

15- Yang, Pu, Wang, Yayong. ” A study on improvement of pushove Analysis”, 12 WCEE

16- J. Skokan, Matthew, and C. Hart, Gary, “Reliability of Nonlinear static methods for seismic per formance prediction of steel frame buildings”, 12 WCEE

مقدمه

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

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

شکل 1-1- ایده آل سازی سازه با جرم گسترده

استفاده از بردارهای ویژه، برای کاهش اندازه سیستمهای سازه‌ای یا ارائه رفتار سازه به وسیله تعداد کمی از مختصات های عمومی (تعمیم یافته) – در فرمول بندی سنتی – احتیاج به حل بسیار گرانقیمت مقدار ویژه دارد

یک روش جدید از تحلیل دینامیکی که نیاز به برآورد دقیق فرکانس ارتعاش آزاد و اشکال مدی ندارد توسط ویلسون Wilson یوان (Yuan) و دیکنز (Dickens) (1.17) ارائه شده است

روش کاهش، بردارهای ریتز وابسته به بار WYD Ritz vectors) که D, Y, W (حروف اختصاری نویسندگان)( بر مبنای بر هم نهی مستقیم بردارهای ریتز حاصل از توزیع مکانی و  بارهای مشخص دینامیکی می‌باشد. این بردارها در کسری از زمان لازم برای محاسبه اشکال دقیق مدی، توسط یک الگوریتم بازگشتی ساده بدست می‌آیند. ارزیابی‌های اولیه و کاربرد الگوریتم در تحلیل تاریخچه زمانی زلزله نشان داده است که استفاده از بردارهای ریتز وابسته به بار منجر به نتایج قابل مقایسه یا حتی بهتری نسبت به حل دقیق مقدار ویژه شده است

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

1-1- اصول اولیه تحلیل دینامیکی

تمام سازه های واقعی هنگام بارگذاری یا اعمال تغییرمکان به صورت دینامیکی رفتار می کنند. نیروهای اینرسی اضافی، با استفاده از قانون دوم نیوتن، برابر نیرو در شتاب می‌باشند. اگر نیروها و یا تغییر مکانها بسیار آرام اعمال شوند نیروهای اینرسی قابل صرفنظر کردن می باشند و یک تحلیل استاتیکی قابل انجام است. بنابراین می توان گفت، تحلیل دینامیکی بسط ساده ای از تحلیل استاتیکی می‌باشد

بعلاوه تمام سازه های حقیقی بالقوه دارای درجات آزادی نامحدودی می باشند. بنابراین بحرانی ترین قسمت در تحلیل سازه ایجاد مدلی با تعداد درجات آزادی محدود می باشد که دارای تعدادی اعضای تقریباً بدون جرم و تعدادی گره باشد، که بتواند رفتار سازه را به طور مناسبی تخمین بزند. جرم سازه را می توان درگره ها متمرکز نمود. نیز برای یک سیستم الاستیک خطی خصوصیات سختی اعضاء را می توان باصحت بسیار خوبی تخمین زد- باتوجه به داده های تجربی- هرچند تخمین بارگذاری  دینامیکی، اتلاف انرژی و شرایط مرزی می تواند بسیار مشکل باشد

با در نظر گیری موارد گفته شده برای کاهش خطاهای موجود لازم است تحلیل های دینامیکی متعدد با استفاده از مدلهای مختلف دینامیکی، بارگذاری و شرایط مرزی به کار گرفته شود و انجام حتی 20 آنالیز کامپیوتری برای طراحی یک سازه جدید و یا برآورد یک سازه موجود ممکن است لازم شود

 با توجه به تعداد زیادی آنالیزهای کامپیوتری که برای یک تحلیل دینامیکی نمونه لازم است  باید در کامپیوترها روشهای عددی مناسبی برای محاسبات به کار رود

2-1- تعادل دینامیکی

تعادل نیرویی برای یک سیستم چند درجه آزادی با جرم متمرکز شده، به صورت تابع زمان را می توان این گونه نوشت

F(t)I + F(t)D + F(t)S = F(t)                                                                                  (1-2-1)

F(t)I : بردار نیروهای اینرسی عمل کننده بروی جرم

F(t)D : بردار نیروی میرایی لزج، یا اتلاف انرژی می باشد

F(t)S : بردار نیروهای داخلی تحمل شده توسط سازه

F(t) : بردار بارهای اعمالی

معادله (121) برمبنای قوانین فیزیکی قرار دارد و برای هر دو دسته سیستمهای خطی و غیرخطی معتبر می باشد

برای بسیاری از سیستمهای سازه ای تخمین رفتار خطی برای سازه انجام می گردد تا معادله فیزیکی
(121) تبدیل به گروهی از معادلات دیفرانسیل مرتبه دوم خطی گردد

                                          (2-2-1)

که M ماتریس جرم، C ماتریس میرایی، K ماتریس سختی می باشند. بردارهای وابسته به زمان, ,, مقادیر مطلق تغییر مکان، سرعت و شتاب می باشند

برای بارگذاری زلزله F(t) نیروی خارجی برابر صفر می باشد. حرکت اساسی لرزه‌ای سه مؤلفه u(t)ig می باشند که در نقطه ای زیر پی ساختمان در نظر گرفته می شوند. بنابراین می توانیم معادله (122) را با توجه به, ,,که کمیاتی نسبی (نسبت به مؤلفه‌های زلزله) می باشند بنویسیم

بنابراین مقادیر مطلق تغییر مکان، سرعت و شتاب را می توان از معادله‌ (122) حذف نمود

u(t)a = u(t) + {rx} u(t)xg + {ry} u(t)yg + {rz} u(t)zg

(t)a = (t) + {rx}  (t)xg + {ry} (t)yg + {rz} (t)zg                                       (3-2-1)

ü(t)a= ü(t) + {rx} ü(t)xg + {ry} ü(t)yg + {rz} ü(t)zg

که {ri} برداری است که در درجات آزادی جهتی 1 می باشد و بقیه عناصر آن صفرند

با قرار دادن این معادله (3-2-1) در (2-2-1) داریم

Mü(t) + C(t) + Ku(t) = -Mx ü(t)xg – My ü(t)yg – Mz ü(t)zg                               (4-2-1)

که

Mi = M{ri}

روشهای کلاسیک گوناگونی برای حل معادله (1-4) وجود دارد که هرکدام دارای محاسن و معایب خاص خود می باشند که آنها را به صورت خلاصه بیان می کنیم

3-1- روش حل گام به گام

عمومی ترین روش تحلیل دینامیکی روش افزایشی است که معادلات تعادل در زمانهای Dt, 2Dt, 3Dt , …  حل می شوند. که تعداد زیادی از اینگونه روشهای افزاینده برای حل وجود دارد. در حالت عمومی این روشها شامل حل گروه کاملی از معادلات تعادل در هر افزایش زمان می باشند. در صورت انجام تحلیلی غیرخطی ممکن است لازم باشد تا ماتریس سختی سازه را شکل دهی مجدد نماییم

نیز امکان دارد در هر گام زمانی برای رسیدن به تعادل نیاز به تکرار داشته باشیم. از دیدگاه محاسباتی ممکن است حل یک سیستم با چند صد درجه آزادی زمان بسیاری طلب نماید

بعلاوه ممکن است نیاز داشته باشیم تا میرایی عددی یا مجازی را به دسته زیادی از این راه حلهای افزایشی برای بدست آوردن راه حلی پایدار اضافه کنیم. برای تعدادی از سازه های غیرخطی که تحت تأثیر حرکت زمین قرار گرفته اند، روشهای حل عددی افزایشی لازم می باشد

برای سیستمهای سازه ای بسیار بزرگ ترکیبی از برهم نهی مودی و روشهای افزایشی می توانند بسیار مؤثر باشند. (برای سیستمهای با تعداد کمی المانهای غیرخطی)

4-1- روش برهم نهی مودی

معمول ترین و مؤثرترین رهیافت برای آنالیز لرزه ای سازه های خطی روش برهم‌نهی‌مودی می باشد. پس از آنکه گروهی از بردارهای متعامد برآورد شدند این روش دسته بزرگ معادلات تعادل را به تعداد نسبتاً کمتری از معادلات دیفرانسیل مرتبه دوم تبدیل می کند که این باعث کاهش قابل توجهی در زمان محاسبات می‌شود

نشان داده شده است که حرکات لرزه ای زمین تنها فرکانسهای پایین سازه را تحریک می نماید.به صورت معمول حرکات زلزله در فواصل زمانی 200 نقطه در ثانیه ثبت می گردند. بنا بر این داده های بارگذاری پایه شامل اطلاعات بالای 50 دور در ثانیه نمی باشند.با توجه به این مطلب صرف نظر از مودها و فرکانسهای بالاتر معمولاَ باعث ایجاد خطا نمی شوند

5-1- تحلیل طیف پاسخ

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

این روش حجم خروجی بالایی ایجاد می کند که این امر سبب زیاد شدن عملیات طراحی به خصوص هنگامی که بخواهیم نتایج را برای کنترل طراحی به کار بریم می‌گردد

تحلیل باید برای چندین زلزله دیگر هم تکرار شود تا اطمینان حاصل گرد که تمام مدها تحریک شده اند

مزایای محاسباتی قابل توجهی در استفاده از تحلیل طیف پاسخ برای پیش بینی تغییر مکانها و نیروهای اعضاء در سیستمهای سازه ای وجود دارد. این روش فقط شامل محاسبه حداکثر مقدار تغییر مکانها و نیروهای اعضاء با استفاده از طیفی هموار شده است که میانگین چندین زلزله است، می باشد. سپس لازم است برای بدست آوردن متحمل‌ترین مقدار اوج تغییر مکان یا نیرو از روشهای CQC ، SRSS و یا CQC3 استفاده  گردد

6-1- حل در حوزه فرکانس

رهیافت پایه استفاده شده در حل معادلات تعادل دینامیکی در دامنه فرکانس بسط نیروهای خارجیF(t) در قالب عبارات سری های فوریه یا انتگرالهای فوریه می باشد

حل شامل عبارات مختلط است که محدوده زمانی¥+ تا ¥- را پوشش می دهد. بنابراین روشی بسیار کارا برای گونه‌های بارهای تکرارای مانند: ارتعاشات مکانیکی، آکوستیک، امواج دریا و باد می باشد. هرچند استفاده از حل در حوزه فرکانس برای تحلیل سازه‌هایی که تحت تأثیر زلزله قرار می گیرند دارای معایب چندی نیز می باشد

فهم ریاضیات به کار رفته برای دسته زیادی از مهندسان سازه بسیار مشکل می باشد. بنابراین مطمئن شدن از صحت حل بسیار مشکل است

برای نوع بارگذاری لرزه ای  این روش از نظر عددی کارا نمی باشد. انتقال نتایج از حوزه فرکانس به حوزه زمان حتی با استفاده از روشهای FFT مقدار محاسبات عددی قابل توجهی را لازم دارد

روش محدود به سیستمهای ساختمانی خطی می باشد

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

7-1- حل معادلات خطی

حل گام به گام معادلات دینامیکی، حل در حوزه فرکانس و برآورد بردارهای ویژه و بردارهای ریتز تماماً احتیاج به حل معادلات خطی دارند که به صورت زیر بیان می‌شود

AX=B                                                                                                                (1-7-1)

که در اینجا A یک ماتریس N×N متقارن است که تعداد زیادی جمله صفر دارد. ماتریسهای B و X که
“N × M”هستند بیانگر این مطلب است که بیشتر از یک حالت بارگذاری در یک زمان قابل حل می باشد. که روشهای متعددی برای کاهش حافظه مصرفی توسط A وحل دستگاه همزمان وجود دارد. (روش حذفی گوس,حل اسکای لاین  و روشهای بسیار متنوع دیگر که برای معکوس سازی ماتریسها به کار می روند از جمله روشهای:افراز کردن,سه قطری کردن,کاهش ماتریس,روش جوردن و;)

1-2- روش جستجوی د ترمینانی (Determinant search method)

معادله تعادل که بر ارتعاش آزاد یک مد نمونه نامیرا حاکم است به صورت زیر نوشته می‌شود

      یا                                                                 (1-1-2)

این معادله را می توان با فرض i و فاکتورگیری به صورت زیر مستقیماً حل کرد

                                                                                                (2-1-2)

می توان نشان داد

                                                                                      (3-1-2)

می توان با تکرار این عمل نموداری از دترمینان در مقابلl رسم نمود. (شکل (1-1-2) این روش کلاسیک برای بدست آوردن فرکانسهای طبیعی سیستم روش جستجوی دترمینانی نام دارد

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

شکل1-1-

2-2- کنترل ترتیب استورم (Sturm sequence check)

شکل (1-1-2) خاصیت بسیار مهمی از دنباله عبارات قطری ماتریس فاکتورگیری شده را نشان می دهد. متوجه می شویم برای مقدار مشخصی از i ، می توان تعداد عبارت منفی در ماتریس قطری را شمرد که برابر تعداد فرکانسهای کمتر از آن مقدار می باشد. بنابراین، این روش می تواند برای کنترل روشی که نتوانسته تمام فرکانسهای طبیعی کمتر از مقدار مشخصی را حساب کند به کار رود

نیز کاربرد مهم دیگر این روش برآورد تعداد فرکانسهای موجود در بازه خاص فرکانسی می باشد. که این مطلب در مسائل ارتعاش ماشین کارآمد می باشد

تکرار معکوس

معادله (1-1-2) را می توان به فرمی مناسب برای روش حل تکراری نوشت  داریم

      یا                                                    (1-2-2)

گامهای محاسباتی برای محاسبه یک بردار ویژه یا مقدار ویژه به صورت زیر خلاصه می‌شود

ماتریس سختی را مثلثی می کنیم به فرم  LDLT. (در فاز حل بار استاتیکی)

برای اولین سعی فرض کنیم R(1) برداری حاوی اعداد تصادفی باشد و برای بردار اولیه  حل کنیم

برای i=1,2,… سعی می کنیم

(a  بردار را نرمال می کنیم به گونه ای که

(b مقدار ویژه را تخمین می زنیم که

(c کنترل  برای همگرایی اگر همگرا شد تمام

(d i=i+1  و محاسبه

(e حل برای بردار جدید

(f گام 3 را تکرار کنید

می توان دید که این روش به سمت کوچکترین مقدار منحصربه فرد مقدار ویژه همگرا می باشد

3-2- متعامدسازی گرام ـ اشمیت

بردارهای ویژه دیگر در روش تکرار معکوس قابل محاسبه اند به شرط آنکه بعد از هر چرخه تکرار، بردار تکرار نسبت به تمامی بردارهای محاسبه شده قبل متعامد شود. برای نشان دادن این روش فرض کنید بردار مفروض `Vموجود می باشد که می خواهیم نسبت به بردار محاسبه شده قبلی Vn متعامد شود. یا بردار جدید می تواند از رابطه زیر حساب شود

V=`V-aVn                                                                                                        (1-3-2)

اگر این معادله را در  پیش ضرب کنیم بدست می آوریم

                                                                      (2-3-2)

بنابراین شرایط تعامد در صورت برآورده شدن شرط زیر ارضا می‌شود

                                                                                       (3-3-2)

اگر این متعامد سازی بعد از گام 3e در تکرار معکوس قرار گیرد، مقادیر ویژه و بردارهای ویژه اضافی قابل محاسبه اند

4-2- تکرار زیرفضای بلوکی ((Block subspace iteration

روش تکرار معکوس با یک بردار در صورت وجود مقادیر ویژه و بردارهای ویژه مشابه ممکن است همگرا نگردد. این حالت برای بسیار از سازه های سه بعدی واقعی با جرم و سختی مشابه در هر دو جهت اصلی ممکن است اتفاق بیافتد

این مشکل را می توان با تکرار بوسیله گروهی (بلوکی) از بردارهای متعامد برطرف ساخت تجربه نشان داده است که اندازه بلوک (b) باید برابر جذر «پهنای متوسط باند ماتریس سختی» قرار داده شود ولی کمتر از 6 نگردد. این الگوریتم (روش) نسبتاً کند می‌باشد هرچند بسیار دقیق می‌باشد

در حالت کلی بعد از آنکه برداری به  بلوک اضافه شد احتیاج به 5 تا 10 کاهش به جلو و جاگذاری به عقب می باشد تا این روش به بردار ویژه دقیق همگرا شود

5-2- حل سیستمهای منفرد

برای انواع کمی از سازه ها، مانند شاتل ها، امکان ندارد که روش تکرار معکوس یا زیرفضا را به طور مستقیم برای بدست آوردن فرکانسهای طبیعی و اشکال مدی به کار برد. دلیل این امر وجود حداقل شش مد صلب با فرکانس صفر می‌باشد و ماتریس سختی منفرد است و قابل مثلثی کردن نیست. برای حل این مشکل تنها لازم است که جابجایی زیر در مقادیر ویژه انجام شود، یا تغییر متغیر بدهیم

ln=`ln-r                                                                                                            (1-5-2)

بنابراین مسئله مقادیر ویژه تکراری را می توانیم به صورت زیر بنویسیم

LDLT `Vn(i) = R(i)  یا  `K`Vn(i) = `ln(i-1)MVn(i-1)                                                             (2-5-2)

(3-5-2) ماتریس سختی جابجا شده به صورت `K=K+rM می‌باشد که دیگر منفرد نمی‌باشد

بردارهای  ویژه با  انتقال دلخواه r  دستخوش تغییر نمی شوند  بنابراین از رابطه (1-5-2)بردارهای درستی حاصل می گردد

برای دریافت اینجا کلیک کنید

سوالات و نظرات شما

برچسب ها

سایت پروژه word, دانلود پروژه word, سایت پروژه, پروژه دات کام,
Copyright © 2014 cpro.ir
 
Clicky