Computation of stresses in jammed packings modeled with Tresca friction

तकनीकी सारांश: ट्रेस्का घर्षण (Tresca Friction) के साथ मॉडल किए गए जैम्ड पैकिंग्स (Jammed Packings) में प्रतिबल (Stresses) की गणना

समस्या विवरण
यह शोध पत्र कठोर बहुभुजीय कोशिकाओं (rigid polygonal cells), जैसे कि चिनाई वाली संरचनाओं या सघन दानेदार सामग्रियों (dense granular materials) के जैम्ड पैकिंग्स के भीतर स्ट्रेस फील्ड के कंप्यूटेशन को संबोधित करता है। पारंपरिक डिस्क्रीट एलीमेंट मेथड्स (DEM) के विपरीत, जो अक्सर अंतर्वेशन (interpenetration) को दंडित करने के लिए लोचदार विरूपण (elastic deformations) पर निर्भर करते हैं, यह कार्य उन प्रणालियों पर केंद्रित है जहाँ कोशिकाएं पूरी तरह से कठोर होती हैं और लोचदार रूप से विरूपित नहीं होती हैं। ऐसी "जैम्ड" (jammed) विन्यासों में, यांत्रिक संतुलन केवल संपर्क बलों और घर्षण के माध्यम से बनाए रखा जाता है, बिना किसी महत्वपूर्ण लोचदार विकृति (elastic strain) के। विशिष्ट चुनौती इन हाइपरस्टेटिक (hyperstatic) प्रणालियों में आंतरिक स्ट्रेस वितरण को निर्धारित करने में निहित है, जहाँ बल चरों (force variables) की संख्या संतुलन बाधाओं (equilibrium constraints) की संख्या से अधिक होती है, जिससे एक गैर-अद्वितीय समाधान स्थान (non-unique solution space) उत्पन्न होता है। लेखक संपर्क भौतिकी को एक ट्रेस्का घर्षण नियम (Tresca friction law) का उपयोग करके मॉडल करते हैं, जो एक निश्चित स्पर्शरेखीय फिसलन (tangential slip) की दहलीज (threshold) को प्रस्तावित करता है, जो तनाव-निर्भर कूलम्ब घर्षण (Coulomb friction) से भिन्न है।

कार्यप्रणाली (Methodology)
प्रस्तावित कार्यप्रणाली एक तीन-चरणीय कम्प्यूटेशनल ढांचे का अनुसरण करती है:

1. प्राइमल ऊर्जा न्यूनीकरण (Primal Energy Minimization):
   लेखक एक बाधित न्यूनीकरण समस्या (constrained minimization problem) तैयार करते हैं जहाँ कुल घर्षण ऊर्जा को न्यूनतम किया जाता है। ऊर्जा फलन (energy functional) सभी आंतरिक इंटरफेस पर ट्रेस्का घर्षण कार्य का योग है, जिसे $s_T |[u]_e \cdot t_e|$ के रूप में परिभाषित किया गया है, जहाँ $s_T$ घर्षण गुणांक है, $[u]_e$ विस्थापन जंप (displacement jump) है, और $t_e$ स्पर्शरेखीय वेक्टर (tangent vector) है। यह न्यूनीकरण निम्नलिखित के अधीन है:

• गैर-अंतर्वेशन बाधाएं (Non-interpenetration constraints): विस्थापन जंप का अभिलंब घटक गैर-ऋणात्मक होना चाहिए ($[u]_e \cdot n_e \geq 0$)।
• परिरोध (Confinement): एक न्यूमैन सीमा स्थिति ($g$) लागू की जाती है ताकि कोशिकाओं को अनंत तक फिसलने से रोका जा सके, क्योंकि ऊर्जा में लोचदार पुनर्स्थापना बल (elastic restoring forces) का अभाव है।
• स्थानांतरण अपरिवर्तनीयता (Translation invariance): द्रव्यमान केंद्र (center of mass) को स्थिर करने के लिए एक वैश्विक बाधा जोड़ी गई है।
  कोशिकाओं को केवल स्थानांतरीय स्वतंत्रता की डिग्री (translational degrees of freedom) वाले कठोर पिंडों के रूप में मॉडल किया गया है; घूर्णी स्वतंत्रता की डिग्री (rotational degrees of freedom) को स्पष्ट रूप से उपेक्षित किया गया है। यह सरलीकरण इस अवलोकन से प्रेरित है कि जैम्ड विन्यासों में, घूर्णन के प्रतिरोध (जिसके लिए डायलेटेंसी/dilatancy की आवश्यकता होती है) की तुलना में स्लाइडिंग के प्रतिरोध की दर कई गुना कम होती है।

