Caveats on formulating finite elasto-plasticity in curvilinear coordinates

यहाँ Prett et al. के शोध पत्र "Caveats on formulating finite elasto-plasticity in curvilinear coordinates" का विस्तृत तकनीकी सारांश दिया गया है।

1. समस्या विवरण (Problem Statement)

यह शोध पत्र वक्राकार निर्देशांक प्रणालियों (विशेष रूप से अक्षीय सममिति/axisymmetry) के भीतर परिमित-विकृति प्रत्यास्थ-प्लास्टिसिटी (finite-strain elasto-plasticity) के संख्यात्मक कार्यान्वयन में एक महत्वपूर्ण अंतराल को संबोधित करता है। जबकि निरंतरता यांत्रिकी (continuum mechanics) मौलिक रूप से टेंसर विश्लेषण पर आधारित है, जो फ्रेम-इनवेरिएंट (frame-invariant) है, परिमित तत्व (Finite Element - FE) कार्यान्वयन उपयोगकर्ता द्वारा चयनित आधारों (bases) में मैट्रिक्स निरूपण पर निर्भर करते हैं।

• मुख्य मुद्दा: जब कार्टेशियन (ऑर्थोनॉर्मल) आधारों से वक्राकार (गैर-ऑर्थोनॉर्मल, जैसे बेलनाकार) आधारों में संक्रमण किया जाता है, तो मानक सूत्रण अक्सर विफल हो जाते हैं या त्रुटियां उत्पन्न करते हैं। विशेष रूप से, विकृति प्रवणता (deformation gradient), जैकोबियन (Jacobian), और तनाव मापों (stress measures) के लिए कार्टेशियन परिभाषाओं का सहज विस्तार आयतन रूपांतरण और गतिक विवरणों में अशुद्धियों का कारण बनता है।
• चुनौती: अक्षीय सममिति वाली समस्याओं में (जो भू-तकनीकी इंजीनियरिंग, धातु निर्माण और बायोमैकेनिक्स में सामान्य हैं), ये त्रुटियां बड़े विरूपण और अनैलिक प्रभाव (anelastic effects) (प्लास्टिसिटी) के समावेश के कारण और भी बढ़ जाती हैं। शोध पत्र इस बात पर प्रकाश डालता है कि परिमित-विकृति प्लास्टिसिटी में कॉची स्ट्रेस (Cauchy stress), वर्तमान इलास्टिक अवस्था से परे विन्यास परिवर्तनों पर निर्भर करता है, जिसके लिए जैकबियन और विकृति प्रवणता के इलास्टिक और प्लास्टिक घटकों के बीच सावधानीपूर्वक अंतर करना आवश्यक है।
• लक्ष्य: एक सुव्यवस्थित, चरण-दर-चरण प्रक्रिया प्रदान करना जो अक्षीय सममित परिमित तत्व विश्लेषण को सही ढंग से तैयार कर सके, जो स्पष्ट आधार परिवर्तन (explicit basis changes) का उपयोग करके मानक कार्टेशियन ढांचे के भीतर इन ज्यामितीय जटिलताओं को संभाल सके, बिना अमूर्त रिमानियन मैनिफोल्ड सिद्धांत (Riemannian manifold theory) का सहारा लिए।

2. कार्यप्रणाली (Methodology)

लेखक ने गतिक विसंगतियों को उजागर करने और उन्हें हल करने के लिए एक कैविटी एक्सपेंशन समस्या (एक खोखले सिलेंडर का विस्तार) को एक बेंचमार्क मॉडल के रूप में उपयोग किया है। कार्यप्रणाली में शामिल हैं:

