2025-11-14T17:43:11.826412

Macroscopic stress, couple stress and flux tensors derived through energetic equivalence from microscopic continuous and discrete heterogeneous finite representative volumes

Eliáš, Cusatis
This paper presents a rigorous derivation of equations to evaluate the macroscopic stress tensor, the couple stress tensor, and the flux vector equivalent to underlying microscopic fields in continuous and discrete heterogeneous systems with independent displacements and rotations. Contrary to the classical asymptotic expansion homogenization, finite size representative volume is considered. First, the macroscopic quantities are derived for a heterogeneous Cosserat continuum. The resulting continuum equations are discretized to provide macroscopic quantities in discrete heterogeneous systems. Finally, the expressions for discrete system are derived once again, this time considering the discrete nature directly. The formulations are presented in two variants, considering either internal or external forces, couples, and fluxes. The derivation is based on the virtual work equivalence and elucidates the fundamental significance of the couple stress tensor in the context of balance equations and admissible virtual deformation modes. Notably, an additional term in the couple stress tensor formula emerges, explaining its dependence on the reference system and position of the macroscopic point. The resulting equations are verified by comparing their predictions with known analytical solutions and results of other numerical models under both steady state and transient conditions.
academic

सूक्ष्म निरंतर और असतत विषमांगी परिमित प्रतिनिधि आयतन से ऊर्जा समतुल्यता के माध्यम से व्युत्पन्न स्थूल प्रतिबल, युग्म प्रतिबल और प्रवाह प्रदिश

मूल जानकारी

  • पत्र ID: 2412.12616
  • शीर्षक: Macroscopic stress, couple stress and flux tensors derived through energetic equivalence from microscopic continuous and discrete heterogeneous finite representative volumes
  • लेखक: Jan Eliáš (ब्रनो प्रौद्योगिकी विश्वविद्यालय), Gianluca Cusatis (नॉर्थवेस्टर्न विश्वविद्यालय)
  • वर्गीकरण: cs.CE (कम्प्यूटेशनल इंजीनियरिंग), cond-mat.other
  • प्रस्तुति समय: दिसंबर 2024
  • पत्र लिंक: https://arxiv.org/abs/2412.12616

सारांश

यह पत्र निरंतर और असतत विषमांगी प्रणालियों में स्थूल प्रतिबल प्रदिश, युग्म प्रतिबल प्रदिश और प्रवाह सदिश का मूल्यांकन करने के लिए एक कठोर व्युत्पत्ति विधि प्रस्तावित करता है, जिनमें स्वतंत्र विस्थापन और घूर्णन हैं। शास्त्रीय अनंतक विस्तार समांगीकरण के विपरीत, यह पत्र परिमित आकार के प्रतिनिधि आयतन तत्व (RVE) पर विचार करता है। पहले विषमांगी Cosserat निरंतर पदार्थ के लिए स्थूल राशियों को व्युत्पन्न किया जाता है, फिर निरंतर पदार्थ समीकरणों को असतत किया जाता है ताकि असतत विषमांगी प्रणालियों में स्थूल राशियां प्रदान की जा सकें, और अंत में असतत गुणों से सीधे असतत प्रणाली के लिए व्यंजकों को पुनः व्युत्पन्न किया जाता है।

अनुसंधान पृष्ठभूमि और प्रेरणा

1. समाधान की जाने वाली समस्याएं

यह पत्र मुख्य रूप से सूक्ष्म पैमाने की विषमांगी सामग्री प्रणालियों से स्थूल यांत्रिकी और परिवहन विशेषताओं को सटीक रूप से निकालने की समस्या को हल करता है, विशेष रूप से:

  • परिमित आकार के प्रतिनिधि आयतन तत्व से स्थूल प्रतिबल प्रदिश को कैसे व्युत्पन्न किया जाए
  • स्वतंत्र घूर्णन स्वतंत्रता वाली प्रणालियों में युग्म प्रतिबल प्रदिश को कैसे संभाला जाए
  • सूक्ष्म और स्थूल पैमानों के बीच ऊर्जा समतुल्यता संबंध को कैसे स्थापित किया जाए

2. समस्या की महत्ता

