Stretching and bending of (really) thick elastic plates

तकनीकी सारांश: (वास्तव में) मोटे इलास्टिक प्लेटों का खिंचाव और झुकना

समस्या विवरण
एक इलास्टिक प्लेट की यांत्रिक ऊर्जा को पारंपरिक रूप से खिंचाव (stretching) और झुकने (bending) के अलग-अलग घटकों में विभाजित समझा जाता है। यह पृथक्करण अत्यंत पतली प्लेटों ($h = O(\varepsilon)$ जहाँ $\varepsilon \ll 1$) के लिए एक कठोर परिणाम है, जिसे उन स्केलिंग धारणाओं से प्राप्त किया गया है जहाँ मिडलाइन खिंचाव छोटा ($E = O(\varepsilon)$) और झुकना परिमित ($K = O(1)$) होता है। हालांकि, यह "पतली-प्लेट चित्र" (thin-plate picture) मोटी प्लेटों के लिए भी आश्चर्यजनक रूप से सटीक सन्निकटन प्रदान करता है, यहाँ तक कि उन मामलों में भी जहाँ मोटाई $h = O(1)$ (क्रम इकाई) है और किरचॉफ परिकल्पना (कि अभिलंब (normals) अभिलंब ही रहते हैं) विफल हो जाती है। यह कार्य जिस केंद्रीय पहेली को संबोधित करता है वह यह है कि क्या छोटे विरूपणों (deformations) के तहत "वास्तव में मोटी" प्लेटों के लिए ऊर्जा का यह विभाजन (खिंचाव और झुकने के मोड में) औपचारिक रूप से उभरता है, और यदि ऐसा है, तो प्रभावी मापांक (effective moduli) कैसे व्यवहार करते हैं।

कार्यप्रणाली
लेखक एक दो-आयामी, अघुलनशील (incompressible), नियो-हुकियन (neo-Hookean) प्लेट के लिए इलास्टिसिटी के समीकरणों को हल करके इस समस्या की जांच करते हैं, जिसकी मोटाई $h = O(1)$ है। वे प्लेट के मिडलाइन पर एक अत्यंत लघु विरूपण लागू करते हैं, जिसे एक खिंचाव स्ट्रेन $E = O(\varepsilon)$ और एक वक्रता (curvature) $K = O(\varepsilon)$ द्वारा परिभाषित किया गया है।

1. काइनेमैटिक्स और कॉन्स्टिट्यूटिव संबंध: प्लेट को कार्तीय निर्देशांकों $(x, z)$ का उपयोग करके वर्णित किया गया है जहाँ $z \in [-h/2, h/2]$ है। विरूपण को मिडलाइन के विस्थापन और मोटाई के माध्यम से बिंदुओं के विस्थापन द्वारा निर्धारित किया जाता है। सामग्री को अघुलनशील नियो-हुकियन के रूप में मॉडल किया गया है जिसका इलास्टिक ऊर्जा घनत्व $\text{tr}(F^T F) - 2$ के समानुपाती है।
2. एसिम्प्टोटिक विस्तार (Asymptotic Expansion): लेखक छोटे पैरामीटर $\varepsilon$ के संदर्भ में एक नियमित क्रमबद्ध विस्तार (regular perturbation expansion) करते हैं। वे विस्थापन क्षेत्रों ($\tilde{u}, \tilde{w}$) और दबाव $p$ को $\varepsilon$ की घातों (powers) में विस्तारित करते हैं।
3. सीमा स्थितियाँ (Boundary Conditions): लेखक ऊपरी और निचले सतहों ($z = \pm h/2$) पर बल-मुक्त सीमा स्थितियाँ लागू करते हैं और मिडलाइन ($z=0$) पर एक निर्धारित विरूपण लागू करते हैं। पतली-प्लेट सिद्धांतों के विपरीत, जहाँ मिडलाइन निरंतरता (continuity) मानी जाती है, लेखक सतहों पर सीमा स्थितियों को संतुष्ट करने के लिए मिडलाइन ($z=0$) के पार समाधान में स्पष्ट रूप से विच्छेदन (discontinuities) की अनुमति देते हैं, यह स्वीकार करते हुए कि किरचॉफ परिकल्पना मोटी प्लेटों के लिए मान्य नहीं है।
4. समाधान रणनीति: अग्रणी-क्रम (leading-order) की रैखिककृत समस्या को फूरियर मोड अपघटन ($e^{ikx}$) का उपयोग करके हल किया जाता है। यह आंशिक अंतर समीकरणों (PDEs) को $z$ में साधारण अंतर समीकरणों (ODEs) के एक निकाय में बदल देता है, जिन्हें विश्लेष的に हल किया जाता है। परिणामों को FEniCSx में किए गए परिमित-तत्व सिमुलेशन (finite-element simulations) के विरुद्ध सत्यापित किया जाता है।

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