• मिश्रित आधारों में टेंसर विश्लेषण (Tensor Analysis in Mixed Bases): शोध पत्र मानक कार्टेशियन निरूपणों के भीतर कार्य करता है लेकिन वक्राकार (बेलनाकार) कोवेरिएंट (covariant) और कंट्रावेरिएंट (contravariant) आधारों में रूपांतरण नियमों को स्पष्ट रूप से व्युत्पन्न करता है। यह स्पष्ट आधार परिवर्तन मैट्रिसेस का उपयोग करके विभेदक ज्यामिति मैनिफोल्ड के बजाय प्रत्यक्ष परिवर्तन पर ध्यान केंद्रित करता है।
• गतिक परिभाषाएं (Kinematic Definitions):
  • विकृति प्रवणता ($F$): इसे सामग्री वक्रों (material curves) को मैप करने वाले एक टू-पॉइंट टेंसर के रूप में कड़ाई से परिभाषित किया गया है। शोध पत्र स्पष्ट करता है कि $F$ केवल वक्राकार निर्देशांकों में एक वेक्टर फील्ड का ग्रेडिएंट नहीं है और इसके लिए अलग-अलग विन्यासों के बीच मानचित्रण करने के लिए शिफ्टर टेंसर (shifter tensor, $S$) की आवश्यकता होती है।
  • जैकोबियन ($J$): यह दिखाया गया है कि आयतन अनुपात मेट्रिक गुणांकों के सारणिक (determinant) पर निर्भर करता है। लेखक $J = \det(F) \sqrt{\det(g)}/\sqrt{\det(G)}$ के संबंध को व्युत्पन्न करते हैं, यह सिद्ध करते हुए कि $\det(F_{curvilinear}) \neq J$ सीधे तौर पर।
  • इलास्टो-प्लास्टिक अपघटन (Elasto-Plastic Decomposition): विकृति प्रवणता का गुणात्मक अपघटन ($F = F^e F^p$) किया जाता है। शोध पत्र मध्यवर्ती तनाव-मुक्त विन्यास के मेट्रिक को परिभाषित करने की कठिनाई को संबोधित करता है। वे इलास्टिक जैकबियन ($J_e$) के लिए एक व्यावहारिक रैखिकीकरण प्रस्तावित करते हैं जो इलास्टिक लॉगरिदमिक स्ट्रेन के ट्रेस ($J_e = \exp(\text{tr}(\epsilon^e))$) का उपयोग करता है, जिससे सुसंगत रैखिकीकरण सरल हो जाता है।
• तनाव माप (Stress Measures): शोध पत्र एस्लेबी-जीटा स्ट्रेस (Eshelby-zeta stress, $\zeta = J_e \sigma$) को एक पसंदीदा तनाव माप के रूप में पेश करता है, जो संबद्ध प्रवाह नियमों (associated flow rules) के साथ परिमित-विकृति प्लास्टिसिटी के लिए थर्मोडायनामिक स्थिरता और तनाव-मुक्त मध्यवर्ती विन्यास को सुनिश्चित करता है।
• वीक फॉर्म्स और विविक्तकरण (Weak Forms & Discretization): मोमेंटम बैलेंस के स्ट्रॉन्ग फॉर्म को टोटल लैग्रेंजियन (TL) और अपडेटेड लैग्रेंजियन (UL) दोनों के लिए वीक फॉर्म में परिवर्तित किया जाता है। इन्हें $2\pi r$ आयतन/सतह पदों को सावधानीपूर्वक संभालते हुए, परिधि निर्देशांक को एकीकृत करके 2D अक्षीय सममित रूपों में कम किया जाता है।
• रैखिकीकरण (Linearization): एक सुसंगत न्यूटन-राफसन रैखिकीकरण व्युत्पन्न किया गया है, जो स्पष्ट रूप से कोवेरिएंट डेरिवेटिव और क्रिस्टोफेल प्रतीकों (Christoffel symbols) को ध्यान में रखता है, जो वक्राकार निर्देशांकों में गैर-शून्य होते हैं।

3. मुख्य योगदान (Key Contributions)

1. शिफ्टर टेंसर का स्पष्टीकरण: शोध पत्र संदर्भ और वर्तमान वक्राकार आधारों के बीच विस्थापन घटकों को मैप करने के लिए शिफ्टर टेंसर की भूमिका को स्पष्ट रूप से परिभाषित करता है, जिससे उस सामान्य त्रुटि को सुधारा जा सके जहाँ विस्थापनों को केवल साधारण वेक्टर अंतर के रूप में माना जाता है।
2. सही जैकबियन सूत्र: यह प्रदर्शित करता है कि वक्राकार घटकों में विकृति प्रgrenता का सारणिक आयतन अनुपात ($J$) के बराबर नहीं होता है, जब तक कि इसे मेट्रिक सारणिकों द्वारा सुधारा न जाए। यह $J$ और इसके इलास्टिक/प्लास्टिक अपघटन के लिए सही सूत्र प्रदान करता है।
3. व्यावहारिक रैखिकीकरण रणनीति: लेखक एक विशिष्ट रैखिकीकरण ($J_e = \exp(\text{tr}(\epsilon^e))$) का प्रस्ताव करते हैं जो अतिरिक्त पदों को न्यूनतम करता है और मानक कॉन्स्टिट्यूटिव सबरूटीन के साथ सहजता से एकीकृत होता है, जिससे प्लास्टिक विकृति प्रवणता के जटिल डेरिवेटिव की गणना करने की आवश्यकता समाप्त हो जाती है।
4. एल्गोरिदम वर्कफ़्लो: शोध पत्र UL और TL दोनों अक्षीय सममित FE कार्यान्वयन के लिए विस्तृत एल्गोरिदम (एल्गोरिदम 1 और 2) प्रदान करता है, जो यह बताता है कि मानक कार्टेशियन कोड की तुलना में वक्राकार पदों (क्रिस्टोफेल प्रतीक, मेट्रिक अपडेट) को कहाँ डाला जाना चाहिए।
5. सत्यापन: यह एक पूर्ण 3D कार्टेशियन कार्यान्वयन के समान समस्या के लिए प्रस्तावित 2D अक्षीय सममित कोड के विरुद्ध अपने सूत्र का सत्यापन करता है।