कई इंजीनियरिंग सामग्रियां (जैसे कंक्रीट, सिरेमिक, दानेदार सामग्री आदि) सूक्ष्म या सूक्ष्मदर्शी पैमाने पर महत्वपूर्ण विषमांगता प्रदर्शित करती हैं। सटीक बहु-पैमाने मॉडलिंग के लिए महत्वपूर्ण है:

  • सामग्री के समग्र यांत्रिक व्यवहार की भविष्यवाणी
  • विश्वसनीय संख्यात्मक मॉडल का विकास
  • घटनात्मक संवैधानिक संबंध स्थापित करना

3. मौजूदा विधियों की सीमाएं

  • अनंतक विस्तार समांगीकरण: प्रतिनिधि आयतन तत्व को शून्य की ओर जाने के लिए मानता है, परिमित आकार के प्रभावों को नजरअंदाज करता है
  • शास्त्रीय Love-Weber सूत्र: मुख्य रूप से प्रतिबल प्रदिश पर लागू होता है, युग्म प्रतिबल प्रदिश के संभालने में अधूरा है
  • सीमा शर्त संभालना: मौजूदा विधियों में सीमा प्रभाव और संदर्भ बिंदु चयन को संभालने में कमी है

4. अनुसंधान प्रेरणा

यह पत्र एक अधिक संपूर्ण और कठोर सैद्धांतिक ढांचा स्थापित करने का लक्ष्य रखता है जो सक्षम हो:

  • परिमित आकार के प्रतिनिधि आयतन तत्वों को संभालना
  • आंतरिक बल और बाहरी बल दोनों सूत्रों पर एक साथ विचार करना
  • युग्म प्रतिबल प्रदिश में संदर्भ बिंदु स्थिति के प्रभाव को स्पष्ट करना

मूल योगदान

  1. आभासी कार्य समतुल्यता पर आधारित कठोर व्युत्पत्ति ढांचा स्थापित किया, जो निरंतर और असतत विषमांगी प्रणालियों पर लागू होता है
  2. परिमित RVE आकार पर विचार करने वाली स्थूल राशियों के व्यंजकों को व्युत्पन्न किया, जिसमें प्रतिबल प्रदिश, युग्म प्रतिबल प्रदिश और प्रवाह सदिश शामिल हैं
  3. दो प्रकार के रूप प्रस्तावित किए: आंतरिक बल और बाहरी बल आधारित सूत्रण
  4. युग्म प्रतिबल प्रदिश में अतिरिक्त पद खोजे, जो संदर्भ तंत्र और स्थूल बिंदु स्थिति पर इसकी निर्भरता को समझाता है
  5. संपूर्ण सत्यापन प्रदान किया, जिसमें स्थिर और क्षणिक स्थितियों के तहत विश्लेषणात्मक समाधान और संख्यात्मक मॉडल तुलना शामिल है

विधि विवरण

कार्य परिभाषा

इस पत्र का मूल कार्य सूक्ष्म विषमांगी सामग्री और स्थूल समांगी निरंतर पदार्थ के बीच समतुल्यता संबंध स्थापित करना है, विशेष रूप से:

  • इनपुट: सूक्ष्म पैमाने का प्रतिबल क्षेत्र σ^mic, युग्म प्रतिबल क्षेत्र μ^mic, प्रवाह क्षेत्र a^mic
  • आउटपुट: स्थूल प्रतिबल प्रदिश σ^mac, युग्म प्रतिबल प्रदिश μ^mac, प्रवाह सदिश a^mac
  • बाधाएं: आभासी कार्य समतुल्यता सिद्धांत और संतुलन समीकरणों को संतुष्ट करना

सैद्धांतिक ढांचा

1. Cosserat निरंतर पदार्थ की नींव

Cosserat निरंतर पदार्थ के मुख्य चर शामिल हैं:

  • विस्थापन सदिश u और स्वतंत्र घूर्णन सदिश θ
  • विकृति प्रदिश γ = ε - E·θ^m और वक्रता प्रदिश κ = ∇⊗θ
  • प्रतिबल प्रदिश σ (असममित) और युग्म प्रतिबल प्रदिश μ

