ANSYS'DE SONUÇ DOĞRULAMA VE HATA KAYNAKLARINI AZALTMAK

620 elemanla çözülen bir radyasyon modelinde yüzey net ışınımı 0.2015 W çıkar. Aynı geometri 1600 elemana çıkarıldığında değer 0.20344 W'a oturur. Fark binde bir. Bu, mesh'in yakınsadığının kanıtıdır. Şimdi aynı modeli düşünün: 480 MPa okunan gerilme, mesh sıklaştığında 720 MPa'ya, sonra 1100 MPa'ya tırmanıyor. Bu yakınsama değil; o nokta keskin bir köşe ve modelde matematiksel bir singülarite var. İki sahne de aynı yöntemi kullanıyor, sonuçlar zıt anlam taşıyor.

Doğrulama (verification) ile geçerleme (validation) ANSYS çıktısını rapor seviyesine çıkaran iki ayrı kapıdır. Verification "denklemleri doğru çözdüm mü" sorusudur; mesh, sayısal yöntem, çözücü ayarlarıyla ilgilenir. Validation "doğru denklemleri çözdüm mü" sorusudur; fiziksel deneyle veya analitik formülle karşılaştırma ister. ANSYS Workbench'te Calculate Stress çıktı kontrolünü açtığınızda görünen Convergence aracı sadece ilk kapının bir kısmını otomatikleştirir; geri kalanı modelleyenin disiplinine kalır.

V&V çerçevesinde ne doğrulanır?

Bir TÜBİTAK MAM kalibrasyon laboratuvarında bir transducer test edilirken üç farklı seviyeden geçer: matematiksel model doğruluğu, sayısal çözüm doğruluğu ve fiziksel gerçeklikle eşleşme. ANSYS analizi de aynı üç katmandan geçer. İlk katmanda yanlış denklem seçimi var mı bakılır — küçük deformasyon varsayımı geçerli mi, lineer malzeme yeterli mi. İkinci katmanda mesh ve çözücü ayarları sorgulanır. Üçüncü katmanda sonuç bir kıyas değeriyle karşılaştırılır: analitik formül, daha önce yapılmış deney, üreticinin yayınladığı katalog değeri.

Üçüncü katman üniversitelerin FEM dersinde en çok atlanan yerdir. Öğrenci konsol kirişi çözer, ANSYS bir sayı düşer, raporu yazar. Ucundaki deplasmanın δ = FL³/(3EI) formülüyle iki saniyede tahmin edilebileceğini kimse hatırlamaz. Sayı 5 katı yanlışsa bile fark edilmez. Endüstride ise bu refleks bir alışkanlıktır: çözmeden önce mertebe tahmini, çözdükten sonra karşılaştırma.

Yakınsama çalışması: h ve p yöntemleri

Sonlu elemanlar yöntemi yaklaşım üretir. Mesh sıklaştıkça çözüm gerçek değere yakınsamalı. Bu yakınsamayı iki yoldan elde edebilirsiniz:

• h-refinement: Eleman boyutunu küçültür, sayısını artırırsınız. ANSYS Mechanical'ın Adaptive Convergence aracı bunu otomatikleştirir.
• p-refinement: Eleman sayısını sabit tutup interpolasyon polinomunun derecesini yükseltirsiniz. Lineer elemandan kuadratik elemana geçmek tipik p-refinement adımıdır.
• hp-refinement: İkisini birlikte kullanır. Gerilme yığılması olan bölgede hem eleman küçülür hem polinom derecesi yükselir.

Pratik bir yakınsama testi şöyle kurulur: aynı modeli mesh yoğunluk faktörü 1, 2 ve 4 olacak şekilde üç kez çözün. İlgilendiğiniz büyüklüğü (maksimum gerilme, ucun deplasmanı, reaksiyon kuvveti) her seferinde kayda alın. İki ardışık değer arasındaki bağıl fark %5'in altına düştüğünde mesh yakınsadı sayılır. %30 fark devam ediyorsa ya yakınsama gerçekleşmedi ya da gerilme yığılması yapay bir singülaritedir.

Richardson ekstrapolasyonu üç sonucu kullanarak teorik "sonsuz sıklıkta mesh" değerini tahmin eder. Kuadratik elemanlarla 2x sıklaştırmada ekstrapolasyon faktörü yaklaşık 1.14'tür. Yöntem matematiksel olarak güçlü, ancak iki ölçüm yerine üç noktayla ampirik trend görmek daha güvenlidir. Akademik makalelerde her ikisi sıkça birlikte raporlanır.