4. परिणाम (Results)

लेखकों ने एक मोटी दीवार वाले सिलेंडर (thick-walled cylinder) को ऊर्ध्वाधर ऊंचाई के साथ रैखिक रूप से बदलते विस्थापन के तहत आंतरिक रेडियल विस्तार के अधीन रखकर अपने सूत्र को सत्यापित किया।

• तुलना: प्रस्तावित 2D अक्षीय सममित सूत्र के परिणामों की तुलना 20-नोड वाले हेक्साहेड्रल तत्वों का उपयोग करके किए गए 3D परिमित तत्व विश्लेषण से की गई।
• सटीकता: प्रस्तावित 2D अक्षीय सममित सूत्र के परिणाम प्रमुख इनवेरिएंट्स (invariants) के लिए 3D विश्लेषण के समान थे:
  • नकारात्मक कॉची दबाव (Negative Cauchy pressure)।
  • डेविएटरिक कॉची स्ट्रेस (Deviatoric Cauchy stress)।
  • इलास्टिक और प्लास्टिक जैकबियन।
• अभिसरण (Convergence): न्यूटन-राफसन पुनरावृत्ति योजना ने मजबूत अभिसरण प्रदर्शित किया (आमतौर पर प्रति समय चरण 3-4 पुनरावृत्तियाँ) और दोनों 3D और 2D स्वरूपों के लिए समान अभिसरण दर देखी गई।
• कंप्यूटेशनल दक्षता: अक्षीय सममित सूत्र ने 3D जैसी सटीकता के साथ काफी कम कंप्यूटेशनल लागत (2D विविक्तकरण बनाम 3D) प्राप्त की, जो प्रस्तावित विधि की व्यवहार्यता की पुष्टि करती है।

5. महत्व (Significance)

यह शोध पत्र कम्प्यूटेशनल मैकेनिक्स समुदाय के लिए कई कारणों से महत्वपूर्ण है:

• त्रुटि निवारण: यह उन सूक्ष्म लेकिन मौलिक त्रुटियों की पहचान करता और सुधारता है जो वक्राकार निर्देशांकों में परिमित-विकृति प्लास्टिसिटी लागू करने के दौरान उत्पन्न होती हैं। गतिक मात्राओं की गलत परिभाषाएँ असंगत टेंगेंट स्टिफनेस मैट्रिक्स, अभिसरण की कमी, या भौतिक रूप से गलत तनाव अवस्थाओं का कारण बन सकती हैं।
• सिद्धांत और अभ्यास के बीच सेतु: अमूर्त मैनिफोल्ड सिद्धांत के बजाय स्पष्ट आधार रूपांतरणों पर ध्यान केंद्रित करके, यह शोध पत्र एक "कुकबुक" दृष्टिकोण प्रदान करता है जो FE डेवलपर्स को गहन विभेदक ज्यामिति विशेषज्ञता की आवश्यकता के बिना मजबूत अक्षीय सममित सॉल्वर लागू करने की अनुमति देता है।
• जटिल सिमुलेशन को सक्षम बनाना: जैकबियन और तनाव मापों का कठोर उपचार बड़े विरूपण वाले जटिल अक्षीय सममित समस्याओं (जैसे मृदा प्रवेश/CPT, टनलिंग और धातु निर्माण) के विश्वसनीय सिमुलेशन को सक्षम बनाता है, जहाँ आयतन परिवर्तन और एनिसोट्रोपिक प्रभाव महत्वपूर्ण होते हैं।
• मानकीकरण: यह अक्षीय सममित FE कोड में "शिफ्टर", मेट्रिक गुणांक और कोवेरिएंट डेरिवेटिव को संभालने के लिए एक स्पष्ट, सुसंगत नोटेशन और प्रक्रिया स्थापित करता है, जिससे भविष्य के अनुसंधान में पुनरुत्पादकता और सटीकता को बढ़ावा मिलता है।

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