ANSYS'DE MESH KALİTESİ KRİTERLERİNİ DOĞRU KURGULAMAK

Bir savunma sanayi yan sanayisinde çalışan genç bir analiz mühendisi, çelik bir tutucu parçanın statik çözümünü kuruyor. İlk koşumda maksimum gerilme 240 MPa görünüyor; aynı parçayı bir gün sonra daha sık mesh ile çözdüğünde sonuç 410 MPa. Aradaki fark malzeme, sınır koşulu veya yük değil — sadece mesh. ANSYS çözücüsünün matematiksel gücü ne kadar yüksekse, ona verdiğiniz elemanların geometrik kalitesi de o kadar belirleyici hale gelir.

Mesh, sürekli bir cismi sonlu elemanlar yöntemine tercüme eden köprüdür. Köprü çarpık veya gerilmiş elemanlarla kurulduğunda çözücü bir sonuç üretmeye devam eder; ama o sonucun fiziksel gerçekle bağı zayıflar. Türk üniversitelerinde Makine Bölümü dördüncü sınıfta verilen FEM derslerinin pratik kısmı bu yüzden mesh kontrolüne yoğunlaşır; ANSYS eğitimini ciddi anlamda tamamlamış mühendisler için mesh kalite metrikleri raporu çözüm başlatmadan önce kontrol edilen ilk listedir.

Hangi element tipi hangi geometriye uyar?

ANSYS Meshing dört temel eleman tipi üretir: tetrahedral (tetra), hexahedral (hex), wedge (prizmatik) ve pyramid. Her birinin matematiksel uyum penceresi farklıdır. Aynı düğüm sayısıyla karşılaştırıldığında hexahedral eleman, doğru yönlendirildiğinde tetrahedrale göre belirgin biçimde daha doğru sonuç verir; ama karmaşık bir bağlantı parçası veya dökme gövde için otomatik hex üretimi pratik değildir.

Tipik tercih mantığı şu şekilde işler:

• Hex (tuğla): Düzgün gövdeler, kanallar, prizmatik kesitler için Sweep veya MultiZone yöntemiyle üretilir. CFD'de boundary layer çözünürlüğü için tercih edilir.
• Tetra: Karmaşık dökme parçalar, organik geometriler. Otomatik üretimi en hızlı, geometrik esnekliği en yüksek tipdir.
• Wedge: İnflasyon katmanı uygulamasında duvar yakınında tetra üzerine kurulan prizmatik tabaka. Sınır tabakası gradyanlarını yakalar.
• Pyramid: Hex ve tetra bölgeleri birbirine bağlamak için geçiş elemanı olarak kullanılır.

Bir TÜBİTAK MAM ar-ge projesinde tipik karar şudur: gövdenin ana bölümü hex ile, bağlantı detayları tetra ile, ikisinin arası pyramid geçişiyle çözülür. Bu hibrit yaklaşım hem hesap maliyetini düşürür hem de detay bölgelerin geometrik uyumunu korur.

Mesh kalite metriklerinin eşik değerleri nedir?

ANSYS, mesh kalitesini ölçen birden fazla metrik sunar. Hepsi aynı şeyi söylemez: biri elemanın açı bozulmasına, biri kenar oranına, biri hacim-yüzey ilişkisine bakar. Mesh detay panelindeki Mesh Metric açılır listesinden seçilen metrik histogram olarak ekrana düşer. Hangi metriğe ne zaman bakacağını bilmek pratik bir alandır; matematiksel arka planı tek tek anlamak için ANSYS'in kendi eğitim merkezi kaynakları başlangıç düzeyini hızla aşmanızı sağlar.

Üretim ortamında bakılan beş ana metrik:

1. Element Quality: 0-1 arası, 1 mükemmel. Tek metrik bakılacaksa bu. Statik yapısalda 0.1 altı eleman çözücü uyarısı tetikler; 0.2 üstü hedeflenir.
2. Skewness: Elemanın ideal şekilden açısal sapması. Tetrahedral mesh için maksimum 0.95, ortalama 0.33 altı yaygın bir CFD eşiğidir. Hexahedral için 0.5 altı çok iyi, 0.8 altı kabul edilebilir.
3. Aspect Ratio: En uzun kenar / en kısa kenar. Yapısal analiz için 20-40 üst sınırdır; CFD'de inflasyon katmanları nedeniyle yüksek aspect ratio normal sayılır (5-10 tipik).
4. Jacobian Ratio: Elemanın gerçek şeklinin ideal şekle dönüşüm jakobiyenine oranı. 1 ideal; %90 elemanın 1-10 bandında kalması iyi mesh kuralıdır. Negatife düşerse eleman ters dönmüştür, çözücü durur.
5. Orthogonal Quality: 0-1 arası, 1 ideal. CFD'de kritik metrik; 0.15 altı eleman Fluent veya CFX'te yakınsamayı bozar. Yapısal analizde 0.1 minimum eşiktir.

