A robust mixed finite element formulation for third medium contact

तकनीकी सारांश: तृतीय माध्यम संपर्क (Third Medium Contact) के लिए एक सुदृढ़ मिश्रित परिमित तत्व (Finite Element) सूत्रीकरण

समस्या विवरण
शास्त्रीय परिमित तत्व संपर्क सूत्रीकरण (finite element contact formulations) संपर्क इंटरफेस के स्पष्ट विविक्तीकरण (explicit discretization) पर निर्भर करते हैं, जो दंड (penalty), लाग्रेंज मल्टीप्लायर, या बैरियर विधियों के माध्यम से बाधाओं को लागू करते हैं। हालांकि ये परिपक्व हैं, लेकिन ये बड़े स्लाइडिंग, स्व-संपर्क (self-contact), गंभीर मेश विरूपण (mesh distortion), या टोपोलॉजी अनुकूलन (topology optimization) जैसी स्थितियों में, जहाँ संपर्क सीमाएँ पूर्व-निर्धारित नहीं होती हैं, महत्वपूर्ण एल्गोरिदम संबंधी बोझ का सामना करते हैं। "तृतीय माध्यम संपर्क" (TMC) विधि इस समस्या को हल करने के लिए एक अत्यधिक लचीले काल्पनिक माध्यम (fictitious medium) का उपयोग करती है जो संभावित रूप से संपर्क करने वाली निकायों के बीच के स्थान को भर देता है। यह माध्यम संपर्क से पहले नरम रहता है लेकिन संपीड़न (compression) के तहत अत्यधिक कठोर हो जाता है, जिससे संपर्क बल निरंतरता (continuum) के माध्यम से प्रसारित होते हैं।

हालाँकि, TMC में एक केंद्रीय चुनौती स्थिरीकरण (stabilization) है। माध्यम संपर्क से पहले इतना नरम होना चाहिए कि स्पूरियस बलों (spurious forces) से बचा जा सके, और इतना स्थिर भी होना चाहिए कि मेश विरूपण के बिना अत्यधिक संपीड़न और कतरनी (shear) को सहन कर सके। मौजूदा स्थिरीकरण रणनीतियोंों को अक्सर उच्च-क्रम के परिमित तत्वों (द्वितीय विस्थापन डेरिवेटिव की गणना करने के लिए) की आवश्यकता होती है, या वे वैश्विक सहायक क्षेत्रों (global auxiliary fields) को पेश करती हैं जो सिस्टम के आकार को बढ़ा देते हैं, या विशिष्ट विविक्तीकरणों पर निर्भर करती हैं। एक सुदृढ़ स्थिरीकरण तंत्र की आवश्यकता है जो सरल, गणनात्मक रूप से कुशल प्रथम-क्रम के परिमित तत्वों के अनुकूल हो और विस्थापन के दूसरे डेरिवेटिव के प्रत्यक्ष मूल्यांकन से बच सके।

कार्यप्रणाली (Methodology)
लेखक एक मिश्रित परिमित तत्व सूत्रीकरण प्रस्तावित करते हैं जो तृतीय माध्यम के लिए एक सहायक-क्षेत्र स्थिरीकरण (auxiliary-field stabilization) पेश करता है। विधि का मुख्य भाग निम्नलिखित है:

1. सहायक क्षेत्र का परिचय: एक विरूपण-ग्रेडिएंट-समान क्षेत्र, जिसे $\Theta$ द्वारा दर्शाया गया है, तृतीय माध्यम में पेश किया जाता है और इसे भौतिक विस्थापन क्षेत्र $u$ से स्वतंत्र एक अज्ञात चर के रूप में माना जाता है।
2. दुर्बल युग्मन (Weak Coupling): सहायक क्षेत्र $\Theta$ को एक दंड पद ($W^p_{tm} \propto |\Theta - F|^2$) के माध्यम से भौतिक विरूपण ग्रेडिएंट $F$ के साथ दुर्बल रूप से जोड़ा जाता है। यह सुनिश्चित करता है कि $\Theta$ भौतिक विरूपण का अनुसरण करे।
3. ग्रेडिएंट नियमितीकरण (Gradient Regularization): भौतिक विरूपण ग्रेडिएंट ($\nabla F$) के ग्रेडिएंट को दंडित करने के बजाय, जो विस्थापन के दूसरे डेरिवेटिव की आवश्यकता रखता है, यह सूत्रीकरण सहायक क्षेत्र ($\nabla \Theta$) के ग्रेडिएंट को दंडित करता है। नियमितीकरण ऊर्जा को $W^r_{tm} \propto |\nabla \Theta|^2$ के रूप में परिभाषित किया गया है।
4. मिश्रित विविक्तीकरण (Mixed Discretization): यह सूत्रीकरण $u$ और $\Theta$ के स्वतंत्र प्रक्षेप (interpolation) की अनुमति देता है। लेखक $\Theta$ के लिए सतत (continuous) और एलिमेंट-वार असंतत (element-wise discontinuous) सन्निकटन दोनों की जांच करते हैं।
   • सतत $\Theta$: नोडल मान तत्वों के बीच साझा किए जाते हैं, जिससे अंतर-एलिमेंट युग्मन (inter-element coupling) बनता है।
   • असंतत $\Theta$: स्वतंत्रता की कोटियाँ (degrees of freedom) प्रत्येक तत्व के लिए स्थानीय होती हैं। इन्हें स्टैटिक कंडेंसेशन (static condensation) के माध्यम से समाप्त किया जा सकता है, जिससे वैश्विक सिस्टम केवल विस्थापन स्वतंत्रता की कोटियों को बनाए रखता है।

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

• निम्न-क्रम अनुकूलता (Low-Order Compatibility): यह सूत्रीकरण विस्थापन के दूसरे डेरिवेटिव की गणना किए बिना ग्रेडिएंट-प्रकार का नियमितीकरण प्राप्त करता है, जिससे यह मानक प्रथम-क्रम (रैखिक) और द्वितीय-क्रम के परिमित तत्वों के अनुकूल बन जाता है।
• स्थिरीकरण का तंत्र: पेपर स्पष्ट करता है कि निम्न-क्रम के विस्थापन तत्वों के लिए (जहाँ $F$ एलिमेंट-वाइज कांस्टेंट है), $\Theta$ का सतत प्रक्षेप आवश्यक है। यह $F$ के एलिमेंट-वाइज जंप्स को $\Theta$ के गैर-शून्य ग्रेडिएंट में बदल देता है, जिससे नियमितीकरण सक्रिय होता है। इसके विपरीत, असंतत $\Theta$ क्षेत्र अंतर-एलिमेंट ग्रेडिएंट उत्पन्न नहीं कर सकते, जिससे निम्न-क्रम तत्वों के लिए स्थिरीकरण प्रभाव कमजोर हो जाता है।
• स्टैटिक कंडेंसेशन: असंतत सहायक क्षेत्रों के लिए, विधि यह प्रदर्शित करती है कि अतिरिक्त अज्ञात चरों को स्थानीय रूप से कंडेंस किया जा सकता है, जिससे वैश्विक सिस्टम के आकार में वृद्धि को रोका जा सकता है जबकि स्थानीय स्थिरीकरण बना रहता है।

संख्यात्मक परिणाम
इस सूत्रीकरण का मूल्यांकन कई बेंचमार्क समस्याओं के माध्यम से किया गया:

• दो-ब्लॉक संपीड़न (Two-Block Compression): प्रदर्शित किया गया कि तृतीय-माध्यम स्केलिंग पैरामीटर $\gamma$ अवशिष्ट अंतराल (residual gap/संपर्क सटीकता) और गैररेखीय अभिसरण मजबूती के बीच संतुलन को नियंत्रित करता है। छोटा $\gamma$ अधिक सटीक संपर्क प्रदान करता है लेकिन अधिक न्यूटन पुनरावृत्तियों (Newton iterations) की आवश्यकता होती है। विधि संयुक्त संपीड़न और कतरनी के तहत स्थिर रही।
• बॉक्स के भीतर स्व-संपर्क (Self-Contact within a Box):
  • प्रक्षेप तुलना: सतत सहायक क्षेत्रों ($T_u^1 T_\Theta^1$, $T_u^2 T_\Theta^2$) ने मजबूत संपर्क क्लोजर प्रदान किया। मिश्रित युग्म $T_u^2 T_\Theta^1$ (क्वाड्रेटिक विस्थापन, लीनियर सहायक) ने $T_u^2 T_\Theta^2$ के समान परिणाम दिए, जिसमें कम सहायक स्वतंत्रता की कोटियों का उपयोग हुआ।
  • असंतत बनाम सतत: असंतत सहायक क्षेत्रों ने कम प्रभावी नियमितीकरण दिखाया, विशेष रूप से रैखिक तत्वों के लिए जहाँ ग्रेडिएंट पद लुप्त हो गया या तत्वों को जोड़ने में विफल रहा।
  • विरूपण-ग्रेडिएंट औसत के साथ तुलना: Faltus et al. [18] की विधि की तुलना में, प्रस्तावित मिश्रित सूत्रीकरण को काफी कम लोड स्टेप्स और न्यूटन पुनरावृत्तियों की आवश्यकता पड़ी, जो बेहतर गैररेखीय मजबूती प्रदर्शित करता है। हालांकि मिश्रित सूत्रीकरण में बड़ा बीजगणितीय सिस्टम था, लेकिन बेहतर अभिसरण ने लागत की भरपाई की, हालांकि अत्यंत सूक्ष्म मेश ($>10^6$ अज्ञात) के लिए सहायक DOFs की लागत प्रभावी हो जाती है।
• C-आकार का बॉक्स: सूत्रीकरण ने बड़े रोटेशन, प्रगतिशील स्व-संपर्क और गंभीर तृतीय-माध्यम संपीड़न को सफलतापूर्वक कैप्चर किया। बल-विस्थापन प्रतिक्रिया शास्त्रीय नोड-टू-सेगमेंट संपर्क सूत्रीकरण के साथ निकटता से मेल खाती है, जिसमें संपर्क में सुचारू संक्रमण होता है।
• 3D विस्तार: तीन-आयामी स्व-संपर्क उदाहरण ने प्लेन स्ट्रेन सेटिंग्स से परे इसकी प्रयोज्यता की पुष्टि की।

महत्व और दावे
पेपर का दावा है कि प्रस्तावित सहायक-क्षेत्र स्थिरीकरण शास्त्रीय संपर्क एल्गोरिदम और मौजूदा TMC नियमितीकरणों के लिए एक सुदृढ़, निम्न-क्रम-अनुकूल विकल्प प्रदान करता है। विस्थापन के दूसरे डेरिवेटिव के प्रत्यक्ष मूल्यांकन से बचकर, यह जटिल संपर्क समस्याओं वाले जटिल विरूपण और स्व-संपर्क के लिए सरल, गणनात्मक रूप से सस्ते परिमित तत्वों (रैखिक और द्विघात) के उपयोग को सक्षम बनाता है।

लेखक इस बात पर जोर देते हैं कि विधि की प्रभावशीलता सहायक क्षेत्र की निरंतरता से आंतरिक रूप से जुड़ी हुई है; निम्न-क्रम विस्थापन तत्वों के लिए ग्रेडिएंट-आधारित स्थिरीकरण तंत्र उत्पन्न करने के लिए सतत प्रक्षेप आवश्यक हैं। यह सूत्रीकरण उन अनुप्रयोगों के लिए विशेष रूप से उपयुक्त है जहाँ संपर्क सीमाएँ अज्ञात या विकसित होती हैं, जैसे कि संपर्क-सहायता प्राप्त तंत्र (contact-aided mechanisms), सॉफ्ट रोबोटिक्स, और टोपोलॉजी अनुकूलन, जहाँ निरंतरता-आधारित दृष्टिकोण स्पष्ट संपर्क खोज एल्गोरिदम और सक्रिय सेट रणनीतियों की आवश्यकता को समाप्त कर देता है। कार्य यह निष्कर्ष निकालता है कि जबकि यह विधि वैश्विक मजबूती के लिए एक बड़े स्थानीय बीजगणितीय सिस्टम का व्यापार करती है, यह उच्च-क्रम के तत्वों या विशिष्ट मेश प्रौद्योगिकियों की जटिलता के बिना गंभीर संपर्क परिदृश्यों के अनुकरण के लिए एक व्यवहार्य मार्ग प्रदान करती है।