💻 computer science

Cohesive phase-field fracture with an explicit strength surface: an eigenstrain-based return-mapping formulation

यह शोध पत्र कोहेसिव फेज़-फील्ड फ्रैक्चर (cohesive phase-field fracture) के लिए एक आइगेनस्ट्रेन-आधारित रिटर्न-मैपिंग फॉर्मूलेशन प्रस्तुत करता है जो फ्रैक्चर आइगेनस्ट्रेन को स्थानीय कंस्टिट्यूटिव वेरिएबल्स के रूप में मानता है, जिससे अतिरिक्त वैश्विक डिग्री ऑफ फ्रीडम के बिना स्पष्ट स्ट्रेंथ सर्फेस और मानक परिमित-तत्व (finite-element) प्रक्रियाओं का एकीकरण सक्षम होता है और साथ ही मेश-स्वतंत्र परिणाम और सुदृढ़ अभिसरण सुनिश्चित होता है।

मूल लेखक: Tim Hageman

प्रकाशित 2026-05-27

📖 6 मिनट में पढ़ें🧠 गहराई से पढ़ें

CC BY 4.0

मूल लेखक: Tim Hageman

मूल पेपर CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/) के तहत लाइसेंस किया गया है। ✨ नीचे दिए गए पेपर की यह व्याख्या AI से तैयार की गई है। इसे लेखकों ने न तो लिखा है, न इसका समर्थन किया है। तकनीकी सटीकता के लिए मूल पेपर देखें। पूरा डिस्क्लेमर पढ़ें

कल्पना कीजिए कि आप यह अनुमान लगाने की कोशिश कर रहे हैं कि कांच या कंक्रीट का एक टुकड़ा कैसे टूटेगा। लंबे समय तक, कंप्यूटर वैज्ञानिकों के पास उनके सिमुलेशन टूल्स के साथ दो मुख्य समस्याएं थीं:

"कब" की समस्या: वे एक बार शुरू होने के बाद दरार के फैलने का सिमुलेशन आसानी से कर सकते थे, लेकिन वे यह अनुमान नहीं लगा सकते थे कि दरार कब और कहाँ पहली बार दिखाई देगी। यह एक जंगल की आग के रास्ते के मानचित्र को रखने जैसा था, लेकिन यह जानने का कोई तरीका नहीं था कि चिंगारी कहाँ गिरेगी।
"कैसे" की समस्या: उन्हें उन सामग्रियों के सिमुलेशन करने में कठिनाई होती थी जो केवल कांच की तरह टूटती (भंगुर/brittle) नहीं हैं, बल्कि च्युइंग गम या गीली मिट्टी की तरह खिंचती और फटती (संसंजक/cohesive या लचीली/ductive) हैं।

यह पेपर एक नया "स्मार्ट रेसिपी" पेश करता है जो इन दोनों समस्याओं को एक साथ हल करता है। यह कैसे काम करता है, इसके लिए रोजमर्रा के उदाहरणों का उपयोग किया गया है।

पुराना तरीका बनाम नया तरीका

पुराना तरीका (केवल भंगुर/Brittle):
मानक सिमुलेशन को एक सख्त रूलर (पैमाने) की तरह समझें। यदि आप इसे बहुत अधिक मोड़ते हैं, तो यह तुरंत टूट जाता है। कंप्यूटर जानता है कि टूटने के लिए कितनी ऊर्जा चाहिए, लेकिन वह यह नहीं जानता कि इसे तोड़ने के लिए कितने बल की आवश्यकता है। इसे तोड़ने के लिए, इंजीनियरों को अक्सर यादृच्छिक "कमजोर बिंदु" जोड़कर या अनुमान लगाकर धोखाधड़ी करनी पड़ती थी।

नया तरीका ("आइजनस्ट्रेन" का कमाल):
लेखक फ्रैक्चर आइजनस्ट्रेन (fracture eigenstrains) की एक अवधारणा पेश करते हैं। कल्पना करें कि सामग्री केवल एक ठोस ब्लॉक नहीं है, बल्कि एक स्पंज की तरह है जो वास्तव में टूटने से पहले आंतरिक रूप से खुद को "पिचक" (squish) सकती है।

