راهنمای جامع انجام پروژه با OpenSees برای تحلیل لرزه‌ای سازه‌ها

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

۱. چرا OpenSees انتخاب اول پژوهشگران مهندسی زلزله است؟

OpenSees به دلیل ساختار متن‌باز، کتابخانه غنی از مصالح غیرخطی پیشرفته و توانایی بالا در تحلیل‌های دینامیکی تاریخچه زمانی غیرخطی (NDTHA)، استاندارد طلایی در پژوهش‌های مهندسی زلزله است. این نرم‌افزار به محققان اجازه می‌دهد تا فراتر از محدودیت‌های نرم‌افزارهای تجاری معمولی، رفتار پسماند دقیق مصالح و مفاصل پلاستیک را شبیه‌سازی کنند.

بر خلاف نرم‌افزارهای گرافیکی مانند SAP2000 یا ETABS، نرم‌افزار OpenSees به شما امکان می‌دهد کنترل کاملی روی الگوریتم‌های حل، معیارهای همگرایی و فرآیند ماتریس سختی سازه داشته باشید. این ویژگی به ویژه در پایان‌نامه‌هایی که به دنبال ارزیابی آسیب، فروریزش سازه‌ها، یا اندرکنش خاک و سازه (SSI) هستند، اهمیت دوچندانی پیدا می‌کند.

۲. مراحل گام‌به‌گام انجام پروژه در OpenSees

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

در ادامه، مراحل اصلی که ساختار اصلی فایل ورودی شما را تشکیل می‌دهند بررسی می‌کنیم:

۱. تعریف هندسه و گره‌ها (Nodes)

ابتدا باید سیستم مختصات و گره‌های سازه را تعریف کنید. در انتخاب واحدها بسیار دقیق عمل کنید؛ چرا که این نرم‌افزار فاقد سیستم خودکار کنترل واحد است. پیشنهاد می‌شود از سیستم استاندارد (نیوتن، متر، کیلوگرم) یا (کیلونیوتن، میلی‌متر، تن) استفاده کنید.

# Define model builder: 2D with 3 Degrees of Freedom per node
model BasicBuilder -ndm 2 -ndf 3

# Define nodes: node tag, X-coordinate, Y-coordinate
node 1   0.0   0.0
node 2   4.0   0.0
node 3   0.0   3.2
node 4   4.0   3.2

۲. تعریف مصالح غیرخطی (Materials)

برای مدل‌سازی رفتار بتن و فولاد تحت بارهای چرخه‌ای زلزله، از مصالح یک‌محوره (uniaxialMaterial) استفاده می‌شود. پرکاربردترین آن‌ها Concrete02 (برای بتن محصور و غیرمحصور با مقاومت کششی خطی) و Steel02 (مدل جیامپ-پینتو برای فولاد تقویت‌کننده) هستند.

# Steel02: tag, Fy, E0, b, R0, cR1, cR2
uniaxialMaterial Steel02 1 400.0e6 200.0e9 0.01 18.5 0.925 0.15

۳. تعریف مقاطع و المان‌های غیرخطی (Sections & Elements)

برای شبیه‌سازی دقیق رفتار غیرخطی تیرها و ستون‌ها، از روش توزیع پلاستیسیته با المان‌های فایبر (الیافی) استفاده می‌شود. المان forceBeamColumn یا dispBeamColumn گزینه‌های استاندارد برای این کار هستند.

۴. اعمال بارهای ثقلی و تحلیل استاتیکی اولیه

پیش از اعمال زلزله، سازه باید تحت بارهای ثقلی (مرده و زنده) قرار گیرد. این تحلیل استاتیکی باید به صورت کنترل نیرو انجام شده و در انتهای آن، حالت بارگذاری قفل (Load Const) شود تا اثرات p-delta و نیروهای محوری اولیه در ستون‌ها به درستی در گام دینامیکی لحاظ شوند.

۵. تعریف بارهای دینامیکی و حلگر (Transient Analysis)

در این گام، شتاب‌نگاشت زلزله با دستور pattern UniformExcitation به پای سازه اعمال می‌شود. همچنین برای پایداری حل دینامیکی، استفاده از میرایی رایلی (Rayleigh Damping) متناسب با فرکانس‌های طبیعی مود اول و دوم سازه ضروری است.

۳. مقایسه روش‌های مدل‌سازی: فایبر در مقابل پلاستیسیته متمرکز

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

۴. استانداردهای مدل‌سازی و اعتبارسنجی (Validation) در پایان‌نامه‌ها

برای پذیرش نتایج یک پایان‌نامه یا مقاله علمی مبتنی بر OpenSees، اعتبارسنجی عددی مدل در برابر نتایج آزمایشگاهی معتبر (مانند پایگاه داده دانشگاه برکلی یا PEER) یا کارهای تحقیقاتی مرجع الزامی است.

برای این کار، معمولاً روند زیر پیشنهاد می‌شود:

• مقایسه پریودهای نوسان: ابتدا مقادیر پریود و فرکانس‌های مودهای اولیه سازه را از تحلیل مقدار ویژه (Eigenvalue Analysis) استخراج کرده و با روابط تجربی آیین‌نامه یا نرم‌افزار تجاری همسان (مانند ETABS) مقایسه کنید.
• شبیه‌سازی یک نمونه آزمایشگاهی: یک ستون یا قاب بتنی/فولادی تک‌دهانه که داده‌های آزمایشگاهی چرخه بارگذاری (Cyclic Loading) آن موجود است را شبیه‌سازی کرده و نمودار هیسترزیس خروجی را بر روی داده‌های واقعی منطبق کنید.
• تعیین میزان خطا: اختلاف بین سختی اولیه، مقاومت نهایی و میزان استهلاک انرژی در مدل عددی و آزمایشگاهی نباید بیش از ۵ تا ۱۰ درصد باشد تا مدل شما از نظر علمی معتبر شناخته شود.

۵. اشتباهات رایج و راه‌حل سریع

در جدول و توضیحات زیر، مهم‌ترین دلایل توقف تحلیل و بروز خطا در کدهای OpenSees به همراه راه‌حل‌های سریع مهندسی آورده شده است:

• خطای عدم همگرایی (Convergence Failure):
  
  راه‌حل: الگوریتم حل را پله‌پله تغییر دهید (مثلاً از Newton به NewtonLineSearch یا KrylovNewton). گام زمانی تحلیل دینامیکی (dt) را به نصف یا یک‌دهم کاهش دهید. همچنین مقدار تلرانس خطای همگرایی را موقتاً با دستور test EnergyIncr تغییر دهید.
• خطای ماتریس منفرد (Singular Matrix / Zero Diagonals):
  
  راه‌حل: این خطا نشان‌دهنده ناپایداری سازه‌ای یا رهاشدگی بیش از حد درجات آزادی است. شرایط مرزی تکیه‌گاه‌ها (دستور fix) و اتصالات صلب تیر به ستون را دوباره بررسی کنید.
• در نظر نگرفتن جرم در گره‌ها:
  
  راه‌حل: برای تحلیل‌های دینامیکی، حتماً باید جرم لرزه‌ای سازه را با دستور mass به گره‌های انتقالی اختصاص دهید. در غیر این صورت، شتاب زلزله هیچ نیروی اینرسی به سازه وارد نخواهد کرد.

۶. پرسش‌های متداول دانشجویی