• ऊर्जा का औपचारिक पृथक्करण: गणनाएँ प्रदर्शित करती हैं कि एक मोटी प्लेट ($h=O(1)$) के लिए, जो लघु मिडलाइन विरूपणों से गुजर रही है, कुल इलास्टिक ऊर्जा घनत्व अभी भी एक खिंचाव पद और एक झुकने वाले पद में विभाजित होता है। ऊर्जा का रूप इस प्रकार है:
  $$E \propto \sum_k \left( h A(kh) E_k E_{-k} + h^3 B(kh) K_k K_{-k} \right)$$
  जहाँ $E_k$ और $K_k$ क्रमशः खिंचाव और झुकने के स्ट्रेन के फूरियर आयाम हैं।
• प्रभावी मापांक (Effective Moduli): लेखक प्रभावी खिंचाव ($A$) और झुकने ($B$) मापांक प्राप्त करते हैं जो आयामहीन पैरामीटर $\eta = kh$ (विकृति की तरंग दैर्ध्य के अनुपात में प्लेट की मोटाई) पर निर्भर करते हैं।
  • जैसे-जैसे $\eta \to 0$ (पतली-प्लेट सीमा), $A \to 2$ और $B \to 1/6$ होता है, जो शास्त्रीय पतली-प्लेट परिणामों को पुनः प्राप्त करता है।
  • जैसे-जैसे $\eta$ बढ़ता है (मोटी प्लेटें या छोटी तरंग दैर्ध्य), $A$ और $B$ दोनों कम हो जाते हैं, जो यह दर्शाता है कि प्लेट खिंचाव और झुकने दोनों में नरम (soften) हो जाती है।
• मापांक का स्थिर अनुपात: एक उल्लेखनीय परिणाम यह है कि खिंचाव और झुकने के योगदानों का अनुपात, जिसे "स्ट्रैचिंग नंबर" $S = A/12B$ द्वारा क्वांटिफाई किया जाता है, $\eta \lesssim 1$ के लिए लगभग स्थिर रहता है। विशेष रूप से, $S - 1 \approx \eta^2/30$। इसका तात्पर्य यह है कि जबकि मोटाई के साथ प्लेट की पूर्ण कठोरता (absolute stiffness) कम हो जाती है, खिंचाव बनाम झुकने की सापेक्ष ऊर्जा लागत एक पतली प्लेट के लगभग समान रहती है।
• सत्यापन: एसिम्प्टोटिक समाधान परिमित-तत्व सिमुलेशन के साथ उत्कृष्ट सहमति दिखाता है, जो प्राप्त प्रभावी मापांक की वैधता और मिडलाइन के पार विच्छेदन की उपस्थिति की पुष्टि करता है।

महत्व
यह शोध पत्र पतली-प्लेट सिद्धांतों द्वारा मोटी प्लेटों (जैसे जैविक कोशिका शीट) के यांत्रिकी को सटीक रूप से वर्णित करने के अनुभवजन्य अवलोकन को एक सैद्धांतिक आधार प्रदान करता है। लेखक तर्क देते हैं कि इन सिद्धांतों की आश्चर्यजनक सटीकता इस तथ्य से आती है कि मोटा होने पर भी ऊर्जा का खिंचाव और झुकने में विभाजन बना रहता है। इसके अलावा, प्लेट की मोटाई के तुल्य लंबाई पैमानों तक खिंचाव बनाम झुकने का सापेक्ष महत्व संरक्षित रहता है।

परिणामस्वरूप, मोटी इलास्टिक प्लेटों के यांत्रिक व्यवहार को पतली-प्लेट सिद्धांतों द्वारा पकड़ा जा सकता है, बशर्ते कि प्रभावी मापांक के नरम होने को ध्यान में रखते हुए बाहरी बलों को पुन: स्केल (rescale) किया जाए। यह निष्कर्ष मॉर्फोइलास्टिसिटी (morphoelasticity) और जैविक प्रणालियों (जैसे Volvox कोशिका शीट) के अध्ययन के लिए विशेष रूप से प्रासंगिक है जहाँ सामग्री मोटी होती है और विकास-प्रेरित विरूपण होते हैं। यह कार्य सुझाव देता है कि आंतरिक खिंचाव और झुकने का वैचारिक ढांचा, जो व्यापक रूप से पतली-शेल सिद्धांतों में उपयोग किया जाता है, मोटी जैविक संरचनाओं के लिए भी लागू रहता है।