۱۰ خطای رایج در مدلسازی گودبرداری عمیق با PLAXIS 2D

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

در این مقاله، ۱۰ اشتباه رایج و مهم در مدل‌سازی گودبرداری‌های عمیق با PLAXIS 2D را از سطح مفاهیم کلی تا جزئیات تخصصی بررسی می‌کنیم.

۱. انتخاب مدل رفتاری نامناسب

یکی از رایج‌ترین خطاهای مهندسان، استفاده از مدل رفتاری موهر-کولمب (Mohr-Coulomb) برای تمامی خاک‌ها بدون توجه به شرایط پروژه و اهداف تحلیل است. مدل موهر-کولمب رفتار غیرخطی خاک را به دو بخش ساده (الاستیک–پلاستیک کامل) تقلیل می‌دهد و سختی خاک (E₅₀) را تا لحظه گسیختگی کاملاً ثابت فرض می‌کند. مشکل بزرگ‌تر این است که این مدل تفاوتی بین سختی بارگذاری اولیه و سختی باربرداری/بارگذاری مجدد (Eur) قائل نمی‌شود. در حالی که در واقعیت، خاک در حالت باربرداری (مانند کف گودبرداری) معمولاً رفتار سخت‌تری از خود نشان می‌دهد و مدول Eur اغلب چندین برابر E₅₀ است.

استفاده از مدل موهر-کولمب می‌تواند باعث شود میزان بالاآمدگی (Heave) کف گودال و همچنین نشست‌های پشت دیوار بیش از مقدار واقعی پیش‌بینی شود. برای دستیابی به نتایج واقع‌بینانه‌تر، معمولاً از مدل‌های پیشرفته‌تری مانند Hardening Soil (HS) استفاده می‌شود که وابستگی سختی به سطح تنش را در نظر گرفته و پارامترهای سه‌گانه E₅₀، E_ur  و E_oed را در مدل‌سازی لحاظ می‌کنند.

۲. تعیین نامناسب شرایط مرزی (Boundary Conditions)

یکی از خطاهای مهم در مدل‌سازی اجزای محدود، انتخاب نامناسب فاصله و نوع شرایط مرزی (Boundary Conditions) است. در PLAXIS، مرزهای مدل باید به اندازه‌ای از ناحیه مورد مطالعه فاصله داشته باشند که بر توزیع تنش‌ها و تغییرشکل‌های ناشی از گودبرداری تأثیر نگذارند.

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

۳. مش‌بندی (Meshing) نامناسب

کیفیت شبکه اجزای محدود (Mesh) مستقیماً بر دقت نتایج و زمان محاسبات تأثیر می‌گذارد. استفاده از یک مش یکنواخت برای کل مدل اشتباه است. پلکسیس امکان استفاده از مش‌های ریز (Refined Mesh) را فراهم کرده است که باید حتماً در نواحی حساس و پرتنش مانند اطراف دیوارهای حائل، فونداسیون‌ها، تونل‌ها و محل اتصال انکرها استفاده شود. در مقابل، برای مرزهای دورترِ مدل که تغییرات تنش در آن‌ها حداقل است، باید از مش‌های درشت‌تر (General Mesh) استفاده کرد تا از افزایش بی‌دلیل و شدید زمان محاسبات جلوگیری شود.

۴. تعیین اشتباه فشار آب حفره‌ای و شرایط زهکشی

در مدل‌سازی خاک‌های رسی طی گودبرداری‌های سریع در نرم‌افزار PLAXIS، انتخاب صحیح نوع زهکشی (Drainage Type) اهمیت ویژه‌ای دارد. در این شرایط، بسیاری از متخصصان استفاده از حالت Undrained A همراه با مدل‌های رفتاری پیشرفته نظیر Hardening Soil را توصیه می‌کنند؛ زیرا در این رویکرد، مقاومت و سختی خاک بر مبنای پارامترهای تنش مؤثر تعریف شده و اضافه فشار آب حفره‌ای به‌صورت طبیعی در فرآیند تحلیل تولید می‌شود.

با این حال، استفاده از حالت Undrained A  در کنار مدل ساده Mohr-Coulomb  همواره نتایج رضایت‌بخشی به همراه ندارد. به دلیل محدودیت‌های ذاتی این مدل در شبیه‌سازی رفتار غیرخطی و وابستگی سختی به تنش، ممکن است مقاومت زهکشی‌نشده بسیج‌شده و تغییرشکل‌های خاک به‌طور واقع‌بینانه بازتولید نشوند. به همین دلیل معمولاً توصیه می‌شود حالت Undrained A همراه با مدل‌های پیشرفته‌تری نظیر Hardening Soil یا Hardening Soil Small به‌کار گرفته شود.

