التطبيقات
التنبؤ بقوة سقالات β-TCP
يربط النمذجة معلمات البنية المسامية بالاستجابة الميكانيكية للسقالة الحيوية؛ وهي أداة بحث وليست استنتاجًا طبيًا.
أُعدت هذه المذكرة التحريرية من KUNWU استنادًا إلى المصدر التقني المقدم. وتخص النتائج الدراسة أو العينة الموصوفة وليست مواصفات لمعدات KUNWU.
لماذا يصعب التنبؤ بمقاومة السقالات؟
β-TCP (β-tricalcium phosphate) عبارة عن سيراميك حيوي قابل للتحلل وهو قريب من معادن العظام البشرية ويستخدم على نطاق واسع في سقالات إصلاح العظام وزراعة الأسنان وإعادة بناء الجمجمة. يمكن للطباعة ثلاثية الأبعاد DLP (معالجة الضوء الرقمية) أن تشكل بنية مسام مترابطة عالية الدقة في قطعة واحدة، والتي تناسب تمامًا طوبولوجيا المسام المطلوبة لنمو الأنسجة العظمية. هذا المزيج من "مرونة التصميم + الدقة العالية" يجعل DLP أكثر تنافسية على مسار سقالة العظام.
ومع ذلك، بمجرد دخول الدعامة إلى مجال إصلاح العظام الحاملة، فإن أول ما يواجهه المهندسون ليس "هل يمكن إصلاحها؟" ولكن "هل يمكن حملها بعد إصلاحها؟" β-TCP نفسه هش للغاية، ويختلف توزيع الضغط الناتج عن هياكل TPMS المختلفة (السطح الأدنى ثلاثي الفترات) اختلافًا كبيرًا. بالإضافة إلى ذلك، فإن عملية DLP نفسها تجلب اضطرابات "غير تصميمية" مثل عيوب الطبقات البينية، والانكماش غير المتساوي، والثقوب الداخلية، ويصبح من الصعب للغاية التنبؤ بالكمية بقوة الضغط. الطريقة التقليدية هي "الضرب على قطعة واحدة، والضغط على قطعة واحدة، والضبط مرة أخرى". يمكن أن تستغرق الجولة الواحدة بضعة أشهر أو قد تستغرق سنة كاملة.
المشكلة ليست فقط "بطيئة". تتكون عائلة TPMS نفسها من نوعين: الورقة (النوع ذو الجدران الرقيقة) والهيكل العظمي (النوع الهيكلي). وهي مقسمة إلى ثمانية طبولوجيا رئيسية مثل Diamond وGyroid وI-WP وFischer-KochS وPrimitive وSplit-P وNeovius وPMY. كل نوع له تغيرات في المسامية ناتجة عن سمك الجدار القابل للتعديل. مساحة التصميم بأكملها ضخمة، ويجب اجتياز كل منها عن طريق الطباعة + اختبارات الضغط، وهو أمر غير واقعي ومكلف.
ما يريد فريق Shuo Yin فعله واضح جدًا: بدلاً من الطباعة مرارًا وتكرارًا، هل يمكننا أولاً استخدام الإحصائيات + العناصر المحدودة "لحساب" القوة؛ هل يمكن استقراءه حتى عندما تتغير البنية والمسامية؟ وهذا هو أحدث أعمالهم المنشورة في المجلة الدولية للعلوم الميكانيكية.
إحصاء وايبل والعناصر المحدودة والضرر التراكمي
قام الفريق بتقسيم المشكلة إلى ثلاثة أجزاء واستخدم ثلاثة أنواع من الأدوات لحلها.
القطعة الأولى هي القوة المادية الجوهرية. من الطبيعي أن يؤدي الكسر الهش لـ β-TCP إلى تشتت البيانات. للحديث عن "القوة"، يجب علينا أولا أن يكون لدينا توزيع إحصائي. استخدم الفريق عملية DLP (Lithoz CeraFab 7500، والملاط التجاري LithaBone TCP 300، بسماكة طبقة 25 ميكرومتر، وتعريض 168 مللي جول/سم²، وتلبد عند 1200 درجة مئوية) لتحضير 86 كتلة كثيفة بحجم 1.2×1.2×2.4 ملم مكعب. معامل Weibull المقاس m = 2.34 وقوة الضغط المميزة σ C = 365 MPa؛ معامل يونغ يتوافق مع m=6.74، EC=10.8 GPa. توفر هاتان المجموعتان من الأرقام نقطة ربط المعايرة "لقوة المادة المحلية" لجميع التنبؤات المعمارية اللاحقة.
القطعة الثانية هي تأثير تضخيم الإجهاد الهيكلي. نفس المادة الموضوعة في وحدات TPMS المختلفة سيكون لها توزيعات ضغط مختلفة بشكل كبير. أجرى الفريق عمليات محاكاة ضغط أحادي المحور للعناصر المحدودة على 14 نوعًا من وحدات TPMS مع مسامية مستهدفة تبلغ 50% (باستثناء Neovius وPMY، اللذين يحتفظان فقط بالنوع الهيكلي، والأنواع الستة الأخرى من الهيكل تظهر جميعها في أزواج من الصفائح/الهيكل العظمي، بإجمالي 14 نوعًا)، واستخرجوا الحد الأقصى لتوزيع حجم الضغط الرئيسي (توزيع حجم SCF) لكل بنية. وهذا يعادل قياس "تركيزات الإجهاد المستحثة هندسيًا" في بصمة قابلة للحساب والمقارنة.
الكتلة الثالثة هي معيار الفشل. قارن الفريق معايير الكسور الهشة المختلفة مثل الحد الأقصى للإجهاد الأساسي، وإجهاد دراكر-براجر، وإجهاد ميزس، وإجهاد القص، وتوصل أخيرًا إلى أن الحد الأقصى للإجهاد الرئيسي هو معيار الفشل الأكثر دقة وثباتًا، وهو ما يتوافق مع الصورة المادية للكسر الجانبي الناتج عن الشد الخزفي. إلى جانب معيار الضرر التراكمي لـ "جزء حجم الفشل الحرج" - عندما يصل حجم الفشل في الهيكل إلى قيمة حرجة، يتم أخذ الفشل الإجمالي في الاعتبار - يمكن عكس "البصمة الهندسية" إلى "تنبؤ بالقوة الكلية".
قُدّر عدم اليقين أيضًا بأخذ عينات عشوائية من توزيعات معامل تركيز الإجهاد والحجم SCF–Volume للهندسات الأربع عشرة، وإرفاق كل تنبؤ بفاصل تنبؤ 95%. بذلك لا يعطي النموذج رقم مقاومة واحدًا فقط، بل نطاقًا يساعد المصمم على تقييم الحالات الأقل ملاءمة؛ ولا يمثل الفاصل ضمانًا مطلقًا لأسوأ حالة.
معايرة واحدة لـ14 هندسة TPMS ومسامية متغيرة
إطار التنبّؤ هذا "مميّز" و "مجاز" في خطوتين.
في الخطوة الأولى، يتم استخدام 14 هيكل TPMS مع مسامية 50% للمعايرة: يتم إجراء 7 عينات ضغط قياسية بحجم 4.0×4.0×8.0 مم مكعب لكل هيكل، ويتم إجراء ضغط أحادي المحور لقياس القوة المقاسة؛ وفي الوقت نفسه، يتم تشغيل العناصر المحدودة للحصول على القوة المتوقعة. تظهر مقارنة التنبؤات مقابل القياسات أن معيار الضغط الرئيسي الأقصى أفضل بكثير من معايير الفشل الأخرى من حيث اتساق التصنيف.
والخطوة الثانية هي التحقق على مستوى المسامية الذي لا يشارك في المعايرة (تم تحديد طبولوجيا ورقتين، D وG، كل منهما بمستويين 40% و60%، وهما D40/D60/G40/G60). وبدون إعادة ضبط المعلمات، تظل النتائج المتوقعة تقع ضمن فترة 95% من القياس الفعلي. وهذا يعني أنه يمكن استقراء هذا الإطار في كلا البعدين الهيكلي + المسامية، وليست هناك حاجة لإجراء تجربة عينة كبيرة في كل مرة يتم فيها تعديل التصميم.
| هندسة TPMS بمسامية 50% | مقاومة الضغط (MPa) | معامل يونغ (GPa) |
|---|---|---|
| I50-sh | 57.5 ± 3.5 | 3.8 ± 0.2 |
| D50-sh | 53.1 ± 6.5 | 3.8 ± 0.2 |
| P50-sk | 30.0 ± 3.8 | 2.8 ± 0.2 |
| G50-sk | 29.2 ± 5.3 | 3.0 ± 0.2 |
| F50-sh | 27.3 ± 1.9 | 3.5 ± 0.1 |
| N50-sk | 26.0 ± 10.7 | 3.0 ± 0.4 |
| S50-sh | 22.0 ± 3.7 | 3.2 ± 0.3 |
| G50-sh | 20.7 ± 2.3 | 3.2 ± 0.3 |
| S50-sk | 17.9 ± 0.9 | 2.8 ± 0.2 |
| PMY50-sk | 16.3 ± 4.7 | 2.5 ± 0.3 |
| P50-sh | 11.6 ± 2.0 | 2.7 ± 0.6 |
| I50-sk / D50-sk / F50-sk | ≈ 11.2–11.4 | 2.4–2.8 |
الجدول 1 · قوة الضغط ومعامل يونغ لـ 14 هيكل TPMS عند مسامية 50% (sh=sheet، sk=skeletal؛ الهياكل الثلاثة في الصف الأخير لها قيم متشابهة وتم إدراجها معًا)
بمجرد النظر إلى الجدول، يمكنك العثور على ظاهرة غير بديهية: عند نفس المسامية البالغة 50%، تتجاوز قوة صفائح I-WP والصفائح الماسية 50 ميجا باسكال، ويصل معامل يونج إلى 3.8 جيجا باسكال، وهو ضعف ارتفاع معظم الهياكل الهيكلية. بالنسبة لنفس الطوبولوجيا، قد تنخفض كثافة الطبقة إلى الهيكل العظمي مباشرة من 50 ميجا باسكال إلى 11 ميجا باسكال - إذا تم تغيير الشكل الهندسي، ستكون الكثافة مختلفة 5 مرات. من الصعب على المصممين الحكم مسبقًا إذا كانوا يعتمدون فقط على الخبرة. وهذا يسلط الضوء فقط على قيمة تنبؤات نوع "البصمة الهندسية".
تتبع العيوب: كيف تؤثر الهندسة في الشقوق؟
توقعات الهندسة المثالية ليست مثالية. وأظهرت إعادة البناء بالأشعة المقطعية الدقيقة أن الانحرافات المتبقية تأتي بشكل رئيسي من العيوب العيانية التي أدخلتها العملية، وأن هذه العيوب ليست موزعة بشكل عشوائي - حيث "يختار" الهيكل أنواع الشقوق التي من المرجح أن تتطور.
وقد لخص الفريق هذه الظاهرة في وصفين معماريين:
① نسبة الكثافة: الجزء الحجمي لوحدات البكسل الكثيفة تمامًا في "الكثافة المتوقعة" في اتجاه الطباعة. كلما زادت نسبة الكثافة، كانت الاستمرارية الرأسية أفضل وكلما زادت قدرة الحمل؛ ولكن في الوقت نفسه، يكون تركيز إجهاد الواجهة أكثر وضوحًا أثناء إزالة الشحوم والتلبيد، ويزداد خطر إزالة الترابط بين الطبقات. نسب الكثافة P50-sk وG50-sk مرتفعة على حد سواء، وشوهدت بالفعل شقوق بينية على الأشعة المقطعية الدقيقة؛ على الرغم من أن نسب كثافة I50-sk وD50-sh وPMY50-sk ليست منخفضة، إلا أن مناطق الكثافة أكثر تشتتًا ولا تشكل "أعمدة داخل الطبقة" مستمرة، لذلك يتحول شكل العيب إلى مسام داخلية.
② نسبة التعبئة المستوية: الحد الأقصى لجزء الطور الصلب داخل المستوى في جميع طبقات الطباعة. القيمة العالية تعني أن حجم الطبقة المفردة يتقلص بشكل أكثر عنفًا أثناء إزالة الشحوم، ويمكن أن يتراكم الضغط الداخلي الناجم عن الانكماش غير المتساوي بسهولة في الطبقة ويتطور بشكل أكبر إلى ثقوب داخلية.
إن استخدام هذين الخيارين "كمحددات للعيوب" يمكن أن يفسر سبب اختلاف أوضاع الفشل التي تظهرها TPMS المختلفة بشكل كبير. استخدم الفريق أيضًا P50-sk لإجراء التحقق العكسي: تم تغذية الهندسة "المحتوية على العيوب" التي أعيد بناؤها بواسطة التصوير المقطعي المحوسب مباشرة إلى العنصر المحدود، وتم تحسين دقة التنبؤ بشكل كبير. وهذا بدوره يثبت أن السبب وراء انحياز تنبؤات الهندسة المثالية لا يرجع في المقام الأول إلى مشكلة النموذج نفسه، بل "الجروح" التي خلفتها الحرفية المخفية في الهندسة.
الدلالات لترميم العظام
وبوضع هذا الإطار في المشهد الفعلي لإصلاح العظام الحاملة، هناك ثلاث نقاط قيمة:
أولاً، قم بفحص البنية بسرعة. في مرحلة تصميم الطباعة ثلاثية الأبعاد للسيراميك DLP، لم يعد المهندسون بحاجة إلى إجراء جولة من الطباعة + الضغط لكل طوبولوجيا TPMS + مجموعة سمك الجدار؛ بدلاً من ذلك، يستخدمون إطار عمل للتنبؤ للقيام أولاً بجولة من "فرز القوة" على الكمبيوتر لتضييق البنى المرشحة إلى 2-3 قبل تشغيلها على الكمبيوتر. يمكن تقصير دورة التطوير بشكل كبير.
ثانيا، يتم إعطاء تقرير مع عدم اليقين. ويعني الفاصل الزمني للتنبؤ بنسبة 95% أن المصممين يمكنهم الحصول على "تقديرات متفائلة" و"تقديرات متحفظة" في نفس الوقت، والحصول على قيمة مرجعية مباشرة لمعايرة أسوأ الحالات للأجزاء الحاملة - وهو أمر لم تتمكن تقارير "قوة النقطة الواحدة" السابقة من توفيره.
ثالثًا، دمج قوانين العيوب في قواعد التصميم. يعمل واصفا نسبة الكثافة ونسبة ملء المستوى على تحويل "أي معمارية يسهل فكها وأيها يسهل الحصول على ثقوب داخلية" إلى مؤشر هندسي قابل للقياس الكمي. يمكن للمصممين بدورهم تقييد اختيار البنية بناءً على متطلبات الحاملة. على سبيل المثال، يتم التحكم في منطقة نسبة الكثافة العالية تحت عتبة معينة لتجنب خطر حدوث تشققات بين الطبقات.
وبطبيعة الحال، يعترف الفريق أيضًا بالقيود المفروضة على الإطار الحالي: حيث يتم استخدام 7 عينات فقط للتحقق لكل بنية، وحجم العينة الإحصائية محدود؛ لم يتم بعد دمج التباين الناتج عن اتجاه البناء في النموذج. هذه كلها اتجاهات يمكن استكمالها في المستقبل. كشف فريق Yin عن جميع البيانات الخاصة بـ Zenodo، ويمكن للقراء المهتمين إعادة استخدامها مباشرة لإجراء مقارنة مرجعية.
بالنسبة للزملاء المشاركين في الطباعة ثلاثية الأبعاد للسيراميك DLP والبحث والتطوير في سقالات العظام، فإن هذه الطريقة سهلة الحمل للغاية: تغيير مادة السيراميك ومجموعة من بنية TPMS، تظل العملية الأساسية دون تغيير. فهو لا يحل مشكلة "التنبؤ بقوة β-TCP" فحسب، بل يدفع أيضًا مسألة "الموثوقية الميكانيكية للأجزاء المسامية المطبوعة ثلاثية الأبعاد من السيراميك" إلى موضع يمكن من خلاله إصدار حكم على المستوى الهندسي.
المراجع
Zhao Y، Liu K، Lupoi R، Yin S. التنبؤ بالقوة مع تحليل العيوب في سقالات β-TCP المطبوعة بواسطة DLP. المجلة الدولية للعلوم الميكانيكية، 2026، 326: 111873.