🔬 physics

Size-dependent elastostatic behavior of non-centrosymmetric solids via a mixed variational finite element method

यह शोध पत्र गैर-केंद्रसममित (non-centrosymmetric) ठोसों के आकार-निर्भर इलास्टोस्टैटिक व्यवहार की जांच करने के लिए सुसंगत कपल स्ट्रेस सिद्धांत (consistent couple stress theory) पर आधारित एक नवीन मिश्रित विवरणात्मक परिमित तत्व विधि (mixed variational finite element method) प्रस्तुत करता है, जो सूक्ष्म और नैनो-पैमाने की प्रतिक्रियाओं में नई अंतर्दृष्टि प्रदान करता है और कपल स्ट्रेस मापदंडों को मापने के लिए संभावित प्रयोगों की पहचान करता है।

मूल लेखक: G. F. Dargush, B. T. Darrall, T. K. Scully, O. Lavan

प्रकाशित 2026-08-10

📖 6 मिनट में पढ़ें🧠 गहराई से पढ़ें

CC BY 4.0

मूल लेखक: G. F. Dargush, B. T. Darrall, T. K. Scully, O. Lavan

मूल पेपर CC BY 4.0 (https://creativecommons.org/licenses/by/4.0/) के तहत लाइसेंस किया गया है। ✨ नीचे दिए गए पेपर की यह व्याख्या AI से तैयार की गई है। इसे लेखकों ने न तो लिखा है, न इसका समर्थन किया है। तकनीकी सटीकता के लिए मूल पेपर देखें। पूरा डिस्क्लेमर पढ़ें

सामग्री विज्ञान (materials science) की दुनिया की कल्पना एक विशाल, अदृश्य लेगो (Lego) सेट के रूप में करें। लंबे समय तक, वैज्ञानिकों ने इन ब्लॉकों के व्यवहार के अपने मॉडलों को "क्लासिकल इलास्टिसिटी" (classical elasticity) का उपयोग करके बनाया, जो नियमों का एक ऐसा समूह है जो हर पदार्थ के टुकड़े को एक चिकने, समान पिंड (uniform blob) की तरह मानता है। इस पुराने दृष्टिकोण में, यदि आप मिट्टी के एक ब्लॉक को दबाते हैं, तो वह बस दब जाता है; यदि आप उसे खींचते हैं, तो वह खिंच जाता है। नियम कहते हैं कि आकार मायने नहीं रखता—मिट्टी का एक छोटा सा कण भी उसी सामग्री के एक विशाल पत्थर की तरह व्यवहार करना चाहिए, बस उसका पैमाना छोटा होता है।

लेकिन, इसमें एक पेंच है। कुछ सामग्रियां, जैसे आपके फोन के अंदर के क्रिस्टल या कुछ प्रकार की रेत, पूरी तरह से सममित (symmetrical) नहीं होती हैं। उनमें एक "हाथों कापन" (handedness) या एक विशिष्ट दिशा होती है, ठीक वैसे ही जैसे आपका बायां हाथ दाएं हाथ से अलग होता है। पुराने नियमों में, ये असममित सामग्रियां मुख्य रूप से एक ही ट्रिक के लिए प्रसिद्ध थीं: पीज़ोइलेक्ट्रिसिटी (piezoelectricity)। यह वह जादू है जहाँ एक क्रिस्टल को दबाने से बिजली पैदा होती है (जैसे लाइटर या माइक्रोफ़ोन में)। वैज्ञानिक जानते थे कि ये सामग्रियां विशेष हैं, लेकिन वे एक दूसरी, छिपी हुई ट्रिक को नहीं जानते थे। उन्हें संदेह था कि जब आप इन असममित क्रिस्टलों को दबाते हैं, तो वे इस तरह से भी मुड़ या घूम सकते हैं जो पूरी तरह से इस बात पर निर्भर करता है कि वे कितने छोटे हैं। यह ऐसा है जैसे यह महसूस करना कि एक छोटा, असममित लेगो ब्रिक केवल दबता नहीं है; बल्कि अगर आप उसे दबाते हैं, तो वह एक पत्ती की तरह मुड़ भी सकता है, लेकिन केवल तभी जब वह ब्रिक सूक्ष्म (microscopic) हो। इस "आकार-निर्भर" (size-dependent) व्यवहार को समझना इसलिए महत्वपूर्ण है क्योंकि जैसे-जैसे हम छोटी और छोटी मशीनें (नैनोटेक्नोलॉजी) बना रहे हैं, ये सूक्ष्म विचित्रताएं हमारे डिजाइनों को बना या बिगाड़ सकती हैं।

यह शोध पत्र "कंसिस्टेंट कपल स्ट्रेस थ्योरी" (consistent couple stress theory) नामक एक नए गणितीय लेंस का उपयोग करके उस लापता ट्रिक की गहराई में जाता है। इस थ्योरी को खेल के नियमों को अपग्रेड करने के रूप में सोचें। केवल यह ट्रैक करने के बजाय कि सामग्री कितनी खिंचती है (strain), नए नियम यह भी ट्रैक करते हैं कि वह कितनी मुड़ती या मुड़ती है (mean curvature)। लेखकों ने, जो इंजीनियरों और गणितज्ञों की एक टीम है, महसूस किया कि असममित क्रिस्टलों के लिए, खिंचाव और मुड़ने के बीच एक गुप्त हाथ मिलाप (secret handshake) होता है। जिस तरह से एक पीज़ोइलेक्ट्रिक क्रिस्टल को दबाने से बिजली पैदा होती है, उसी तरह एक असममित क्रिस्टल को दबाने से एक "कपल स्ट्रेस" (couple stress) पैदा हो सकता है—एक सूक्ष्म आंतरिक टॉर्क (torque) जो सामग्री को घुमाने की कोशिश करता है।

टीम ने केवल अनुमान नहीं लगाया; उन्होंने एक परिष्कृत कंप्यूटर सिमुलेशन, एक "मिक्स्ड वेरिएशनल फाइनाइट एलीमेंट मेथड" (mixed variational finite element method) बनाया, जो एक आभासी प्रयोगशाला के रूप में कार्य करता है। इसे एक उच्च-तकनीकी वीडियो गेम की तरह समझें जहाँ वे विभिन्न क्रिस्टलों (जैसे सिलिकॉन कार्बाइड और क्वार्ट्ज) से बनी छोटी बीम, बार और घुमावदार पुल बना सकते हैं और उन्हें धक्का दे सकते हैं, खींच सकते हैं और घुमा सकते हैं ताकि देखा जा सके कि क्या होता है। उन्होंने गणित के उपकरणों का एक विशेष सेट बनाया जिसने उन्हें सामान्य कम्प्यूटेशनल सिरदर्द में फंसे बिना इन समस्याओं को हल करने की अनुमति दी, जो विशेष रूप से इस बात पर केंद्रित था कि कैसे सामग्री का आंतरिक "हाथों कापन" (handedness) उसकी कठोरता (stiffness) को बदल देता है।

उन्होंने जो पाया वह बेहद दिलचस्प है। सबसे पहले, उन्होंने पुष्टि की कि जैसे-जैसे ये सामग्रियां छोटी होती जाती हैं, वे अधिक कठोर होती जाती हैं, जैसा कि पिछले सिद्धांतों ने सुझाव दिया था। लेकिन असली आश्चर्य विषमता (asymmetry) से आया। उनके सिमुलेशन में, जब उन्होंने हेक्सागोनल सिलिकॉन कार्बाइड से बनी एक छोटी बार पर दबाव डाला, तो वह केवल छोटी नहीं हुई; बल्कि वह बगल की ओर झुक गई! यह एक विशुद्ध रूप से रैखिक प्रभाव (linear effect) है, जिसका अर्थ है कि यह भार के तहत तुरंत होता है, बिना किसी क्रिटिकल "टिपिंग पॉइंट" के जैसे कि एक बकलिंग कॉलम में होता है। यह ऐसा है जैसे एक छोटे, असममित स्प्रिंग को नीचे से दबाने पर वह धक्का देने की दिशा के आधार पर बाईं या दाईं ओर मुड़ जाता है।

इस प्रभाव का आकार सामग्री के भीतर दो अदृश्य "लंबाई पैमानों" (length scales) पर निर्भर करता है। एक पैमाना ( $l$ ) सामग्री की प्राकृतिक मोटाई की तरह है, और दूसरा ( $l_p$ ) एक नया, रहस्यमय पैरामीटर है जो विशेष रूप से असममित क्रिस्टलों के लिए है। सिमुलेशन ने दिखाया कि जैसे-जैसे इन लंबाई के अनुपातों में बदलाव आता है, सामग्री की कठोरता नाटकीय रूप से बढ़ सकती है—कभी-कभी शास्त्रीय भौतिकी की भविष्यवाणी से 40 गुना अधिक कठोर। हालांकि, यह कोई सार्वभौमिक नियम नहीं है; यह वस्तु के आकार और बल की दिशा पर बहुत अधिक निर्भर करता है। उदाहरण के लिए, एक गहरी बीम जिसे शीयर (sideways slide) किया जा रहा था, उसने इन प्रभावों को शायद ही ध्यान में लिया, जबकि एक पतली, मुड़ने वाली बीम अत्यधिक कठोरता के साथ व्यवहार करती है।

लेखक सावधानीपूर्वक नोट करते हैं कि ये कंप्यूटर सिमुलेशन हैं, वास्तविक भौतिक माप नहीं। उन्होंने अभी तक वास्तविक लैब में इन नए "कपल स्ट्रेस" मापदंडों के सटीक मानों को नहीं मापा है क्योंकि ऐसा करने के लिए उपकरण बनाना अविश्वसनीय रूप से कठिन है। इसके बजाय, उन्होंने अपने सिमुलेशन का उपयोग भविष्य के प्रयोगों के लिए एक रोडमैप तैयार करने के लिए किया। वे सुझाव देते हैं कि यदि वैज्ञानिक इन क्रिस्टलों से बनी छोटी कैंटिलीवर बीम (जैसे डाइविंग बोर्ड) बनाते हैं और उन्हें दबाते हैं, तो उन्हें यह पार्श्व झुकाव और भारी कठोरता दिखाई देनी चाहिए। यदि वे ऐसा देखते हैं, तो यह सिद्ध कर देगा कि इन असममित क्रिस्टलों का एक छिपा हुआ, आकार-निर्भर व्यक्तित्व है जिसे शास्त्रीय भौतिकी ने मिस कर दिया था। शोध पत्र का निष्कर्ष है कि हालांकि गणित जटिल है, लेकिन इसके द्वारा की गई भौतिक वास्तविकता सूक्ष्म दुनिया के दबाव में मुड़ने और घूमने की एक नई सीमा को समझने का एक नया मोर्चा है।

तकनीकी सारांश: एक मिश्रित परिवर्तनशील परिमित तत्व विधि के माध्यम से गैर-केंद्रसममित ठोस पदार्थों का आकार-निर्भर इलास्टोस्टैटिक व्यवहार

समस्या विवरण
जबकि शास्त्रीय प्रत्यास्थता (elasticity) और सुसंगत कपल स्ट्रेस थ्योरी (C-CST) को सममित (centrosymmetric) पदार्थों के लिए व्यापक रूप से लागू किया गया है, C-CST ढांचे के भीतर गैर-केंद्रसममित (non-centrosymmetric) ठोस पदार्थों के आकार-निर्भर इलास्टोस्टैटिक प्रतिक्रिया का अन्वेषण काफी हद तक अनछुआ रहा है। गैर-केंद्रसममित क्रिस्टल संरचनाएं, जो पीज़ोइलेक्ट्रिसिटी जैसी घटनाओं के लिए जानी जाती हैं, एक तीसरे क्रम के निर्माण संबंधी टेंसर (constitutive tensor) को धारण करती हैं जो बल-तनाव (force-stress) को विकृति (strain) के साथ और कपल-तनाव (couple-stress) को माध्य वक्रता (mean curvature) के साथ जोड़ती है। यह शोध पत्र उन कम्प्यूटेशनल उपकरणों और सैद्धांतिक सूत्रीकरणों की कमी को संबोधित करता है जो सूक्ष्म और नैनो-पैमानों पर इन पदार्थों के विश्लेषण के लिए आवश्यक हैं, जहाँ आकार-निर्भर प्रभाव महत्वपूर्ण हो जाते हैं। प्राथमिक चुनौती एक ऐसा सूत्रीकरण विकसित करने में निहित है जो विस्थापन क्षेत्र (displacement field) के लिए उच्च-क्रम निरंतरता (high-order continuity) की आवश्यकता के बिना बल-तनाव, कपल-तनाव, विकृति और इंजीनियरिंग माध्य वक्रता के बीच ऊष्मप्रवैगिकी (thermodynamic) और सातत्य यांत्रिकी (continuum mechanical) अंतःक्रियाओं को पूर्णतः समाहित कर सके।

कार्यप्रणाली
लेखक सुसंगत कपल स्ट्रेस थ्योरी (C-CST) पर आधारित एक नवीन मिश्रित परिवर्तनशील परिमित तत्व विधि (mixed variational finite element method - FEM) विकसित करते हैं। कार्यप्रणाली निम्नलिखित चरणों के माध्यम से आगे बढ़ती है:

सैद्धांतिक सूत्रीकरण: प्रयुक्त समीकरणों को C-CST से व्युत्पन्न किया गया है, जिसमें अर्ध-स्थैतिक (quasistatic) स्थितियों और लघु विरूपण (infinitesimal deformations) को माना गया है। निर्माण संबंधी संबंध एक सामान्यीकृत हेकमैन आरेख (Heckmann diagram) दृष्टिकोण का उपयोग करके स्थापित किए गए हैं, जो पीज़ोइलेक्ट्रिसिटी के सादृश्य पर आधारित है। यह एक तीसरे क्रम के निर्माण संबंधी टेंसर ( $Q$ या $P$ ) को पेश करता है जो विकृति और माध्य वक्रता को क्रमशः बल-तनाव और कपल-तनाव से जोड़ता है। यह सूत्रीकरण स्टिफनेस-आधारित और कंप्लायंस-आधारित टेंसरों के बीच अंतर करता है, जो विभिन्न प्रयोगात्मक बाधाओं (जैसे स्थिर विकृति बनाम स्थिर तनाव) के तहत उनके संबंधों को परिभाषित करते हैं।
मिश्रित परिवर्तनशील सिद्धांत (Mixed Variational Principle): मानक विस्थापन-आधारित सूत्रीकरणों द्वारा आवश्यक $C^1$ निरंतरता (जिसमें विस्थापन के दूसरे डेरिवेटिव शामिल होते हैं) की कम्प्यूटेशनल कठिनाई से बचने के लिए, लेखक एक मिश्रित परिवर्तनशील सिद्धांत विकसित करते हैं। यह सिद्धांत विस्थापनों ( $u_i$ ) और कपल-तनावों ( $\mu_m$ ) को प्राथमिक चरों के रूप में मानता है। गतिज बाधाओं (kinematic constraints) को संतुष्ट करने के लिए लैग्रेंज मल्टीप्लायर्स को पेश करके, लेखक एक "स्थिर मिश्रित क्षमता सिद्धांत" (Principle of Stationary Mixed Potential Energy) व्युत्पन्न करते हैं। यह विस्थापन और कपल-तनाव दोनों क्षेत्रों के लिए निरंतरता की आवश्यकता को $C^0$ तक कम कर देता है।
परिमित तत्व कार्यान्वयन (Finite Element Implementation): दो-आयामी प्लेन स्ट्रेन विविक्तीकरण (discretization) को निम्न-क्रम त्रिकोणीय तत्वों का उपयोग करके लागू किया गया है।
- विस्थापन: मानक तीन-नोड रैखिक आकार कार्यों (shape functions) का उपयोग किया गया है।
- कपल-तनाव: कपल-तनाव वेक्टर का प्रतिनिधित्व करने के लिए, लेखक एज-आधारित आकार कार्यों का उपयोग करते हैं। केंद्रसममित मामलों के विपरीत, जहाँ कपल-तनाव वेक्टर विचलन-मुक्त (divergence-free) होता है, गैर-केंद्रसममित पदार्थों को आम तौर पर तत्व किनारों (element edges) के साथ स्पर्शरेखीय (tangential) और अभिलंब (normal) दोनों घटकों की आवश्यकता होती है। इसके परिणामस्वरूप प्रति किनारा दो स्वतंत्रता स्तर (degrees of freedom) प्राप्त होते हैं।
- संख्यात्मक समाधान: कमजोर रूप (weak form) को सममित और विषम-सममित (skew-symmetric) तनाव घटकों वाले रैखिक बीजगणितीय समीकरणों के एक निकाय में विविक्त किया गया है। इस प्रणाली को MATLAB में स्पार्स मैट्रिक्स तकनीकों का उपयोग करके हल किया जाता है।

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

गैर-केंद्रसममित ठोस पदार्थों का प्रथम अनुप्रयोग: यह कार्य C-CST का उपयोग करके गैर-केंद्रसममित ठोस पदार्थों में आकार-निर्भर इलास्टोस्टैटिक व्यवहार का पहला कम्प्यूटेशनल अन्वेषण प्रस्तुत करता है।
नवीन मिश्रित परिवर्तनशील ढांचा: यह कार्य (जो पहले केंद्रसममित पदार्थों के लिए उपयोग किया गया था) को बाय-एनिसोट्रोपिक (गैर-केंद्रसममित) मामले के लिए विस्तारित करता है, जो बल/कपल-तनाव और विकृति/वक्रता के बीच पूर्ण युग्मन (coupling) को सफलतापूर्वक संभालता है।
निर्माण संबंधी मॉडलिंग (Constitutive Modeling): लेखक गैर-केंद्रसममित पदार्थों के लिए ऊष्मप्रवैगिकी संबंधों और निर्माण संबंधी टेंसरों को स्पष्ट रूप से परिभाषित करते हैं, जिससे यह स्पष्ट होता है कि परीक्षण के दौरान स्थिर रखे गए चरों पर सामग्री के पैरामीटर कैसे निर्भर करते हैं।
एलिमेंट विकास: यह सूत्रीकरण एक विशिष्ट त्रिकोणीय एज-एलिमेंट पेश करता है जो कपल-तनाव के स्पर्शरेखीय और अभिलंब दोनों घटकों को पकड़ने में सक्षम है, जो सामान्य गैर-केंद्रसममित सममिति के लिए आवश्यक है।

परिणाम
कार्यप्रणाली को तीन गैर-केंद्रसममित पदार्थों (क्यूबिक सिलिकॉन कार्बाइड (SiC), हेक्सागोनल SiC, और ट्रिगोनल क्वार्ट्ज) से जुड़े चार कम्प्यूटेशनल उदाहरणों के माध्यम से मान्य किया गया। इन समस्याओं में बेंडिंग के तहत एक प्रिजमैटिक कैंटिलीवर, अक्षीय भार के तहत एक प्रिजमैटिक बार, शियर के तहत एक डीप बीम, और संयुक्त लोडिंग के तहत एक कर्व्ड बीम शामिल थे।

अभिसरण (Convergence): मिश्रित FEM समाधानों ने उत्कृष्ट अभिसरण प्रदर्शित किया। सापेक्ष ऊर्जा नॉर्म (energy norm) बेंडिंग समस्याओं के लिए इकाई से थोड़ा कम दर से और अक्षीय लोडिंग के लिए लगभग इकाई के करीब अभिसरित हुआ, जो रैखिक आकार कार्यों के उपयोग के अनुरूप है।
आकार-निर्भर कठोरता (Size-Dependent Stiffness): जैसा कि केंद्रसममित पदार्थों में देखा गया है, जैसे-जैसे गैर-विमिश्रित आंतरिक लंबाई पैमाने ( $l/h$ ) में वृद्धि होती है, संरचनाएं अधिक कठोर हो जाती हैं और अंततः एक संतृप्ति शासन (saturation regime) तक पहुँच जाती हैं। बेंडिंग-प्रधान समस्याओं (जैसे कैंटिलीवर) में महत्वपूर्ण कठोरता वृद्धि (लगभग ~43 गुना) देखी गई, जबकि शियर-प्रधान समस्याओं (डीप बीम) में न्यूनतम वृद्धि (<2 गुना) देखी गई।
गैर-केंद्रसममित प्रभाव:
- कपलिंग पैरामीटर ( $l_p$ ): दूसरा लंबाई पैमाना पैरामीटर ( $l_p$ ), जो गैर-केंद्रसममित कपलिंग टेंसर से जुड़ा है, पाया गया कि कठोरता को और बढ़ाता है, विशेष रूप से जब $l_p$ , $2l$ के करीब पहुँच जाता है।
- विषमदैशिकता (Anisotropy): प्रतिक्रिया क्रिस्टल ओरिएंटेशन और लोडिंग प्लेन पर अत्यधिक निर्भर थी। उदाहरण के लिए, हेक्सागोनल SiC में, $x-y$ तल में अक्षीय भार और अनुप्रस्थ वक्रता के बीच युग्मन महत्वपूर्ण था, लेकिन $x-z$ तल में नगण्य था।
- रैखिक अनुप्रस्थ विक्षेपण (Linear Transverse Deflection): एक महत्वपूर्ण निष्कर्ष यह था कि अक्षीय संपीड़न भार (axial compressive loading) के तहत, गैर-केंद्रसममित पदार्थ बिना किसी क्रिटिकल बकलिंग लोड के रैखिक अनुप्रस्थ विक्षेपण (बेंडिंग) प्रदर्शित करते हैं। यह एक शुद्ध रैखिक यांत्रिक प्रभाव है जो 'कन्वर्स पीज़ोइलेक्ट्रिक प्रभाव' के समान है, जहाँ विक्षेपण का चिह्न कपलिंग पैरामीटर के चिह्न पर निर्भर करता है।
- तनाव विषमता (Stress Asymmetry): बल-तनाव टेंसर का विषम-सममित भाग गैर-शून्य था, जिससे असमान शियर स्ट्रेस घटक ( $\sigma_{xy} \neq \sigma_{yx}$ ) और शास्त्रीय प्रत्यास्थता की तुलना में विशिष्ट तनाव वितरण उत्पन्न हुआ।

महत्व और दावे
यह शोध पत्र गैर-केंद्रसममित प्रत्यास्थ ठोस पदार्थों में आकार-निर्भर व्यवहार का विश्लेषण करने के लिए पहला कम्प्यूटेशनल ढांचा प्रदान करने का दावा करता है। लेखक इस बात पर जोर देते हैं कि अक्षीय लोडिंग और अनुप्रस्थ वक्रता के बीच देखा गया रैखिक युग्मन, मायावी कपल-स्ट्रेस निर्माण संबंधी मापदंडों ( $l_p$ और $P$ टेंसर घटकों) को मापने के लिए भौतिक प्रयोगों के लिए एक संभावित तंत्र प्रदान करता है। अध्ययन यह सुझाव देता है कि उथले कैंटिलीवर बेंडिंग परीक्षण $l$ के आंतरिक लंबाई पैमाने का अनुमान लगाने के लिए प्रभावी हैं, जबकि अक्षीय लोडिंग परीक्षण गैर-केंद्रसममित कपलिंग मापदंडों को अलग कर सकते हैं। वे निष्कर्ष निकालते हैं कि जबकि वर्तमान कार्य रैखिक इलास्टोस्टैटिक्स पर केंद्रित है, विकसित ढांचा भविष्य के अनुसंधान (कंपन, क्षणिक गतिकी, त्रि-आयामी समस्याएं और बड़े विरूपण सहित) के लिए आधार तैयार करता है, जिसमें केंद्रसममित क्रिस्टल में संभावित बड़े-विकृति प्रेरित सममिति भंग (symmetry breaking) भी शामिल है।

इसी तरह के और पेपर