Summary
ويصف البروتوكول الإجراءات اللازمة للحصول علي صور التصوير المقطعي (CT) عاليه الاستبانة المكانية للتربة الحبيبية اثناء الضغط ثلاثي المحاور ، ولتطبيق تقنيات معالجه الصور علي هذه الصور المقطعية لاستكشاف السلوك الميكانيكي علي نطاق الحبوب التربة تحت التحميل.
Abstract
وقد مكن التطور السريع لتقنيات التصوير بالاشعه السينية مع مهارات معالجه الصور وتحليلها من اقتناء صور مقطعيه للتربة الحبيبية ذات الاستبانات المكانية العالية. واستنادا إلى هذه الصور المقطعية ، يمكن التحقيق كميا في السلوك الميكانيكي علي نطاق الحبوب مثل الجسيمات الجزيئية (اي ترجمات الجسيمات وتناوب الجسيمات) وتوطين السلالات وتطور الاتصال بين الجسيمات من التربة الحبيبية. ومع ذلك ، لا يمكن الوصول اليها باستخدام الأساليب التجريبية التقليدية. وتبين هذه الدراسة استكشاف السلوك الميكانيكي علي نطاق الحبوب من عينه التربة الحبيبية تحت ضغط ثلاثي المحاور باستخدام السنكروتروني الاشعه السينية المجهرية التصوير المقطعي (μct). مع هذا الأسلوب ، يتم استخدام جهاز التحميل مصغره ملفقه خصيصا لتطبيق الضغوط الحصر والمحورية للعينه خلال اختبار ثلاثي المحاور. يتم تركيب الجهاز في التصوير المقطعي الاشعه السينية السنكروتروني بحيث يمكن جمع الصور المقطعية الدقة العالية المكانية للعينه في مراحل تحميل مختلفه من الاختبار دون اي إزعاج للعينه. مع القدرة علي استخراج المعلومات علي النطاق الكلي (علي سبيل المثال ، والضغوط عينه الحدود والسلالات من الاعداد الجهاز ثلاثي المحاور) ومقياس الحبوب (علي سبيل المثال ، حركات الحبوب والتفاعلات الاتصال من الصور المقطعية) ، وهذا الاجراء يوفر منهجيه فعاله للتحقيق في الميكانيكا متعددة المستويات من التربة الحبيبية.
Introduction
ومن المسلم به علي نطاق واسع ان الخواص الميكانيكية علي النطاق الكلي للتربة الحبيبية ، مثل الصلابة وقوه القص والنفاذيه ، هي ذات اهميه حاسمه بالنسبة للعديد من الهياكل الجيوتقنيه ، علي سبيل المثال ، الأسس والمنحدرات والسدود المليءه بالصخور. لسنوات عديده ، والاختبارات في الموقع والاختبارات المختبرية التقليدية (علي سبيل المثال ، اختبارات ضغط أحاديه البعد ، اختبارات ضغط ثلاثي المحاور واختبارات نفاذيه) وقد استخدمت لتقييم هذه الخصائص في التربة المختلفة. كما وضعت رموز ومعايير لاختبار الخواص الميكانيكية للتربة لأغراض هندسية. في حين ان هذه الخصائص الميكانيكية علي نطاق واسع قد درست بشكل مكثف ، والسلوك الميكانيكي علي نطاق الحبوب (علي سبيل المثال ، الجسيمات الجزيئية ، والتفاعل الاتصال وتوطين سلاله) التي تحكم هذه الخصائص قد اجتذبت اهتمام اقل بكثير من المهندسين والباحثين. أحد الأسباب هو عدم وجود طرق تجريبية فعاله متاحه لاستكشاف السلوك الميكانيكي علي نطاق الحبوب للتربة.
وحتى الآن ، فان معظم فهم السلوك الميكانيكي علي نطاق الحبوب من التربة الحبيبية قد جاء من النمذجة عنصر منفصل1 (DEM) ، وذلك بسبب قدرته علي استخراج المعلومات علي نطاق الجسيمات (علي سبيل المثال ، الجسيمات والاتصال الجسيمات قوات). في الدراسات السابقة لاستخدام تقنيات DEM لنمذجة السلوكيات الميكانيكية التربة الحبيبية ، وكان كل الجسيمات الفردية تمثل ببساطه من قبل دائره واحده أو المجال في النموذج. وقد ادي استخدام مثل هذه الاشكال الجسيمات الإفراط في التبسيط إلى الإفراط في دوران الجزيئات التالي انخفاض قوه الذروة السلوك2. لتحقيق أداء أفضل النمذجة ، وقد استخدمت العديد من المحققين نموذج مقاومه المتداول3،4،5،6 أو الجسيمات غير النظامية الاشكال7،8، 9،10،11،12 في المحاكاة الخاصة بهم DEM. ونتيجة لذلك ، تم الحصول علي فهم أكثر واقعيه للسلوك الحركي الجسيمات. وبصرف النظر عن الجسيمات الجزيئية ، فقد استخدمت DEM بشكل متزايد للتحقيق في تفاعل الحبوب وتطوير النماذج النظرية. ومع ذلك ، وبسبب الحاجة إلى أعاده إنتاج اشكال الجسيمات الحقيقية واستخدام نماذج الاتصال المتطورة ، DEM يتطلب قدره حسابيه عاليه للغاية في نمذجة التربة الحبيبية مع الاشكال غير النظامية.
وفي الاونه الاخيره ، أتاح تطوير المعدات البصرية وتقنيات التصوير (علي سبيل المثال ، المجهر ، والاشعه المقطعية بالليزر ، والصور المقطعية المحوسبة للاشعه السينية ، والاشعه المقطعية الدقيقة) العديد من الفرص للفحص التجريبي السلوك الميكانيكي علي نطاق الحبوب من التربة الحبيبية. عن طريق الحصول علي وتحليل الصور عينه التربة قبل وبعد اختبار ثلاثي المحاور ، وقد استخدمت هذه المعدات والتقنيات في التحقيق في الهياكل المجهرية التربة13،14،15،16 و17و18و19. في الاونه الاخيره ، في الموقع الاختبارات مع الاشعه السينية ct أو μct وقد استخدمت بشكل متزايد للتحقيق في تطور نسبه الفراغ20، وتوزيع سلاله21،22،23،24، حركه الجسيمات25،26،27،28، الاتصال بين الجسيمات29،30،31 والجسيمات سحق32 من التربة الحبيبية. هنا ، "في الموقع" يعني المسح بالاشعه السينية التي أجريت في نفس الوقت التحميل. وعلي النقيض من المسح الضوئي العام للاشعه السينية ، فان اختبارات المسح بالاشعه السينية في الموقع تتطلب جهاز تحميل ملفق خصيصا لتسليم الإجهاد إلى عينات التربة. مع الاستخدام المشترك للجهاز التحميل والاشعه السينية CT أو جهاز μCT ، يمكن الحصول علي صور CT من العينات في مراحل تحميل مختلفه من الاختبارات غير مدمره. واستنادا إلى هذه الصور المقطعية ، يمكن الحصول علي ملاحظات علي مقياس الجسيمات لسلوك التربة الحبيبية. هذه الملاحظات علي مستوي الجسيمات المستندة إلى صوره CT مفيده للغاية للتحقق من النتائج العددية واكتساب رؤى جديده في السلوك الميكانيكي علي نطاق الحبوب من التربة الحبيبية.
تهدف هذه المقالة إلى تبادل التفاصيل حول كيفيه اجراء فحص الاشعه السينية في الموقع من عينه التربة ، وذلك باستخدام تجربه مثاليه التي تلاحظ الجسيمات الجزيئية ، وتوطين سلاله وتطور الاتصال بين الجسيمات داخل عينه التربة. وتبين النتائج ان الاشعه السينية في الموقع اختبارات المسح لديها إمكانات كبيره لاستكشاف السلوك علي مستوي الحبوب من التربة الحبيبية. ويغطي البروتوكول اختيار جهاز الاشعه السينية μCT واعداد أجهزه التحميل ثلاثية المحاور المصغرة ، ويتم توفير إجراءات مفصله لتنفيذ الاختبار. الاضافه إلى ذلك ، فان الخطوات التقنية لاستخدام معالجه الصور وتحليلها لتحديد كميات الجسيمات الجزيئية (اي ترجمه الجسيمات ودوران الجسيمات) ، وتوطين السلالات ، وتطور الاتصال بين الجسيمات (اي كسب الاتصال ، وفقدان الاتصال اتصلت حركه) من التربة وصفت.
Subscription Required. Please recommend JoVE to your librarian.
Protocol
1. تصميم التجربة مقدما بشكل جيد
- تحديد المواد الاختبار ، وحجم الجسيمات ، وحجم العينة وعينه المساميه الاوليه.
ملاحظه: يستخدم الرمل لايتون بونار مع قطر 0.15 ~ 0.30 mm وحجم عينه من 8 × 16 ملم (قطر x الارتفاع) كمثال لإثبات بروتوكول هذه الدراسة. ويمكن أيضا ان تستخدم الرمال الأخرى مثل الرمل فوجيان ، هيوستن الرمل ، اوتاوا الرمل وكايكوس ooids ، الخ واحجام عينه مماثله. - اختيار كاشف المناسبة (الشكل 1ا) وفقا لدقه المكانية المطلوبة ومنطقه المسح الضوئي ، والتي يتم تحديدها وفقا لحجم الجسيمات سلفا وحجم العينة. علي سبيل المثال ، يتم استخدام كاشف مع دقه المكانية من 6.5 μm في هذه الدراسة. لديها مساحة مسح فعاله من 2048 x 860 بكسل (اي 13.3 × 5.6 مم).
ملاحظه: اثناء اختبار ضغط ثلاثي المحاور ، يجب ان تبقي العينة المشوهة في منطقه المسح الضوئي للكشف. وينبغي استخدام كاشف الاستبانة المكانية العالية بحيث تحتوي الجزيئات الفردية علي voxels كافيه لاستخراج خصائص الجسيمات بالشكل المناسب. - تحديد الطاقة المطلوبة من مصدر الاشعه السينية (الشكل 1ا) ووقت التعرض وفقا لماده الاختبار وحجم العينة. عموما ، ينبغي استخدام طاقة اعلي لعينه أكبر تتالف من ماده أكثر كثافة. استخدام طاقة الاشعه السينية من 25 كيلو ووقت التعرض 0.05 s لعينات الرمل في هذه الدراسة.
ملاحظه: يمكن تحديد طاقة الاشعه السينية المطلوبة ووقت التعرض بواسطة التجربة والخطا باستخدام إسقاط ممسوح ضوئيا للعينه. وينبغي الا تقل نسبه كثافة الحجم الرمادي الأدنى للإسقاط إلى قيمته القصوى عن 0.2. والا ، ينبغي استخدام طاقة الاشعه السينية الأعلى أو وقت التعرض الأطول. - تحديد سرعه الدوران المطلوبة ω (درجات في الثانية) لمرحله الدوران (الشكل 1ا) من جهاز الاشعه السينية. يتم احتساب سرعه دوران ω وفقا للعدد المطلوب من الإسقاطات N (علي سبيل المثال ، N = 1,080) لأعاده بناء شريحة CT.
ملاحظه: ω = 180 Vs/N. هنا ، Vs هو سرعه المسح الضوئي للجهاز الاشعه السينية ، اي عدد الاشعه الممسوحة ضوئيا وسجلت في الثانية الواحدة. يتاثر Vs بشكل رئيسي بأداء الكاشف والاجهزه المرتبطة بالكاشف مثل الكمبيوتر. - افتعال جهاز تحميل ثلاثي المحاور (الأرقام 1B ، C، انظر أيضا المرجع 33) لاستخدامها بالاقتران مع جهاز الاشعه السينية μct. يجب ان يكون الجهاز نفس الوظائف الرئيسية كجهاز ضغط ثلاثي محوري التقليدية. يجب ان التصميم النظر في متطلبات حجم العينة ، ومجموعه من الضغوط حصر ومعدلات التحميل.
ملاحظه: يجب ان يكون الجهاز قادرا علي احتواء جهاز الاشعه السينية μCT ويكون خفيفا لتسهيل دورانه باستخدام مرحله الدوران. يجب ان تكون الخلية ثلاثية المحاور شفافة للاشعه السينية. النظر إلى متطلبات الشفافية ، يمكن استخدام الأكريليك والبولي لافتعال الخلية ثلاثية المحاور. - اجراء اختبار مع نفس الضغط الحصر ، وسرعه التحميل وخصائص العينة (اي المواد ، وحجم العينة والمساميه الاولي) خارج الاشعه السينية CT الماسح الضوئي للتخطيط عند إيقاف التحميل للمسح المقطعي المحوسب.
2. اجراء اختبار ضغط ثلاثي المحاور في الموقع
- وضع معدات التحميل ثلاثي المحاور ومواد الاختبار في الموقع.
ملاحظه: يتم وضع جهاز التحميل وحصر الضغط التي تقدم الاجهزه (انظر جدول المواد) في غرفه المسح المقطعي الاشعه السينية ، في حين ان الحصول علي البيانات والسيطرة علي الاجهزه تقع خارج. ثم يتم تشغيل التحميل ثلاثي المحاور والمسح المقطعي المحوسب للعينه خارج غرفه المسح الضوئي. - إصلاح مرحله رفع علي متن جهاز الاشعه السينية الصغرى CT (الشكل 1ب). إصلاح مرحله أماله علي مرحله رفع ومرحله التناوب علي مرحله أماله ، علي التوالي (الشكل 1ب).
ملاحظه: يجب ان يكون مرحله رفع وأماله المرحلة قدره التحميل كافيه لنقل المعدات ذات الصلة وضعت عليها. - اضبط موضع واتجاه مرحله الاستدارة عبر مرحله الاماله بحيث يمر اي أشعه سينية واحده عبر نفس النقاط داخل العينة عندما يتم تدويرها عبر 180 درجه حول محور مرحله التدوير.
ملاحظه: الخطوات 2.2 إلى 2.3 قابله للتطبيق علي جهاز CT الاشعه السينية الصغرى في مركز الإشعاع شنغهاي Synchrotron (SSRF). بالنسبة للاجهزه المقطعية الدقيقة للاشعه السينية المستخدمة خصيصا للاختبار ثلاثي المحاور في الموقع ، يمكن حذف هذه الخطوات بعد التحديد الدقيق للموقع وتثبيت مرحله التدوير. - اعداد عينه التربة علي متن السفينة وفقا للإجراءات التالية.
- أضافه كميه صغيره من الشحوم سيليكون حول السطح الجانبي من الطرف العلوي من لوحه قاعده ووضع حجر مساميه علي سطحه العلوي. وضع غشاء حول السطح الجانبي للنهاية العليا (الشكل 2ا).
- أضافه كميه صغيره من الشحوم سيليكون علي أسطح الاتصال بين جزاين من صانع عينه وقفل عليه. وضع صانع عينه علي لوحه قاعده والسماح لغشاء لتمرير من خلال ذلك (الشكل 2ب).
- إنشاء شفط (علي سبيل المثال ، 25 كيلو باسكال) داخل صانع عينه من خلال فوهه لها باستخدام مضخة فراغ. إصلاح الغشاء إلى السطح الجانبي من نهايته العليا. تاكد من إرفاق الغشاء بالسطح الداخلي لصانع العينة (الشكل 2ج).
- إسقاط المواد الحبيبية اختبار من ارتفاع معين في صانع عينه باستخدام قمع حتى يتم ملؤها تماما. وينبغي ان يكون السطح العلوي لعينه التربة هو نفس مستوي الحافة العليا لصانع العينة (الشكل 2د).
- وضع حجر مساميه آخر علي راس عينه التربة ، ولوحه وساده الفولاذ المقاوم للصدا علي راس الحجر مساميه. تطبيق بعض الشحوم سيليكون حول السطح الجانبي للوحه وساده. أزاله الجانب العلوي من الغشاء من صانع عينه وإصلاحه إلى لوحه وساده (الشكل 2ه).
- أزاله شفط داخل فوهه صانع عينه وخلق شفط داخل صمام علي لوحه قاعده. وأخيرا ، قم بازاله صانع العينة. ويتم إنتاج عينه جافه مصغره ، كما يتبين من الشكل 2واو.
ملاحظه: توضح هذه الخطوة الاجراء الخاص بإنتاج عينه تربه مصغره باستخدام أسلوب pluviation الهواء. ويمكن أيضا استخدام طريقه الضغط الجاف التقليدية لإنتاج العينة.
- إصلاح الخلية حصر علي لوحه قاعده وإصلاح لوحه اعلي الغرفة في الأعلى من الخلية حصر (الشكل 1ج).
- إصلاح رمح المكبس من الخلية علي لوحه اعلي الغرفة (الشكل 1C).
- ضع اللوحة الاساسيه مع الخلية التي تم حصرها واللوحة العلوية للغرفة في مرحله التدوير. ويستخدم اطار لضبط ارتفاع العينة للمسح المقطعي المحوسب (الشكل 1ب).
ملاحظه: يتم استخدام هذا الإطار نظرا لنطاق حركه محدوده من مرحله رفع في SSRF. ليست هناك حاجه لاستخدام الإطار إذا تم استخدام مرحله رفع مع نطاق حركه كبيره. - افيكس بقية جهاز التحميل علي لوحه اعلي الغرفة.
- تثبيت المحولات التفاضلية المتغيرة الخطية (LVDT) ، خليه التحميل والمحرك يخطو وتفعيلها (الشكل 1ج).
- ملء الخلية مع أزاله المياه من خلال صمام ضغط الخلية (CP) (انظر الشكل 1ج) باستخدام المياه الموردة من جهاز عرض الضغط الحصر (انظر جدول المواد). إغلاق صمام الخروج (WE) المياه (انظر الشكل 1ج) عندما يبدا الماء في التدفق من الصمام.
ملاحظه: اضبط الضغط الذي يوفره الجهاز علي وضع الضغط المستمر مع قيمه ضغط ثابت منخفضه جدا (علي سبيل المثال ، 10 كيلو باسكال). - أضافه ضغط ثابت الحصر من 25 كيلو باسكال إلى العينة وأزاله الشفط داخل العينة.
- زيادة تدريجيا الضغط الحصر إلى قيمه محدده مسبقا باستخدام الجهاز الضغط الحصر العرض.
- اجراء المسح الضوئي الأول من العينة. بالنسبة لماسح الاشعه المقطعية عالي المساحة (مثلا ، بحجم البكسل 6.5 ميكرومتر) ، فان الفحص الكامل للعينه (علي سبيل المثال ، بارتفاع 16 مم) يتطلب عاده مسح العينة علي عده ارتفاعات مختلفه (اي ان المسح ينقسم إلى عده أقسام).
ملاحظه: إذا تم استخدام كاشف دقه مكانيه منخفضه وعينه صغيره الحجم ، قد تكون منطقه المسح كافيه للحصول علي فحص حقل كامل من العينة باستخدام مقطع واحد.- مسح مقطع من العينة. تعيين الماسح الضوئي CT إلى وضع التقاط الصورة ومن ثم بدء مرحله الدوران لتدوير الجهاز بأكمله عبر 180 درجه في معدل دوران ثابت محدد مسبقا (علي سبيل المثال ، 3.33 درجه/ثانيه) للتقاط إسقاطات CT للعينه في زوايا مختلفه.
ملاحظه: ويقترح ان يتم مسح العينة من أسفل إلى اعلي (اي ، القسم الأول يحتوي علي جميع الجزيئات الموجودة في الجزء السفلي من العينة). - قم بإيقاف تشغيل وضع التقاط الصورة عند انتهاء الاستدارة. تدوير الجهاز مره أخرى إلى الموقف الاولي.
- رفع العينة مع الجهاز بأكمله حتى باستخدام مرحله رفع (الشكل 1ب) بارتفاع معين (علي سبيل المثال ، 4 مم) لمسح المقطع التالي من العينة.
ملاحظه: ينبغي ان يضمن الرفع وجود تداخل بين القسم الحالي والقسم الأخير (اي ان هناك تداخلا بين المقطعين المتتاليين). يجب ان يكون التداخل 10 بكسل علي الأقل لتسهيل خياطه منها. - كرر الخطوات 2.13.1-2.13.3 حتى يتم مسح المقطع الأخير من العينة.
- مسح مقطع من العينة. تعيين الماسح الضوئي CT إلى وضع التقاط الصورة ومن ثم بدء مرحله الدوران لتدوير الجهاز بأكمله عبر 180 درجه في معدل دوران ثابت محدد مسبقا (علي سبيل المثال ، 3.33 درجه/ثانيه) للتقاط إسقاطات CT للعينه في زوايا مختلفه.
- تطبيق تحميل محوري علي العينة مع معدل تحميل ثابت. هنا ، يتم استخدام معدل تحميل 0.2 ٪/دقيقه في هذه الدراسة. يمكن للمستخدمين تعيين معدل تحميل مختلف وفقا لمتطلبات التجربة.
- وقفه التحميل المحوري في سلاله محوريه محدده مسبقا. انتظر حتى تصل القوه المحورية المقاسه إلى قيمه ثابته (عموما في غضون 2 دقيقه) وتنفيذ الفحص التالي. إجراءات المسح الضوئي هي نفسها كما هو موضح في الخطوة 2.13.
- كرر الخطوات 2.14 و 2.15 حتى نهاية التحميل.
- تفريغ الاختبار وأزاله العينة من جهاز triaxial.
- قم بتركيب اللوحة الاساسيه والخلية التي تم حصرها في مرحله الدوران للحصول علي عده إسقاطات مسطحه (عموما 10 إسقاطات) من الكاشف. اغلق مصدر الاشعه السينية للحصول علي نفس العدد من الإسقاطات المظلمة من الكاشف.
ملاحظه: يتم استخدام الإسقاطات المسطحة والداكنة لاستعاده المرحلة من إسقاطات CT الخام. ويعزز تنفيذ التصحيح المسطح والغامق التباين بين العينة والخلفية المحيطة بها في الشرائح المقطعية المعاد بناؤها. كما انه يساعد علي تخفيف القطع الاثريه الحلقة الناتجة عن بكسل معيبه للكشف.
3-معالجه الصور وتحليلها
- معالجه الصور
- تنفيذ استرجاع المرحلة (الشكل 3ب) من إسقاطات CT الخام (الشكل 3ا) من العينة باستخدام البرمجيات الحرة pitre34. إسقاطات التحميل (بما في ذلك الإسقاطات المسطحة والداكنة) في PITRE من القائمة تحميل الصورة. انقر فوق الرمز Ppci. ادخل معلمات المسح الضوئي ذات الصلة وانقر فوق Single لتنفيذ استرداد المرحلة.
ملاحظه: يوفر تنفيذ استرجاع المرحلة تعزيز الوصلات البينية بين المراحل المختلفة (اي مرحله الفراغ والمرحلة الصلبة) في الشرائح المقطعية المعاد بناؤها ، والتي تكتسي اهميه كبيره بالنسبة للتحليل اللاحق المستند إلى الصور اتصالات بين الجسيمات. - أعاده بناء شرائح CT من العينة باستخدام PITRE استنادا إلى إسقاطات الاشعه المقطعية بعد استرجاع المرحلة (الشكل 3ج). تحميل الإسقاطات في PITRE من القائمة تحميل الصورة. انقر فوق الرمز بروجسينو. ادخل المعلمات ذات الصلة في نافذه ظهرت وانقر فوق واحد لأعاده بناء شريحة CT.
ملاحظه: تحقق من الشرائح الافقيه للتاكد من عدم وجود شعاع ثقيل تصلب التحف أو التحف الدائرية. والا فان تغيير معلمات المسح الضوئي الحالية وأعاده تفحص العينة مطلوبه. تحقق من الشرائح العمودية. إذا كانت العينة مائله بشده قبل القص ، يعتبر الاختبار غير ناجح. - تنفيذ تصفيه الصور علي شرائح CT. يستخدم فلتر نشر متباين الخواص لاجراء تصفيه الصور (الشكل 3د).
- قم باجراء عمليه التقطير بالصور علي الشرائح المقطعية المصفاة. تنفيذ الصورة التقطير (الشكل 3ه) عن طريق تطبيق عتبه قيمه الكثافة إلى شرائح ct ، والذي يتحدد وفقا لكثافة الرسم البياني للشرائح ct باستخدام أسلوب اوتسو35.
ملاحظه: بالنسبة لشرائح CT مع الرسم البياني كثافة رمادية النطاق التي تظهر تداخل كبير من الكثافة بين المرحلة الصلبة ومرحله الفراغ ، والتحقق من صحة صوره التقطير مطلوب باستخدام كتله المرحلة الصلبة36. - افصل الجزيئات الفردية من الشرائح المقطعية المحوسبة باستخدام خوارزميه مستجمعات المياه القائمة علي العلامات وخزن النتائج في صوره مسماه ثلاثية الابعاد (الشكل 3و). التحقق من صحة النتائج عن طريق مقارنه توزيع حجم الجسيمات المحسوبة من الصورة CT إلى تلك من اختبار الغربلة الميكانيكية.
ملاحظه: الوحدة النمطية كائنات منفصلة من البرنامج Avizo النار يمكن استخدامها لتنفيذ هذه الخوارزميه. أزاله الحجارة مساميه من شرائح CT بيناريزيد باستخدام وحده الحدود قتل افيزو النار. للحصول علي نتائج فصل الجسيمات يمكن الاعتماد عليها ، ويقترح القراء لمحاولة مختلفه الجسيمات تجزئه خوارزميات37،38،39.
- تنفيذ استرجاع المرحلة (الشكل 3ب) من إسقاطات CT الخام (الشكل 3ا) من العينة باستخدام البرمجيات الحرة pitre34. إسقاطات التحميل (بما في ذلك الإسقاطات المسطحة والداكنة) في PITRE من القائمة تحميل الصورة. انقر فوق الرمز Ppci. ادخل معلمات المسح الضوئي ذات الصلة وانقر فوق Single لتنفيذ استرداد المرحلة.
- تحليل الصور
- استخراج خصائص الجسيمات من الصورة المسمية. ويستخدم البرنامج النصي matlab لاستخراج خصائص الجسيمات بما في ذلك حجم الجسيمات, منطقه سطح الجسيمات, اتجاه الجسيمات وإحداثيات الجسيمات متمركز.
ملاحظه: يتم استخدام وظائف MATLAB المضمنة الاقليميه، bwprim و pca للحصول علي هذه الخصائص لكل الجسيمات. ويمكن الاطلاع علي وصف أكثر تفصيلا لهذه الإجراءات في عمل تشنغ وانغ28. - استخراج voxels الاتصال من الشرائح CT بواسطة تنفيذ عمليه منطقيه وبين الصورة الثنائية من شرائح ct (الشكل 4) وصوره ثنائيه من خطوط مستجمعات المياه المكتسبة من تنفيذ القائم علي علامة خوارزميه مستجمعات المياه31.
ملاحظه: يمكن ان يحدث الإفراط في الكشف عن voxels الاتصال بسبب تاثير الصوت الجزئي والضوضاء العشوائية من الصور المقطعية40،41. ومع ذلك ، فان الكشف المفرط الطفيف للاتصالات بين الجسيمات لن يكون له اثار كبيره علي الاتجاه العام للسلوك الاتصال بين الجسيمات42.
- استخراج خصائص الجسيمات من الصورة المسمية. ويستخدم البرنامج النصي matlab لاستخراج خصائص الجسيمات بما في ذلك حجم الجسيمات, منطقه سطح الجسيمات, اتجاه الجسيمات وإحداثيات الجسيمات متمركز.
4. CT القائم علي الصورة التنقيب عن السلوك الميكانيكي علي نطاق الحبوب من التربة
ملاحظه: لا ينطبق التحليل التالي المستند إلى الصورة علي الجزيئات الكروية المثالية أو العينات ذات نطاقات التقدير الضيقة جدا (اي عينات أحاديه التشتت). ومع ذلك ، بالنسبة للجسيمات ذات الاستدارة العالية والدرجات الضعيفة (علي سبيل المثال ، 0.3 ~ 0.6 مم الخرز الزجاجي) ، والمنهجية تعطي نتائج جيده (انظر تشنغ وانغ31).
- قياس الجسيمات الجزيئية للعينه. استخدام طريقه تتبع الجسيمات لتتبع الجزيئات الفردية داخل العينة في الفحص المختلفة استنادا إلى حجم الجسيمات أو منطقه سطح الجسيمات. ويرد وصف مفصل لهذه الطريقة في تشنغ وانغ28.
- حساب ترجمه كل الجسيمات خلال اي اثنين من المسح المتتالية. ويحسب الفرق في إحداثيات الجسيمات المركزية بين عمليتي المسح.
- تحديد زاوية دوران كل الجسيمات وفقا للفرق في توجات محورها الرئيسي الرئيسي بين عمليتي المسح.
- قياس حقل السلالة من العينة. استخدم طريقه قائمه علي الشبكة لحساب حقل السلالة خلال اي عمليتي مسح متتاليتين استنادا إلى ترجمه الجسيمات ودوران الجسيمات.
ملاحظه: يتطلب الأسلوب الصور المسمية من العينة من كل من المسح الضوئي ونتائج الجسيمات. يشار إلى القراء بالعمل السابق24 للحصول علي وصف مفصل. - تحليل تطور الاتصال بين الجسيمات من العينة. بناء علي voxels الاتصال المستخرجة ، والصور المسمية من الجزيئات ونتائج تتبع الجسيمات ، وتحليل الاتجاه ناقلات فرع من الاتصالات المفقودة والاتصالات المكتسبة داخل العينة خلال كل زيادة القص.
ملاحظه: يتم إعطاء وصف كامل لهذه الطريقة في تشنغ وانغ31.
Subscription Required. Please recommend JoVE to your librarian.
Representative Results
ويصور الشكل 5 نتائج الجسيمات الجزيئية لعينه من الرمل لايتون بونار (رطل) في شريحة ثنائيه الابعاد خلال فترتين من زيادات القص النموذجية ، الأول والثاني. ويتم تعقب معظم الجزيئات بنجاح ويتم تحديد حجم ترجماتها وتناوبها بعد البروتوكول المذكور أعلاه. واثناء الزيادة الاولي في القص ، لا تظهر عمليات أزاحه الجسيمات ولا تناوب الجسيمات تعريبا واضحا. ومع ذلك ، يتم تطوير النطاق الموضعي في كل من خريطة أزاحه الجسيمات وخريطة دوران الجسيمات اثناء زيادة القص الثانية. ويبين الشكل 6 خرائط السلالة الثمانية والحجمية للعينه خلال الزيادتين في القص. ويلاحظ وجود منطقه تعريب واضحة في خرائط السلالة في زيادة القص الثانية ، مما يدل علي قدره الأسلوب علي تصور فشل الرمل تحت القص ثلاثي المحاور. ويصور الشكل 7 التردد التوجيهي الموحد لمتجات الفروع للاتصالات المكتسبة والاتصالات المفقودة في العينة اثناء زيادات القص الاثنين. تظهر الاتصالات المفقودة تفضيل اتجاهي واضح نحو اتجاه الإجهاد الرئيسي البسيط (اي الاتجاه الأفقي) اثناء كل من زيادات القص.
الشكل 1: الاشعه السينية الصغرى الاعداد CT وجهاز التحميل ثلاثي المحاور. (ا) جهاز ثلاثي المحاور يستخدم بالاقتران مع اعداد الاشعه المقطعية الدقيقة للاشعه السينية. (ب) رؤية موسعه لتركيب الجهاز ثلاثي المحاور اثناء الاختبار ثلاثي المحاور. (ج) جهاز ثلاثي محوري من زاوية مختلفه. وقد تم تعديل هذا الرقم من تشنغ وانغ28. يرجى النقر هنا لعرض نسخه أكبر من هذا الرقم.
الشكل 2: عمليه صنع عينه. (ا) تركيب حجر مساميه وغشاء علي الصفيحة القاعدية ، (ب) تركيب صانع عينه ، (ج) إنشاء شفط داخل صانع العينة ، (د) إسقاط جزيئات الرمل في صانع العينات ، (ه) تركيب الحجر مساميه آخر ولوحه وساده علي راس عينه الرمل ، و (و) أزاله صانع عينه من لوحه قاعده. يرجى النقر هنا لعرض نسخه أكبر من هذا الرقم.
الشكل 3: معالجه صوره الصور المقطعية. (ا) الإسقاط المقطعي الخام ، (ب) الإسقاط المقطعي المحوسب بعد استرجاع المرحلة ، (ج) شريحة ct أفقيه أعيد بناؤها ، (د) الشريحة ct بعد تصفيه الصورة ، (ه) الشريحة المقطعية بعد التقطير بالصورة ، و (و) CT شريحة بعد فصل الجسيمات. يرجى النقر هنا لعرض نسخه أكبر من هذا الرقم.
الشكل 4: توضيح لاستخراج الاتصالات بين الجسيمات من رطل في شرائح 2D. (ا) تنفيذ عمليه منطقيه وبين الصورة الثنائية لشريحة التصوير المقطعي المحوسب والصورة الثنائية لخطوط مستجمعات المياه ، و (ب) الاتصال النموذجي لجزيئي رطل في الفضاء الثلاثي الابعاد (تظهر الجزيئات باللونين الأخضر والأزرق ويكون التلامس يظهر باللون الأحمر). يرجى النقر هنا لعرض نسخه أكبر من هذا الرقم.
الشكل 5: نتائج الجسيمات النمطية النموذجية لعينه رطل خلال زيادات القص اثنين. (ا) الإجهاد-منحني سلاله العينة تحت ضغط ثلاثي المحاور ، (ب) عمليات نزوح الجسيمات وتناوب الجسيمات في العينة اثناء زيادة القص الأول ، و (ج) عمليات نزوح الجسيمات وتناوب الجسيمات في العينة خلال القص الزيادة الثانية. وقد تم تعديل هذا الرقم من تشنغ ووانغ24. يرجى النقر هنا لعرض نسخه أكبر من هذا الرقم.
الشكل 6: حقول السلالة النموذجية من رطل خلال زيادات القص اثنين.
(ا) السلالة الثمانية للقص والسلالة الحجمية للعينه اثناء زيادة القص أولا-(ب) السلالة الثمانية للقص والسلالة الحجمية للعينه اثناء الزيادة الثانية للقص. وقد تم تعديل هذا الرقم من تشنغ ووانغ24. يرجى النقر هنا لعرض نسخه أكبر من هذا الرقم.
الشكل 7: نتائج الاتصال النموذجية بين الجسيمات من رطل خلال زيادات القص اثنين. (ا) التردد التوجيهي الموحد لمتجات الفروع من الاتصالات المكتسبة والاتصالات المفقودة من رطل اثناء زيادة القص أولا. (ب) التردد التوجيهي الموحد لمتجات الفروع من الاتصالات المكتسبة والاتصالات المفقودة من رطل اثناء القص الزيادة الثانية. يرجى النقر هنا لعرض نسخه أكبر من هذا الرقم.
Subscription Required. Please recommend JoVE to your librarian.
Discussion
وقد مكنت تقنيات الاستبانة العالية المكانية للاشعه السينية الدقيقة والمعالجة المتقدمة للصور وتحليلها من اجراء التحقيق التجريبي في السلوك الميكانيكي للتربة الحبيبية تحت القص علي مستويات متعددة النطاقات (اي علي النطاق الكلي والمتوسط مستويات حجم الحبوب). ومع ذلك ، فان التحقيقات التي تستند إلى صور الاشعه المقطعية المتوسطة والحبوب تتطلب اقتناء صور مقطعيه عاليه الاستبانة لعينات التربة اثناء التحميل. ولعل الجانب الأكثر تحديا في هذه العملية هو تلفيق جهاز التحميل ثلاثي المحاور المصغرة التي يمكن استخدامها بالاقتران مع جهاز CT الاشعه السينية الصغرى. وينبغي للمرء ان يجعل النظر الشامل لحجم العينة المطلوبة ، والضغوط التحميل والمعدلات ، بالاضافه إلى القيود المفروضة علي أجهزه الاشعه السينية الصغرى CT مثل الدقة المكانية ، ومنطقه المسح الضوئي وقدره الحمولة من مرحله الدوران.
ويمكن ان يستغرق تحديد الطاقة المثلي للاشعه السينية وزمن التعرض للوقت وقتا طويلا ، ولكنه أمر حاسم للحصول علي صور مقطعيه عاليه الجودة. من المستحسن ان المستخدمين محاولة الطاقات المختلفة وأوقات التعرض خلال المسح الضوئي الأول وتحديد الطاقة المناسبة والتعرض للوقت وفقا لنوعيه الشرائح التي أعيد بناؤها. إلى جانب ذلك ، يمكن ان تنتج العينات مع الخنازير الاوليه المختلفة اثناء اعداد العينة عن طريق إسقاط جزيئات الرمل في قالب العينة من مرتفعات مختلفه. ومع ذلك ، بسبب حجم العينة الصغيرة ، وإنتاج عينه مع المساميه الاوليه محدده أكثر صعوبة بالمقارنة مع الاختبارات ثلاثية المحاور التقليدية. لإنتاج عينه مع المساميه الاوليه التي هي قريبه من قيمه محدده للاختبار ثلاثي المحاور مع المسح المقطعي المحوسب ، ينصح المستخدمين لممارسه إنتاج العينات مقدما.
بالمقارنة مع اختبار ثلاثي المحاور التقليدية ، مصغره في الموقع اختبار ثلاثي المحاور لديه ميزه ان تكون قادره علي استكشاف السلوك الميكانيكي علي نطاق الحبوب من التربة الحبيبية ، بما في ذلك الحبوب كينيماتيك ، وتوطين سلاله والاتصال بين الجسيمات التفاعل ، الخ. حاليا ، طريقه شعبيه بديله للتحقيق في السلوك الميكانيكي علي نطاق الحبوب من التربة الحبيبية هو DEM. علي الرغم من ان هذه التقنية تمكن من نمذجة السلوك الميكانيكي الرمل تحت ظروف التحميل المعقدة ، والاشكال الحبوب ونماذج الاتصال عموما الإفراط في تبسيط لتحقيق كفاءه عاليه في الحوسبة في معظم الدراسات DEM. وفي هذه الحالة ، فان المعلومات المستقية من الحبوب المستخرجة من الرمال الحقيقية باستخدام هذا البروتوكول ضرورية لتحسين التحقق من نماذج الماركات المانيه علي المستويات المتعددة النطاقات. ومن المزايا الأخرى للطريقة المقدمة لحساب السلالة المستندة إلى الصورة المقطعية المحوسبة إدماج دوران الجسيمات في حساب السلالة. وتبين ان طريقه حساب السلالة تنتج نتائج أكثر موثوقيه من طريقه الشبكة القاعدية دون النظر في اثار تناوب الجسيمات24.
وحتى مع مزاياه العديدة ، فان استخدام الاشعه السينية الدقيقة لدراسة تطور الاتصال بين الجسيمات من التربة الحبيبية قد يعاني من الإفراط في الكشف عن الاتصالات بين الجسيمات. وتعتمد دقه نتائج الكشف فيما بين الجسيمات بشده علي الاستبانة المكانية للاشعه السينية الدقيقة. ويرجع ذلك إلى تاثير الحجم الجزئي للاشعه السينية المجهرية ، التي يمكن فيها تحديد جزيئين معزولتين لهما مسافة أصغر من حجم فوكسيل باعتبارهما جزيئين للاتصال. ولحسن الحظ ، تبين ان الاتجاه العام لتطور الاتصال بين الجسيمات داخل التربة الحبيبية لا يتاثر بالكشف المفرط للاتصالات بين الجسيمات. وفي الوقت نفسه ، فان عدم القدرة علي استخراج قوات الاتصال بين الجسيمات داخل التربة الحبيبية هو عيب آخر من الاشعه السينية الصغرى CT مقارنه بالدراسات43،44،45،46،47 و [فوتو-مرن] دراسات48,49. وعلاوة علي ذلك ، وبسبب التحقيق المذكور أعلاه علي نطاق الحبوب المحوسبة المستندة إلى الصورة المطلوبة لتحديد واستخراج الجسيمات الفردية بشكل صحيح من الصور المقطعية ، وتطبيق هذه الطريقة علي التربة مع اشكال الجسيمات غير النظامية للغاية أو عاليه التربة الصليبية التي تحتوي علي فراغات داخل الجسيمات غير النظامية صعبه للغاية.
في المستقبل ، في الموقع اختبار ثلاثي المحاور توفير بيانات وافره علي شكل الحبوب والحبوب والبذور وسوف تسهل دمج الاشكال الجسيمات الحقيقية في النمذجة مارك ألماني. وفي وقت لاحق ، فان النمذجة المحوسبة للصور المقطعية المعتمدة علي الصور ستوفر فهما أفضل للسلوك الميكانيكي علي نطاق الحبوب للتربة الحبيبية قيد التحميل. وفي الوقت نفسه ، ونظرا للقدرة علي استخراج قوات الاتصال بين الجسيمات50، فان الجمع بين انكسار الاشعه السينية والاشعه السينية الدقيقة لاختبار ثلاثي المحاور في الموقع سيكون مفيدا لاستخراج المعلومات الكاملة علي نطاق الحبوب (اي الحبوب [كينميتيكس] وحبه قوه اتصال) من ترب محببه تحت يقص.
Subscription Required. Please recommend JoVE to your librarian.
Disclosures
وليس لدي المؤلفين ما يفصحون عنه.
Acknowledgments
وقد حظيت هذه الدراسة بدعم صندوق البحوث العامة رقم CityU 11213517 من مجلس منح البحوث في منطقه هونغ كونغ الاداريه الخاصة ، ومنحه البحوث رقم 51779213 من المؤسسة الوطنية للعلوم في الصين ، وBL13W beamline من مرفق الإشعاع شنغهاي Synchrotron (S.SF).
Materials
| Name | Company | Catalog Number | Comments |
| Confining pressure offering device | GDS | STDDPC | |
| De-aired water | N/A | N/A | Water de-aired in the lab |
| Leighton Buzzard sand | Artificial Grass Cambridge | Drained Industrial Sand 25 kg | Can be replaced with different soils |
| Miniature triaxial loading device | N/A | N/A | The miniature loading device is specially fabricated by the authors |
| Silicon grease | RS company | RS 494-124 | |
| Synchrotron radiation X-ray micro CT setup | Shanghai Synchrotron Radiation Facility Center (SSRF) | 13W1 | The triaxial testing is carried out at the BL13W beam-line of the SSRF |
| Vacuum pump | Hong Kong Labware Co., ltd. | Rocker 300 |
References
- Cundall, P. A., Strack, O. D. A discrete numerical model for granular assemblies. Géotechnique. 29 (1), 47-65 (1979).
- Rothenburg, L., Bathurst, R. J. Micromechanical features of granular assemblies with planar elliptical particles. Géotechnique. 42 (1), 79-95 (1992).
- Iwashita, K., Oda, M. Rolling resistance at contacts in simulation of shear band development by DEM. Journal of Engineering Mechanics. 124 (3), 285-292 (1998).
- Jiang, M. J., Yu, H. S., Harris, D. A novel discrete model for granular material incorporating rolling resistance. Computers and Geotechnics. 32 (5), 340-357 (2005).
- Ai, J., Chen, J. F., Rotter, J. M., Ooi, J. Y. Assessment of rolling resistance models in discrete element simulations. Powder Technology. 206 (3), 269-282 (2011).
- Zhou, B., Huang, R., Wang, H., Wang, J. DEM investigation of particle anti-rotation effects on the micromechanical response of granular materials. Granular Matter. 15 (3), 315-326 (2013).
- Matsushima, T., Saomoto, H. Discrete element modeling for irregularly-shaped sand grains. Proc. NUMGE2002: Numerical Methods in Geotechnical Engineering. , 239-246 (2002).
- Price, M., Murariu, V., Morrison, G. Sphere clump generation and trajectory comparison for real particles. Proceedings of Discrete Element Modelling. , (2007).
- Ferellec, J., McDowell, G. Modelling realistic shape and particle inertia in DEM. Géotechnique. 60 (3), 227-232 (2010).
- Wiącek, J., Molenda, M., Horabik, J., Ooi, J. Y. Influence of grain shape and intergranular friction on material behavior in uniaxial compression: Experimental and DEM modeling. Powder Technology. , 435-442 (2012).
- Ng, T. T. Fabric study of granular materials after compaction. Journal of Engineering Mechanics. 125 (12), 1390-1394 (1999).
- Cleary, P. W. The effect of particle shape on simple shear flows. Powder Technology. 179 (3), 144-163 (2008).
- Oda, M. Initial fabrics and their relations to mechanical properties of granular material. Soils and Foundations. 12 (1), 17-36 (1972).
- Konagai, K., Tamura, C., Rangelow, P., Matsushima, T. Laser-aided tomography: a tool for visualization of changes in the fabric of granular assemblage. Structural Engineering/Earthquake Engineering. 9 (3), 193-201 (1992).
- Johns, R. A., Steude, J. S., Castanier, L. M., Roberts, P. V. Nondestructive measurements of fracture aperture in crystalline rock cores using X ray computed tomography. Journal of Geophysical Research: Solid Earth. 98 (2), 1889-1900 (1993).
- Ohtani, T., Nakano, T., Nakashima, Y., Muraoka, H. Three-dimensional shape analysis of miarolitic cavities and enclaves in the Kakkonda granite by X-ray computed tomography. Journal of Structural Geology. 23 (11), 1741-1751 (2001).
- Oda, M., Takemura, T., Takahashi, M. Microstructure in shear band observed by microfocus X-ray computed tomography. Géotechnique. 54 (8), 539-542 (2004).
- Fonseca, J., O'Sullivan, C., Coop, M. R., Lee, P. D. Quantifying the evolution of soil fabric during shearing using directional parameters. Géotechnique. 63 (6), 487-499 (2013).
- Fonseca, J., O'Sullivan, C., Coop, M. R., Lee, P. D. Quantifying the evolution of soil fabric during shearing using scalar parameters. Géotechnique. 63 (10), 818-829 (2013).
- Desrues, J., Chambon, R., Mokni, M., Mazerolle, F. Void ratio evolution inside shear bands in triaxial sand specimens studied by computed tomography. Géotechnique. 46 (3), 529-546 (1996).
- Lenoir, N., Bornert, M., Desrues, J., Bésuelle, P., Viggiani, G. Volumetric digital image correlation applied to X-ray microtomography images from triaxial compression tests on argillaceous rock. Strain. 43 (3), 193-205 (2007).
- Higo, Y., Oka, F., Sato, T., Matsushima, Y., Kimoto, S. Investigation of localized deformation in partially saturated sand under triaxial compression using microfocus X-ray CT with digital image correlation. Soils and Foundations. 53 (2), 181-198 (2013).
- Alikarami, R., Andò, E., Gkiousas-Kapnisis, M., Torabi, A., Viggiani, G. Strain localisation and grain breakage in sand under shearing at high mean stress: insights from in situ X-ray tomography. Acta Geotechnica. 10 (1), 15-30 (2015).
- Cheng, Z., Wang, J. Quantification of the strain field of sands based on X-ray micro-tomography: A comparison between a grid-based method and a mesh-based method. Powder Technology. , 314-334 (2019).
- Hall, S. A., Bornert, M., Desrues, J., Pannier, Y., Lenoir, N., Viggiani, G., Bésuelle, P. Discrete and continuum analysis of localised deformation in sand using X-ray μCT and volumetric digital image correlation. Géotechnique. 60 (5), 315-322 (2010).
- Andò, E., Hall, S. A., Viggiani, G., Desrues, J., Bésuelle, P. Grain-scale experimental investigation of localised deformation in sand: a discrete particle tracking approach. Acta Geotechnica. 7 (1), 1-13 (2012).
- Watanabe, Y., Lenoir, N., Otani, J., Nakai, T. Displacement in sand under triaxial compression by tracking soil particles on X-ray CT data. Soils and Foundations. 52 (2), 312-320 (2012).
- Cheng, Z., Wang, J. A particle-tracking method for experimental investigation of kinematics of sand particles under triaxial compression. Powder Technology. 328, 436-451 (2018).
- Matsushima, T., Katagiri, J., Uesugi, K., Nakano, T., Tsuchiyama, A. Micro X-ray CT at Spring-8 for granular mechanics. Soil Stress-Strain Behavior: Measurement, Modeling and Analysis. A Collection of Papers of the Geotechnical Symposium in Rome. 146, 225-234 (2006).
- Andò, E., Viggiani, G., Hall, S. A., Desrues, J. Experimental micro-mechanics of granular media studied by X-ray tomography: recent results and challenges. Géotechnique Letters. 3, July-September 142-146 (2013).
- Cheng, Z., Wang, J. Experimental investigation of inter-particle contact evolution of sheared granular materials using X-ray micro-tomography. Soils and Foundations. 58 (6), 1492-1510 (2018).
- Karatza, Z., Andò, E., Papanicolopulos, S. A., Ooi, J. Y., Viggiani, G. Evolution of deformation and breakage in sand studied using X-ray tomography. Géotechnique. , 1-11 (2017).
- Cheng, Z., Wang, J. F., Coop, M. R., Ye, G. L. A miniature triaxial apparatus for investigating the micromechanics of granular soils with in-situ X-ray micro-tomography scanning. Frontiers of Structural and Civil Engineering. , (2019).
- Chen, R. C., Dreossi, D., Mancini, L., Menk, R., Rigon, L., Xiao, T. Q., Longo, R. PITRE: software for phase-sensitive X-ray image processing and tomography reconstruction. Journal of Synchrotron Radiation. 19 (5), 836-845 (2012).
- Otsu, N. A threshold selection method from gray-level histograms. IEEE Trans. Systems, Man Cybernet. 9 (1), 62-66 (1979).
- Karatza, Z. Study of temporal and spatial evolution of deformation and breakage of dry granular materials using X-ray computed tomography and the discrete element method. , University of Edinburgh. PhD Thesis (2017).
- Shi, Y., Yan, W. M. Segmentation of irregular porous particles of various sizes from X-ray microfocus computer tomography images using a novel adaptive watershed approach. Géotechnique Letters. 5 (4), 299-305 (2015).
- Zheng, J., Hryciw, R. D. Segmentation of contacting soil particles in images by modified watershed analysis. Computers and Geotechnics. 73, 142-152 (2016).
- Lai, Z., Chen, Q. Reconstructing granular particles from X-ray computed tomography using the TWS machine learning tool and the level set method. Acta Geotechnica. 14 (1), 1-18 (2019).
- Wiebicke, M., Andò, E., Herle, I., Viggiani, G. On the metrology of interparticle contacts in sand from x-ray tomography images. Measurement Science and Technology. 28 (12), 1-14 (2017).
- Karatza, Z., Andò, E., Papanicolopulos, S. A., Viggiani, G., Ooi, J. Y. Effect of particle morphology and contacts on particle breakage in a granular assembly studied using X-ray tomography. Granular Matter. 21 (3), 1-13 (2019).
- Cheng, Z. Investigation of the grain-scale mechanical behavior of granular soils under shear using X-ray micro-tomography. , City University of Hong Kong. PhD Thesis (2018).
- Antony, S. J. Evolution of force distribution in three-dimensional granular media. Physical Review E. , (2000).
- Kruyt, N. P., Rothenburg, L. Probability density functions of contact forces for cohesionless frictional granular materials. International Journal of Solids and Structures. 39 (3), 571-583 (2002).
- Marketos, G., Bolton, M. D. Quantifying the extent of crushing in granular materials: a probability-based predictive method. Journal of the Mechanics and Physics of Solids. 55 (10), 2142-2156 (2007).
- Cheng, Z., Wang, J. Quantification of particle crushing in consideration of grading evolution of granular soils in biaxial shearing: A probability-based model. International Journal for Numerical and Analytical Methods in Geomechanics. 42 (3), 488-515 (2018).
- Zhou, B., Wang, J., Wang, H. A new probabilistic approach for predicting particle crushing in one-dimensional compression of granular soil. Soils and Foundations. 54 (4), 833-844 (2014).
- Geng, J., Reydellet, G., Clément, E., Behringer, R. P. Green's function measurements of force transmission in 2D granular materials. Physica D: Nonlinear Phenomena. 182, 274-303 (2003).
- Majmudar, T. S., Behringer, R. P. Contact force measurements and stress-induced anisotropy in granular materials. Nature. 435, 1079 (2005).
- Hurley, R. C., Hall, S. A., Andrade, J. E., Wright, J. Quantifying interparticle forces and heterogeneity in 3D granular materials. Physical Review Letters. 117 (9), 098005 (2016).