गुप्त पिचक (The Secret Squish): जब आप सामग्री को खींचते हैं, तो यह "गुप्त पिचक" (आइजनस्ट्रेन) बढ़ना शुरू हो जाता है। यह एक छिपे हुए स्प्रिंग की तरह काम करता है जो तनाव (stress) को सोख लेता है।
शक्ति सीमा (The Strength Limit): सामग्री की एक "शक्ति सीमा" होती है (जैसे एक पुल कितना वजन सह सकता है)। जब तक तनाव इस सीमा से नीचे रहता है, गुप्त पिचक छोटा रहता है।
टूटना (The Breaking Point): जैसे ही तनाव इस सीमा तक पहुँचता है, गुप्त पिचक तेजी से बढ़ता है। यह कंप्यूटर को बताता है, "ठीक है, यहाँ एक दरार शुरू हो रही है!" यह "कब" की समस्या को हल करता है।
फटना (The Tear): जैसे-जैसे पिचक बढ़ता है, सामग्री केवल गायब नहीं होती; यह धीरे-धीरे अपनी ताकत खो देती है। यह सिमुलेशन को यह दिखाने की अनुमति देता है कि सामग्री केवल टूटती नहीं है बल्कि खिंचती और फटती (cohesive behavior) है। यह "कैसे" की समस्या को हल करता है।

"लोकल" जादू (अतिरिक्त कंप्यूटरों की आवश्यकता नहीं)

आमतौर पर, सिमुलेशन में इन "गुप्त पिचक" वेरिएबल्स को जोड़ने के लिए कंप्यूटर को एक साथ पूरे ऑब्जेक्ट के लिए एक विशाल, जटिल पहेली को हल करने की आवश्यकता होती है। यह एक सुडोकू पहेली को हलने जैसा है जहाँ हर संख्या एक साथ बोर्ड के अन्य सभी नंबरों को प्रभावित करती है। यह धीमा है और इसके लिए विशेष, महंगे सॉफ्टवेयर की आवश्यकता होती है।

पेपर का नवाचार:
लेखकों ने महसूस किया कि इस "गुप्त पिचक" को अपने पड़ोसियों से बात करने की आवश्यकता नहीं है। इसे केवल अपने विशिष्ट स्थान (सिमुलेशन में एक सिंगल पिक्सेल) पर क्या हो रहा है, इसकी परवाह है।

उपमा: एक स्टेडियम में बैठे लोगों की कल्पना करें। पुराने तरीके में, किसी के भी बैठने से पहले सभी को हाथ पकड़कर एक चाल पर सहमत होना पड़ता था। इस नए तरीके में, प्रत्येक व्यक्ति केवल अपने स्वयं के आराम के स्तर के आधार पर बैठने का निर्णय लेता है।
परिणाम: कंप्यूटर हर एक छोटे स्थान के लिए स्वतंत्र रूप से टूटने के बिंदु की गणना कर सकता है, ठीक वैसे ही जैसे वह प्लास्टिक के मुड़ने की गणना करता है। इसका मतलब है कि इस नए तरीके को विशेष सुपरकंप्यूटर या कोड के पूर्ण पुनर्लेखन की आवश्यकता के बिना मानक इंजीनियरिंग सॉफ्टवेयर में डाला जा सकता है।

उन्होंने क्या परीक्षण किया (द "ड्राइविंग टेस्ट्स")

यह सिद्ध करने के लिए कि यह रेसिपी काम करती है, उन्होंने अपने कंप्यूटर मॉडल पर तीन विशिष्ट "ड्राइविंग टेस्ट" चलाए:

