מודל פרידינמיקה מבוסס-מצב לא-אורתודוקסי להתפשטות סדקים בסלעים: שיטה משולבת מבוססת-מתח ואנרגיה לשבר

למה חשובה שבירת סלעים

מתחת לכבישים וערים ועד למאגרי אחסון של פחמן דו-חמצני בעומק האדמה, פרויקטים מודרניים רבים תלויים באופן שבו סלעים סדוקים ושברים מתפתחים. לצפייה בזמן אמת בצמיחת סדקים בתוך בלוק סלע מוצק דרוש ציוד יקר ומורכב ולעתים בלתי אפשרי. מאמר זה מציג שיטת סימולציה חדשה במחשב לדימות האופן שבו סדקים מתחילים, גדלים, סוטים ומתחברים בתוך הסלע, במיוחד תחת התערובת המורכבת של לחיצה והחלקה שמתרחשת בתנאי הנדסה אמיתיים.

דרך חדשה להבנת הסדיקה

רוב המודלים הממוחשבים המסורתיים מתייחסים לסלע כאל רשת רציפה של נקודות ש"מרגישות" רק את שכנותיהן המיידיות. זו גישה עובדתית עד שמופיע סדק, שכן סדק הוא מהותו שבר פתאומי ברציפות החומר. השיטה שבמחקר זה, המכונה פרידינמיקה, מתחילה מתמונה שונה: כל קטע זעיר של החומר מתקשר ישירות עם רבים אחרים במרחק מסוים, מחוברים בעזרת קישורים בלתי נראים. כאשר הסלע עומס, הקישורים הללו נמתחים, נדחסים או מחליקים; אם הם נלחצים מעבר לגבול, הם נשברים, וכתוצאה מכך נוצרים סדקים באופן טבעי בלי חוקים חיצוניים או רימשינג מחדש.

תיקון ליקויים נסתריים במודלים קודמים

למרות שכיוּן זה של קישורים יעיל, גרסאות מתקדמות קודמות סבלו מבעיה נומרית עדינה. בשל האופן שבו המודל מממץ את הדפורמציה, עשוּתם לעתים מצב של "מצבי אפס-אנרגיה"—תבניות תנועה רועדות שעלות האנרגיה שלהן כמעט אפס ואינן פיזיקליות. אלו הופיעו כהזזות משומרות ויכלו לשבש את החיזוי של מסלולי הסדקים. המחברים מסירים את חולשה זו על ידי שיוך לכל קשר מדד דפורמציה מקומי שנבנה במדויק מתוך ממוצע של שתי קצותיו ותוקן כך שמיקומי ההתחלה והמידע של אותו קשר יהיו תואמים בהחלט. תיאור ברמת הקשר הזה משיב תמונה מקומית עקבית של מתיחה וגזירה ומקטין באופן ניכר את התנודות הלא-פיזיקליות.

לימוד המודל להבחין בין גזירה למתיחה

סדקים בסלע אינם זהים זה לזה. תחת מתיחה טהורה, סדקים נוטים להיפתח בכיוון ישר, אך תחת דחיסה בשילוב החלקה לצד, הם מתעקמים, מסתעפים ויוצרים דפוסים מורכבים של שברים מתיחה, גזירה ומעורבים. קריטריונים הנדסיים מקובלים כמו חוק מוהר–קולומב לא תמיד היו אמינים כאשר הוטמעו ישירות בפרידינמיקה, במיוחד כי הם מתעלמים במידה רבה מרכיב המתוח הביניים. בעבודה זו, המחברים משלבים קריטריון מדויק יותר של "אנרגיית גזירה משולשת" ברמת הקשר. במקום להשוות רק בין מתחים מקסימליים ומינימליים, גישה זו מודדת את אנרגיית הגזירה על שלוש מישורים פנימיים אפשריים ומשכללת את השפעת המתוח האמצעי. קשר נכשל אם או שהמתח המתיחה שלו חורג מעמידת המתיחה של הסלע או שאנרגיית הגזירה המצטברת שלו עוברת סף הקשור לקוהזיה של הסלע. הדבר מאפשר לסימולציה להבחין בין כשלי פתיחה לכשלי החלקה בצורה שתואמת יותר לתצפיות במעבדה.

בדיקת השיטה

כדי להראות שהמודל המשופר אינו רק מתמטי אלא מועיל בפועל, המחברים בוחנים אותו כנגד מספר בעיות מדד. תחילה מושכים לוח עם חור עגול ומשווים את ההסטות לתוצאות אלמנט סופי סטנדרטיות, ומוצאים שדות חלקים כמעט זהים ללא הארטיפקטים של מצבי האפס-אנרגיה הקודמים. לאחר מכן מדמים קורה עם שני ליקויים בקצה תחת ניסוח גזירה בארבע נקודות, ניסוי קלאסי בפגיעת סדיקה מעורבת. מסלולי הסדק המתעקמים, העומס השיא ועקומת העומס–ההסטה הכללית חזקים בהתאמה לניסויים ולמחקרים מספריים איכותיים אחרים. בהמשך, דגימות כפיפה חצי-מעגליות עם חריצים מוטים נבדקות כך שמצב השבר נשנה בהדרגה מפתיחה טהורה לשבר מעורב-פתיחה/גזירה. המודל משחזר את זוויות ייזום הסדק ואת המסלולים הסופיים במגוון זוויות חריצים. לבסוף, הוא מטפל בבלוקים סלעיים ריאליסטיים יותר עם פגמים קיימים תחת דחיסה וגזירה. הסימולציות תופסות דפוסים ידועים כגון סדקי כנף היורים קדימה מטיפי הפגמים, סרטי גזירה מזוגגים, וקישורים מורכבים על פני גשרים סלעיים, כל זאת בהתאמה לתמונות מעבדתיות קודמות.

מה זה אומר להנדסת סלעים

בסך הכל, המחקר מדגים כי תיאור מדוקדק של דפורמציה וכשל ברמת הקשר הבודד משנה באופן מהותי את היכולת לחזות כיצד סלעים שבורים באמת. על ידי הסרת מצבי אנרגיה מדומים ושימוש בכלל כשל מבוסס אנרגיית גזירה המכבד את מצב המתחים התלת-ממדי המלא, המודל יכול לעקוב היכן וכיצד סדקים מתחילים, כיצד הם מתעקמים תחת מתחים משתנים וכיצד סדקים נפרדים מתחברים בסופו של דבר. אם כי עדיין תובעני חישובית ומבוסס על תכונות סלע אידיאליות, מסגרת הפרידינמיקה המשופרת הזו מציעה "מעבדה" נומרית אמינה יותר לחקירת התפתחות סדקים בחומרים שבירים, עם פוטנציאל ברור לשיפור בבטיחות תכנון בכרייה, בבנייה תת-קרקעית ובמערכות אנרגיה גיאוטכניות.

ציטוט: Gong, B., Song, Y., Zhang, L. et al. Non-ordinary state-based peridynamics model for rock crack propagation: a combined stress-energy fracture method. Sci Rep 16, 11386 (2026). https://doi.org/10.1038/s41598-026-40833-8

מילות מפתח: דגום שברים בסלעים, פרידינמיקה, התפשטות סדקי גזירה, סימולציה נומרית, חומרים שבירים