संतुलन समीकरण हैं:

-ρü + ∇·σ + f = 0                    (रैखिक गति)
-ρJ_ρθ̈ + ∇·μ + E:σ + z = 0         (कोणीय गति)

2. आभासी कार्य समतुल्यता सिद्धांत

मुख्य धारणा यह है कि आभासी विस्थापन और आभासी घूर्णन को स्थान निर्देशांक की श्रृंखला के रूप में व्यक्त किया जा सकता है:

δu_j(x) = α_j^(0) + x_iα_ij^(1) + x_kx_iα_ikj^(2) + ...
δθ_j(x) = β_j^(0) + x_iβ_ij^(1) + ...

आभासी कार्य समतुल्यता शर्त के माध्यम से:

δW_int^mac = δW̄_int^mic = δW̄_ext^mic

3. स्थूल प्रतिबल प्रदिश व्युत्पत्ति

खिंचाव और कतरनी विकृति मोड के लिए (α^(1) पद), प्राप्त होता है:

σ_ij^mac = ⟨σ_ij^mic⟩ = (1/V)∫_Ω σ_ij^mic dV        (आंतरिक बल रूप)
σ_ij^mac = (1/V)[∫_Γ x_it_j^mic dΓ + ∫_Ω x_if_j^mic dV]  (बाहरी बल रूप)

4. स्थूल युग्म प्रतिबल प्रदिश व्युत्पत्ति

मोड़ और ट्विस्ट विकृति मोड के लिए (β^(1) पद), प्राप्त होता है:

μ_ij^mac = ⟨μ_ij^mic + E_jklx_kσ_il^mic⟩ - E_jklx_k^macσ_il^mac

मुख्य खोज यह है कि युग्म प्रतिबल प्रदिश में एक अतिरिक्त पद होता है जो स्थूल संदर्भ बिंदु स्थिति पर निर्भर करता है।

असतत प्रणाली का संभालना

1. बाहरी बल सूत्रण का असतत करण

निरंतर पदार्थ व्यंजकों को सीधे असतत करना:

σ^mac = (1/V)∑_f x_f ⊗ F_f                    (प्रतिबल प्रदिश)
μ^mac = σ^mac ⊗ x^mac : E + (1/V)∑_z x_z ⊗ Z_z  (युग्म प्रतिबल प्रदिश)
a^mac = -(1/V)∑_q x_qQ_q                       (प्रवाह सदिश)

2. आंतरिक बल सूत्रण (Love-Weber प्रकार)

बल और बल आघूर्ण संतुलन के माध्यम से, प्रसिद्ध Love-Weber सूत्र का सामान्यीकरण प्राप्त होता है:

σ^mac ≈ (1/V)∑_e A_el_ee_N^e ⊗ t^e
μ^mac ≈ σ^mac ⊗ x^mac : E + (1/V)∑_e A_el_ee_N^e ⊗ [m^e + E:(x_c^e ⊗ t^e)]
a^mac ≈ (1/V)∑_d S_dh_de_λ^d j_d

3. सीमा-त्रिज्या अंतराल पद

आंतरिक बल सूत्रण बाहरी बलों की वास्तविक स्थिति को नजरअंदाज करता है, "सीमा-त्रिज्या अंतराल" त्रुटि पद का परिचय देता है, जो बड़ी संख्या में कणों वाली प्रणालियों के लिए नजरअंदाज किया जा सकता है।

प्रायोगिक सेटअप

सत्यापन रणनीति

पत्र दो प्रकार की सत्यापन विधियों का उपयोग करता है:

1. Poisson समस्या सत्यापन

  • ज्यामिति: 1×1 m² वर्गाकार डोमेन
  • असतत करण: लगभग 48,000 नोड्स, 72,000 नलिका तत्व
  • परीक्षण मामले:
    • स्थिर प्रवाह पैच परीक्षण: p = 2(2x-1) + 2(2y-1), a = (-4,-4)^T
    • अरैखिक विश्लेषणात्मक समाधान: p = 16x₁x₂(1-x₁)(1-x₂)

