منحنی های شکنندگی ابزارهای قدرتمندی در ارزیابی ریسک لرزه ای هستند که احتمال رسیدن یا تجاوز یک سازه (یا المان سازهای) از یک حالت آسیب مشخص را در برابر سطوح مختلف شدت زلزله (مانند شتاب طیفی) نشان می دهند. بسیاری از دانشجویان و مهندسین به دنبال یک دستور مستقیم در ETABS برای تولید این منحنی ها هستند، اما باید توجه داشت که ETABS به طور مستقیم قابلیت رسم منحنی شکنندگی را در حال حاضر ندارد, امید است در آینده با توسعه این نرم افزار امکانات لازم جهت استخراج منحنی های شکنندگی به صورت مستقیم در این نرم افزار فراهم شود. با این حال، این نرم افزار نقش اساسی در تولید داده های خام مورد نیاز برای ساخت این منحنی ها ایفا می کند. در این مقاله آموزشی، فرآیند دقیق و گام به گام استفاده از ETABS به عنوان یک نرم افزار کاربردی و شناخته شده را برای دستیابی به منحنی های شکنندگی را شرح می دهیم.
خلاصه فرآیند کلی
1. در ETABS:
آنالیز تاریخچه زمانی غیرخطی (Nonlinear Response History Analysis) را روی مدلی که المان های غیرخطی (مفاصل پلاستیک) دارد، برای یک مجموعه از رکورد های زلزله که به سطوح مختلفی مقیاس شده اند، انجام می دهیم. این روش “آنالیز دینامیکی افزایشی (IDA)” نامیده می شود.
قبلا یک در مطلب آموزشی مراحل انجام تحلیل IDA را شرح داده بودیم که می توانید به آن مراجعه کنید (لینک زیر):
مراحل انجام تحلیل IDA در Etabs | مرجع آموزش Etabs , SAP2000, SAFE
2. خارج از ETABS (در Excel، Python و غیره):
نتایج حاصل از ETABS (مانند حداکثر تغییرمکان نسبی طبقات) را استخراج کرده و با استفاده از روش های آماری، منحنی شکنندگی را رسم می کنیم.
مراحل گام به گام در ETABS
گام ۱: ایجاد و کالیبراسیون مدل غیرخطی
دقت منحنی شکنندگی شما مستقیماً به دقت مدل غیرخطی بستگی دارد.
* مدل سازه ای خود را با جزئیات کامل در ETABS بسازید.(Geometry Model)
* غیرخطی بودن مصالح (Material Nonlinearity): ویژگی های مفاصل پلاستیک (Hinge) را برای المان های بحرانی تعریف کنید.
* تیرها و ستون ها: از مفاصل خمشی (M3) استفاده کنید. می توانید از ویژگی های پیش فرض مطابق با آیین نامه هایی مانند FEMA 356 یا ASCE 41 استفاده کرده یا ویژگی های سفارشی تعریف کنید.
* دیوارهای برشی: از المان های پوسته غیرخطی (Nonlinear Shell) یا مفاصل فیبری (P-M2-M3) استفاده کنید.
* بادبندها: از مفاصل محوری (P) برای مدل کردن کمانش استفاده کنید.
* بارگذاری (Load Cases): بارهای مرده و زنده را تعریف کنید.
گام ۲: تعریف حالت بار برای آنالیز غیرخطی
برای انجام “آنالیز تاریخچه زمانی غیرخطی” یک حالت بار جدید ایجاد کنید.
* نوع آنالیز (Analysis Type): از منوی مربوطه، Nonlinear Load Case جدید ایجاد کنید.
* نوع آنالیز: گزینه *Response History* را انتخاب کنید.
* اعمال بار (Load Application): گزینه *Applied Loads* را با یک ضریب مقیاس (برای اعمال بارهای ثقلی) انتخاب کنید. سپس برای فاز غیرخطی، گزینه *Multi-Step Static* را انتخاب کنید (این بخش، زلزله را اعمال می کند).
* غیرخطی هندسی(Geometric Nonlinearity): برای ساختمان های بلند، استفاده از گزینه `P-Delta` توصیه می شود.
گام ۳: انتخاب و مقیاس کردن رکوردهای زلزله
به مجموعه ای از شتابنگاشت های واقعی نیاز دارید.
* انتخاب رکوردها (Record Selection): مجموعه ای متشکل از ۱۰ تا ۴۰ شتابنگاشت را انتخاب کنید که نماینده شرایط لرزه خیزی سایت پروژه شما باشند.
* مقیاس کردن (Scaling): این رکوردها را به سطوح مختلف شدت مقیاس کنید. این هسته اصلی روش IDA است.
* معیار شدت (Intensity Measure – IM): یک معیار شدت انتخاب کنید. رایج ترین معیار، شتاب طیفی در دوره تناوب اصلی سازه، Sa(T1, 5%) است.
* سطوح مقیاس (Scale Levels): یک بازه برای عوامل مقیاس تعیین کنید. برای مثال، می توانید هر رکورد را طوری مقیاس کنید که Sa(T1) از ۰.۱g تا ۲.۰g با گام های ۰.۱g تغییر کند. این کار صدها سناریوی آنالیز ایجاد می کند.
گام ۴: اجرای آنالیزهای متعدد (آنالیز دینامیکی افزایشی – IDA)
این بخش نیازمند محاسبات سنگین است. شما باید ETABS را بارها اجرا کنید.
1. برای هر شتابنگاشت در مجموعه خود.
2. برای هر ضریب مقیاس (سطح شدت) برای آن شتابنگاشت.
3. آنالیز تاریخچه زمانی غیرخطی را در ETABS اجرا کنید.
نکته بسیار مهم: انجام این کار به صورت دستی غیرممکن است. باید از ETABS API (رابط برنامه نویسی) با یک زبان اسکریپت نویسی مانند Python یا C استفاده کنید. یک اسکریپت می تواند به طور خودکار در بین رکوردها و عوامل مقیاس حلقه بزند، مدل ETABS را تغییر دهد، آنالیز را اجرا کرده و نتایج را استخراج کند.
گام ۵: استخراج نتایج (پارامترهای Demand)
برای هر اجرای آنالیز، باید مقادیر حداکثر پاسخ را استخراج کنید.
* پارامترهای سراسری (Global EDPs): حداکثر تغییرمکان نسبی طبقات (MIDR)، تغییرمکن بام.
* پارامترهای محلی (Local EDPs): میزان چرخش مفاصل پلاستیک، کرنش های برشی در دیوارها و غیره.
این استخراج داده نیز باید از طریق API و به صورت خودکار انجام شود.
گام ۶: توسعه منحنی شکنندگی (در خارج از ETABS)
اکنون شما یک جدول از داده های `IM` (شتاب طیفی) در مقابل `EDP` (مثلاً تغییرمکان نسبی) دارید.
۱. تعریف حالت های آسیب (Damage States – DS): پارامترهای پاسخ (EDP) را به سطوح کیفی آسیب ارتباط دهید.
* DS1 (آسیب جزئی): تغییرمکان نسبی = ۰.۷%
* DS2 (آسیب متوسط): تغییرمکان نسبی = ۱.۵%
* DS3 (آسیب شدید): تغییرمکان نسبی = ۲.۵%
* DS4 (فروریزی): تغییرمکان نسبی = ۵.۰%
(اعداد فوق مثال هستند. مقادیر دقیق را باید از منابع معتبری مانند FEMA P-58 استخراج کنید.)
۲. تعیین ظرفیت برای هر رکوردزلزله:
برای هر شتابنگاشت و هر حالت آسیب، سطح شدت (Sa) ای را پیدا کنید که در آن، پارامتر پاسخ سازه برای اولین بار از آستانه آن حالت آسیب تجاوز کرد. این کار با درون یابی بین نقاط داده حاصل از آنالیز IDA انجام می شود.
۳. انجام آنالیز آماری:
* اکنون شما یک مجموعه از مقادیر ظرفیت (یک مقدار برای هر شتابنگاشت) برای یک حالت آسیب مشخص دارید.
* لگاریتم طبیعی (`ln`) این مقادیر ظرفیت را محاسبه کنید.
* میانگین لگاریتمی (θ) و انحراف معیار لگاریتمی (β) این مقادیر `ln(IM)` را محاسبه کنید.
* `θ = exp( average( ln(IM_capacity) ) )`
* `β = standard deviation( ln(IM_capacity) )`
۴. رسم منحنی شکنندگی:
با استفاده از مقادیر میانگین (`θ`) و پراکندگی (`β`)، می توانید تابع شکنندگی را در یک نرم افزار مانند Excel رسم کنید. این تابع یک تابع توزیع تجمعی نرمال لگاریتمی است:
`P(DS | IM = x) = Φ( (ln(x) – ln(θ)) / β )`
این تابع را برای یک بازه از مقادیر `x` (مقادیر شتاب طیفی) رسم کنید. این کار را برای هر حالت آسیب تکرار کنید.
سخن پایانی
همانطور که ملاحظه کردید، ETABS به عنوان هسته محاسباتی این فرآیند عمل می کند، اما خروجی نهایی حاصل تلفیق این محاسبات با آمار و احتمالات است. توصیه می شود برای انجام پروژه های تحقیقاتی و اجرایی ، حتماً از API استفاده کنید تا بتوانید حجم عظیم آنالیزهای مورد نیاز را به صورت کارآمد و بدون خطای دستی مدیریت نمایید.
