2025-11-19T11:34:14.421060

Isogeometric Topology Optimization Based on Topological Derivatives

Teixeira, Krenn, Gangl et al.
Topology optimization is a valuable tool in engineering, facilitating the design of optimized structures. However, topological changes often require a remeshing step, which can become challenging. In this work, we propose an isogeometric approach to topology optimization driven by topological derivatives. The combination of a level-set method together with an immersed isogeometric framework allows seamless geometry updates without the necessity of remeshing. At the same time, topological derivatives provide topological modifications without the need to define initial holes [7]. We investigate the influence of higher-degree basis functions in both the level-set representation and the approximation of the solution. Two numerical examples demonstrate the proposed approach, showing that employing higher-degree basis functions for approximating the solution improves accuracy, while linear basis functions remain sufficient for the level-set function representation.
academic

تحسين الطوبولوجيا المتساوي الهندسي بناءً على المشتقات الطوبولوجية

المعلومات الأساسية

  • معرّف الورقة: 2509.09236
  • العنوان: Isogeometric Topology Optimization Based on Topological Derivatives
  • المؤلفون: Guilherme Henrique TEIXEIRA (جامعة غراتس للتكنولوجيا)، Nepomuk KRENN (معهد يوهان رادون)، Peter GANGL (معهد يوهان رادون)، Benjamin MARUSSIG (جامعة غراتس للتكنولوجيا)
  • التصنيف: math.NA cs.CE cs.NA math.OC
  • تاريخ النشر: 15 أكتوبر 2025 (الإصدار الثالث من arXiv)
  • رابط الورقة: https://arxiv.org/abs/2509.09236

الملخص

يعتبر تحسين الطوبولوجيا أداة مهمة في الهندسة لتصميم الهياكل المحسّنة. ومع ذلك، فإن التغييرات الطوبولوجية عادة ما تتطلب إعادة تشبيك، وهو ما قد يصبح تحديًا كبيرًا. تقترح هذه الورقة طريقة تحسين طوبولوجيا متساوية الهندسة مدفوعة بالمشتقات الطوبولوجية. يسمح الجمع بين طريقة مجموعة المستويات والإطار الهندسي المتساوي المغمور بتحديثات هندسية سلسة دون الحاجة إلى إعادة تشبيك. في الوقت ذاته، توفر المشتقات الطوبولوجية تعديلات طوبولوجية دون الحاجة إلى تحديد الثقوب الأولية. تم دراسة تأثير الدوال الأساسية من الرتب العليا في تمثيل مجموعة المستويات وتقريب الحل. يتحقق مثالان عدديان من الطريقة المقترحة، مما يدل على أن تقريب الحل باستخدام دوال أساسية من رتب عليا يحسّن الدقة، بينما تظل الدوال الأساسية الخطية كافية لتمثيل دالة مجموعة المستويات.

الخلفية البحثية والدافع

تعريف المشكلة

تتمحور المشكلة الأساسية التي يسعى هذا البحث إلى حلها حول صعوبة إعادة بناء الشبكة في تحسين الطوبولوجيا التقليدي. أثناء عملية تحسين الهيكل، التغييرات الطوبولوجية في مجال التصميم (مثل ظهور الثقوب أو دمجها أو اختفاؤها) عادة ما تتطلب إعادة إنشاء شبكة حسابية، وهو ما لا يرفع تكاليف الحساب فحسب، بل قد يؤدي أيضًا إلى عدم استقرار عددي.

الأهمية

يعتبر تحسين الطوبولوجيا أداة مهمة في تصميم الهياكل، وتطبيقاته واسعة في مجالات الفضاء والملاحة والتصنيع السيارات والهندسة المعمارية. يحمل حل مشكلة إعادة بناء الشبكة أهمية كبيرة لتحسين كفاءة التحسين وضمان الاستقرار العددي.

قيود الطرق الموجودة

  1. طريقة مجموعة المستويات التقليدية: بينما يمكنها تجنب إعادة بناء الشبكة الصريحة، إلا أنها تواجه صعوبة في إنشاء ثقوب جديدة تلقائيًا، وتعتمد على تكوين الثقوب في الهندسة الأولية
  2. طريقة الكثافة: بينما تتجنب مشاكل الشبكة، قد تنتج عناصر رمادية وأنماط رقعة الشطرنج
  3. طريقة المشتقات الشكلية: يمكنها التعامل فقط مع التغييرات الحدودية، وليس التغييرات الطوبولوجية

الدافع البحثي

