مقدمه به عنوان قسمتی از هر برآورد لرزه ای یا طراحی لرزه ای، مهندس طراح باید تحلیلی از سازه با در نظر گیری خطر لرزه ای در محل ساختمان، برای برآورد کمیت‌های پاسخ سازه انجام دهد.
این پاسخ ها اگر در حدود پاسخ مجاز سازه قرار گیرد، قبول می گردند.
در حالت کلی، تحلیل سازه شامل اثر دادن توزیع جانبی نیروهای زلزله به علاوه نیروهای ثقلی بر یک مدل ریاضی از سازه می باشد.
روشهای تحلیل سازه با توجه به این مطلب که مدل ریاضی خطی یا غیرخطی می باشند و یا اینکه نیروها، دینامیکی و استاتیکی می باشند قابل تمایز می باشند.
فرض اصلی در یک مدل ساختمان خطی آنست که المانهای ساختمان برای مثال تیرها و ستونها، دارای قدرت نامحدود و سختی ثابت در حین تحلیل می باشند.
از طرف دیگر یک مدل غیرخطی سعی می کند کاهش مقاومت و سختی اعضا را در مدل ساختمان در هنگام خرابی در نظر بگیرد.
در تحلیل استاتیکی، یک توزیع فرضی نیروهای زلزله به ساختمان اعمال می‌گردد.
بنابراین، این نوع تحلیل ،پاسخ های متغیر با زمان سازه را که در تحلیل دینامیکی منظور می شوند، در نظر نمی گیرد.
روش های تحلیل سازه به روشهای زیر تقسیم بندی می شوند (ATC،1997)، استاتیکی خطی (LS) دینامیکی خطی (LD)، استاتیکی غیرخطی (NS)، دینامیکی غیرخطی (ND).
به خاطر فرضیات موجود در به کارگیری روش خطی فرض می شود (ND) روش دقیق است و سایر روشها، روشهای تقریبی می باشند.
تا این اواخر آیین نامه های طراحی ساختمانهای جدید بیشتر بر روی روشهای تحلیل خطی تمرکز کرده اند.
دلیل آنکه به روشهای غیرخطی کمتر پرداخته شده آنستکه روشهای غیرخطی برای امر طراحی بسیار پیچیده می باشند.
به علاوه این روشها نیازمند به کارگیری کامپیوترهای بسیار قدرتمند می باشد و نیز عدم قطعیتهایی نیز دارند.
ATC-40 بر روش NS برای برآورد لرزه ای تکیه می کند.
FEMA 273/274 روشهای NS و ND را تحت پوشش قرار می دهد.
در زمان حاضرآنالیز استاتیکی فزاینده غیرخطی، POA به عنوان یک روش کامپیوتری ساده شده برای برآورد عملکرد ساختمانها روز به روز اقبال بیشتری می یابد.
این روش فرض می کند پاسخ سازه تنها با مد اول کنترل می شود و این مد در کل تاریخچه زمانی ثابت می‌ماند.
مثالهای فراوانی نشان داده اند که حداکثر پاسخ سازه با این روش، بدون در نظر گیری مدهای بالاتر کمتر تخمین زده می شود و خطای قابل توجهی به خصوص برای سازه های با زمان تناوب بلند یا هنگامی که مکانیزمهای محلی تشکیل گردد (به علت تغییر خواص دینامیکی سازه) ایجاد می گردد.
که محققین زیادی در این موردد تحقیق نموده اند(، Fajfar & Fishinger ،Faella 1981، Saiid & Sozen 1981 ،Fajfar & Gaspersic ، Rutenberg & Destefano 1997 Fajfar & Kilar) و نتایج برخی از این تحقیقات در ادامه ارائه خواهد شد.
لازم به توضیح است کارآیی قابل قبول و سطح خطر زلزله به صورت زوج در آیین نامه های مربوطه تعریف شده‌اند.
در سالهای اخیر، روش NS نسبت به ND با توجه بیشتری روبرو بوده است.
که این لزوم بررسی صحت این روش را موجب می گردد.
2-1- پیدایش روش غیر خطی استاتیکی درسال 1975 فریمن (Freeman) وهمکارانش یک روش ارزیابی سریع بدست آوردند که درواقع مانند روش طیف ظرفیت امروزی می باشد.درسال 1981سیدوسوزن (Saiid & Sozen) یک روش آنالیز دینامیکی غیر خطی برروی یک سیستم یک درجه آزادی معادلSDOF پیشنهاد کردند.براساس این ایده فایفر وفیشینگر (Fajfar & Fischinger)در اواسط سال 1980 اولین نسخه از روش N2 بدست آوردند.(منظور ازN تحلیل غیر خطی و عدد 2 مدلهای ریاضی یک درجه ازادیFSDO وچنددرجه آزادی MDOF می‌باشد).
به هرحال جامعه مهندسی زلزله توجه زیادی به روشهای ساده غیرخطی تا اواسط سال1990 نداشتند .در روشهای ساده شده غیرخطی تحت عنوان روش طیف ظرفیت در آیین نامه های ATC 40 و ترای سرویسزTriservices، استاندارد ساختمان ژاپن ، روش آنالیز استاتیکی غیر خطی اعمال شده درFEMA273 وFEMA356، روش N2اجرا شده در پیش نویس Eurocode 8 وآنالیز Pushover مودال تفاوتهای بسیاری دیده می شود و نیزچندین روش غیر خطی در کتاب آبیSEAOC اجرا شده است.تمام روشهای موجود آنالیزPushover یک مدل چند درجه‌آزادی را با آنالیز طیف پاسخ سیستم یک درجه آزادی معادل SDOF،ترکیب می کنندوطیف الاستیک وغیر الاستیک بر اساس یک میرایی معادل وپریود اعمال می شوند که در واقع طیف غیر الاستیک نقطه جاری شدن (تسلیم)را حاصل می کند.روشهای ارائه شده در آیین نامه FEMA روشهایی هستند که براساس فرمت شتاب– تغییر مکان (A-D) فرمول بندی شده اند.در این شکل ظرفیت سازه به طور مستقیم با تقاضای حرکت زمین در زلزله روی یک سازه مقایسه می شود.منحنی نیرو-تغییرمکان که توسط آنالیز استاتیکی غیر خطی به دست آمده است ، ظرفیت سازه را نشان می دهد.
نیروهای برش پایه وتغییر مکانهای بام به طیفهای شتاب وطیفهای تغییر مکان یک سیستم یک درجه آزادی معادل SDOF برگردانده می شوند که مقادیر این طیفها نمودار ظرفیت را تعریف می کند.مقادیر طیف لرزه ای تقاضا در بین روشهای متفاوت دارای محتوای مختلفی است .درتمام حالات، تقاطع منحنی ظرفیت وطیف تقاضا یک تخمینی ازشتاب غیر الاستیک(مقاومت) وتغییر مکان تقاضا را بدست می دهد.بنابراین بسیاری از روشهای ساده شده به سمت طراحی براساس کنترل تغییرشکل (که طراحی براساس تغیرمکان نامیده میشود)سوق یافته ودرحال گسترش می باشند.تمام روشها براساس سازه های صفحه ای محدود می شوند که اخیرا فرضیاتی جهت توسعه وکاربردی کردن روشهای موجود برای سازه های نامتقارن با استفاده ازتحلیلهای سه بعدی3D انجام گرفته است وروشهای تحلیل استاتیکی غیر خطی جهت تعیین وتخمین پاسخ لرزه ای سازه های نامتقارن توسعه داده شده است.روشهای پیشنهاد شده در واقع شامل پاسخهای پیچشی ناشی از تاثیرات آنالیزهای سه بعدی این ساختمانهای نامتقارن می باشند که درفرم جدید این روشها با استفاده از آنالیزهای دینامیکی الاستیک ساختمانها ،جهت تعیین تغییر مکانهای هدف وشکل توزیع بار جانبی برای آنالیز استاتیکی غیر خطی استفاده می کنند.
3-1- فرضیات اساسی تمام آنالیزها توسط برنامه هایی تحت عنوان برنامه های آنالیز اجزای محدود FEAP(Finite Element Analysis Program)به صورت استاتیکی ودینامیکی غیر خطی انجام پذیر است.آنالیزهای فوق قابل تعمیم از حالت دوبعدی به حالت سه بعدی می باشند.نتایج آنالیز نسبت نیرو - تغییر شکل را برای تمام مقاطع در طول المان ونسبت تنش - کرنش را برای تمام نقاط مقطع المان نشان میدهد.فرضیات بنیادین در آنالیز استاتیکی غیر خطی فزاینده شامل موارد زیر می باشد (بحث تکمیلی موارد 1 تا 4 در بخشهای بعدی ارایه خواهد شد): 1-3-1- کنترل بر اساس نیرو یا تغییرمکان Displacment Control Force vs.
آنالیز مورد نظر قابل انجام در دوحالت زیر می باشد: الف- یک توزیع تغییر مکان افقی ب- یک توزیع نیروی افقی درحالتی که یک نوع توزیع تغییر مکان افقی وجود دارد، سازه متحمل یک نوع تغییر فرم از پیش تعیین شده می گردد که باید تفاوت مناسبی با فرم طبیعی سازه داشته باشد،به خصوص زمانی که اجزای سازه منجمله قابها وارد ناحیه پلاستیک می شود.
2-3-1- الگوهای بارگذاری Pattern Shape Load بارهای جانبی اعمال شده به هر کف در هنگام آنالیز استاتیکی غیرخطی می تواند تاثیر قابل ملاحظه ای در پیش بینی سیستم و نیروهای المانها داشته باشد.
این الگوی بارگذاری جانبی باید به گونه ای باشد که نشان دهنده نیروهای اینرسی موجود بر سازه در هنگام زلزله باشد.
الگوهای بارگذاری جانبی به دو دسته ثابت و متغیر تقسیم می‌گردند.
در مورد الگوی بارگذاری جانبی ثابت، فرض می شود نیروهای اینرسی در هنگام پاسخ سازه به زلزله ثابت می مانند.
این نوع بارگذاری قابلیت به حساب آوردن تغییرات توزیع نیروی جانبی سازه را به علت تغییر سختی، ناشی از رفتار غیر خطی را ندارد.
از آنجایی که توزیع رفتار غیرخطی بسته به الگوی بار ممکن است تغییر کند، پیشنهاد شده است که چند توزیع در هنگام انجام آنالیز در نظر گرفته شود.
برای بارگذاری می توان روشهای مختلفی استفاده نمود سه نوع روش معمول عبارتند از: اشکال مختلف توزیع بار الف- شکل توزیع نیرو به صورت مثلثی معکوس ب- شکل توزیع یکنواخت ج- شکل توزیع نیرو با در نظرگیری مدهای بالاتر با توجه به نتایج حاصل از آنالیزهای مختلف استاتیکی غیر خطی افزاینده مشخص شده است این بارگذاریها برای انواع خاصی از سازه ها جوابهای مناسبی بدست می دهند و قابل کاربرد برای تمامی انواع سازه ها نمی باشند.
همانطوریکه اشاره شد برای بارگذاری حالات گوناگون و متعددی وجود دارد و بنا بر فلسفه اولیه تحلیل فزاینده استاتیکی غیرخطی (در نظر گیری مد اول) باید برای رسیدن به درک درستی از رفتار لرزه ای سازه حداقل 2 توزیع متفاوت در نظر گرفته شود برای مثال مطابق آیین نامه FEMA 356 داریم: I) الگوی مدی الف) در صورتی که مشارکت جرمی در مد اول بیشتر از 75% باشد از توزیع Cvx به صورت زیر می توان استفاده نمود.
Kبین این دو عدد به صورت خطی تغییر می کند.
ب) یک توزیع قائم متناسب با شکل مد پایه (مد اول) در جهت مورد بررسی (تنها در هنگامی که مشارکت جرمی بیشتر از 75% باشد).
ج) یک توزیع قائم با در نظر گیری ترکیب برش مدی حاصل از تحلیل طیف پاسخ که تعداد مد به گونه ای انتخاب شده باشد که حداقل 90% مشارکت جرمی حاصل شده باشد.
این توزیع زمانی استفاده می گردد که تناوب مد اول از Sec1 بیشتر باشد.
II) الگوی دوم (یکی از دو صورت زیر استفاده شود) الف)یک توزیع یکنواخت که شامل نیرویی متناسب با جرم در هر تراز می باشد.
ب) توزیع سازگار شونده که با توجه به تغییر شکل سازه تغییر می یابد.
این توزیع بار با در نظرگیری خواص تسلیم سازه باید اصلاح شود.
انواع توزیع بار پیشنهادی: الگوی A: الگوی B: الگوی C:(الگوی مدی) الگوی D: توضیح الگوی D: بعد از انجام آنالیز مد اول برش طبقات به واسطه روش SRSS قابل محاسبه می‌باشد اگر فرض کنیم Fij و Qij به عنوان نیروی جانبی و برش طبقه i در مد jام باشند و Qi برش طبقه i حاصل از SRSS، N مد می‌توانیم برای محاسبه Pi (نیروی جانبی معادل) روابط عنوان شده را به کار بریم.
ضریب تاثیر افقی زلزله طبق آیین نامه چین می باشد.
ضریب مشارکت مدی می باشد.
Xij تغییر مکان نسبی افقی طبقه نامتناظر با مد j است.
Wi نیز وزن طبقه iام می باشد.
بنا بر مطالعات انجام شده این روش نسبت به روشهای قبلی جوابهای بهتری ارایه می دهد.
لازم به توضیح است علی رغم تمامی تلاشهای صورت گرفته و بهبودهای نسبی هنوز هم مشکلات این روش به جای خود باقی می باشند.
اخیراً توسط Chopra و Goel روشی ارائه شده است که به نام MPA یا Modal Pushover Analysisشناخته میشود در ادامه به آن خواهیم پرداخت.
3-3-1- تبدیل سازه چنددرجه آزادیMDOF به سازه یک درجه آزادی معادلSDOF تبدیل سازه چنددرجه آزادی به سازه یک درجه آزادی معادل فرمول بندی سیستم تک درجه آزادی منحصر به فرد نمی باشد اما روشهای زیادی بر مبنای بردار شکل و معادله تعادل دینامیکی شکل گرفته اند.
معادله دیفرانسیل حاکم بر سازه چند درجه آزادی عبارتست از: (1) که [M] ماتریس جرم، [C] ماتریس میرایی، {X} بردار تغییر مکان نسبی، {Q} نیروی طبقه و تاریخچه شتاب حرکت زمین می باشد.
اگر بردار شکل {} را نسبت به جابجایی بام xt نرمال کنیم داریم.
(2) با جاگذاری در معادله 1 (3) اگر تغییر مکان مرجع تک درجه آزادی را xr بنامیم.
(4) با پیش ضرب معادله (3) را در و با جایگزینی xt از معادله (4) معادله حاکم عبارت خواهد بوداز: (5) که داریم.
[FEMA (1998)] (6) (7) (8) تناوب اولیه سیستم تک درجه آزادی (Teq) به صورت زیر قابل محاسبه می باشد.
(9) که Kr سختی الاستیک سیستم تک درجه آزادی معادل می باشد.
ارتباط نیرو- تغییر مکان در سیستم تک درجه آزادی معادل از نتایج NSP یک ساختمان چند درجه آزادی با توجه به بردار شکل محاسبه شده در بالا قابل محاسبه می باشد.
برای تشخیص قدرت و تغییر مکان کلی یک رابطه چند خطی به صورت رابطه ای دو خطی ارائه می گردد.
به عنوان راهنما روابط میان نیرو در مقابل تغییر مکان برای سیستم MDF و سیستم معادل SDF در شکل ارائه شده اند.
ایده آل سازی دوخطی FEMA 273 FEMA 356 نیروی برش پایه در تسلیم (Vy) و تغییر مکان معادل بام (xt,y) از شکل فوق به همراه معادلات 3و4 به منظور محاسبه رابطه نیرو- تغییر مکان سیستم SDF معادل به صورت زیر به کار می روند.
تناوب اولیه سیستم تک درجه آزادی معادل (Teq) به صورت زیر بیان می شود.
(10) در جایی که برای سیستم SDF مرجع به صورت زیر محاسبه می گردد.
(11) و نیروی تسلیم سیستم مرجع SDF به صورت زیر محاسبه می گردد.
(12) که {Qy} نیروی طبقه در هنگام تسلیم می باشد.
، برابر در سیستم چند درجه آزادی فرض می شود.
4-3-1- تغییر مکان هدف همانطوریکه عنوان شد صحت روش POA وابستگی زیادی به توزیع بار و تغییر مکان هدف معرفی شده دارد در این قسمت به روشهای پیشنهادی برای بدست آوردن تغییر مکان هدف می پردازیم.
با استفاده از اطلاعات بالا می توان تغییر مکان هدف را بدست آورد.
در حقیقت استفاده از خصوصیات سیستم تک درجه آزادی به همراه اطلاعات طیف مربوط به سیستمهای یک درجه آزادی غیرالاستیک، اطلاعات لازم برای تخمین تغییر مکان هدف را بدست می دهد.
برای سیستمهای الاستیک تغییر مکان طیفی از رابطه زیر قابل محاسبه می باشد.
اما تغییر مکان سیستمهای غیرخطی متفاوت است (شکل زیر).
در محدوده تناوب کم، نسبت تغییر مکان غیرالاستیک به الاستیک قویاً به نیاز تغییر مکان غیرالاستیک سیستم بستگی دارد.
رابطه بین R و C1 رابطه میان نسبت شکل پذیری و نسبت تقاضای قدرت الاستیک به غیرالاستیک توسط روابط ارائه شده توسط میراندا و برترو 1994 قابل بیان می باشد.
بنابراین برای محاسبه تغییر مکان هدف تقاضای شکل پذیری سیستم SDF باید محاسبه گردد.
داریم: (14) رابطه 14 بیانگر پاسخ یک سیستم SDF با جرم واحد و تناوب Teq و قدرت تسلیم Fy,eq که به صورت زیر بیان می گردد، می باشد.
(15) اگر طیف الاستیک معلوم باشد داریم: (16) پس می توانیم فاکتور کاهش مقاومت R را به صورت زیر محاسبه کنیم.
(17) حال می توانیم با توجه به روابط منتشر شده میان نیاز شکل پذیری سیستم را بدست آوریم.
به این نکته توجه شود که روابط منتشر شده حاوی مقادیر میانگین می باشند و برای بدست آوردن برای سیستمهای مهم بهتر است از میانگین به علاوه یک انحراف معیار استاندارد استفاده نماییم.
به علت آنکه برای سیستم SDF و MDF یکسان فرض می شود مقدار تغییر مکان هدف MDF برابر خواهد بود با (xt,t): ممکن است اصلاحات بیشتری نیز به علت در نظر گیری اثر خاک محلی، تاثیرات کاهش سختی و مقاومت، تاثیرات مرتبه دوم و… برای تغییر مکان هدف لازم باشد.
دو کمیت کلیدی برای محاسبه تغییر مکان هدف عبارتند از Fy,Teq این کمیات وابسته به بردار شکل و نیرو و توزیع جرم در ارتفاع می باشند.
اما مطالعات نشان داده است که فرضیات زیر صحت مناسبی دارند.
Teq=T1 صحت این فرضیات برای سازه ای با توزیع بار مثلثی و جرم مساوی طبقات و بردارهای شکل متفاوت نشان داده شده است.
این مطالعات نشان داده اند که Teq و Fy,eq به انتخاب بردار شکل حساس نمی باشد.
بردارهای شکل استفاده شده در مطالعه حساسیت حساسیت به بردار شکل انتخاب شده در روش دیگر منحنیهای ظرفیت حاصل از سیستمهای یک درجه آزادی را می توان توسط منحنیهای دو خطی تقریب زد .منحنی ظرفیت دوخطی سیستم یک درجه آزادی توسط یک خط که در تراز85.0برابر نیروی ماکزیمم به صورت افقی رسم شده قابل محاسبه خواهد بود.
در مجموع چون در اینجا هدف ما بررسی این مطلب نمی باشد لذا برای سایرروشهامی‌توان به دستورالعمل بهسازی لرزه ای و تفسیر آن، همچنین به استاندارد FEMA273 و تفسیر آن مراجعه نمود.
5-3-1- حداکثر شتاب زمین بر اساس منحنی ظرفیت دوخطی حاصل از پاسخ سازه (تغییر مکان – برش پایه)،GAP قابل ارزیابی خواهدبود.
که دررابطه فوق Fmax حداکثر برش پایه سازه ، eff m جرم موثر سازه ، شکل پذیری ارزیابی شده براساس منحنی دوخطی حاصل از اولین مکانیزم شکست (طبقه نرم،کاهش ظرفیت شکل پذیر در ستون ، شکست برشی ستون یا شکست اتصال تیر به ستون )، (T) a S شتاب طیفی مطابق با پریود بنیادیT تعیین شده براساس سختی الاستیک اولیه .
شتاب ظرفیتCapacity Accelerationدر رابطه فوق باید با شتاب تقاضاDemand Accelerationتعیین شده براساس مطالعات خطر(تحلیل ریسک سایت مورد نظر)Hazard Analysis درمحل مقایسه شود.
درکل با توجه به مطالب عنوان شده می توان گفت POA پایه تئوری سختی ندارد فرضیات اصلی آن عبارتند از: 1) پاسخ سازه را می توان به پاسخ سیستم یک درجه آزادی معادل مربوط نمود.
که این بدان مفهوم است که پاسخ توسط یک مد کنترل می شود.
2) بردار شکل در کل پاسخ تاریخچه زمانی ثابت می ماند.
به وضوح هر دو فرض ذکر شده نادرست می باشند.
اما مطالعاتی که توسط چندین محقق انجام شده است نشان می دهد که این فرضیات،ما را به سوی تخمین خوبی از پاسخ حداکثر لرزه ای سیستم چند درجه آزادی راهنمایی می کنند.
4-1- روش آنالیز استاتیکی غیرخطی انجام یک تحلیل استاتیکی غیرخطی نیازمند آنست که یک مدل ریاضی غیرخطی تحت بارهای افزاینده جانبی قرار گیرد تا به یک «تغییر مکان هدف» برسیم.
که این تغییر مکان حداکثر تغییر مکانی است که انتظار داریم در زمان زلزله در بام ساختمان رخ دهد.
بار جانبی به صورت افزاینده اعمال می شوند و بعد از هر افزایش بار، سختی و مقاومت المانهای سازه به روز آوری می گردند.
برش پایه در مقابل تغییر مکان بام که به نام (منحنی ظرفیت) شناخته می شود، نتیجه اصلی تحلیل NS می باشد زیرا نشان دهنده کارایی کلی ساختمان می باشد.
از تحلیل فزاینده استاتیکی غیرخطی، نیازهای نیرویی و تغییر مکانی در هنگام تغییر مکان هدف برآورد می شوند و در مقابل نیرو و تغییر مکان مجاز سنجیده می شوند.
5-1- روش گام به گام در محاسبه منحنی ظرفیت با استفاده از آنالیز فزاینده استاتیکی غیرخطی قبل از آنکه روش گام به گام را ذکر کنیم لازم به توضیح است که ما در طراحی براساس عملکرد با دو عنصر کلیدی ظرفیت و تقاضا سر و کار داریم.
ظرفیت قابلیت ایستادگی سازه در مقابل تقاضای لرزه ای می‌باشد.
روش های ساده شده ای که از تحلیل استاتیکی فزاینده غیرخطی استفاده می‌کنند مانند روش طیف ظرفیت، و روش ضریب تغییر مکان نیازمند مشخص بودن سه عنصر اصلی می باشند.
1) ظرفیت 2) تقاضا 3) عملکرد که در این قسمت هدف ما بدست آوردن ظرفیت با روش آنالیز استاتیکی فزاینده غیرخطی می باشد.
1-5-1- روش گام به گام محاسبه منحنی ظرفیت ظرفیت سازه توسط یک منحنی فزاینده بدست می آید.
مناسب ترین روش برای رسم نمودار نیرو- تغییر مکان ردیابی برش پایه در مقابل جابجایی بام می باشد.
بعضی برنامه های کامپیوتری مانند DRAIN-2DX و DRAIN-3D-X و SAP2000 و ETABS2000 قابلیت انجام مستقیم این تحلیل را دارند اما وقتی از برنامه ای خطی استفاده می کنیم ممکن است روندی مثل روند زیر در نظر گرفته شود.
توجه: منحنی ظرفیت معمولاً به گونه ای ساخته می شود که نشانگر پاسخ مد اول سازه باشد (با این فرض که مد پایه پاسخ غالب سازه می باشد).
این در حالت کلی برای سازه هایی که تناوب پایه تا حدود Sec1 دارند درست می باشد.
در سازه های با تناوب طبیعی بزرگتر از Sec1 باید اثر مدهای بالاتر هم دیده شود.
1) با توجه به قواعد مدلسازی مدلی از سازه (با پی یا بدون پی) بسازید.تمامی اجزایی که دارای تاثیر قابل توجه در سختی ، مقاومت،وزن و پایداری می باشند باید در مدل منظور شوند.
2) المانهای موجود در مدل را دسته بندی کنید (اصلی و فرعی).
3) بار جانبی را به سازه اعمال نمایید.
لازم است بارهای ثقلی نیز در نظر گرفته شوند.
روش استاتیکی غیر خطی به فرمهای گوناگونی با متدولوژی های گوناگون ارائه شده است (Seneviranta & Krawinkler 1994,Moehle 1992).
همانگونه که از نام آن بر می آید روندی است که در طی آن سازه با حالت از پیش تعریف شده بارگذاری به صورت افقی هول داده می شود، تا به تدریج سازه به وضعیتی حدی برسد.
تنوع زیادی در روش های بارگذاری جانبی وجود دارد که در زیر به تعدادی از آنها اشاره می گردد.
الف) تنها یک بار متمرکز افقی به بالای سازه اعمال شود (معمولاً برای سازه‌های یک درجه آزادی) ب) بارگذاری بر روی هر طبقه با توجه به روند موجود در آیین نامه ها انجام گردد.
یا ج) اعمال بارهای جانبی به صورت سازگار با حالت مد اول به صورت حاصل ضرب جرم و شکل مدی مد اول الاستیک سازه .
معمولاً منحنی ظرفیت به گونه ای بدست می آید که اثر این مد (مد اول) را ارائه نماید.
این در حالت کلی برای سازه های با تناوب پایه کمتر از Sec1 قابل اعمال است.
(روش پایه) د) در این حالت مانند حالت (ج) عمل می کنیم برای هر افزایشی در ورای تسلیم، نیرو به گونه ای تنظیم شود که سازگار با شکل تغییر یافته باشد.
(معمولاً برای سازه‌های با طبقات ضعیف) هـ) مانند (ج‌ود) اما به منظور نمودن تاثیر مدهای بالاتر (بارها را با توجه به اشکال مدی مدهای بالاتر بر سازه اعمال می کنیم) برای مدهای بالاتر، سازه همزمان هول داده می شود و یا کشیده می گردد.
لازم است حداقل دو توزیع جانبی برای برآورد سازه به کار گرفته شود.
4) محاسبات نیروهای اجزا برای ترکیبات مورد نیاز بارهای قایم و جانبی.
5) تنظیم نیروی جانبی که گونه ای که تنش داخل المان (یا گروه المانها) تا حدود 10% مقاومتش برسد.
6) ثبت برش پایه و جابجایی بام.
(بهتر است نیروهای المانها و چرخشها نیز ثبت شوند تا در کنترل عملکرد استفاده گردند).
7) بازنگری در مدل با قرار دادن سختی صفر(یا بسیار کوچک) برای المانهای تسلیم شده.
8) اعمال افزایش جدید در بار جانبی و بازنگری در سازه تا آنکه المان یا گروه المان دیگری تسلیم شوند.
9) افزایش تغییر مکان بام و برش پایه را به مقادیر متناظر قبلی اضافه نمایید تا مقدار تغییر مکان بام و برش پایه جدید را بدست آورید.
10) گامهای 7و8و9 را تا رسیدن سازه به حد تغییر مکان نهایی ادامه دهید.
مانند ناپایداری به علت اثرات .
اغتشاش قابل توجه فراتر از سطح عملکردی مورد نظر، المان یاگروهی از المانها به سطحی از تغییر مکان جانبی برسند که در آن کاهش قابل توجه قدرت شروع شود یا المان یا گروهی از المانها به حدی از تغییر مکان جانبی برسند که قابلیت حمل بارهای ثقلی را از دست بدهند.
استثناء: در وضعیت های ویژه ای در جایی که المانها تمام یا قسمت عمده ای از قابلیت باربری جانبی خود را از دست می دهند اما می توانند بدون هیچ گونه اثرات غیرقابل قبول دیگری به تغییر شکل ادامه دهند ادامه تحلیل ممکن است موجه باشد.
برای مثال می توانیم به تیرهای رابط یا هم بند اشاره نماییم.از بین رفتن این تیرها باعث باز توزیع نیروی جانبی می شود برای مدلسازی این مطلب از مدلسازی صریح استفاده می گردد.
که باعث رسیدن به شکل دندانه اره‌ای می‌گردد.
6-1- محدودیتهای POA گامهای اساسی در POA عبارتند از: 1) فرض رابطه غیرخطی نیرو- تغییر مکان برای المانهای سازه (مقاومت تسلیم، سخت شدگی مجدد، کاهش سختی و…) 2) محاسبه تغییر مکان هدف سازه 3) انتخاب بارگذاری جانبی معقول و هل دادن سازه تحت این بارگذاری که به صورت یکنواخت در هر گام افزایش می یابد.
اگر یک عضو سازه‌ای تسلیم شود، سختی آن اصلاح می گردد تا زمانی که به تغییر مکان هدف برسیم یا سازه فرو بریزد.
در این لحظه برآورد لرزه ای سازه انجام شده است.
از گامهای بالا مشخص می شود که تغییر مکان هدف و بارگذاری جانبی برای برآورد کارایی لرزه‌ای سازه ها توسط POA بسیار مهم و حیاتی می باشد.
فرض میشود که تغییر مکان هدف نشانگر تغییر مکان حداکثری است که سازه در اثر زلزله تجربه خواهد کرد و الگوی بارگذاری جانبی نشان دهنده توزیع نیروی اینرسی در زلزله می باشد.
رفتار واقعی سازه هر زمانی که تحت تاثیر جنبشی نیرومند قرار گیرد توسط انجام تحلیل دینامیکی غیرخطی قابل تخمین زدن می باشد.
نیز می توان ترتیب ترک و یا تسلیم اعضا را نیز مشخص نمود و در نتیجه رفتار سازه را درک کرد.
بنابراین تحلیل دینامیکی غیرخطی به عنوان روشی قوی و منطقی برای طراحی لرزه ای سازه ها فرض می شود.
در حال حاضر بسیاری از سازه های ویژه، پیچیده و مهم توسط روشهای دینامیکی غیرخطی تحلیل می شوند.
با این وجود، این روش بسیار وقت گیر می باشد و نیز عدم اطمینان های زیادی موجود می باشند (برای مثال لزوم در نظرگیری چندین رکورد زلزله، رفتار هیترسیس اعضا و…).
بنابراین روش تحلیل دینامیکی غیرخطی به طور معمول در مطالعات تئوری استفاده می شود و در طراحی های معمول به کار نمی رود.
پس کار ارزنده ای می باشد که روش ساده و کامپیوتری برای برآورد لرزه ای سازه ها پیشنهاد شود.
برای این منظور POA معرفی شده است.
برای بسیاری از موارد، می توانیم اطلاعات ارزنده تری از POA نسبت به تحلیل دینامیکی غیر خطی بدست آوریم، این روش، روش ساده و اقتصادی می باشد و امروزه توجه بیشتری را به خود جلب نموده است.
بخش دوم: روش MPA: Modal Pushover Analysis روند استاتیکی غیرخطی NSP یا POA همانطوریکه در FEMA273 و FEMA356 بیان شده است امروزه به صورت گسترده ای در میان مهندسان سازه استفاده می گردد.
در طی سالیان گذشته تحقیقات زیادی بر روی فرضیات و محدودیتهای این روش انجام شده است.
برای مثال بارگذاری مطابقت کننده که سعی می کند توزیع وابسته به زمان نیروهای اینرسی را مدل کند و نیز در نظر گیری مدهایی به جز مد اصلی نوسان از جمله راه‌حلها بوده است.
با توجه به اصول دینامیک سازه روش MPA پیشنهاد شده است که امکان در نظرگیری تمامی مدهایی که تاثیر قابل ملاحظه ای بر پاسخ دارند را به استفاده کننده میدهد.
این روش همان سادگی روندهای استاندارد تحلیل فزاینده استاتیکی غیرخطی را حفظ می کند و در ضمن توزیع بار نیز مستقل از زمان می باشد.
از سال 1997، محققین متعددی روش POA را برای ساختمانهای غیر متقارن توسعه دادند.
با اعمال توزیع متناسب با ارتفاع (height-wise) نیروی جانبی و انجام POA صفحه ای بر مرکز جرم طبقات یک روش تقریبی بسط یافت [18](kilar،fajfar) اما نویسندگان این گونه عنوان نمودند که «به نظر نمی رسد این روش صحت زیادی داشته باشد».
روش دیگری که شامل (i) یک آنالیز سه بعدی طیف پاسخ الاستیک برای بدست آوردن تغییر مکان بام و توزیع نیروی ارتفاعی برای المان های مقاوم و (ii) تحلیل صفحه‌ای برای هر المان مقام بود نیز توسط (Moghadam، TSO) ارائه شد.
مطالعات چندی نیز بر روی تاثیر اندرکنش قاب- دیوار در سازه های قاب دیوار انجام شد (Dstefano،Rutenberg).
مقالاتی نیز در مورد اعمال نیروی جانبی در محلهای گوناگون پلان ارائه شد (Faella، Kilar).
مقایسات نسبتاً کمی که بین‌روش POA و RHA (Response History Analysis) انجام شد نشان دهنده موفقیتهای محدودی بودند.
بنابراین ارائه روشی که جوابهای تقریبی نسبتاً قابل قبولی ارائه دهد یک نیاز برای مهندسین می‌باشد.
در اینجا هدف اصلی ما بسط MPA برای تخمین تقاضای لرزه ای ساختمانهای متقارن و غیرمتقارن در پلان بوده است.
برای بیان MPA ابتدا لازم است مفهوم جدیدی به عنوان (UMRHA) Uncoupled Modal Response History Analysis بسط داده شود، که نشان داده می شود برای سیستمهای الاستیک خطی معادل (RHA) می باشد و برای سیستمهای غیرخطی روشی تقریبی می باشد.
سپس MPA ارائه می گردد و نشان داده می شود که برای سیستمهای الاستیک معادل تحلیل طیف پاسخ می باشد(RSA) و فرضیات آن برای سیستمهای غیرخطی بررسی می گردد.
1-2- معادلات حرکت ترکیبی از قابهای مقاوم در مقابل نیروی جانبی را در نظر بگیرید.
در شکل (1) پلان نشان داده شده حول x و یا حول y متقارن نمی باشد این بدان معناست که توزیع جرم و یا سختی متقارن نمی باشند.
با فرض کف صلب می توانیم سه درجه آزادی را در مرکز جرم CM تعریف نماییم.
شکل1- (الف) پلان (ب) قابها معادله حرکت بردار u دارای 1*3N مولفه می باشد (سه مولفه 1*N) یعنیux،uy، معادله حرکت به صورت زیر می باشد.
(1) و ماتریس M یک ماتریس قطری می باشد که دارای مرتبه 3N می باشد و سه زیر ماتریس قطری IO،m،m (ممان اینرسی قطبی) دارا می باشد.
بردارهای تاثیر به صورت زیر می باشند.
(2) طرف راست معادله را میتوان به عنوان نیروی موثر زلزله درنظر گرفت.
(3) 2-2- معرفی سیستمهای مورد بررسی وحرکت زمین سیستمهای در نظر گرفته شده گونه های مختلف ساختمان فلزی 9 طبقه می باشد که بیانگر سازه ای با ارتفاع متوسط است و مطابق سازه های در نظر گرفته شده برای پروژه فولاد SACمی باشد.
شکل 2- پلان انتخاب شده برای ساختمانهای نامتقارن در پلان مقدار خروج از مرکز برای نامتقارن ساختن سازه حول محور y برابر 10% بعد انتخاب شد.
با اعمال ضریب بر نسبت Ioj/m در سیستم متقارن سه سیستم زیر حاصل شدند: U1 : با اعمال ضریب 1 سیستمی حاصل شد که از نظر پیچشی سخت بود.که آنرا سیستم Torsionally-Stiff می نامیم.
در این سیستم حرکات انتقالی در مد اول وحرکات پیچشی در مد دوم غالب بودند.
U2 : با اعمال ضریب 2.95 سیستمی حاصل شد که دارای تناوبهای مودال نزدیک به هم بود که آنرا سیستم Torsionally-Similarly-Stiff می نامیم.در این سیستم حرکات انتقالی و پیچشی به شدت به یکدیگر وابسته اند.
U3 : با اعمال ضریب 6 سیستمی حاصل شد که از نظر پیچشی انعطاف پذیر بود.که آنرا سیستم Torsionally-Flexible می نامیم.
در این سیستم حرکات پیچشی در مد اول وحرکات انتقالی در مد دوم غالب بودند.
شکل 3-مدها و تناوبهای طبیعی برای سیستمهای انتخاب شده (a) U1 (b) U2 (c) U3 حرکت زمین انتخاب شده حرکت زمین انتخاب شده LA25 می باشد که در شکل 4 نشان داده شده است.این یکی از 20 حرکت زمینی است که برای پروژه فلزیSAC تشکیل شده است و دارای تناوب بازگشت 2475 سال است.این حرکت مشابه حرکات حوزه نزدیک می باشد و از حرکات ثبت شده در رینالدی به هنگام زلزله Northridge مشتق شده است.
شکل 4- حرکت زمین مفروض 3-2- روند تقریبی تحلیل 1-3-2- بسط مدی نیروهای موثر توزیع مکانی نیروهای موثر ( معادله 3 ) به صورت زیر نوشته میشود: (4) که ، nامین مد ارتعاش طبیعی سیستم می باشد که سه زیر بردار دارد.
(5) بنابراین می توانیم بنویسیم: (6) و مشارکت مد nام در Peff (7) که (8) هرچند به مولفه حرکت زمین وابسته می باشد معادله (5).در شکل زیر یک بسط مدی نشان داده شده است.
شکل 5- نمونه ای از بسط مدی برای سیستم U2: 2-3-2- ایده اساسی دو روش برای تحلیل تقریبی ساختمانهای غیر الاستیک بررسی خواهند شد.
تحلیل تاریخچه زمانی غیر درگیر (UMRHA) و MPA.
در روش UMRHA تاریخچه زمانی پاسخ یک سازه به Peff،n(t)، nامین مد ارتعاش توسط یک RHA غیرخطی بر روی یک سازه یک درجه آزادی بدست می آید و برهم نهی این پاسخ های مودال، کل پاسخ را بدست می دهد.
در روش MPA پاسخ حداکثر به Peff،n(t) توسط یک آنالیز استاتیکی غیرخطی محاسبه می شود و حداکثر پاسخ های مودال توسط قوانین ترکیب مدی ترکیب می گردند تا پاسخ نهایی حاصل شود.
4-2- روش UMRHA 1-4-1- سیستمهای خطی با توجه به آنالیز مودال کلاسیک داریم.