2. यांत्रिकी समस्या सत्यापन

  • ज्यामिति: कैंटिलीवर बीम, गहराई D=3m, लंबाई S=12m
  • असतत करण: 1,034,196 स्वतंत्रता की डिग्री
  • सामग्री पैरामीटर: E₀=40GPa, α=0.3, β=0~10⁵
  • तुलना मॉडल:
    • असतत मॉडल
    • समांगीकृत Cosserat मॉडल
    • सरलीकृत Cosserat मॉडल
    • Cauchy निरंतर पदार्थ मॉडल

मूल्यांकन संकेतक

  • सटीकता मूल्यांकन: विश्लेषणात्मक समाधान से विचलन
  • अभिसरण: विभिन्न जाल आकारों के तहत परिणाम
  • सामंजस्य: आंतरिक बल और बाहरी बल सूत्रों की तुलना
  • सीमा प्रभाव: सीमा-त्रिज्या अंतराल पद का प्रभाव

प्रायोगिक परिणाम

मुख्य परिणाम

1. Poisson समस्या सत्यापन

  • रूप I (सीमा-त्रिज्या अंतराल पद के साथ): हमेशा सटीक समाधान प्रदान करता है
  • रूप II (सीमा-त्रिज्या अंतराल पद को नजरअंदाज करना): जाल आकार कम होने पर सटीकता में गिरावट
  • अभिसरण: बड़े आकार के जाल के लिए दोनों रूपों के परिणाम समान हैं

2. यांत्रिकी समस्या सत्यापन

स्थिर विश्लेषण:

  • असतत मॉडल समांगीकृत मॉडल के साथ अच्छी तरह मेल खाता है
  • मोड़ कठोरता β बढ़ने के साथ, प्रतिबल कम होता है जबकि युग्म प्रतिबल बढ़ता है
  • β<10³ के लिए Cosserat प्रभाव नजरअंदाज किया जा सकता है, Cauchy निरंतर पदार्थ का उपयोग किया जा सकता है

क्षणिक विश्लेषण:

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

मुख्य खोजें

1. सीमा-त्रिज्या अंतराल प्रभाव

प्रयोग सीमा-त्रिज्या अंतराल पद की महत्ता को सत्यापित करता है:

  • छोटे आकार के नियंत्रण आयतन के लिए, यह पद नजरअंदाज नहीं किया जा सकता
  • बड़ी संख्या में कणों वाली प्रणालियों में यह पद का प्रभाव कम है
  • मात्रात्मक त्रुटि मूल्यांकन प्रदान करता है

2. संदर्भ बिंदु स्थिति का प्रभाव

युग्म प्रतिबल प्रदिश में संदर्भ बिंदु स्थिति पद के सैद्धांतिक पूर्वानुमान को सत्यापित करता है:

  • जब स्थूल बिंदु डोमेन के द्रव्यमान केंद्र पर होता है, तो अतिरिक्त पद लुप्त हो जाता है
  • अन्य स्थितियों में अतिरिक्त युग्म प्रतिबल योगदान उत्पन्न होता है

3. सामग्री पैरामीटर का प्रभाव

  • β पैरामीटर Cosserat विशेषता लंबाई ℓ को महत्वपूर्ण रूप से प्रभावित करता है
  • β=0 के लिए ℓ≈4.6mm, β=10⁵ के लिए ℓ≈0.7m
  • Lamé पैरामीटर λ,μ अपेक्षाकृत स्थिर हैं

संबंधित कार्य

शास्त्रीय सैद्धांतिक आधार

  • Love-Weber सूत्र: असतत प्रणाली प्रतिबल प्रदिश का शास्त्रीय सूत्रण
  • अनंतक विस्तार समांगीकरण: Forest आदि का अनंतक RVE सिद्धांत
  • आभासी कार्य सिद्धांत: Chang और Kuhn की ऊर्जा समतुल्यता विधि

इस पत्र की नवीनता

  1. परिमित आकार RVE: अनंतक सिद्धांत की तुलना में आकार प्रभावों पर विचार करता है
  2. संपूर्ण युग्म प्रतिबल संभालना: सूक्ष्म प्रतिबल के दूरी पर योगदान को शामिल करता है
  3. एकीकृत ढांचा: निरंतर और असतत प्रणालियों को एक साथ संभालता है
  4. दो प्रकार के सूत्रण: आंतरिक बल और बाहरी बल की संपूर्ण पत्राचार