يسعى المؤلفون إلى دمج التحليل الهندسي المتساوي (IGA) وطريقة مجموعة المستويات والمشتقات الطوبولوجية، بهدف:

  • القضاء على متطلبات إعادة بناء الشبكة
  • إنشاء ثقوب جديدة تلقائيًا دون الحاجة إلى تعريفها مسبقًا
  • الاستفادة من النعومة من الرتب العليا للمنحنيات B-spline لتحسين الدقة
  • تبسيط عملية التحديث الهندسي

المساهمات الأساسية

  1. إطار عمل مبتكر: دمج طريقة المشتقات الطوبولوجية البحتة لـ Amstutz و Andrä مع التحليل الهندسي المتساوي للمرة الأولى
  2. تحسين بدون إعادة بناء شبكة: تحقيق تحديثات هندسية سلسة أثناء التغييرات الطوبولوجية من خلال الإطار الهندسي المتساوي المغمور
  3. دراسة الدوال الأساسية من الرتب العليا: دراسة منهجية لتأثير رتب متعددة الحدود المختلفة على تمثيل مجموعة المستويات وتقريب الحل
  4. التحقق العددي: التحقق من فعالية الطريقة من خلال أمثلة كلاسيكية، مما يثبت أن الدوال الأساسية الخطية كافية لتمثيل مجموعة المستويات بينما تحسّن الدوال الأساسية من الرتب العليا دقة الحل

شرح التفاصيل الطريقة

تعريف المهمة

الإدخال: مجال التصميم D، الشروط الحدودية، الأحمال، خصائص المواد الإخراج: توزيع المواد الأمثل Ω ⊂ D الهدف: تقليل المرونة مع الامتثال لقيود الحجم القيود: معادلات المرونة الخطية، الشروط الحدودية

معمارية النموذج

1. صيغة مشكلة المرونة الخطية

حل على مجال ثابت D:

-∇ · (αΩσ(u)) = 0 في D
u = 0 على ΓD
σ(u) · n = τ على ΓN
σ(u) · n = 0 على Γ0

حيث يتم تعريف معامل المادة αΩ كما يلي:

αΩ = {
  αin   في Ω (منطقة المادة)
  αout  في D\Ω (منطقة الفراغ، αout ≪ 1)
}

2. تمثيل مجموعة المستويات

يتم تمثيل المجال Ω بشكل ضمني من خلال دالة مجموعة المستويات ϕ:

ϕ(x) < 0 ⟺ x ∈ Ω
ϕ(x) = 0 ⟺ x ∈ ∂Ω  
ϕ(x) > 0 ⟺ x ∈ D\Ω

يتم تقسيم دالة مجموعة المستويات باستخدام منحنيات B-spline:

ϕ(ξ,η) = Σᵢⱼ Bᵢ,d(ξ)Bⱼ,d(η)cᵢⱼ

3. حساب المشتقات الطوبولوجية

بالنسبة لمشكلة المرونة الخطية، يكون التعبير التحليلي للمشتقة الطوبولوجية:

dℐ(Ω)(x₀) = {
  -3αin((αout-αin)/(2αout+αin))σ(u):ε(u) - l  إذا كان x₀ ∈ Ω
  -3αout((αin-αout)/(2αin+αout))σ(u):ε(u) + l  إذا كان x₀ ∈ D\Ω
}

4. خوارزمية تحديث مجموعة المستويات

استخدام الاستيفاء الخطي الكروي لتحديث مجموعة المستويات:

ϕⁱ⁺¹ = (1/sin θᵢ)[sin((1-κᵢ)θᵢ)ϕⁱ + sin(κᵢθᵢ)gΩᵢ]

حيث θᵢ هي الزاوية L² بين مجموعة المستويات الحالية والمشتقة الطوبولوجية المعممة gΩᵢ.

نقاط الابتكار التقني

1. الطريقة الهندسية المتساوية المغمورة

  • استخدام شبكة خلفية B-spline ثابتة، تجنب إعادة بناء الشبكة
  • التمييز بين مناطق المادة والفراغ من خلال معامل المادة αΩ
  • الحفاظ على النعومة من الرتب العليا للتمثيل الهندسي

2. المشتقات الطوبولوجية البحتة المدفوعة

  • عدم حل معادلة Hamilton-Jacobi، الاعتماد فقط على المشتقات الطوبولوجية
  • إنشاء ثقوب جديدة تلقائيًا دون الحاجة إلى تعريفها مسبقًا
  • تبسيط عملية التحسين وتحسين الكفاءة الحسابية

3. معالجة العناصر المقطوعة

  • استخدام مكتبة التكامل الهندسي الضمني لحساب نسب القطع بدقة
  • استخدام التكامل الغاوسي القياسي لتجنب عدم الاستقرار في التكامل على المناطق الصغيرة
  • تطبيق التصفية المتوسطة لتنعيم قفزات خصائص المادة

