مقاله بررسی بردارهای ریتز وابسته به بار و روش MPA


دنلود مقاله و پروژه و پایان نامه دانشجوئی

مقاله بررسی بردارهای ریتز وابسته به بار و روش MPA   مربوطه  به صورت فایل ورد  word و قابل ویرایش می باشد و دارای ۱۶۳  صفحه است . بلافاصله بعد از پرداخت و خرید لینک دنلود مقاله بررسی بردارهای ریتز وابسته به بار و روش MPA  نمایش داده می شود، علاوه بر آن لینک مقاله مربوطه به ایمیل شما نیز ارسال می گردد

فهرست مطالب

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

 

مراجع:

 ۱- 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-464.

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

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

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-582.

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

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

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-11.

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.

مقدمه

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

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

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

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

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

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

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

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

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

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

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

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

۲-۱- تعادل دینامیکی

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

F(t)I + F(t)D + F(t)S = F(t)                                                                                  (۱-۲-۱)

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

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

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

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

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

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

                                          (۲-۲-۱)

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

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

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

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                                       (۳-۲-۱)

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

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

با قرار دادن این معادله (۳-۲-۱) در (۲-۲-۱) داریم:

(t) + C(t) + Ku(t) = -Mx ü(t)xgMy ü(t)ygMz ü(t)zg                               (۴-۲-۱)

که

Mi = M{ri}

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

۳-۱- روش حل گام به گام

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

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

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

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

۴-۱- روش برهم نهی مودی

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

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

۵-۱- تحلیل طیف پاسخ

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

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

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

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

۶-۱- حل در حوزة فرکانس

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

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

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

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

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

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

۷-۱- حل معادلات خطی

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

AX=B                                                                                                                (۱-۷-۱)

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

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

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

      یا                                                                 (۱-۱-۲)

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

                                                                                                (۲-۱-۲)

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

                                                                                      (۳-۱-۲)

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

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

شکل۱-۱-۲

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

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

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

تکرار معکوس

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

      یا                                                    (۱-۲-۲)

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

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

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

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

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

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

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

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

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

(f گام ۳ را تکرار کنید.

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

۳-۲- متعامدسازی گرام ـ اشمیت

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

V=`V-aVn                                                                                                        (۱-۳-۲)

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

                                                                      (۲-۳-۲)

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

                                                                                       (۳-۳-۲)

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

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

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

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

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

۵-۲- حل سیستمهای منفرد

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

ln=`ln-r                                                                                                            (۱-۵-۲)

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

LDLT `Vn(i) = R(i)  یا  `K`Vn(i) = `ln(i-1)MVn(i-1)                                                             (۲-۵-۲)

(۳-۵-۲) ماتریس سختی جابجا شده به صورت `K=K+rM می‌باشد که دیگر منفرد نمی‌باشد.

بردارهای  ویژه با  انتقال دلخواه r  دستخوش تغییر نمی شوند  بنابراین از رابطة (۱-۵-۲)بردارهای درستی حاصل می گردد.

120,000 ریال – خرید

تمام مقالات و پایان نامه و پروژه ها به صورت فایل دنلودی می باشند و شما به محض پرداخت آنلاین مبلغ همان لحظه قادر به دریافت فایل خواهید بود. این عملیات کاملاً خودکار بوده و توسط سیستم انجام می پذیرد.

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

مطالب پیشنهادی: برای ثبت نظر خود کلیک کنید ...

براي قرار دادن بنر خود در اين مکان کليک کنيد
به راهنمایی نیاز دارید؟ کلیک کنید


جستجو پیشرفته مقالات و پروژه

سبد خرید

  • سبد خریدتان خالی است.

دسته ها

آخرین بروز رسانی

    چهارشنبه, ۱۷ آذر , ۱۳۹۵

اولین پایگاه اینترنتی اشتراک و فروش فایلهای دیجیتال ایران
wpdesign Group طراحی و پشتیبانی سایت توسط دیجیتال ایران digitaliran.ir صورت گرفته است
تمامی حقوق برایdjkalaa.irمحفوظ می باشد.