Pratik bir yaklaşım: tek bir min/max değer yerine dağılım histogramına bak. 50 bin tetrahedrali olan bir modelde 20 elemanın skewness'i 0.93 olsa bile geri kalanının %95'i 0.5 altındaysa mesh kullanılabilir. Asıl soru şu: yüksek skewness'li elemanlar yük yolu üzerinde mi, yoksa yüksüz bir bölgede mi? Yük yolundaysa o bölge yeniden meshlenir; değilse kalabilir.

Geometri hazırlığı: temiz mesh temiz geometriden çıkar

Mesh hatalarının önemli bir kısmı aslında geometriden kaynaklanır. Otomotiv yan sanayisinde NVH analizi yapan bir mühendis, bir tutucu parçanın çözümünde yüksek skewness'li yüzlerce mikro tetra ile karşılaşıyorsa, bunların çoğu büyük olasılıkla CAD'den gelen 0.3-0.5 mm fillet'lerin yarattığı problemdir. ANSYS bu detayları temizlemek için iki ortam sunar: SpaceClaim hızlı tamir, DesignModeler parametrik temizleme.

Tipik temizleme adımları:

• Merge Faces: 1 mm altı yüzeyleri birleştir; aksi halde o bölgede mikro elemanlar üreyip aspect ratio'yu patlatır.
• Fillet/pah simplification: Yük yolundan uzak 0.5-1 mm radyuslar hesaba katkı vermez ama eleman sayısını iki katına çıkarır.
• Imprint Faces: Yük ve mesnet uygulanacak yüzeyler partition edilir; sınır koşulu o yüzeye temiz oturur.
• Beam-shell-solid ayrımı: Profil kesitler için beam, ince saclar için shell, kalın gövde için solid. Yanlış idealizasyon her metrik kontrolünü çürütür.
• Named Selection: Yüzey ve gövdeler isimlendirilir; sizing ve sınır koşulu atamaları kopukluk olmadan ilerler.

Türk savunma sanayisinde tipik bir senaryo: bir radar montaj braketinin statik analizinde brakete giren cıvata deliklerinin çevresinde countersink fillet'leri vardır. Bu fillet'ler yük yolundan uzaksa SpaceClaim'de tek tıkla kaldırılır; üzerinde dururlarsa Sphere of Influence ile yerel mesh refinement uygulanır.

Refinement ve inflasyon katmanı nasıl kurulur?

Mesh sıklığını her yerde aynı tutmak hem hesap maliyetini hem zamanı boşa harcar. Akıllı yaklaşım kritik bölgede yerel refinement uygulamaktır. ANSYS bunun için birkaç araç sunar: Face Sizing, Edge Sizing, Body Sizing, Sphere of Influence ve Refinement (1-3 seviye).

CFD analizinde duvar yakınında inflasyon katmanı şarttır. Sınır tabakası gradyanları çözümün en hassas bölgesidir ve burada hexahedral wedge elemanlar gerekir. Tipik inflasyon parametreleri:

• İlk katman yüksekliği: Hedeflenen y+ değerine göre belirlenir. 1 m/s akış hızında 0.005 m ilk katman y+ ~ 50 değerini verir; düşük Re modeli için y+ < 1 hedeflenirse ilk katman mikrometre mertebesine düşer.
• Büyüme oranı: 1.1-1.2 arası standart; ani geçişler eleman aspect ratio'sunu bozar.
• Katman sayısı: 8-15 katman tipik. Reynolds sayısı ve duvar yasası yaklaşımına göre değişir.
• Transition ratio: İnflasyon katmanından dış mesh'e geçişin yumuşaklığı. 0.272 default değer çoğu durumda iyi çalışır.

Yapısal analizde ise refinement kararı genelde gerilme yığılması olan bölgeye verilir. Bir bağlantı deliğinin çevresinde, kaynak dikişinin başlangıç noktasında veya keskin köşe yakınında 2-3 seviye Refinement uygulamak doğal bir reflekstir.

