DOI: https://doi.org/10.5194/tc-20-595-2026
تاريخ النشر: 2026-01-23
المؤلف: Faranak Sahragard وآخرون
الموضوع الرئيسي: دراسات وملاحظات الكريوسفير
نظرة عامة
تقدم هذه الدراسة نموذجًا شاملًا ثلاثي الأبعاد للخصائص الميكانيكية للجليد، يتناول التداخلات المعقدة بين درجة الحرارة، الضغط، ومعدل التشوه، والتي تُعرف مجتمعة باسم الديناميكا الحرارية للثلج. يعزز النموذج الإطار الحالي للماكسويل المرن الهش من خلال تقديم تعبير موحد من نوع أرهينيوس الذي يحكم كل من اللزوجة وعتبات بدء التلف. يسمح هذا النهج بتمثيل العمليات الحرارية الميكانيكية الرئيسية باستخدام مجموعة بسيطة من المعلمات القابلة للتفسير الفيزيائي. تُظهر معايرة النموذج والتحقق من صحته ضد بيانات اختبارات الضغط الثلاثي والاسترخاء على الجليد متعدد البلورات قدرته على التقاط السلوكيات الميكانيكية الحرجة بدقة، بما في ذلك الانتقالات من الهش إلى المرن، وقوة تعتمد على معدل التشوه، والتليين الحراري.
تشير النتائج الرئيسية إلى أن النموذج يعيد إنتاج استجابة الإجهاد-التشوه للجليد بفعالية، بما في ذلك أحجام الإجهاد القصوى والمتبقية، بينما يعكس التأثيرات المثبطة للضغط المحصور العالي على تطور التلف وحساسية قوة الجليد ومرونته لتغيرات درجة الحرارة. بالإضافة إلى ذلك، يتضمن النموذج آلية شفاء جديدة مستوحاة من التلبيد اللزج، تربط بين معدل استعادة التلف وتفريغ الطاقة اللزجة وتعديلها بواسطة الضغط ودرجة الحرارة. تشير المحاكاة الأولية إلى أن هذه العملية الشفائية تؤدي إلى استعادة تدريجية للصلابة أثناء الاسترخاء، مما يوفر رؤى حول السلوك طويل الأمد واستعادة التلف لهياكل الجليد.
مقدمة
تؤكد مقدمة هذه الورقة البحثية على الحاجة الملحة لنمذجة دقيقة لميكانيكا الجليد عبر مجالات مختلفة، بما في ذلك علم الجليد والهندسة الجيوتقنية. تنشأ تعقيدات سلوك الجليد من التفاعل المعقد بين الضغط، معدل التشوه، ودرجة الحرارة، مما يعقد تطوير نماذج خصائص ميكانيكية قوية. على عكس المواد الصلبة الطبيعية الأخرى، تؤدي التفاعلات غير الخطية بين هذه المعلمات في الجليد إلى تحديات في عزل الآليات الفردية للتحليل. أظهرت الدراسات التجريبية أن معدل التشوه يؤثر ليس فقط على الزحف والاسترخاء ولكن أيضًا على إجهاد الخضوع وأنماط الفشل، بينما يؤثر الضغط ودرجة الحرارة بشكل كبير على الخصائص الميكانيكية مثل الصلابة واللزوجة.
تستعرض الورقة النماذج الحالية للخصائص الميكانيكية للجليد، مشيرة إلى أن العديد منها تم تكييفه من الأطر المستخدمة في المواد السائلة، الحبيبية، والبلورية. تشمل النماذج البارزة نموذج جسم بورجر، الذي يلتقط السلوك اللزج المرن، وتكييفات مختلفة تتضمن آليات التلف. ومع ذلك، غالبًا ما تكافح هذه النماذج لدمج تأثيرات الإجهاد، معدل التشوه، ودرجة الحرارة بشكل متماسك، مما يحد من قابليتها للتطبيق عبر ظروف متنوعة. يقترح المؤلفون نموذجًا جديدًا ثلاثي الأبعاد للخصائص الميكانيكية، يُعرف باسم الديناميكا الحرارية للجليد، والذي يبني على إطار الماكسويل المرن الهش. يهدف هذا النموذج إلى معالجة العيوب السابقة من خلال دمج اللزوجة المعتمدة على الإجهاد ودرجة الحرارة، وتلف مدفوع بالطاقة، وحدود بدء تلف تعتمد على درجة الحرارة ومعدل التشوه، مع الحفاظ على معلمات قابلة للإدارة. تشير المقارنات الأولية مع البيانات التجريبية إلى أن هذا النموذج يلتقط بفعالية الميزات الرئيسية لسلوك الجليد عبر مجموعة من الظروف.
مناقشة
في هذا القسم، يناقش المؤلفون الإطار الرياضي وصياغة النموذج لتكوين ماكسويل اللزج المطبق على الجليد البحري، مع دمج آليات التلف. يعرفون الرموز الرياضية الرئيسية، مثل استخدام الحروف الغامقة للموترات من الدرجة الثانية والحروف المجوفة للموترات من الدرجة الرابعة، ويقدمون مفاهيم التشوهات والإجهادات الحجمية والانحرافية. يبني النموذج على الأعمال السابقة لـ Dansereau وآخرين (2016) ويتضمن معيار بدء التلف بناءً على معيار فشل موهر-كولوم، الذي يتأثر بدرجة الحرارة ومعدل التشوه. يؤكد المؤلفون على الطبيعة الديناميكية الحرارية لبدء التلف والعلاقة بين الآليات اللزجة والمرنة، مشيرين إلى أن التشوه اللزج بلا أثر وانحرافي.
تتناول المناقشة أيضًا اللزوجة المعتمدة على الإجهاد للجليد، معتمدة على قانون تدفق جلين-ناي، الذي تم تعديله ليأخذ في الاعتبار تدهور التلف. يقدم المؤلفون معيار فشل بيضاوي يتضمن اعتماديات على درجة الحرارة ومعدل التشوه، مما يسمح بتنبؤ أكثر دقة بقوة الجليد تحت ظروف متغيرة. يتحققون من صحة النموذج ضد البيانات التجريبية، مما يوضح قدرته على التقاط الانتقال من الفشل الهش إلى الفشل المرن وتأثيرات الضغط المحصور ومعدل التشوه على تطور التلف. تختتم هذه القسم باستكشاف بارامتري لتنبؤات النموذج تحت درجات حرارة ومعدلات تشوه مختلفة، موضحة قوته في محاكاة مسارات التحميل المعقدة وآلية الشفاء الكامنة في الجليد.
القيود
تسلط القيود في الدراسة الضوء على عدة افتراضات مبسطة تقيد قابلية تطبيق الصياغة المطورة. بشكل أساسي، يفترض النموذج التماثل، متجاهلاً خصائص الصلابة والتدفق المعتمدة على الاتجاه للجليد متعدد البلورات، والتي تتأثر بالنسيج البلوري ومحاذاة الشقوق الدقيقة. بالإضافة إلى ذلك، من المتوقع أن تتناقص دقة النموذج التنبؤية بالقرب من نقطة انصهار الضغط بسبب الآليات المنشطة حراريًا مثل عمليات حدود الحبيبات وإعادة التبلور، التي يمكن أن تسرع الزحف بما يتجاوز سلوك أرهينيوس. كما أن العوامل الميكروهيكلية الأخرى، بما في ذلك الملوحة، محتوى المحلول الملحي، المسامية، وتطور حجم الحبيبات، لم تؤخذ في الاعتبار، على الرغم من تأثيرها الكبير على الصلابة الفعالة وآليات الزحف.
علاوة على ذلك، فإن تكوين ماكسويل المستخدم في الدراسة لا يمثل بشكل كافٍ المرحلة الأولية للزحف، التي ترتبط بالتشوه القابل للاستعادة جزئيًا وحقول الإجهاد الداخلية المتطورة بدلاً من التدفق اللزج فقط. تظل آلية الشفاء المدفوعة باللزوجة غير مضبوطة وتحتاج إلى تحقق تجريبي. تم تقييد مقارنة السلوك الثلاثي المحاور بمستويات ضغط متوسطة عالية (حوالي 10 ميجا باسكال)، ذات صلة بالضغوط الأساسية تحت الأنهار الجليدية، بسبب ندرة البيانات التجريبية. يجب أن تركز الأبحاث المستقبلية على تجارب مستهدفة، مثل اختبارات التحميل-الاحتفاظ-إعادة التحميل واختبارات ثلاثية المحاور بتردد منخفض، للتحقق من الآليات المقترحة وتعزيز قابلية تطبيق النموذج عبر سياقات مختلفة في علم الجليد، والهندسة الجيوتقنية، خاصة تحت ظروف التحميل طويل الأمد والتشوه الدوري.
DOI: https://doi.org/10.5194/tc-20-595-2026
Publication Date: 2026-01-23
Author(s): Faranak Sahragard et al.
Primary Topic: Cryospheric studies and observations
Overview
This research presents a comprehensive three-dimensional constitutive model for ice, addressing the complex interdependencies of temperature, pressure, and strain rate, collectively termed the thermobarokinetics of ice. The model enhances the existing Maxwell elasto-brittle framework by introducing a unified Arrhenius-type expression that governs both viscosity and damage initiation thresholds. This approach allows for the representation of key thermomechanical processes using a minimal set of physically interpretable parameters. The model’s calibration and validation against triaxial compression and relaxation test data on polycrystalline ice demonstrate its capability to accurately capture critical mechanical behaviors, including ductile to brittle transitions, strain-rate-dependent strength, and thermal softening.
Key findings indicate that the model effectively reproduces the stress-strain response of ice, including peak and residual stress magnitudes, while reflecting the suppressive effects of high confining pressure on damage development and the sensitivity of ice strength and ductility to temperature variations. Additionally, the model incorporates a novel healing mechanism inspired by viscous sintering, linking the rate of damage recovery to viscous energy dissipation and modulated by pressure and temperature. Preliminary simulations suggest that this healing process leads to a gradual restoration of stiffness during relaxation, offering insights into the long-term behavior and damage recovery of ice structures.
Introduction
The introduction of this research paper emphasizes the critical need for accurate modeling of ice mechanics across various fields, including glaciology and geotechnical engineering. The complexity of ice behavior arises from the intricate interplay of pressure, strain rate, and temperature, which complicates the development of robust constitutive models. Unlike other natural solids, the nonlinear interactions among these parameters in ice lead to challenges in isolating individual mechanisms for analysis. Experimental studies have demonstrated that strain rate affects not only creep and relaxation but also yield stress and failure modes, while pressure and temperature significantly influence mechanical properties such as stiffness and viscosity.
The paper reviews existing constitutive models for ice, noting that many have been adapted from frameworks used in fluid, granular, and crystalline materials. Notable models include the Burgers body model, which captures viscoelastic behavior, and various adaptations that incorporate damage mechanisms. However, these models often struggle to cohesively integrate the effects of stress, strain rate, and temperature, limiting their applicability across diverse conditions. The authors propose a new three-dimensional constitutive model, termed the thermobarokinetics of ice, which builds on the Maxwell elasto-brittle framework. This model aims to address previous shortcomings by incorporating stress and temperature-dependent viscosity, energy-driven damage, and a temperature and rate-dependent damage initiation limit, while maintaining manageable parameters. Preliminary comparisons with experimental data suggest that this model effectively captures key features of ice behavior across a range of conditions.
Discussion
In this section, the authors discuss the mathematical framework and model formulation for a viscoelastic Maxwell configuration applied to sea ice, incorporating damage mechanisms. They define key mathematical notations, such as the use of bold characters for second-order tensors and hollow characters for fourth-order tensors, and introduce the concepts of volumetric and deviatoric strains and stresses. The model builds upon previous work by Dansereau et al. (2016) and includes a damage initiation criterion based on the Mohr-Coulomb failure criterion, which is influenced by temperature and strain rate. The authors emphasize the thermokinetic nature of the damage initiation and the relationship between viscous and elastic mechanisms, highlighting that the viscous strain is traceless and deviatoric.
The discussion also covers the stress-dependent viscosity of ice, adopting the Glen-Nye flow law, which is modified to account for damage degradation. The authors present an elliptical failure criterion that incorporates temperature and strain rate dependencies, allowing for a more accurate prediction of ice strength under varying conditions. They validate the model against experimental data, demonstrating its ability to capture the transition from brittle to ductile failure and the effects of confining pressure and strain rate on damage evolution. The section concludes with a parametric exploration of the model’s predictions under different temperatures and strain rates, illustrating its robustness in simulating complex loading paths and the healing mechanism inherent in ice.
Limitations
The limitations of the study highlight several simplifying assumptions that restrict the applicability of the developed formulation. Primarily, the model assumes isotropy, neglecting the direction-dependent stiffness and flow characteristics of polycrystalline ice, which are influenced by crystallographic fabric and microcrack alignment. Additionally, the model’s predictive accuracy is expected to diminish near the pressure melting point due to thermally activated mechanisms such as grain boundary processes and recrystallization, which can accelerate creep beyond Arrhenius behavior. Other microstructural factors, including salinity, brine content, porosity, and grain size evolution, are also not accounted for, despite their significant impact on effective stiffness and creep mechanisms.
Furthermore, the Maxwell configuration employed in the study does not adequately represent the primary creep phase, which is associated with partially recoverable deformation and evolving internal stress fields rather than solely viscous flow. The proposed viscosity-driven healing mechanism remains uncalibrated and requires experimental validation. The comparison of triaxial behavior has been limited to high mean stress levels (approximately 10 MPa), relevant to basal overburden pressures beneath glaciers, due to a scarcity of experimental data. Future research should focus on targeted experiments, such as load-hold-reload and low-frequency cyclic triaxial tests, to verify the proposed mechanisms and enhance the model’s applicability across various cryospheric, glaciological, and geotechnical contexts, particularly under long-term loading and cyclic deformation scenarios.