⚡ electrical engineering

A Physically Constrained Energy-Based Excess Pore Pressure Model for Saturated Sands with an Ordinary Differential Equation Interpretation

यह अध्ययन संतृप्त स्वच्छ रेत (saturated clean sands) के लिए ऊर्जा-आधारित अतिरिक्त छिद्र दबाव मॉडल को एक कॉची-यूलर साधारण अवकल समीकरण (Cauchy–Euler ordinary differential equation) के समाधान के रूप में व्याख्यायित करके पुनर्गठित करता है, जिससे भौतिक सीमाबद्धता सुनिश्चित होती है, अंशांकन को एकल पैरामीटर तक कम किया जाता है, और विभिन्न सापेक्ष घनत्वों में उच्च सटीकता प्राप्त की जाती है।

मूल लेखक: Ferruh Gündoğan, Ayfer Erken

प्रकाशित 2026-06-30

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

CC BY 4.0

मूल लेखक: Ferruh Gündoğan, Ayfer Erken

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

एक बड़ी तस्वीर: रेत के महल क्यों ढह जाते हैं?

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

अंततः, पानी का दबाव इतना बढ़ जाता है कि वह रेत के कणों को पूरी तरह से अलग कर देता है। रेत एक ठोस की तरह व्यवहार करना बंद कर देती है और एक गाढ़े तरल (जैसे मिल्कशेक) की तरह व्यवहार करने लगती है। इसे द्रवीकरण (Liquefaction) कहा जाता है, और यही भूकंप के दौरान इमारतों के धंसने या झुकने का कारण बनता है।

इंजीनियरों को यह अनुमान लगाने के लिए एक तरीके की आवश्यकता होती है कि यह ठीक कब होता है। वे अपव्यस्त ऊर्जा (dissipated energy) नामक एक अवधारणा का उपयोग करते हैं—इसे "हिलने-डुलने" या काम की कुल मात्रा के रूप में समझें जो रेत को खुद को पुनर्गठित करने के लिए करनी पड़ती है। आप इसे जितना अधिक हिलाएंगे, यह उतनी ही अधिक ऊर्जा अवशोषित करेगा, और पानी का दबाव उतना ही अधिक बढ़ेगा।

पुरानी समस्या: एक टूटा हुआ पैमाना

द दशकों से, इंजीनियर इस पानी के दबाव का अनुमान लगाने के लिए एक प्रसिद्ध सूत्र (जो बेरिल और डेविस द्वारा 1985 में बनाया गया था) का उपयोग करते आए हैं। यह एक साधारण पैमाने की तरह है जो कहता है: "आप जितनी अधिक ऊर्जा डालेंगे, दबाव उतना ही अधिक बढ़ेगा।"

हालाँकि, इस पुराने पैमाने में दो प्रमुख कमियां हैं:

इसकी कोई सीमा नहीं है: यह सूत्र अनंत तक बढ़ता रहता है। वास्तविक दुनिया में, पानी का दबाव मिट्टी के ऊपर के वजन से अधिक नहीं हो सकता (यह 100% से ऊपर नहीं जा सकता)। पुराना सूत्र कभी-कभी असंभव संख्याएं बताता है, जैसे 150% दबाव, जो भौतिक विज्ञान के नियमों के विरुद्ध है।
इसके 'नॉब्स' उलझे हुए हैं: इस सूत्र में दो सेटिंग्स (पैरामीटर) हैं जिन्हें स्वतंत्र होना चाहिए, लेकिन वे वास्तव में आपस में जुड़े हुए हैं। यदि आप एक को बदलते हैं, तो दूसरा भी बदल जाता है। यह विभिन्न प्रकार की रेत के लिए मॉडल को सटीक बनाना बहुत कठिन बना देता है।

नया समाधान: एक स्मार्ट, सीमित मॉडल

इस शोध पत्र के लेखकों (गुंडोगान और एर्केन) ने इस पैमाने को ठीक करने का निर्णय लिया। उन्होंने दो प्रकार की रेत पर वास्तविक लैब परीक्षणों के डेटा का अध्ययन किया: ओटावा F65 (एक मानक, साफ रेत जिसका उपयोग कई अध्ययनों में किया जाता है) और अकपिनार सैंड (इस्तांबुल की एक स्थानीय रेत)।

केवल एक ग्राफ पर सीधी रेखा खींचने के बजाय, उन्होंने महसूस किया कि हिलने वाली ऊर्जा और पानी के दबाव के बीच का संबंध एक तीन-चरणीय यात्रा की तरह है:

शुरुआत: एक धीमी, हल्की बढ़त।
मध्य: एक तीव्र, खड़ी चढ़ाई।
अंत: एक चिकना वक्र (curve) जो ठीक सीमा (100% दबाव) पर आकर सपाट हो जाता है।