إعداد التجارب

الأمثلة العددية

تم اختيار مثالين معياريين كلاسيكيين لتحسين الطوبولوجيا:

  1. مشكلة الكابولي: مجال مستطيل 2×1، مثبت من الطرف الأيسر، حمل نقطي على الطرف الأيمن
  2. مشكلة الربع حلقة: ربع حلقة بنصف قطر داخلي 1 ونصف قطر خارجي 2، مثبت من الأسفل، حمل نقطي من الأعلى

إعدادات المعاملات

  • الشبكة: شبكة خلفية 128×128
  • معاملات المادة: αin = 1, αout = 10⁻⁴
  • التحكم في الحجم: l = 5
  • معامل التصفية: γ = 10⁻⁴
  • خصائص المادة: E = 1, ν = 1/3
  • شروط التقارب: εθ = 1، الحد الأقصى للتكرارات 200

إعدادات المقارنة

دراسة تكوينين:

  1. نفس الرتبة: مجموعة المستويات والحل يستخدمان نفس رتبة متعددة الحدود (d = p)
  2. رتب مختلطة: مجموعة مستويات خطية (d = 1) مع حل من رتب عليا (p > 1)

مؤشرات التقييم

  • قيمة دالة الهدف J (المرونة)
  • زاوية التقارب θ
  • المساحة النسبية Ai/A₀
  • عدد تكرارات التقارب

نتائج التجارب

النتائج الرئيسية لمشكلة الكابولي

أداء التقارب

  • ميزة الرتب العليا: p ≥ 2 يتقارب أسرع من p = 1، مع انخفاض أسرع لدالة الهدف في المراحل الوسيطة
  • التقارب المتسق: جميع التكوينات تتقارب إلى حل أمثل متشابه (J ≈ 6.97-7.00)
  • كفاءة التكرار: p = 2 و p = 4 يتطلبان نفس عدد التكرارات، بينما p = 3 يتطلب 29 خطوة مع مجموعة مستويات من رتبة عليا مقابل 24 خطوة مع مجموعة مستويات خطية

تأثير رتبة مجموعة المستويات

  • مجموعة المستويات الخطية (d = 1) تنتج أشكالًا نهائية مشابهة لمجموعات المستويات من الرتب العليا (d = p)
  • الفروقات بين الرتب الفردية والزوجية: الرتب الفردية تشارك نقاط Greville مما يسبب تأثير توسيط، ينتج عنه اختلافات طفيفة

النتائج الرئيسية لمشكلة الربع حلقة

ظاهرة الحد الأدنى المحلي

  • اختلافات التقارب: التكوينات من الرتب المنخفضة p = 1,2 تتوقف عند زوايا أعلى (θ = 4.72, 9.33)
  • الحل الأمثل: p = 3(d = 1) يحقق أقل قيمة لدالة الهدف 4.12، بينما p = 3(d = 3) يحقق 4.13
  • توازن الكفاءة: p = 2(d = 2) يحقق 4.15 في 129 تكرار، بكفاءة أعلى

تنوع الأشكال

تنتج التكوينات المختلفة هياكل طوبولوجية مختلفة، مما يعكس التحديات الكامنة للحد الأدنى المحلي في تحسين الطوبولوجيا.

النتائج الرئيسية

  1. فائدة تقريب الحل من الرتب العليا: تحسين سرعة التقارب ودقة الحل
  2. كفاية مجموعة المستويات الخطية: d = 1 كافية لتمثيل الحدود الهندسية بدقة
  3. عدم الاعتماد على الشبكة: مع شبكة دقيقة 128×128، ميزة الرتب العليا محدودة
  4. أهمية التصفية: التصفية المتوسطة وتصفية PDE حاسمة لجودة النتائج

الأعمال ذات الصلة

تصنيف طرق تحسين الطوبولوجيا

  1. طرق الكثافة: SIMP و RAMP وغيرها، قد تنتج عناصر رمادية
  2. طرق مجموعة المستويات: تتجنب العناصر الرمادية لكن تتطلب ثقوبًا أولية
  3. طرق المشتقات الطوبولوجية: يمكنها إنشاء ثقوب تلقائيًا، وهي الطريقة المستخدمة في هذه الورقة

تحسين الطوبولوجيا الهندسي المتساوي

  • الأعمال المبكرة: نموذج المرحلة لـ Dedè وآخرون (2012)، الدوال الأساسية الشعاعية لـ Shojaee وآخرون (2012)
  • طرق B-spline: Jahangiry و Tavakkoli (2017)، Khatibinia وآخرون (2020)
  • ابتكار هذه الورقة: دمج المشتقات الطوبولوجية البحتة مع الطريقة الهندسية المتساوية للمرة الأولى