از سوی دیگر، بسیاری از کاربران از ترکیب Undrained B  و Mohr-Coulomb  استفاده می‌کنند که در آن مقاومت برشی زهکشی‌نشده (Su یا Cu) مستقیماً به مدل وارد می‌شود. اگرچه این روش برای برخی تحلیل‌های پایداری کوتاه‌مدت کاربرد دارد، مطالعات مختلف نشان داده‌اند که در مسائل تغییرشکل‌محور نظیر گودبرداری‌های عمیق، این رویکرد می‌تواند منجر به کم‌برآوردی نیروهای سازه‌ای و تغییرشکل‌ها نسبت به رفتار واقعی شود. از این رو، برای پروژه‌هایی که کنترل نشست زمین، جابه‌جایی دیوار و عملکرد سازه نگهبان اهمیت بالایی دارد، استفاده از مدل‌های مبتنی بر تنش مؤثر معمولاً ترجیح داده می‌شود.

نکته مهم دیگر آن است که تنظیمات Drainage Type تنها در تحلیل‌های Plastic و Safety مورد استفاده قرار می‌گیرند. در مقابل، در تحلیل‌های وابسته به زمان نظیر Consolidation  و Fully-coupled Flow-Deformation ، تولید و استهلاک فشار آب حفره‌ای مستقیماً بر اساس نفوذپذیری خاک و معادلات جریان آب زیرزمینی محاسبه می‌شود و انتخاب حالت‌های Undrained A  یا Undrained B دیگر نقشی در فرآیند تحلیل نخواهد داشت.

۵. تولید نادرست تنش‌های اولیه (Initial Stresses) در مدل

تولید تنش‌های اولیه، گام صفر هر تحلیل ژئوتکنیکی است. بسیاری از مهندسان عادت دارند این کار را با استفاده از روش ضریب فشار خاک در حالت سکون (K₀ Procedure) انجام دهند که بر پایه رابطه σ′h0 = K₀ σ′v0 استوار است.

با این حال، باید توجه داشت که این روش تنها در شرایطی معتبر است که سطح زمین، لایه‌بندی خاک و سطح آب زیرزمینی کاملاً افقی باشند. در صورت وجود شیب طبیعی زمین، لایه‌بندی‌های نامنظم یا هندسه‌های پیچیده، استفاده از روش K₀ می‌تواند منجر به ایجاد نیروهای نامتعادل (Unbalanced Forces) و تنش‌های برشی غیرواقعی در مدل شود.

در چنین شرایطی، استفاده از روش بارگذاری ثقلی (Gravity Loading) توصیه می‌شود. در این روش، تنش‌های اولیه از طریق اعمال تدریجی وزن خاک ایجاد شده و میدان تنش به‌صورت طبیعی در مدل شکل می‌گیرد. در نتیجه، تعادل تنش‌ها با دقت بیشتری برقرار شده و احتمال ایجاد تنش‌های مصنوعی ناشی از فرضیات ساده‌کننده کاهش می‌یابد.

۶. انتخاب رفتار الاستیک به جای الاستوپلاستیک برای المان‌های ورق (Plate)

برای شبیه‌سازی دیوارهای حائل، شمع‌ها و دیوارهای دیافراگمی از المان ورق (Plate) استفاده می‌شود. اگر متریالِ این المان روی حالت الاستیک (Elastic) تنظیم شود، نرم‌افزار فرض می‌کند که مقطع سازه ظرفیت تحمل بی‌نهایت لنگر خمشی و نیروی محوری را دارد و هیچ‌گاه تسلیم نمی‌شود. این فرض ایمنی طرح را نقض می‌کند. همواره توصیه می‌شود رفتار المان ورق به صورت الاستوپلاستیک (Elastoplastic) تعریف شود و پارامترهای مقاومت نهایی مقطع شامل حداکثر لنگر پلاستیک (M_p) و حداکثر نیروی محوری (N_p) به درستی وارد شوند تا ظرفیت سازه محدود گردد.

۷. وارد کردن اشتباهِ وزن المان دیوار

یکی از خطاهای نسبتاً رایج در مدل‌سازی دیوارهای نگهبان، نحوه تعریف وزن المان‌های Plate در PLAXIS است. المان Plate در نرم‌افزار به‌صورت یک خط با ضخامت مجازی روی محیط خاک قرار می‌گیرد و جایگزین هندسه واقعی دیوار نمی‌شود. بنابراین خاک موجود در محل دیوار همچنان در محاسبات مدل حضور دارد.

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

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

W_plate = (γ_concrete − γ_soil) × d_real

که در آن:

• γ_concrete  وزن مخصوص مصالح دیوار (بتن یا فولاد معادل)
• γ_soil  وزن مخصوص خاک جایگزین‌شده
• d_real  ضخامت واقعی دیوار

عدم توجه به این موضوع می‌تواند به ویژه در دیوارهای دیافراگمی، شمع‌های درجا و سایر سازه‌های نگهبان سنگین، بر وضعیت تنش اولیه و نتایج تحلیل تأثیر قابل توجهی بگذارد.

