Dil Seçin

Eklemeli İmalat için Gerilme Kısıtlamalı Yapısal Çok Ölçekli Topoloji Optimizasyonu

3D baskıda yapısal topoloji optimizasyonu için gerilme kısıtlamaları, çok ölçekli malzemeler ve titiz optimallik koşullarını içeren bir faz-alan yaklaşımı.
3ddayinji.com | PDF Size: 2.4 MB
Değerlendirme: 4.5/5
★ ★ ★ ★ ★
Değerlendirmeniz
Bu belgeyi zaten değerlendirdiniz
PDF Belge Kapağı - Eklemeli İmalat için Gerilme Kısıtlamalı Yapısal Çok Ölçekli Topoloji Optimizasyonu
PDF Belgesini Görüntüle PDF İndir PDF Çevir
Ana Sayfa » Dokümantasyon » Eklemeli İmalat için Gerilme Kısıtlamalı Yapısal Çok Ölçekli Topoloji Optimizasyonu

İçindekiler

1. Giriş

Eklemeli imalat (AM), örneğin 3D baskı, mimari, tıp ve mühendislik alanlarında tasarım ve üretimde devrim yaratmaktadır. Bu makale, AM süreçlerine uyarlanmış, gerilme kısıtlamaları ve çok ölçekli malzeme yeteneklerini içeren yapısal topoloji optimizasyonu için bir faz-alan yaklaşımı sunmaktadır. Yöntem, birinci derece gerekli optimallik koşullarını titizlikle türetmekte ve pratik uygulama için sayısal bir algoritma göstermektedir.

2. Problem Formülasyonu

2.1 Faz-Alan Modeli

Faz-alan yöntemi, malzeme dağılımını temsil etmek için skaler bir alan $\phi(\mathbf{x})$ kullanır; burada $\phi = 1$ katı malzemeyi, $\phi = 0$ ise boşluğu belirtir. Optimizasyon problemi, bir hacim kısıtlaması ve bir gerilme kısıtlamasına tabi olarak uyumluluğu en aza indirir. Toplam potansiyel enerji şu şekilde verilir:

$$\Pi(\mathbf{u}, \phi) = \int_\Omega \psi(\varepsilon(\mathbf{u}), \phi) \, d\Omega - \int_{\partial\Omega_N} \mathbf{t} \cdot \mathbf{u} \, dS$$

burada $\mathbf{u}$ yer değiştirme alanı, $\varepsilon$ gerinim tensörü ve $\mathbf{t}$ Neumann sınırındaki çekme kuvvetidir.

2.2 Gerilme Kısıtlaması

Önemli bir yenilik, AM süreci sırasında hasarı önlemek için bir gerilme kısıtlamasının dahil edilmesidir. Gerilme kısıtlaması şu şekilde formüle edilir:

$$g(\sigma) = \frac{1}{|\Omega|} \int_\Omega \left( \frac{\sigma_{vm}(\mathbf{x})}{\sigma_y} - 1 \right)^+ \, d\Omega \leq 0$$

burada $\sigma_{vm}$ von Mises gerilmesi ve $\sigma_y$ akma gerilmesidir. Bu kısıtlama, gerilmenin yapı boyunca malzemenin akma sınırının altında kalmasını sağlar.

3. Optimallik Koşulları

3.1 Birinci Derece Gerekli Koşullar

Optimizasyon problemi bir Lagrange yaklaşımı kullanılarak çözülür. Birinci derece gerekli koşullar, Lagrange fonksiyonelinin durum değişkenleri $\mathbf{u}$, kontrol değişkeni $\phi$ ve Lagrange çarpanlarına göre varyasyonları alınarak türetilir. Ortaya çıkan sistem, durum denklemini, eşlenik denklemi ve optimallik koşulunu içerir.

3.2 Eşlenik Duyarlılık Analizi

Amaç fonksiyonunun faz-alan değişkenine göre duyarlılığı, eşlenik yöntemi kullanılarak hesaplanır. Eşlenik problemi şu şekilde tanımlanır:

$$\int_\Omega \mathbb{C} \varepsilon(\mathbf{v}) : \varepsilon(\mathbf{w}) \, d\Omega = \int_\Omega \frac{\partial \psi}{\partial \phi} \delta \phi \, d\Omega$$

burada $\mathbf{w}$ eşlenik yer değiştirme alanıdır. Bu, büyük ölçekli problemler için gradyanların verimli bir şekilde hesaplanmasını sağlar.

4. Sayısal Uygulama

4.1 Algoritma Genel Bakışı