इसे पकड़ने के लिए, उन्होंने गणित को एक नुस्खे (recipe) की तरह माना। उन्होंने सूत्र को दो भागों में विभाजित किया:

ग्रोथ इंजन (Growth Engine): एक बुनियादी शब्द जो तीव्र उछाल को संभालता है।
मध्यस्थ (Mediator): एक "ट्रैफिक कंट्रोलर" जो सीमा के करीब पहुँचते ही विकास को धीमा कर देता है, जिससे यह सुनिश्चित होता है कि यह कभी भी 100% की छत को पार न करे।

"सिंथेटिक" रहस्य: एक गणितीय कहानी

यहाँ उनकी विधि का चतुर हिस्सा है। उन्होंने देखा कि उनके नए वक्र (curve) का आकार बिल्कुल एक विशिष्ट प्रकार के भौतिक समीकरण के समाधान जैसा दिखता है जिसे साधारण अवकल समीकरण (Ordinary Differential Equation - ODE) कहा जाता है।

एक ODE को एक ऐसी कहानी के रूप में सोचें कि कैसे एक सिस्टम समय के साथ बदलता है।

समीकरण का बायां हिस्सा रेत के आंतरिक व्यवहार (कैसे वह प्रतिरोध करती है) का प्रतिनिधित्व करता है।
दायां हिस्सा बाहरी कंपन (भूकंप) का प्रतिनिधित्व करता है।

इस समीकरण के रूप में अपने मॉडल को फिर से लिखकर, वे सख्त "सड़क के नियम" (भौतिक बाधाएं) लागू कर सके:

नियम 1: दबाव हमेशा बढ़ना चाहिए (आप रेत को "अन-शेक" नहीं कर सकते)।
नियम 2: दबाव ठीक तब रुक जाना चाहिए जब यह 100% (द्रवीकरण) पर पहुँचता है।
नियम 3: जिस क्षण यह 100% पर पहुँचता है, बढ़ने की गति शून्य होनी चाहिए (यह सहजता से समतल हो जाता है, न कि किसी नुकीले उछाल के साथ)।

गणित को इन नियमों का पालन करने के लिए मजबूर करके, वे "उलझे हुए नॉब्स" को खत्म करने में सक्षम हुए। उन्होंने पूरे जटिल मॉडल को केवल एक संख्या (जिसे b कहा जाता है) तक सीमित कर दिया जिसे उन्हें खोजना है। बाकी गणित भौतिकी के नियमों द्वारा निर्धारित है।

उन्होंने क्या पाया

उन्होंने अपने नए "वन-नंबर मॉडल" का लैब डेटा के साथ परीक्षण किया:

सटीकता: यह वास्तविक लैब परिणामों के साथ अविश्वसनीय रूप से मेल खाता है (99% सटीकता)।
सुरक्षा: यह कभी भी असंभव संख्या का अनुमान नहीं लगाता (यह 0% और 100% के बीच रहता है)।
सरलता: यह पुराने मॉडल की तुलना में उपयोग करने में बहुत आसान है क्योंकि आपको केवल एक सेटिंग को ट्यून करने की आवश्यकता होती है।

"अहा!" क्षण: विफलता के दो चरण

गणित ने एक दिलचस्प बात का खुलासा किया कि रेत कैसे विफल होती है।

प्रारंभिक चरण: रेत की अपनी आंतरिक संरचना (कण कितने कसकर भरे हुए हैं) दबाव वृद्धि को नियंत्रित करती है। यह समीकरण का "होमोजेनियस" (Homogeneous) हिस्सा है।
अंतिम चरण: एक बार जब रेत द्रवीकृत होने के बहुत करीब पहुँच जाती है (लग लगभग 75% दबाव), तो वह अपनी "याददाश्त" खो देती है। उसकी आंतरिक संरचना टूट जाती है, और दबाव का बढ़ना लगभग पूरी तरह से बाहरी कंपन ऊर्जा द्वारा संचालित होता है। यह समीकरण का "पार्टिक्युलर" (Particular) हिस्सा है।

सारांश

लेखकों ने रेत के द्रवीकरण की भविष्यवाणी करने वाले एक अस्त-व्यस्त, असीमित और जटिल सूत्र को फिर से बनाया। उन्होंने एक गणितीय "कहानी" (एक समीकरण) का उपयोग किया ताकि मॉडल को वास्तविक भौतिकी की तरह व्यवहार करने के लिए मजबूर किया जा सके: सुचारू रूप से ऊपर उठना और ठीक सीमा पर रुकना। परिणाम एक सरल, सुरक्षित और अधिक सटीक उपकरण है जिससे इंजीनियर यह समझ सकते हैं कि भूकंप के दौरान रेत कब तरल में बदल जाएगी।