प्लेट में छेद (The Hole in the Plate): उन्होंने एक धातु की प्लेट का सिमुलेशन किया जिसके बीच में एक छेद था, जिसे खींचकर (तनाव/tension) और दबाकर (संपीड़न/compression) अलग किया गया।
- परिणाम: मॉडल ने सही ढंग से भविष्यवाणी की कि बल बहुत अधिक होने पर छेद के किनारों पर दरारें शुरू होंगी। इसने यह भी दिखाया कि संपीड़न (compression) के तहत, सामग्री बिना कुचले हुए फिसल सकती है (जैसे दो किताबें एक-दूसरे के बगल से फिसलती हैं), जिसे पुराने मॉडल अक्सर गलत कर देते थे।
नॉच्ड प्लेट (Shear): उन्होंने एक प्लेट ली जिसके किनारे पर एक कट था और उसे मरोड़ा (twist) गया।
- परिणाम: उन्होंने दिखाया कि "शक्ति" सेटिंग्स को ट्यून करने के आधार पर, दरार या तो सीधे आर-पार (shear) फिसल सकती है या नीचे की ओर फट (tension) सकती है। यह साबित करता है कि मॉडल जटिल, मिश्रित ब्रेकिंग पैटर्न को संभाल सकता है।
डायनेमिक क्रैश (The Dynamic Crash): उन्होंने एक प्लेट का सिमुलेशन किया जिस पर अचानक, तेज़ बल (जैसे विस्फोट या दुर्घटना) लगा।
- परिणाम: जब सामग्री "भंगुर" (कम ऊर्जा टूटने वाली) थी, तो दरार तेजी से फैली और शाखाओं में बंट गई (जैसे बिजली की चमक)। जब उन्होंने सामग्री को "मजबूत" (उच्च ऊर्जा टूटने वाली) बनाया, तो दरार धीमी हो गई और एक एकल रेखा के रूप में बनी रही। इसने साबित किया कि मॉडल स्वाभाविक रूप से "ब्रांचिंग" (शाखाओं में बंटना) के लिए विशिष्ट नियम प्रोग्राम किए बिना जटिल, अराजक टूटने को कैप्चर कर सकता है।

मुख्य निष्कर्ष

पेपर का दावा है कि यह नया तरीका मौजूदा टूल्स के लिए एक "ड्रॉप-इन रिप्लेसमेंट" है। यह इंजीनियरों को निम्न की अनुमति देता है:

दरार कब शुरू होती है (सामग्री की शक्ति के आधार पर)।
यह कैसे टूटती है (एक साफ टूटने से लेकर एक धीमी फटने तक)।
जटिल व्यवहार जैसे कि दरार का शाखाओं में बंटना या फिसलना, वह भी बिना किसी अतिरिक्त कंप्यूटर पावर या विशेष सॉफ्टवेयर के।

लेखक इस बात पर जोर देते हैं कि परिणाम इस बात पर निर्भर नहीं करते कि कंप्यूटर सिमुलेशन कितना "ज़ूम इन" है (mesh-independent), जिसका अर्थ है कि उत्तर भौतिक रूप से वास्तविक हैं, न कि कंप्यूटर ग्रिड का कोई भ्रम। उनके द्वारा उपयोग किया गया सारा कोड किसी के भी उपयोग और निर्माण के लिए खुला है।

तकनीकी सारांश: एक स्पष्ट स्ट्रेंथ सरफेस के साथ सुसंगत फेज़-फील्ड फ्रैक्चर (Cohesive Phase-Field Fracture)