2. ड्यूल मैक्सिमाइजेशन (Dual Maximization):
   चूंकि प्राइमल समस्या विस्थापनों (जो स्थिर जैम्ड अवस्था में प्रभावी रूप से शून्य होते हैं) की खोज करती है, इसलिए लेखक एक ड्यूल मैक्सिमाइजेशन समस्या का निष्कर्ष निकालते हैं। ड्यूल चर आंतरिक किनारों पर कार्य करने वाले बलों के अनुरूप होते हैं। इस ड्यूल समस्या को हल करने से इंटरफेस पर अभिलंब और स्पर्शरेखीय स्ट्रेस प्राप्त होते हैं। परिरोध शब्द के संबंध में उपयुक्त परिकल्पनाओं के तहत प्राइमल ऊर्जा के न्यूनीकरणकर्ताओं (minimizers) और ड्यूल समस्या के समाधानों का अस्तित्व स्थापित किया गया है। ड्यूल फॉर्मूलेशन प्रभावी रूप से हाइपरस्टेटिक मैनिफोल्ड के संभावित स्ट्रेस स्टेट्स में से एक विशिष्ट समाधान का चयन करता है।

3. स्ट्रेस पुनर्निर्माण (Stress Reconstruction):
   ड्यूल समस्या का आउटपुट किनारों पर असतत स्ट्रेस मान प्रदान करता है। प्रत्येक कोशिका के भीतर एक निरंतर स्ट्रेस फील्ड $\sigma$ को पुनर्गठित करने के लिए, लेखक लो-ऑर्डर रेविएर-थॉमस ($RT_0$) परिमित तत्वों (finite elements) का उपयोग करते हैं। ये तत्व $H(\text{div})$-कन्फॉर्मिंग हैं, जिसका अर्थ है कि उनके स्वतंत्रता की डिग्री किनारों के बैरीसेंटर (barycenters) पर अभिलंब स्ट्रेस के अनुरूप होती है। पुनर्निर्माण में प्रत्येक कोशिका में एक स्थानीय मिश्रित समस्या (local mixed problem) को हल करना शामिल है ताकि एक स्ट्रेस टेंसर $\sigma_c$ प्राप्त किया जा सके जो:

• $\text{div}(\sigma_c) = 0$ (आयतन भार की अनुपस्थिति में संतुलन) सुनिश्चित करे।
• सीमा पर अभिलंब स्ट्रेस, ड्यूल समाधान से प्राप्त मानों से मेल खाता हो।
• स्ट्रेस टेंसर की समरूपता ($\sigma = \sigma^T$) को एक लैग्रेंज मल्टीप्लायर के माध्यम से दुर्बल रूप से (weakly) लागू किया जाए। यह आवश्यक है क्योंकि अंतर्निहित असतत मॉडल मोमेंट बैलेंस समीकरणों की उपेक्षा करता है, जिसका अर्थ है कि पुनर्निर्मित स्ट्रेस में एक स्क्यू-सिमेट्रिक (skew-symmetric) घटक हो सकता है जो अवशिष्ट टॉर्क (residual torques) को दर्शाता है।

प्रमुख योगदान (Key Contributions)

• रिजिड फ्रिक्शन मॉडल: यह शोध पत्र एक मॉडल पेश करता है जो लोचदार मापदंडों के बिना केवल रिजिड फ्रिक्शन (ट्रेस्का मानदंड) पर निर्भर करता है, जो इसे मानक DEM दृष्टिकोणों से अलग करता है।
• हाइपरस्टैटिसिटी हैंडलिंग: यह कार्य जैम्ड पैकिंग की स्टैटिक अनिश्चितता (static indeterminacy) को स्पष्ट रूप से संबोधित करता है। यह प्रदर्शित करता है कि अनुकूलन दृष्टिकोण घर्षण ऊर्जा फलन को न्यूनतम करके स्व-तनाव (self-stress) अवस्थाओं के स्थान से एक अद्वितीय समाधान चुनता है।
• सुसंगत स्ट्रेस पुनर्निर्माण: एक नवीन पुनर्निर्माण तकनीक प्रस्तुत की गई है जो सेल के भीतर स्थानीय बल संतुलन सुनिश्चित करने के लिए $RT_0$ तत्वों का उपयोग करती है।
• घूर्णन की उपेक्षा: लेखक एक सरलीकृत मॉडल का प्रस्ताव करते हैं जो घूर्णी स्वतंत्रता की डिग्री को अनदेखा करता है, यह तर्क देते हुए कि यह संकुचित पैकिंग के लिए भौतिक रूप से उचित है जहाँ स्थानांतरीय स्लाइडिंग हावी होती है, हालांकि वे स्वीकार करते हैं कि इससे गैर-सममित स्ट्रेस टेंसर प्राप्त होते हैं।