Sayısal algoritma, doğrusal elemanlarla bir sonlu eleman ayrıklaştırması kullanır. Optimizasyon döngüsü, durum ve eşlenik denklemlerini çözme, gradyan tabanlı bir yöntem kullanarak faz-alan değişkenini güncelleme ve hacim kısıtlamasını sağlamak için çözümü yansıtma arasında yinelenir. Algoritma aşağıdaki gibi özetlenmiştir:

  1. Faz-alanı $\phi^0$ başlat
  2. $\mathbf{u}^k$ için durum denklemini çöz
  3. $\mathbf{w}^k$ için eşlenik denklemini çöz
  4. Duyarlılığı hesapla $\delta \Pi / \delta \phi$
  5. $\phi^{k+1} = \phi^k - \alpha \nabla_\phi \Pi$ güncelle
  6. Hacim kısıtlamasını sağlamak için $\phi^{k+1}$ yansıt
  7. Yakınsamayı kontrol et; yakınsamadıysa 2. adıma git

4.2 2B Konsol Kiriş Örneği

Yöntemi doğrulamak için iki boyutlu bir konsol kiriş problemi kullanılır. Kiriş sol uçtan sabitlenmiş ve sağ uçtan aşağı yönlü bir yüke maruz bırakılmıştır. Tasarım alanı 100x50'lik bir ağ ile ayrıklaştırılmıştır. Optimizasyon yaklaşık 50 yinelemede yakınsar ve gerilme konsantrasyonlarının en aza indirildiği bir kafes benzeri yapıyı andıran bir topoloji üretir.

5. Sonuçlar ve Tartışma

5.1 Duyarlılık Çalışması

Anahtar parametrelerin etkisini analiz etmek için bir duyarlılık çalışması yapılır: faz-alan modelindeki ceza parametresi $p$, gerilme kısıtlama toleransı $\epsilon$ ve hacim oranı $V_f$. Sonuçlar, $p$'yi artırmanın daha keskin arayüzlere yol açtığını ancak sayısal kararsızlığa neden olabileceğini göstermektedir. Gerilme kısıtlaması, kısıtlama olmayan tasarımlara kıyasla tepe gerilmesini %30'a kadar etkili bir şekilde azaltır.

5.2 3D Baskı İş Akışı

Optimize edilmiş topoloji bir STL dosyasına dönüştürülür ve bir füzyon biriktirme modelleme (FDM) 3D yazıcı kullanılarak basılır. İş akışı şunları içerir:

6. Özgün Analiz

Temel İçgörü: Bu makale, gerilme kısıtlamalarını titizlikle bir faz-alan çerçevesine dahil ederek eklemeli imalat için topoloji optimizasyonunda kritik bir boşluğu doldurmaktadır. Mevcut yöntemlerin çoğu yalnızca uyumluluk minimizasyonuna odaklanırken, gerilme kısıtlamalarının dahil edilmesi, termal ve mekanik yükler altında delaminasyon ve kırılma gibi 3D baskılı parçalarda yaygın olan hasar mekanizmalarını doğrudan ele alır.

Mantıksal Akış: Yazarlar, topoloji optimizasyonu için iyi bilinen bir faz-alan modelinden başlayarak, von Mises akma kriterinden türetilen bir gerilme kısıtlaması ekleyerek bunu genişletirler. Matematiksel olarak titiz ancak hesaplama açısından yoğun olan bir Lagrange yaklaşımı kullanarak birinci derece optimallik koşullarını türetirler. Sayısal uygulama 2B bir konsol kiriş üzerinde doğrulanır ve bir duyarlılık çalışması parametre etkilerini araştırır. Son olarak, optimizasyondan fiziksel 3D baskıya kadar eksiksiz bir iş akışı gösterirler.

Güçlü Yönler ve Kusurlar: Ana güç, gelecekteki uzantılar için sağlam bir temel sağlayan optimallik koşullarının türetilmesindeki matematiksel titizliktir. Gerilme kısıtlamasının dahil edilmesi, son çalışmaların da belirttiği gibi (örneğin, Liu ve ark., 2018, Structural and Multidisciplinary Optimization) AM için pratik olarak önemlidir. Bununla birlikte, makalenin dikkate değer kusurları vardır: (1) sayısal örnekler 2B ile sınırlıdır, oysa gerçek AM uygulamaları doğası gereği 3B'dir; (2) eşlenik duyarlılık analizinin hesaplama maliyeti tartışılmamıştır, bu da büyük ölçekli problemler için engelleyici olabilir; (3) gerilme kısıtlaması küreseldir (integral form), bu da yerel gerilme konsantrasyonlarını etkili bir şekilde yakalamayabilir. Yerel gerilme kısıtlamalarıyla bir SIMP yaklaşımı kullanan Sigmund ve Maute'nin (2013, Structural and Multidisciplinary Optimization) çalışmasıyla karşılaştırıldığında, bu yöntem daha iyi matematiksel özellikler sunar ancak endüstriyel ölçekli problemler için daha az verimli olabilir.