निष्कर्ष और चर्चा

मुख्य निष्कर्ष

  1. सैद्धांतिक योगदान: परिमित RVE का कठोर समांगीकरण सिद्धांत स्थापित किया, शास्त्रीय Love-Weber सूत्र को सामान्यीकृत किया
  2. विधि नवीनता: युग्म प्रतिबल प्रदिश संदर्भ बिंदु निर्भरता को संभालने की संपूर्ण विधि प्रस्तावित की
  3. व्यावहारिक मूल्य: बहु-पैमाने मॉडलिंग के लिए अधिक सटीक उपकरण प्रदान करता है

सीमाएं

  1. आभासी विकृति मोड का चयन: कुछ उच्च-क्रम आभासी कार्य संतुलन को नजरअंदाज किया जाता है
  2. सीमा शर्त सीमा: सभी सीमाओं के लिए स्थिर बल सीमा शर्त की आवश्यकता होती है
  3. कम्प्यूटेशनल जटिलता: सीमा-त्रिज्या अंतराल पद कम्प्यूटेशनल बोझ बढ़ाता है

भविष्य की दिशाएं

  1. अरैखिक विस्तार: अरैखिक संवैधानिक संबंधों तक सामान्यीकरण
  2. युग्मित समस्याएं: यांत्रिकी-परिवहन युग्मित समस्याओं तक विस्तार
  3. उच्च-क्रम सिद्धांत: उच्च-क्रम आभासी विकृति मोड पर विचार करना

गहन मूल्यांकन

शक्तियां

  1. सैद्धांतिक कठोरता: आभासी कार्य सिद्धांत पर आधारित कठोर गणितीय व्युत्पत्ति
  2. व्यावहारिकता: मौजूदा असतत और निरंतर मॉडलों पर सीधे लागू होता है
  3. पर्याप्त सत्यापन: विश्लेषणात्मक समाधान, संख्यात्मक तुलना और पैरामीटर अध्ययन शामिल
  4. स्पष्ट लेखन: तर्क स्पष्ट, गणितीय व्यंजक सटीक

कमियां

  1. लागू क्षेत्र: मुख्य रूप से रैखिक लोचदार समस्याओं के लिए, अरैखिक विस्तार को आगे के अनुसंधान की आवश्यकता है
  2. कम्प्यूटेशनल दक्षता: संपूर्ण सूत्र की कम्प्यूटेशनल लागत अधिक है
  3. पैरामीटर संवेदनशीलता: कुछ परिणाम सामग्री पैरामीटर β के प्रति अत्यधिक संवेदनशील हैं

प्रभाव

  1. शैक्षणिक मूल्य: बहु-पैमाने यांत्रिकी के लिए महत्वपूर्ण सैद्धांतिक आधार प्रदान करता है
  2. इंजीनियरिंग अनुप्रयोग: कंक्रीट, मिट्टी आदि विषमांगी सामग्रियों के मॉडलिंग के लिए सीधे उपयोग किया जा सकता है
  3. सॉफ्टवेयर कार्यान्वयन: खुला स्रोत कोड प्रदान किया गया है, प्रसार को सुविधाजनक बनाता है

लागू परिस्थितियां

  • दानेदार सामग्री का बहु-पैमाने मॉडलिंग
  • विषमांगी समग्र सामग्री का समांगीकरण
  • सूक्ष्म घूर्णन वाली सामग्री प्रणालियां
  • आकार प्रभावों पर विचार करने की आवश्यकता वाली इंजीनियरिंग समस्याएं

संदर्भ

पत्र 63 संबंधित संदर्भों का हवाला देता है, मुख्य रूप से शामिल हैं:

  • Cosserat निरंतर पदार्थ सिद्धांत के शास्त्रीय साहित्य
  • असतत तत्व विधि और समांगीकरण सिद्धांत
  • बहु-पैमाने मॉडलिंग के संबंधित कार्य
  • सत्यापन के लिए उपयोग की गई संख्यात्मक विधि साहित्य

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