A stochastic model for elastoplastic contact of rough surfaces incorporating scale-dependent hardness

तकनीकी सारांश: पैमाने-निर्भर कठोरता (Scale-Dependent Hardness) को समाहित करने वाले खुरदरे सतहों के इलास्टोप्लास्टिक संपर्क के लिए स्टोकेस्टिक मॉडल

समस्या विवरण
प्राकृतिक और निर्मित सतहों में अंतर्निहित खुरदरापन होता है जो संपर्क असपेरिटीज़ (asperities) पर तनाव सांद्रता (stress concentrations) की ओर ले जाता है, जिससे अक्सर प्लास्टिक विरूपण (plastic deformation) प्रेरित होता है। इस इलास्टोप्लास्टिक संपर्क का मॉडलिंग करना दो प्रतिस्पर्धी कारकों के कारण जटिल हो जाता है: कई पैमानों पर स्व-समान (self-affine) सतह स्थलाकृति की प्रकृति और प्लास्टिक विरूपण का "आकार प्रभाव" (size effect)। दूसरा कारक यह बताता है कि सामग्री की कठोरता एक अंतर्निहित स्थिरांक नहीं है, बल्कि इंडेंटेशन आकार घटने (या आवर्धन बढ़ने) के साथ ज्यामितीय रूप से आवश्यक डिस्लोकेशन्स (Geometrically Necessary Dislocations - GNDs) के संचय के कारण बढ़ती है। मौजूदा मॉडल, जैसे कि पर्सन के डिफ्यूजन समीकरण या मल्टी-एस्पेरिटी दृष्टिकोण पर आधारित मॉडल, अक्सर स्थिर कठोरता मानते हैं या उच्च रिज़ॉल्यूशन बनाए रखते हुए पैमाने-निर्भर कठोरता को कुशलतापूर्वक शामिल करने में संघर्ष करते हैं। इसके अलावा, पिछले फॉर्मूलेशन अक्सर रैखिक समीकरणों (जैसे कि डिफ्यूजन-आधारित दृष्टिकोणों में ट्राई-डायगोनल मैट्रिसेस) के जटिल सिस्टम को हल करने पर निर्भर करते हैं, जो गणनात्मक रूप से गहन हो सकता है।

कार्यप्रणाली (Methodology)
लेखक पैमाने-निर्भर कठोरता वाले इलास्टोप्लास्टिक संपर्क समस्याओं को हल करने के लिए कंपाउंडेड चैपमैन-कोलमोगोरोव (C-K) समीकरण पर आधारित एक नवीन स्टोकेस्टिक फ्रेमवर्क विकसित करते हैं। मुख्य धारणाएं और चरण निम्नलिखित हैं:

1. मार्कोव प्रक्रिया धारणा: आवर्धन ($\zeta$) के सापेक्ष यादृच्छिक संपर्क दबाव (random contact pressure) के परिवर्तन को एक मार्कोव प्रक्रिया के रूप में माना गया है।
2. इंटीग्रल समीकरण फॉर्मूलेशन: डिफ्यूजन समीकरण का उपयोग करने के बजाय, लेखक तीन इंटीग्रल समीकरणों की एक प्रणाली विकसित करते हैं जो निम्नलिखित के विकास का वर्णन करती है:
   • इलास्टिक संपर्क दबाव का संभाव्यता घनत्व फलन (PDF), $P_0(p, \zeta)$।
   • सापेक्ष प्लास्टिक संपर्क क्षेत्र, $A^*_{pl}(\zeta)$।
   • सापेक्ष गैर-संपर्क क्षेत्र, $A^*_{nc}(\zeta)$।