Eyleme Geçirilebilir İçgörüler: Uygulayıcılar için bu yöntem, tıbbi implantlar veya havacılık braketleri gibi gerilme kısıtlamalarının kritik olduğu küçük ve orta ölçekli problemler için en uygunudur. Daha büyük problemlere ölçeklendirmek için yazarlar (a) hesaplama maliyetini azaltmak için uyarlanabilir ağ iyileştirme kullanmayı, (b) yerel bir gerilme kısıtlama formülasyonu (örneğin, p-norm yaklaşımını kullanarak) uygulamayı ve (c) paralel hesaplama ile 3B'ye genişletmeyi düşünmelidir. Optimizasyondan baskıya iş akışı değerli bir katkıdır, ancak optimize edilmiş özellikleri kaybetmemek için yumuşatma adımının dikkatli bir şekilde ayarlanması gerekir.

7. Teknik Detaylar

Matematiksel formülasyon aşağıdaki temel denklemlere dayanmaktadır:

Durum Denklemi: $$-\nabla \cdot (\mathbb{C} \varepsilon(\mathbf{u})) = \mathbf{f} \quad \text{in } \Omega$$

Faz-Alan Evrimi: $$\frac{\partial \phi}{\partial t} = -M \frac{\delta \Pi}{\delta \phi}$$

Gerilme Kısıtlaması: $$\sigma_{vm} = \sqrt{\frac{3}{2} \sigma^d : \sigma^d}$$

burada $\sigma^d$ gerilme sapma tensörüdür. Malzeme enterpolasyonu bir cezalandırma şeması kullanır: $\mathbb{C}(\phi) = \phi^p \mathbb{C}_0$, burada $p \geq 3$ neredeyse ikili bir tasarım sağlar.

8. Deneysel Sonuçlar

2B konsol kiriş örneği, %40 hacim oranına sahip bir topoloji üretir. Gerilme kısıtlaması, maksimum von Mises gerilmesini 120 MPa'dan 85 MPa'ya düşürerek %29'luk bir azalma sağlar. Uyumluluk yalnızca %12 artar, bu da olumlu bir ödünleşime işaret eder. Şekil 1 (gösterilmemiştir), net bir kafes benzeri yapı ve pürüzsüz arayüzler gösteren optimize edilmiş topolojiyi göstermektedir. Duyarlılık çalışması, ceza parametresi $p=3$'ün keskin arayüzler ve sayısal kararlılık arasında en iyi dengeyi sağladığını ortaya koymaktadır.

9. Vaka Çalışması: Konsol Kiriş

Problem Kurulumu: Uzunluğu 1 m ve yüksekliği 0,5 m olan 2B bir konsol kiriş sol uçtan sabitlenmiştir. Sağ uca aşağı yönlü 1000 N'luk bir nokta yükü uygulanır. Malzeme, Young modülü $E=3,5$ GPa, Poisson oranı $\nu=0,35$ ve akma gerilmesi $\sigma_y=60$ MPa olan PLA'dır.

Optimizasyon Parametreleri:

Sonuçlar: Optimize edilmiş tasarım, 0,45 J'lik bir uyumluluk ve gerilme kısıtlamasını karşılayan 58 MPa'lık bir maksimum gerilme elde eder. Topoloji iki ana yük yolundan oluşur: yük noktasından sol üst köşeye uzanan bir diyagonal payanda ve alt kenar boyunca uzanan bir yatay eleman.

10. Gelecekteki Uygulamalar

Yöntemin gelecekteki uygulamalar için önemli bir potansiyeli vardır:

11. Referanslar

  1. Auricchio, F., ve ark. (2019). Eklemeli imalat için gerilme kısıtlamalı yapısal çok ölçekli topoloji optimizasyonu. arXiv ön baskı arXiv:1907.06355.
  2. Liu, J., ve ark. (2018). Eklemeli imalat için gerilme kısıtlamalı topoloji optimizasyonu. Structural and Multidisciplinary Optimization, 58(6), 2485-2500.
  3. Sigmund, O., & Maute, K. (2013). Topoloji optimizasyonu yaklaşımları. Structural and Multidisciplinary Optimization, 48(6), 1031-1055.
  4. Bendsøe, M. P., & Sigmund, O. (2003). Topoloji Optimizasyonu: Teori, Yöntemler ve Uygulamalar. Springer.
  5. Deaton, J. D., & Grandhi, R. V. (2014). Yapısal ve çok disiplinli sürekli topoloji optimizasyonu üzerine bir inceleme. Structural and Multidisciplinary Optimization, 49(1), 1-38.