الطرق المغمورة

تطبيقات واسعة في ديناميكا الموائع والميكانيكا الصلبة وغيرها، وتقدم هذه الورقة تطبيقها في تحسين الطوبولوجيا لتحقيق تحسين بدون إعادة بناء شبكة.

الخلاصة والنقاش

الاستنتاجات الرئيسية

  1. تم تطوير إطار عمل تحسين طوبولوجيا هندسي متساوي مغمور مدفوع بالمشتقات الطوبولوجية بنجاح
  2. الدوال الأساسية من الرتب العليا تحسّن تقريب الحل لكن الدوال الأساسية الخطية كافية لتمثيل مجموعة المستويات
  3. الطريقة تتجنب إعادة بناء الشبكة ويمكنها إنشاء ثقوب تلقائيًا
  4. تم التحقق من فعالية الطريقة على مشاكل معيارية كلاسيكية

القيود

  1. مشكلة الحد الأدنى المحلي: قد تتقارب تكوينات معاملات مختلفة إلى حدود محلية مختلفة
  2. معالجة قفزات المادة: لا تظهر الدوال الأساسية من الرتب العليا ميزة واضحة في تمثيل قفزات واجهة المادة-الفراغ
  3. حساسية المعاملات: معاملات التصفية مثل γ لها تأثير مهم على النتائج
  4. قيود ثنائية الأبعاد: تم التحقق حاليًا فقط في مشاكل ثنائية الأبعاد

الاتجاهات المستقبلية

  1. دراسة متعمقة لتمثيل قفزات واجهة المادة من الرتب العليا
  2. التوسع إلى مشاكل ثلاثية الأبعاد
  3. تحسين الطوبولوجيا متعدد المواد
  4. شروط أكثر تعقيدًا مثل قيود الإجهاد

التقييم المتعمق

المميزات

  1. ابتكار الطريقة قوي: دمج عضوي لثلاث تقنيات متقدمة لحل مشكلة هندسية عملية
  2. أساس نظري متين: اشتقاق المشتقات الطوبولوجية دقيق، الصيغ الرياضية كاملة
  3. تصميم تجريبي معقول: دراسة منهجية لتكوينات رتب مختلفة، مقارنة شاملة
  4. قيمة هندسية عملية عالية: تجنب إعادة بناء الشبكة، مناسبة للتطبيقات العملية

أوجه القصور

  1. غياب التوسع ثلاثي الأبعاد: التحقق فقط في ثنائي الأبعاد، معظم التطبيقات الهندسية ثلاثية الأبعاد
  2. إرشادات معاملات غير كافية: نقص التحليل المتعمق لاختيار معاملات حاسمة مثل معاملات التصفية
  3. غياب تحليل الكفاءة الحسابية: عدم مقارنة تكاليف الحساب مع الطرق التقليدية
  4. التحقق من الهندسة المعقدة غير كافٍ: الأمثلة نسبيًا بسيطة، نقص التحقق من الهياكل الهندسية المعقدة

التأثير

  1. المساهمة الأكاديمية: توفير مسار تقني جديد لمجال تحسين الطوبولوجيا
  2. القيمة الهندسية: من المتوقع تطبيقها في تصميم الهياكل في مجالات الفضاء والملاحة والسيارات وغيرها
  3. قابلية الاستنساخ: بناءً على كود GeoPDEs مفتوح المصدر، يسهل الاستنساخ والتوسع

السيناريوهات المناسبة

  1. تحسين الهياكل المستمرة: مناسبة بشكل خاص لمشاكل التصميم التي تتطلب تغييرات طوبولوجية متكررة
  2. مرحلة التصميم المفاهيمي: استكشاف سريع لتكوينات طوبولوجية مختلفة
  3. تطبيقات متطلبات دقة عالية: الاستفادة من النعومة من الرتب العليا لمنحنيات B-spline
  4. التصميم الآلي: تحسين تلقائي دون الحاجة إلى تعريف يدوي للثقوب الأولية

المراجع

تستشهد الورقة بـ 34 مرجعًا مهمًا، تغطي الأعمال الكلاسيكية والمتقدمة في المجالات الأساسية مثل تحسين الطوبولوجيا والتحليل الهندسي المتساوي وطريقة مجموعة المستويات، مما يوفر أساسًا نظريًا متينًا للبحث. تشمل المراجع الرئيسية طريقة المشتقات الطوبولوجية لـ Amstutz و Andrä (2006) والنظرية الأساسية للتحليل الهندسي المتساوي لـ Hughes وآخرون (2005).