संख्यात्मक परिणाम (Numerical Results)
कार्यप्रणाली का परीक्षण विभिन्न 2D बहुभुजीय मेश (विशेष रूप से वोरोनोई टेसेलेशन) पर किया गया, जिसके निम्नलिखित निष्कर्ष निकले:

• सत्यापन (Verification): समान सीमा स्थितियों (uniform boundary conditions) के परीक्षणों में, पुनर्निर्मित स्ट्रेस फील्ड अपेक्षित समान स्ट्रेस अवस्था से मेल खाया, जो विधि की निरंतरता की पुष्टि करता है।
• शियर टेस्ट (Shear Test): एक शियर लोडिंग परिदृश्य में, पुनर्निर्मित स्ट्रेस टेंसर असममित ($\sigma_{12} \neq \sigma_{21}$) था। लेखक इस विषमता का श्रेय असतत मॉडल में मोमेंट बैलेंस को लागू न करने को देते हैं, यह नोट करते हुए कि स्क्यू-सिमेट्रिक भाग अवशिष्ट टॉर्क को दर्शाता है।
• समांगी संपीड़न (Homogeneous Compression): समान संपीड़न सीमा स्थितियों के बावजूद, सॉल्वर ने एक गैर-समांगी स्ट्रेस फील्ड उत्पन्न किया। लेखक इसे हाइपरस्टैटिसिटी के परिणाम के रूप में समझाते हैं; अनुकूलन एक विशिष्ट ऊर्जा-न्यूनतम अवस्था का चयन करता है, जो स्वाभाविक रूप से "फोर्स चेन्स" (उच्च स्ट्रेस के विरल नेटवर्क) और निम्न या शून्य स्ट्रेस वाले क्षेत्रों के निर्माण की ओर ले जाता है।
• संपीड़न में तनाव (Tension in Compression): परिणामों ने वैश्विक संपीड़न के तहत भी कोशिकाओं में तनाव (धनात्मक मुख्य स्ट्रेस) दिखाया। लेखक इसे एक भौतिक घटना के रूप में पहचानते हैं जो घर्षण इंटरलॉकिंग और ब्रिजिंग प्रभावों से उत्पन्न होती है, न कि किसी संख्यात्मक त्रुटि (numerical artifact) से।

महत्व और दावे (Significance and Claims)
शोध पत्र का दावा है कि इसकी प्रस्तावित कार्यप्रणाली रिजिड फ्रिक्शन के भौतिकी और ऊर्जा न्यूनीकरण का लाभ उठाकर जैम्ड पैकिंग में स्ट्रेस की गणना करने का एक मजबूत और सुसंगत तरीका प्रदान करती है। लेखक तर्क देते हैं कि उनका दृष्टिकोण लोचदार-आधारित DEM की तुलना में जैम्ड सिस्टम के अंतर्निहित भौतिकी को बेहतर ढंग से पकड़ता है, विशेष रूप से उन क्षेत्रों में जहाँ लोचदार विरूपण नगण्य है।

परिणामों का महत्व इस तथ्य के प्रदर्शन में निहित है कि घर्षण ऊर्जा को न्यूनतम करना (एक $L_1$-प्रकार का ऑब्जेक्टिव) स्वाभाविक रूप से विरल बल नेटवर्क (sparse force networks) की ओर ले जाता है, जो ग्रेनुलर सामग्री में फोर्स चेन्स के प्रायोगिक अवलोकनों के अनुरूप है। शोध पत्र अपनी सीमाओं को विनम्रतापूर्वक स्वीकार करता है: घूर्नात्मक संतुलन की उपेक्षा के कारण गैर-सममित स्ट्रेस टेंसर प्राप्त होते हैं, और ट्रेस्का घर्षण (निश्चित दहलीज) का उपयोग कूलम्ब घर्षण के लिए आवश्यक अधिक जटिल, गैर-कन्वेक्स क्वासी-वैरिएशनल इनइक्वेलिटीज (Quasi-Variational Inequalities) की तुलना में समस्या को सरल बनाता है। लेखक सुझाव देते हैं कि भविष्य के कार्य में पूर्ण यांत्रिक निरंतरता प्राप्त करने के लिए घूर्णी स्वतंत्रता की डिग्री और अधिक यथार्थवादी घर्षण नियमों को शामिल करने के लिए मॉडल का विस्तार किया जा सकता है।