तकनीकी सारांश: संतृप्त रेत के लिए एक भौतिक रूप से बाधित ऊर्जा-आधारित अतिरिक्त पोर प्रेशर मॉडल

समस्या विवरण
मिट्टी का द्रवीकरण (Liquefaction) एक महत्वपूर्ण भू-तकनीकी खतरा बना हुआ है, जिसका पारंपरिक रूप से स्ट्रेस-आधारित, स्ट्रेन-आधारित, पोर-प्रेशर-आधारित और ऊर्जा-आधारित विधियों के माध्यम से मूल्यांकन किया जाता है। इनमें से, ऊर्जा-आधारित विधियाँ (EBM) एक विशिष्ट लाभ प्रदान करती हैं क्योंकि वे संचयी विसरित ऊर्जा (cumulative dissipated energy) के साथ अतिरिक्त पोर प्रेशर जनरेशन को जोड़ती हैं, जिससे तुल्य-चक्र रूपांतरणों (equivalent-cycle conversions) की अनिश्चितता के बिना अनियमित भूकंपीय लोडिंग को समायोजित किया जा सके। सबसे प्रभावशाली EBM सूत्रीकरण बेरिल-डेविस (Berrill–Davis, 1985) पावर लॉ ( $r_u = \alpha W^\beta$ ) है। हालांकि, यह अध्ययन बेरिल-डेविस मॉडल में दो मौलिक सीमाओं की पहचान करता है:

अबाधित प्रतिक्रिया (Unbounded Response): पावर लॉ गणितीय रूप से अबाधित है, जो अतिरिक्त पोर प्रेशर अनुपात ( $r_u$ ) को इकाई (unity) से अधिक होने की अनुमति देता है, जो संतृप्त रेत के लिए भौतिक रूप से असंभव है जहाँ $r_u$ द्रवीकरण पर 1.0 के रूप में स्पर्शोन्मुखी (asymptotically) होना चाहिए।
पैरामीटर अन्योन्याश्रयता (Parameter Interdependence): अंशांकन (calibration) पैरामीटर $\alpha$ और $\beta$ को अक्सर स्वतंत्र माना जाता है, फिर भी अनुभवजन्य अध्ययन दर्शाते हैं कि वे परस्पर सह-संबंधित हैं। यह निर्भरता अंशांकन को जटिल बनाती है, भौतिक व्याख्यात्मकता को कम करती है, और एक्सट्रपलेशन (extrapolation) में बाधा डालती है।

कार्यप्रणाली
लेखक एक पुनर्गठित मॉडल प्रस्तावित करते हैं जो समान, स्वच्छ, संतृप्त रेत के लिए है और जो इन सीमाओं को एक बहु-चरणीय व्युत्पन्न प्रक्रिया के माध्यम से संबोधित करता है:

डेटा विश्लेषण और अपघटन: ओटावा F65 रेत (LEAP-2017) और अक्पिंर (Akpınar) रेत के चक्रीय त्रिकक्षीय (cyclic triaxial) डेटा का उपयोग करते हुए, सामान्यीकृत पोर प्रेशर अनुपात ( $r_u$ ) और सामान्यीकृत विसरित ऊर्जा ( $W_N$ ) के बीच संबंध का लॉग-लॉग स्केल पर विश्लेषण किया गया है। लेखक देखते हैं कि डेटा एक एकल सीधी रेखा का पालन नहीं करता है बल्कि कम से कम तीन रैखिक खंडों से बना है।
बहुपद सन्निकटन (Polynomial Approximation): डेटा को बिना किसी स्थिरांक पद के छठे क्रम के बहुपद द्वारा अनुमानित किया गया है। इस फलन को एक प्राथमिक विकास पद ( $W_N^b$ ) और एक मध्यस्थ फलन (mediation function) $g(W_N)$ में विघटित किया गया है।
टेलर विस्तार और सरलीकरण: मध्यस्थ फलन का आगे विश्लेषण किया गया है। एक्सपोनेंशियल घटक का टेलर विस्तार यह प्रकट करता है कि एक तृतीय-क्रम का बहुपद डेटा प्रवृत्तियों का पर्याप्त प्रतिनिधित्व करता है। इससे एक संभावित मॉडल रूप प्राप्त होता है: $r_u = W_N^b (a_0 + a_1 W_N + a_2 W_N^2 + a_3 W_N^3)$ ।
साधारण अवकल समीकरण (ODE) व्याख्या: गणितीय संरचना को पहचानते हुए, मॉडल को एक गैर-समांगी कॉची-यूलर (non-homogeneous Cauchy–Euler) साधारण अवकल समीकरण के समाधान के रूप में पुनर्व्याख्यायित किया गया है। बायां पक्ष प्रणाली की आंतरिक प्रतिक्रिया का प्रतिनिधित्व करता है, जबकि दायां पक्ष बाहरी लोडिंग स्थितियों का प्रतिनिधित्व करता है।
भौतिक बाधाएं: अनभौतिक उतार-चढ़ाव और पैरामीटर अन्योन्याश्रयता को समाप्त करने के लिए, मॉडल को द्रवीकरण बिंदु ( $W_N = 1, r_u = 1$ $W_{N} = 1, r_{u} = 1$ ) पर विशिष्ट सीमा स्थितियों के अधीन किया गया है:
- $W_N = 1$ पर $r_u = 1$ ।
- $W_N = 1$ पर एक क्षैतिज स्पर्शोन्मुख (horizontal tangent, $r_u' = 0}$ ), जो यह सुनिश्चित करता है कि वक्र सीमा तक पहुँचते समय समतल हो जाए।
- लोडिंग प्रक्रिया के दौरान निरंतर वृद्धि (monotonic increase, $r_u' > 0$ )।
- गुणांकों के बीच एक विशिष्ट संबंध ( $a_2 = -3a_3$ ), जो द्रवीकरण पर गैर-समांगी पद के शून्य होने से व्युत्पन्न है।

ये बाधाएं मध्यस्थ बहुपद के गुणांकों को $b$ के एकल पैरामीटर के संदर्भ में विश्लेषणात्मक रूप से निर्धारित करती हैं, जिससे मॉडल केवल $b$ द्वारा नियंत्रित एक रूप में सिमट जाता है।

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

भौतिक रूप से बाधित सूत्रीकरण: प्रस्तावित मॉडल यह सुनिश्चित करने के लिए निर्मित किया गया है कि $0 \le r_u \le 1$ बना रहे, जो बेरिल-डेविस मॉडल में सामान्य अनभौतिक पोर प्रेशर ओवर-प्रेडिक्शन को रोकता है।
एकल अंशांकन पैरामीटर: मध्यस्थ गुणांकों को विकास घातांक $b$ से विश्लेषणात्मक रूप से जोड़कर, मॉडल कई स्वतंत्र पैरामीटर के अंशांकन की आवश्यकता को समाप्त करता है, जो पिछले सूत्रीकरणों में पाई जाने वाली अन्योन्याश्रयता की समस्या को हल करता है।
ODE ढांचा: यह अध्ययन पोर प्रेशर संबंध की एक गैर-समांगी कॉची-यूलर ODE के समाधान के रूप में एक नवीन व्याख्या प्रदान करता है। यह ढांचे को आंतरिक सामग्री विशेषताओं (होमोजेनियस सॉल्यूशन, जो प्रारंभिक चरणों में प्रभावी है) और बाहरी लोडिंग स्थितियों (पार्टिकुलर सॉल्यूशन) में विभाजित करता है।
सत्यापन: मॉडल को ओटावा F65 और अक्पिंर रेत के उच्च-गुणवत्ता वाले चक्रीय त्रिकक्षीय डेटा के विरुद्ध विभिन्न सापेक्ष घनत्वों (ढीली से बहुत घनी रेत तक) के लिए मान्य किया गया है, जिसमें $R^2 \ge 0.99$ प्राप्त हुआ है।

परिणाम

सरल मॉडल (समीकरण 20) ओटावा F65 और अक्पिंर रेत दोनों के मापे गए पोर-प्रेशर इतिहास को सफलतापूर्वक पुनरुत्पादित करता है।
पाया गया कि पैरामीटर $b$ सापेक्ष घनत्व के साथ बदलता है (जैसे, 51% सापेक्ष घनत्व के लिए $b \approx 0.36$ , 80–91% घनत्व के लिए $b \approx 0.43$ )।
मॉडल पोर प्रेशर के गैर-रेखीय विकास को कैप्चर करता है, जिसमें प्रारंभिक तीव्र वृद्धि और पूर्ण द्रवीकरण के करीब पहुँचने पर बाद का पठार (plateauing) शामिल है।
मूल बेरिल-डेविस मॉडल के साथ तुलना दर्शाती है कि प्रस्तावित सूत्रीकरण पोर प्रेशर विकास का अधिक सुसंगत विवरण प्रदान करता है, विशेष रूप से पूर्ण द्रवीकरण की अंतिम स्थिति के पास, जहाँ मूल मॉडल महत्वपूर्ण बिखराव और अबाधित वृद्धि प्रदर्शित करता है।
शेष LEAP-2017 डेटासेट पर सरल मॉडल का अनुप्रयोग विभिन्न सापेक्ष घनत्वों और चक्रीय स्ट्रेस रेश्यो (CSR) के लिए पोर प्रेशर जनरेशन को सिम्युलेट करने की इसकी क्षमता की पुष्टि करता है, बशर्ते कि स्थानीय अधिकतम (local maxima) से बचने के लिए $b \ge 0.27$ हो।

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