समस्या विवरण
मानक फेज़-फील्ड फ्रैक्चर विधियाँ मौलिक रूप से ग्रिफिथ के भंगुर फ्रैक्चर सिद्धांत (Griffith's brittle fracture theory) पर आधारित हैं। फलस्वरूप, उनमें सामग्री की शक्ति (material strength) के लिए कोई अंतर्निहित निर्देश का अभाव होता है, जिससे दरार के न्यूक्लिएशन (crack nucleation) की निर्भरता एड-हॉक मानदंडों या फेज़-फील्ड लंबाई पैमाने (regularization parameter) पर हो जाती है, न कि किसी भौतिक यील्ड सरफेस पर। इसके अतिरिक्त, ये विधियाँ कोहेसिव फ्रैक्चर व्यवहार और डक्टाइल प्रतिक्रियाओं को कैप्चर करने में संघर्ष करती हैं, क्योंकि डिग्रेडेशन फंक्शन आमतौर पर यूनिएक्सियल टेंशन के लिए कैलिब्रेट किए जाते हैं और स्वाभाविक रूप से मल्टी-एक्सियल लोडिंग या प्रेशर-डिपेंडेंट स्ट्रेंथनिंग तक विस्तारित नहीं होते हैं। हालांकि हालिया आइजनस्ट्रेन-आधारित फॉर्मुलेशन (जैसे विन्सेन्टिनी एट अल., बोरडिन एट अल.) ने फ्रैक्चर ऊर्जा से सामग्री की शक्ति को अलग करने के लिए फ्रैक्चर आइजनस्ट्रेन पेश करके सफलतापूर्वक अलग किया है, लेकिन उनका कार्यान्वयन ग्लोबल एनर्जी-मिनिमाइजेशन फ्रेमवर्क पर निर्भर करता है। यह दृष्टिकोण आइजनस्ट्रेन को अतिरिक्त वैश्विक डिग्री ऑफ फ्रीडम (global degrees of freedom) के रूप में मानता है, जिसके लिए विशेष सॉफ़्टवेयर की आवश्यकता होती है, जो मानक परिमित तत्व (finite-element - FE) कोड के एकीकरण में बाधा डालता है जो पॉइंटवाइज इंटीग्रेशन पॉइंट्स पर कंस्ट्रिट्यूटिव मॉडल का मूल्यांकन करते हैं।

कार्यप्रणाली (Methodology)
यह कार्य आइजनस्ट्रेन-आधारित कोहेसिव फेज़-फील्ड मॉडल को एक स्थानीय कंस्ट्रिट्यूटिव मॉडल के रूप में पुनर्गठित करता है, जो प्लास्टिक रिटर्न-मैपिंग स्कीम्स के समान है, ताकि इसे अतिरिक्त वैश्विक डिग्री ऑफ फ्रीडम के बिना मानक FE कोड में एकीकृत किया जा सके।

शासी समीकरण (Governing Equations): फ्रेमवर्क एक कुल ऊर्जा फंक्शनल का उपयोग करता है जिसमें इलास्टिक ऊर्जा (विस्थापन और फ्रैक्चर आइजनस्ट्रेन पर आधारित), फ्रैक्चर ऊर्जा (फेज़-फील्ड पैरामीटर $\phi$ के माध्यम से वितरित), और एक स्ट्रेंथ पोटेंशियल शामिल है। इलास्टिक स्ट्रेन को $\varepsilon^{el} = \varepsilon - \eta$ के रूप में परिभाषित किया गया है, जहाँ $\eta$ फ्रैक्चर आइजनस्ट्रेन का प्रतिनिधित्व करता है।
स्थानीय कंस्ट्रिट्यूटिव फॉर्मुलेशन: मुख्य नवाचार आइजनस्ट्रेन के विकास को पूरी तरह से इंटीग्रेशन-पॉइंट स्तर पर हल करना है। आइजनस्ट्रेन को $\eta = \lambda G(\varepsilon)$ के रूप में परिभाषित किया गया है, जहाँ $\lambda$ एक मल्टीप्लायर है और $G$ स्ट्रेन अवस्था पर निर्भर दिशा फलन (direction function) है। यह अनुमति देता है कि कंस्ट्रिट्यूटिव अपडेट को स्थानीय रूप से न्यूटन-र्राफसन स्कीम का उपयोग करके हल किया जाए, जैसा कि प्लास्टिसिटी में होता है, न कि वैश्विक स्तर पर।
स्ट्रेंथ मानदंड (Strength Criteria): दो अलग-अलग स्ट्रेंथ पोटेंशियल लागू किए गए हैं:
1. नॉन-स्मूथ $r_1$ मानदंड: इसमें स्वतंत्र टेंसाइल ( $f_t$ ) और शियर ( $f_s$ ) स्ट्रेंथ के साथ उच्च संपीड़न शक्ति (compressive strength) शामिल है ताकि नॉन-पेनेट्रेशन को लागू किया जा सके।
2. स्मूथ ड्रकर-प्रापर-लाइक (Drucker-Prager-like) मानदंड: यह संपीड़न के तहत प्रेशर-डिपेंडेंट स्ट्रेंथनिंग को कैप्चर करता है। यह एक संदर्भ स्ट्रेन पैरामीटर ( $\varepsilon_{ref}$ ) का उपयोग करता है जो टेंसाइल और कंप्रेसिव रेजीम के बीच संक्रमण को निर्देशित करता है, जिससे एक स्मूथ फेलियर सरफेस सुनिश्चित होता है।
एल्गोरिदम कार्यान्वयन: लेखक दोनों स्मूथ और नॉन-स्मूथ मानदंडों के लिए सुसंगत टेंगेंट ऑपरेटर्स (consistent tangent operators) व्युत्पन्न करते हैं ताकि ग्लोबल न्यूटन-र्राफसन सॉल्वर की मजबूती सुनिश्चित हो सके। एक स्टेगर्ड समाधान योजना का उपयोग किया जाता है, जिसमें पहले विस्थापन और आइजनस्ट्रेन को हल किया जाता है, उसके बाद फेज़-फील्ड वेरिएबल को, और अभिसरण (convergence) तक पुनरावृत्ति की जाती है। कार्यान्वयन में क्वाड्रेटिक लैग्रेंज तत्वों का उपयोग किया गया है और इसे FEniCSx में कोड किया गया है।

प्रमुख योगदान

लोकल रिटर्न-मैपिंग फॉर्मुलेशन: प्राथमिक योगदान आइजनस्ट्रेन के विकास को एक स्थानीय कंस्ट्रिट्यूटिव मॉडल में पुनर्गठित करना है। यह अतिरिक्त वैश्विक डिग्री ऑफ फ्रीडम की आवश्यकता को समाप्त करता है, जिससे कोहेसिव फेज़-फील्ड मॉडल को मौजूदा FE कोड में एक "ड्रॉप-इन" रिप्लेसमेंट के रूप में लागू करना संभव हो जाता है।
स्पष्ट स्ट्रेंथ सरफेस (Explicit Strength Surface): विधि दरार न्यूक्लिएशन के लिए एक स्पष्ट सामग्री स्ट्रेंथ सरफेस निर्धारित करती है, जो फ्रैक्चर ऊर्जा रिलीज रेट ( $G_c$ ) से अलग है। यह न्यूक्लिएशन स्ट्रेस और प्रोपेगेशन एनर्जी के स्वतंत्र नियंत्रण की अनुमति देता है।
सुसंगत टेंगेंट ऑपरेटर्स: स्मूथ और नॉन-स्मूथ स्ट्रेंथ सर्फेस के लिए क्लोज्ड-फॉर्म सुसंगत टेंगेंट ऑपरेटर्स का व्युत्पन्न ग्लोबल सॉल्वर में संख्यात्मक स्थिरता और क्वाड्रेटिक कन्वर्जेंस रेट सुनिश्चित करता है।
प्रेशर-डिपेंडेंट मानदंड: एक स्मूथ ड्रकर-प्रापर-लाइक मानदंड का परिचय संपीड़न के तहत प्रेशर-डिपेंडेंट स्ट्रेंथनिंग को मॉडल करने में सक्षम बनाता है, जो एक ऐसी विशेषता है जो अक्सर मानक फेज़-फील्ड फॉर्मुलेशन में अनुपस्थित होती है।

परिणाम और सत्यापन
फ्रेमवर्क को तीन बेंचमार्क समस्याओं के विरुद्ध सत्यापित किया गया है:

छेद के साथ प्लेट (टेंशन/कंप्रेशन): मॉडल मेष-स्वतंत्र व्यवहार प्रदर्शित करता है यदि तत्व का आकार फेज़-फील्ड लंबाई पैमाने ( $\ell/dx \geq 2$ ) के सापेक्ष पर्याप्त छोटा है। यह ब्रिटल से कोहेसिव व्यवहार के संक्रमण को सही ढंग से कैप्चर करता है क्योंकि $G_c$ बढ़ता है। संपीड़न के तहत, मॉडल नॉन-पेनेट्रेशन को लागू करता है और शियर बैंड गठन की सही भविष्यवाणी करता है, हालांकि कम $G_c$ पर हॉरिजॉन्टल टेंसाइल क्रैकिंग और डायगोनल शियर क्रैकिंग के बीच प्रतिस्पर्धा देखी जाती है।
सिंगल-एज नॉच्ड प्लेट (शियर): मॉडल सफलतापूर्वक मिक्स-मोड फ्रैक्चर को कैप्चर करता है। कम संदर्भ स्ट्रेन (उच्च संपीड़न शक्ति) के साथ ड्रकर-प्रापर-लाइक मानदंड का उपयोग करके, दरार का पथ शियर-ड्रिवन से टेंसाइल-ड्रिवन में बदल जाता है, जो जटिल स्ट्रेस स्टेट्स को बिना किसी एड-हॉक स्प्लिटिंग स्कीम के संभालने की मॉडल की क्षमता को प्रदर्शित करता है।
डायनेमिक फ्रैक्चर: डायनेमिक लोडिंग परिदृश्यों में, मॉडल दरार शाखाओं (crack branching) की घटनाओं को स्वाभाविक रूप से कैप्चर करता है। परिणाम दिखाते हैं कि $G_c$ को बढ़ाने से क्रैक ब्रांचिंग दब जाती है और प्रसार गति कम हो जाती है, जबकि स्ट्रेंथ सरफेस प्रारंभिक न्यूक्लिएशन को नियंत्रित करती है। मॉडल पीक लोड और दरार पथों के संबंध में मेष-स्वतंत्र और लेंथ-स्केल-स्वतंत्र परिणाम उत्पन्न करता है।

महत्व और दावे
पेपर का दावा है कि यह फ्रेमवर्क कोहेसिव फ्रैक्चर प्रक्रियाओं को सिम्युलेट करने के लिए एक शक्तिशाली और लचीला उपकरण प्रदान करता है। स्ट्रेंथ सरफेस को फ्रैक्चर ऊर्जा और फेज़-फील्ड लंबाई पैमाने से अलग करके, यह विधि निम्नलिखित की अनुमति देती है:

एक निर्धारित सामग्री स्ट्रेंथ के आधार पर दरार न्यूक्लिएशन की भविष्यवाणी करना।
फ्रैक्चर ऊर्जा रिलीज रेट ( $G_c$ ) के माध्यम से दरार प्रसार और ब्रिटल-टू-कोहेसिव संक्रमण का नियंत्रण।
एड-हॉक मानदंडों के बिना क्रैक ब्रांचिंग और मिक्स-मोड प्रोपेगेशन जैसी जटिल घटनाओं का सिमुलेशन।
वैश्विक सॉल्वर या तत्व फॉर्मुलेशन को संशोधित किए बिना मानक परिमित तत्व कोड में निर्बाव एकीकरण।

लेखक इस बात पर जोर देते हैं कि लंबाई पैमाना $\ell$ केवल फेज़-फील्ड वितरण चौड़ाई के लिए एक संख्यात्मक रेगुलराइजेशन पैरामीटर के रूप में कार्य करता है और यदि मेष पर्याप्त रूप से रिफाइंड है, तो यह भौतिक सामग्री की शक्ति या ऊर्जा अपव्यय को प्रभावित नहीं करता है। आगे के अनुप्रयोग और सत्यापन को सुगम बनाने के लिए सोर्स कोड को ओपनली उपलब्ध कराया गया है।

पुराना तरीका बनाम नया तरीका

"लोकल" जादू (अतिरिक्त कंप्यूटरों की आवश्यकता नहीं)

उन्होंने क्या परीक्षण किया (द "ड्राइविंग टेस्ट्स")

मुख्य निष्कर्ष

इसी तरह के और पेपर