۸. استفاده نادرست از ضرایب سطح تماس (R_inter)

المان‌های سطح تماس (Interface) برای شبیه‌سازی افت اصطکاک و لغزش بین المان‌های سازه‌ای و خاک استفاده می‌شوند. مقدار R_inter = 1.0  به معنای عدم کاهش مقاومت برشی در سطح تماس است، در حالی که برای بسیاری از مسائل عملی، مقدار این ضریب کمتر از یک در نظر گرفته می‌شود. به‌طور معمول، برای تماس ماسه و فولاد مقادیری در حدود 6/0 تا 7/0، برای رس و فولاد حدود 5/0 و برای تماس خاک و بتن بین 8/0 تا 0/1 توصیه می‌شود.

با این حال، یکی از خطاهای رایج، اعمال ضریب کاهش مقاومت در بخش تزریق‌شده انکرهای تحت فشار (Grouted Ground Anchors) است. در این حالت، دوغاب تزریق‌شده تحت فشار معمولاً پیوستگی مناسبی با خاک اطراف ایجاد می‌کند و گسیختگی غالباً در داخل توده خاک رخ می‌دهد، نه در مرز تماس دوغاب و خاک. بنابراین، اعمال کاهش اصطکاک در این ناحیه معمولاً توجیه فنی نداشته و مقدار Rinter = 1.0 انتخاب مناسب‌تری خواهد بود. همچنین باید توجه داشت که در مدل‌سازی متداول انکرها با استفاده از المان Embedded Beam Row در PLAXIS، اساساً نیازی به تعریف Interface وجود ندارد؛ زیرا اندرکنش بین دوغاب و خاک از طریق پارامترهای تعبیه‌شده در خود المان، از جمله Tskin (حداکثر مقاومت اصطکاکی واحد طول)، مدل می‌شود. در نتیجه، ظرفیت انتقال بار طول گیردار انکر مستقیماً توسط پارامترهای Embedded Beam کنترل شده و تعریف Interface مجزا در اطراف آن نه‌تنها ضروری نیست، بلکه می‌تواند موجب برداشت نادرست از مکانیزم واقعی انتقال بار در سیستم شود.

۹. عدم رعایت مراحل دقیق اجرایی (Staged Construction)

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

۱۰. خطاهای رایج در مدل‌سازی طول آزاد و گیردار انکرها

دل‌سازی مهارها (Ground Anchors) یکی از حساس‌ترین بخش‌های تحلیل در PLAXIS است و بخش قابل توجهی از خطاهای مدل‌سازی در این قسمت رخ می‌دهد. انکر باید به صورت فیزیکی به دو المان مجزا تفکیک شود:

• بخش آزاد (Free Length): باید با المان Node-to-Node Anchor مدل شود که رفتار فنر محوری داشته و هیچ اندرکنشی با خاک پیرامون ندارد.
• بخش گیردار یا تزریق‌شده (Bonded Length): باید با المان Embedded Beam Row مدل شود تا انتقال بار اصطکاکی بین دوغاب و خاک از طریق پارامترهایی نظیر T_skin شبیه‌سازی شود.

یکی از خطاهای رایج هنگام تعریف المان Node-to-Node Anchor، تقسیم دستی سختی محوری (EA) و ظرفیت کششی (Fmax)  بر فاصله انکرها (Spacing) است. در PLAXIS، مقادیر واردشده به‌صورت خودکار و بر اساس مقدار تعریف‌شده در فیلد L_spacing  به شرایط کرنش مسطح (Plane Strain) تبدیل می‌شوند. بنابراین اگر کاربر پیش از ورود داده‌ها، مقادیر EA و Fmax را نیز بر فاصله انکرها تقسیم کند، این کاهش دوبار اعمال شده و ظرفیت و سختی سیستم کمتر از مقدار واقعی در نظر گرفته خواهد شد. از این رو، پارامترهای سختی و مقاومت باید همواره بر اساس مشخصات واقعی یک انکر منفرد وارد نرم‌افزار شوند.

نتیجه‌گیری

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

در گودبرداری‌های عمیق، عواملی نظیر انتخاب مدل رفتاری مناسب، تعریف صحیح تنش‌های اولیه، در نظر گرفتن رفتار زهکشی‌نشده خاک، مدل‌سازی واقع‌بینانه سیستم‌های نگهبان و رعایت توالی واقعی اجرا، تأثیر مستقیمی بر پیش‌بینی جابه‌جایی دیوارها، نشست زمین و نیروهای وارد بر المان‌های سازه‌ای دارند. از این رو، استفاده از PLAXIS نباید صرفاً به عنوان یک ابزار نرم‌افزاری تلقی شود، بلکه باید در کنار دانش مکانیک خاک، تجربه مهندسی و قضاوت فنی به کار گرفته شود.

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