Hybrid Dynamics of Rocking Blocks Beyond Overturning: Saltation Analysis, Bifurcations, and Stability Characterization

तकनीकी सारांश: शास्त्रीय और वैकल्पिक प्रतिस्थापन मॉडलों के साथ रॉकिंग ब्लॉक व्यवहारों का अभिलक्षण

समस्या विवरण
यह अध्ययन रॉकिंग ब्लॉकों में प्रयोगात्मक रूप से देखे गए प्रभाव व्यवहारों और शास्त्रीय हाउसनर (Housner) प्रत्यास्थता गुणांक (restitution coefficient) की भविष्यवाणियों के बीच विसंगति को संबोधित करता है। जबकि हाउसनर मॉडल, जो नए धुरी (pivot) के चारों ओर कोणीय संवेग के संरक्षण को मानता है, कठोर-पिंड रॉकिंग गतिकी के विश्लेषण के लिए एक आधार रहा है, यह अक्सर प्रभाव प्रतिक्रियाओं को सटीक रूप से पकड़ने में विफल रहता है, विशेष रूप से कम स्लींडनेस अनुपात (slenderness ratios) वाले ब्लॉकों के लिए। यद्यपि माओ एट अल. (Mao et al.) द्वारा प्रस्तावित जैसे वैकल्पिक प्रत्यास्थता सूत्र, प्रयोगात्मक डेटा के साथ बेहतर तालमेल दिखाते हैं, लेकिन रॉकिंग प्रणालियों की वैश्विक गैर-रैखिक गतिकी—विशेष रूप से जटिल दोलनों, अट्रैक्टर स्थिरता और पलटने की संवेदनशीलता (overturning susceptibility) के संबंध में—पर उनके निहितार्थ अभी भी काफी हद तक अनछुए हैं।

कार्यप्रणाली
लेखक एक रॉकिंग ब्लॉक को एक हाइब्रिड नॉन-स्मूथ डायनेमिकल सिस्टम के रूप में तैयार करते हैं। भौतिक मॉडल द्रव्यमान $M$, चौड़ाई $2b$, और ऊँचाई $2h$ वाला एक समांग आयताकार ब्लॉक है, जो हार्मोनिक बेस एक्सिटेशन (harmonic base excitation) के तहत एक कठोर आधार पर बिना फिसले रॉक कर रहा है। गति के समीकरण सक्रिय धुरी के चारों ओर मोमेंट के संतुलन से व्युत्पन्न किए गए हैं, जिसमें आयाम ($\alpha$), आवृत्ति ($\omega$), और श्यान अवमंदन (viscous damping) ($\xi$) के लिए आयामहीन पैरामीटर शामिल हैं।

प्रणाली का विश्लेषण करने के लिए, लेखक दो प्रत्यास्थता मॉडलों की तुलना करते हैं:

1. शास्त्रीय हाउसनर मॉडल ($e_H$): कोणीय संवेग संरक्षण से व्युत्पन्न, जो केवल ज्यामिति ($e_H = 1 - \frac{3}{2}\cos^2(\theta_c)$) पर निर्भर है।
2. वैकल्पिक माओ एट अल. मॉडल ($e_M$): गतिज ऊर्जा पुनर्वितरण और वेक्टर-आधारित वेग क्षेत्र विश्लेषण से व्युत्पन्न, जो आम तौर पर कम ऊर्जा अपव्यय और उच्च प्रत्यास्थता मानों की भविष्यवाणी करता है, विशेष रूप से कम-स्लींडनेस वाले ब्लॉकों के लिए।

गतिकीय अभिलक्षणन (Dynamical characterization) के लिए तीन प्राथमिक उपकरणों का उपयोग किया गया है:

• बाइफरकेशन आरेख (Bifurcation Diagrams): बल आयाम $\alpha$ के परिवर्तन के साथ गुणात्मक परिवर्तनों (आवधिक, बहु-आवधिक, अराजक) की पहचान करने के लिए।
• अधिकतम लियापुनोव घातांक ($\lambda_{max}$): प्रभाव विच्छिन्नता (impact discontinuities) को संभालने के लिए साल्टेशन मैट्रिसेस (saltation matrices) के साथ बेनेटिन-वुल्फ एल्गोरिदम का उपयोग करके गणना की गई, जो स्थिर आवधिक प्रतिक्रियाओं ($\lambda_{max} < 0$) और अराजक गतिकी ($\lambda_{max} > 0$) के बीच अंतर करती है।
• अट्रैक्शन के बेसिन (Basins of Attraction): दीर्घकालिक गतिकी को प्रारंभिक स्थितियों $(\theta_0, \theta'_0)$ पर निर्भरता को मैप करने के लिए, जिसमें परिणाम (डैम्प्ड मोशन, निरंतर दोलन [सममित, असममित, बहु-आवधिक], और पलटना [$|\theta| \ge \pi/2$]) को वर्गीकृत किया गया है।

सिमुलेशन पायथन 3 के साथ इम्प्लिसिट BDF इंटीग्रेशन का उपयोग करके दो स्लींडनेस अनुपात ($h/b = 1$ और $h/b = 2$) के लिए आयोजित किए गए थे।

प्रमुख योगदान और परिणाम
अध्ययन यह प्रदर्शित करता है कि प्रत्यास्थता मॉडल का चयन अनुमानित प्रणाली प्रतिक्रिया को महत्वपूर्ण रूप से प्रभावित करता है, विशेष रूप से कम-स्लींडनेस ब्लॉकों के लिए।

1. गतिकीय व्यवस्थाओं में विचलन: $h/b = 1$ के लिए, माओ एट अल. मॉडल हाउसनर मॉडल की तुलना में जटिल दोलन व्यवहारों की जल्दी शुरुआत और अधिक व्यापकता की भविष्यवाणी करता है। वैकल्पिक सूत्र अनियमित व्यवहार की ओर संक्रमण को कम बल आयामों की ओर स्थानांतरित कर देता है, जिसमें बिखरे हुए बाइफरकेशन संरचनाओं और मुख्य रूप से सकारात्मक लियापुनोव घातांकों के साथ व्यापक पैरामीटर अंतराल प्रदर्शित होते हैं।
2. अट्रैक्शन के बेसिन की संरचना: प्रत्यास्थता मॉडल बंधित गतियों की वैश्विक सुलभता को बदल देता है। $h/b = 1$ के लिए, हाउसनर मॉडल सममित दोलनों के लिए एक विस्तृत, सघन बेसिन प्रदान करता है, जो प्रारंभिक स्थिति गड़बड़ी के विरुद्ध मजबूती को दर्शाता है। इसके विपरीत, माओ एट अल. मॉडल एक काफी संकीर्ण दोलनी बेसिन उत्पन्न करता है जिसकी अनियमित सीमाएं एक प्रमुख पलटने वाले क्षेत्र (overturning region) के भीतर अंतर्निहित होती हैं। यह विखंडन प्रारंभिक स्थितियों के प्रति बढ़ती संवेदनशीलता को दर्शाता है।
3. मात्रात्मक असहमति: $h/b = 1$ के लिए, 10,000 प्रारंभिक स्थितियों के संबंध में दोनों मॉडलों के बीच अंतिम गतिक परिणाम के बारे में सहमति केवल 64.09% है। माओ एट अल. का सूत्र पलटने वाली घटनाओं के अनुपात को काफी बढ़ा देता है (80.36% बनाम हाउसनर के लिए 44.45%) और नियमित सममित दोलनों के डोमेन को कम कर देता है।
4. स्लींडनेस के साथ अभिसरण: जैसे-जैसे स्लींडनेस अनुपात बढ़कर $h/b = 2$ हो जाता है, दोनों फॉर्मुलेशन द्वारा उत्पन्न गतिक विशेषताएं प्रगतिशील रूप से अभिसरित होती हैं। गतिक परिणामों में सहमति बढ़कर 86.93% हो जाती है, और बाइफरकेशन संरचनाओं और लियापुनोव घातांक वितरणों में अंतर कम स्पष्ट हो जाता है। यह $\theta_c$ के बढ़ने के साथ $e_H$ और $e_M$ के अभिसरण के अनुरूप है।

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