Hata hangi kaynaklardan geliyor?

Mesh sıklaştırarak gidermediğiniz dört ana hata kategorisi vardır. Yakınsama testi sadece birinciyi kapatır; geri kalan üçü kalıcı şekilde sonuçta kalır ve sayfalarca rapor sonrası fark edilir.

1. Diskretizasyon hatası: Mesh kaynaklı. Sıklaştırma ve eleman tipini doğru seçmekle azalır. Kuadratik elemanda 2 eleman kalınlık bile lineer elemanlı 4 elemana göre daha az hata verir; literatürde lineer için %6.5 mertebesinde hatalar tipiktir.
2. Model hatası: Yanlış fizik seçimi. Küçük deformasyon varsayımı altında büyük rotasyon çözmek, viskoelastik davranışı lineer elastik modellemek, dinamik problemi statik çözmek bu kategoride. Mesh ne kadar sıklaşsa da çözüm gerçeğe yaklaşmaz.
3. Sınır koşulu hatası: Fixed Support'un abartılı kullanımı en yaygın olanı. Tüm yüzeyi sabit yapmak gerilmeyi 2-3 katına şişirir; cıvatalı bağlantı gerçekte hiçbir yönde mutlak sıfır değildir. Yerine cıvata bölgesini ayrı tutarak pretension veya Remote Displacement ile gerçek esnekliği yakalayın.
4. Temas hatası: Connections panelinde otomatik tanınan kontaklar her zaman doğru tipte değildir. Bonded sandığınız bir temas aslında frictional olabilir; iki parça arasında olmaması gereken yapışma da olabilir. Otomatik tanıma çıktısını her zaman görsel olarak doğrulayın.

Türk savunma sanayi simülasyon ekipleri raporlama disiplininde bu dört kategoriyi ayrı ayrı kontrol listesi yapar. ASELSAN tedarikçi geçerleme süreçlerinde mesh yakınsama belgesi tek başına yeterli sayılmaz; sınır koşulu doğrulama çıktısı, temas listesi denetim raporu ve analitik karşılaştırma birlikte istenir. Bir ANSYS eğitim programı kapsamında ilerlerken bu dört kategoriyi ayrı bloklar olarak düşünmek alışkanlığa dönüşür.

Mesh kalite metrikleri: gözle değil sayıyla

Güzel görünen mesh her zaman iyi mesh değildir. ANSYS Solution Information panelinde otomatik geçilen kalite uyarıları farklı önemde olur. Hangisi neyi söylüyor:

• Element Quality: 0-1 aralığı. 0.7 üstü tercih, 0.2 altı kırmızı bayrak.
• Skewness: 0 mükemmel, 0.95 üstü diverjans riski. Akış analizinde belirleyici.
• Aspect Ratio: Yapısalda 10 üstü problemli; akışta inflation katmanları için 20'ye kadar tolere edilir.
• Orthogonal Quality: 0.15 altı ciddi sorun göstergesi, akışta öncelikli metrik.
• Jacobian Ratio: Yüksek dereceli elemanlarda 30 üstü deformasyon işareti.

İlgilendiğiniz bölgeye Sphere of Influence ile yerel sıklaştırma uygulayın; tüm modeli aşırı sıklaştırmak hem süreyi uzatır hem Jacobian'ı bozar. Akış analizinde inflation layer büyüme oranı 1.15-1.25 bandında tutulur; ani büyüme y+ değerini kaçırır.

Sınır koşulu ve yük nasıl doğrulanır?

Doğru meshli bir modelin yanlış sınır koşuluyla verdiği sonuç, kötü meshli doğru kurulumdan daha tehlikelidir. Çünkü gözle bakılınca şıktır. Sık karşılaşılan tuzaklar şunlar:

Kuvvet yönü kontrolü atlanır. Workbench koordinatı global iken parça döndürülmüş olabilir. Yükü her zaman görsel okla doğrulayın; sayı tek başına aldatıcı. Basınç işaret kuralı: dışa doğru pozitif, içe doğru negatif. Bu işaret hatası genç mühendislerin en çok düzeltilen yanlışıdır.

Tepki kuvveti testi en hızlı denetleyicidir. Yapıya 1000 N uyguladıysanız mesnetlerdeki toplam Force Reaction 1000 N olmalıdır. 850 N veya 1200 N geliyorsa bir yerde temas tanımı kaçırılmış, yük doğrultusu yanlış girilmiş veya model rijit bir parça üzerinden boşa zıplıyor demektir. Bu test, başka türlü sayfalarca rapor sonrası fark edilecek temel modelleme hatalarını saniyeler içinde yakalar.