Sınır koşulu ve sonuç yorumlama

Mesh ne kadar iyi olursa olsun, sınır koşulu yanlış kuruluysa sonuç da yanlıştır. Tipik hata: mesnet yüzeyini tüm alt yüzeyle değil, gerçek bağlantı bölgesiyle tanımlamak. Cıvata bağlantılı bir taban Imprint Faces ile cıvata deliklerinin etrafındaki ring'lere ayrılıp Cylindrical Support uygulandığında çok daha gerçekçi davranır.

Yük tanımında yön ve referans koordinat sistemi netleştirilmelidir. Pressure yüzeye dik etkir, Force yön vektörüne göre uygulanır. Remote Force, geometrik olarak modelde olmayan bir noktadan yük aktarmak için kullanılır — bir titreşim analizinde kütle merkezi tipik kullanım örneğidir.

Çözüm bittiğinde ilk bakılan genelde Equivalent (von Mises) Stress'tir. Bu metrik sünek malzemeler için doğrudur; dökme demir veya seramik gibi gevrek malzemelerde Maximum Principal Stress daha uygundur. Tek bir kırmızı nokta görüp parçanın kırılacağına hükmetmek yaygın bir hatadır; o noktanın geometrik singülarite olup olmadığı mesh sıklaştırma testi ile kontrol edilmelidir.

Mesh yakınsama çalışması

Profesyonel bir analizin kapanışı mesh independence study'dir. Aynı problem üç farklı sıklıkta meshlenip ilgilenilen değerin (maksimum gerilme, kritik nokta deplasmanı, basınç düşüşü) nasıl değiştiği izlenir. Son iki koşum arasındaki fark %5 altına düştüyse mesh yakınsamış sayılır. Bir TÜBİTAK projesinin kalite raporunda yakınsama tablosu standart bir kapanıştır.

Tipik bir tablo şöyle görünür:

• Mesh 1 (kaba): 50 bin eleman, max gerilme 380 MPa, çözüm süresi 4 dk
• Mesh 2 (orta): 150 bin eleman, max gerilme 410 MPa, çözüm süresi 12 dk
• Mesh 3 (sık): 400 bin eleman, max gerilme 418 MPa, çözüm süresi 38 dk

Mesh 2 ile Mesh 3 arasındaki fark %2; yakınsama gerçekleşmiş, Mesh 2 raporlanabilir. MATLAB tarafında otomasyon yapmak isterseniz workbench journal dosyalarını MATLAB eğitimini tamamlamış mühendisler tipik olarak parametrik tarama scriptleriyle entegre eder.

Otomotiv ve savunma sanayisinde pratik notlar

Otomotiv yan sanayisinde NVH analizleri için mesh kalitesi özellikle önemlidir çünkü modal analiz mesh'e gerilme analizinden daha duyarlıdır. Düşük kaliteli elemanlar yapay frekanslar üretebilir; bu da NVH raporunu çürütür. Bursa-Demirtaş OSB'deki bir motor montaj braketi tedarikçisi, ana sanayi tarafından beklenen titreşim spektrumu için mesh global boyutunu 2 mm'ye, kritik bölgede 0.5 mm'ye kadar indirir.

Savunma sanayisinde durum farklıdır. Ankara-İvedik OSB'deki bir radar muhafazası üreticisi, şok yükü altındaki bir kapağın simülasyonunda eksplisit dinamik kullanır ve burada eleman boyutu zaman adımını belirler (CFL koşulu). Çok küçük elemanlar çözüm süresini patlatır; bu yüzden mesh refinement çok daha cerrah bir karardır. Aynı geometri NVH için 1 milyon eleman gerektirebilirken patlama yüklemesinde 200 bin elemanla daha akıllıca çözülür.

Türk üniversitelerinin makine ve havacılık bölümlerinde lisansüstü öğrencilere verilen ANSYS pratik dersi tam da bu farkı oturtmaya çalışır: "iyi mesh" tek bir sayı değil, problemin tipine göre değişen bir karar zinciridir. ANSYS'in verdiği metrikler ham veridir; o veriyi yorumlamak mühendisin işidir.

Mesh kalitesi bir analizin tek başına kazandıracağı bir şey değildir, ama atlandığında her şeyi kaybettirir. Geometri temizliğinden başlayıp metrik histogramına, refinement kararından yakınsama testine kadar uzanan disiplini kuran mühendis, ANSYS'in verdiği sayılara güvenle imza atabilir.