3. पैमाने-निर्भर कठोरता: मॉडल एक कठोरता फलन $H(\zeta) = H_0 \zeta^n$ को शामिल करता है, जहाँ $n$ कठोरता घातांक (hardness exponent) है। यह यह अनुमति देता है कि अवलोकन पैमाने के घटने (आवर्धन बढ़ने) के साथ इलास्टिक दबाव की ऊपरी सीमा बढ़ जाए।
4. बैकफ्लो तंत्र (Backflow Mechanism): एक महत्वपूर्ण सैद्धांतिक प्रगति "नो री-एंट्री" (no re-entry) धारणा का त्याग है जो स्थिर-कठोरता वाले मॉडलों में पाई जाती है। जब कठोरता पैमाने के साथ बढ़ती है, तो पहले से प्लास्टिक संपर्क ($p = H(\zeta')$) में मौजूद सतह के बिंदु उच्च आवर्धन पर वापस इलास्टिक संपर्क क्षेत्र या यहाँ तक कि गैर-संपर्क क्षेत्र में भी परिवर्तित हो सकते हैं। इस "संभाव्यता के बैकफ्लो" को इंटीग्रल समीकरणों में गणितीय रूप से कैप्चर किया गया है, जो इस मॉडल को पिछले स्थिर-कठोरता वाले फॉर्मूलेशन से अलग करता है।
5. संख्यात्मक कार्यान्वयन (Numerical Implementation): आवर्धन अक्ष को लॉगरिदमिक स्केल का उपयोग करके विविक्त (discretize) किया गया है। समाधान आवर्धन चरणों के माध्यम से क्रमिक रूप से आगे बढ़ता है, जिसमें एब्जॉर्बिंग (absorbing) और नॉन-एब्जॉर्बिंग बाउंड्री कंडीशंस के साथ एक ट्रांजिशन प्रोबेबिलिटी कर्नेल का उपयोग किया जाता है।

प्रमुख परिणाम
अध्ययन विभिन्न स्थितियों के तहत मॉडल के व्यवहार को प्रदर्शित करने के लिए संख्यात्मक समाधान और अभिसरण परीक्षण (convergence tests) प्रस्तुत करता है:

• अभिसरण (Convergence): मॉडल विविक्त बिंदुओं ($n_\zeta$) के बढ़ने के साथ उत्कृष्ट अभिसरण दिखाता है, जिससे रैखिक इलास्टिक संपर्क के लिए पर्सन के सिद्धांत और स्थिर कठोरता के लिए मौजूदा क्लोज्ड-फॉर्म समाधानों के अनुरूप परिणाम प्राप्त होते हैं।
• संपर्क स्थिति का विकास:
  • कठोरता घातांक ($n$): $n$ (मजबूत आकार प्रभाव) बढ़ाने से संपर्क व्यवहार पूरी तरह से प्लास्टिक से रैखिक इलास्टिक की ओर स्थानांतरित हो जाता है। पर्याप्त बड़े $n$ के लिए, संपर्क पूर्ण नाममात्र संपर्क (full nominal contact) प्राप्त होने तक इलास्टिक बना रहता है।
  • हर्स्ट घांक (Hurst Exponent - $H$): उच्च $H$ मान (चिकनी सतहें) सतह की ऊंचाई के विचरण और दबाव पुनर्वितरण की दर को कम करके इलास्टिक व्यवहार को बढ़ावा देते हैं।
  • सिमेट्री ब्रेकिंग (Symmetry Breaking): स्थिर-कठोरता वाले मॉडलों के विपरीत, जहाँ इलास्टिक प्रेशर PDF $H_0/2$ के आसपास सममित होता है, पैमाने-निर्भर मॉडल इस समरूपता को तोड़ देता है। बैकफ्लो तंत्र के कारण वर्तमान कठोरता सीमा $H(\zeta)$ के पास संभाव्यता जमा हो जाती है।
• क्षेत्र-भार संबंध (Area-to-Load Relationship): मॉडल रैखिक इलास्टिक से इलास्टोप्लास्टिक व्यवहार और अंततः पूर्ण प्लास्टिसिटी तक एक सुचारू संक्रमण की भविष्यवाणी करता है। कम भार ($\kappa$) पर क्षेत्र-भार संबंध का ढलान $n$ के बढ़ने के साथ घटता है, जो पूर्णतः प्लास्टिक सीमा से रैखिक इलास्टिक सीमा की ओर संक्रमण करता है।
• टोपोग्राफिक यील्ड पैरामीटर: लेखक एक नया टोपोग्राफिक यील्ड पैरामीटर, $w(\zeta)$, प्रस्तावित करते हैं, जो इलास्टिक दबाव के पुनर्वितरण और कठोरता में वृद्धि के बीच प्रतिस्पर्धा को मापता है। पिछले मापदंडों के विपरीत, यह नया फॉर्मूलेशन स्पष्ट रूप से यंग्स मॉड्यूलस, कट-ऑफ वेवनंबर और PSD स्थिरांकों को शामिल करता है।
• संपर्क स्थिति आरेख (Contact Status Diagram): मिसफिट मापदंडों ($M_{le}$ और $M_{fp}$) का उपयोग करके, जो शुद्ध इलास्टिक और शुद्ध प्लास्टिक सीमाओं से विचलन को मात्रात्मक रूप से दर्शाते हैं, वे $(H, n)$ स्पेस में एक संपर्क स्थिति आरेख उत्पन्न करते हैं। यह आरेख रैखिक इलास्टिक, इलास्टोप्लास्टिक और पूर्णतः प्लास्टिक शासन (regimes) को स्पष्ट करता है। परिणाम बताते हैं कि सामान्य खुरदरी सतहों ($H \approx 0.8$) के लिए, एक छोटा आकार प्रभाव ($n > 0$) भी प्रमुख रूप से इलास्टिक व्यवहार की ओर ले जा सकता है, जो स्थिर-कठोरता वाले मॉडलों की भविष्यवाणियों के विपरीत है।

महत्व और दावे
यह शोध संपर्क यांत्रिकी (contact mechanics) के क्षेत्र में निम्नलिखित महत्वपूर्ण योगदान का दावा करता है:

1. नवीन ढांचा (Novel Framework): यह खुरदरी सतह संपर्क विश्लेषण के लिए कंपाउंडेड चैपमैन-कोल्मोगोरोव समीकरण के अनुप्रयोग का नेतृत्व करता है, जो एक ऐसा स्टोकेस्टिक प्रोसेस फ्रेमवर्क है जो डिफ्यूजन-आधारित या मल्टी-एस्पेरिटी समेशन विधियों से भिन्न है।
2. कंप्यूटेशनल दक्षता: इंटीग्रल समीकरण दृष्टिकोण डिफ्यूजन-आधारित मॉडलों के लिए आवश्यक बड़े रैखिक समीकरणों (ट्राई-डायगोनल मैट्रिसेस) को हल करने की आवश्यकता को समाप्त करता है, जो संभावित रूप से अधिक कुशल समाधान प्रक्रिया प्रदान करता है।
3. आकार प्रभाव में भौतिक अंतर्दृष्टि: मॉडल यह समझने के लिए एक तंत्र प्रदान करता है कि कैसे पैमाने-निर्भर कठोरता संपर्क यांत्रिकी को बदलती है, विशेष रूप से, प्लास्टिक से इलास्टिक अवस्थाओं में संभाव्यता के "बैकफ्लो" के माध्यम से, जो बड़े पैमाने (कम आवर्धन) पर अधिक इलास्टिक व्यवहार की ओर ले जाता है।
4. व्यावहारिक उपयोगिता: व्युत्पन्न संपर्क स्थिति आरेख सीमित भौतिक और ज्यामितीय गुणों के आधार पर संपर्क शासन (इलास्टिक बनाम प्लास्टिक) की तीव्र पहचान के लिए एक उपकरण प्रदान करता है, जो ट्राइबोलॉजिकल सिस्टम के डिजाइन और विश्लेषण में सहायता करता है।
5. व्यापक प्रयोज्यता: लेखक सुझाव देते हैं कि इंटरफेशियल गुणों के पैमाने के साथ विकास को वर्णित करने वाले इंटीग्रल समीकरण अन्य बहु-विषयक क्षेत्रों के लिए मूल्यवान अंतर्दृष्टि प्रदान कर सकते हैं, जिसमें भूकंप विज्ञान (earthquake mechanics), विद्युत संपर्क और संपर्क विद्युतीकरण (contact electrification) शामिल हैं।

अध्ययन निष्कर्ष निकालता है कि हालांकि इस मॉडल का आधार "प्रोट्रूबेंस-ऑन-प्रोट्रूबेंस" (protuberance-on-protuberance) अवधारणा (अर्चर्ड, गाओ-बावर, जैक्सन-स्ट्रिएटर) के समान है, लेकिन विशिष्ट ज्यामितीय प्रतिनिधित्व (स्केल-डिपेंडेंट वेरिएंस के साथ गॉसियन रफ सर्फेस) और C-K समीकरण के माध्यम से पैमाने-निर्भर कठोरता का समावेश इसे विशिष्ट गवर्निंग समीकरण और भौतिक भविष्यवाणियों की ओर ले जाता है।