Benchmark karşılaştırması ve sanity check

Doğrulamanın en güçlü adımı dış referansla karşılaştırmadır. Üç tip benchmark kullanılır:

• Analitik formül: Konsol kiriş için δ = FL³/(3EI), kalın silindirik basınç kabı için Lamé denklemleri, dönen disk için Yoshimoto. Mertebe kontrolüne yeter.
• Standart test case: NAFEMS benchmark setleri, ASME V&V dokümanları, ANSYS Verification Manual problemleri. Bunlar bilinen tam çözümlü modeldir.
• Fiziksel deney: Strain gauge ölçümü, ivmeölçer testi, kalibrasyon laboratuvar raporu. Validation'ın altın standardı.

Sanity check (akıl sağlığı kontrolü) çözüm öncesinde yapılır: 100 kg yük altında çelik bir kirişin mertebe deplasmanı milimetrenin onda biri civarındaysa, ANSYS 50 cm gösteriyorsa bir yerde temel bir hata var. Birim sistemi yanlış, kuvvet 1000 katı büyük girilmiş veya mesnet eksik. Bu refleks tek bir alışkanlıkla kurulur: çözmeden önce zihninde bir tahmin oluştur. Tahmin doğru çıkarsa güven artar, yanlış çıkarsa hangi sebepten saptığını sorarsın.

Sonuç yorumlarken hangi tuzaklar var?

Equivalent (von-Mises) Stress sünek malzemeler için doğru kıyas büyüklüğüdür; çelik, alüminyum, paslanmaz. Dökme demir, cam, beton, seramik gibi gevrek malzemelerde Maximum Principal Stress kullanılır çünkü çekme yönündeki ana gerilme kırılmayı belirler. İkisini karıştırmak yorum hatasının en pahalısıdır.

Total Deformation sıfır kontrolünü kolaylaştırır: 1 nm seviyesinde değer çıkıyorsa muhtemelen modeliniz fix'ten kilitlenmiş, kuvvet eksik veya birim yanlış. Lokal bir bölgede 800 MPa görünüyor ve etrafı 60 MPa ise yüksek olasılıkla tekil gerilme noktası vardır. Kontrol için o bölgenin etrafına yerel sıklaştırma yapın, mesh'i bir kademe daha sık çözün. Sayı %5'in altında değişiyorsa çözüm yakınsamış; %30'un üstünde değişiyorsa orada anlamsız bir matematiksel singülarite var, fiziksel gerçek değil.

Doğrulanmış kontrol listesi

Çözümü "raporlamaya hazır" sayılan modelin geçtiği son kontrol listesi:

1. Birim sistemi tek ve tutarlı (mm-N-s veya m-kg-s; yarı yolda değişmemiş)
2. Eleman tipi soyutlamaya uygun (kirişe beam, ince saca shell, kalın yapıya solid)
3. Mesh metrikleri eşik içinde (Element Quality > 0.7, Skewness < 0.95)
4. Yakınsama testi üç sıklıkta çalıştırıldı, fark < %5
5. Sphere of Influence ile yerel sıklaştırma yapıldı, singülarite ayıklandı
6. Sınır koşulu görsel okla doğrulandı, Fixed Support gereksiz kullanılmadı
7. Temas listesi gözden geçirildi, otomatik bonded'lar manuel onaylandı
8. Force Reaction = Uygulanan Yük (denge sağlandı)
9. Sonuç tipi malzemeye uygun (sünek için von-Mises, gevrek için Maximum Principal)
10. Analitik formül veya benchmark ile aynı mertebede

Bu liste tek başına teknik bir kontrolden öte, sonuca güvenmenin yapısal halidir. Yıllık iki yüz analiz yapan bir mühendis kendini bu adımları otomatik geçer halde bulur; yeni başlayan biri için aynı liste rapor şablonunun arkasına yazılmış bir hatırlatma kâğıdıdır. Disiplin tekrarla içselleşir.

ANSYS güçlü bir araç ama doğrulanmamış bir sayı sadece kâğıttaki bir karaktirler dizisidir. Güveni kuran araç değil, kullanıcının sayıyı sorgulama refleksidir. Mesh yakınsama testini standartlaştıran, sınır koşulunu görsel doğrulayan, tepki kuvvetini her seferinde kontrol eden ve sonucu analitik bir tahminle yan yana koyan ekip — sonuçların güvenilirliğinde belirgin bir fark yaratır.